We compared median runoff (R) and precipitation (P) relationships over 25 years from 20 mesoscale (50 to 5,000 km2) catchments on the Boreal Plains, Alberta, Canada, to understand controls on water sink and source dynamics in water‐limited, low‐relief northern environments. Long‐term catchment R and runoff efficiency (RP−1) were low and varied spatially by over an order of magnitude (3 to 119 mm/year, 1 to 27%). Intercatchment differences were not associated with small variations in climate. The partitioning of P into evapotranspiration (ET) and R instead reflected the interplay between underlying glacial deposit texture, overlying soil‐vegetation land cover, and regional slope. Correlation and principal component analyses results show that peatland‐swamp wetlands were the major source areas of water. The lowest estimates of median annual catchment ET (321 to 395 mm) and greatest R (60 to 119 mm, 13 to 27% of P) were observed in low‐relief, peatland‐swamp dominated catchments, within both fine‐textured clay‐plain and coarse‐textured glacial deposits. In contrast, open‐water wetlands and deciduous‐mixedwood forest land covers acted as water sinks, and less catchment R was observed with increases in proportional coverage of these land covers. In catchments dominated by hummocky moraines, long‐term runoff was restricted to 10 mm/year, or 2% of P. This reflects the poor surface‐drainage networks and slightly greater regional slope of the fine‐textured glacial deposit, coupled with the large soil‐water and depression storage and higher actual ET of associated shallow open‐water marsh wetland and deciduous‐forest land covers. This intercatchment study enhances current conceptual frameworks for predicting water yield in the Boreal Plains based on the sink and source functions of glacial landforms and soil‐vegetation land covers. It offers the capability within this hydro‐geoclimatic region to design reclaimed catchments with desired hydrological functionality and associated tolerances to climate or land‐use changes and inform land management decisions based on effective catchment‐scale conceptual understanding.

Abstract. Preferential flow paths have been found to be important for runoff generation, solute transport, and slope stability in many areas around the world. Although many studies have identified the particular characteristics of individual features and measured the runoff generation and solute transport within hillslopes, very few studies have determined how individual features are hydraulically connected at a hillslope scale. In this study, we used dye staining and excavation to determine the morphology and spatial pattern of a preferential flow network over a large scale (30 m). We explore the feasibility of extending small-scale dye staining techniques to the hillslope scale. We determine the lateral preferential flow paths that are active during the steady-state flow conditions and their interaction with the surrounding soil matrix. We also calculate the velocities of the flow through each cross-section of the hillslope and compare them to hillslope scale applied tracer measurements. Finally, we investigate the relationship between the contributing area and the characteristics of the preferential flow paths. The experiment revealed that larger contributing areas coincided with highly developed and hydraulically connected preferential flow paths that had flow with little interaction with the surrounding soil matrix. We found evidence of subsurface erosion and deposition of soil and organic material laterally and vertically within the soil. These results are important because they add to the understanding of the runoff generation, solute transport, and slope stability of preferential flow-dominated hillslopes.

[1] Our understanding of hillslope subsurface flow relies on assumptions about how storm characteristics affect the hillslope runoff response. Experiments in hillslopes dominated by preferential flow features often show that runoff is dynamic and is affected by antecedent conditions, rainfall conditions, and position of the slope. We applied tracers to a hillslope under natural and steady state flow boundary conditions to determine the relationship between lateral tracer velocities and various hillslope lengths and storm indicators. Tracer velocities were similar to the fastest velocities measured in other similar experiments. The velocities were dependent on the boundary conditions and slope length, and the subsurface flow velocity was most closely related to the 1-h rainfall intensity. Unlike some studies, there was little correlation between our measured flow velocities and storm volume or antecedent conditions. We attributed this to the hillslope characteristics and the relatively consistent wet antecedent conditions during the experiments. This experiment showed that the connectivity of the hillslope preferential flow network is an important factor governing the average subsurface flow velocity.

