Columbia Mountain landslide: late-glacial emplacement and indications of future failure, Northwestern Montana, USA

The well preserved and undissected Columbia Mountain landslide, which is undergoing suburban development, was studied to estimate the timing and processes of emplacement. The landslide moved westward from a bedrock interfluve of the northern Swan Range in Montana, USA onto the deglaciated floor of the Flathead Valley. The landslide covers an area of about 2 km², has a toe-to-crown height of 1100 m, a total length of 3430 m, a thickness of between 3 and 75 m, and an approximate volume of 40 million m³. Deposits and landforms define three portions of the landslide; from the toe to the head they are: (i) clast-rich diamictons made up of gravel-sized angular rock fragments with arcuate transverse ridges at the surface; (ii) silty and sandy deposits resting on diamictons in an internally drained depression behind the ridges; and (iii) diamictons containing angular and subangular pebble-to block-sized clasts (some of which are glacially striated) in an area of lumpy topography between the depression and the head of the landslide. Drilling data suggest the diamictons cover block-to-slab-sized bedrock clasts that resulted from an initial stage of the failure.The landslide moved along a surface that developed at a high angle to the NE-dipping, thinly bedded metasediments of the Proterozoic Belt Supergroup. The exposed slope of the main scarp dips 30–37°W. A hypothetical initial rotational failure of the lower part of a bedrock interfluve may have transported bedrock clasts into the valley. The morphology and deposits at the surface of the landslide indicate deposition by a rock avalanche (sturzstrom) derived from a second stage of failure along the upper part of the scarp.The toe of the Columbia Mountain landslide is convex-west in planview, except where it was deflected around areas now occupied by glacial kettles on the north and south margins. Landsliding, therefore, occurred during deglaciation of the valley while ice still filled the present-day kettles. Available chronostratigraphy suggests that the ˜1-km thick glacier in the region melted before 12,000 ¹⁴C years BP—within 3000 years of the last glacial maximum. Deglaciation and hillslope failure are likely causally linked. Failure of the faceted interfluve was likely due tensile fracturing of bedrock along a bedding-normal joint set shortly after glacial retreat from the hillslope.Open surficial tension fractures and grabens in the Swan Range are limited to an area above the crown of the landslide. Movement across these features suggests that extensional flow of bedrock (sackung) is occurring in what remains of the ridge that failed in the Columbia Mountain landslide. The fractures and grabens likely were initiated during failure, but their morphologies suggest active extension across some grabens. Continued movement of bedrock above the crown may result in future mass movements from above the previous landslide scarp. Landslides sourced from bedrock above the scarp of the late-glacial Columbia Mountain landslide, which could potentially be triggered by earthquakes, are geologic hazards in the region.     

Late glacial and Holocene marine records from the Independence Fjord and Wandel Sea regions, North Greenland

The first marine sediment cores from the unexplored Independence Fjord system and the Wandel Sea, North Greenland, have been investigated to reveal the glacial marine history of the region. Two key sites in the Independence Fjord system, and an earlier analysed site from the Wandel Sea continental slope, off the mouth of Independence Fjord, are presented. The Independence Fjord sites reveal an early Holocene record (10.0–8.9 Kya) of fine-grained reddish muds with calcareous microfossils, dominated by the benthic foraminifera Cassidulina neoteretis . We suggest that a semi-permanent fast ice cover characterized the region in the early Holocene, and that the deeper troughs in the mouth region of the Independence Fjord system were intruded by subsurface Atlantic water. A stiff diamicton, at least 1.3 m thick, with coal and sandstone clasts of mainly local origin, and a 0.5-m-thick Holocene cover, are found in one of the sites. The diamicton is assumed to represent a subglacial till predating the early Holocene sediments (>10 Kya). Shallow seismic records off the mouth of Independence Fjord reveal kilometre-sized troughs with signs of glacial erosion, till deposition and a Holocene glaciomarine deposition. These features could indicate that glacial ice debouching from the Independence Fjord system at some time during the last glacial period extended to the mid-outer Wandel Sea shelf. Data from a high-resolution sediment core previously retrieved from the adjacent Wandel Sea slope indicate that the maximum ice sheet advance in this area culminated about 25–20 Kya.     

An Alpine Lacustrine Record of Early Holocene North American Monsoon Dynamics from Dry Lake, Southern California (USA)

