Environment and Climate Change Canada has initiated the production of a 1980–2018, 10 km, North American precipitation and surface reanalysis. ERA-Interim is used to initialize the Global Deterministic Reforecast System (GDRS) at a 39 km resolution. Its output is then dynamically downscaled to 10 km by the Regional Deterministic Reforecast System (RDRS). Coupled with the RDRS, the Canadian Land Data Assimilation System (CaLDAS) and Precipitation Analysis (CaPA) are used to produce surface and precipitation analyses. All systems used are close to operational model versions and configurations. In this study, a 7-year sample of the reanalysis (2011–2017) is evaluated. Verification results show that the skill of the RDRS is stable over time and equivalent to that of the current operational system. The impact of the coupling between RDRS and CaLDAS is explored using an early version of the reanalysis system which was run at 15 km resolution for the period 2010–2014, with and without the use of CaLDAS. Significant improvements are observed with CaLDAS in the lower troposphere and surface layer, especially for the 850 hPa dew point and absolute temperatures in summer. Precipitation is further improved through an offline precipitation analysis which allows the assimilation of additional observations of 24 h precipitation totals. The final dataset should be of particular interest for hydrological applications focusing on transboundary and northern watersheds, where existing products often show discontinuities at the border and assimilate very few – if any – precipitation observations.

From 19 to 22 June 2013, intense rainfall and concurrent snowmelt led to devastating floods in the Canadian Rockies, foothills and downstream areas of southern Alberta and southeastern British Columbia, Canada. Such an event is typical of late-spring floods in cold-region mountain headwater, combining intense precipitation with rapid melting of late-lying snowpack, and represents a challenge for hydrological forecasting systems. This study investigated the factors governing the ability to predict such an event. Three sources of uncertainty, other than the hydrological model processes and parameters, were considered: (i) the resolution of the atmospheric forcings, (ii) the snow and soil moisture initial conditions (ICs) and (iii) the representation of the soil texture. The Global Environmental Multiscale hydrological modeling platform (GEM-Hydro), running at a 1 km grid spacing, was used to simulate hydrometeorological conditions in the main headwater basins of southern Alberta during this event. The GEM atmospheric model and the Canadian Precipitation Analysis (CaPA) system were combined to generate atmospheric forcing at 10, 2.5 and 1 km over southern Alberta. Gridded estimates of snow water equivalent (SWE) from the Snow Data Assimilation System (SNODAS) were used to replace the model SWE at peak snow accumulation and generate alternative snow and soil moisture ICs before the event. Two global soil texture datasets were also used. Overall 12 simulations of the flooding event were carried out. Results show that the resolution of the atmospheric forcing affected primarily the flood volume and peak flow in all river basins due to a more accurate estimation of intensity and total amount of precipitation during the flooding event provided by CaPA analysis at convection-permitting scales (2.5 and 1 km). Basin-averaged snowmelt also changed with the resolution due to changes in near-surface wind and resulting turbulent fluxes contributing to snowmelt. Snow ICs were the main sources of uncertainty for half of the headwater basins. Finally, the soil texture had less impact and only affected peak flow magnitude and timing for some stations. These results highlight the need to combine atmospheric forcing at convection-permitting scales with high-quality snow ICs to provide accurate streamflow predictions during late-spring floods in cold-region mountain river basins. The predictive improvement by inclusion of high-elevation weather stations in the precipitation analysis and the need for accurate mountain snow information suggest the necessity of integrated observation and prediction systems for forecasting extreme events in mountain river basins.

This work explores the potential of the distributed GEM-Hydro runoff modeling platform, developed at Environment and Climate Change Canada (ECCC) over the last decade. More precisely, the aim is to develop a robust implementation methodology to perform reliable streamflow simulations with a distributed model over large and partly ungauged basins, in an efficient manner. The latest version of GEM-Hydro combines the SVS (Soil, Vegetation and Snow) land-surface scheme and the WATROUTE routing scheme. SVS has never been evaluated from a hydrological point of view, which is done here for all major rivers flowing into Lake Ontario. Two established hydrological models are confronted to GEM-Hydro, namely MESH and WATFLOOD, which share the same routing scheme (WATROUTE) but rely on different land-surface schemes. All models are calibrated using the same meteorological forcings, objective function, calibration algorithm, and basin delineation. GEM-Hydro is shown to be competitive with MESH and WATFLOOD: the NSE √ (Nash–Sutcliffe criterion computed on the square root of the flows) is for example equal to 0.83 for MESH and GEM-Hydro in validation on the Moira River basin, and to 0.68 for WATFLOOD. A computationally efficient strategy is proposed to calibrate SVS: a simple unit hydrograph is used for routing instead of WATROUTE. Global and local calibration strategies are compared in order to estimate runoff for ungauged portions of the Lake Ontario basin. Overall, streamflow predictions obtained using a global calibration strategy, in which a single parameter set is identified for the whole basin of Lake Ontario, show accuracy comparable to the predictions based on local calibration: the average NSE √ in validation and over seven subbasins is 0.73 and 0.61, respectively for local and global calibrations. Hence, global calibration provides spatially consistent parameter values, robust performance at gauged locations, and reduces the complexity and computation burden of the calibration procedure. This work contributes to the Great Lakes Runoff Inter-comparison Project for Lake Ontario (GRIP-O), which aims at improving Lake Ontario basin runoff simulations by comparing different models using the same input forcings. The main outcome of this study consists in a new generalizable methodology for implementing a distributed hydrologic model with a high computation cost in an efficient and reliable manner, over a large area with ungauged portions, using global calibration and a unit hydrograph to replace the routing component.