Abstract. Preferential flow features have been found to be important for runoff generation, solute transport, and slope stability in many areas around the world. Although many studies have identified the particular characteristics of individual features and measured the runoff generation and solute transport within hillslopes, no studies have determined how individual features are hydraulically connected at a hillslope scale. In this study, we used dye staining and excavation to determine the morphology and spatial pattern of a preferential flow network over a large scale (30 m). We explore the feasibility of extending small-scale dye staining techniques to the hillslope scale. We determine the lateral preferential flow features that are active during the steady state flow conditions and their interaction with the surrounding soil matrix. We also calculate the velocities of the flow through each cross-section of the hillslope and compare them to hillslope scale applied tracer measurements. Finally, we investigate the relationship between the contributing area and the characteristics of the preferential features. The experiment revealed that larger contributing areas coincided with highly developed and hydraulically connected preferential features that had flow with little interaction with the surrounding soil matrix. We found evidence of subsurface erosion and deposition of soil and organic material laterally and vertically within the soil. These results are important because they add to the understanding of the runoff generation, solute transport, and slope stability of these types of hillslopes.

Abstract:Characterizing zones of a watershed based on the water table is used to understand and predict internal watershed processes. In watersheds dominated by lateral preferential flow, the water table response typically shows a distinct hydraulically limited pattern. This response is characterized by a capping of the rising water table when the lateral preferential flow features are activated and subsurface flow still increases. We expected that this response would be related to the contributing area since nearby hillslope excavations showed that the development of preferential flow network was positively correlated with the contributing area. The watershed was stratified into three predetermined zones and installed 25 piezometers to measure the water table dynamics. The objectives were (1) to characterize the water table-runoff relationship, (2) to prove preferential flow by observable characteristics and (3) to test the feasibility of identifying areas within a watershed that are dominated by lateral preferential flow. Watershed zones were not well defined and there was no strong relationship between the hydraulically limited response and observable watershed characteristics. Although zones might still be useful for grouping the hillslope processes, the piezometric response may not be an appropriate indicator for mapping the watershed into areas with runoff dominated by lateral preferential flow.

Manually designing road networks for planning purposes is labour-intensive. As an alternative, we have developed a computer algorithm to generate road networks under a variety of assumptions related to road design standards. This method does not create an optimized road network, but rather mimics the procedure a professional might use when projecting roads by hand. Because many feasible road networks are possible, sensitivity analysis is required to choose the best ones. Such analysis gives forest planners additional information with which to assess the long-term consequences of road density and road standards common in forest management decisions. The procedures used to create road networks are presented in this paper, along with a sensitivity analysis of assumptions on total network length, percentage of landings connected, grades, and horizontal and vertical alignment for a case study. We also include a sensitivity analysis of spatial detail such as node density and link characteristics. Although the road network generation algorithm requires manipulation of many input parameters to create desired road networks, and variation between outputs is a concern, the method still offers considerable improvement over manual methods, especially for applications in strategic planning, and appears to be suitable for all types of topography and road standards.

No abstract

Summary Stream crossing structures are an increasingly prevalent anthropogenic feature on North American riverscapes, particularly in watersheds affected by industrial resource development in sensitive boreal environments. If improperly managed, stream crossings have the potential to alter fish habitat and impede fish movement. This study assessed instream habitat characteristics and fish communities from 33 culverted, bridged and reference streams in an industrialising region of the boreal forest in west‐central Alberta. Mixed‐effects modelling and multivariate analysis were used to determine impacts of stream crossings at three scales: whole‐stream scale, within‐stream scale and the interaction of scales. Instream habitat characteristics such as mean depth, water velocity, percent fines, turbidity, water temperature and dissolved oxygen showed significant between‐stream as well as within‐stream differences among stream crossings. The majority of fish species exhibited significantly lower densities (n m−2) in upstream habitats as compared to downstream habitats, including a significant reduction in Slimy Sculpin densities in culverted streams. Multivariate tests corroborated these results, showing that fish assemblages differ as a function of stream type. This study suggests industrial stream crossings influence abiotic habitat characteristics in freshwater ecosystems, restrict biotic connectivity and impact fish community structure at the whole‐stream and within‐stream scales. Alterations to stream ecosystems associated with stream crossings may be driving large‐scale changes in stream fish communities in the boreal forest. With expanded development expected in much of North America's boreal region, mitigation measures which limit impacts from stream crossings are needed to ensure proper ecosystem function in freshwater systems.