Dry Lake (2763 m), located in the San Bernardino Mountains of southern California, USA, provides a high-resolution climate record from the coastal southwest depicting early Holocene terrestrial climate. 27 AMS ¹⁴C dates and multi-proxy analyses, including magnetic susceptibility, total organic mater, microfossil counts, and grain size, suggest the early Holocene was significantly wetter then present, due to an enhanced North American Monsoon (NAM). Elevated insolation at 9000 cal year B.P., raised summer sea surface temperatures in the Gulf of California and the eastern tropical Pacific, as well as land surface temperatures, extending the NAM into southern California. The data also provide evidence of the 8.2 ka event, which registers as a 300-year cool period characterized by reduced monsoonal precipitation, depressed basin productivity, and increased erosion. We suggest this event is the most likely period for the early to middle Holocene (9000–5000 cal year B.P.) glacial advance in the San Bernardino Mountains proposed by Owen et al. (2003, Geology 31: 729–732).     

云南千湖山第四纪冰川发育特点与环境变化

千湖山(4249 m) 是横断山脉中段保存确切第四纪冰川遗迹的山地,受西南季风影响强烈。对于研究青藏高原边缘山地冰川发育与气候和构造之间的耦合关系具有十分重要的科学意义。在千湖山海拔3500 m以上保存着古冰川侵蚀与堆积地貌,冰川发育依托海拔4000~4200 m的夷平面及其支谷地形。冰川形态类型为小型的冰帽以及由冰帽边缘溢流进入山谷的山谷冰川。应用相对地貌法,光释光(OSL) 年代测试,本文确定千湖山地区的冰进系列:末次冰盛期(LGM,22.2±1.9 ka BP)、末次冰期中期(MIS3b,37.3±3.7 ka BP、45.6±4.3 ka BP45.6±4.3 ka BP)、末次冰期早期(MIS4)。千湖山冰川前进规模是MIS3b 阶段大于末次冰盛期,主要原因是末次冰期中期(MIS3b) 时本区气候相对湿润,而在末次冰盛期(MIS2) 时气候条件比较干燥。在总体相似的气候背景下,与横断山其它存在多期次冰川作用的山地相比,千湖山只发育末次冰期的冰川作用,其差异性说明该地区冰川发育主要受山体构造抬升控制。     

哈巴雪山第四纪冰川发育特点的初步研究

哈巴雪山(5 396 m)位于横断山脉中北段,与玉龙雪山(5 596 m)以金沙江相隔,受西南季风影响强烈。在哈巴雪山3 100 m以上保存第四纪冰川侵蚀与堆积地貌,主要沿哈巴雪山西北脊两侧分布。应用相对地貌法,对比邻近山地尤其是玉龙雪山,将哈巴雪山冰期系列初步划分为:倒数第二次冰期、末次冰期早期和末次冰盛期(LGM)。以哈巴雪山哈巴河谷的冰碛物为研究对象,其末端海拔高度在倒数第二次冰期、末次冰期早期、末次冰盛期分别约为3 100 m、3 500 m、3 900 m。应用TSAM法、MEIM法、CF法计算得出哈巴雪山古雪线高度在倒数第二次冰期、末次冰期早期、末次冰盛期分别为3 675 m、4 000 m、4 200 m。     

The Southern Margin of the Late Cainozoic Ice Cap on the Central Plateau of Tasmania

The easternmost extremity of the ice cap that developed in the Tasmanian Central Highlands during the time of most extensive Late Cainozoic glaciation lay on the doleritecapped Central Plateau east and north-east of Lake St Clair. During the Last Glacial Maximum (LGM), the more restricted ice cover included a small discrete ice cap (probably less than 250-300 m thick) that formed on the Central Plateau. The LGM ice limits on the southern part of the Central Plateau, including all five southern outlet valleys, are reported here. Earlier ice limits have been identified in two of these valleys, but on the plateau proper earlier glacial deposits have been generally extensively reworked beyond the LGM limit, such that confirmation of a glacial origin for diamictons on slopes is difficult. South of the plateau, the oldest deposits flooring lower reaches of two outlet valleys indicate that ice flowed southwards directly from the plateau, but later deposits indicate diffluent flow from the Derwent Glacier.     

Snow-Push Processes in Pronival (Protalus) Rampart Formation: Geomorphological Evidence from Smørbotn, Romsdalsalpane, Southern Norway

