Votre recherche

Summary Projected climate change effects on streamflow are investigated for the 2041–2070 horizon following the SRES A2 emissions scenario over two snowmelt-dominated catchments in Canada. A 16-member ensemble of SWAT hydrological model (HM) simulations, based on a comprehensive ensemble of the Canadian Regional Climate Model (CRCM) simulations driven by two global climate models (GCMs), with five realizations of the Canadian CGCM3 and three realizations of the German ECHAM5 is established per catchment. This study aims to evaluate, once model bias has been removed by statistical post-processing (SP), how the RCM-simulated climate changes differ from those of the parent GCMs, and how they affect the assessment of climate change-induced hydrological impacts at the catchment scale. The variability of streamflow caused by the use of different SP methods (mean-based versus distribution-based) within each statistical post-processing pathway of climate model outputs (bias correction versus perturbation) is also evaluated, as well as the uncertainty of natural climate variability. The simulations cover 1971–2000 in the reference period and 2041–2070 in the future period. For a set of criteria, results based on raw and statistically post-processed model outputs for the reference climate are compared with observations. This process demonstrates that SP is important not only for GCMs outputs, but also for CRCM outputs. SP leads to a high level of agreement between the CRCM and the driving GCMs in reproducing patterns of observed climate. The ensemble spread of the climate change signal on streamflow is large and varies with catchments and hydrological periods (winter/summer flows). The results of various hydrological indicators show that most of the uncertainty arises from the natural climate variability followed by the statistical post-processing. The uncertainty linked to the choice of statistical pathway is much larger than that associated with the choice of the method in quantifying the hydrological impacts. We find that the incorporation of dynamical downscaling of global models through the CRCM as an intermediate step in the GCM–RCM–SP–HM model chain does not lead to considerable differences in the assessment of the climate change impacts on streamflow for the study catchments.

Canada has experienced some of the most rapid warming on Earth over the past few decades with a warming rate about twice that of the global mean temperature since 1948. Long-term warming is observed in Canada’s annual, winter and summer mean temperatures, and in the annual coldest and hottest daytime and nighttime temperatures. The causes of these changes are assessed by comparing observed changes with climate model simulated responses to anthropogenic and natural (solar and volcanic) external forcings. Most of the observed warming of 1.7°C increase in annual mean temperature during 1948–2012 [90% confidence interval (1.1°, 2.2°C)] can only be explained by external forcing on the climate system, with anthropogenic influence being the dominant factor. It is estimated that anthropogenic forcing has contributed 1.0°C (0.6°, 1.5°C) and natural external forcing has contributed 0.2°C (0.1°, 0.3°C) to the observed warming. Up to 0.5°C of the observed warming trend may be associated with low frequency variability of the climate such as that represented by the Pacific decadal oscillation (PDO) and North Atlantic oscillation (NAO). Overall, the influence of both anthropogenic and natural external forcing is clearly evident in Canada-wide mean and extreme temperatures, and can also be detected regionally over much of the country.

The impact of snow-atmosphere coupling on climate variability and extremes over North America is investigated using modeling experiments with the fifth generation Canadian Regional Climate Model (CRCM5). To this end, two CRCM5 simulations driven by ERA-Interim reanalysis for the 1981–2010 period are performed, where snow cover and depth are prescribed (uncoupled) in one simulation while they evolve interactively (coupled) during model integration in the second one. Results indicate systematic influence of snow cover and snow depth variability on the inter-annual variability of soil and air temperatures during winter and spring seasons. Inter-annual variability of air temperature is larger in the coupled simulation, with snow cover and depth variability accounting for 40–60% of winter temperature variability over the Mid-west, Northern Great Plains and over the Canadian Prairies. The contribution of snow variability reaches even more than 70% during spring and the regions of high snow-temperature coupling extend north of the boreal forests. The dominant process contributing to the snow-atmosphere coupling is the albedo effect in winter, while the hydrological effect controls the coupling in spring. Snow cover/depth variability at different locations is also found to affect extremes. For instance, variability of cold-spell characteristics is sensitive to snow cover/depth variation over the Mid-west and Northern Great Plains, whereas, warm-spell variability is sensitive to snow variation primarily in regions with climatologically extensive snow cover such as northeast Canada and the Rockies. Furthermore, snow-atmosphere interactions appear to have contributed to enhancing the number of cold spell days during the 2002 spring, which is the coldest recorded during the study period, by over 50%, over western North America. Additional results also provide useful information on the importance of the interactions of snow with large-scale mode of variability in modulating temperature extreme characteristics.

