Younes Alila


Relevant Degree Programs


Graduate Student Supervision

Doctoral Student Supervision (Jan 2008 - Mar 2019)
An investigation into the Flow Duration Curve in eastern United States : environmental controls and prediction at ungauged basins (2018)

The Flow Duration Curve (FDC) is a probabilistic flow representation relevant to streamflow investigation and physical understanding given its use in wide range of hydrological and ecological applications. A regional study that investigates FDCs and their prediction at ungauged catchments is important to develop causal models and provide insights to solve issues of water resources planning and management of aquatic ecosystems and habitats. Hydrological modelling and model parametrization in gauged and ungauged catchments are fundamental steps preceding the regional investigation of flow response using FDC. By means of Sacramento rainfall-runoff model (SAC-SMA), in 73 catchments from the eastern United States, I investigated the effect of SAC-SMA a priori parameters in the constrained calibration of the model. This analysis revealed and discussed limitations of a priori parameters that are intensively used to facilitate model calibration and make predictions at ungauged basins (PUB). The PUB using parameter transfer within homogeneous regions of similar climate and flow characteristics outweighed in performance the a priori parameters. The FDC was advantageous in revealing the effect of lack of efficiency and bias. A parameter regionalization approach is more efficient for PUB than a priori parameters. The ultimate limitations of the within-region parameter transfer are recognized and discussed. The interaction between climate and landscape properties was central to develop the physical understanding of the FDC. The high precipitation variability does not necessarily lead to FDC of a steep slope. Characteristics of the catchment— equivalent to a precipitation filter— interplayed with the precipitation and affected FDC shapes. The analysis highlighted the role of soil infiltration rates in specific conditions of soil moisture storage capacity and predominant runoff generation mechanism. The study revealed that processes underlying FDC and the FDC shapes are by far more complex than being characterized in the wider literature using characteristics of landscape or climate. In conditions of humid climate and perennial flow, a meta-analysis utilizing a process-based investigation showed that mean monthly runoff FDC —readily available at ungauged catchments—predicts only FDC middle third (exceedance probabilities between 33% and 66%). The method is constrained by the value of flow variability (slope of the FDC).

View record

Modelling sediment dynamics at the basin scale : implications of changes in climate and hydrological regimes (2018)

Basin-wide sediment dynamics are closely linked to hydrological processes and landscape and therefore expected to be susceptible to climate change. Simulating sediment transport through large basins presents a challenging problem to modellers; the relationship between water flux and sediment load is complex and non-linear, and significant sediment generation can occur over small spatial and time scales. To date, most studies have employed lumped empirical models that predict annual load at the outlet of a study basin, but do not consider variability across the basin or sub-annually. In this study, we adapt a physically-based, distributed suspended sediment transport model for large-scale use. The sediment model is integrated into the Terrestrial Hydrology Model with Biochemistry (THMB) as a routine to make use of THMB’s dynamic water routing. The coupled model is applied to the 230,000 km² Fraser River Basin (FRB) in British Columbia, Canada using 1) historical hydrological input to test the model and 2) synthetic input derived from Intergovernmental Panel on Climate Change (IPCC) scenarios A1B, A2, and B1 to study potential impacts of climate change. In both cases the input data is provided by the Pacific Climate Impacts Consortium (PCIC) and comes from simulations using the Variable Infiltration Capacity (VIC) model. Simulation results using historical inputs are compared with observations at five stations using the coefficient of determination (R²), Nash-Sutcliffe coefficient of efficiency (NSE), and percent bias (PBIAS) metrics. Overall, simulated load values match well with observed values, with the monthly simulations at the station nearest the outlet scoring R² = 0.78, NSE = 0.77, and PBIAS = -20%. Simulation results using climate scenario-driven inputs are studied for potential future changes in sediment dynamics. Results re- veal a general shift in hillslope erosion and sediment yield towards larger values from autumn to spring, reduced summer values, and an overall annual increase, with hillslope generation growing 35-45% from baseline levels and yield at the basin’s outlet increasing 10-15%. These physically-based results offer unique insights into the impacts of climate change on sediment processes within a large basin and their potential implications.

View record

Forests, floods and channel processes : illuminating links between forest harvesting, the flood regime and channel response in snowmelt headwater streams (2014)

