UBC Theses and Dissertations

UBC Theses Logo

UBC Theses and Dissertations

Hydrologic control on proglacial suspended sediment dynamics Richards, George 2001

You don't seem to have a PDF reader installed, try download the pdf

Item Metadata


ubc_2001-0503.pdf [ 5.97MB ]
JSON: 1.0090147.json
JSON-LD: 1.0090147+ld.json
RDF/XML (Pretty): 1.0090147.xml
RDF/JSON: 1.0090147+rdf.json
Turtle: 1.0090147+rdf-turtle.txt
N-Triples: 1.0090147+rdf-ntriples.txt
Original Record: 1.0090147 +original-record.json
Full Text

Full Text

HYDROLOGIC CONTROL ON PROGLACIAL SUSPENDED SEDIMENT DYNAMICS By GEORGE RICHARDS B.A. (Hons), U n i v e r s i t y of.Oxford,  1999  THESIS SUBMITTED IN PARTIAL FULFILLMENT THE REQUIREMENTS FOR THE DEGREE OF MASTER OF SCIENCE in THE  FACULTY OF GRADUATE STUDIES (Department of Geography)  We accept t h i s t h e s i s as conforming to the r e q u i r e d standard  THE UNIVERSITY OF BRITISH COLUMBIA September  2001  © George Richards,  2001  In  presenting this  degree at the  thesis  University  in  partial  fulfilment  of  British Columbia,  of  the  requirements  I agree that the  freely available for reference and study. I further agree that copying  of  department  this thesis for scholarly or  publication of  by  his  her  or  representatives.  fe°6&AP 7 K  The University of British Columbia Vancouver, Canada  Date K  DE-6 (2/88)  It  is  an  fcP<r£r^&&L  f  "2001  advanced  Library shall make it  permission for extensive  granted by the understood  head of  that  this thesis for financial gain shall not be allowed without  permission.  Department of  purposes may be  for  copying  my or  my written  Abstract Measuring suspended sediment concentration (SSC) is both costly and labour intensive. Temporal records of S S C are, however, of paramount importance in elucidating issues relating to geomorphology, ecology and water quality.  Rating curves, that relate S S C and  discharge by a simple linear regression function, are frequently employed by workers to address the problems of recording S S C . Such functions, however, rarely account for more than 50% of the variability in observed S S C . The aim of this thesis is to formulate subseasonal  predictive S S C models  suspended  sediment  and to investigate hydrologic controls on proglacial  dynamics using data collected  from a glacier-fed  stream, Coast  Mountains, British Columbia. In order to model proglacial S S C , the hydrologic season was initially divided into nival, nival-glacial, glacial and autumn recession periods, according to sudden shifts in the ratio of stream discharge between the glacierised and a neighbouring unglacierised catchment of similar size and aspect.  Multiple regression functions, to predict S S C , were then  developed for each period. These regression models incorporate a suite of easily measured variables and are shown heteroskedasticity  to reduce  significantly,  initial problems  of autocorrelation,  and non-linearity of the SSC-discharge relationship.  Analysis of the  significant parameters in the multiple regression models, the hysteretic relationship between S S C and discharge, and downstream changes in S S C reveal that short-term, within channel, storage of fine sediment may be an important control on proglacial suspended sediment dynamics in this complex glaciofluvial lacustrine system.  ii  Table of Contents Abstract  ii  Table of Contents  iii  List of Tables  vi  List of Figures  vii  Acknowledgements  ix  Chapter 1 Introduction  1  1.1  Alpine hydrology and suspended sediment transport  2  1.2  Modelling proglacial S S C  6  1.3  Objectives  8  1.4  Thesis outline  9  Chapter 2 Study Area and Methods 2.1  10  Regional context  10  2.1.1 Climate  10  2.1.2 Geology  10  2.1.3 Vegetation  10  2.1.4 Glacial influences on hydrology 2.2  Study basins  11 11  2.2.1 Place Creek  11  2.2.2 Eight Mile Creek  18  2.3  Streamflow measurement  20  2.3.1 Stage measurement  20  2.3.2 Discharge measurement  20  2.3.3 Rating curves  21  2.4  Suspended sediment concentration  23  2.5  Water chemistry  24  2.5.1 Electrical conductivity  24  iii  25  2.5.2 Cations and anions 2.6  Meteorological variables  26  Chapter 3 Results 3.1  27  Overview of study period  27  3.2 Data quality  27  3.2.1 Discharge measurement  27  3.2.2 SSC measurement  30  3.3 Streamflow and S S C variations  32  3.3.1 Place Creek discharge and SSC time series  32  3.3.2 Eight Mile Creek discharge and SSC time series  34  3.3.3 Standardised discharge series 3.4  Seasonal hydrograph division  37 37  3.4.1 CUSUMplot of Q ratio  37  3.4.2 Relations between electrical conductivity and discharge 3.5  S S C - Q relations for Place Creek  39 41  3.5.1 Seasonal SSC-Q relation  41  3.5.2 Sub-season SSC-Q relations  46  3.5.3 Multiple regression models  50  3.5.4 SSC-Q hysteresis relations  55  3.6 Upstream-downstream comparisons of S S C  61  3.7  63  Chapter summary  Chapter 4 Discussion 4.1  65  Hydrologic control on S S C  65  4.1.1 Seasonal scale  65  4.1.2 Synoptic scale  66  4.1.3 Diurnal scale - glacial period  68  4.1.4 Integrated temporal scales  69  4.2 Modelling S S C  72  Chapter 5 Conclusions  77  5.1 Summary of key findings  77  iv  5.1.1 SSC-Q relations in a glacier-fed stream  77  5.1.2 Multiple regression models  77  5.1.3 Within-channel transient sediment storage 5.2 Future research  78 78  References  80  Appendix A  85  v  List of Tables Table 3.1: Results of ISCO - DH48 Paired Sample T-test  30  Table 3.2: Descriptive Statistics for Place Creek Sub-seasons  39  Table 3.3: Summary of OLS SSC-Q Seasonal Regression Model (values in brackets are the p-values for bo and bi)  43  Table 3.4: Summary of SSC-Q OLS Regression Models (values in brackets are the pvalues for bo and bi)  47  Table 3.5: Results of Sub-season MLR Analysis  53  Table 3.6: Summary of MLR Assumption Test Statistics (PRESS is the Prediction Sum of Squares)  54  vi  List of Figures Figure 2.1:  Location  map of study area (provided  courtesy of E . Leinberger,  Department of Geography, University of British Columbia)  12  Figure 2.2: M a p showing drainage areas (based on data from 1987) for Place and Eight M i l e Creeks  13  Figure 2.3: Vegetation in drainage basins for Place and Eight M i l e Creeks (based on 2001 data from Ministry of Forests, Government of British Columbia)  14  Figure 2.4: Geology of Place and Eight M i l e basins (based on data courtesy of D r Journey, Geological Survey of Canada)  15  Figure 2.5: Place Lake. A braided stream flows from a small lake at the terminus of Place Glacier into the Place Lake on the right side. Place Lake outlet is on the far side  16  Figure 2.6: Place Creek in the valley bottom during spring snowmelt  17  Figure 2.7: Eight M i l e Creek at low flow  19  Figure 2.8: (a) Place Creek and (b) Eight M i l e Creek stage-discharge rating relations. Eight M i l e has two plots (old and new) because of a shift in the original rating curve resulting from movement of the stilling well during high flows in early July 2000  22  Figure 3.1: Normal (1984 to 2000) versus 1998,  1999 and 2000 (a) monthly snow  course, (b) summer monthly air temperature ( T ) and (c) summer monthly A  precipitation (P).  Data provided by (a) B . C . Environment Historical Snow  Survey Data and ((b) and (c)) D . Hutchinson, Hydrology Applications Group, Environment Canada  28  Figure 3.2: Manual stage-mV scatter plots for (a) Place Creek and (b) Eight M i l e Creek  :  29  Figure 3.3: (a) I S C O versus D H 4 8 values of S S C for all samples and (b) left versus right bank D H 4 8 values of S S C  31  Figure 3.4: Time series of discharge, S S C , air temperature and precipitation for Place Creek (2000 season)  33  Figure 3.5: Time series of discharge, S S C , air temperature and precipitation for Eight M i l e Creek (2000 season)  35  vii  Figure 3.6: Standardised discharge time series for (a) Place Creek and (b) Eight M i l e Creek.  Arrows on the (c) C U S U M S plot indicate sub-seasonal transition  dates  38  Figure 3.7: Total cationic concentration-EC relation for Place Creek  40  Figure 3.8: E C - Q relations for Place Creek. Data points are grouped by sub-seasonal division dates defined in Figure 3.6c  42  Figure 3.9: (a) simple S S C - Q rating relation; (b) S S C - Q relation residual plot; (c) seasonal and diurnal scale autocorrelation in the S S C - Q relation; and (d) log transformed S S C - Q relation (note the lack of fit for low values of S S C )  44  Figure 3.10: (a) Seasonal versus (b) sub-seasonal S S C - Q relation residual plots  48  Figure 3.11: Sub-seasonal S S C - Q rating relations (note that S S C scale has been reduced so that the grouping of points can be more readily distinguished)  49  Figure 3.12: Examples of S S C - Q hysteresis for: ((a) and (b)) storm runoff events; (c) a snowmelt event; and ((d) and (e)) more complex S S C - Q relations.  Arrows  on S S C - Q plots indicate direction of hysteresis (filled and open arrows divide discharge cycles/events).  Time plots (right) correspond to hysteresis loops  56  Figure 3.13: Diurnal and synoptic scale S S C - Q hysteresis during the main.ablation period (arrows indicate direction of synoptic scale hysteresis).  A l l individual  loops have a clockwise direction  59  Figure 3.14: Time series of organic S S C as a percentage of total S S C for four different discharge events in Place Creek  60  Figure 3.15: Upstream-downstream comparison of lake outlet S S C , and valley discharge and S S C for (a) early August and (b) late August 2000 Figure 4.1: Flow chart model for prediction of S S C from sub-seasonal M L R models  viii  62 73  Acknowledgements I am grateful for the support and encouragement given to me by Dan Moore, who supervised this work throughout.  Olav Slaymaker and Michael Church made valuable suggestions  during the formulation of this research. I am grateful to Brian Menounos for providing invaluable advice on field work and data analysis and for helping with map making. determine their chemistry.  Giles Brown analysed water samples to  Ken Hall provided I S C O water samplers.  Geoff H i l l was a  friendly and dedicated field assistant, who tolerated a season of filtering and long nights of hourly sampling at Place Glacier. Matt Gellis, Dion Whyte and Carolyne Albert also helped in the field.  Catherine Griffiths gave several word processing tutorials.  Vincent Kujala  assisted with Fortran programming. I am also grateful to Jo for providing accommodation at the Birkenhead Resort. This work was completed while the author was recipient of a B U N A C Educational Trust Scholarship and the John S. Ball Award. Funding for field work was provided by Dan Moore ( N S E R C grant). I owe great thanks to Elaine Cho and the office staff, who made my time at U . B . C . much easier, and to my family and friends for their continual support.  ix  Chapter 1 Introduction Temporal variations in suspended sediment transport have received considerable attention. Such studies have traditionally been spurred by geomorphic interests in, for example, global sediment yield to oceanic basins (Milliman and Syvitsky, 1992), estimates of landscape denudation rates (Collins, 1998), models of sediment transport (Syvitski et al, 2000; Picouet et al. 2001) and relative rates of glacial and non-glacial erosion (Hicks et al., 1990; Harbor and Warburton, 1993; Hallet et al. 1996). Suspended sediment is, however, also important from an ecological and water quality perspective because it can alter water chemistry and cause temperature decreases and turbidity increases (Ryan, 1991). In particular, Richardson et al. (2001) stressed that high turbidity levels decrease upstream fish migration and thus species' populations.  Most studies of alpine proglacial suspended sediment dynamics have concentrated on glacial erosion processes inferred from measurements of suspended sediment concentration (SSC), discharge (Q) and electrical conductivity (EC) made within close proximity of the glacier terminus. From initial studies of proglacial SSC in the 1950s, when "almost nothing was known of the seasonal variation in sediment concentration of glacial streams" (Matthews, 1956), the literature has grown extensively, incorporating North American, European, Asian and Arctic based studies (e.g. Collins, 1979; Hammer and Smith, 1983; Gurnell and Fenn, 1984; Gurnell et al, 1994; Willis et al, 1996; Collins, 1998). Relatively few investigations, however, have focused directly on downstream transport and storage of glacier-derived sediment. Intuitively, it would seem that the next step is to conduct a study that provides an integral understanding of the temporal and spatial link between glaciers and 1  downstream hydrology and suspended sediment transport.  This thesis is an attempt to  achieve this link by presenting results of an investigation of temporal and spatial changes in the hydrology and suspended sediment transport of a retreating cirque glacier system in the Coast Mountains, British Columbia.  Under conditions of sustained glacier retreat, the  delivery of water, sediment and dissolved substances  by proglacial stream systems to  downstream reaches will undoubtedly change, and will likely have ecological impacts. Thus the study of these processes is of great environmental and academic importance.  Sections 1.1 and 1.2 provide a focused literature review and identify knowledge gaps that this thesis will aim to fill through the objectives stated in Section 1.3.  Section 1.1  concentrates on the delivery of water and sediment from alpine headwaters to downstream reaches.  Section 1.2 focuses on attempts to model proglacial glacier-derived water and  suspended sediment dynamics.  1.1 Alpine hydrology and suspended sediment transport A n appreciation of hydrology is required to understand suspended sediment transport, because sediment input to a stream will occur where hydrologic contributing areas and sediment source areas overlap. In addition, the mode of runoff generation can control the erosion process. For example, the kinetic energy of rainfall on snow-free ground can detach soil particles, rendering them vulnerable to transport, while the presence of snow cover protects the soil surface.  Alpine seasonal hydrographs are typically characterised by spring snowmelt and, where basins are glacierised, late spring and summer glacier melt sources  2  of water.  Throughout the season, flows may also be punctuated by short-term rainfall-driven storm runoff events. Net radiation and sensible heat transfer provide the dominant mechanisms for snow and ice melt (Hannah and McGregor, 1997). During the nival season, meltwater slowly percolates to the ground through the porous snowpack and then travels as overland or subsurface flow to the stream channel. Discharge typically peaks late in the evening or early in the morning because of routing delays through the snowpack. The rising limb of a diurnal snowmelt hydrograph is thus controlled by melt rates and the falling limb by storage characteristics of the basin (Singh et al., 2000).  In mountain catchments, hydrologic  contributing areas tend to be dominated by snowmelt in the lower elevation zones at the beginning of the spring. Indeed, the higher zones may experience snowfall while lower zones experience rain and melt. As air temperatures and energy supply increase through the season, the zone of active snowmelt progresses to higher elevation zones.  The lower  boundary of the active contributing area is dictated by the snowline elevation, which also tends to progress upslope.  Thus, the active hydrologic contributing area progresses  systematically upslope through the nival melt season.  Glaciers tend to be located in the  higher elevation zones, and thus become hydrologically active in mid- to late summer. Once most of the non-glacial snowcover has melted, streamflow would be dominated by glacier melt contributions. During the glacial melt season, the daily pattern of discharge is driven by diurnal variations in energy input. Discharge consists of baseflow (subglacial meltwater, water stored in cavities, meltwater percolating through snow and groundwater (Rbthlisberger and Lang, 1987)), superimposed by a daily cycle consisting of a rapid subglacial meltwater routing component and a supraglacial component (Benn and Evans, 1998).  3  The hydrologic effect of rain events may depend on storm characteristics (particularly air temperatures). During cool storms, precipitation may fall as snow at higher elevations, so that only lower elevation bands may generate rainfall runoff. The quantity of rainfall runoff generated will also depend on soil moisture conditions. Where soils are wet (e.g. areas that are experiencing active snowmelt or that have only recently become snow-free), generation of rainfall runoff may be relatively efficient. However, if an area has been snow-free for a period of warm, dry weather, the resulting soil moisture deficit may effectively limit runoff generation. elevation.  Added to these factors is the fact that precipitation tends to increase with Thus, the hydrologically active contributing area for rain events may vary in  complex ways from storm to storm, as opposed to the seasonally progressive nature of contributing areas for snowmelt.  During the ablation season, fine sediment evacuated from the glacier subsole is transported as suspended load in proglacial melt streams.  Several studies have also  highlighted the significance of the bed (0strem, 1975; Hammer and Smith, 1983; Gurnell et al., 1999) and dissolved (Barsch and Caine, 1984) components of sediment load in high mountain rivers, but these components of sediment load were not measured in this thesis and the focus will be on suspended sediment.  Temporal variability of proglacial SSC is  frequently claimed to result from the interaction between discharge fluctuations and accumulation of sediment at the glacier base, as the subglacial drainage system develops through the season (Collins, 1998). Proglacial lakes are, however, also important because they act both as sediment traps and buffers to daily SSC fluctuations. Weirich (1985), for example, formulated a sediment budget for a high-energy glacial lake in southeastern British Columbia. The study stressed that SSC was much lower and varied less in the outlet than  inlet stream. Several anomalous spikes of S S C were, however, recorded and concluded to be the product of overflow events.  The spatial variability of S S C beyond glacier margins has, however, not been studied greatly.  Gurnell (1982) investigated downstream changes in S S C in order to highlight the  importance of sampling location, but the study focused on tributary sediment sources.  More  recently, Lenzi and Marchi (2000) argued that, during flood events in a steep alpine stream, within-channel suspended sediment sources provided a significant proportion of sediment loads during storms. N o data were provided to substantiate this claim and the study was only located at one downstream sampling site.  There is a need, therefore, for more intense  seasonal monitoring programmes at both up and downstream S S C sampling locations.  Attempts to link alpine runoff sources to suspended sediment transport in mountain streams have been made in the literature.  In a Norwegian study comparing suspended  sediment transport during glacial melt and rainfall-dominated runoff, Richards (1984) concluded by proposing that different hydrologic source areas determined the dynamics of sediment quantities and particle sizes. Chikita (1993) examined seasonal changes in S S C in a snowmelt fed river in western Hokkaido, Japan. In that study it was discovered that as the snowmelt season progressed, S S C decreased for a given value of discharge because of an upslope decrease in the supply of bank sediment.  Finally, Lenzi and Marchi (2000)  suggested that for unglaciated alpine basins, suspended sediment transport takes place at low discharge during snowmelt because erosion processes consist of the removal of loose, finegrained material from slopes by surface runoff. Except for large floods when bank failure may contribute much sediment to the channel, low discharges during summer and autumn rainfall storms were seldom found to transport suspended sediment because of snowmelt  5  depletion of sediment loosened during the previous winter.  If diurnal through to seasonal  scale downstream S S C dynamics are to be understood, further investigations of the link between suspended sediment transport and hydrologic source area must be made.  1.2 Modelling proglacial SSC High values of S S C can occur under the following conditions: (1) at the start of the ablation season, when there is an extensive supply of glacially eroded basal sediment that has accumulated over the winter; (2) during large melt events when the capacity of the subglacial drainage system (Richards, 1984).  to transport sediment  increases;  and (3)  during rain-induced events  Periods of relatively low S S C are caused by exhaustion of suspended  sediment supply on diurnal (Collins, 1979), synoptic (Gurnell, 1987) and seasonal timescales (e.g. 0strem, 1975). In order to predict S S C these factors must be accounted for.  Stream S S C is both a labour intensive and time-consuming variable to measure.  In  the absence of a continuous (hourly) sampling programme, the simplest and most widespread procedure adopted is to construct  a suspended  sediment  rating curve derived from  observations of S S C and Q . A rating curve consists of a graph relating S S C and Q . From the equation of the resulting curve, which is typically a power function, S S C can be calculated from a continuous record of Q.  In this way rating curves can be used to interpolate and  extrapolate measured S S C time series.  Applicability and accuracy of suspended sediment rating curves are, however, highly questionable.  Although errors may stem from inaccurate field and laboratory procedures,  much of the scatter in the S S C - Q rating plot can be explained by: (1) seasonal variations in sediment availability, associated with early season flushing and late season exhaustion of  sediment supply (0strem, 1975), (2) hysteresis effects related to variations in sediment supply on the rising and falling limbs of diurnal flow cycles (Collins, 1979, Bogen, 1980), and (3) sediment supply dynamics associated with rainfall induced flow events (Richards, 1984, Stott and Grove, 2001).  Various methods used to overcome problems with suspended sediment rating curves have  been  adopted  extensively  in  the  literature.  Principally, these  problems  are  autocorrelation in the residual series and the non-normal distribution of error terms. Gurnell and  Fenn  (1984)  and Fenn  (1989)  accounted  for  the  former by  using  univariate  AutoRegressive Integrated Moving Average ( A R I M A ) stochastic time series models of Q and S S C .  Although predictions based on A R I M A models were superior to the original rating  functions, A R I M A modelling has not been adopted extensively in the literature because a detailed series of uniformly spaced observations is needed.  In order to normalise the error  distribution, logarithmic transformation of the Q and S S C series is common.  Ferguson  (1986, 1987), however, stressed that such transformations, when back-transformed to a rating curve, lead to underestimation of S S C , by up to 10 to 50% (Asselman, 2000), because the residual term is additive in a power function non-linear regression model and multiplicative in a log-transformed model. One of the properties of a regression function is that it must pass through the arithmetic and geometric means for non-linear and log-transformed models, respectively.  Because geometric means are systematically lower than arithmetic means, the  log-transformed equation underestimates true values of S S C because of the bias afforded to lower values of S S C (Walling, 1977; Jansson, 1985).  Miller (1984) proposed a bias  correction factor to reduce the transformation bias that arises when regressions of natural logarithms are detransformed.  7  Separate rating curves can also be fitted to different ablation seasons, for specific subperiods within ablation seasons and for conditions of rising and falling discharge only (e.g. Fenn et al., 1985; Gurnell et al., 1992a).  This approach allows for changing sediment  availability and supply dynamics to be partly accounted for. Methods used to sub-divide seasonal discharge hydrographs are, however, subjective and have been a result of 'hydrological judgement' (Hannah et ah, 2000). Hodgkins (1996), for example, described his method of division as "arbitary" and provided no quantitative justification. In another case, Hodson et al. (1998) based their sub-seasonal divisions on changes in the amount of energy delivered to the basin and on the position of the transient snowline. Hannah et al. (2000) attempted a more statistical approach to the classification of the shape and magnitude of diurnal hydrographs using a combination of principal components analysis and cluster analysis.  Nine different classes of hydrograph classes were developed using these  techniques, but the division of discharge series into their relevant class was subjective.  Recently, workers have developed multivariate rating equations that incorporate additional explanatory variables such as the rate of change of discharge, cumulative discharge and hourly precipitation (e.g. Hodgkins, 1999; Hodson and Ferguson, 1999). Generally, though, such complex rating equations contain variables derived from measured SSC (e.g. lagged values of SSC) and so independent model validation is uncommon because of the time and cost involved in measuring the predictor variables.  1.3 Objectives Drawing from the knowledge gaps identified in the preceding literature review, the objectives of this thesis are: (1) to investigate whether and how the SSC-Q relation varies with changing hydrologic conditions in a glacier-fed mountain catchment, (2) to evaluate the 8  use of multiple regression models for predicting S S C , and (3) to investigate the significance of channel storage on the downstream delivery of suspended sediment. Data from a glacierfed stream in the Coast Mountains, British Columbia, Canada, will be used to fulfil these objectives. A neighbouring creek draining an unglacierised basin will be used as a reference basin with which to aid interpretation of discharge and S S C series in the glacierised basin.  1.4 Thesis outline Details of the study area and methods are presented in Chapter 2. In Chapter 3 the results of the investigation are presented.  Firstly, discharge and S S C time series for both basins are  analysed over the 2000 season and then attention is drawn to the glacierised basin.  The  seasonal hydrograph is separated and for each sub-season, multiple regression models are developed to predict S S C . Hysteresis in the S S C - Q relation is examined on diurnal, synoptic and seasonal timescales. Finally an upstream-downstream comparison of S S C is undertaken for two diurnal data series.  Chapter 4 discusses the hydrologic control on S S C within Place  Creek on diurnal, synoptic and seasonal scales, and addresses the importance of adopting an integrated approach to such timescales.  Multiple regression models for predicting S S C are  then assessed. Chapter 5 summarises the key findings of the thesis and makes suggestions for future research.  9  Chapter 2 Study Area and Methods 2.1 Regional context 2.1.1 Climate The winter climate of coastal British Columbia is dominated by successions of midlatitude cyclones which move onshore from the Pacific Ocean. Associated frontal activity provides over 60% of annual precipitation from November-March in the south coastal region (Moore and McKendry, 1996).  Mean annual precipitation averages more than 3500 m m at high  elevations within the Coast Mountains, but in the valleys ranges from 1500 to 2500 mm. The study basins in this investigation are, however, located within a transitional zone from a strong maritime to dry interior influence, where annual precipitation ranges from 750 to 1000 m m at valley sites. In the study area, mean monthly air temperature ranges from 0 to - 5 ° C in January and from 18 to 20 ° C in July (Farley, 1979).  2.1.2 Geology The Coast Mountains lie along a northwest band of Cretaceous and Tertiary plutonic and metamorphic geology. The topography is dominated by glacial landforms. The study basins lie on granodiorite and dacite-pyroclastic rocks, and so are typical of the regional geology.  2.1.3 Vegetation From the coast to the study area, the dominant vegetation changes from Coastal Douglas Fir to Coastal Western Hemlock in the valleys and from Subalpine Mountain Hemlock to Subalpine Engelmann Spruce at higher elevations. Within the Coast Mountains the tree-line ranges from 1500 to 2000 m (Ryder and Thomson, 1986).  10  2.1.4 Glacial influences on hydrology The majority of glaciers within the Coast Mountains reached their maximum Holocene extent during the 1800s and have since experienced overall recession punctuated by short periods of advance (Ryder and Thomson, 1986).  Few studies have examined the contribution of  glaciers to streamflow in the Coast Mountains. Moore (1993), however, used a conceptual model to estimate that for the Lillooet River in the southern Coast Mountains, over 50% of total basin runoff in August is derived from glacial melt.  The glacial cover of the Lillooet  River basin is approximately 17%.  2.2 Study basins Figures 2.2 to 2.4 illustrate the topography, vegetation and geology of both study basins.  2.2.1 Place Creek Place Creek (PC) is located approximately 140 km north of Vancouver in the southern Coast Mountains, British Columbia, Canada (Figure 2.1). Its basin faces northwest and has an area of 13 k m (Figure 2.2). In 1965 Place Glacier had an area of 3.98 k m , when a mass balance 2  2  study was instigated as part of the International Hydrological Decade.  In 1987 the glacier  was 3.4 k m . Evidence from air photographs indicates that from 1982 to present, a series of 2  three proglacial lakes has formed within 200 m of the glacier terminus.  Place Creek  originates at the outlet of the largest of these lakes (Figure 2.5), Place Lake, and flows along a steep channel for 4 km to the valley below.  Large boulders dominate the bed morphology  of Place Creek, except over a short braided section about halfway down the creek where the channel gradient declines. Channel width ranges from approximately 4 to 5 m (Figure 2.6) at the gauging site.  11  Figure 2.1: Location map of study area (provided courtesy of E. Leinberger, Department of Geography, University of British Columbia).  12  Kilometres  N  A  Lake 100 m contour Watershed Ice River  Figure 2.2: M a p showing drainage areas (based on data from 1987) for Place and Eight Mile Creeks.  13  Kilometres  4  Vegetation H H Alpine Tundra and Ice I  N  I Coastal Western Hemlock  lllllllllHI Engelmann Spruce - SubAlpine Fir £441 Interior Douglas-Fir Lake I  I Watershed divide  A  / V River  Figure 2.3: Vegetation in drainage basins for Place and Eight M i l e Creeks (based on data from Ministry of Forests, Government of British Columbia).  14  2001  N  Kilometres  A  Geology HHI basalt-andesite flows, breccia, tuff, carbonate M I  blotlte-homblende granodiorite | calcareous sandstone, slltstone, shale, conglomerate, siliceous argillite, silty limestone  I ' I daclte, rhyolite, andesite, pyroclastlcs, tuff, argillite, conglomerate, calcareous slltstone, shale, sandstone volcanic sandstone, calcarenlte, turbldite, shale, limestone • 1 Lake [~| Watershed boundary  m  Ice  / V River  Figure 2.4: Geology of Place and Eight M i l e basins (based on data courtesy of D r . Journey, Geological Survey of Canada).  15  16  Figure 2.6: Place Creek in the valley bottom during spring snowmelt.  17  The tree-line within the basin is at approximately 1800 m, above which Alpine Tundra vegetation dominates.  Below the tree-line, land is forested with Engelmann Spruce  and Sub-Alpine F i r from 1000 to 1800 m, Coastal Western Hemlock from 700 to 1000 m and Douglas Fir below 700 m (Figure 2.3).  The geology of the basin is dominated by a dacite-pyroclastic complex above 2000 m and granodiorite below 2000 m (Figure 2.4).  Most of the sediment produced from glacial  erosion is thus derived from the former, while down-basin sediment sources originate from the latter, more resistant rock.  2.2.2 Eight Mile Creek The basin containing Eight M i l e Creek ( E M ) neighbours Place and also faces northwest (Figure 2.2).  It has an area of 25 k m and is essentially unglacierised, except for a small 2  portion of residual ice (0.18 km ) in the southwestern corner (Figure 2.2). 2  Figure 2.2 was  constructed using T R I M data from 1987 and it is probable that the area of ice cover is currently even smaller. Therefore, this ice likely has an insignificant hydrological influence at present.  From its headwaters to the valley below, Eight M i l e Creek also flows along a high gradient boulder-dominated channel that contains frequent large woody debris jams.  Near  the gauging station the channel width varies from 4 to 6 m (Figure 2.7). The basin is actively logged at present and a 5 km logging road that has three minor tributary roads stretches from the mouth to approximately half way along the length of the basin. The effect of logging on streamflow characteristics was not, however, investigated.  18  Below the measurement site, the  Figure 2.7: Eight Mile Creek at low flow  19  creek joins Gates Creek, which flows northeast from Gates Lake. Vegetation for the basin is as described in Sub-section 2.2.1.  The geology of Eight M i l e is approximately 75% dacite-pyroclastics.  The lowest  portion of the basin is a less-resistant and more permeable sandstone and conglomerate, while the highest portion is mostly tuff and volcanic sandstone (Figure 2.4).  Thus, although  Place and Eight M i l e are neighbouring basins, their dominant geology is different.  2.3 Streamflow measurement 2.3.1 Stage measurement Stage was recorded at stable gauging sites (established at Place Creek by R D Moore in 1999 and at Eight M i l e Creek by R D Moore and the author in 2000) with a Campbell Scientific C R 1 0 X datalogger attached to a submerged Keller PSI pressure transducer that was housed in a stilling well secured firmly in the stream. averaged every 10 minutes.  Stage was measured every 10 seconds and  Stage was also measured manually, from the top of the stilling  well to the water level within it. A n insulated wire attached along the length of a ruled steel rod was connected to a multi-meter and stage was recorded by eye when resistance changed from an off-scale reading as the exposed end of the wire touched the water and completed the circuit. These measurements provided a check on the stability of the data logger-transducer system.  2.3.2 Discharge measurement For small creeks with steep channel gradients and boulder-bed morphologies, salt dilution gauging has become the preferred gauging technique over the traditional velocity-area method (Elder et al., 1990).  For Place Creek and Eight M i l e Creek, 5 L of salt solution,  20  approximately 100 g/L, was injected as a "slug."  Stream electrical conductivity was  measured approximately 75 m downstream from the injection point to allow complete mixing of the solution. From measurement of the conductivity wave as it travelled downstream, and background conductivity, discharge was calculated from  Q =  fV  —  ,  (EC -EC )D ave  (2.1)  b  where Q is stream discharge, /  is the slope of the relation between E C and relative  concentration, determined by field calibration, V is the volume of the injected salt solution, s  ECave is the average conductivity during the salt wave passage, ECb is the background conductivity and D is the duration of the salt wave at the measurement point.  2.3.3  Rating curves  Stage-discharge rating curves were established for Place Creek (Equation 2.2) using 20 salt dilution gauges from 2000 and 5 from 1999 and for Eight M i l e Creek using 19 salt dilution gauges from 2000. A shift in the Eight M i l e rating curve required two rating relations to be established (Equation 2.3a for the original relation and 2.3b for the adjusted relation). A l l rating curves were fitted using the power function non-linear trendline tool in Excel2000.  Q = 1.383S ' 1  Q = 0.2845  Q = 0.688S  (2.2)  856  (2.3a)  2485  1 8 7 8  ,  (2.3b)  where Q is discharge (m /s) and S is stage (mV). The rating relation scatter plots are shown 3  in Figure 2.8.  21  0.5  1  1.5  2  2.5  3  3.5  Stage (mV)  (b) Eight Mile Creek (new) Q = 0.688S""" ^ = 0.982  •  /  •  Stage (mV)  Figure 2.8: (a) Place Creek and (b) Eight M i l e Creek stage-discharge rating relations. Eight M i l e has two plots (old and new) because of a shift in the original rating curve resulting from movement of the stilling well during high flows in early July 2000. 22  2.4 Suspended sediment concentration The objectives of this thesis did not require sampling S S C at a fixed interval. Rather, a range of sample intervals, over a number of field campaigns from mid-May through to the beginning of October, was employed to capture a variety of hydrologic conditions.  Sample  intervals thus ranged from 1 to 12 hr, during various periods of the field season.  Stream  water samples of 150 to 200 m L were collected at the gauging sites, at Place Creek and Eight M i l e Creek, using an I S C O (model 2900) automatic pumping water sampler which was located on the channel bank. The intake of the 7.6 m of Teflon 1 I S C O hose was fixed in the stream cross-section, approximately 1 m from the bank. Sample volumes were measured to ± 1 m L using a 250 m L cylinder before being vacuum filtered by a hand pump through individually pre-weighed numbered 0.7 um glass-fiber filters. Care was taken to ensure that no sediment remained in either the cylinder or sample bottle by rinsing both with the filtrate and re-filtering using the same filter. marked sample bags.  Sediment-laden filters were stored in individually  Filters were oven-dried in the laboratory at 1 0 5 ° C for 1.5 hr, both  before and after being used in the field.  Mass of dry sediment was obtained using an  electronic balance, by subtraction of the original mass from that of the dried sediment laden filter.  Organic proportions of S S C were determined by ashing the dry sediment laden filters at 5 5 0 ° C for 1.5 hr and subtracting the mass of the ashed from the pre-ashed filter. Adjustment for loss in water of hydration when ashing filters at 5 5 0 ° C was calculated from the average loss in mass of 30 randomly selected ashed sample filters, from each site, that had been used in the field and recorded a S S C of 0 mg/L.  23  It is noted widely in the literature that proglacial streams are highly turbulent and that fine suspended sediment is, therefore, well-mixed in the vertical and horizontal water profiles (Gurnell et al., 1992b; Wren et al., 2000).  Such claims support the use of fixed depth  sampling in this thesis. T o substantiate such inferences, 32 D H 4 8 depth-integrated samples were taken at various times within the season on either side of the cross-section of Place Creek in which the I S C O intake hose was fixed.  From the right bank, 24 samples were  taken, but only eight were taken from the left bank because of difficulties in crossing the stream at high flows. A l l filtering and laboratory procedures followed those described earlier in this section.  Samples from the outlet of Place Lake were taken using an improvised depthintegrated grab technique, where 500 m L polyethylene bottles were lowered and raised by hand, at an arm's length, through approximately 80% of the water depth with the bottle neck facing upstream. Sample volumes were approximately 300 m L . Again, subsequent filtering procedures were conducted as described for Place Creek.  2.5 Water chemistry 2.5.1  Electrical  conductivity  Electrical conductivity (EC) of the I S C O water samples was measured using a W T W L F 340 E C probe, which corrected E C to 2 5 ° C using a non-linear function.  The E C of sampled  stream water containing fine sediment may, however, continue to change because of solute acquisition from water-rock flour chemical reactions (Brown et al., 1994; Brown et al., 1996). Periodic spot measurements of stream E C were thus conducted at the. same time as samples were being collected by the I S C O sampler. The relation between E C of the I S C O samples and stream water followed the functions 24  EC  = (1.015 • EC  EC  = (1.056 • EC  PC  EM  )-0.438  ISCOpc  )  ISC0EU  - 4.563  (r = 0.996, S E E = 0.3, n = 17)  (2.4)  (r = 0.992, S E E = 0.5, n = 6),  (2.5)  2  2  where ECPC and ECEM, and ECISCOPC and ECISCOEM are the E C of Place Creek and Eight Mile Creek, and of the I S C O samples at Place and Eight M i l e Creek, respectively. Although time between I S C O stream sampling and ECisco measurement varied from 5 minutes to 1 week, Equations 2.4 and 2.5 provide calibration factors to adjust for any post-sampling geochemical reactions.  2.5.2 Cations  and  anions  Stream samples for chemical analyses were taken throughout the season at Place and Eight M i l e Creek. Sampling was conducted as close as possible to maximum and minimum daily flows, to capture chemical conditions close to high and low flows. In total, 66 samples were taken at these two sites. A l l 125 m L polyethylene sample bottles had been soaked in warm soapy water and were pre-washed with a 30% hydrochloric acid solution and then with distilled water, before being taken into the field.  Having rinsed all filtering equipment with stream water,  approximately 30 m L of water was scooped from the stream into the body of the filterer, swilled around and immediately vacuum filtered by a hand pump through a 0.45  |im  cellulose nitrate membrane that had been placed in the filter unit using a clean knife to avoid hand-contamination. The filtrate was used to rinse the inside of the unit and then discarded. This procedure was repeated three times.  Approximately 100 m L of stream water was then  filtered and part of this filtrate was used to rinse the sample bottle three times before filling the bottle completely, and sealing and labelling it. 25  A l l samples were kept out of direct  sunlight and refrigerated in the dark within 2 hr of collection. Sample storage before analysis ranged from 6 to 9 months.  Samples were analysed by Dr. G . Brown at the Centre for Glaciology, University of Aberystwyth, Wales. Major cation concentrations ( C a , K , M g 2+  +  2 +  and N a ) were determined +  using Atomic Absorption Spectrophotometry and Inductively Coupled Plasma analysis, and major anion concentrations (CI", N O 3 " and SO4 ") using Flow Injection Analysis. During air 2  2  and road transportation to Wales, samples would have been exposed to variations in ambient temperature, but this was unavoidable given the distance between sampling and analysis locations.  2.6 Meteorological variables A tipping bucket rain gauge, attached to a H O B O event recorder, was installed in a clearing opposite the gauging site at Place Creek. The H O B O logger recorded the time of each tip, which represented 0.25 m m rainfall.  After downloading the event data, number of tips per  hour were calculated using a Fortran programme written by Vincent Kujala, Department of Geography, University of British Columbia. A bulk precipitation gauge was also positioned at the same site and emptied at the beginning of each field campaign in case the tippingbucket gauge failed. A t the end of the season no known problems with the tipping bucket gauge were encountered.  Hourly air temperature data from Phelix Creek, located approximately 30 km northwest of Place and Eight M i l e Creeks, was provided by Brian Menounos, Department of Geography, University of British Columbia.  26  Chapter 3 Results 3.1 Overview of study period The year 2000 was characterised by near-normal snowpack conditions, based on data from the Tenquille Lake snow course located about 20 km N W of Place Glacier at an elevation of 1680 m (Figure 3.1a). Based on data from Pemberton, located 18 km S W of Place Glacier at an elevation of about 200 m, summer precipitation and summer air temperature conditions were also near-normal in 2000 (Figure 3.1a, b).  3.2 Data quality 3.2.1 Discharge measurement For Place Creek, the values of r  2  and Root Mean Squared Error ( R M S E ) for the stage-  discharge rating curve (Figure 2.8a) were 0.978 and 0.11 m /s respectively, indicating a 3  reasonable precision.  Figure 3.2a also indicates that no shift in the pressure transducer  occurred during the season.  At Eight M i l e Creek, high flows during early July induced movement of the stilling well.  Although the stilling well was re-fastened to the bank in approximately the same  position (indicated by the stability of the manual stage-mV scatter plot (Figure 3.2b)), a shift in the stage-Q rating curve was experienced (r = 0.995 and 0.982, and R M S E = 0.08 and 2  0.10 m /s for old and new rating curves, respectively) (Figure 2.8b). 3  27  -Normal - » - 2 0 0 0  1999  1998  500 -i 450 -  A  A.,  400 -  ? 0 (A 3  O  350 •  "A.  300 • " " ^ ^ ^  ''A  250 •  O  3 200 o  c V)  " " • A * \  A'\^-  150 100 50 0• 20-Dec  9-Jan  29-Jan  18-Feb  10-Mar 30-Mar  19-Apr  9-May  29-May 18-Jun  8-Jul  25 •  May  Figure 3.1: Normal (1984 to 2000) versus 1998, 1999 and 2000 (a) monthly snow course, (b) summer monthly air temperature (TA) and (c) summer monthly precipitation (P).  Data  provided by (a) B . C . Environment Historical Snow Survey Data and ((b) and (c)) D . Hutchinson, Hydrology Applications Group, Environment Canada. 28  Figure 3.2: Manual stage-mV scatter plots for (a) Place Creek and (b) Eight M i l e Creek.  29  3.2.2 SSC measurement Error analyses were conducted on a range of S S C measurements.  Sources of error were the  precision of the electronic balance ( ± 0 . 1 mg) and volume measurement ( ± 1 m L ) . For low values of S S C (< 5 mg/L) error was around 10% (0.5 mg/L); for mid-range values (40-60 mg/L), about 1.5% (0.8 mg/L); and for high values (>120 mg/L), approximately 0.7% (0.7 mg/L). Thus measurement errors were small.  A n important question is whether the I S C O samples were representative of the stream cross section. Table 3.1 summarises results of a paired sample T-test that was used to test the hypothesis:  Ho :  p: = 0 D  H i : P-D  ^  0,  where p, is the mean of the differences between paired I S C O and right bank D H 4 8 samples. D  Table 3.1: Results of I S C O - D H 4 8 Paired Sample T-test  T-test information  Value  n  24  Degrees of Freedom (df)  23  T-statistic (T )  -1.206  t  T-critical for T w o Tailed Test ( T )  2.069  P (1T I < T )  0.240  c  T  c  The T-statistic is less than the critical T-value (Table 3.1) and so there is no evidence to reject Ho at a significance level of 5%, thus indicating that the I S C O sampler provided unbiased estimates of S S C in comparison with the D H 4 8 sampler.  Figure 3.3a shows a -  scatter plot of I S C O versus D H 4 8 values of S S C and the proximity of the points to the 1:1 line.  30  35  Figure 3.3: (a) I S C O versus D H 4 8 values of S S C for all samples and (b) left versus right bank D H 4 8 values of S S C .  31  T o test the hypothesis that there was a difference in S S C between the I S C O sampler and the eight paired D H 4 8 samples, S S C values were analysed by a two-way Analysis O f Variance ( A N O V A ) using the General Linear Model function in S Y S T A T 9.0. The factors were time of sampling (8 levels) and type of sample (3 levels: I S C O , D H 4 8 right bank, D H 4 8 left bank). Results indicated that time of sampling was significant (p<0.001) and type of sample was insignificant (p=0.502).  There is, therefore, evidence that sampling location does not  have a significant influence on S S C . Again, Figure 3.3b shows a plot of left versus right bank S S C , indicating that all points plot close to the 1:1 line.  3.3 Streamflow and SSC variations 3.3.1 Place Creek discharge and SSC time series A continuous record of discharge was obtained from 12 M a y to 31 N o v 2000, after which streamflow was minimal and thus not of interest for the purpose of this thesis (Figure 3.4). The S S C series stretched from 18 M a y to 8 Oct 2000, but was discontinuous.  Hourly air  temperature (TA) and daily precipitation (P) records spanned the periods 19 M a y to 26 Sept and 19 M a y to 30 Nov, respectively.  Through M a y and June, background discharge rose and was punctuated by rainfall events on 20 to 21 M a y (11.4 mm), 26 to 27 M a y (18.5 mm), 5 June (13.7 mm), 18 to 19 June (7.4 mm) and 1 to 2 July (12.4 mm), and a snowmelt event on 28 to 29 June (Figure 3.4).  Following a lull in discharge during early July, discharge rose again and was  characterised by an increasing diurnal oscillatory amplitude. From 27 to 28 July, 31.3 m m of rainfall produced the largest recorded discharge of the season (3.83 m /s). 3  During early  August discharge fluctuated daily, from around 2.0 to 2.8 m Is. Cooler air temperatures then  32  25  4-May  jli!  3-Jun  J_J  UL 3-Jul  2-Aug  1-Sep  1-Oct  iL  31-Oct  30-Nov  Figure 3.4: Time series of discharge, S S C , air temperature and precipitation for Place Creek  (2000 season).  33  prevailed from 10 to 21 A u g and average daily discharge fell from approximately 2.4 to 1.3 m /s. After a period of renewed glacial melt from 22 to 26 A u g , induced by warm air 3  temperatures, background discharge followed a seasonal recession trend and was again punctuated by rainfall and late glacial melt events. From 10 to 30 Nov, discharge remained minimal at <0.1 m /s. 3  A tentative relation between S S C and discharge can be noted from Figure 3.4, in which peaks in S S C generally corresponded to peaks in discharge. Following the first two S S C events of the season, during mid-May and at the beginning of June, S S C declined from around 40 to 10 mg/L, despite similar discharge event magnitudes.  During July, S S C rose  and peaked at 123 m g / L during the storm event on 27 to 28 July.  For the first 10 days in  August, S S C mimicked diurnal flow cycles, but S S C then decreased in mid-August during a dry period of cooler air temperatures.  N o values of S S C were recorded from 3 to 25 Sep  because of time constraints on access to the field site. The last field campaign of the season recorded S S C during the rain event of 28 to 30 Sep, during which similar peak concentrations (43 mg/L) to those at the same discharge during the rain storm on 29 A u g , were measured (44 mg/L).  3.3.2  Eight Mile  Creek discharge  and SSC time  series  The streamflow record for Eight M i l e Creek covered the period from 14 M a y to 31 Nov 2000, again beyond which discharge was minimal investigation (Figure 3.5).  and thus not of interest for this  Following a period of high discharge during late June and  subsequent shift in the stage-discharge rating curve (Figure 2.8b), a break in the discharge series was experienced from 2 July 00:00 to 6 July 15:50, after which the stilling well was  34  4-May  3-Jun  3-Jul  2-Aug  1-Sep  1-Oct  31-Oct  30-Nov  Figure 3.5: Time series of discharge, S S C , air temperature and precipitation for Eight M i l e  Creek (2000 season).  35  re-secured. The record of S S C ranged from 19 M a y to 3 Sep 2000 and as for Place Creek, the series was sporadic.  A t the start of the season, discharge rose from 0.6 m /s to around 5 m /s at the end of 3  3  June. During this period streamflow responded to the same rainfall and snowmelt events as •a  Place Creek did. Following a drop in discharge to just above 2 m /s at the beginning of July, background flow was approximately 2.5 m /s, but displayed similar synoptic scale increases 3  as for Place Creek. Despite maximum daily air temperatures of above 30 ° C during 24 to 30 June, streamflow dropped rapidly over this period. Background levels of discharge remained at about 2.5 m /s through July. 3  The storm event on 27 to 28 July, which produced the  maximum seasonal discharge in Place Creek, did not have such a significant impact in Eight M i l e Creek.  Discharge then declined rapidly, from 6 to 22 A u g , and then more steadily  through to the end of the season.  Again, during this latter portion of the season, synoptic  scale storm events induced responses in discharge similar to those in Place Creek. Background S S C in Eight M i l e Creek was generally <10 m g / L , except during periods of elevated discharge.  The maximum recorded S S C over the measurement season was 56  mg/L, which occurred during the first rain event of the season (20 to 23 May). The following three events (same dates as for Place Creek), all of similar peak discharges (around 5 m /s), were associated with progressively declining S S C , from >51 m g / L during the second rain event, to 23 and then 8 m g / L for the third and fourth events, respectively.  The only other  significant periods of suspended sediment transport occurred during the second of two discharge events, closely spaced over time, on 23 and 28 July.  For the remainder of the  season, minor fluctuations of S S C above 0 m g / L can be attributed to small rain events and possibly measurement error. 36  3.3.3  Standardised  discharge  series  During the snowmelt and autumn portion of the season, streamflow was greater in Eight Mile than in Place Creek because the basin area of Eight M i l e is almost double that of Place Creek and thus collects more snowmelt and rainfall runoff. Time series of standardised discharge (Q ), that is discharge per unit drainage area (m s" km ), reveal a good correspondence 3  1  2  s  during the major period of snowmelt and spring storm events, from 14 M a y to June 30 (Figure 3.6). For this period, the linear relation between the two time series is  Q (EM) = 0.02 + 0.82Q (PC) S  S  (r = 0.949, SEE=0.01), 2  (3.1)  where S E E is the standard error of the estimate. From the beginning of July, melt from Place Glacier began to make a significant contribution to streamflow in Place Creek as inferred from the increase in diurnal discharge cycles. In Eight M i l e Creek, however, flows generally decreased from July to October. The differences in standardised discharge for this period probably reflected the effects of glacier melt in Place Creek. From the end of October until the end of the season, Eight M i l e yielded more discharge per k m than Place Creek, which was probably because there are greater volumes of groundwater in the former basin.  3.4 Seasonal hydrograph division 3.4.1  CUSUMplot  of Q  ratio  The use of cumulative sums ( C U S U M S ) provides a powerful technique for identifying regime shifts in time series (Mac Nally and Hart, 1997). For this sub-section, the discharge series for Place Creek will be sub-divided using the time-ordered E M / P C discharge ratio series. The C U S U M S are thus defined as the running tally of the deviations of each value of the E M / P C discharge ratio from the mean of the series, such that  37  0.00 0.30 0.25+  |i  ii  (b)  ^0.20 0.15  S  !  go.io 111  0.05  -400 4-May  3-Jun  3-Jul  2-Aug  1-Sep  1-Oct  '  31-Oct  30-Nov  Figure 3.6: Standardised discharge time series for (a) Place Creek and (b) Eight M i l e Creek. Arrows on the (c) C U S U M S plot indicate sub-seasonal transition dates.  38  Sr = S _ r  l  (3.2)  +  Qr  where S is the r* partial sum, S .i is the previous partial sum, and Q is the r value of and th  r  r  r  Q is the mean of the EM/PC discharge ratio series. The CUSUMS can then be plotted as a r  time series (Figure 3.6c). Abrupt shifts in the slope of the plot indicate that the ratio between Eight Mile and Place Creek discharges has changed. According to the greatest changes in gradient of Figure 3.6c, the season can be sub-divided into a 'nival' (N) period (start of the season to 17 June), a 'nival-glacial' (NG) transition period (18 June to 23 July), a 'glacial' (G) period (24 July to 2 Oct) and an 'autumn recession' (AR) period (3 Oct to 30 Nov). Discharge and SSC varied amongst the 4 sub-seasons (Table 3.2). Table 3.2: Descriptive Statistics for Place Creek Sub-seasons. SSC (mg/L) Mean Standard Deviation Minimum Maximum No. Data Points  3.4.2  Relations  Q (mVs)  N  NG  G  AR  N  NG  G  AR  9.9 11.7 0.0 70.6 104  16.3 7.6 1.5 42.5 233  27.7 12.9 8.8 122.6 388  10.6 2.9 7.4 16.8 11  0.84 0.50 0.00 2.55 5238  1.71 0.39 0.97 2.54 5184  1.45 0.61 0.52 3.83 10224  0.23 0.19 0.02 1.33 8496  between  electrical  conductivity  and  discharge  Unfortunately, laboratory restrictions did not permit HCO3" to be determined and so a charge balance calculation for measuring chemical analysis accuracy could not be performed. Given that total anion concentrations ( 2 " ) should equal total cation concentrations (Z ), however, the +  strong relation between EC and S (Figure 3.7) indicates that EC was a valid measure of +  Total Dissolved Solids (TDS) for Place Creek.  39  40  T o support the validity of the C U S U M S method for defining sub-seasons of the Place Creek discharge series, the relation between E C and discharge for Place Creek, divided into the 4 sub-seasons, can be examined (Figure 3.8). In the literature, hydrograph separation using E C has been successful in separating hydrograph components (Teti, 1979; Laudon, 1995). From Figure 3.8 it appears that there are three distinct groupings of points: a concave-up relation for the nival period, a concave-up relation for the combined glacial and autumn recession periods, and an unstructured group for the nival-glacial transition. These groupings indicate that water sources and/or flow paths varied amongst the hydrologic sub-seasons defined by the C U S U M S technique, supporting their validity.  3.5 SSC-Q relations for Place Creek 3.5.1 Seasonal SSC-Q relation The time series of S S C and Q (Figure 3.4) indicated that there was a tentative relation between S S C and Q . Accordingly, the simplest way to predict S S C is from a S S C - Q rating curve.  A common technique for constructing the rating relation is to use Ordinary Least  Squares (OLS) regression, following the function  SSC, = b +b Q,, 0  (3.3)  1  where SSC, is the estimated value of S S C , subscript i represents a particular value at the time of sampling, and bo is the intercept and bi the slope of the function. The parameters for the fitted O L S equation for S S C are shown in Table 3.3.  41  30  • • 26  H  22  A  • Nival A NivalGlacial o Glacial o Autumn Recession  *A  A A A^  *A« *A ^  A ^ A  ^  &A* * A  AA  18H A  A O  14H o floo  0.0  0.5  1.0  1.5  2.0  %>  2.5  §  go  3.0  3.5  4.0  Q (m7s)  Figure 3.8: E C - Q relations for Place Creek. division dates defined in Figure 3.6c.  42  Data points are grouped by sub-seasonal  Table 3.3: Summary of O L S S S C - Q Seasonal Regression Model (values in brackets are the p-values for bo and bi).  n 686  r  2  0.445  SEE (mg/L)  F-ratio  9.8  548.94  bo (mg/L)  bi mg L V ( m s ) 3  -2.36  14.86  (0.026)  (<0.001)  1  Consistent with the relatively low r value, there is considerable scatter about the regression 2  line (Figure 3.9a).  Such low model accuracy is widely noted in the literature (reasons  outlined in Section 1.2).  The significance of the seasonal S S C - Q model and its predictive ability cannot, however, be correctly interpreted unless the assumptions of regression are met.  These  assumptions are: normality of the error terms, equal variances of the errors, independence between error terms, and linearity.  Plots of the residuals of the fitted model (Table 3.3)  reveal violations of the assumptions (Figure 3.9b, c, d). A s values of  values of e, (that is  SSC increase, absolute t  SSCr SSC ) increase, indicating heteroskedasticity (Figure 3.9b). Beyond t  SSC, values of 35 m g / L there are no negative values of e „ which provides evidence for a lack of fit in the model (Figure 3.9b). dependence) of the residuals.  Figure 3.9c illustrates autocorrelation (i.e., time  Over the entire season a trend towards higher positive  residuals, and thus S S C underestimation, can be seen. O n a shorter timescale, this trend is also punctuated by marked short-term S S C underestimation, particularly during the rainfall induced discharge events described in Sub-section 3.3.1.  Violations of the assumptions were tested statistically.  T o test for normality, the  Shapiro-Wilks (W) test (Neter et al., 1996) was used under the hypothesis:  43  44  Ho: error terms normal H i : error terms not normal  The probability value for the computed W-statistic was less than 0.001, indicating that at a significance level of 0.05, Ho can be rejected. Thus, it can be concluded the error terms are not distributed normally.  Homoskedasticity was tested for by generating a 'pseudo-White's' test statistic (x ) using S A S 6.12. The hypothesis tested was:  Ho: equal variances H i : unequal variances  The probability value for the % -statistic was less than 0.001. A t a 5% significance level % is 2  2  thus significant and the null hypothesis can be rejected.  It can be concluded that there is  evidence for heteroskedasticity.  For independence, a Durbin-Watson (DW) test was performed under the hypothesis:  Ho:p = 0 H  i :  p*0,  where p is the lag-1 autocorrelation amongst the error terms.  The D W test statistic was  0.269. A t a 5% significance level, the critical bounds for the test statistic were 1.65 for the lower and 1.69 for the upper bound (for >100 samples and 1 parameter (QJ).  The test  statistic is below the lower bound and so Ho can be rejected, thus suggesting the presence of residual autocorrelation.  45  Lastly, for lack of fit, the residual plot was assessed visually and concluded to be nonlinear at the upper end. attempt  to  reduce  The linearity assumption must, therefore, also be rejected.  heteroskedasticity  and normalise  the  variable distributions  transformation of the discharge and S S C data was performed.  The log  e  In an a  log  e  transformation  linearised the S S C - Q relation for high values of discharge, but a large number of low values of S S C then exhibited a non-linear relation with streamflow (Figure 3.9d). Additionally, the variances clearly became heteroskedastic.  Thus, transforming the variables did not appear to  provide a solution to the problems of non-linearity and heteroskedasticity.  3.5.2 Sub-season SSC-Q relations This section tests the hypothesis that the nature of the S S C - Q relation for Place Creek differs among the hydrologic periods defined in Sub-section 3.4.1. Given that the autumn recession period only has 11 data points, which are all spaced between 3 to 8 Oct, this sub-season will be ignored. Further analyses and reference to the "Entire season" thus refer to the other three periods.  Given that the data set has been sub-divided into three groups it is not possible to compare r values, either between sub-seasons or against the whole season, because the data 2  in each group are different.  B y introducing a dummy variable into a S S C - Q O L S model for  the first three sub-seasons, however, a single value of r can be generated and compared to the value for the Entire season.  Beginning with an O L S model for the Entire season (Table 3.4), a dummy variable for each sub-season was introduced, leading to the expression  SSC, = b + bi Qt +b S + b S + b QiS + b QiS , 0  2  u  3  i2  A  u  46  5  i2  (3.4)  .  where Su = 1 if nival period, otherwise 0; and 5,2 = 1 if nival-glacial period, otherwise 0. Equation 3.4 includes interaction terms for both the slope and intercept of the original Entire season rating function (Table 3.4).  The R value for the dummy-season model is 0.648, an  increase of 0.202 from the Entire season model, and the S E E has also decreased from 9.8 to 7.8 m g / L (Table 3.4).  Thus the dummy-season model has improved the accuracy of S S C  predictions. Statistical analysis revealed significantly different slopes and intercepts amongst subseasons (details provided in Appendix A ) .  Table 3.4: Summary of S S C - Q O L S Regression Models (values in brackets are the p-values for bo and bi).  Sub-season  n  r  1  SEE F-ratio (mg/L)  675  0.446  9.8  540.83  675  0.684  7.8  246.70  Nival  104  0.583  7.6  142.70  Nival-glacial  233  0.393  5.9  149.57  Entire season Dummy model  Glacial  338  0.516  9.0  358.43  b (mg/L) 0  bi mg I/V(m s ) 3  -3.20  15.31  (0.004)  (<0.001)  -  -  -18.10  31.32  (<0.001) -4.38 (0.012)  (<0.001) 11.86  -1.33  (<0.001) 15.67  (0.367)  (<0.001)  1  Comparing the residual plot for the sub-seasonal model against the Entire season model (Figure 3.10), it appears that the residuals are more homoskedastic and linear in the former, which indicates that seasonal division has also improved the validity of the regression assumptions.  From Figure 3.11, however, it is clear that although the different slopes and  intercepts of the sub-season models account for general seasonal differences in hydrologic control on proglacial suspended sediment  dynamics, over 30% scatter remains to be  explained. In order to account for remaining scatter, a number of other explanatory variables were introduced into each sub-season rating function. 47  E OT 0) DC  E, 75 3  •g 'in  v rr  10  15  20  25  30  35  40  45  50  55  60  Predicted SSC (mg/L)  Figure 3.10: (a) Seasonal versus (b) sub-seasonal S S C - Q relation residual plots.  48  x Nival A Nival-Glacial o Glacial  Glacial  £  40 H  Nival-Glacial  30 H  20 H  10 H  4.0  0.0  Q (m7s)  Figure 3.11: Sub-seasonal S S C - Q rating relations (note that S S C scale has been reduced so that the grouping of points can be more readily distinguished).  49  3.5.3  Multiple  regression  models  A suite of variables from which to select for Multiple Linear Regression ( M L R ) was derived from published literature and consideration of local factors that may have proglacial suspended sediment transport in Place Creek. variables are outlined below.  Reasons  influenced  for selecting these  Each variable has been drawn from relevant literature unless  otherwise indicated.  1. Q  2  was selected to account for the non-linearity present in Figure 3.11.  Having fitted a  non-linear function to the data set, in the form S S C = a Q , using the non-linear b  regression tool in S Y S T A T 9.0, an associated r  2  value of 0.454 did not show great  improvement over the value of 0.446 for the Entire season model (Table 3.4).  2. Q M . B y lagging Q at different intervals (current time (t) - 0.5 hr, t-1 hr, t-1.5 hr ... t-250 hr) and then separately adding each variable into regression function 3.3, the optimum lag for Q was selected as t-1 hr, based on R values. 2  The reason for lagging discharge  was to account for the residual autocorrelation noted in Figure 3.9c (Gurnell at al., 1994).  A s stated in Section 2.4, the S S C sampling interval was irregular.  irregularity complicates A R L M A modelling).  attempts to correct for residual autocorrelation (e.g.  Such using  B y lagging the diurnal phasing of discharge to the best-match fit  with the diurnal phase of S S C , it is possible, however, to account for some degree of the residual autocorrelation in the Entire season O L S function (Table 3.4).  3. Z Q . Cumulative discharge over the monitoring period (m xlO" ) was incorporated into the suite of predictor variables because it accounts for seasonal changes in suspended sediment supply. The assumption of E Q is, therefore, that processes governing seasonal  50  changes in suspended snowmelt,  sediment  supply are represented by the increasing flux of  glacial melt and storm runoff that passed through the drainage basin.  Specifically, Z Q acts as a surrogate for suspended sediment supply exhaustion (Hodson and Ferguson, 1999).  4. A Q 3 represents the rate of change of discharge over 3 consecutive hours and was obtained by subtracting from current values discharge, values recorded 3 hr previously.  It was  thus positive during the rising hydrograph limb and negative during the falling limb. Recent literature suggests that S S C is more dependent on A Q , rather than the current discharge, Q (Willis et al., 1996).  The optimum time-interval selection for A Q was  obtained using the same method outlined for Q . i , with AQ0.5, A Q i , AQ1.5 . . . AQ250t  5.  SSC( ef)r  Many S S C M L R models published in contemporary literature adopt various  lagged values of S S C as autoregressive predictors (e.g. Hodson and Ferguson, 1999). Such variables, however, limit the predictive capacity of the model because an intense sampling programme must be conducted to obtain lagged values of S S C .  A n alternative  approach would be to include a measured, reference value of S S C for each day.  If  successful, this approach could provide an effective trade-off between minimising costs of monitoring while maintaining accuracy of predicted S S C values.  Inspection of  diurnal-scale S S C time series for Place Creek indicated that peak values  generally  occurred at approximately 18:00. B y using the daily value of S S C , measured at 15:00, to predict values between 9:00 that day and 8:00 the following morning (approximate timing of beginning and end of diurnal phase of S S C ) , a low maintenance sampling programme could be used to 'tune' the overall model to the physical state of the sediment  51  supply system for each day.  This is a variable that has not, to the knowledge of the  author, been used in previous studies similar to this one.  6.  P24 was calculated by summing precipitation (mm) over the 24 hour period previous to the time of each S S C sample. Evidence from the literature suggests that precipitation is an important determinant of S S C during high discharge storm events (Hodgkins, 1999).  7.  T  A  was used as a surrogate for glacial melt and thus accounted for diurnal fluctuations in  discharge during the dominant summer ablation period (Aizen et al., 1996).  Values of variables 1 to 7 were entered into a S Y S T A T 9.0 spreadsheet for each of the 675 values of S S C . T o construct M L R models, the data set was again divided into its respective  sub-seasons and variables were selected by backwards stepwise interactive  selection. The selection was interactive because SSC( f) and A Q 3 were 'forced' into each re  model. After analysing similar M L R models from the literature, it was expected that SSC( f) re  and A Q 3 should, physically, provide the most explanation for the scatter in Figure 3.11. Selection criteria were an Alpha-to-Enter level of 0.049 and an Alpha-to-Remove level of 0.050. The resulting sub-season M L R models are given in Table 3.5.  The number of data  points used to derive each model is less than those used for corresponding sub-season models (Table 3.4) because on occasions S S C was not measured at 15:00 and thus a value of SSC( f) re  could not be used with other values of S S C that same day.  52  O  1—1  <  9 §  9 § in o d  cn  CO  CN O  9 §  CN  i  d  o  o o o  CO  ^  d  oo  > n o O  cn  o  cn  <d  | > CN  i  o  o o m o  in  •  ON  o o il n o -H cn  d  NO  O  cn  o  ^ r-NO  d  00 ON  00 CN  ©  O  d  d  O  ^O  o  CN  o  O  00  P  O  -H  ON  /—V  CN  (000  1^  O o  (000  CN  u  CN  CN  055)  Sf  O o o  c  CO  1  CN  d  CO  CO  NO  8  "i °  CN  o  OO CO  ON <^ CM o  cn  O CN CN  d  vq o oo o  NO  O CN  NO  cn  NO  ON  o  o  d  oo  oo  cn  o  O  ON  fi o vi  CN CN  cn  _c3  O cd  US CU CO I  .a fi  o  cn  C3 >  I  2 'o  >  3  53  03  In Sub-section 3.5.1 it was shown that the original O L S model (Table 3.3) failed to meet all of the regression assumptions. Under the same hypotheses and using the same tests as those described in Sub-section 3.5.1, Table 3.6 summarises the sub-seasonal M L R model assumption test statistics.  Table 3.6: Summary of M L R Assumption Test Statistics (PRESS is the Prediction Sum of Squares).  Sub-season  Pr<W  Pr>%  DW  PRESS  Nival  0.1126  0.2536  1.176  1839  Nival-glacial  0.1325  0.0257  1.139  2779  Glacial  <0.001  0.0304  0.863  11860  2  For the nival period, the normality and equal variances tests were passed.  For the  nival-glacial sub-season, the data were normally distributed, but slightly heteroskedastic. For the glacial period, the data were again slightly heteroskedastic, but not normally distributed. The upper and lower bounds of the Durbin Watson test statistic were, for each sub-season, 1.78 and 1.57, respectively.  There is, therefore, evidence of minor autocorrelation in each  model.  The validity of the models can be assessed from the values of the P R E S S statistic, which are calculated by summing the squares of the differences between each observed value of S S C and its predicted value based on a regression function developed from the remaining n - 1 values of S S C . The values of P R E S S are thus a useful measure of how well the use of the fitted values for a subset of each sub-seasonal M L R model can predict the observed values of S S C (Neter et al,  1996).  statistic (not shown) was 65242.  For the Entire season model (Table 3.4), the P R E S S Although P R E S S values cannot be directly compared  between models, since they are all based on different samples, the values for the M L R nival 54  and nival-glacial sub-seasons (Table 3.6) are much lower than for the Entire season model, indicating that both M L R models have small prediction errors. For the glacial period the prediction errors are larger than for the other two periods, but still represent a great improvement in the predictive capacity over the O L S Entire season model.  3.5.4 SSC-Q hysteresis relations Several types, strengths and scales of S S C - Q hysteresis can be noted according to seasonal, synoptic and diurnal changes in dominant hydrologic sources. T o allow process inferences to be made from the hysteresis loops in Place Creek, examples from different sub-seasons and synoptic-scale events will be compared.  It must be noted that quantitative measures of  hysteresis have not been developed and so qualitative analysis must suffice.  T o ensure  minimum subjectivity, the direction of the loop will describe its type and relative enclosed area its strength. Figure 3.12 gives examples of hysteresis loops with their associated discharge and S S C time series plots.  For storm runoff events, hysteresis was strong and clockwise early in the season (Figure 3.12a), but later in the season was weaker and anticlockwise (Figure 3.12b).  For  snowmelt events, both strong and weak clockwise hysteresis occurred within sequential days during late June (Figure 3.12c). O n the first day, discharge and S S C peaked at approximately the same time (27 June at 00:00). Although they again peaked at similar times the following evening (27 June at 21:00), discharge did not begin to fall for another 6 hr, therefore inducing stronger clockwise hysteresis.  Further complexities can be noted in Figure 3.12d, in which  the latter part of the loop is vertical and anticlockwise, indicating a negative S S C - Q relation. Another interesting loop on 20 to 21 July (Figure 3.12e) depicts a 'figure 8' formation  55  12:00  0:00 12:00 0:00 12:00  0:00  20-22 May 2.0  1.8 1.6  if  S, o  1.4  V 12:00  1.2 12:00  0:00  29-30 Aug 2.2  12:00  0:00  12:00  12:00 0:00  26-28 June  Figure 3.12: Examples of S S C - Q hysteresis for: ((a) and (b)) storm runoff events; (c) a snowmelt event; and ((d) and (e)) more complex S S C - Q relations. SSC-Q  plots indicate  direction  of  hysteresis  (filled  and open  Arrows on  arrows  discharge cycles/events). Time plots (right) correspond to hysteresis loops.  56  divide  1.8  2.0  2.2  2.4  2.6  12:00  18:00  Q (m /s)  0:00  20-21 July  3  — Q SSC  Figure 3.12 (continued)  57  6:00  12:00  (Williams, 1992), combining initially clockwise and then anticlockwise sub-loops.  The  former occurred during high discharge and the latter at lower discharge.  During the glacial season, the dominant diurnal S S C - Q hysteresis was strong and clockwise (Figure 3.13), indicating that S S C peaked before discharge. O n a synoptic scale, hysteresis was also clockwise during mid to late August, as indicated by the clockwise progression of the loops themselves.  T o aid the comparison of hysteresis loops in Figure 3.12, Figure 3.14 illustrates the pattern of S S C and its proportion that was organic material, for four discharge events. A l l four plots indicate that organic S S C was approximately constant throughout the events, but from event a through d, a progressive decrease in organic S S C occurred.  Given that total  S S C rose and fell in each event, while the percentage organics remained approximately constant, the latter had a negative relation with S S C (Figure 3.14).  This percentage was  much lower (10-20%) during the fourth event than for the previous three events (20-100%). Since this last event was dominated by snowmelt, while the others were principally rainfall induced, Figure 3.14d confirms that sediment sources differed for different sources of runoff to the main channel.  58  Figure 3.13: Diurnal and synoptic scale S S C - Q hysteresis during the main ablation period (arrows indicate direction of synoptic scale hysteresis). clockwise direction.  59  A l l individual loops have a  60  3.6 Upstream-downstream comparisons of SSC This section focuses on changes in S S C between Place Lake outlet and the valley site during two glacial melt cycles, one at the beginning and one at the end of August (Figure 3.15). A s a result of time constraints, an adequate stage-discharge rating curve could not be established at the lake outlet.  Data from 4 to 6 A u g clearly show that S S C at the lake outlet was  approximately constant over the sampling period, at around 55 mg/L. In the valley, however, S S C followed a diurnal increase-decrease pattern, with a peak value of 60 m g / L at 20:00 and minimum value of about 30 m g / L at 10:00 on 5 A u g . Respectively, these times are also similar to those of the maximum and minimum values of valley discharge.  A t the end of August the S S C data series were complicated by a rain event, which appears to have caused an increase in S S C at the lake outlet, from 59 to 75 m g / L , during the period of rain (01:00 to 19:00 on 29 Aug). During this period, sediment runoff was observed from morainal slopes at the lakeside.  In the valley, minimum S S C values were similar to  values recorded in early August, but maximum values were approximately 40 m g / L and thus less than upstream values.  For the case of approximately constant discharge at the lake outlet (i.e., early August upstream-downstream field campaign), an approximate suspended sediment budget was constructed, assuming that the total outflow over a diurnal cycle is the same at Place Lake outlet and Place Creek at the valley site. The sediment fluxes can thus be expressed as  (3.5)  (3.6)  61  Figure 3.15: Upstream-downstream comparison of lake outlet S S C , and valley discharge and S S C for (a) early August and (b) late August 2000.  62  where SSF o  and SSF  PL  V  are the suspended sediment fluxes over a daily discharge cycle at  Place Lake outlet and the valley bottom, respectively, SSCPLO  is the average hourly S S C at  the lake outlet over a daily discharge cycle, and Q(PC) is average hourly discharge at the valley bottom. Equations 3.5 and 3.6 can be approximated as  SSFPLO = SSCPLO • £  SSFv =  (3.7)  (<2 • At)  ^(Q-SSC-At),  (3.8)  where At is the change in time (seconds).  For the lake outlet, the total daily sediment flux for 5 A u g 12:00 to 6 A u g 11:59 was 12300 kg (12.3 t) and for the valley bottom was 9300 kg (9.3 t).  Approximately, the  suspended sediment budget at the basin outlet was -3000 ± 2 0 0 kg (-3 ± 0 . 2 t).  Errors for  suspended sediment fluxes were calculated considering errors in S S C ( ± 2 % ) and discharge ( ± 5 % ) measurement.  3.7 Chapter summary Sub-division of the seasonal discharge hydrograph for Place Creek was comparison of its standardised discharge with that for Eight M i l e Creek.  achieved by  The C U S U M S  technique revealed four clear sub-seasons (nival, nival-glacial transition, glacial and autumn recession).  For each sub-season, four distinct groups of points on an EC-discharge plot  supported the method of hydrograph separation. Relations between S S C and discharge were highly scattered, but showed improvement when data were divided into nival, nival-glacial and glacial sub-seasons.  Multiple regression models further improved the accuracy of S S C  63  prediction and the ability of the models to meet the regression assumptions. Hysteresis loop characteristics for diurnal scale SSC-Q relations provided additional evidence of the influence of different sources of runoff on suspended sediment transport within glacier-fed rivers. Finally, comparison of upstream-downstream SSC suggested that SSC at the lake outlet was approximately constant, but that at the valley bottom SSC fluctuated with diurnal discharge in the glacial period. Further evidence of the significance of in-channel sediment storage on the downstream delivery of sediment in Place Creek was provided by a negative daily sediment budget at the basin outlet.  64  Chapter 4 Discussion Discussion in this chapter will be focused on Place Creek and Eight M i l e Creek will only be used as a reference basin.  4.1 Hydrologic control on SSC 4.1.1  Seasonal scale  From the groups of points in the E C - Q relation (Figure 3.8) it can be hypothesized that there are three main hydrologic sources of streamflow within Place Creek: (1) water from Place Lake with an E C of about 10 to 12 uS/cm (measured at the lake outlet during early and late August 2000 by the author); (2) water from a baseflow reservoir with values of E C around 40 uS/cm (as measured on 16 Jan 2000); and (3) high-EC water that is generated during the nival period (25 to 30 u,S/cm). Unfortunately, this hypothesis cannot be addressed with the current data.  The distinct groupings of points in Figure 3.8, however, suggest that the  relative contributions of these sources vary amongst sub-seasons.  The significance of the  dummy variables and their interaction with Q in function 3.4 indicates that seasonal changes in the hydrologic sources of discharge (i.e., snow-melt, glacial-melt and a nival-glacial combination) alter the simple S S C - Q rating relation. Different areas of a catchment deliver different quantities and quality of sediment to a stream channel. In the initial part of the nival period, the slope of the nival S S C - Q regression function was steeper than for the nival-glacial and glacial functions (Figure 3.11), indicating that sediment entrainment efficiency greatest early in the season, as previously found by Chikita (1993).  was  The decrease in the  regression slope during the nival-glacial transition period may reflect the retreat of the  65  snowpack to the upper portion of the basin above the tree-line (Figure 2.3), where soils were likely thinner and provided less sediment for entrainment.  Additionally, during this  transition period, the glacier was only making a small contribution to streamflow,  as  indicated by the similarity of standardised discharges between Place and Eight M i l e Creeks (Figure 3.6), and thus suspended sediment supply from Place Lake was also minor. The regression slope for the glacial period was similar to the slope for the nival-glacial subseason, but the intercept was greater for the former period. It could be argued, therefore, that the sediment discharge mechanisms during these two periods were analogous, but that during the latter, the increase in glacier-melt discharge from Place Lake supplied more suspended sediment than for the same discharge during the former sub-season.  A s a result of differences in the S S C - Q relation through sub-seasonal changes in water sources, the need to sub-divide seasonal hydrographs in glacierised basins is of fundamental importance to any investigation of suspended sediment transport within such catchments. Previous studies were highlighted in Section 1.2 as using subjective methods for hydrograph division. In this study, the comparison of Place Creek with a similar basin that had no current glacial effect provided a better means of dividing the seasonal hydrograph for Place Creek.  4.1.2 Synoptic scale During the first two discharge events of the season (20 to 23 M a y and 4 to 7 June respectively), S S C peaked at 47 and >71 m g / L in Place Creek. The third event, on 18 to 19 June, produced discharge values similar to those for the second event, but S S C peaked at only 43 mg/L. Although streamflow during both the second and third events was similar, rainfall during the former was almost double that during the latter. It is possible, therefore, 66  that the latter was a rain-on-snow event because air temperature was higher on 18 to 19 June than on 5 June. Intuitively, it would seem that snowmelt runoff should carry less sediment to the channel than does storm runoff because the snowpack acts as a protective layer, thus reducing runoff access to sediment sources, whilst raindrops provide an impact mechanism to loosen soil particles. Further evidence of this claim is provided by the fourth, snowmelt, event (28 to 30 June), during which little rain fell. Again, peak discharge during this event was similar to that reached during the second and third events. Assuming from the minimal rainfall that the bulk of discharge was induced by snowmelt runoff, the associated peak recorded SSC of 23 mg/L suggests that snowmelt runoff only carried a relatively minor sediment load to the stream channel.  Figure 3.14 also showed that snowmelt events  transported a smaller percentage of organic sediment as the nival season progressed. This finding supports the speculation made in Sub-section 4.1.1 that snowmelt runoff entrains less sediment towards the end of the nival period because of the rising elevation of snowmelt to areas where soils and vegetation cover are thinner.  Another explanation for the progressive decrease in event-discharge SSC from the second to fourth events would be synoptic scale sediment supply exhaustion, in which, for example, most of the small grains, loosened by freeze-thaw cycles during the winter, or deposited during the previous autumn, were washed into the stream by the first and second events.  Additionally, on daily timescales within the fourth synoptic period, increasing  diurnal snowmelt discharge cycles were associated with progressively smaller corresponding diurnal SSC cycles. Such a negative relation between SSC and discharge could be evidence of intra-event sediment exhaustion.  67  4.1.3 Diurnal scale - glacial period In order to discuss the hydrologic control on S S C on a diurnal timescale, focus will be made on the glacial period when predominantly fine-grained material, the product of glacial erosion, was transported down Place Creek.  Figure 3.15a (4 to 6 Aug) suggests that there  was a constant supply of fine sediment evacuated from Place Lake during diurnal glacial melt cycles. Since S S C in the valley mimicked the pattern of discharge, it is unlikely that inputs of water between the two sites induced the diurnal fluctuation of S S C in early August. Assuming no significant measurement errors, a more plausible explanation is that short-term storage and re-mobilisation of fine sediment between Place Lake and the valley bottom occurred during the 24 hr cycle.  O n the falling hydrograph limb, large within-channel  boulders may have provided a surface onto which clay and silts were deposited.  During the  rising hydrograph limb, some of this fine sediment could have been re-mobilised as the water level rose up around the boulders.  Estimates made in Section 3.6 revealed that between 5  A u g 12:00 and 6 A u g 11:59, approximately 12,000 kg of sediment left Place Lake, but only 9000 kg passed by the valley-bottom site.  Thus the value of -3000 ± 2 0 0 kg for the  approximate daily suspended sediment budget at the mouth of the basin may provide evidence for a net loss of sediment to channel storage during diurnal glacier melt cycles.  It  must be noted, however, that the sediment budget is an estimate and to validate the difference in upstream-downstream fluxes would require lake outlet discharge to be measured.  A t the end of August (Figure 3.15b), suspended sediment supply from the Place Lake was complicated by rainfall runoff, as evidenced by the undulating S S C series at the lake outlet,  which was probably induced by periodic inputs of sediment  from the slopes  surrounding the lake. The only time at which valley S S C approached values measured at the  68  lake outlet was at the glacial melt discharge peak on 30 A u g . It can be deduced, therefore, that either transient stores of fine glacial sediment were not re-mobilised during the storm event, or that the rainfall component of discharge, generated downstream of the lake outlet, had a relatively low S S C and thus served to 'dilute' the glacial component of total S S C . The former is possible because streamflow was approximately 1 m /s less in the late August measurement period than in early August. Thus maximum discharge on 29 A u g would not have risen to heights on within-channel boulders that minimum streamflow reached on 5 Aug.  The latter is plausible because at the synoptic scale it was shown in Section 4.1.2 that  storm event sediment supplies were progressively exhausted from event to event, unless discharge rose beyond values previously reached in the season (i.e., during the largest event of the season, on 27 to 28 July).  4.1.4  Integrated  temporal  scales  The separate S S C - Q rating relations for Place Creek (Figure 3.11) indicated that seasonal suspended sediment transport is principally controlled by two sources of streamflow, glacial and nival, combined with a nival-glacial mixture and rainfall runoff. The points within the three sub-seasons plot in different areas of the S S C - Q plot. Initially, the slope of the relation was steep during the nival period as early season discharge events (both rainfall and snowmelt derived) flushed out sediment that had accumulated over the preceding winter. In the nival-glacial period, the slope decreased because only sand-sized material was being transported and large discharge events did not entrain as much sediment as earlier in the season. Finally, during the glacial period, the average suspended sediment particle size likely decreased from sand- to silt- and clay-sized material. Accordingly, the intercept of the S S C Q rating function increased as more sediment could be transported at relatively lower  69  discharges, than over the nival-glacial transition period.  Thus from a seasonal to sub-  seasonal scale, the groups of points in Figure 3.11 moved in an anticlockwise direction, producing anticlockwise hysteresis.  At the diurnal scale, Sub-section 3.5.4 highlighted different types and strengths of diurnal scale hysteresis. Changes in diurnal hysteresis imply that different sediment and water delivery processes occur through the season. For storm runoff events, the strength and clockwise direction of hysteresis during the early season suggest that the bulk of suspended sediment was flushed through the basin during the rising limb of storm hydrographs (Figure 3.12a). Later in the season, however, the loops were both tighter and anticlockwise during rain events (Figure 3.12b), indicating that the sediment did not reach the channel until the falling hydrograph limb. The seasonal change in the character of storm-hysteresis probably reflects the decrease in supply of near-bank sediment, which would have been washed into the channel during early season storm events. Thus sediment that is supplied to the channel during late season storms may take a longer time to be transported to the channel by rainfall runoff.  During the nival period, complex hysteresis was noted in Figure 3.12d, in which the latter part of the loop was vertical and anticlockwise.  This irregular loop illustrates the  importance of preceding events, which in this case was a rain event. O n the falling limb of the storm hydrograph, maximum daily air temperature rose from 2 0 ° C on 2 July to 2 8 ° C on 3 July.  It is plausible, therefore, that although discharge was declining, the storm runoff  proportion of streamflow was also declining and being replaced by nival-glacial meltwater which may have initially been carrying higher values of S S C . Another interesting loop on 21-22 July (Figure 3.12e) depicts a 'figure 8' formation. The 'figure 8' loop in Figure 3.12e 70  combined initially clockwise and then anticlockwise sub-loops. The former occurred during higher and the latter at lower discharge, implying that following peak flow,  sediment  availability was great enough that S S C decreased more slowly than discharge. This figure 8 loop also occurred at the end of the nival-glacial season and so again it may be that diurnal changes in hydrologic sources, from nival to glacial, are reflected in the pattern of S S C - Q hysteresis.  During the glacial season, clockwise hysteresis at both diurnal and synoptic scales (Figure 3.13) provided evidence of sediment flushing prior to daily and synoptic peaks in discharge. It is probable that such flushing is of sediment stored within the channel. This sediment may have been derived from its deposition during falling diurnal hydrograph limbs on a daily timescale and, on a synoptic scale, from its deposition during the receding weekly hydrograph limb following high flows during early August.  The integration of diurnal, synoptic and seasonal timescales  is of fundamental  importance in the attempt to understand the proglacial suspended sediment transport system in Place Creek. It has been argued in this thesis that different water sources exert different hydrologic controls on downstream values of S S C in Place Creek.  However, to predict  accurately how S S C will respond to different hydrologic conditions, effects of processes at different timescales must be integrated.  The aim of the next section is, therefore, to show  how the M L R models presented in Sub-section 3.4.3 incorporate the ideas of integrated timescales presented in this Sub-section and to analyse the process inferences that can be made from the significant variables in each sub-season M L R model.  71  4.2 Modelling SSC Figure 4.1 summarises the stages followed in this study in order to obtain predictions of S S C , firstly using the C U S U M S technique to sub-divide the seasonal discharge series and then using the M L R models in Table 3.5.  Temporal scales are shown above each stage and the  final predicted value of S S C represents the integration of the timescales within which each stage operates.  Within the sub-seasonal stage, the nival-glacial and then glacial periods are  indirectly influenced by the legacies of previous sub-seasons within the entire hydrological season. The diurnal and synoptic scales are combined because each of the predictor variables operates within both scales.  Considerable improvement in the accuracy of the sub-season M L R models, over the sub-season O L S models (Table 3.4), is shown by the high adjusted R 3.5).  2  (R dj) values (Table 2  a  The R dj for the nival period may be slightly higher because n was smaller and only 2  a  spanned the period 18 M a y to 17 June. The prevailing hydrologic conditions and sources of suspended sediment within this time were likely to be more similar than for the entire nival period. For the nival-glacial period, however, the hydrology of the basin was transitional and for the glacial season the period encompassed over two months.  During these latter sub-  seasons, therefore, it could be expected that model accuracy would be lower than for the nival season because as time increases the stochastic element in the data set may also increase.  Such stochastic model elements can, in part, be attributed to short-term S S C  autocorrelation and be accounted for by adding a lagged S S C predictor variable to each of the models in Table 3.5 (Hodson and Ferguson, 1999). Having added S S C , lagged by sample interval, as an additional predictor, R dj values improved by up to 0.1. 2  a  72  The M L R models  SEASONAL  SUB-SEASONAL  !  SYNOPTIC / DIURNAL  NIVAL Q Q IQ 2  Q  w  NIVAL-GLACIAL  TA  AQ  3  SSC(ref)  SSC  GLACIAL  Figure 4.1: Flow chart model for prediction of S S C from sub-seasonal M L R models.  73  would, however, become redundant for cost-effective S S C prediction because a closeinterval sampling scheme would be required to calculate values of such lagged S S C variables.  In all of the sub-season M L R models, Q _i was not significant because its contribution t  to the model is provided through the combination of the other significant predictor variables. Its insignificance in, for example, the nival season may reflect low mean S S C (9.9 mg/L) for the nival period (Table 3.2) and that when S S C was highest (during rain events), most of the variability in S S C was explained by P24.  For the nival period, the predictors are (in order of decreasing significance) Q , P24, 2  A Q 3 , E Q and SSC( f>. A l l these variables have positive coefficients, indicating that as they re  increase they cause a collective increase in S S C . From what was reported earlier in this section, it appears that for the nival period, E Q does not, therefore, represent a decreasing supply of sediment.  Rather, its positive coefficient possibly reflects the relatively high  values of S S C recorded during the rainfall events over the nival period rather than an increasing area of suspended sediment supply. included in the model because it was  Even though it is insignificant, SSC( f) is re  'forced' in through the selection process.  Its  contribution to the overall model may be relatively small because mean S S C for the nival season was low compared to the nival-glacial and glacial seasons (Table 3.2).  The M L R model for the nival-glacial period incorporates SSC( f), A Q 3 , Q , E Q , T A , 2  re  P24 and Q in order of decreasing significance.  For this model, the increasing importance of  SSC( f) is reflected in the increased variability of S S C through the nival-glacial season and re  thus the need for a reference value. A s glacial melt became a more established component of  74  stream discharge during July, A Q 3 also became more important because the major portion of daily suspended sediment transfer occurred on the rising hydrograph limb in Place Creek. The coefficient of P24 was negative because rainy conditions correlated with subdued snow and glacial melt production. Discharge and T became significant, but also had negative A  coefficients, possibly because the individual effects of the variables are confounded by the other predictors.  The negative coefficient of E Q implies that seasonal sediment supply  exhaustion played a significant role in the temporal control on SSC through the nival-glacial period.  During the glacial sub-season, the M L R model contains Q, SSC f), EQ, A Q 3 , P24 and (re  T A in order of decreasing significance. Air temperature is the only variable with a negative coefficient and this is likely to be a function of the passage of storm systems with associated pulses of high SSC. A positive coefficient is produced for EQ, indicating that during the ablation season, the supply of suspended sediment (dominantly glacially derived) was not exhausted as commonly noted in the literature. As for the nival-glacial model, SSC( f) and re  AQ3 were significant. The former shows that the hydrologic control on proglacial suspended sediment dynamics changed on short timescales.  The latter again provides evidence that  most sediment is transported during the rising limb of the diurnal glacial melt hydrograph.  Overall, the M L R models in Table 3.5 can be relied upon to give accurate predictions of SSC because the PRESS statistics were low (Table 3.6). assumptions are generally met for these models.  Additionally, the regression  The need for log transformations and e  associated problems of bias is thus avoided. Although the glacial M L R model did not meet the regression assumptions, most M L R SSC prediction models proposed in the literature  75  acknowledge that such assumptions are difficult to meet, given white noise in the data set (e.g. Hodgkins, 1999).  Frequently, measurements are discarded in the literature without  adequate reason, besides being outliers.  In this study, no data were discarded because no  obvious measurement error was noted for any of the samples. F r o m discussions of transient sediment storage in Section 4.1, it may be that during the glacial period, the relation between S S C and its suite of glacial period predictors is complicated by short term pulses of sediment derived from within channel stores such as the surface of large boulders.  76  Chapter 5 Conclusions 5.1 Summary of key findings 5.1.1  SSC-Q  relations  in a glacier-fed  stream  Place Creek is dominated by a nival- and glacial-melt discharge regime, which is punctuated by storm runoff events. A simple S S C - Q rating relation was shown in Sub-section 3.5.1 to be an inaccurate S S C prediction tool.  In order to investigate the hypothesis that different  sources of water exert different controls on the S S C - Q relation, the seasonal Place Creek discharge series was divided into sub-seasons.  Corresponding sub-seasonal O L S rating  curves showed great improvement. Coupled with E C data that supported the hypothesis that streamflow within the sub-seasons was derived from different sources (snowmelt, glaciermelt and groundwater), it can be concluded that a fundamental source of explanation for the scatter in the seasonal S S C - Q relation is provided by the transition from nival to nival-glacial and finally to glacial sources of runoff, combined with periodic inputs of storm runoff through the season.  5.1.2  Multiple  regression  models  The use of M L R models to predict S S C from a suite of variables has arguably been extended in the literature to the point at which they are no longer time or cost efficient to evaluate over subsequent seasons.  The simple sub-seasonal M L R models presented in Sub-section 3.5.3  can be evaluated from easily measured predictors and thus provide a cost efficient tool with which to predict accurately S S C in Place Creek over future seasons.  The M L R models  themselves represent an integration of timescales, from seasonal to diurnal. Their improved accuracy from  the  sub-seasonal  O L S models  77  indicates  that they  provide additional  explanation for the scatter in the S S C - Q relation that cannot be accounted for at the subseasonal scale.  The detailed data set that was used to derive the M L R model parameters in  this thesis represents an extensive and accurate record of S S C , something that has rarely been achieved in the literature.  5.1.3  Within-channel transient sediment storage  Place Creek appears to be an example of a sediment-transfer system that is dominated by a transient diurnal storage mechanism within the ablation season, possibly provided by the surface of large within-channel boulders. Such storage has not been widely acknowledged in the literature, largely because of a lack of intensive studies of upstream-downstream changes in S S C in streams with minimal tributary or groundwater inflow. be a regulating sink for fine glacially-derived sediment,  Place Lake was shown to  while diurnal fluctuations  of  downstream S S C and a negative suspended sediment budget at the basin mouth suggested that the sediment was deposited and re-mobilised within the channel.  These data are  compatible with the boulder storage hypothesis.  5.2 Future research Clearly this thesis presents a number of claims that require further evidence to validate. Future research projects could (a) generate daily sediment budgets for the entire glaciofluvial lacustrine system; (b) investigate the trapping efficiency of Place Lake and how it varies with discharge; (c) undertake a more extensive upstream-downstream monitoring campaign of S S C and discharge; (d) analyse the caliber of the suspended sediment transported down Place Creek in order to identify spatial and temporal changes in sediment-source areas; (e) collect geochemical  data from water samples  for diurnal, synoptic  and seasonal hydrograph  separation with the aim of providing more detailed divisions of water sources that affect the  78  SSC-Q relation, particularly to investigate the hypothesis that the changes in the EC-Q relation reflect seasonal changes in the location of water sources.  79  References Aizen V B , Aizen E M , Melack J M . 1996. Precipitation, melt and runoff in the northern Tien  Shan. Journal of Hydrology 186: 229-251. Asselman N E M . 2000. Fitting and interpretation of sediment rating curves. Journal of  Hydrology 234: 228-248. Barsch D and Caine N . 1984. The nature of mountain geomorphology.  Mountain Research  and Development 4(4): 287-298. Benn D I , Evans J A . 1998.  Glaciers and Glaciation. London: Arnold: 734 pp.  Bogen J. 1980. The hysteresis effect of sediment transport systems.  Norsk Geogr. Tidsskr 34:  45-54. Brown G H , Sharp M J , Tranter M , Gurnell A M and Nienow P W . 1994. impact of postmixing chemical reactions on the major ion chemistry of bulk meltwaters draining the Haut Glacier D ' A r o l l a , Valais, Switzerland.  Hydrological Processes 8: 465-480.  Brown G H , Tranter M and Sharp M J . 1996. Experimental investigations of the weathering of suspended sediment by alpine glacial meltwater.  Hydrological Processes 10: 579-  597. Chikita K . 1993. A comparative study on suspended sediment discharge initiated by snow- or glacier-melting. Publication 217, International Association of Hydrological Sciences: Wallingford; 257-267. Collins D N . 1979. Sediment concentration in meltwaters as an indicator of erosion processes beneath an Alpine glacier.  Journal of Glaciology 23: 247-257.  Collins D N . 1998. Suspended sediment flux in meltwaters draining from Batura glacier as an indicator of the rate of glacial erosion in the Karakoram mountains.  Quaternary  Proceedings 6: 1-10. Elder K , Kattelmann R, Ferguson R. 1990. Refinements in dilution gauging for mountain streams.  Publication  193,  International  Association  of  Hydrological  Sciences:  Wallingford: 247-254.  Farley A L . 1979. Atlas of British Columbia. Vancouver: U B C Press: 136pp. Fenn C R . 1989. Quantifying the errors involved in transferring suspended sediment rating equations across ablation seasons. Annals  of Glaciology 13: 64-68.  Fenn C R , Gurnell A M , Beecroft IR. 1985. A n evaluation of the use of suspended sediment rating curves for the prediction of suspended sediment concentration in a proglacial stream.  Geografiska Annaler 61A{\-2): 71-82. 80  Ferguson R L 1986. River loads underestimated by rating curves. Water Resources Research 22: 74-76. Ferguson RI 1987. Accuracy and precision of methods for estimating river loads. Earth  Surface Processes and Landforms 12: 95-104. Gurnell A M . 1982. The dynamics of suspended sediment concentration in an Alpine proglacial stream network. Publication 138, International Association of Hydrological Sciences: Wallingford: 319-330. Gurnell A M . 1987. Suspended sediment. In Glacio-fluvial Sediment Transfer: an Alpine  Perspective, Gurnell A M , Clark M J (eds). Wiley: Chichester: 305-353. Gurnell A M , Fenn C R . 1984. Box-Jenkins transfer function models applied to suspended sediment concentration-discharge relationships in a proglacial stream. Arctic and Alpine Research 16: 93-106. Gurnell A M , Clark M J , H i l l C T . 1992a. Analysis and interpretation of patterns within and between hydroclimatological time series in an alpine glacier basin. Earth Surface  Processes and Landforms 17: 821-839. Gurnell A M , Clark M J , Hill C T and Greenhalgh J. 1992b. Reliability and representativeness of a suspended sediment concentration monitoring programme for a remote alpine proglacial river. International Association of Hydrological Sciences Publications 210: 191-200. Gurnell A M , Hodson A , Clark M J , Bogen J , Hagen J O , Tranter M . 1994. Water and sediment flux from glacier basins: an arctic and alpine comparison. Publication 224, International Association of Hydrological Sciences: Wallingford: 325-334. Gurnell A M , Edwards PJ, Petts G E , Ward J V . 1999. A conceptual model for alpine proglacial channel evolution under changing climatic conditions. Catena 38: 223-242. Hallet B , Hunter L , Bogen J. 1996. Rates of erosion and sediment evacuation by glaciers: a review of field data and their implications. Global and Planetary Change 12: 213235. Hammer K M , Smith N D . 1983. Sediment production and transport in a proglacial stream: Hilda Glacier, Alberta, Canada. Boreas 12: 91-106. Hannah D M , McGregor G R . 1997. Evaluating the impact of climate on snow- and ice-melt dynamics in the Taillon Basin, French Pyrenees. Journal of Glaciology 43(145): 563568. Hannah D M , Gurnell A M . 2001. A conceptual, linear reservoir runoff model to investigate melt season changes in cirque glacier hydrology. Journal of Hydrology 246: 123-141.  81  Hannah D M , Smith B P G , Gurnell A M , McGregor G R . 2000. A n approach to hydrograph  classification. Hydrological Processes 14: 317-338. Harbor J , Warburton J. 1993. Relative rates of glacial and non-glacial erosion in alpine  environments. Arctic and Alpine Research 25(1): 1-7. Hicks D M , McSaveney M J , Chinn T J H . 1990. Sedimentation in proglacial Ivory lake, Southern Alps, New Zealand. Arctic and Alpine Research 22(1): 26-42. Hodgkins R. 1996. Seasonal trend in suspended-sediment transport from an Arctic glacier, and implications for drainage-system structure. Annals ofGlaciology 22: 147-151. Hodgkins R. 1999. Controls on suspended-sediment transfer at a high-arctic glacier, determined from statistical modelling. Earth Surface Processes and Landforms 24: 121. Hodson A , Gurnell A , Tranter M , Bogen J , Hagen J O , Clark M . 1998. Suspended sediment yield and transfer processes in a small High-Arctic glacier basin, Svalbard. Hydrological Processes 12: 73-86. Hodson A J , Ferguson R L 1999. Fluvial suspended sediment transport from cold and warmbased glaciers in Svalbard. Earth Surface Processes and Landforms 24: 957-974. Jansson M . 1985. A comparison of detransformed logarithmic regressions function regressions. Geografiska Annaler 67A(l-2): 61-70.  and power  Laudon H . 1995. Spatial and temporal variability of the stream water chemistry of an alpine/sub-alpine catchment, in the Coast Mountains of British Columbia, Unpublished MSc. Thesis, Department of Geography, University of British Columbia. Lenzi M A , Marchi L . 2000. Suspended sediment load during floods in a small stream of the Dolomites (northeastern Italy). Catena 39: 267-282. Mac Nally R, Hart B T . 1997. Use of C U S U M methods for water-quality monitoring in storages. Environmental Science and Technology 31(7): 2114-2119. Matthews W H . 1956. Physical limnology and sedimentation in a glacial lake. Bulletin of the  Geological Society of America 67: 537-552. Miller D M . 1984. Reducing transformation bias in curve fitting. The American Statistician 38(2): 124-126. Milliman J D , Syvitski J P M . 1992. Geomorphic/tectonic control of sediment discharge to the Ocean: the importance of small mountainous rivers. Journal of Geology 100: 522544. Moore R D . 1993. Application of a conceptual streamflow model in a glacierized drainage basin. Journal of Hydrology 150: 151-168.  82  Moore R D and McKendry IG. 1996. Spring snowpack anomaly patterns and winter climate variability, British Columbia, Canada.  Water Resources Research 32(3): 623-632.  Neter J , Wasserman M , Kutner W H , Nachtsheim C J . 1996.  Applied Linear Statistical  Models. Chicago: McGraw-Hill: 1408pp.  Glaciofluvial and Glaciolacustrine Sedimentation, Jopling A V , McDonald B C (eds). 1975. Society of  0strem G . 1975. Sediment transport in glacial meltwater streams. In  Economic Palaeontologists and Mineralogists Special Publication 23: 101-122. Picouet C , Hingray B , Olivry  J C . 2001. Empirical and conceptual modelling of the  suspended sediment dynamics in a large tropical African river: the Upper Niger river  basin. Journal of Hydrology 250(1-4): 19-39. Richards K S . 1984. Some observations on suspended sediment dynamics in Storbregrova,  Jotunheim. Earth Surface Processes and Landforms 9: 101-112. Richardson J , Rowe D K , Smith JP. 2001. Effects of turbidity on the migration of juvenile banded kokopu (Galaxias fasciatus) in a natural stream. The New Zealand Journal of  Marine and Freshwater Research 35(1): 191-196. Rothlisberger H , Lang H . 1987. Glacial Hydrology. In Gurnall A M , Clark M J (eds), Glaciofluvial Sediment Transfer. New York: Wiley: 207-284. Ryan P A . 1991. Environmental-effects of sediment on New-Zealand streams - a review. The  New Zealand Journal of Marine and Freshwater Research 25(2): 207-221. Ryder J M , Thomson B . 1986. Neoglaciation in the southern Coast Mountains of British Columbia: chronology prior to the late Neoglacial maximum.  Canadian Journal of  Earth Sciences 23(3): 273-287. Singh P, Huebl H , Weinmeister H W . 2000. Use of recession characteristics of snowmelt hydrographs in the assessment of snow water storage in a basin.  Hydrological  Processes 14: 91-101. Stott T A , Grove JR. 2001. Short-term discharge and suspended sediment fluctuations in the proglacial Skeldal River, north-east Greenland.  Hydrological Processes 15: 407-423.  Syvitski JP, Morehead M D , Bahr D B , Mulder T . 2000. Estimating fluvial sediment transport: The rating parameters.  Water Resources Research 36(9): 2747-2760.  Teti P A . 1979. The variability of stream chemistry in a Coast Mountain watershed, British Columbia,  Unpublished MSc.Thesis, Department of Geography, University of British  Columbia. Walling D E . 1977. Assessing the accuracy of suspended sediment rating curves for a small basin.  Water Resources Research 13(3): 531-538.  83  Weirich F H . 1985. Sediment budget for a high energy glacial lake.  Geografiska Annaler 67A  (1-2): 83-99. Williams G P . 1989. Sediment concentration versus discharge during single hydrologic events in rivers. Journal  of Hydrology 111: 89-106.  Willis IC, Richards K S , Sharp M J . 1996. Links between proglacial stream suspended sediment dynamics, glacier hydrology and glacier motion at Midtdalsbreen, Norway.  Hydrological Processes 10: 629-648. Wren D G , Barkdoll B D , Kuhnle R A , Derrow R W . 2000. Field techniques for suspendedsediment  measurement.  Journal of Hydraulic Engineering 126(2): 97-104.  84  Appendix A Dummy c o d i n g used f o r c a t e g o r i c a l v a r i a b l e s i n model. C a t e g o r i c a l v a l u e s encountered d u r i n g p r o c e s s i n g a r e : SEASON (3 l e v e l s ) 0, 1, 2 Step # 0 R = 0.805 R-Square = 0.648 Effect Coefficient Std Error S t d Coef In 1 2 3 4 Out  Constant Q SEASON SEAS0N*Q  15.673 . . Part. Corr.  0.722 . .  T o l . df  0.683 0.53046 . 0.08622 . 0.09959  F  'P'  1 471.126 2 23.614 2 21.165  0.000 0.000 0.000  none  Estimates of e f f e c t s SSC CONSTANT Q SEASON SEASON SEASON Q SEASON Q  Source Regression Residual  0 1 0  1.332 15.673 -19.43 6 -5.716 15.642  1 -3.811 A n a l y s i s of Variance Sum-of-Squares df Mean-Square 75547.344 40977.537  5 669  85  15109.469 61.252  F-ratio  P  246.677  0.000  


Citation Scheme:


Usage Statistics

Country Views Downloads
China 8 14
United States 7 0
France 3 0
Canada 1 0
Argentina 1 0
City Views Downloads
Beijing 5 0
Unknown 4 33
Shenzhen 3 14
Ashburn 3 0
Vancouver 1 0
Wilmington 1 0
Los Angeles 1 0
Sunnyvale 1 0
El Calafate 1 0

{[{ mDataHeader[type] }]} {[{ month[type] }]} {[{ tData[type] }]}
Download Stats



Customize your widget with the following options, then copy and paste the code below into the HTML of your page to embed this item in your website.
                            <div id="ubcOpenCollectionsWidgetDisplay">
                            <script id="ubcOpenCollectionsWidget"
                            async >
IIIF logo Our image viewer uses the IIIF 2.0 standard. To load this item in other compatible viewers, use this url:


Related Items