This study focuses on the evaluation of daily precipitation and temperature climate indices and extremes simulated by an ensemble of 12 Regional Climate Model (RCM) simulations from the ARCTIC-CORDEX experiment with surface observations in the Canadian Arctic from the Adjusted Historical Canadian Climate Dataset. Five global reanalyses products (ERA-Interim, JRA55, MERRA, CFSR and GMFD) are also included in the evaluation to assess their potential for RCM evaluation in data sparse regions. The study evaluated the means and annual anomaly distributions of indices over the 1980–2004 dataset overlap period. The results showed that RCM and reanalysis performance varied with the climate variables being evaluated. Most RCMs and reanalyses were able to simulate well climate indices related to mean air temperature and hot extremes over most of the Canadian Arctic, with the exception of the Yukon region where models displayed the largest biases related to topographic effects. Overall performance was generally poor for indices related to cold extremes. Likewise, only a few RCM simulations and reanalyses were able to provide realistic simulations of precipitation extreme indicators. The multi-reanalysis ensemble provided superior results to individual datasets for climate indicators related to mean air temperature and hot extremes, but not for other indicators. These results support the use of reanalyses as reference datasets for the evaluation of RCM mean air temperature and hot extremes over northern Canada, but not for cold extremes and precipitation indices.

AbstractIn this time of a changing climate, it is important to know whether lake levels will rise, potentially causing flooding, or river flows will dry up during abnormally dry weather. The Great Lakes region is the largest freshwater lake system in the world. Moreover, agriculture, industry, commerce, and shipping are active in this densely populated region. Environment and Climate Change Canada (ECCC) recently implemented the Water Cycle Prediction System (WCPS) over the Great Lakes and St. Lawrence River watershed (WCPS-GLS version 1.0) following a decade of research and development. WCPS, a network of linked models, simulates the complete water cycle, following water as it moves from the atmosphere to the surface, through the river network and into lakes, and back to the atmosphere. Information concerning the water cycle is passed between the models. WCPS is the first short-to-medium-range prediction system of the complete water cycle to be run on an operational basis anywhere. It currently produces ...

AbstractSoil moisture–atmosphere interactions play a key role in modulating climate variability and extremes. This study investigates how soil moisture–atmosphere coupling may affect future extreme events, particularly the role of projected soil moisture in modulating the frequency and maximum duration of hot spells over North America, using the fifth-generation Canadian Regional Climate Model (CRCM5). With this objective, CRCM5 simulations, driven by two coupled general circulation models (MPI-ESM and CanESM2), are performed with and without soil moisture–atmosphere interactions for current (1981–2010) and future (2071–2100) climates over North America, for representative concentration pathways (RCPs) 4.5 and 8.5. Analysis indicates that, in future climate, the soil moisture–temperature coupling regions, located over the Great Plains in the current climate, will expand farther north, including large parts of central Canada. Results also indicate that soil moisture–atmosphere interactions will play an imp...

