Dropstones in Lacustrine Sediments as a Record of Snow Avalanches—A Validation of the Proxy by Combining Satellite Imagery and Varve Chronology at Kenai Lake (South-Central Alaska)

Snow avalanches cause many fatalities every year and damage local economies worldwide. The present-day climate change affects the snowpack and, thus, the properties and frequency of snow avalanches. Reconstructing snow avalanche records can help us understand past variations in avalanche frequency and their relationship to climate change. Previous avalanche records have primarily been reconstructed using dendrochronology. Here, we investigate the potential of lake sediments to record snow avalanches by studying 27 < 30-cm-long sediment cores from Kenai Lake, south-central Alaska. We use X-ray computed tomography (CT) to image post-1964 varves and to identify dropstones. We use two newly identified cryptotephras to update the existing varve chronology. Satellite imagery is used to understand the redistribution of sediments by ice floes over the lake, which helps to explain why some avalanches are not recorded. Finally, we compare the dropstone record with climate data to show that snow avalanche activity is related to high amounts of snowfall in periods of relatively warm or variable temperature conditions. We show, for the first time, a direct link between historical snow avalanches and dropstones preserved in lake sediments. Although the lacustrine varve record does not allow for the development of a complete annual reconstruction of the snow avalanche history in the Kenai Lake valley, our results suggest that it can be used for long-term decadal reconstructions of the snow-avalanche history, ideally in combination with similar records from lakes elsewhere in the region.


Introduction
Over the period 2006 to 2016, an average of 27 people died each year in the United States of America due to snow avalanches [1]. In the European Alps, the number of fatalities is even higher, with a yearly average of 100 casualties over the past 40 years [2]. Furthermore, snow avalanches also have socio-economic consequences, such as road blockages, and power plant and power line damage [3]. Ongoing climate change affects the air temperature in mountainous regions and influences River and Snow River, along the east side of the lake, are the main inlets and sources of water and sediment. A smaller inlet is situated at the north side of the lake [13].
The Kenai Peninsula is situated above the Alaska-Aleutian subduction zone, where numerous active volcanoes are present, and earthquakes occur frequently [16]. These earthquakes and volcanic eruptions can result in the deposition of turbidites and tephras in lakes that can then be used as time markers. In Kenai Lake, the 1964 Great Alaska Earthquake (M W 9.2) triggered deposition of a turbidite [17], which formed the base of the interval that is studied here. Within this interval, Boes et al. [13] suggested that one of the peaks in the magnetic susceptibility curve, taken along a sediment core, is related to the 1989/1990 Redoubt Volcano eruption.
Kenai Lake sediments are composed of clastic varves, representing a couplet of a coarse-grained spring-summer lamina and a fine-grained winter lamina [13,18]. Boes et al. [13] observed internal laminations in some varves that could be mistaken for a new varve year. These laminations were interpreted to be related to flood events in the late summer or to spring melting of large amounts of snow/ice in a short period. The Kenai Peninsula is situated above the Alaska-Aleutian subduction zone, where numerous active volcanoes are present, and earthquakes occur frequently [16]. These earthquakes and volcanic eruptions can result in the deposition of turbidites and tephras in lakes that can then be used as time markers. In Kenai Lake, the 1964 Great Alaska Earthquake (MW 9.2) triggered deposition of a turbidite [17], which formed the base of the interval that is studied here. Within this interval, Boes et al. [13] suggested that one of the peaks in the magnetic susceptibility curve, taken along a sediment core, is related to the 1989/1990 Redoubt Volcano eruption.
Kenai Lake sediments are composed of clastic varves, representing a couplet of a coarse-grained spring-summer lamina and a fine-grained winter lamina [13,18]. Boes et al. [13] observed internal laminations in some varves that could be mistaken for a new varve year. These laminations were interpreted to be related to flood events in the late summer or to spring melting of large amounts of snow/ice in a short period.

Multibeam Data, Core Collection and Core Analysis
Multibeam bathymetry data for Kenai Lake were collected in the summer of 2015 using an ELAC SeaBeam 1050 system (ELAC Nautik GmbH, Kiel, Germany). This sonar system has two 50 kHz transducer arrays and a 120° swath angle. Roll, pitch, and heave were corrected by an IXSEA Octans 3000 motion sensor (iXBlue, Cedex, France). CTD measurements (conductivity, temperature and depth) at the lake surface were carried out with a Valeport Modus CTD. A CastAway CTD and a Sea and Sun Technology CTD 48M were used for sound velocity profiles, measured several times each day. The echosounder was controlled by ELAC HydroStar software, which also served as an interface for external sensors that acquired hydrographic information. Processing of the bathymetric data was done with HDPedit (ELAC Nautik GmbH, Kiel, Germany)(outlier cleaning), HDPpost (ELAC Nautik GmbH, Kiel, Germany) (correcting of sound velocity profiles) and IVS 3D Fledermaus™ (QPS B.V., Zeist, The Netherlands) (detailed ping editing). The final high-resolution bathymetry map has a 9-mgrid spacing and a vertical accuracy of a few decimeters.
A total of 27 sediment cores were retrieved during the summer from Kenai Lake using a gravity corer with hammering extension; 9 in 2012 and 18 in 2015. After shipping to Belgium, the cores were split at Ghent University, the working halves were described macroscopically, and geophysical