A meta-analysis of four snowmelt catchments with moderate harvest levels (30% to 40%) utilizing a frequency-based approach demonstrates, despite a century’s-worth of studies to the contrary, how harvesting increases the magnitude and frequency of all floods on record including the largest floods (R.I. = 1:50 yrs.) and how such effects increase unchecked with increasing return period as a consequence of changes to both the mean and standard deviation of the flood frequency distribution. Additionally, meta-analysis outcomes reveal up to three-fold increases in the number and duration of peak flows with return periods ranging from 0.6*Q₁.₅ to Q₁₀, which includes floods capable of mobilizing bedload and altering the form of gravel-bed streams. A frequency-based meta-analysis provides new insights concerning the physical processes governing the relation between forests and floods in snowmelt environments that were previously unrecognized using traditional chronological-pairing methods. The dominant process responsible for flood regime changes following harvesting is the increase in basin-average snowmelt rates that are amplified or mitigated by physical characteristics such as aspect distribution, elevation range, slope gradient and amount of alpine area. The outcomes of a high-resolution, nested-monitoring investigation of hydrologic and geomorphic controls on bedload mobility indicate that flow regime changes from harvesting in snowmelt streams can alter rates of bedload transport, a first-order determinant of channel form in fluvial systems. Regression analysis shows that annual sediment yield is controlled by the number of peaks-over-threshold discharge. During peak events, repeated destabilization of channel armor and re-mobilization of sediment temporarily stored behind LWD generates bedload transport across the entire snowmelt season. However, the potential for channel response to changes in the flow regime depends on characteristics of channel morphology. Study results indicate that patterns of bedload entrainment and mobility are influenced by flow resistance associated with channel form, grain size and LWD while the value of the critical dimensionless shear stress varies with channel gradient, relative and absolute grain effects and flow resistance. Differences in bed texture, gradient and channel form contributes to variations in rates and characteristics of bed mobility, hence differences in potential for alteration due to changes in the flow regime.

View record

Modeling deep groundwater flow through fractured bedrock in a mountainous headwater catchment using a coupled surface water - groundwater model, Okanagan Basin, British Columbia (2013)

Quantifying recharge to the mountain block from headwater catchments in snowmelt dominated upland mountainous regions is an important aspect of hydrologic studies. This study contributes to understanding of the interaction between surface water, soil water and deep groundwater flow in headwater catchments. A novel approach was developed for estimating the bedrock hydraulic conductivity of a regional-scale fractured bedrock aquifer using discrete fracture network (DFN) modeling. The methodology was tested in the mountainous Okanagan Basin, British Columbia, Canada. Discrete fractures were mapped in outcrops, and larger-scale fracture zones (corresponding to lineaments) were mapped from orthophotos and LANDSAT imagery. Outcrop fracture data were used to generate DFN models for estimating hydraulic conductivity for the fractured matrix (Km). The mountain block hydraulic conductivity (Kmb) was estimated using larger-scale DFN models. Simulated Km and Kmb values range from 10⁻⁸ to 10⁻⁷ m/s, are consistent with estimates from regional modeling studies, and are greatest in a N-S direction, coinciding with the main strike direction of Okanagan Valley Fault Zone. Kmb values also decrease away from the fault, consistent with the decrease in lineament density. Simulated hydraulic conductivity values also compare well with those estimated from pumping tests. The estimates of Kmb were then used to represent the deep bedrock in a coupled surface water - groundwater model using MIKE SHE for the Upper Penticton Creek 241 headwater catchment in the Okanagan Basin. Although highly uncertain due to parameter uncertainty and calibration error, recharge to deep groundwater was ~4% of the annual water budget. An specified outward flux from the catchment boundary, representing ~6% of annual water budget, did not significantly impact streamflow calibration, indicating that such deep groundwater losses from the catchment can be accommodated in a model. This outflow may contribute to cross-catchment flow and, ultimately, to groundwater inflow to lower elevation catchments in the mountain block. The modeling exercise is one of the first in catchment hydrology modeling within steep mountainous terrain in which the lower boundary of the model is not treated as impermeable, and in which recharge to the deep bedrock and discharge to the surrounding mountain block were estimated.

View record

Bankfull and effective discharge in small mountain streams of British Columbia (2012)

In this thesis, the concept of channel-forming discharge developed for large lowland rivers is critically re-evaluated for small mountain streams. In large lowland rivers under equilibrium conditions, the effective discharge (the discharge interval that transports the greatest proportion of sediment) approximates the bankfull discharge. The effective and bankfull discharges therefore provide measurable analogues for the theoretical channel-forming discharge, responsible for the form and dimensions of the stream channel. In small mountain streams this concordance of bankfull and effective discharges has been suspected to break down, in part because the channel form and dimensions are determined by the non-fluvial sediment supply as well as by fluvial processes. This interaction of non-fluvial sediment supply with fluvial processes makes the application of models and concepts developed for lowland rivers to mountain streams difficult. The history of past glaciation in mountain environments adds complexity in the form of increased sediment supply and changes in streamflow regime over time.Bankfull and effective discharges are measured and estimated for small mountain streams in three distinct hydroclimatic zones in southern British Columbia. These discharges are determined using surveyed stream cross-sections, photogrammetric analysis of in-channel sediment, long-term Water Survey of Canada gauging records, and a two-fraction sediment transport model.Results indicate multiple orders of magnitude of variability in the frequency and magnitude of bankfull and effective discharge for the studied streams, and little correlation between the bankfull and effective discharges. Bankfull frequency can be related to hydraulic and hydroclimatic variables while effective discharge frequency cannot. Most of the studied streams are incised, with bankfull discharges one or more orders of magnitude greater than their effective discharge, a condition attributable to the legacy of Pleistocene glaciation. In these streams the effective discharge is not a channel-forming discharge; at best, it is a channel-maintaining flow, and at worst it is geomorphically meaningless. In a few streams, the bankfull and effective discharges are both large, rare events; these are the only streams in the study in which the bankfull and effective discharges approximate a truly channel-forming discharge.