Abstract. River ice is a common occurrence in cold climate hydrological systems. The annual cycle of river ice formation, growth, decay and clearance can include low flows and ice jams, as well as mid-winter and spring break-up events. Reports and associated data on river ice occurrence are often limited to site and season-specific studies. Within Canada, the National Hydrometric Program (NHP) operates a network of gauging stations with water level as the primary measured variable to derive discharge. In the late 1990s, the Water Science and Technology Directorate of Environment and Climate Change Canada initiated a long-term effort to compile, archive and extract river ice related information from NHP hydrometric records. This data article describes the original research data set produced by this near 20-year effort: the Canadian River Ice Database (CRID). The CRID holds almost 73,000 variables from a network of 196 NHP stations throughout Canada that were in operation within the period 1894 to 2015. Over 100,000 paper and digital files were reviewed representing 10,378 station-years of active operation. The task of compiling this database involved manual extraction and input of more than 460,000 data entries on water level, discharge, date, time and data quality rating. Guidelines on the data extraction, rating procedure and challenges are provided. At each location, a time series of up to 15 variables specific to the occurrence of freeze-up and winter-low events, mid-winter break-up, ice thickness, spring break-up and maximum open-water level were compiled. This database follows up on several earlier efforts to compile information on river ice, which are summarized herein, and expands the scope and detail for use in Canadian river ice research and applications. Following the Government of Canada Open Data initiative, this original river ice data set is available at: https://doi.org/10.18164/c21e1852-ba8e-44af-bc13-48eeedfcf2f4 (de Rham et al., 2020).

Abstract Flood quantile estimation at sites with little or no data is important for the adequate planning and management of water resources. Regional Hydrological Frequency Analysis (RFA) deals with the estimation of hydrological variables at ungauged sites. Random Forest (RF) is an ensemble learning technique which uses multiple Classification and Regression Trees (CART) for classification, regression, and other tasks. The RF technique is gaining popularity in a number of fields because of its powerful non-linear and non-parametric nature. In the present study, we investigate the use of Random Forest Regression (RFR) in the estimation step of RFA based on a case study represented by data collected from 151 hydrometric stations from the province of Quebec, Canada. RFR is applied to the whole data set and to homogeneous regions of stations delineated by canonical correlation analysis (CCA). Using the Out-of-bag error rate feature of RF, the optimal number of trees for the dataset is calculated. The results of the application of the CCA based RFR model (CCA-RFR) are compared to results obtained with a number of other linear and non-linear RFA models. CCA-RFR leads to the best performance in terms of root mean squared error. The use of CCA to delineate neighborhoods improves considerably the performance of RFR. RFR is found to be simple to apply and more efficient than more complex models such as Artificial Neural Network-based models.

Study region Hudson Bay Lowlands watersheds, Ontario, Canada. Study Focus The rivers in the Hudson Bay Lowlands are a major source of freshwater entering the Arctic Ocean and they also cause major floods. In recent decades, this region has been affected by major changes in hydroclimatic processes attributed to climate change and natural climate variability. In this study, we used ERA5 reanalysis data, hydrometric observations, and the hydrological model MESH, to investigate the impact of atmospheric circulation on the inter-decadal variability of streamflow between 1979 and 2018 in the Hudson Bay Lowlands. The natural climate variability was assessed using a weather regimes approach based on the discretization of daily geopotential height anomalies (Z500) from ERA5 reanalysis, as well as large scale oceanic and atmospheric variability modes. New hydrological insights The results showed an anomalous convergence of atmospheric moisture flux between 1995–2008 that enhanced precipitation and increased streamflow in the western part of the region. This moisture convergence was likely driven by the combination of (i) low pressure anomalies in the East Coast of North America and (ii) low pressure anomalies in western regions of Canada, associated with the cold phase of the pacific decadal oscillation (PDO). Since 2009, streamflow remains high, likely due to more groundwater discharge associated with the degradation of permafrost.

