Channel morphology and bed-load yield in fluvial, formerly-glaciated headwater streams of the Columbia Mountains, Canada

Abstract

This study examines channel-reach morphology and bedload yield dynamics in relation to landscape structure and snowmelt hydrology in headwater streams of the Columbia Mountains, Canada. Data collection relies on field surveys and geographic information systems analysis in conjunction with a nested monitoring network of water discharge and bedload transfer. The landscape is characterized by subdued, formerly-glaciated upland topography in which the geomorphic significance of landslides and debris flows is negligible and fluvial processes prevail. While the spatial organization of channel morphology is chiefly controlled by glacially imposed local slope in conjunction with wood abundance and availability of glacigenic deposits, downstream patterns of the coarse grain-size fraction, bankfull width, bankfull depth, and stream power are all insensitive to systematic changes of local slope along the typically stepped long profiles. This is an indication that these alluvial systems have adjusted to the contemporary snowmelt-driven water and sediment transport regimes, and as such are able to compensate for the glacially-imposed boundary conditions. Bedload specific yield increases with drainage area suggesting that fluvial re-mobilization of glacial and paraglacial deposits dominate the sedimentary dynamics of basins as small as 2km2. Stepwise multiple regression analysis shows that annual rates of sediment transfer are mainly controlled by the number of peak events over threshold discharge. During such events, repeated destabilization of channel bed armoring and re-mobilization of sediment temporarily stored behind LWD structures can generate bedload transport across the entire snowmelt season. In particular, channel morphology controls the variability of bedload response to hydrologic forcing. In the present case studies, we show that the observed spatial variability in annual bedload yield appears to be modulated by inter-basin differences in morphometric characteristics, among which slope aspect plays a critical part.