View record

Master's Student Supervision (2010-2017)
Revealing forest harvesting effects on large peakflows in rain-on-snow environment with new stochastic physics (2017)

The full abstract for this thesis is available in the body of the thesis, and will be available when the embargo expires.

View record

Hydrogeomorphic disturbance, landscape development and riparian vegetation dynamics of an alluvial, temperate rainforest in the Carmanah River valley, British Columbia, Canada (2012)

The alluvial forest of the Carmanah River valley on the west coast of Vancouver Island, British Columbia, was studied to examine the role of hydrogeomorphic disturbance in perpetuating the shifting-mosaic of habitats within this diverse ecosystem. Field-based research was complemented by a landscape-scale analysis that examined changes in the extent of specific forest types using a 70-year aerial photographic record. Thirty-eight plots containing 4509 trees were sampled for forest structure, composition, age, understory composition, and elevation above the contemporary channel. These field data, including a vegetation chronosequence spanning over 500 years, were used to examine vegetation dynamics. Over the past century, Carmanah River has eroded nearly 30% of the alluvial forest in this study area – 65% over the past 500 years. High magnitude floods result in diminished floodplain forest area by converting forests to channel. This results in a subsequent course of vegetation succession and geomorphic development. Fluvial deposits are colonized by a high density of Alnus rubra accompanied by a subcanopy of Picea sitchensis individuals. As Alnus die off after 60-100 years, Picea increasingly dominates the canopy while Tsuga heterophylla regenerate within the understory. The original cohort of Picea dies off after 300-500 years, which allows Tsuga to dominate old growth terrace forests. Picea or Alnus do not tend to regenerate under these dense canopies and without disturbance Tsuga may remain dominant indefinitely. Understory composition was related to landform age, however species distributions at low elevation floodplain sites were also driven by elevation above thalweg and flood frequency. Light availability was also a significant factor in driving community composition. It appears that understory dynamics were linked to overstory succession and geomorphic development processes, which alter environmental conditions at the understory level. That is, species distributions are driven by dynamic environmental filters, which change as a result of biogeomorphic succession. Mature forest patches tended to persist longer than young forests. The landscape composition reflects a balance between episodes of hydrogeomorphic disturbance and periods of successional development. Increased hydrogeomorphic disturbance rates due to climate change have the capacity to alter the landscape composition resulting in diminished mature forest area.

View record

Stage discharge estimation using a 1d river hydraulic model and spatially-variable roughness (2010)

Stage-discharge relations (rating curves) are integral to stream gauging, yet the existing empirical calibration methods are expensive, particularly in remote areas, and are limited to low flows. Numerical modelling can provide stage-discharge relations from a single site survey, reducing the overall cost, and can be fit to changing surface conditions. This study explores a one-dimensional model to calculate theoretical stage-discharge relations for four field sites in British Columbia that range in bed stability, bed structure, hydrology and sediment supply. However, due to the non-linear relation between flow and roughness we do not assume the conventional reach-averaged roughness and instead employ a spatially-distributed roughness model. Furthermore, based on local grain size distribution and refined field survey technique, new formulae for wetted perimeter, flow area, and flow depth were developed that eliminate commonly held modelling assumptions and reduce topographic error. The results show (1) good agreement with Water Survey of Canada measurements, (2) distributed roughness provided an improvement over spatially-averaged roughness, (3) spatial variability of the geomorphology within the channel reach leads to shifts in the stage-discharge relations and high sediment amplifies those shifts, and (4) the relations must be re-evaluated following events that mobilize the bed. The method can be used to estimate high flows and flows in remote locations and it does not require calibration.

View record

Recent Tri-Agency Grants

The following is a selection of grants for which the faculty member was principal investigator or co-investigator. Currently, the list only covers Canadian Tri-Agency grants from years 2013/14-2016/17 and excludes grants from any other agencies.

  • From chronology to frequency pairing: a paradigm shift in the science of forests and floods - Natural Sciences and Engineering Research Council of Canada (NSERC) - Discovery Grants Program - Individual (2013/2014)

Current Students & Alumni

This is a small sample of students and/or alumni that have been supervised by this researcher. It is not meant as a comprehensive list.

Membership Status

Member of G+PS

Program Affiliations


If this is your researcher profile you can log in to the Faculty & Staff portal to update your details and provide recruitment preferences.