Fluvial systems in southern Ontario are regularly affected by widespread early-spring flood events primarily caused by rain-on-snow events. Recent studies have shown an increase in winter floods in this region due to increasing winter temperature and precipitation. Streamflow simulations are associated with uncertainties mainly due to the different scenarios of greenhouse gas emissions, global climate models (GCMs) or the choice of the hydrological model. The internal variability of climate, defined as the chaotic variability of atmospheric circulation due to natural internal processes within the climate system, is also a source of uncertainties to consider. Uncertainties of internal variability can be assessed using hydrological models fed by downscaled data of a global climate model large ensemble (GCM-LE), but GCM outputs have too coarse of a scale to be used in hydrological modeling. The Canadian Regional Climate Model Large Ensemble (CRCM5-LE), a 50-member ensemble downscaled from the Canadian Earth System Model version 2 Large Ensemble (CanESM2-LE), was developed to simulate local climate variability over northeastern North America under different future climate scenarios. In this study, CRCM5-LE temperature and precipitation projections under an RCP8.5 scenario were used as input in the Precipitation Runoff Modeling System (PRMS) to simulate streamflow at a near-future horizon (2026–2055) for four watersheds in southern Ontario. To investigate the role of the internal variability of climate in the modulation of streamflow, the 50 members were first grouped in classes of similar projected change in January–February streamflow and temperature and precipitation between 1961–1990 and 2026–2055. Then, the regional change in geopotential height (Z500) from CanESM2-LE was calculated for each class. Model simulations showed an average January–February increase in streamflow of 18 % (±8.7) in Big Creek, 30.5 % (±10.8) in Grand River, 29.8 % (±10.4) in Thames River and 31.2 % (±13.3) in Credit River. A total of 14 % of all ensemble members projected positive Z500 anomalies in North America's eastern coast enhancing rain, snowmelt and streamflow volume in January–February. For these members the increase of streamflow is expected to be as high as 31.6 % (±8.1) in Big Creek, 48.3 % (±11.1) in Grand River, 47 % (±9.6) in Thames River and 53.7 % (±15) in Credit River. Conversely, 14 % of the ensemble projected negative Z500 anomalies in North America's eastern coast and were associated with a much lower increase in streamflow: 8.3 % (±7.8) in Big Creek, 18.8 % (±5.8) in Grand River, 17.8 % (±6.4) in Thames River and 18.6 % (±6.5) in Credit River. These results provide important information to researchers, managers, policymakers and society about the expected ranges of increase in winter streamflow in a highly populated region of Canada, and they will help to explain how the internal variability of climate is expected to modulate the future streamflow in this region.

AbstractA snow model forced by temperature and precipitation is used to simulate the spatial distribution of snow water equivalent (SWE) over a 600,000 km2 portion of the province of Quebec, Canada. We propose to improve model simulations by assimilating SWE data from sporadic manual snow surveys with a particle filter. A temporally and spatially correlated perturbation of the meteorological forcing is used to generate the set of particles. The magnitude of the perturbations is fixed objectively. First, the particle filter and direct insertion were both applied on 88 sites for which measured SWE consist of more or less five values per year over a period of 17 years. The temporal correlation of perturbations enables to improve the accuracy and the ensemble dispersion of the particle filter, while the spatial correlation lead to a spatial coherence in the particle weights. The spatial estimates of SWE obtained with the particle filter are compared with those obtained through optimal interpolation of the sno...

AbstractThe Canadian Land Data Assimilation System (CaLDAS) has been developed at the Meteorological Research Division of Environment Canada (EC) to better represent the land surface initial states in environmental prediction and assimilation systems. CaLDAS is built around an external land surface modeling system and uses the ensemble Kalman filter (EnKF) methodology. A unique feature of CaLDAS is the use of improved precipitation forcing through the assimilation of precipitation observations. An ensemble of precipitation analyses is generated by combining numerical weather prediction (NWP) model precipitation forecasts with precipitation observations. Spatial phasing errors to the NWP first-guess precipitation forecasts are more effective than perturbations to the precipitation observations in decreasing (increasing) the exceedance ratio (uncertainty ratio) scores and generating flatter, more reliable ranked histograms. CaLDAS has been configured to assimilate L-band microwave brightness temperature TB ...