Multibeam Data, Core Collection and Core Analysis
Multibeam bathymetry data for Kenai Lake were collected in the summer of 2015 using an ELAC SeaBeam 1050 system (ELAC Nautik GmbH, Kiel, Germany). This sonar system has two 50 kHz transducer arrays and a 120 • swath angle. Roll, pitch, and heave were corrected by an IXSEA Octans 3000 motion sensor (iXBlue, Cedex, France). CTD measurements (conductivity, temperature and depth) at the lake surface were carried out with a Valeport Modus CTD. A CastAway CTD and a Sea and Sun Technology CTD 48M were used for sound velocity profiles, measured several times each day. The echosounder was controlled by ELAC HydroStar software, which also served as an interface for external sensors that acquired hydrographic information. Processing of the bathymetric data was done with HDPedit (ELAC Nautik GmbH, Kiel, Germany)(outlier cleaning), HDPpost (ELAC Nautik GmbH, Kiel, Germany) (correcting of sound velocity profiles) and IVS 3D Fledermaus™ (QPS B.V., Zeist, The Netherlands) (detailed ping editing). The final high-resolution bathymetry map has a 9-m-grid spacing and a vertical accuracy of a few decimeters.
A total of 27 sediment cores were retrieved during the summer from Kenai Lake using a gravity corer with hammering extension; 9 in 2012 and 18 in 2015. After shipping to Belgium, the cores were split at Ghent University, the working halves were described macroscopically, and geophysical parameters were logged at 0.2 cm resolution using the Geotek Multi-Sensor Core Logger (Geotek Ltd, Northamptonshire, United Kingdom) (MSCL). The MSCL is equipped with a Cs gamma source and detector to calculate gamma density, a Konica Minolta CM-2600d spectrophotometer (Konica Minolta Sensing Americas Inc, Ramsey, New Jersey) to analyze colour spectra and a Bartington MS2E point sensor (Barington Instruments Ltd, Oxon, United Kingdom) to measure magnetic susceptibility. Photographs of the split core surfaces were acquired using the line scan camera mounted on the MSCL.