As part of the National Hydrological Services Transformation Initiative, Environment and Climate Change Canada (ECCC) designed and implemented the National Surface and River Prediction System (NSRPS) in order to provide surface and river flow analysis and forecast products across Canada. Within NSRPS, the Canadian Land Data Assimilation System (CaLDAS) produces snow analyses that are used to initialise the land surface model, which in turn is used to force the river routing component. Originally, CaLDAS was designed to improve atmospheric forecasts with less focus on hydrological processes. When snow data assimilation occurs, the related increments remove/add water from/to the system, which can sometimes be problematic for streamflow forecasting, in particular during the snowmelt period. In this study, a new snow analysis method introduces multiple innovations that respond to the need for higher quality snow analyses for hydrological purposes, including the use of IMS snow cover extent data instead of in situ snow depth observations. The results show that the new snow assimilation methodology brings an overall improvement to snow analyses and substantially enhances water conservation, which is reflected in the generally improved streamflow simulations. This work represents a first step towards a new snow data assimilation process in CaLDAS, with the final objective of producing a reliable snow analysis to initialise and improve NWP as well as environmental predictions, including flood and drought forecasts.

Hydrologic model intercomparison studies help to evaluate the agility of models to simulate variables such as streamflow, evaporation, and soil moisture. This study is the third in a sequence of the Great Lakes Runoff Intercomparison Projects. The densely populated Lake Erie watershed studied here is an important international lake that has experienced recent flooding and shoreline erosion alongside excessive nutrient loads that have contributed to lake eutrophication. Understanding the sources and pathways of flows is critical to solve the complex issues facing this watershed. Seventeen hydrologic and land-surface models of different complexity are set up over this domain using the same meteorological forcings, and their simulated streamflows at 46 calibration and seven independent validation stations are compared. Results show that: (1) the good performance of Machine Learning models during calibration decreases significantly in validation due to the limited amount of training data; (2) models calibrated at individual stations perform equally well in validation; and (3) most distributed models calibrated over the entire domain have problems in simulating urban areas but outperform the other models in validation.

Lakes and reservoirs have critical impacts on hydrological, biogeochemical, and ecological processes, and they should be an essential component of regional-scale hydrological and eco-hydrological models. This is particularly important in Canada with its tens of thousands of lakes. Past large-scale hydrologic modelling efforts tend to either ignore the impacts of all lakes or explicitly simulate the behaviour of only the largest lakes in a watershed. This research derives a suite of Pan-Canadian subwatershed-based lake and river routing GIS products at multiple spatial resolutions (average catchment size ranges from 60 to 306 km² while the number of lakes explicitly represented ranges from 11,000 to 87,000). These publicly available data products supply all the necessary hydrologic routing model inputs, including network topology, subwatershed geometry, channel characteristics (slope, length, roughness, and geometry), and lake characteristics (area, volume, and outlet description), and were derived primarily from the HydroLAKES, HydroBASINS and HydroSHEDS databases. All Water Survey of Canada streamflow gauging stations are used to define subwatershed outlets in the products. The routing product is used to inform a hydrologic routing model in the Raven hydrologic modelling framework. This is the first demonstration of Raven in routing-only mode. As a case study, the Hudson Bay drainage basin (∼40% of Canada) is simulated using GEM-Surf land surface model gridded runoff and recharge as inputs and includes more than 20,000 river reaches and more than 10,000 lakes explicitly represented at an hourly timestep. Uncalibrated streamflows compare reasonably well to measured streamflows at select locations. The sensitivity of the routing model prediction quality to the discretization level for represented lakes is evaluated and shows, for example, that ignoring multiple smaller lakes can have a significant impact on predictions.