Over the last decades, different methods have been used by hydrologists to extend observed hydro-climatic time series, based on other data sources, such as tree rings or sedimentological datasets. For example, tree ring multi-proxies have been studied for the Caniapiscau Reservoir in northern Quebec (Canada), leading to the reconstruction of flow time series for the last 150 years. In this paper, we applied a new hydro-climatic reconstruction method on the Caniapiscau Reservoir and compare the obtained streamflow time series against time series derived from dendrohydrology by other authors on the same catchment and study the natural streamflow variability over the 1881–2011 period in that region. This new reconstruction is based not on natural proxies but on a historical reanalysis of global geopotential height fields, and aims firstly to produce daily climatic time series, which are then used as inputs to a rainfall–runoff model in order to obtain daily streamflow time series. The performances of the hydro-climatic reconstruction were quantified over the observed period, and showed good performances, in terms of both monthly regimes and interannual variability. The streamflow reconstructions were then compared to two different reconstructions performed on the same catchment by using tree ring data series, one being focused on mean annual flows and the other on spring floods. In terms of mean annual flows, the interannual variability in the reconstructed flows was similar (except for the 1930–1940 decade), with noteworthy changes seen in wetter and drier years. For spring floods, the reconstructed interannual variabilities were quite similar for the 1955–2011 period, but strongly different between 1880 and 1940. The results emphasize the need to apply different reconstruction methods on the same catchments. Indeed, comparisons such as those above highlight potential differences between available reconstructions and, finally, allow a retrospective analysis of the proposed reconstructions of past hydro-climatological variabilities.

Abstract River confluences are characterized by a complex mixing zone with three‐dimensional (3D) turbulent structures which have been described as both streamwise‐oriented structures and Kelvin–Helmholtz (KH) vertical‐oriented structures. The latter are visible where there is a turbidity difference between the two tributaries, whereas the former are usually derived from mean velocity measurements or numerical simulations. Few field studies recorded turbulent velocity fluctuations at high frequency to investigate these structures, particularly at medium‐sized confluences where logistical constraints make it difficult to use devices such as acoustic doppler velocimeter (ADV). This study uses the ice cover present at the confluence of the Mitis and Neigette Rivers in Quebec (Canada) to obtain long‐duration, fixed measurements along the mixing zone. The confluence is also characterized by a marked turbidity difference which allows to investigate the mixing zone dynamics from drone imagery during ice‐free conditions. The aim of the study is to characterize and compare the flow structure in the mixing zone at a medium‐sized (~40 m) river confluence with and without an ice cover. Detailed 3D turbulent velocity measurements were taken under the ice along the mixing plane with an ADV through eight holes at around 20 positions on the vertical. For ice‐free conditions, drone imagery results indicate that large (KH) coherent structures are present, occupying up to 50% of the width of the parent channel. During winter, the ice cover affects velocity profiles by moving the highest velocities towards the centre of the profiles. Large turbulent structures are visible in both the streamwise and lateral velocity components. The strong correlation between these velocity components indicates that KH vortices are the dominating coherent structures in the mixing zone. A spatio‐temporal conceptual model is presented to illustrate the main differences on the 3D flow structure at the river confluence with and without the ice cover. © 2019 John Wiley & Sons, Ltd.

AbstractA new land surface scheme has been developed at Environment and Climate Change Canada (ECCC) to provide surface fluxes of momentum, heat, and moisture for the Global Environmental Multiscale (GEM) atmospheric model. In this study, the performance of the Soil, Vegetation, and Snow (SVS) scheme in estimating the surface and root-zone soil moisture is evaluated against the Interactions between Soil, Biosphere, and Atmosphere (ISBA) scheme currently used operationally at ECCC within GEM for numerical weather prediction. In addition, the sensitivity of SVS soil moisture results to soil texture and vegetation data sources (type and fractional coverage) has been explored. The performance of SVS and ISBA was assessed against a large set of in situ observations as well as the brightness temperature data from the Soil Moisture Ocean Salinity (SMOS) satellite over North America. The results indicate that SVS estimates the time evolution of soil moisture more accurately, and compared to ISBA, results in highe...

In time series of essential climatological variables, many discontinuities are created not by climate factors but changes in the measuring system, including relocations, changes in instrumentation, exposure or even observation practices. Some of these changes occur due to reorganization, cost-efficiency or innovation. In the last few decades, station movements have often been accompanied by the introduction of an automatic weather station (AWS). Our study identifies the biases in daily maximum and minimum temperatures using parallel records of manual and automated observations. They are selected to minimize the differences in surrounding environment, exposition, distance and difference in elevation. Therefore, the type of instrumentation is the most important biasing factor between both measurements. The pairs of weather stations are located in Piedmont, a region of Italy, and in Gaspe Peninsula, a region of Canada. They have 6years of overlapping period on average, and 5110 daily values. The approach implemented for the comparison is divided in four main parts: a statistical characterization of the daily temperature series; a comparison between the daily series; a comparison between the types of events, heat wave, cold wave and normal events; and a verification of the homogeneity of the difference series. Our results show a higher frequency of warm (+10%) and extremely warm (+35%) days in the automated system, compared with the parallel manual record. Consequently, the use of a composite record could significantly bias the calculation of extreme events.