It is demonstrated that pronival (protalus) ramparts can be formed by a snow-push mechanism and need not accumulate in the conventional manner as a result of supranival processes. Ridges in pronival positions up to 1.2 m high and of unequivocal snow-push origin are described from two sites in Smørbotn cirque, Romsdalsalpane, southern Norway. The seven lines of evidence are: (1) parallel abrasion tracks on large boulders embedded in the substrate; (2) displaced surface and embedded clasts with proximal furrows; (3) corrugated (flute-like) substrate surfaces; (4) the sickle-shaped plan-form of the ridges; (5) generally asymmetrical ridge cross-profiles (shallow, concave proximal; steep, convex distal); (6) strong preferred orientations and dips of surface-embedded clasts on ridge proximal slopes; and (7) a subnival/pronival ridge comprising loosely packed diamicton forming along the contact zone between the snowbed and substrate. This evidence indicates ridge formation by snow sliding involving bulldozing of the substrate. Factors considered important in favouring snow push producing distinct pronival ramparts at the sites include: a maritime periglacial climate with heavy winter snowfall and rapid snow-firn conversion producing snow densities of up to 900 kg m^–3; a deformable substrate with relatively small inputs of rockfall or avalanche debris; and a steep headwall susceptible to snow avalanching and hence enhanced snow supply. Consideration is given to the prospect that larger pronival ramparts can form incrementally by a snow-push mechanism.     

Reconstructing a cultural fire regime in the Pennsylvania Anthracite Region

Abstract

Fire history reconstructions provide fire managers with valuable information regarding historical fire regime dynamics. Yet, in the Central Hardwood Forest Region (CHF), there is an absence of dendroecological studies where frequent wildfires have persisted into the twenty-first century. This study presents the first documented tree-ring reconstruction of fire history in the Pennsylvania anthracite coal-producing region. An approximately 150-year fire chronology was developed from fire-scarred pitch pine (Pinus rigida) and sassafras (Sassafras albidum) collected on Spring Mountain in Schuylkill County. Principal component analysis and k-means cluster analysis were used to identify four fire regime clusters based on interacting fire and non-fire parameters. Results indicate that fires occurred under a variety of moisture conditions and that a mid-twentieth century increase in fire activity on Spring Mountain was associated with the collapse of the regional anthracite industry and possibly amplified by the 1960s drought. Fires continue to burn on Spring Mountain in the twenty-first century, contrary to the dominant pattern of fire exclusion across the CHF. Future research needs are highlighted regarding historical fire regime dynamics and cultural ties to landscape burning in the Pennsylvania Anthracite Region.     

Origin of well-rounded gravels in glacial deposits from Brøggerhalvøya, northwest Spitsbergen: potential problems caused by sediment reworking in the glacial environment

Well-rounded gravels are described from moraine-mound complexes, diamicton forefields and modern englacial thrusts at the margins of four glaciers on the northern side of Brøggerhalvøya, northwest Spitsbergen. Their shape charcteristics are compared with modern and fossil glacigenic, modern beach and Early Weichselian beach gravels from this peninsula. The best discriminators of the well-rounded gravels have been found to be the percentage-frequency roundness histograms, the roundness mid-point and roundness range diagrams and the sphericity-roundness plots. It is concluded that the gravels have been derived by englacial thrusting from Early Weichselian or last interglacial beaches in the inner parts of the fjord and in the low level cirques when sea level reached at least 50m a.s.l. and deposited the beach gravels. The discrimination between gravel in basal diamictons, proglacial outwash and modern beaches is difficult as the reworking has resulted in little particle shape change. The potential major problem caused by reworking in the glacial environment is emphasised. especially when clast shape comparisons from modern environments to older sediments are used.     

Permafrost distribution in the Posets massif,Central Pyrenees

The Posets massif is located in the Central Pyrenees and reaches a height of 3363 m a.s.l. at the Posets peak, the second highest massif in the Pyrenees. Geomorphological maps of scales 1:25000 and 1:10000, BTS (bottom temperature of winter snow), ground measurements and snow poles were used to observe the more representative periglacial active landform association, ground thermal regime, the winter snow cover evolution and basal temperatures of snow. The main active periglacial landforms and processes related to the ground thermal regime and snow cover were studied. Mountain permafrost up to 2700 m a.s.l. on northexposed slopes and up to 2900 m a.s.l. on south-exposed slopes were detected. Three permafrost belts were differentiated: sporadic permafrost between 2700 and 2800 m a.s.l. and between 2850 and 3000 m a.s.l., discontinuous permafrost between 2800 and 2950 and between 2950 and 3050, and continuous permafrost up to 2900 m a.s.l. and up to 3050 m a.s.l. on northern and southern slopes, respectively.     