Core Scanning and Debris Detection
All cores were imaged with a medical X-ray computed tomography scanner (Siemens SOMATOM Definition Flash (Siemens Healthineers, Erlangen, Germany) at the Ghent University Hospital. The scanner was employed at 120 kV, an effective mAs of 200 and a pitch of 0.45. The helical scans were reconstructed using two different protocols. A "soft tissue" reconstruction produced smoother images that have a relatively strong contrast between sediments with low X-ray attenuation. They were used to image the varves in the cores. The "bone" reconstruction produces sharper (but more noisy) images, creating sharp boundaries between objects that have a highly contrasting X-ray attenuation. They were used to image and analyse the medium to coarse sand and gravel in the muddy matrix.
The 2012 cores were split before scanning; while those collected in 2015 were scanned as whole cores. The reconstructed DICOM (Digital Imaging and Communications in Medicine) images have a resolution of 0.15 mm/pixel and a thickness of 0.6 mm and were stacked to create a 3D volume. The step size between consecutive images was 0.6 mm for the 2012 cores, but 0.3 mm for the 2015 cores. Consecutive images of the latter, thus, overlap 50%. The resulting voxel volumes are~0.014 and~0.007 mm 3 for 2012 and 2015 cores, respectively. The stacks were imported in VG Studio (2.2) for visualisation and were converted to tiff images in ImageJ (64-bit Java 1.6.0_24; [20]). The tiff images of each core were uploaded in Octopus Analysis (formerly known as Morpho+) [21]. This software allows a 3D-analysis of the cores and the segmentation and isolation of groups of voxels, which represent sand or gravel, based on their specific Hounsfield Units (HU), i.e., characteristic grey values which are related to their specific X-ray attenuation [12,22]. A dual threshold of 2200 to 2250 HU in the 2015 cores and 2350 to 2400 HU in the 2012 cores was used to select and isolate sand or gravel. The higher threshold values for the 2012 cores were determined by comparing the HU histograms and are the result of the drier state of these split cores compared to the closed 2015 cores. For each isolated group of voxels (i.e., sand or gravel grains), the equivalent spherical diameter (ESD) was calculated. The ESD represents the diameter of an irregular shaped object calculated from a sphere with the same volume as the aspheric particle [23]. The minimal grain diameter that theoretically can be distinguished (i.e., 0.24 mm) is affected by the minimal voxel volume of the CT-images (i.e.,~0.007 mm 3 ). Following Nesje et al. and Vasskog et al. [7,10], we refer to the particles with an ESD larger than 0.25 mm (medium sand) as debris particles. For each core, we analyzed both the number of particles present and the ESD of each individual isolated particle ( Figure 2). The ESD of the debris particles was used to study their average grain size. The number of particles found in each layer represents a minimum, as some grains are too small or are located too close to the core liner, which makes it very difficult or impossible to isolate this group of voxels. The stratigraphic position of each grain was linked to a particular varve. Based on the CT-images and the results of the grain-size analysis of the debris, two cores (i.e., KE15-07 and KE15-08) were selected for sampling of grains of different sizes ( Figure 2). These fragments were visually examined to identify the rock type, which can help to assign a possible source region to the debris. Frontal cross -section of the "soft tissue" reconstructed image. X-ray computed tomography (CT) "sum-along-ray" (equivalent to inverted radiograph) image of the "bone" reconstructed image stack of the same core; laminated (varved) mud is shown in grey scale (intermediate HU), debris in blue (high HU), and organic material in green (low HU). Each layer that contains debris in this core is labelled by DL (debris lamina) and the layer number. The number of grains per equivalent spherical diameter (ESD) class is shown on the right. Pictures of selected grains from three layers.

Core Correlation, Varve and Tephra Chronology
The cores containing debris were correlated by marker horizons and the sequences of varves present in each core. Varves containing debris were marked as debris laminae (DL) followed by the number of the layer, starting at the sediment-water interface and counting downwards to the top of the turbidite deposited in response to strong shaking during the 1964 earthquake [17]. Dating of the different debris laminae was based on the varve chronology, corroborated by tephra marker beds. The different varves are formed in varve years (VY), which start at the beginning of the melting season (spring) and end at the end of winter [13,18], e.g., VY2014 starts approximately in May 2014 and ends in April 2015.
The varve chronology of the different cores is based on varve counting and chronology by Boes et al. [13]. This original chronology is correlated to historical volcanic eruptions based on newly acquired glass geochemistry data from two cryptotephras indicated by sharp MS peaks ( Figure 2). The uncertainty of the varve chronology between marker horizons was calculated using the presence of uncertain varves that were retained for the chronology (a potential "erroneous" varve [17]) and uncertain varves that were not used in the chronology (a potential "missing" varve [17]). From each marker horizon downward, minimum and maximum ages were calculated by summing "erroneous" and "missing" varves, respectively. From each marker horizon upward, minimum and maximum ages were calculated by summing "missing" and "erroneous" varves, respectively. As a result, each lamina receives two minimum and two maximum ages (i.e., one relative to both the nearest marker horizon above and below). From those, the minimum and maximum ages representing the lowest uncertainty were retained.
Cryptotephra samples were treated with 35% hydrogen peroxide to remove organic material and washed through a 20 μm sieve. Glass was separated using LST (lithium heteropolytungstate) heavy liquid at ~2.45 g/cm 3 and mounted in acrylic pucks using SpeciFix-20 epoxy. All analyses were by wavelength-dispersive spectrometry (WDS) on a JEOL 8900 SuperProbe (JEOL USA Inc, Peabody, Frontal cross -section of the "soft tissue" reconstructed image. X-ray computed tomography (CT) "sum-along-ray" (equivalent to inverted radiograph) image of the "bone" reconstructed image stack of the same core; laminated (varved) mud is shown in grey scale (intermediate HU), debris in blue (high HU), and organic material in green (low HU). Each layer that contains debris in this core is labelled by DL (debris lamina) and the layer number. The number of grains per equivalent spherical diameter (ESD) class is shown on the right. Pictures of selected grains from three layers.

Core Correlation, Varve and Tephra Chronology
The cores containing debris were correlated by marker horizons and the sequences of varves present in each core. Varves containing debris were marked as debris laminae (DL) followed by the number of the layer, starting at the sediment-water interface and counting downwards to the top of the turbidite deposited in response to strong shaking during the 1964 earthquake [17]. Dating of the different debris laminae was based on the varve chronology, corroborated by tephra marker beds. The different varves are formed in varve years (VY), which start at the beginning of the melting season (spring) and end at the end of winter [13,18], e.g., VY2014 starts approximately in May 2014 and ends in April 2015.
The varve chronology of the different cores is based on varve counting and chronology by Boes et al. [13]. This original chronology is correlated to historical volcanic eruptions based on newly acquired glass geochemistry data from two cryptotephras indicated by sharp MS peaks ( Figure 2). The uncertainty of the varve chronology between marker horizons was calculated using the presence of uncertain varves that were retained for the chronology (a potential "erroneous" varve [17]) and uncertain varves that were not used in the chronology (a potential "missing" varve [17]). From each marker horizon downward, minimum and maximum ages were calculated by summing "erroneous" and "missing" varves, respectively. From each marker horizon upward, minimum and maximum ages were calculated by summing "missing" and "erroneous" varves, respectively. As a result, each lamina receives two minimum and two maximum ages (i.e., one relative to both the nearest marker horizon above and below). From those, the minimum and maximum ages representing the lowest uncertainty were retained.
Cryptotephra samples were treated with 35% hydrogen peroxide to remove organic material and washed through a 20 µm sieve. Glass was separated using LST (lithium heteropolytungstate) heavy liquid at~2.45 g/cm 3 and mounted in acrylic pucks using SpeciFix-20 epoxy. All analyses were by wavelength-dispersive spectrometry (WDS) on a JEOL 8900 SuperProbe (JEOL USA Inc, Peabody, Massachusetts) at the University of Alberta using 15 KeV, 6 nA current and 5 µm beam. The software 'Probe for EPMA' (Probe Software Inc, Eugene, Oregon) was used to apply time-dependent intensity corrections for Na loss induced by the smaller beam size [24]. Two natural glasses, Lipari obsidian (ID3506) and Old Crow tephra, were analyzed before and after the samples to track calibration and the quality of analyses.

Satellite Image Analysis
Satellite images were used to visually verify the occurrence of snow avalanches in winters over the period 2004 to 2014. The location and timing of these avalanches were compared to the ages of debris laminae found in the lake sediments. For other years, we were not able to locate satellite images with sufficient resolution.

Climate Data
The Alaska Climate Research Center provided climate data for the period 1963 to 2013 (for winter temperature and snowfall) collected at two stations located near Kenai Lake. These two stations were Cooper Lake Project (60 • 23'N, 149 • 41'W, 182 m above sea level) and Cooper Landing (60 • 29'N, 149 • 58'W, 105 m above sea level) ( Figure 1).

Avalanches on Satellite Images
We observed snow-avalanche deposits on the five available satellite images over the period 2004 to 2014 along the south shore of the central basin of Kenai Lake (Figures 3 and 4). The majority of these can be found between the Hydroelectric Power Station (HPS) and the Ship Creek Delta (SCD) (Figure 1), while a few are situated east of Ship Creek Delta. We refer to these places as the western part and the eastern part of the central basin, respectively. In the captions of Figures 3 and 4, we describe the distribution and transport processes of observed avalanche deposits based on satellite images.

Bathymetry
The bathymetry of Kenai Lake can be characterized as a flat basin plain bordered by steep slopes that are typically 30 to 45 • , and locally, as steep as 70 • (Figure 1). Deltas have average slopes of 20 to 30 • . At the foot of major deltas, the hummocky mounds are probably deposits from delta failures.

Core Sedimentology
Most sediments in Kenai Lake are clastic varves (see 2. Setting; [13]). The varves have a grain-size distribution ranging from clay to silt [13] and form the matrix in which the coarse debris was deposited. The post-1964 sedimentary sequence thickens from the west to the east. The considerably thicker sedimentary sequence in the central and eastern basins (21 to 30 cm) compared to the western basin (11 to 16 cm) ( Figure 5) is a result of the river inlets, which are situated in the eastern and central part of the lake. These rivers bring sediment into the lake, which is deposited primarily near the inlets. Within the central basin, the thickness is relatively constant, with slightly thicker (max. 5 cm) sequences on delta-proximal locations.  Creek Delta (SCD) is indicated with a red star as a reference point. Two snow-avalanche deposits occurred in the western part of the central basin with little to no debris visible. Bottom: Magnified images of the observed snow-avalanche deposits and ice melting along the north shore. Melting was followed by the growth of cracks in the ice and finally cracked ice pans that floated around in the lake. Such cracks were also observed along the border of the eastern avalanche deposits.
Twenty-four debris laminae were identified in the 14 debris-containing cores; most of the debris laminae were found in cores from the southern part of the western basin and in the western part of the central basin ( Figures 5 and 6). Core KE15-07, located at the western end of the central basin, records 20 of the 24 debris laminae ( Figure 5). The equivalent spherical diameters (ESD) of the debris grains range from 0.35 mm to 12.30 mm. In general, debris consists of very coarse sand to fine gravel. An overview of the physical properties (mean, median, maximum and minimum grain size) of all the grains found in each debris lamina is presented in Table 1.   Overview of all cores from Kenai Lake where snow-avalanche debris was found. For each core there is from left to right: photograph of the core with enhanced colour histogram, frontal crosssection of the "soft tissue" reconstructed image. X-ray CT "sum-along-ray" (equivalent to inverted radiograph) image of the "bone" reconstructed image stack and a plot with the number of grains found per debris lamina for different classes of equivalent spherical diameter (ESD).

Figure 5.
Overview of all cores from Kenai Lake where snow-avalanche debris was found. For each core there is from left to right: photograph of the core with enhanced colour histogram, frontal cross-section of the "soft tissue" reconstructed image. X-ray CT "sum-along-ray" (equivalent to inverted radiograph) image of the "bone" reconstructed image stack and a plot with the number of grains found per debris lamina for different classes of equivalent spherical diameter (ESD). Debris was only sampled in cores KE15-07 and KE15-08. The debris in both cores generally consists of angular, often flattened, grains with a dark gray colour and shiny surface, sometimes accompanied by organic detritus. In debris lamina 4 from core KE15-08 (at 3.5 to 4 cm from the top of the core) one dark grey, angular pebble, 0.6 cm long with a flattened shape and a flaky surface texture was present, which was accompanied by smaller sand grains (<1 mm) (Figure 2). In debris lamina 5 (at 5 to 6 cm from the top) we found sand grains with a diameter up to 1 mm. Whereas most of the grains are angular and dark grey, some grains were more rounded and vary in colour between yellow and dark grey. Organic material, such as wood chips, was found in the matrix, together with debris particles. In debris lamina 16 at a depth of 14.5 to 15.5 cm from the top, sand grains similar to Debris was only sampled in cores KE15-07 and KE15-08. The debris in both cores generally consists of angular, often flattened, grains with a dark gray colour and shiny surface, sometimes accompanied by organic detritus. In debris lamina 4 from core KE15-08 (at 3.5 to 4 cm from the top of the core) one dark grey, angular pebble, 0.6 cm long with a flattened shape and a flaky surface texture was present, which was accompanied by smaller sand grains (<1 mm) (Figure 2). In debris lamina 5 (at 5 to 6 cm from the top) we found sand grains with a diameter up to 1 mm. Whereas most of the grains are angular and dark grey, some grains were more rounded and vary in colour between yellow and dark grey. Organic material, such as wood chips, was found in the matrix, together with debris particles. In debris lamina 16 at a depth of 14.5 to 15.5 cm from the top, sand grains similar to those in debris lamina 4 and 5 were found. In addition to the small grains, we found a 1.2-cm-long, dark grey, angular pebble.

Tephra Geochemistry
The two cryptotephras, which were identified by their magnetic susceptibility peaks, occur throughout the lake ( Figure 5). They were sampled in core KE12-05 at a depth of 10.0 to 10.4 cm (UA3196) and 18.7 to 19.0 cm (UA3197). Glass geochemical analyses are shown in Table 2. Both samples have high SiO 2 values. However, sample UA3197 is characterized by relatively high CaO and FeOt, and low K 2 O values compared to sample UA3196. All individual glass shard analyses, including standards, are available as supplementary materials (Table S1).  1 Fe is presented as total FeO (FeOt); 2 All analyses are normalized to 100% on a water-free basis; 3 n = number of analyses.

Chronology
Glass geochemical analyses show that the upper tephra (UA 3196) correlates to the 1989/90 Redoubt Volcano eruption (Figure 7). This tephra is a regional marker and present as a cryptotephra across a large portion of Yukon and Alaska [25]. The second tephra (UA 3197) was, based on its stratigraphic position, presumed to be sourced from the large 1976 eruption of Augustine Volcano with a Volcanic Explosivity Index (VEI) of 4. While we do not have access to reference material from this particular event for comparison, existing data support this correlation. Overall, the glass geochemical characteristics of this sample fall within the variability displayed by historic and pre-historic Augustine tephras, with distinctive traits, such as low K 2 O (<2.5 wt%) and higher Cl (~0.25-0.35) contents (e.g., [26,27]). The most similar material are several pumice analyses presented by Roman et al. [27] from the 1986 eruption of Augustine Volcano (Figure 7). The average of their analyses plot with, or very near, UA 3197.  [25]. Data from Roman et al. [27] of the Augustine eruption are the average of five analyses. Unfortunately, no standard deviation or standards are provided in the paper. However, the analyses for this population identified in the 1986 eruption are similar to UA 3197. All analyses, including standards, are available as supplementary materials (Table S1).
The 1976 eruption of Augustine Volcano had several phases over four months from January until April, including a large event on January 24 where column heights reached ~14,000 m a.s.l., and ashfall was reported on Kenai Peninsula and as far as Anchorage (e.g., [28,29]). The Alaska Volcano Observatory's comprehensive database of eruptions shows that only four other volcanoes had eruptions around this time: Shishalden, Pavlof, Akutan, and the Ukinrek Maars. It is very unlikely that ash from Akutan, Shishalden, and Pavlof would have reached this region; the most proximal of these three volcanoes (Pavlof) is 750 km from Kenai Peninsula, all eruptions had a VEI ≤2 with column heights not reaching above 2500 m a.s.l., and no ash was reported beyond the immediate region. Additionally, the eruptions appear to have been largely strombolian involving lava flows, suggesting that products would not have been rhyolitic as erupted at Augustine. The only other large event around this time were the eruptions of the Ukinrek Maars in 1977 (VEI 3) that ejected ash over a wide region, reaching Kodiak Island, south of Kenai Peninsula. However, these were basaltic eruptions (e.g., [30,31]).
We use the two tephra ages to revise the original varve chronology of Kenai Lake by Boes et al. [13] (Figure 8). From the top of the core until the appearance of the first tephra layer (UA3196) the original chronology is unchanged. Only in the lower part of the studied interval, from UA3196 to the 1964-turbidite, the original chronology underwent minor modification. The early 1976 Augustine provenance of tephra UA3197 allows us to attribute the varve with the tephra to varve year 1975. We Figure 7. Plots of glass geochemistry of unknowns and reference data. UA 3196 plots with all elements with the reference data for the 1989/90 eruption of Redoubt [25]. Data from Roman et al. [27] of the Augustine eruption are the average of five analyses. Unfortunately, no standard deviation or standards are provided in the paper. However, the analyses for this population identified in the 1986 eruption are similar to UA 3197. All analyses, including standards, are available as supplementary materials (Table S1).
The 1976 eruption of Augustine Volcano had several phases over four months from January until April, including a large event on January 24 where column heights reached~14,000 m a.s.l., and ashfall was reported on Kenai Peninsula and as far as Anchorage (e.g., [28,29]). The Alaska Volcano Observatory's comprehensive database of eruptions shows that only four other volcanoes had eruptions around this time: Shishalden, Pavlof, Akutan, and the Ukinrek Maars. It is very unlikely that ash from Akutan, Shishalden, and Pavlof would have reached this region; the most proximal of these three volcanoes (Pavlof) is 750 km from Kenai Peninsula, all eruptions had a VEI ≤2 with column heights not reaching above 2500 m a.s.l., and no ash was reported beyond the immediate region. Additionally, the eruptions appear to have been largely strombolian involving lava flows, suggesting that products would not have been rhyolitic as erupted at Augustine. The only other large event around this time were the eruptions of the Ukinrek Maars in 1977 (VEI 3) that ejected ash over a wide region, reaching Kodiak Island, south of Kenai Peninsula. However, these were basaltic eruptions (e.g., [30,31]).
We use the two tephra ages to revise the original varve chronology of Kenai Lake by Boes et al. [13] (Figure 8). From the top of the core until the appearance of the first tephra layer (UA3196) the original chronology is unchanged. Only in the lower part of the studied interval, from UA3196 to the 1964-turbidite, the original chronology underwent minor modification. The early 1976 Augustine provenance of tephra UA3197 allows us to attribute the varve with the tephra to varve year 1975. We also know that the position of the 1964-turbidite in the cores was triggered in March 1964 by the M w 9.2 Great Alaska earthquake [16]. As it likely takes several months, or even longer, to deposit the fine-grained cap of the turbidite [32], the first varve following the turbidite would have been deposited in varve year 1965. The original varve chronology was lacking one varve in between the two absolute time markers (i.e., the 1964 turbidite and the 1976 Augustine tephra), and we, therefore, reinterpret as two varves the multipulsed varve in which debris layer 23 is located (Figure 8).  [13]. Updated varve chronology based on two cryptotephra layers (lower blue lines). Age-depth model of core KE15-09. For each debris layer (DL; green lines) (y-axis) the age can be found on the x-axis.

Evidence for Snow Avalanches as the Origin of Coarse Debris
The mountains surrounding the lake are the source area of the debris, and the angularity of the fragments indicates that snow avalanches are likely the initial mechanism responsible for the erosion and transport of sand and gravel onto the lake ice. The metasedimentary rocks in the catchment match the rock type of the large clasts from cores KE15-07 and KE15-08. We infer that during spring melting of the lake ice, sediment transported by snow avalanches is then distributed across much of the lake by drifting ice pans.
Satellite images of the Kenai Lake area provide more evidence for snow avalanches as the primary mechanism to bring debris to the lake. Images taken at the end of the winter season (March-April) show snow-avalanche deposits as clearly visible on the frozen lake. However, they do not reach the middle of the lake, where the cores were taken. We explain the deposition of debris in the middle of the lake by the melting processes of the ice in the spring. On the 2014 satellite image ( Figure  3), ice along the north coast of the lake starts melting first, and the ice starts cracking across the lake. Cracking of the ice results in ice floes on the lake water, which carry debris initially deposited on the ice near the shore and deliver it to the centre of the lake [9]. The outlet of the lake, the Kenai River, is situated at the northwest end of the lake, and in combination with the predominant easterly wind direction in the months April and May [33], we infer that ice floes would drift westward towards the outlet. Evidence of this westward drift is found in core KE15-07, which contains the most debris laminae of all the cores and is located west of the zone with most snow-avalanche deposits ( Figure  5). The other cores with debris were all taken near to, or west of, the source areas, while no debris was found in the cores east of the source areas. Other potential transport mechanisms, such as  [13]. Updated varve chronology based on two cryptotephra layers (lower blue lines). Age-depth model of core KE15-09. For each debris layer (DL; green lines) (y-axis) the age can be found on the x-axis.

Evidence for Snow Avalanches as the Origin of Coarse Debris
The mountains surrounding the lake are the source area of the debris, and the angularity of the fragments indicates that snow avalanches are likely the initial mechanism responsible for the erosion and transport of sand and gravel onto the lake ice. The metasedimentary rocks in the catchment match the rock type of the large clasts from cores KE15-07 and KE15-08. We infer that during spring melting of the lake ice, sediment transported by snow avalanches is then distributed across much of the lake by drifting ice pans.
Satellite images of the Kenai Lake area provide more evidence for snow avalanches as the primary mechanism to bring debris to the lake. Images taken at the end of the winter season (March-April) show snow-avalanche deposits as clearly visible on the frozen lake. However, they do not reach the middle of the lake, where the cores were taken. We explain the deposition of debris in the middle of the lake by the melting processes of the ice in the spring. On the 2014 satellite image (Figure 3), ice along the north coast of the lake starts melting first, and the ice starts cracking across the lake. Cracking of the ice results in ice floes on the lake water, which carry debris initially deposited on the ice near the shore and deliver it to the centre of the lake [9]. The outlet of the lake, the Kenai River, is situated at the northwest end of the lake, and in combination with the predominant easterly wind direction in the months April and May [33], we infer that ice floes would drift westward towards the outlet. Evidence of this westward drift is found in core KE15-07, which contains the most debris laminae of all the cores and is located west of the zone with most snow-avalanche deposits ( Figure 5). The other cores with debris were all taken near to, or west of, the source areas, while no debris was found in the cores east of the source areas. Other potential transport mechanisms, such as transport by birds or floating logs, would all be from initially fluvially transported sediment, and the grains would likely be rounded. One would expect a more random distribution of the debris both in time and in space if such a process was responsible for the observed debris. Turbidity currents resulting from onshore rock failure can transport angular grains. We rule out this transport mechanism because the debris particles in the Kenai Lake sediment cores occur typically in the muddy matrix of the varved sediments and not in graded turbidites.
A similar process was proposed to explain the presence of debris in Norwegian [11], Canadian [9], and French lakes [34]. Seierstad et al. [11] observed a dirty snow avalanche, indicating that a significant amount of debris and organic material together with snow and ice was transported by avalanches onto the frozen lake. The process was first described by Luckman [9] who observed it on Canadian lakes in the 1970s. Later, Seierstad et al., Nesje et al., Vasskog et al.,and Fouinat et al. [7,10,11,34] explained the distribution of snow-avalanche debris across a lake by the same process. The snow-avalanche debris gets deposited in the spring melting season, at the transition to a new varve year. The coarse particles sink into the fine-grained top of the previous varve. This can be observed in the X-ray CT images (Figures 2 and 5), where all debris particles are situated in the black winter lamina.

Variable Avalanche Extent on Satellite Images
The satellite images from Kenai Lake in the month of April over several years show that snow avalanches were deposited on the lake almost every year. However, not every varve contains debris. In years with large snow-avalanche deposits, such as 2004, a debris deposit is present in several cores throughout the lake, whereas smaller avalanche deposits, like in the years 2012 and 2014, are not always represented by debris in one of the sediment cores ( Figure 4). The satellite images of these years show that the avalanche deposits are not reaching as far onto the lake ice as they do in other years, and it is possible that during cracking of the ice along the shore, part of the deposit already ends up near the shore [10]. As a result, these snow avalanches are not recorded in the cores taken in the middle of the lake. The 2014 satellite image ( Figure 3) shows that ice cracks form around some of the avalanche deposits along the south shore. It is likely that these deposits will not detach from the shoreline, resulting in deposition of the debris on the slopes. Furthermore, depending on which shore of a lake (i.e., shady or sunny-the northern shoreline of Kenai Lake receives more sunlight every day then the south shore, thus lake ice will melt first along the north shore) a snow avalanche is deposited, the debris may or may not get redistributed by drifting ice across the lake.

Natural Factors Affecting the Lacustrine Snow-Avalanche Records
The chance that a snow avalanche is recorded in the lacustrine sediments depends on several factors. The avalanche needs to be large enough and have a sufficiently high velocity so that it spreads far enough onto the ice when it reaches the lake. When an avalanche is too large, it could possibly break the ice when it hits the lake [34], even though this process was not observed on any of the satellite images of Kenai Lake. This process is more likely to happen in the beginning or towards the end of the freezing season as lake ice is still/already thinner and less stable. Only dirty snow avalanches, i.e., avalanches that contain organic material and coarse rock grains, will be recorded in the lake sediments as we can distinguish the coarse avalanche material from the fine lake sediments [9]. In addition to the physical characteristics of the snow avalanche, the melting processes of the lake ice also play an important role in the recording of avalanche deposits in the lake. Depending on where the cracks in the ice form, on wind direction and lake currents [9], the avalanche debris will be deposited at different places across the lake.
Apart from avalanches themselves, the physical characteristics of the lake will influence the record. For example, Lake Lauvitel, studied by Fouinat et al. [6,34] using X-ray CT, is about 130 times smaller than Kenai Lake, which strongly affects the redistribution of avalanche debris. In a smaller lake, such as Lake Lauvitel, there is an increased probability that the snow-avalanche deposits will be spread over most of the lake, thereby increasing the likelihood that a single core will contain a complete record of past avalanches reaching the frozen lake [e.g., 34]. The larger size of Kenai Lake, however, inhibits snow avalanches from reaching the central part of the lake and debris redistribution is, thus, strongly affected by the pathway of drifting ice floes during the melting season. In such large lakes, it is, thus, crucial to conduct a calibration study, such as this one, where multiple short cores are used to map the deposits ( Figure 6) and determine the optimal coring locations for a complete record (it should be noted that the cores in this study were not targeted for a snow-avalanche reconstruction). In Kenai Lake, this is the western part of the central basin, in front and to the west of the main avalanche tracks. Cores KE15-07, KE12-09, and KE15-08 contain 18, 14, and 11 out of 24 debris layers, respectively, and combined they contain all 24 recorded debris layers ( Figure 5). In the eastern part of the central basin, primarily cores close to the shores contain snow-avalanche debris (e.g., KE15-13 and KE15-15; Figure 5). Such relatively near-shore locations may, thus, contain a more complete avalanche record than cores in the center of the deep basin. Finally, while for a small lake it may be easier to reconstruct a more complete snow avalanche history of a few avalanche tracks, a large lake has the advantage to potentially record avalanches from a larger number of tracks, providing a regional avalanche record.

Comparison of Snow Avalanche and Climate Records
Climate variables, such as mean, maximum and minimum monthly temperature and monthly snowfall (data from the Alaska Climate Research Center), were compared to grain-size distributions and the number of grains from all debris laminae. Baggi et al. [35] suggested that a snowpack is weakened by 1) gradual warming of the snowpack to 0 • C, 2) loss of strength by water infiltration into the snowpack, and 3) overload of the (wet) snowpack by precipitation. These weakening factors of the snowpack could trigger a wet-snow avalanche. In Kenai Lake, we can link the presence of a debris lamina to a varve year where there was snow precipitation in a month with a mean temperature above 0 • C ( Figure 9). For example, in 1980 (varve year 1979) a high amount of snow fell in months with a mean temperature above 0 • C and in this varve year, DL17-the debris laminae with most grains of the record-was deposited.
Snow avalanches, specifically slab avalanches, can also initiate from weak bonding of snow above melt-freeze crusts. These crusts typically form during a warmer (melting) period followed by a colder (freezing) period. When fresh snow precipitates on the crust, it will fail to bond to it and can release slab avalanches [36][37][38]. Our climate data show a link between warmer temperatures in January-February and the presence of debris laminae in that varve year. For example: in 2003 (varve year 2002), when DL6 was deposited, the months of December and January had a mean monthly temperature of −3.4 and −3.2 • C, respectively. In February the mean temperature increased to 0.6 • C and decreased again in March to −4.1 • C. In addition to a drop in temperature, March showed a high amount of snow precipitation. The climate data indicate a melt-freeze crust formed in February-March and that, with increased snow in March, the snowpack became unstable and triggered slab-avalanches. A similar series of events occurred in 2004 (varve year 2003), resulting in DL5, during the year in which most avalanche deposits were observed on the satellite imagery (Figures 3 and 4). We did not identify any relationship among the mean annual winter temperature (Figure 9; yellow plot) or the total annual snowfall (Figure 9; blue plot) and the appearance of debris laminae in lake cores (Figure 9; black plot).
Based on the two mechanisms for avalanche formation, it is possible to link most of the debris laminae to the climate data. However, several other factors can play a role in the triggering of snow avalanches, such as daily temperature and snowfall, wind strength and direction, and snowpack characteristics. The climate data show several years with good snow-avalanche conditions, but during which no debris was found in the lake sediments, as well as years in which debris was found, but weather conditions did not match with snow-avalanche conditions. This may be due to insufficient climate data, coring location, ice conditions, or the redistribution process of the avalanche debris by drifting ice floes.
Quaternary 2019, 2, x 18 of 21 Figure 9. Correlation plot between climate data from Cooper Landing and Cooper Lake Project weather stations (Alaska Climate Research Center) and the snow-avalanche record derived from the sediment cores of Kenai Lake. From left to right: Total annual snowfall. Plots with monthly snowfall when (i) the mean monthly temperature is higher than 0 °C (pink) and (ii) winter temperatures are variable: the mean monthly temperature is rising from November through January, or January and/or February mean temperature is higher than December and/or March average temperature (green). In the 50-year period studied , linear fits of temperature and snow/precipitation plots show decreasing trends. The linear fit of the number of debris laminae, and, thus, the recorded snow avalanches in the lake, also shows a minor decreasing trend. This may be the result of fewer snow avalanches as the amount of snowfall and precipitation decreased during this period.