Improving Disaster Preparedness Through Mutual Catastrophe Insurance In “A Mutual Catastrophe Insurance Framework for Horizontal Collaboration in Prepositioning Strategic Reserves,” H. Zbib, B. Balcik, M.-È. Rancourt, and G. Laporte present an innovative approach to collaborative disaster preparedness. The novel framework considers a risk-averse mutual insurer offering multiyear insurance contracts with coverage deductibles and limits to a portfolio of risk-averse policyholders. It is designed to foster horizontal collaboration among policyholders for joint disaster preparedness by effectively integrating operational and financial functions. The problem is modeled as a large-scale nonlinear multistage stochastic program and solved by using an effective Benders decomposition algorithm. The framework is validated with real data from 18 Caribbean countries focusing on hurricane preparedness. Given the predicted impacts of climate change, the proposed multiyear mutual catastrophe insurance framework promises to reshape global disaster preparedness and make a profound societal impact by providing a transparent disaster financing plan to protect vulnerable regions. The study’s findings stress the importance of long-term cooperation, prenegotiation of indemnification policies, and strategic setting of deductibles and limits by taking into account the correlation between policyholders. , We develop a mutual catastrophe insurance framework for the prepositioning of strategic reserves to foster horizontal collaboration in preparedness against low-probability high-impact natural disasters. The framework consists of a risk-averse insurer pooling the risks of a portfolio of risk-averse policyholders. It encompasses the operational functions of planning the prepositioning network in preparedness for incoming insurance claims, in the form of units of strategic reserves, setting coverage deductibles and limits of policyholders, and providing insurance coverage to the claims in the emergency response phase. It also encompasses the financial functions of ensuring the insurer’s solvency by efficiently managing its capital and allocating yearly premiums among policyholders. We model the framework as a very large-scale nonlinear multistage stochastic program, and solve it through a Benders decomposition algorithm. We study the case of Caribbean countries establishing a horizontal collaboration for hurricane preparedness. Our results show that the collaboration is more effective when established over a longer planning horizon, and is more beneficial when outsourcing becomes expensive. Moreover, the correlation of policyholders affected simultaneously under the extreme realizations and the position of their claims in their global claims distribution directly affects which policyholders get deductibles and limits. This underlines the importance of prenegotiating policyholders’ indemnification policies at the onset of collaboration. Funding: G. Laporte and M.-È. Rancourt were funded by the Canadian Natural Sciences and Engineering Research Council (NSERC) [Grants 2015-06189 and 2022-04846]. Funding was also provided by the Institute for Data Valorisation (IVADO) and the Canada Research Chair in Humanitarian Supply Chain Analytics. B. Balcik was partially supported by a grant from the Scientific and Technological Research Council of Turkey (TUBITAK) 2219 program. This support is gratefully acknowledged. Supplemental Material: The online appendix is available at https://doi.org/10.1287/opre.2021.0141 .