GLACIAL GEOMORPHOLOGY OF THE WHITE MOUNTAINS,CALIFORNIA AND NEVADA: ESTABLISHMENT OF A GLACIAL CHRONOLOGY

The White Mountains, astride the California-Nevada stateline, are the highest and westernmost of the Great Basin ranges. This range was extensively glaciated during the Quaternary Period. Glacial landforms and scattered erratics were identified in the field using primarily morpho-and lithostratigraphic criteria, and mapped on aerial photographs and topographic maps. Topographic characteristics of the glacial deposits were analyzed using standard statistical procedures. A sequence of glacial deposits was identified in terms of six glacial stages, these glaciations named according to type site, and relative ages inferred. With the exception of perched Stage I (early) deposits along the range crest, reconstructed equilibrium-line altitudes and elevation of the glacier termini increase to the present, with glacier length and inferred size decreasing through time. Preliminary data suggest that weathering and pedogenesis are also progressive, though environmental gradients mask some of these distinctions between deposits. It is hypothesized that the White Mountains have a similar glacial chronology to the adjacent Sierra Nevada, but that with the intensification of the Sierran rainshadow during the Quaternary, the extent of glaciation in the White Mountains apparently decreased through time owing to regional tectonic uplift. This had led to the preservation of a more complete sequence of glacial deposits than in the Sierra Nevada, making the delineation of multiple mid-Quaternary events possible in many valleys. Further radiometric and chronometric dating of these deposits is in progress. [Key words: Glaciation, Quaternary, glacial geomorphology, White Mountains, California, Nevada.     

Investigation of mountain permafrost in the Kozia Dolinka valley,Tatra Mountains,Poland

The Kozia Dolinka valley lies at an altitude above 1900 m a.s.l. on the northern slope of the main ridge of the High Tatra Mountains. Mountain permafrost occurrences were studied with the use of BTS, infrared imaging, water and ground temperature measurements and DC resistivity soundings. The data suggest the existence of isolated patches of permafrost. The lowest observed bottom temperature of winter snow values was in the order of-10C. DC soundings revealed the existence of a high resistivity layer of limited extent. Permafrost seasonal monitoring was conducted with resistivity soundings. Measurements were carried out in spring-autumn 1999, when a distinct change in permafrost thickness was observed.     

The karst of Sagada Mountain Province,Luzon, Philippines

Speleological investigations during a three month period in late 1982 and early 1983 in the karst of Sagada, Mountain Province, Luzon, Philippines further clarified the initial 1980 study by Deharveng and Orousset. Over 13 km of cave passages were mapped and two main watersheds, Ambasing Spring and Balangagan Spring, were recognized.

The karst of Sagada lies between 1100 and 1600 m elevation in the Oligocene limestones of the Tineg formation. The karst is well developed, both on the surface and in the subsurface. The caves fall into five main types:

1. Potholes (Kitungan Kampus Pit, entrance to Tataya En).

2. Fossil Caves (Lomyang Cave, upper part of many caves).

3. Swallets (Latang Cave, many small swallets, especially during the wet season).

4. Underground rivers (Agoyo-Ige System. Lomyang-Latipan System, Agew Na Nayawakan System, Tataya En System, part of Balangagan Cave).

5. Springs (Ambasing and Balangagan Springs).

The caves are characterized by very wide passages, massive secondary calcite deposits, voluminous collapses, abundant volcanic debris and floods up to 600 times the low flow discharge. Study is not complete, and much work needs to be done, especially to differentiate completely the watersheds of the two main springs.     

Small valley glaciers and the effectiveness of the glacial buzzsaw in the northern Basin and Range, USA

The glacial buzzsaw hypothesis suggests that efficient erosion limits topographic elevations in extensively glaciated orogens. Studies to date have largely focussed on regions where large glaciers (tens of kilometres long) have been active. In light of recent studies emphasising the importance of lateral glacial erosion in lowering peaks and ridgelines, we examine the effectiveness of small glaciers in limiting topography under both relatively slow and rapid rock uplift conditions. Four ranges in the northern Basin and Range, Idaho, Montana, and Wyoming, USA, were chosen for this analysis. Estimates of maximum Pleistocene slip rates along normal faults bounding the Beaverhead–Bitterroot Mountains (~ 0.14 mm y^− 1), Lemhi Range (~ 0.3 mm y^− 1) and Lost River Range (~ 0.3 mm y^− 1) are an order of magnitude lower than those on the Teton Fault (~ 2 mm y^− 1). We compare the distribution of glacial erosion (estimated from cirque floor elevations and last glacial maximum (LGM) equilibrium line altitude (ELA) reconstructions) and fault slip rate with three metrics of topography in each range: the along-strike maximum elevation swath profile, hypsometry, and slope-elevation profiles. In the slowly uplifting Beaverhead–Bitterroot Mountains, and Lemhi and Lost River Ranges, trends in maximum elevation parallel ELAs, independent of variations in fault slip rate. Maximum elevations are offset ~ 500 m from LGM ELAs in the Lost River Range, Lemhi Range, and northern Beaverhead–Bitterroot Mountains, and by ~ 350 m in the southern Beaverhead–Bitterroot Mountains, where glacial extents were less. The offset between maximum topography and mean Quaternary ELAs, inferred from cirque floor elevations, is ~ 350 m in the Lost River and Lemhi Ranges, and 200–250 m in the Beaverhead–Bitterroot Mountains. Additionally, slope-elevation profiles are flattened and hypsometry profiles show a peak in surface areas close to the ELA in the Lemhi Range and Beaverhead–Bitterroot Mountains, suggesting that small glaciers efficiently limit topography. The situation in the Lost River Range is less clear as a glacial signature is not apparent in either slope-elevation profiles or the hypsometry. In the rapidly uplifting Teton Range, the distribution of ELAs appears superficially to correspond to maximum topography, hypsometry, and slope-elevations profiles, with regression lines on maximum elevations offset by ~ 700 and ~ 350 m from the LGM and mean Quaternary ELA respectively. However, Grand Teton and Mt. Moran represent high-elevation “Teflon Peaks” that appear impervious to glacial erosion, formed in the hard crystalline bedrock at the core of the range. Glacier size and drainage density, rock uplift rate, and bedrock lithology are all important considerations when assessing the ability of glaciers to limit mountain range topography. In the northern Basin and Range, it is only under exceptional circumstances in the Teton Range that small glaciers appear to be incapable of imposing a fully efficient glacial buzzsaw, emphasising that high peaks represent an important caveat to the glacial buzzsaw hypothesis.     

长白山现代理论雪线和古雪线高度

According to the glacial landforms and deposits with the optically stimulated luminescence (OSL) dating results, two glacial stages of the last glacial cycle (LGC) and Late Glacial were identified. The Late Glacial stage (Meteorological Station glacier advance) took place about 11 ka (11.3±1.2 ka), and the last glacial maximum (LGM), named Black Wind Mouth glacier advance, occurred at 20 ka (20.0±2.1 ka). Based on the Ohmura’s formula in which there is a relationship between summer (JJA) atmospheric temperature (T) and the annual precipitation (P) at ELA, the present theoretical equilibrium line altitude (ELA_t) in Changbai Mountains was 3380±100 m. Six methods of accumulation-area ratio (AAR), maximum elevation of lateral moraines (MELM), toe-to headwall altitude ratios (THAR), the terminal to summit altitudinal (TSAM), the altitude of cirque floor (CF), and the terminal to average elevation of the catchment area (Hofer) were used for calculation of the former ELAs in different stages. These methods provided the ELA for a range of 2250–2383 m with an average value of 2320±20 m during the LGM, which is 200 m higher than the value of previous investigation. The snowlines during the Late Glacial are 2490 m on northern slope, and 2440 m on western slope. The results show that the snowline on northern slope is 50 m higher than that on western slope during the Late Glacial, and the average snowline is 2465m. The ΔELA values were more than 1000 m during the LGM, and about 920 m lower than now during the Late Glacial stage respectively. Compared with Taiwanese and Japanese mountains in East Asia during the LGM, the effect of the uplift on ELA in Changbai Mountains during the glaciations (i.e. 20 m uplift in the LGM and 11 m in the Late Glacial) is not obvious. Foundation: National Natural Science Foundation of China, No.40571016 Author: Zhang Wei (1969–), Ph.D and Professor, specialized in Quaternary environment and climate geomorphology.     

Factors controlling carbonate dissolution rates quantified in a field test in the Austrian alps

Carbonate dissolution rates were investigated by measuring the mass difference of carbonate tablets exposed to natural dissolution for 1 year. 70 tablets were distributed over 13 test sites on the north slope of the Hochschwab Massif in the Austrian Alps. The influences of altitude, subsoil vs. sub-aerial exposure, vegetation, karst morphology, soil humidity, sample lithology, and sample surface morphology were investigated. The observed dissolution rates varied between 13 and 40 μm/a for subcutaneous samples and about 11 μm/a for sub-aerial exposure. Outstandingly high rates of 48 μm/a were observed in a doline and nearly zero rates were measured at a site influenced by seeping spring water.A mass balance, using high-resolution hydrological data, was calculated for the Kläffer Spring, which has an average outflow of 4.8 m³/s. It indicated a loss of 21×10⁶ kg of carbonate rock per year which gives an average dissolution rate of 95 μm/a for the catchment area of 83 km². The dissolution rates of 10 μm/a from the sub-aerially exposed samples are comparable to values from limestone pedestals, which were protected from dissolution by glacial erratics for the past 15 ka.     

质疑粤西北的冰川遗迹

从地貌、沉积物、第四纪地层关系、古气候、古环境等方面论证了广东封开-怀集-带所谓的冰川谷、冰蚀三角面、悬谷、角峰、刃脊、冰半、鼻山尾、冰碛砾石、蛇形丘以及近百万年来发生过3次冰期，是既不符合地学基础理论也不符合野外实际，所以是不确实的，那里没有第四纪冰川遗迹。     

Landscape evolution and glaciation of the Rwenzori Mountains, Uganda: Insights from numerical modeling

The Rwenzori Mountains are a high alpine mountain chain, about 40 × 80 km in size, just north of the equator in the western branch of the East African Rift System in Africa. The central part of the mountain chain is located in Uganda, and the highest peak, the Margherita Peak with 5119 m, lies on the border to the Democratic Republic of Congo. Topography is very pronounced, with steeply incised valleys and clear glacial landforms in the upper part of the mountain chain. The Rwenzori Mountains are an unusually high mountain chain located in the extensional setting of the East African Rift System, and the large elevation poses a challenging problem for geodynamists to explain.We have used the landscape evolution model ULTIMA THULE, which combines hillslope diffusion, fluvial erosion, and glacial abrasion and is driven by a climate driver, simulating the variations in temperature, precipitation, and relief over several glacial cycles. With a simulation time of 800 ka, we test the hypothesis of climate-tectonic interactions on the uplift of the Rwenzori Mountains.Our results show that a moderate cooling of around 6° causes widespread glaciation of the high mountain regions as observed during the peak glacial phases, and that morphological processes degrading the landscape allow for a tectonic uplift rate of around 0.5 mm a^− 1.     

Late summer glacial meltwater contributions to Bull Lake Creek stream flow and water quality,Wind River Range,Wyoming, USA

ABSTRACT

The Wind River Range in Wyoming contains more glacial ice than any other location within the USA’s Rocky Mountain states of Colorado, Idaho, Montana, and Wyoming. Bull Lake Creek watershed in the southeast portion of the range contains five major (0.6–1.5 km²) glaciers along with numerous smaller glaciers that contribute to the Wind River. Field measurements were made of discharge from the Knife Point and Bull Lake Glaciers to determine the contribution of glacial meltwater to the river system. Water samples were collected and analyzed for stable isotopes, major ions, nutrients, and selected trace elements. Meltwater from the two glaciers contributed 13.9% to Bull Lake Creek streamflow (site BL-3), with all glaciers within the Bull Lake Creek watershed estimated to be contributing 55.6% to the streamflow of Bull Lake Creek (United States Geological Survey gage) during the August 2015 study period. Hydrogen and oxygen stable isotope analysis indicated as much as 80% of late summer discharge in the upper Bull Lake Creek watershed was attributed to glacial meltwater. This study also found that nutrients (NO₃ – NO_2, total P) from glacial meltwater can be a significant source of nutrient loading to Bull Lake Creek.