Conclusions
-X-ray CT scans of the cores from Kenai Lake show varved lake sediments and help identify coarse sediment laminae, which we infer as snow-avalanche debris, in a non-destructive methodology. We identified 24 debris laminae in 14 cores, with the coarsest and highest number of grains occurring in the western part of the lake. It is the first time that this technique was used on such a large scale for the identification of snow-avalanche debris.

-
We made minor corrections to a previously developed varve chronology using cryptotephras from the 1976 Augustine Volcano eruption and the 1989/90 Redoubt Volcano eruption. -Based on the rock type of coarse debris and satellite images of avalanche extent over 10 years, we concluded that snow avalanches emanate from the mountains around the lake, and deposit snow and sediment debris onto the frozen lake surface. In spring, the ice breaks up, and ice floes Figure 9. Correlation plot between climate data from Cooper Landing and Cooper Lake Project weather stations (Alaska Climate Research Center) and the snow-avalanche record derived from the sediment cores of Kenai Lake. From left to right: Total annual snowfall. Plots with monthly snowfall when (i) the mean monthly temperature is higher than 0 • C (pink) and (ii) winter temperatures are variable: the mean monthly temperature is rising from November through January, or January and/or February mean temperature is higher than December and/or March average temperature (green). Total number of grains present in each debris lamina over all the cores, red dots indicate that no debris was found in that year. Mean winter temperature measured from December to March. Total precipitation from November to March, red crosses indicate incomplete records. Linear fits of the different plots are indicated by dashed lines. A few examples are highlighted of years with both debris in the cores and high amounts of snowfall during warm months or when the probability of melt-freeze crusts is high due to variable temperatures.
In the 50-year period studied , linear fits of temperature and snow/precipitation plots show decreasing trends. The linear fit of the number of debris laminae, and, thus, the recorded snow avalanches in the lake, also shows a minor decreasing trend. This may be the result of fewer snow avalanches as the amount of snowfall and precipitation decreased during this period.