Abstract. Climate change affects natural streamflow regimes globally. To assess alterations in streamflow regimes, typically temporal variations in one or a few streamflow characteristics are taken into account. This approach, however, cannot see simultaneous changes in multiple streamflow characteristics, does not utilize all the available information contained in a streamflow hydrograph, and cannot describe how and to what extent streamflow regimes evolve from one to another. To address these gaps, we conceptualize streamflow regimes as intersecting spectrums that are formed by multiple streamflow characteristics. Accordingly, the changes in a streamflow regime should be diagnosed through gradual, yet continuous changes in an ensemble of streamflow characteristics. To incorporate these key considerations, we propose a generic algorithm to first classify streams into a finite set of intersecting fuzzy clusters. Accordingly, by analyzing how the degrees of membership to each cluster change in a given stream, we quantify shifts from one regime to another. We apply this approach to the data, obtained from 105 natural Canadian streams, during the period of 1966 to 2010. We show that natural streamflow in Canada can be categorized into six regime types, with clear hydrological and geographical distinctions. Analyses of trends in membership values show that alterations in natural streamflow regimes vary among different regions. Having said that, we show that in more than 80 % of considered streams, there is a dominant regime shift that can be attributed to simultaneous changes in streamflow characteristics, some of which have remained previously unknown. Our study not only introduces a new globally relevant algorithm for identifying changing streamflow regimes but also provides a fresh look at streamflow alterations in Canada, highlighting complex and multifaceted impacts of climate change on streamflow regimes in cold regions.

Abstract. In northern cold-temperate countries, a large portion of annual streamflow is produced by spring snowmelt, which often triggers floods. It is important to have spatial information about snow parameters such as snow water equivalent (SWE), which can be incorporated into hydrological models, making them more efficient tools for improved decision-making. The future Terrestrial Snow Mass Mission (TSMM) aims to provide high-resolution spatially distributed SWE information; thus, spatial SWE calibration should be considered along with conventional streamflow calibration for model optimization since the overall water balance is often a key objective in the hydrological modelling. The present research implements a unique spatial pattern metric in a multi-objective framework for calibration approach of hydrological models and attempts to determine whether raw SNODAS data can be utilized for hydrological model calibration. The SPAtial Efficiency (SPAEF) metric is explored for spatially calibrating SWE. The HYDROTEL hydrological model is applied to the Au Saumon River Watershed (∽1120 km2) in Eastern Canada using MSWEP precipitation data and ERA-5 land reanalysis temperature data as input to generate high-resolution SWE and streamflow. Different calibration experiments are performed combining Nash-Sutcliffe efficiency (NSE) for streamflow and root-mean-square error (RMSE), and SPAEF for SWE, using the Dynamically Dimensioned Search (DDS) and Pareto Archived Multi-Objective Optimization (PADDS) algorithms. Results of the study demonstrate that multi-objective calibration outperforms sequential calibration in terms of model performance. Traditional model calibration involving only streamflow produced slightly higher NSE values; however, the spatial distribution of SWE could not be adequately maintained. This study indicates that utilizing SPAEF for spatial calibration of snow parameters improved streamflow prediction compared to the conventional practice of using RMSE for calibration. SPAEF is further implied to be a more effective metric than RMSE for both sequential and multi-objective calibration. During validation, the calibration experiment incorporating multi-objective SPAEF exhibits enhanced performance in terms of NSE and Kling-Gupta Efficiency (KGE) compared to calibration experiment solely based on NSE. This observation supports the notion that incorporating SPAEF computed on raw SNODAS data within the calibration framework results in a more robust hydrological model.

Abstract Landslides involving sensitive clays are recurrent events in the world's northern regions and are especially notorious in eastern Canada. The two critical factors that separate sensitive clay landslides from traditional slope stability analysis are the highly brittle behavior in undrained conditions (strain-softening) characteristic of progressive or retrogressive failures and the large deformations associated with them. Conventional limit equilibrium analysis has numerous shortcomings in incorporating these characteristics when assessing landslides in sensitive clays. This paper presents an extensive literature review of the failure mechanics characteristics of landslides in sensitive clays and the existing constitutive models and numerical tools to analyze such slopes' stability and post-failure behavior. The advantages and shortcomings of the different techniques to incorporate strain-softening and large deformation in the numerical modeling of sensitive clay landslides are assessed. The literature review depicts that elastoviscoplastic soil models with non-linear strain-softening laws and rate effects represent the material behavior of sensitive clays. Though several numerical models have been proposed to analyze post-failure runouts, the amount of work performed in line with sensitive clay landslides is very scarce. That creates an urgent need to apply and further develop advanced numerical tools for better understanding and predicting these catastrophic events.