LATE GLACIAL AND HOLOCENE ENVIRONMENTAL CHANGE INFERRED FROM SEDIMENTARY ARCHIVES OF KUSAWA LAKE, BOUNDARY RANGE MOUNTAINS, YUKON TERRITORY, CANADA
by
Nicole Angela Chow
A thesis submitted in conformity with the requirements for the degree of Master of Science Graduate Department of Geography University of Toronto
© Copyright by Nicole Angela Chow (2009)
LATE GLACIAL AND HOLOCENE ENVIRONMENTAL CHANGE INFERRED FROM SEDIMENTARY ARCHIVES OF KUSAWA LAKE, BOUNDARY RANGE MOUNTAINS, YUKON TERRITORY, CANADA
M.Sc., 2009
Nicole Angela Chow
Deptartment of Geography, University of Toronto
Abstract
Modern Kusawa Lake (60° 19' 55” N, 136° 4' 48” W, 142 km 2) of southwestern Yukon
Territory drains a 4290 km 2 catchment, 4.7 % of which is glacier covered. Sediment cores show variability both down-lake and within specific sub-basins of the lake. In Regions II -V of
Kusawa Lake, sediments are mainly clastic with massive to weakly laminated silts and clays interrupted by fine sand units, which reflect distinct runoff events into Region IV from glacier sources. In Region I, massive silts, silt-clay couplets are interrupted by thick sand deposits derived from the Primrose River delta. Further up-lake, the sediment record is further interrupted by modern sediment delivery from the Kusawa Campground alluvial fan.
The relatively small accumulation of lake glacial and Holocene sediment input in Kusawa
Lake is similar to other large lakes of the Canadian Cordillera. These patterns reflect a particular style of deglaciation and Holocene sediment inputs.
ii
Frontispiece: Kusawa Lake, looking south from the outlet. July 2007 (Photograph taken by Joe Desloges)
“Field work can be carried out in all weather conditions provided that you are dressed appropriately for it! ” – Dr. Joe Desloges.
iii
Acknowledgements
First, I would like to thank Dr. Joe Desloges. His passion, enthusiasm and encyclopedic knowledge for geomorphic studies have been inspiring. I would also like to thank Dr. Sarah
Finkelstein for her encouragement and words of wisdom. Both Professors have been a huge support to my academic growth and especially the development of this project. A huge thank you to Mr. Mircea Pilaf, who over looked the splitting of cores, floods, shipments, ceiling collapses, and issues with the photocopier. Life in PGB would be a disaster without him.
I would like to acknowledge the Natural Science and Engineering Research Council of
Canada for their funding. In addition, to Dr. John Westgate and Dr. Mike Gorton in Geology,
Sam Roshdi in Chemical Engineering and Dr. Rolf Beurkens in Physics for access to different laboratory analyses. To Dr. Sharon Cowling, Dr. Tony Davis and Susan Calanza whom have generously provided guidance throughout my time here at U of T.
Muddy times spent in the field with Dr. Bob Gilbert, and Tim Phillpot were certainly fun and memorable. Thank you, Monique Stewart and Dr. Gilbert for the acquisition of the cores and to Krish Chakraborty for the analysis of diatoms.
Fellow graduate PGBers also provided a source of diversion: Carlos Avendaño, Cameron
Balfour, Feng Deng and Lisa Zhang who made lunch hours a multi-cultural (tri-lingual) experience. For times of laughter and advice during those mind-boggling times: Jen Adams,
Nyssa Clubine Jane Devlin, Anastasia Gousseva, J-P Iamonaco, Vito Lam, Kathy Miller, Young-
Lan Shin, Rebecca Snell, Roger Philips, and Jenn Weaver.
Many written pages, figures and tables of thanks to David Pabke who provided a solution to laptop woes six days before this thesis was due.
iv
Finally, my parents Joyce Yeo and Chow Lok Leung have loved and supported me unconditionally. Armed with a B.Ed, and a Dip-Ing in Civil Engineering respectively, they listened to me patiently natter on about the significance of mud and flipped the phone bill regardless. Thank you!
v
Table of Contents
Abstract ii
Frontispiece iii
Acknowledgements iv
Table of Contents vi
List of Figures ix
List of Tables xiii
Chapter 1 – Introduction 1 1 Introduction 1 1.2 Research objective 4
Chapter 2 – Literature review 5 2 Introduction 5 2.1 Glaciations 6 2.1.1 Late Pleistocene 7 2.1.2 Late Wisconsin glaciation of northwestern Canada 8 2.1.3 Early Holocene sedimentary environment of the upper Takhini River drainage basin 11 2.2 The Holocene climate record of northwestern Canada 15 2.3 Pro-glacial fluvial hydrology and sediment transport 18 2.4 Lacustrine processes 19 2.4.1 Thermal stratification 20 2.4.2 Inflow behaviour 21 2.4.3 Turbidity current dynamics 22 2.4.4 The effect of turbidity currents upon lake bottoms 25 2.5 Varves 27 2.5.1 Varve formation 27 2.5.2 Varves as inferences of past climates 28
Chapter 3 – Study Area 30 3.1 Physiography of the Yukon Territory 30 3.2 Upper Takhini River drainage basin morphology 34 3.2.1 Primrose River sub-basin 34 3.2.2 Upper-most Takhini River sub-basin 35 vi
3.2.3 Kusawa River sub-basin 35 3.2.4 Jo-Jo Creek sub-basin 35 3.2.5 Devilhole Creek sub-basin 36 3.2.6 Kusawa Campground alluvial fan-delta complex 36 3.3 Bedrock and surficial geology 37 3.4 Glacial cover 38 3.5 Hydrology 39 3.6 Climate 42 3.7 Vegetation 45
Chapter 4 – Methods 46 4.1 Field Methods 46 4.1.1 Acoustic Profiling 46 4.1.2 Sediment cores 46 4.1.3 CTD’s 46 4.1.4 Geographic Information Systems 47 4.2 Laboratory methods 49 4.2.1 Sediment core properties 49 4.2.2 Loss on Ignition 50 4.2.3 Laser particle analysis 52 4.2.4 Geology – X-ray fluorescence 53 4.2.5 Microprobe tephra glass analysis 54 4.2.6 Radiocarbon-14 Analysis 54 4.3 Secondary data 55 4.3.1 Pb 210 and Cs 137 55
Chapter 5 – Results 56 5 Introduction 56 5.1 CTD’s 56 5.1.1 Temperature 56 5.1.2 Turbidity 58 5.1.3 Conductivity 60 5.2 Aerial imagery 63 5.3 Acoustic records and lake bathymetry 68 5.4 The sedimentology of Kusawa Lake 74 5.4.1 Sediment structure and grain size 74 5.4.2 Stratigraphy and grain size trends 87 5.4.3 Loss on Ignition 93 5.4.4 X-ray fluorescence 98 5.5 Chronology 99 5.5.1 Lead 210 and Caesium 137 99 vii
5.5.2 Radiocarbon-14 100 5.5.3 White River Ash 101 5.5.4 Rate of sediment deposition 103
Chapter 6 – Discussion 105 6.1 The Kusawa lacustrine system 105 6.2 Post-glacial sedimentary environment of Kusawa Lake 106 6.3 Sediment as a proxy for Holocene environmental change in Kusawa Lake 109 6.3.1 Late-glacial, early Holocene (10.5 – 7.0 ka cal. BP) 110 6.3.1.1 Sedimentary environment 110 6.3.1.2 Kusawa Lake diatoms 112 6.3.1.3 Climate 113 6.3.2 The Neo-glacial (7.0 – 2.0 ka cal. BP) 114 6.3.2.1 Sedimentary environment 114 6.3.2.2 Climate 115 6.3.3 Little Ice Age (1200 – 1900 AD) 117 6.3.3.1 Sedimentary environment 117 6.3.3.2 Climate 118 6.3.4 Post-LIA to present (1900 – present) 120 6.4.1 Specific Sediment Yield of Kusawa Lake 121 6.4.2 Sediment trapping 124
Chapter 7 – Conclusion 127 7.1 Spatial conclusions 127 7.2 Temporal conclusions 128 7.3 Future directions 130
List of Citations 131
Appendix A – Mean annual discharge of the lower Takhini River 141
Appendix B – CTD profiles 142
Appendix C – Laser particle size results 144
Appendix D – Loss on Ignition results 187
Appendix E – X-ray Fluorescence results 191
Appendix F – Chronology: Pb 210 and Cs 137 results, Microprobe tephra glass results 192
viii
List of Figures:
Frontispiece Kusawa Lake, July 2007 ii Figure 2.1 Conceptual model of the pro-glacial system generating lacustrine sediment 5 Figure 2.2 Glacial ice limits and ice flow directions of the Cordilleran Ice Sheet of southern Yukon Terriotry 10
Figure 2.3 Southewestern Yukon showing names referred to in the text 14
Figure 2.4 Schematic subdivision of a turbidity current 24
Figure 3.1 Southwestern Yukon showing names referred to in the text 32
Figure 3.2 The upper Takhini River drainage basin of Kusawa Lake and sub-basins 33
Figure 3.3 Mean annual discharge of the lower Takhini River 1950 – 2007 40
Figure 3.4 Annual hydrograph for the lower Takhini River at the outlet of Kusawa Lake 1986 41
Figure 3.5 Annual hydrograph for the lower Takhini River at Highway Bridge 1986 42
Figure 3.6 Temperature precipitation graph of Whitehorse based on 1971 – 2000 averages 44
Figure 3.7 Monthly rainfall normals for the Takhini River Branch based on 1971 – 2000 averages 44
Figure 3.8 Monthly snowfall normals for the Takhini River Ranch based on 1971 – 2000 averages 45
Figure 4.1 Kusawa Lake vibra core sample locations 48
Figure 5.1a Temperature profiles of Kusawa Lake taken on July 19 th , 20 th and 21 st , 2004 (southern region) 58
Figure 5.1b Temperature profiles of Kusawa Lake taken on July 19 th , 20 th and 21 st , 2004 (northern region) 58
Figure 5.2a Turbidity profiles of Kusawa Lake taken on July 19 th , 20 th and 21 st , 2004 (southern region) 60
Figure 5.2b Turbidity profiles of Kusawa Lake taken on July 19 th , 20 th and 21 st , 2004 (northern region) 60
ix
Figure 5.3 Conductivity profiles of Kusawa Lake taken on July 19 th , 20 th and 21 st , 2004 62
Figure 5.4 Aerial photo A27149-12 of Hendon River and surrounding glacial features 64
Figure 5.5a Aerial photo of A27149-100 Upper most Takhini River sub-basin mouth 65
Figure 5.5b Aerial photo A27149-63 Takhini Lake 65
Figure 5.6 Aerial photos A27217-136 and A27217-134 Primrose River sub-basin delta 66
Figure 5.7 Aerial photo A27327-52 Campground Alluvial fan 67
Figure 5.8 Locations of Regions, Lake bathymetry and CHIRP Acoustic transects 69
Figure 5.9 Acoustic section from southern Kusawa Lake near the upper-most Takhini River mouth 70
Figure 5.10 Acoustic section from Region II, south of Primrose River delta 72
Figure 5.11 Acoustic section from Region I, proximal to the Primrose River delta 73
Figure 5.12 Acoustic section from Region I, northern Kusawa Lake 73
Figure 5.13 Core stratigraphy characteristics 75
Figure 5.13a KUS324 76
Figure 5.13b KUS322 76
Figure 5.13c KUS325 76
Figure 5.13d KUS326 77
Figure 5.13e KUS319 77
Figure 5.13f KUS327 78
Figure 5.13g KUS328 78
Figure 5.13h KUS318 78
Figure 5.13i KUS330 78
Figure 5.13j KUS331 79
Figure 5.13k KUS329 79
Figure 5.14 Major sand units with laminations core KUS324 (2.65 – 2.74m) 81
x
Figure 5.15 Silt – clay laminations core KUS322 (0.84 – 0.93 m) 82
Figure 5.16 KUS319 (2.1 – 2.2 m) 83
Figure 5.17 KUS326 (2.24 – 2.31 m) 83
Figure 5.18 KUS327 (2.14 – 2.22 m) 83
Figure 5.19 KUS318 (1.04 – 1.15 m) 85
Figure 5.20 KUS331 (1.03 – 1.14 m) 85
Figure 5.21 KUS330 (0.24 – 0.44 m) 86
Figure 5.22a Percentage grain size of cores taken in Regions II and IV 89
Figure 5.22b Percentage of grain size of cores taken in Region I 89
Figure 5.23a Cross reference of major sand units in cores (0 – 1 m) 90
Figure 5.23b Cross reference of major sand units in cores (1 – 2 m) 91
Figure 5.23c Cross reference of major sand units in cores (2 – 3 m) 92
Figure 5.24 Mean organic matter relative to increasing distance from the upper Takhini River outlet 93
Figure 5.25a Percentage organic matter KUS324 &322 96
Figure 5.25b Percentage organic matter KUS325, 326 & 319 96
Figure 5.25c Percentage organic matter KUS327, 328 & 318 97
Figure 5.25d Percentage organic matter KUS330, 331 & 329 97
Figure 5.26 Major element geochemical composition of Kusawa Lake 98
Figure 5.27 Kusawa Lake down-core profiles of lead-210 and caesium-137 100
Figure 5.28 Ilmentie geochemical plots of KUS324-20 and KUS322-8 in relation to White River Ash – eastern lobe (WRA-E) and the northern lobe (WRA-N) 102
Figure 5.29 Age – depth curve of core KUS324 using radio-carbon -14 and White River eastern-lobe dates 104
Figure 6.1 Stratigraphic log of core KUS324, arbitrary climate record and Jelly Bean Lake 18 O Aleutian Low 116
xi
Figure 6.2 Specific sediment yield comparison for glacier-fed lakes in the Canadian Cordillera, plotted by percent glacier cover 123
Figure 6.3 Specific sediment yield comparison for glacier-fed lakes in the Canadian Cordillera, plotted by drainage basin area 123 Figure B 1 a & b Close-up of top (15 m) of figures 5.2a and 5.2b 142
Figure B.2 Close-up (top 30 m) of figure 5.3 143
Figure F.3 Major Element geochemistry on Kusawa Lake and Duke River White River Ash layers referring to figure 5.28 194
xii
List of Tables:
Table 3.1 Glacial cover distribution within the Kusawa Drainage Basin 38
Table 3.2 Site characteristics of Lower Takhini River discharge stations. 40 Mean annual discharge (m 3/s) figures based on 1955 – 2005 averages
Table 3.3 Climatic characteristics for Whitehorse based on 1971 – 2001 averages Lat / Long: 60° 43' N, 135° 4' W. Elevation: 706.2 m. Climate ID: 2101300. WMO ID: 71964 43
Table 4.1 Site characteristics and core sample lengths taken in July 5 – 8th 2006 49
Table 4.2 Cores and site characteristics of Pb 210 and Cs 137 analyses 55
Table 5.1 Radiocarbon dates and calibrated calendar ages from core KUS324 101
Table A.1 Mean annual discharge of the lower Takhini River 1950 – 2007 referring to Figure 3.3 141
Table C.1 Percentage grain size of cores taken, referring to figures 5.22 a – d 144
Table C.2 Laser particle size results of cores, referring to figures 5.22 a-d 147
Table D.1 Loss on Ignition results referring to figures 5.24 and 5.25 a – d 187
Table E.1 Major element geochemical composition of Kusawa Lake, H 20 calibrated, referring to figure 5.26 191
Table F.1 Kusawa Lake down-core profiles of lead-210 and caesium-137 , referring to figure 5.27 192
Table F.2 Average percentage major element composition of glass shards from two Kusawa Lake cores 193
Table F.3 Average percentage magnetite and ilmenite composition of glass shards from two Kusawa Lake cores referring to figure 5.26 194
xiii
I – Introduction
1.1 Introduction
Within the last 50 years, there has been a rapidly increasing concern regarding global climate change and its impact on water resources. Glaciers in Canada and Alaska cover ~ 90,000 km 2 (13 %) of mountain glacier area on the earth’s surface. More than 50 % of the Canadian landmass is underlain by permafrost. Of that 50 %, approximately 60 % is at temperatures just a few degrees below the melting point (Bockheim, 2006). Environmental scientists have noted that continual increases in mean annual air temperatures in the Arctic latitudes have induced dramatic changes (Arendt et al., 2002; Moore, 2002a). Warmer temperatures can gradually lead to thawing and destabilization of perennially frozen grounds (Camil, 2005; Hinzman et al.,
2006b). Since the 20 th C., glaciers in the Canadian Cordillera have lost 25% of their area
(Luckman, 2000).
Glacial melt in northwestern Canada and Alaska can exert critical controls upon the landscape, most notably, hillslope processes and hydrology including surface and sub-surface ground water regimes and fluvial – lacustrine processes. Glacier and ground ice melt have contributed to a rise in rates of watershed discharge. With an increase in surface glacier melt, a decrease in snowpack accumulation there is a noticeable increase in, surface water discharge increases and, as a result, an elevated hydrologic base-flow in the summer. As glacial ice degrades to the point of becoming discontinuous, surface water infiltration through soils to groundwater can eventually enter the stream drainage network and provide base-flow during both the dry and winter periods (Burn, 2002; Camil, 2005; Hinzman et al. 2006a). Consequently, coastal and lacustrine environments have experienced an increase in the magnitude and frequency of higher-level discharge (Arendt et al. 2002; Lowey, 2002; Gilbert et al. 2006).
1
In order to put contemporary climatic warming and the current state of glaciers in north- west North America into context, past climatic conditions prior to the onset of anthropogenic influences must be understood. Unfortunately, Canada only holds approximately 70 years worth of detailed metrological data and the availability of climatic data for remote locations such as the glacial mountains of northwestern Canada is limited. Consequently, and especially within the last three decades, there has been a surge in the study of paleoenvironmental reconstruction for northwestern North America (Serink 2004). Current research challenges the fundamentals of unveiling historical climatic variation using proxies such as O 18 – 16 isotopes from ice cores, tree rings, marine sediment and clastic sediments from glacier-fed lakes, and diatoms (Last & Smol.
2001; Menounos et al., 2004; Anderson, 2005).
Paleolimnological records from large glacier-fed lakes have contributed significantly to our understanding of Pleistocene and Holocene environments. Studies regarding glacial lake environments within western Canada by Gilbert (1975), Desloges & Gilbert (1994a, 1994b
1998), Lamoureux et al. (1996), Dirszowsky & Desloges (1997, 2004), Menounos et al. (2005),
Menounos (2006), Gilbert et al. (2006) and Hodder et al. (2006) have all shown that lacustrine deposits offer a potentially continuous archive of sediment production from a contributing drainage basin. This in turn has the potential to help infer past climatic variability. Clastic sedimentary deposits can reveal information regarding sediment source characteristics, sediment transport and depositional mechanisms, extreme depositional events and lake conditions. This form of study is not limited to Canada. Similar approaches in reconstructing environmental change in glaciated mountain regions have also been pursued in Sweden (Anderson et al., 1996), in the Swiss Alps (Ohlendorf et al.1997), in Argentina (Strelin & Malagino 2000), and in Central
Finland (Ojala & Saarnisto 1999).
2
Relatively little research has examined the potential to reconstruct paleoenvironments from clastic sediment records in large glaciated lakes within the southwestern Yukon Territory.
Kusawa Lake, located within the Dezadeash Region of the Boundary Ranges, was selected for this purpose. It is 70 km in length and has a total area of 142 km 2. Kusawa Lake is part of the
Takhini River basin, in which the upper part (feeding Kusawa Lake) has a total area of 4292 km 2. Glaciers cover 4.7% of the drainage basin. Glacial history of the Dezadeash Region has played a considerable role in determining the morphology of Kusawa Lake. Once occupied by ice-dammed Glacial Lake Champagne during the Late Pleistocene deglaciation, the onset of climatic warming during the early Holocene caused ice-dam breakage, regional flooding, and the deposition of thick and fast sediment (Eyles, 1990; Gilbert & Desloges, 2005).
Significant shifts in climate were experienced within the region throughout the Holocene.
Research carried out within the northwest of Canada has shown this variability. Paleo proxies from neighbouring environments to Kusawa Lake such as Pyramid Lake and Jellybean Lake
(Anderson et al., 2005; Mazzaucchi et al., 2003) of the Y.T. have shown that the Aleutian Low played a profound role in controlling northwest Pacific Holocene climates. Flood events generated from incursions of warmer and moist Pacific air masses have resulted in numerous synchronous glacial activities in the St. Elias and Coast Mountain Ranges, such as the Little Ice
Age advance, which brought about instability of both sub-aerial and subaqueous terraces in
Kusawa Lake. Given the continual shift in water levels throughout the Holocene, these changes have had a significant effect upon the rate of sediment delivery into Kusawa Lake.
More recently, catastrophic flood experienced at the Kusawa Lake Campground alluvial fan show how rapid climatic changes acts as a factor in determining the morphology of present day glaciolacustrine environments. Throughout the summer of 1982, the persistence of above average precipitation within the Kusawa Lake region caused permafrost thawing and hillslope 3 failure in the basin above the alluvial fan. Approximately 6.7 x 105 m 3 of unconsolidated sediment was transported downslope and deposited on the alluvial fan-delta (Lowey 2002).
Private property and recreational facilities within the area were flooded and destroyed. Along with previous studies using proxies such as tree-ring evidence, Lowey (2002), Gilbert &
Desloges (2005), and Lipovsky (2006) have shown that relatively smaller flood events have occurred more frequently within the last 150 years. Although the torrent system is presently less active, there is still the potential for future flood events and landslide dam failures.
Using the lacustrine sediment accumulation record from Kusawa Lake as a proxy evidence for understanding environmental change may help in understanding shifting climatic patterns throughout Holocene.
1.2 Research Objective
The central focus of this research is to determine how the sedimentary record of glacially- fed Kusawa Lake has varied throughout the Late Glacial and Holocene intervals and subsequently whether these variations are tied to known regional climatic fluctuations.
The more specific research objectives are:
a) Examine how the Cordilleran Ice Sheet and de-glacial events have determined the shape
of the lake basin and, subsequently, the inherent sediment properties of Kusawa Lake
deposits.
b) To characterise the spatial variability of clastic depositional sequences and relate these to
possible sediment sources and delivery processes.
c) To examine Holocene changes in the rate and character of deposition, and relate this to
possible local and regional hydroclimatic forcing.
4
II – Literature Review
2. Introduction
A number of factors influence the relation between glacier activity and lacustrine sedimentation. In some cases, these factors are linked in unidirectional or bidirectional ways
(Fig. 2.1). Fig. 2.1 is a conceptual model, which summarizes the process – network of a typical pro-glacial system generating lacustrine sediment (Hodder et al., 2007). Six inter-related systems are specified leading to the ultimate product of lacustrine sediment deposits.
Fig. 2.1 Conceptual model of the pro-glacial system generating lacustrine sediment.
Dotted lines = connections between lacustrine sediment and the fluvial, climate and / or glacial system
(Hodder et al. 2007) 5
Given the complexities associated with the range of processes that control lacustrine sediment yield, it is important to simplify individual components that act as contributors (Fig.
2.1). Due to the nature of this project, this chapter will focus on the processes and interactions between climate, glacial, fluvial and lacustrine systems in northwestern Canada.
This chapter will begin by briefly examining the Late Pleistocene glacial episode and its impacts, particularly relating to the advance and retreat of the Cordilleran Ice Sheet in northern
British Columbia (B.C.), and southwestern Yukon Territory (Y.T.) during the Late Wisconsinan.
This is followed by a review of the climatic and geomorphic conditions throughout the Holocene up until the present day. Because of heavy modification of the regional topography by the
Cordilleran Ice Sheet, the relevant climatic and geomorphic elements controlling sediment yield will be reviewed to provide context for subsequent inferences from, and discussion of, lacustrine sediments.
2.1 Glaciations
During the Pleistocene epoch (1.65 M.Y – 10 ka cal. BP), a series of warming and cooling events have led to extensive ice sheets that covered up to 35 % of the earth’s high altitude and high latitude regions (Easterbrook, 1999). Causes of these climatic shifts include plate tectonics and intermittent volcanic activity, the Milankovitch insolation cycles and carbon and methane cycling (MacDonald et al., 2000; Hodder et al., 2007). These climatic shifts triggered glacial events, which have had a significant impact upon the landscape in North
America. At a regional scale, glaciations during the Late Pleistocene played a large role in the erosional and depositional landscapes in the southwestern Yukon Territory. Of particular importance to the lacustrine system was the formation of Glacial Lake Champagne during the last deglaciation. 6
2.1.1 Late Pleistocene
The Late Pleistocene (130 – 10 ka BP) began at the Eemian interglacial phase before the final glacial episode of the Pleistocene known as the Holocene (Easterbrook, 1999). The Late
Pleistocene demonstrated distinctive climatic characteristics that included repeated advances and retreats of the continental ice sheets (Flint, 1957; Clague et al. 1987). There have been indications that more than twelve glacial cycles have occurred within northwest Canada alone.
The most notable glacial episode was the Wisconsinan Glaciation, which dominated North
America (Flint, 1957; Fulton, 1991; Easterbrook, 1999). The last glacial stage can be subdivided into the early (~ 75 – 64 ka), middle (64 – 26 ka) and late (26 – 10 ka) Wisconsinan sub-stages
(Easterbrook, 1999; Trenhaile, 2004). During the Late Wisconsinan, the continental ice sheet had three main components. These included the Laurentide ice sheet in central and eastern North
America, the glacial complex in the High Arctic, and the Cordilleran glacial complex (Trenhaile,
2004). Stratigraphic and sedimentological studies supported by C 14 dating show that the
Cordilleran Ice Sheet formed sometime between 29 – 25 ka cal. BP (Ryder & Maynard, 1991;
Easterbrook, 1999). The Cordilleran Ice sheet spanned north – south from the Northwest
Aleutian Islands, to Mt. Adams near the Columbia River in Washington State, and was approximately 900 km wide from the Pacific to the foothills of the Rockies (Fulton, 1991). C 14 dating of organic lacustrine deposits and the characteristic presence of loess deposits has shown that the maximum glacier extent to the south occurred at ~ 18 ka cal. BP and had extended to
Montana, and Washington States (Hobbs, 1947). The extreme northwest Yukon remained ice free.
7
2.1.2 Late Wisconsin glaciation of northwestern Canada
Bostock (1966) identified four advances of the Cordilleran Ice Sheet, within the southern
– central Yukon: the Nansen (oldest), Klaza, Reid, and McConnell (youngest and the latest) glaciations. During the preceding interglacials, upland areas were ice-free. Since the more recent ice advances were the most extensive and most erosive, landforms found throughout southern – central Yukon are predominantly from the Reid and McConnell glaciations (Kerr,
1934; Hughes 1987).
Based on stratigraphic and geomorphic evidence, the Late Wisconsin Reid – McConnell glacial advance began circa 25 – 29 ka cal. BP and covered a vast region from the Pelly and
Selwyn Mountains in the east, to south central Yukon, and from the Cassiar Lobe and to the eastern Coast Mountains of B.C. (Jackson et al. 1991) (Fig. 2.2). Glacial ice from the Cassiar and Coastal Mountains converged and flowed 1600 km in a north and northwesterly through the
Yukon and Teslin River Valleys and ceased north of the Yukon Plateau (Gilbert & Desloges
2005). Westward ice flow from the Pelly and Selwyn Mountains was contiguous with an eastward piedmont glacier complex from the St. Elias Mountains. Ice from the Lowell and
Dusty Glacier lobes, which constituted the main mass of the St Elias complex, flowed east, northeasterly, and north into the valleys of the Slims, Kaskawulch, Dusty, and Alsek rivers (Fig.
2.2) (Kindle, 1953; Clague, 1989; Braher et al., 2008). Glacial ice reached the Shakwak Valley in the Dezadeash Region through the Denali fault line. At Much Lake, ice flow continued in a northerly direction along the Duke Depression through Alder Creek and eastwards towards
Dezadeash Lake where it accumulated (Kindle, 1953; Jackson et al, 1991; Gilbert & Desloges,
2005). This ice subsequently merged with the densely compact glacial ice flowing from the
8
Klukshu and Takhanne River Valleys. The enlarged ice field flowed northeast along the
Dezadeash Valley and Shakwak Trench into the Kusawa Lake area.
Prominent glacial positions of the glacial maximum can be traced through moraines and erratic boulders. Evidence of glacial scouring can be found at elevations up to 2000 m a.s.l.
(Kindle, 1953; Gilbert & Desloges, 2005) suggesting this was the maximum ice thickness.
Kindle (1953) noted Shakwak till exposures 120 km northeast of the present divide in the Ice
Field Ranges between the Kaskawulsh and Hubbard glaciers. Isolated nunataks, which protruded through the ice, suggest ice elevations in the southern Dezadeash, Ruby, and Kluane
Ranges reached 1700 m a.s.l. Based upon the floor elevations of the lowest cirques in the Ruby
Range of the Yukon Territory, which once supported cirques – glaciers, paleo-firn lines fell as low as 1500 m a.s.l. (Jackson et al., 1991). Kame terraces found at elevations of 1050 – 1350 m a.s.l. near Jo-Jo Lake and Dezadeash River indicate the long-term presence of the ice. The movement of ice, in the Shakwak Trench and Takhini River valley caused bedrock scouring that resulted in the deepening of the valley system. Glacial Lake Champagne was formed towards the end of the McConnell glaciation, when valley glaciers still occupied the large valleys of the
Kluane Ranges west of the Shakwak Valley, Takhanne, and Klukhu river Valleys south of
Dezadeash Lake (Kerr, 1934; Kindle; 1953). The Kusawa Valley remained glaciated to the early
Holocene (Jackson et al., 1991; Gilbert & Desloges, 2005).
9
Fig. 2.2 Glacial ice limits and ice flow directions of the Cordilleran Ice Sheet of southern Yukon Territory.
EL = Eastern Coast Ranges Lobe; EC = Eastern Coast Ranges; CL = Cassiar Lobe; C = Cassiar Mountains; LL =
Liard Lobe; P = Pelly Mountains; S = Selwyn Mountains, SL = Selwyn Lobe. Arrows = flow directions. Red dash line = Ice sheet divide. NTS Digital Elevation Model (Hillshade) at 500 m intervals.
Ages documenting the retreat of the Cordilleran Ice Sheet in the Yukon have been limited. The C 14 dating of organic – rich silt in lacustrine sediments derived from the terminus of the St. Elias piedmont lobe complex has provided an estimated age of 13,660 +/- 180 yr BP
(Rampton, 1971 in Jackson et al. 1991). The end of McConnell de-glaciation marked the start of the Holocene, which encompasses approximately the past 10.5 ka BP. A major feature of the
10
McConnell de-glaciation were the large ice-dammed lakes that occupied the interior valleys of the Yukon River drainage basin. Thick sequences of glaciolacustrine sediments were deposited and these have shaped the valley morphology we see today (Gilbert & Desloges 2005).
2.1.3 Early Holocene sedimentary environment of the upper Takhini River drainage basin
Due to the large spatial extent of the strandlines (indicators of water levels) along the
Dezadeash Valley from the Takhini River to the south-east of Haines Junction, and in the valleys of the Dezadeash, Kathleen and Frederick Lakes, it is possible that Glacial Lake Champagne was dammed by glaciers from the west in close proximity to the St. Elias ice complex (Gilbert &
Desloges 2005) during the early Holocene. Glaciers in the Primrose Valley and glaciers near the
Kusawa Lake outlet (Fig. 2.3) formed large deltas in Glacial Lake Champagne (Gilbert &
Desloges, 2005). Because a trunk glacier occupied the southern portion of Kusawa Lake, the majority of sediment accumulation was derived in response to the early Holocene deglaciation and the complete drainage of Glacial Lake Champagne.
An assessment of strandlines by Gilbert & Desloges (2005) showed that water levels once stood at ~ 772 m a.s.l. in the northern (outlet) part of Kusawa Lake during the early
Holocene (contemporary lake level is 671 m a.s.l). This was coincident with a persistent high – level of Glacial Lake Champagne (Gilbert & Desloges, 2005). Less prominent strandlines were found at 853 m a.s.l, which suggests that an early phase of Glacial Lake Champagne once stood a higher level. Subsequently, the lake grew and levels fell. Glacial-lake elevation data suggests that the most prominent strandlines along valley walls occur at 765 and 725 m a.s.l. with lesser strandlines occurring between 725 and 765 m a.s.l and down to 701 m a.s.l. This indicates that water levels had fluctuated significantly throughout the deglaciation process.
11
Throughout this drainage process, Lake levels were controlled by both a spillway floored at 756 m a.s.l at the north of the Nordenskiold River along with a continual down-cutting of the sediment plug at the outlet of Kusawa Lake. This indicates differential isostatic rebound of 0.2 m/km from south to north (Kindle, 1953; Gilbert & Desloges, 2005). Uniformly thick glaciolacustrine fine sand, silt and clay sediment were deposited on the Takhini and Dezadeash valley floors (Fig. 2.3) (Lowey, 2002; Gilbert & Desloges, 2005). Using acoustic records,
Gilbert & Desloges (2005) noted three distinct sediment facies within Kusawa Lake. They propose that the first facies was deposited when the glacier was positioned in the middle of
Kusawa Lake, and another glacier had occupied the Primrose Valley. The second sediment facies was deposited over facies 1 as the lake level continued to decrease and was found to be discontinuous and thickest near Primrose River delta. However, this facies is not continuous.
The third facie is finely layered and found to lie conformably over facie 2 and extends to the modern sediment surface.
However, Gilbert & Desloges (2005) noticed that reflectors are stronger near the surface and the sediment is more focused towards the deeper parts of the lake. Given the acoustic records, they estimate that the total glaciolacustrine sediment thickness is greater than 100 m, with the majority of it related to immediate de-glaciation and glacial lake draining. This thick and fast process of sedimentation conforms to the thick and acoustically stratified sediment found in other large lakes in B.C. such as Okanagan and Lillooet Lakes (Eyles, 1990; Desloges
& Gilbert, 1994)
Repeated glacial advances and retreats have shaped the regional topography, and have a significant influence upon lacustrine sedimentation processes (Fig. 2.1). Unfortunately, due to the remoteness and limited accessibility of glaciated regions, archives of climate change within the Cordillera are hard to develop. Therefore, there has been an emphasis on interpretation of 12 small sites (ponds, bogs etc.) that hopefully provide a regional picture of hydrologic forcing.
Assessing late glacial and Holocene climate impacts on large watersheds have been more problematic, but such an approach has the potential to identify regional wide geomorphic effects of climate change. Kusawa Lake and the upper Takhini River drainage basin provide the opportunity to assess the utility of large lake records in reconstructing Holocene climate change.
In order to achieve this, the next section will discuss the Holocene climate record within northwestern Canada.
13
Fig. 2.3 Southwestern Yukon, showing names referred to in the text. Spot height elevations from Gilbert & Desloges (2005) Arrows
indicate large deltas built into Kusawa Lake
14
2.2 The Holocene climate record of northwestern Canada
Scientists use multi-proxies to reconstruct Pleistocene and Holocene paleoenvironmental conditions at a variety of resolution scales. Most common analyses include: isotopic δ018 derived from snow pack or sediment accumulation, dendrochronology (Alley et al., 1993;
Gedalof & Smith, 2001; Moore, G.W.K et al., 2002; Anderson et al., 2005), lithostratigraphic analyses, pollen records, and biochemical remains of past phototropic communities preserved in lake sediments (Pienitz et al., 2000; Menounos, 2006). However, due to response and lag-time of glacier movement and sudden shifts in atmospheric – oceanic processes, determining climatic conditions during the Holocene can be extremely complex. In northwest Canada, it is possible to recognize a series of minor Holocene glacial advances and retreats that characterize the impact of the Holocene climate change.
During the late Pleistocene and into the early Holocene (ca. 14 – 10.5 ka cal. BP), rapid shifts in the sea surface temperature and temperature gradients in the atmospheric circulation occurred over the northwest Pacific and Arctic Oceans. The Aleutian Low (AL); a semi- permanent low-pressure system located over the Gulf of Alaska, controlled northwest Pacific
Holocene climates (Moore, G.W.K et al., 2002; Anderson et al., 2005) and was a major factor in explaining glacial activity. Regional factors such as topography and proximity to the Pacific
Ocean led to the decay of the Cordilleran Ice Sheet (Jackson et al. 1991). By the end of the
McConnell glaciation, the Cordilleran Ice Sheet decayed rapidly via of down wasting and stagnation and exposed large areas of the continental shelf (MacDonald et al., 2000).
The Holocene (10.5 ka cal. BP – present) is associated with a series of climatic warming and cooling events that are known to affect lacustrine sedimentary records. Multi-proxy analysis of pro-glacial lakes within northwestern Canada showed periods of high temperatures and dry 15
conditions between 10.5 – 7.0 ka cal. BP (Johnson, 1992; Wiles et al., 2002; Anderson, 2005).
Higher rates of sedimentary accumulation and more variable inputs of sediment were recorded in many Cordilleran lakes (Barber & Finney, 2000; Hodder et al., 2006, Gilbert et al., 2006).
Clastic material inputs into several lakes across the interior of B.C. were also high and nutrient levels were poor (MacDonald et al., 2000). Within the Y.T, an analysis of fossil pigments, diatoms and lake stratigraphy on Lake U60 in the Pelly Mountains showed that primary productivity of chlorophytes and cyanobacteria were greatest following the retreat of the glaciers and the occurrence of fine sediment lamination was indicative of meromitic water conditions
(Pienitz et al., 2000). This coincides with a peak in Populus found in late glacial sediments across the Y.T. (Lacourse & Gajewski, 2000). All results obtained suggest that these dry conditions persisted until ~ 7.0 ka cal. BP.
The mid-Holocene (~ 7.0 – 2.0 ka cal. BP.) also known as the Neoglacial is characterized by a slightly cooler and wetter climate due to an increase in effective moisture content from the
Aleutian Low and an increased intensity of the Pacific air mass (Last, 2002; Spooner, 2002;
Mazzucchi et al., 2003). High levels of exotic western hemlock pollen were noticed in
Waterdevil Lake, northwestern B.C. Similarly, a profusion of white spruce, poplar and willow pollen were found in glacial lakes across southwest Yukon and the Mackenzie Mountains
(Lacourse & Gajewski, 2000). This suggests an increase in effective moisture, which encouraged growth and subsequently the distribution of pollen. In addition, this corresponds with sedimentation rates in Lake U60 of the Y.T., which varied by an order of magnitude reflecting an increased variability of sediment and organic matter input from the watershed due to a moister climate (Pienitz et al., 2000). Using isotopic δ 18 0 analyses of sediment cores excavated out of Jellybean Lake, Anderson et al. (2005) noted a significant increase in atmospheric circulation intensity from the southwest Pacific between 4.5 – 3.5 ka cal. BP. The 16
Aleutian Low is theorised to have been intense and is estimated to have been situated eastward over Alaska and the Y.T.. This eastward shift had a profound orographic influence upon precipitation and subsequent lake levels and increased rates of sedimentation (Moore, G.W.K et al., 2002).
During the onset of the Neoglacial warming phase (~ 3.5 – 2.0 ka cal. BP), the Aleutian
Low was signified as weak (Moore, G.W.K. et al., 2002; Anderson et al., 2005). The westward shift in the Beaufort Gyre allowed the counter-flowing Mackenzie Current to reposition, strengthen, and persist until present day (Dyke & Savelle, 2001; Wiles et al., 2002). Once again, this abrupt change in circulation patterns had an impact upon glaciolacustrine environments across northwestern Canada. The slightly warmer temperatures and distinctive varves were recorded in Lake C2 of northern Ellesmere Island between the years 3.3 – 3.1 ka cal. BP.
(Lamoureaux & Bradley, 1996). Last et al. (2002) recorded similar results in Oro Lake, B.C.
Likewise, Brahney et al., (2008) noted that distinct sequences of rhythmic sediment laminae in
Kluane Lake at around 2.8 ka cal. BP corresponds to sediment exhaustion from preceeding glacial advances in the St. Elias Mountains (Denton & Stuiver, 1966).
The end of the Neoglacial was marked by an intensification of the Aleutian Low and shifts in the North Pacific Index ~ 1.2 ka cal. B.P. (MacDonald et al., 2000; Anderson et al.,
2005) which brought about a warmer and moister climate. Subsequently, glaciers retreated
(Denton & Stuiver, 1996). This corresponds with the profusion of western hemlock ( Tsuga heterophylla ), which extended into the eastern ranges near Pyramid Lake ~1.5 ka cal. BP
(Mazzucchi et al., 2003). Likewise, a sudden high influx of Duke River sediment into Kluane
Lake ~ 1.3 ka cal. BP. also coincides with the presence of a warmer climate (Brahney et al.,
2008).
17
The ‘Little Ice Age’ (LIA) is a climatic cooling period, which is thought to have occurred between the 12 th and 13 th Centuries, and the early 18 th to mid 19 th Centuries throughout most of northern B.C and Y.T. Dendrochronology, lichenometry and radiocarbon dates show that at least three cooling periods occurred during the 18 th – 19 th C. (1650, 1700 and 1850 A.D), each separated by slightly warmer intervals (Luckman, 2000; Loso et al., 2006). The Aleutian Low moved eastward and intensified over northwest North America, which brought cooler temperatures during the spring and periods of increased precipitation throughout the summer
(Luckman, 2000; Wiles et al. 2002). The increase in effective moisture triggered episodes of synchronous glacier advances in the Coast and St. Elias Mountain Ranges ~1850 A.D. (Johnson,
1992; Lacourse & Gajewski, 2000; Anderson et al., 2005).
2.3 Pro-glacial fluvial hydrology and sediment transport
Lake sedimentation processes are closely related to the temporal and regional variation in hydrology and the topography of the drainage basin (Fig. 2.1) (Håkanson, 1983). In order to understand the relationships, this section will discuss, with relevant examples, the interaction between pro-glacial fluvial environments, which ultimately leads to the sedimentation of lakes.
Major sources of sediment into a fluvial and lacustrine environment of northwest Canada typicaly include glacier-derived sediment, rockfall / slide material and peri-glaical slumped from lateral valley sides. Studies conducted at Lillooet Lake, B.C. and the surrounding regions indicate that glacier derived debris alone can contribute up to 75 – 83% of the fine sediment load
(Desloges & Gilbert, 1994). However, not all sediment is deposited directly into the lake. In some cases, sediment traps within the upper reaches of the drainage basin (Smith, 1981), and alluvial fans (Parker et al., 1998; Lipowsky, 2006) act as storage sites. The trapping of sediment introduces a lag effect into the sediment cascade (Brierley et al., 2006; Hodder et al., 2007). In a 18
study of Small River Glacier basin in the Canadian Rockies Orwin & Smart (2004) discovered that at least 80% of the downstream sediment load is derived from the storage in the pro-glacial area and not from contemporary glacier bed erosion.
The mobilization of sediment from stored sediments within a pro-glacial environment is controlled by the hydrology during a given year (Marren, 2005). In addition, the pro-glacial area can function as both a source and a sink of sediment during a melt season or period within a melt season (Brierley et al., 2006). However, a transport-limited sediment system occurs when sediment size and sediment supply is greater than the energy available to move it. Within a sediment-rich pro-glacial system, each flood event of a given magnitude has approximately the same ability of mobilizing and thus transporting the same volume of sediment. If the system is supply-limited, then the sediment supply may become exhausted over time and hydrologic changes will not be recorded in the downstream and lake sedimentary records. Supply limitations may occur when the channel sediment recharge rate becomes lower than the rate of sediment mobilization in the channels (Knighton, 1998; Marren, 2005). Pleistocene glaciations have ensured a mostly transport-limited situation during the Holocene for many of the glacially covered terrains of the western Cordillera (Church & Ryder, 1972). A transport-limited situation offers better understanding of the magnitude-frequency relation in comparison to a supply- limited system. Temporal lags in the availability of sediment may mean that floods of similar magnitudes yield different amounts and sizes of sediment due to the exhaustion effect (Leopold et al., 1964; Knighton 1998; Trenhaile, 2004). A number of studies have shown the importance of pro-glacial sediment sources during the onset of a high magnitude-flooding event (Eyles,
1990; Warburton, 1990; Desloges & Gilbert, 1994).
19
2.4 Lacustrine processes
Sedimentation patterns are determined by seasonal and diurnal changes within the lacustrine environment. This section will discuss the theory behind water circulation, its control on the distribution of sediments, and in particular, the formation of turbidity currents (see box in
Fig. 2.1). This section will conclude by looking at the effects of turbidity currents upon lake bottom sediments.
2.4.1 Thermal stratification
Glacier-fed lakes of northwestern Canada are commonly diamictic (Crookshanks, 2008).
Water and sediment circulation patterns are dependent on the thermal structure of the water column and the Coriolis Effect (Hodder et al., 2007). This ultimately influences the characteristics and distribution of sediment throughout the lake (Drewry, 1986). The physical structure of lakes is governed by seasonal changes in temperature. In general, there is often a difference between the density of lake water near the bottom and the top of the water column.
The thermocline separates the warmer water above from the colder and denser water below
(density of water is greatest at ~ 4°C) during a warm season (Hodder et al., 2007). The Coriolis
Effect also has an important control upon internal waves. This occurs only when the Coriolis
Effect is balanced by the force of gravity. As a result, a standing wave rotates around the basin along the thermocline and contributes to significant mixing, which can then influence the distribution of energy and particulates within the water column (Smith & Ashley, 1985; Drewry,
1986; Hodder et al. 2007).
In a pro-glacial environment, direct heating and runoff from surrounding tributaries in the spring and early summer encourages lake water temperatures to rise and gradually become denser, which ensures some mixing. Pro-glacial lakes are commonly stratified by the early 20
summer. During the autumn, the fall in lake temperatures encourages lake turnover. There is also the potential for the entire lake to become isothermal at several intervals of the year. This occurs when the lake surface water cools to approximately 4 oC and the density of the lake water becomes highest at the surface allowing the upwelling of less dense colder water from the bottom. It is unclear if all large glacier-fed, high latitude, lakes are diamictic, but at least one major turnover event probably would occur (Drewry, 1986; Knighton, 1998; Hodder et al.,
2007).
2.4.2 Inflow behaviour
When a river enters a pro-glacial lake, three spatial patterns of inflow may be observed.
These are overflows (hypocycnal), interflows (homopycnal) or underflows (hyperpycnal). The difference between these types of flows is dependent on the hydrostatic pressure of the lake, the density difference between the lake and river water and changes in temperature gradients
(Middleton & Hampton, 1973; Leeder, 1982; Lowe, 1982; Edwards, 1992; Knighton, 1998).
Overflow occurs when inflow density is less than the density of the lake epilimnion.
Wind stress also influences the surface current flow velocity and in some cases, the Coriolis force deflects the dispersal of sediment towards the right (Smith & Ashley, 1985; Drewry, 1986).
Overflows are normally observed during times of low sediment input when thermal properties control the density difference between the two water masses (Smith & Ashley, 1985). In glacial lakes with low rates of fine-sediment input, the process of interflow may occur where the density of inflowing river water is equal to the density of unstratified lake water. In the case of thermally stratified pro-glacial lakes with high sediment input, interflow at the thermocline may alternate with underflow on a seasonal basis. The concentration of suspended silts and clays in overflows and interflows are typically low (5 – 30 gm/l). Consequently, sediment deposition via 21
settling from the water column, accumulations on the lake bed are thin, regular, have sharp contacts and display no current structures (Gilbert, 1975; Smith & Ashley, 1985; Middleton,
1993).
Underflows are widely recognised as the principal mechanism of sediment distribution in pro-glacial lakes. There are two types of underflows which are typical of glacier-fed lakes and these are disintegrated by differences in density and velocity. The first type is quasi-continuous flows, where bottom currents form when high-density river water “underflows” the hypolimnion of the lake (Lowe, 1982; Smith & Ashley, 1983). Sediment deposits are a mixture of fine sands and silts (i.e.: little upward fining) because sediment delivery persists over a long period of time with intermittent pulses. The second type of underflow is from a more catastrophic source such as slumps from mass movement processes that generate surge-type currents (Middleton &
Hampton, 1976; Edwards, 1992). Sediment deposits display upward thinning and fining due to the short-lived nature of sediment input. They become more ‘local’ in their effects. Since high inflow densities are often caused by high-suspended sediment concentrations, both types of underflows are a type of turbidity current (Smith & Ashley, 1983; Middleton, 1993; Serink,
2004).
2.4.3 Turbidity current dynamics
Turbidity currents are particle-laden underflows, which travel along the lake or ocean bed. Because turbidity currents involve the transport of water and sediment, they can be classified as both fluid and sediment gravity flows. In a fluid context, turbidity currents are classified as stratified gravity currents (or density currents). Flow takes place below, above or between ambient fluids because of a difference in hydrostatic pressure-force between the two.
The difference in hydrostatic pressure-force may be due to differences in sediment composition, 22
temperature or salinity (Allen, 1970; Middleton, 1993; Kneller & Buckee, 2000). In a sedimentary context, turbidity currents are classified as a type of sediment gravity flow where grain suspension is mainly supported by the upward component of fluid turbulence due to water circulation in the sediment water mixture. Sediment gravity flows can be distinguished from other types of flow i.e.: debris flows where sediment is dispersed by other mechanisms
(Middleton, 1993).
The vertical concentration structure of the turbidity current is a critical factor in understanding flow and sediment dynamics (Allen, 1971). Depending on the height at which maximum velocity is achieved, the turbidity current can be divided into the inner and outer regions. Both regions have separate flow processes: within the inner region, turbulence is created by bed friction, whereas in the outer region turbulence is created by friction due to the entrainment of ambient fluid (Fig. 2.4). The height of the velocity maximum is a function of the ratio of the drag forces at the upper and lower boundaries (Middleton, 1976; Allen, 1971) and in many experimental studies, this height is typically 0.2 – 0.3 ratio of the current depth (Edwards,
1992; Middleton, 1993; Kneller & Buckee, 2000).
23
Fig. 2.4 Schematic subdivision of a turbidity current flow pattern and structure of the turbidity current head where y is the height from the lake bed and u is the velocity (modified from: Middleton, 1993; Kneller & Buckee, 2000)
In terms of a latitudinal concentration structure of a turbidity current, turbidity currents are found to have three distinct parts: the head, the body and the tail. The head of a density surge is lobate in form and possesses a characteristic ~ 10% overhang (Fig. 2.4) (Allen, 1970;
Middleton, 1993). Unlike the body and the tail, the head is denser and has a higher concentration of coarse material (Bouma, 1964; Middleton & Hampton, 1976). In order for the head to advance, it must have more gravitational potential energy to displace the ambient fluid which produces resistance to the flow (Middleton, 1993). Large eddies within the head are usually found to have velocities that exceed 50% of the maximum mean downstream velocity (Kneller &
Buckee, 2000). Due to the friction experienced at the base of the current, the head is a region in which erosion takes place (Middleton, 1993).
Flow within the body region of the current is found to be uniform in thickness and less dense in terms of sediment concentration. In order to achieve a constant rate of flow, the average velocity in the body must be greater than that of the head. However, this is dependent on the rate 24
of loss of fluid from the head and slope of the bed (Leeder, 1982; Edwards, 1992; Middleton,
1993). In a study carried out in a flume setting, Kneller & Buckee (2000) note that high magnitude inflow velocity maybe up to 40% higher than the maximum mean velocity in the body. This maximum inflow velocity may therefore be equivalent to the velocities in the current head. With increased Reynolds stress (Re>2500) there is the potential for sediment entrainment within the body and subsequently bed erosion.
A gradual thinning of flow depth and sediment concentration is found at the tail of the current (Bouma, 1962; Allen, 1970). The mixing of the current and the water above produces an entrained layer.
Sediment sorting in turbidity currents is highly dependent on clastic properties such as size and weight (Edwards, 1992). The sediment particle entrained within a turbidity current must overcome the gravity component and have the energy to overcome the flow of the fluid drag
(friction) exerted at the bed boundary of the lake. Turbidity currents can travel over long periods
(days) and over a great distance (10's of kilometres) until energy (velocity) is lost (Edwards,
1992; Hånkanson, 1983). The deposition of sediment commences with the coarsest material, typically from the head of the turbidity current, with finer particulates deposited above the initial deposit or further up lake, as the body and tail of the turbidity current passes (Lowe, 1982;
Edwards, 1992; Middleton, 1993).
2.4.4 The effect of turbidity currents upon lake bottoms
In general, the concentration of suspended sediment in the water column of glacier-fed lakes is, on average, much greater than in freshwater or marine environments. The inflow of sediment is largely inorganic with low concentrations of dissolved sediment loads (Gilbert, 1975;
Smith & Ashley, 1983). As discussed in section 2.4.3, any abrupt changes in lake temperature 25
variations and / or inflowing water have the ability to create turbidity currents and induce sedimentation (Serink, 2004; Edwards, 1992).
The effects of turbidity currents in a glaciolacustrine environment can be identified in a number of pro-glacial lakes. When a river reaches peak discharge, turbidity currents entering the lake can create trenches at the bottom. In a study carried out by Lamoureux et al. (2002) on Bear
Lake in Devon Island, large vortex-like underflows eroded and created trenches in the deep basin as seen today. Underflows deposited rhythmic sediments and these were dominant in the proximal reaches. Interflows, associated with a decrease in energy and settling from suspension, were found to be the main sedimentation processes in the more shallow and distal locations of the lake. Turbidity processes were recorded in Atlin Lake, B.C. by Gilbert et al. (2006). They noted that sediment ponding along the thalweg of a lake may be due to turbidity currents following the deepest part of the lake. However, the greater depths of Atlin lake resulted in more isolation from the vigorously circulating interflows which created an evenly distributed and quieter form of deposition. In terms of sedimentary deposition, Desloges & Gilbert (1994) noticed that sediments found in the large basin of Lillooet Lake, B.C. were unconformable, which indicated that the turbidity currents originating off the delta had travelled as underflows and as a result, eroded and flattened the underlying surface. Areas where sediments were found to be more conformable, suggested lower energy conditions and the deposition of sediments was through the process of settling from suspension.
26
2.5 Varves
One of the first analyses of varve chronologies was by de Geer in 1912. de Geer (1912) suggested that annually laminated lacustrine sediment deposits from deep-water lakes have the ability to provide us with detailed and continuous chronological data about the rate of sediment production from a contributing drainage basin (Ohlendorf, 1997; Dirszowsky et al., 2002).
Regional studies of varve thicknesses carried out within Cordilleran lakes by Desloges et al.
(1994, 1998, 2002), Lamoureux et al. (1996), Slaymaker (2003), Dirszowsky et al. (1997, 2004), and Menounos (2004, 2006) have all demonstrated that the composition and deposition of lake sediments in glacial environments are closely related to climatic changes within that region.
Studies have also demonstrated varve formations can be used to infer seasonal or sub-annual climatic changes leading back to the late Holocene (Last, 2002; Hodder et al., 2006). Varve analysis is therefore a potentially important parameter for measuring environmental change in a drainage basin (Desloges et al., 2002; Menounos, 2006).
2.5.1 Varve formation
The formation of varved structures depends greatly on the type and characteristics of the lake such as its tropic state, its morphology, and the physical setting of the lake (Serink, 2004).
Under regular circumstances, a single varve will consist of two distinctive layers forming a couplet. The first is a thick and light coloured layer of silt and sand rich laminae, which forms in the summer through inflow pulses associated with, melt water during the nivial or subsequent rain events. Turbidity currents may deliver part of these higher energy laminae. A thin dark - coloured layer of silty-clay forms during the autumn and winter periods. This reflects the decrease in the sediment supply and quieter water conditions under a frozen lake surface. Low water temperatures delay the settling of clay particles, resulting in the formation of clastic varves 27
which are thinly laminated (Desloges & Gilbert, 1994; Lamoreux, 1996; Slaymaker et al., 2003).
Below average thickness in a varve sequence may be indicative of cooling and / or limited sediment supplies. Extreme flood events from high temperature induced snow and ice melt and high rates of precipitation are noticeable factors that contribute to the formation of thick and above average varves.
Desloges & Gilbert (1994) noted that smaller runoff events, which follow larger flows, might reflect a hysteresis in the supply of sediment. As the melt season progresses, sediment exhaustion leads to thinner, less frequent and finer grained laminations. The thickening of a varve indicates a period of increased sediment delivery rates. Extreme events in particular lead to a reduction in storage and consequently an increase in sediment supply. It is thought that availability and sediment transport capacity would both be at their highest just as glacier retreat begins (Desloges & Gilbert, 1994; Dirszowsky, 2004).
2.5.2 Varves as inferences of past climates
Although numerous studies have demonstrated statistically significant relationships between hydro-climatic records and varve thicknesses, Hodder et al. (2007) question the universal accuracy of this proxy analysis and refer to many “marginal correlations” between varve thicknesses and climate. Using a previously studied glacier-fed lake in the Cordillera of
British Columbia, Hodder et al. (2006), found that, despite the overall physical connection between varve correlations thicknesses and seasonal temperatures and discharge (Fig. 2.1), there is no significant information regarding the processes by which temperature anomalies translate into thicker varves. The “marginal” correlations found at Mud Lake, B.C. demonstrate an incomplete understanding of drainage basin sediment cascades. Hodder et al. (2007) conclude
28
that each varve is not a unique situation as numerous non-linear relationships and feedbacks exist among processes. These situations are either unknown or stochastic.
29
III – Study Area
3. 1 Physiography of the Yukon Territory
The Yukon Territory is within the Canadian Cordilleran region and topography ranges from rolling uplands in the interior to the rugged Coast and St. Elias Mountains (Hughes, 1987).
The Dezadeash region of the southwestern Yukon Territory encompasses three distinct physiographic subdivisions of the Cordilleran region. The Boundary Ranges of the Coast
Mountains occupy the southeast region; the St. Elias Mountains cover the southwest area; and the southeast part of the Kluane Plateau, which is part of the greater Yukon Plateau (Fig. 3.1)
(Kerr, 1934; Kindle, 1953). The Shakwak Trench is a structural trench created by the Denali fault line. The Shakwak Trench extends northwest along the eastern side of the St. Elias
Mountain ranges through to the Kluane Lake drainage basin (Gilbert, 2004) and separates the
Kluane Plateau from the St. Elias Mountains. The St. Elias Mountains are predominantly comprised of the Ice Field Ranges, which lie west of the Alsek River, and the Alsek Ranges, which forms the front of the St. Elias Mountains. The Duke Depression, which occupies the southwest corner of the Dezadeash Range, trends northwestwardly and separates the Kluane
Ranges from the St. Elias Mountains (Fig. 3.1) (Kerr, 1934; Kindle, 1953).
The Takhini River is a tributary of the Yukon River west of Whitehorse. It drains northwest from the Coast Mountains and the Yukon Plateau (Hamilton, 1995). The upper
Takhini River drainage basin straddles the border between the Boundary Range Mountains of southern Yukon and northern British Columbia, with a very small segment in southeast Alaska
(Fig. 3.1, 3.2). The upper Takhini River basin outlet that flows into Kusawa Lake is 60 km east of Haines Junction and 37 km east of the Dezadeash Range. Kusawa Lake (60° 19' 55” N, 136°
4' 48” W), which means “long, narrow lake” in the Tlingit language, is 65 km southwest of
30
Whitehorse. Kusawa Lake occupies the central – southern portions of the structural Shakwak
Trench east of the Boundary Ranges within the upper Takhini River basin. Mountains around
Kusawa Lake are as high as 2000 m a.s.l. (Fig 3.1) (Kindle 1953, Gilbert, 2004).
31
Fig. 3.1 Southwestern Yukon showing names referred to in the text. Spot height elevations determined by Gilbert & Desloges (2005). Arrows = major drainage networks into the Upper Takhini River drainage basin 32
Fig. 3.2 The upper Takhini drainage basin of Kusawa Lake (dark blue bold line) and the sub-basin divides (dark blue dashed line). Glaciated regions are shaded in grey. Contour intervals at 1000 m.
33
3.2 Upper Takhini River drainage basin morphology
Kusawa Lake is 70 km in length and has a total area of 142 km 2. The total drainage of the upper Takhini River basin above Kusawa Lake is 4292 km 2. The basin area to lake area ratio is a modest 30.2, according to the index established by Gilbert (2004). Lake water flows from the south to the north. There are three major sub-drainage basins of the upper Takhini River basin: the Primrose River, Kusawa River, and the upper-most Takhini River, which account for
62.50 % of the total upper Takhini River basin (Gilbert, 2004). Two smaller sub-basins: Jo-Jo
Creek and Devilhole Creek, along with small tributaries found on the west side of the lake, account for the 37.5 % of the entire basin. In proximity to the outlet of Kusawa Lake, is an active alluvial fan complex known as the Kusawa Campground alluvial fan.
3.2.1 Primrose River sub-basin
The Primrose River sub-basin covers an area of 1466 km 2, or 34 % of the upper Takhini
River basin. 3.0 % is glacier covered (Table 3.1) and it has significant sediment traps between the glaciers and outlet. Within the Primrose River basin is Primrose Lake (23 km long, 0.7 km wide) and Rose Lake (8 km long, 1 km wide). In addition, within the upper valleys of the
Primrose River basin, an 8 km long sandur named Silt Lake can be found (Gilbert, 2004). Given the presence of sediment traps, and in comparison to the Kusawa and upper-most Takhini River basins, it was assumed that the Primrose contributes less sediment to Kusawa Lake due to its diminutive glacial cover (Table 3.1). Near the Primrose River outlet are a number of minor debris aprons off the valley walls. The Primrose River fan-delta is approximately 0.8 km 2, is terraced and pro-grades into Kusawa Lake creating lake narrows with the opposite shoreline
(Fig. 3.2)
34
3.2.2 Upper-most Takhini River sub-basin
The upper-most Takhini River sub-basin covers an area of 545 km 2 (12.7 % of the entire basin) and is the smallest of the three major sub-basins (Gilbert 2004). However, it has the greatest proportion of glacial cover of 13.9 % (Table 3.1). In proximity to the mouth of the upper-most Takhini basin is Takhini Lake (9 km long and 0.6 km wide at 800 m a.s.l.). Takhini
Lake probably acts as a significant sediment trap. A very small delta of approximately 0.4 km 2 is found at the mouth of the upper-most Takhini River as it flows into Kusawa Lake (Fig. 3.2).
3.2.3 Kusawa River sub-basin
The Kusawa River sub-basin covers an area of 677 km2 (15.8% of the upper Takhini basin). Glaciers cover an area of 75 km 2 (or 11.2%) of the sub-which is 1.8 % of the entire
Upper Takhini River basin. Unlike the Primrose and the upper-most Takhini River sub-basins, the Kusawa River basin does not appear to have significant sediment traps. However, there is a large (13 km 2) unnamed lake in the headwaters of the Kusawa River. This region of the lake denotes rising valley sides, which are narrow (Gilbert, 2004). The Kusawa River delta is confined to the steep valley sides and is 6 km long. The nature of the delta is braided and the presence of eyots suggests that the Kusawa River acts as a primary sediment source into Kusawa
Lake (Fig. 3.2).
3.2.4 Jo-Jo Creek sub-basin
Jo-Jo Creek sub-basin is the most distal sub-basin and covers an area of 330 km 2 (7.7 % of the entire basin). Jo-Jo Creek sub-basin has no glacial cover. Jo-Jo Lake found at 953 m a.s.l. probably forms a significant sediment trap. The overall length of Jo-Jo Lake is ~12 km but its average width is 1.6 km as it occupies a narrow, steep-walled valley. A massive kame moraine, 35
~ 70 m high found 3.5 km north of the lake, and stratified gravel, sand and silt deposits found south of the lake, suggest that Jo-Jo Lake was once dammed (Fig. 3.2).
3.2.5 Devilhole Creek sub-basin
The Devilhole Creek sub-basin is situated west of Kusawa Lake and covers an area of
244.5 km 2 (5.7% of the entire basin). Like the Kusawa River sub-basin, the region denotes steep valley sides. Glaciers occupy 5 km 2 (Table 3.1) of the headwaters and a number of moraines can be found (Fig. 3.2).
3.2.6 Kusawa Campground alluvial fan-delta complex
The total catchment the alluvial fan-delta complex occupies an area of 32 km 2 (Lipovsky,
2006), or 0.0075 % of the entire Upper Takhini River basin. The fan and delta features act as the primary sediment source into Kusawa Lake (Fig. 3.2).
Recent field studies conducted by Lowery (2002), Yukon Parks and the Yukon Geological
Survey, showed that sediment in the upper valleys and on exposed land at the northern end of the lake is comprised of a complex of kames, kame terraces and kettles. Sedimentary debris deposited in the valley floor of the Kusawa floodplain consists of gravelly muddy sand: 20% gravel, 55 % sand, 19 % silt and 6 % clay (Lowey, 2002). The composition of the terrain around
Kusawa Lake varies with changing elevation. Slopes at high elevation display steep rock outcrops, occasionally covered by thin moraine or colluvial deposits. Moraine deposits, with permafrost present in local areas, cover slopes at mid-elevation. At lower elevations and on valley floors it is comprised of glaciofluvial sand and gravel. Colluvial deposits found in higher elevations consist mainly of cobbley and bouldery diamicton (a finer grained matrix with permafrost common at shallow depths), which is susceptible to creep, solifluction, gullying and 36 active layer detachment flows. The till normally found at mid-elevations is also composed of diamicton, yet with a higher amount of sand, silt and clast content (Lowey, 2002; GeoProcess
NTS 115A, 2002).
3.3 Bedrock and surficial geology
The bedrock geology of Kusawa Valley consists primarily of crystalline rocks of the
Coast Plutonic Complex, essentially comprised of 100 – 55 million year old (M.Y.) granodiorite and granite intrusions. The Yukon – Tanana Terrane includes pre – 570 M.Y quartz mica schist, gneiss, slate, quartzite, limestone, and greenstone. The occurrence of mineral deposits within these rock types include platinum group elements, copper, gold, silver, lead and zinc most commonly found along the region near the Denali Fault line and towards the east of Kusawa
Lake (GeoProcess NTS 115A, 2002).
Glacial, glaciofluvial, colluvial, and fluvial processes have contributed to the formation of the rugged topography found in the Kusawa Valley. Kames, kame terraces, and kettles are located north of the Kusawa Lake outlet. Steep rock outcrops are dominant at high elevations above 1500 m a.s.l. Thin ground moraine or colluvial deposits which are comprised of cobbly, boulder diamicton may also be found throughout the basin. Slopes are susceptible to creep and solifluction from permafrost effects. Gullying, and under extreme circumstances, avalanches and rockslides may occur. Streamlined or crested moraine deposits of sand and gravel cover the mid elevation slopes (1000 – 1500 m a.s.l.) in the region of active glacier cover. Permafrost may also be present at these elevations and is extremely susceptible to creep and thermokarst (Gilbert,
2004). Glaciofluvial silts, sands, and gravel deposits are abundant at the lower elevations (680 –
1000 m a.s.l.) and throughout the valley floors. The relative coarseness of these sedimentary
37 deposits allows for stable well-drained surfaces. Glaciolacustrine deposits of fine sand, silts, and clays are dominant at lake level (GeoProcess NTS115A 2002; Lowey 2002).
3.4 Glacial cover
The distribution of glacial cover within the upper Takhini River drainage basin has been heavily influenced by coastal – atmospheric interactions. The total glacier area in 2005 is ~ 200 km 2, which accounts for 4.7 % of the total basin area (Table 3.1). The majority of the glaciers are small, with the largest being 8 km in length. The southern region is influenced by cool maritime climate conditions, which has a profound influence upon winter accumulation and summer ablation. As a result, it is assumed that the glaciers are relatively active (Gilbert, 2004).
Sub-basin / Location Glacial Cover % of % of Upper Takhini area (km 2) sub-basin River Drainage Basin Primrose River 44 3.0 1.0 Upper-most Takhini River 76 13.9 1.8 Kusawa River 75 11.2 1.8 Devilhole Creek Headwaters 5 0.07 0.1 TOTAL 200 - 4.7
Table 3.1 Glacial cover distribution within the Kusawa Drainage Basin (Gilbert, 2004)
38
3.5 Hydrology
Hydrometric stations associated with the Takhini River watershed are situated at the
Kusawa Lake outlet (09AC004) and along the lower Takhini River, 40 km downstream of the outlet near Whitehorse (09AC004) (Fig. 3.3). The Kusawa Lake outlet station was commissioned in 1952 and recorded discharge between 1955 – 1986. Since then, it has ceased operation. The downstream gauge, known as Highway Bridge, was installed in 1948 and is still in operation (Table 3.2) (WSC Hydat 2008).
Rivers north of the 60 th parallel typically exhibit a nival hydrologic regime: high flows in the spring and summer, which are associated with glacial and nival-melt and rainstorms respectively, and low flows throughout the winter. Throughout the year, the mean annual
o o discharge (Qm) of the lower Takhini River at Highway Bridge (60 51’ 0” N, 135 43’ 47.9” W)
3 is on average 10 m /s greater than the Qm at the Kusawa Lake outlet (Table 3.2, Fig. 3.3, 3.4 &
3.5) as it drains additional sub-basins: Ibex River, Mendenhall River and Arkell Creek. During the winter, discharge is limited to baseflow at the Kusawa Lake outlet due to ice formation.
Moore, R.D et al. (2002) studied winter streamflow variability at both sites and discovered that under extreme winter conditions (i.e.:1965 and 1994) ice can plug the lower Takhini channel for
2 months and discharge levels do not return to a pre-freeze stage until March. Late spring melt in mid-April caused minor increases in both stage and discharge. However, the marked rise in discharge during late April - May corresponded with a sudden decline in stage. The authors believe that this occurrence was consequence of ice break-up and the associated decrease in frictional flow resistance. The rise in discharge levels was a result of channel storage release caused by backwater effect of ice cover. Throughout May – August, snow and ice melt and storm water feed into Kusawa Lake. Up to 220 m3/s of water can be discharged at the lake outlet
(Moore, R.D et al., 2002b). 39
Takhini River Site coordinates (DMS) Station Drainage Qm Number Area (km 2) (m 3/s) Highway Bridge 60 o 51’ 0” N, 135 o 43’ 47.9” W 09AC001 6990 61.9 Kusawa Lake outlet 60 o 36’ 35.9” N, 136 o 7’ 12.1” W 09AC004 4292 51.7
Table 3.2 Lower Takhini River discharge stations. Mean annual discharges (m 3/s) are based on 1955 – 2005 averages (WSC Hydat 2008)
Fig. 3.3 Mean annual discharge of the lower Takhini River 1950 – 2007 (App. 3A) (WCS 2008)
Peaks in mean annual discharge (Q m) of the lower Takhini River are cyclic in nature. At
3 Highway Bridge, Qm peaks greater than 70 m /s occur approximately every 6.1 years and have a
16.4% chance of doing so in any given year. This is typically followed by a significant decrease
3 to < 60 m /s in Q m within the next year. Qm values at the Kusawa Lake outlet resemble those at
40
3 Highway Bridge although values are ~ 10 m /s lower. Qm values above the overall mean of 51.8 m3/s at the Kusawa Lake outlet have a 57.7 % recurrence probability (Fig. 3.3).
Fig. 3.4 Annual hydrograph (mean daily discharge) for the lower Takhini River at the outlet of Kusawa Lake
(09AC004) for 1986 (WSC 2008). 1986 data are shown in red. The maximum and minimum values for the entire record are shown in green and blue respectively.
41
Fig. 3.5 Annual hydrograph for the lower Takhini River at Highway Bridge (09AC001) 1986 (WSC2008).
1986 data are shown in red. The maximum and minimum values for the entire record are shown in green and blue respectively.
3.6 Climate
The Dezadeash region has a sub-arctic climate. Average temperatures rise above 10 oC for approximately three months of the summer (Table 3.3). The warm air mass from the Pacific moderates temperatures within the region. As Kusawa lies within the regions of the Boundary
Range Mountains, temperature and precipitation received is considerably lower due to the orographic effects from the Coastal and St. Elias mountain ranges (Lowey, 2002).
Temperatures in January can range between -22 oC to -13.3 oC with an average of -17 oC. In July, the average temperature is 14 oC although temperatures may dip to 7 oC and peak to 20.5 oC
(Environment Canada, 2008).
The total annual precipitation is low and averages 267 mm. Most precipitation occurs in the winter rather than in the summer. The average depth of winter snow is 22 cm. The months of July and August are considered some of the wettest months with an average rainfall of 26 mm 42
(Table 33, Fig. 3.6, 3.7 & 3.8) (Pienitz, 1997; Environment Canada, 2008). Lakes in the southern part of Yukon have protracted periods of ice cover. This usually extends from
November to May, with the mean number of days between the clearing of ice and initial ice formation being around 150 days in the Whitehorse area (Pienitz et al., 1997).
Average Average Average Total Total Total Average Rainfall Snowfall Precipitation Month Temp ( oC) (mm) (cm) (mm) Jan -17.7 0.2 23.7 16.7 Feb -13.7 0.1 16.8 11.4 March -6.6 0 14.9 10.4 April 0.9 1.3 8 7.0 May 6.9 13 2.4 15.2 June 11.8 29.7 0.5 30.3 July 14.1 41.4 0 41.4 Aug 12.5 38.5 1 39.4 Sept 7.1 29.3 5.1 34.1 Oct 0.6 8.8 18.9 23.8 Nov -9.4 0.7 27.3 19.2 Dec -14.9 0.3 26.4 18.5 Year -0.7 163.1 145 267.4
Table 3.3 Climatic characteristics for Whitehorse based on 1971 – 2001 averages
Lat / Long: 60° 43' N, 135° 4' W. Elevation: 706.2 m. Climate ID: 2101300. WMO ID: 71964
(Environment Canada 2008).
43
Fig. 3.6 Temperature precipitation graph of Whitehorse based on 1971 – 2000 averages.
(Environment Canada 2008).
Fig. 3.7 Monthly rainfall normals for the Takhini River Ranch 1971 – 2000 averages
(Environment Canada 2008).
44
Fig. 3.8 Monthly snowfall normals for the Takhini River Ranch 1971 – 2000 averages
(Environment Canada 2008).
3.7 Vegetation
The outlet of Kusawa Lake is at 750 m a.s.l. Due to the influence of both coastal and arctic air masses, conifers and some mixed forest trees are dominant up to 1300 m a.s.l. In the very compressed sub-alpine zone, white spruce and fir are present although white spruce is the dominant tree species in the valleys. An under-story of feather-moss in the forested areas is also common along the valley floor. Balsam poplar may also be found on the margins of the lake.
Above tree line (~ 1200m a.s.l) herbs, dwarf shrubs, shrub birch, willow, moss, and lichen dominate the alpine tundra terrain (GeoProcess NTS 115A, 2002).
45
IV – Methods
4.1 Field methods
4.1.1 Acoustic profiling
During the summer of 2004, all of Kusawa Lake was surveyed using a Datasonics Dual –
Frequency CHIRP II sub-bottom acoustic profiler. Sixty-eight oblique cross – lake transects from the northern (outlet – Takhini River outlet from Kusawa Lake) to the southern reaches of the lake (inlet of Kusawa River) were completed covering 125 km worth of line work. The profiler can penetrate up to 100 m of thick silt and clay muds and produce sub-bottom imagery with a maximum resolution of 0.5 m. Data derived from the acoustic records allowed for the mapping of lake bathymetry and sediment depths (Gilbert & Desloges 2005).
4.1.2 Sediment cores
Nine cores in the northern region of the lake and two cores in the central and southern regions were extracted from the lake using a Rossfelder submersible vibra corer (Table 4.1).
Core tubes were 6 m in length, 0.75 cm in diameter and were made of aluminium. A core catcher was riveted into the bottom of the tube and the top was inserted into the vibra corer head.
Cores were extracted at mid-points in the lake to avoid local slope effects. The lengths of cores obtained ranged from 0.8 to 5.5 m and were cut into smaller sections for shipment.
4.1.3 Conductivity Turbidity and Density Profiles
CTD profiles of Kusawa Lake were measured using a Hydrolab DataSonde 4a. Sixteen sites at approximately 1.5 to 1.6 km intervals from the source to the outlet of the Lake were surveyed during the summer of 2004 measuring primarily temperature, turbidity and conductivity with depth. Identifying where water columns become unstable is an integral part of 46 understanding how Kusawa Lake mixes vertically (James, 2004). Regions of mixing can be easily identified from density profiles where density inversions occur. How this changes throughout the long profile of the lake will also determine sediment distribution as abrupt changes in lake temperature variations and / or inflowing water have the ability to create turbidity currents and induce sedimentation (Serink, 2004).
4.1.4 Geographic Information Systems
A GARMIN Geographic Positioning System (GPS) recorded acoustic track lines, coring locations and CTD measurement locations. The metadata was transferred onto the Earth Science
Research Institute’s (ESRI) ArcGIS and visually displayed using ArcMap. The ArcMap file contains recorded data indicating the water depth, length of core, maximum thickness of sediment from the sub-bottom acoustic records and surface water temperature of each coring location (Fig. 4.1; Table 4.1)
47
Fig. 4.1 Kusawa Lake vibra core sample locations
48
Vibra Site Coordinates Lake Depth Surface Vibra core core # (D.M.S) (m) Water Temp length (m) (oC) 318 60 o 20’ 15.216” N, 136 o 4’ 39.755” W 30 .0 11.1 2.14 319 60 o 19’ 10.847” N, 136 o 4’ 55.559” W 28.4 11.6 3.34 322 60 o 15’ 35.423” N, 136 o 6’ 29.627” W 47 .0 13.8 1.10 324 60 o 10’ 11.639” N, 136 o 10’ 5412” W 135 .0 11.6 2.78 325 60 o 18’ 21234” N, 136 o 4’ 22.583” W 24.7 12.1 2.87 326 60 o 18’ 32.327” N, 136 o 4’ 42.239” W 30.2 13.3 2.96 327 60 o 19’ 17.759” N, 136 o 4’ 32.663” W 29.6 13.4 2.70 328 60 o 19’ 52.860” N, 136 o 4’ 48.071” W 32.4 11.0 1.70 32 9 60 o 21’ 18.971” N, 136 o 4’ 42.780” W 32.2 10.8 5.48 330 60 o 20’ 54.239” N, 136 o 4’ 35.940” W 34.7 10.8 0.80 331 60 o 21’ 4.7874” N. 136 o 4’ 45.551” W 33.5 10.9 2.90
Table 4.1 Site characteristics and core sample lengths taken in July 5 – 8th 2006.
4.2 Laboratory methods
4.2.1 Sediment core properties
In the laboratory, core sections were stood upright to separate interstitial water. Any surfacing water was removed with the use of a pipette. The cores were split longitudinally by incising the aluminium tube and passing a steel wire through to separate the mud. The lengths of the sediment in the cores were measured and recorded. One-half of the core was wrapped in saran wrap and aluminium foil and stored in the cooler at 4 oC for archival purposes.
The different stages of drying (which were primarily controlled by temperature and moisture content) were logged through digital photography. The process of drying allowed the rhythmic laminations to become clearly visible (Lamoreaux, 1996). Photographs of the cores in a semi- dried state were developed and manually cross-referenced with other cores. The physical characteristic of each core (i.e.: grain size, volcanic ash layers, organics, colour etc.) were logged manually. This was graphically displayed on Corel Draw 12 as a schematic core section for visual interpretation. 49
Couplets of silt overlain by clay were analysed extensively through the measurement of thickness. Counts were conducted from high – quality photographic close-ups. Where the sections of the core showed couplets less than 2 mm in thickness, the use of a high precision microscope with x40 the magnification (the Increment Measurer System) was calibrated to conduct couplet counts. Thickness and couplet counting were carried out as consistently as possible but the possibility of error arises from the possible presence of sub-annual (sub-event) layers.
4.2.2 Loss on Ignition
The Loss on Ignition (LOI) method was used to determine the percentage of organic matter and carbonate content (Dean, 1974; Heiri, et al., 2001). LOI was carried out with approximately 2 g of sediment sample taken from all different sedimentary units within the cores. Samples were oven dried at 105 oC for 24 hours to exonerate excess moisture and re- weighed. The first reaction occurs when the organic matter is oxidized after burning at a temperature of 550 oC to carbon dioxide and ash for one hour. Following a cooling period the samples were then reweighed to determine the percentage loss:
OM = LOI 550 {[MS 105 – MS 550 ] /MS 105 }*100
Where:
OM = Organic Matter %
C = Carbonate %
o LOI 550 = Loss on Ignition at 550 C
o MS 105 = Mass of Dried Sample at 105 C
o MS 550 = Mass of Dried Sample after LOI at 550 C
50
Given the bedrock geology, no major sources of carbonate content were expected.
However, the LOI procedure on determining the carbonate content was followed anyway to test error on the OM content.
The second reaction occurred when samples were burned at a temperature of 1100 oC and here, the carbonate was lost leaving oxide residuals. The sample was reweighed and the percentage loss from the initial sample weight was calculated:
C = LOI 1100 {[MS 550 – MS 1100 ] / MS 550 }*100
Where:
OM = Organic Matter %
C = Carbonate %
o LOI 550 = Loss on Ignition at 550 C
o LOI 1100 = Loss on Ignition at 1100 C
o MS 105 = Mass of Dried Sample at 105 C
o MS 550 = Mass of Dried Sample after LOI at 550 C
o MS 1100 = Mass of Dried Sample after LOI at 1100 C
(Dean, 1974; Håkanson 1995; Heiri et al. 1999)
Results obtained from the LOI procedure were used as a proxy for the clastic content of lake sediments. In lakes draining glacial environments, changes in LOI appear to reflect the changes in the fraction of clastic sedimentation (Menounos et al. 2004) where sediment sequestration occurs. Periods of high OM content may reflect possible periods of high deposition from a flooding and / or outwash event or changes in lake primary productivity levels
(Håkanson 1995; Heiri et al. 1999; Lacourse & Gajewski, 2000).
51
4.2.3 Laser particle analysis
Evaluating the sediment particle size is an important textural parameter of glaciolacustrine sedimentology as it supplies information on the conditions of transport, sorting, and deposition of the sediment. In addition, it provides some information on the history of events such as floods that occurred at the depositional site (Lowey 2002; Schnurrenberger et al.,
2003). Determining the percentage of clay and silt allows for possible inferences of inputs from glaciolacustrine deposits in the Kusawa Lake drainage basin.
Laser Particle Size (LPS) Analysis was conducted using two different instruments from two different laboratories: The Coutler Counter from the Department of Geography, Queen’s
University and the Malvern Laser Sizer from the Department of Chemical Engineering,
University of Toronto. The Coulter Counter was used for cores KUS318, 330, 331, and 329.
The Malvern Laser Sizer was used for cores KUS319, 322, 324, 325, 326, 327 and 332.
Sedimentary units within the core section, which displayed distinctive thick silt and clay deposits, thick sand events or a distinctive amalgamation of sandy fine silt over a sand layer, were selected for grain size distribution analysis. Approximately 20 samples were taken from each core with the exception of cores KUS330 and 332. Ten samples were taken from these cores, as cores were 0.8 and 1.0 m in length respectively.
The best form of data regarding fine particle size distributions, hence the clay – silt fraction, was obtained from samples which do not contain any organic matter. Approximately
250 mg of clay – silt and 150 mg of fine sand fractions were placed in glass vials and 7 pipette drops of 32.5 % concentrated hydrogen peroxide (H 202) were used in each sample to remove any organic matter. Samples were left to evaporate and dry for one week in the fume hood.
The sample was mixed with distilled water and further disaggregated using ultrasonics when passed through the refractor. The Coulter and the Malvern instruments were set to read
52 grain size intervals ranging from 0.375 m to 2000 m and from 0.00582 m to 880 m respectively. The mean, median, mode, standard deviation, variance, skewness in a normal distribution, kurtosis, D 10 , D 50 and D 90 and the specific surface area of a particle were all recorded. A replicate of three readings were conducted for each sample. The standard cut off for fluvial geomorphic classification of particle size was used: clay (4 m), Silt (63 m), Sand (2000