Conclusions
-X-ray CT scans of the cores from Kenai Lake show varved lake sediments and help identify coarse sediment laminae, which we infer as snow-avalanche debris, in a non-destructive methodology. We identified 24 debris laminae in 14 cores, with the coarsest and highest number of grains occurring in the western part of the lake. It is the first time that this technique was used on such a large scale for the identification of snow-avalanche debris.
-We made minor corrections to a previously developed varve chronology using cryptotephras from the 1976 Augustine Volcano eruption and the 1989/90 Redoubt Volcano eruption. -Based on the rock type of coarse debris and satellite images of avalanche extent over 10 years, we concluded that snow avalanches emanate from the mountains around the lake, and deposit snow and sediment debris onto the frozen lake surface. In spring, the ice breaks up, and ice floes and pans with snow-avalanche deposits drift westward on the lake (forced by river outflow and prevailing wind directions). When the ice melts, debris is deposited on the lake bottom predominantly in locations westward of where the avalanches reached the lake ice. Debris deposition in the varved sediments allows us to date snow-avalanche events with annual resolution. -Comparison of climate data with the debris laminae record shows that snow avalanches tend to occur in months with a mean positive temperature and increased snow fall. Snow avalanches also occur when there is a sudden temperature rise in winter followed by colder temperatures and snow fall. -Although the annual debris record is not complete enough to fully reconstruct the regional snow avalanche history, it can be used for decadal reconstructions of past snow-avalanche activity. -From 1966 to 2014, less snowfall has resulted in a lowering frequency of snow avalanches.