This catalog contains annual 30 m spatial resolution snow dynamics metrics for each snow-year from 2018-2019 to 2023-2024 for all of Canada. We gather all Landsat and Sentinel-2 images collected over Canada and identify the status of each pixel observation on the image collection date: snow (and ice), non-snow (i.e., land, water), unclear (i.e., clouds, shadows). We built an algorithm to calculate snow cover metrics for each pixel during each winter: start date of the first (and biggest) snow period [startF, startB], end date of the last (and biggest) snow period [endL, endB], number of days with snow cover in total (or in the biggest snow period) [lengthT, lengthB], number of snow periods (i.e., separated times with multiple confirmed snow observations) [periods], and a status classification (e.g., continuous snow, snow free) [status]. We do not obtain a clear observation every day because of satellite orbit frequencies and clouds. This means that timing-based metrics are identified by the middle date between two clear observations, with uncertainty quantified as half the length of the gap (i.e., ± days) [startF_u, startB_u, endL_u, endB_u, lengthT_u, lengthB_u].

This dataset corresponds to daily snow cover percentage at 1km resolution grid over land areas of Canada from 2006-2010. The data are subsampled by 4km to reduce data volumes and considering the geolocation uncertainty of the input satellite imagery. The daily maps are generated by assimilation of daily cloud screened NOAA AVHRR satellite imagery and Canadian Meteorological Centre (CMC) snow depth analysis snow depth and density fields within an off-line version of the CMC daily snow depth model. The snow depth model is modified to include snowpack reflectance model and a surface radiative transfer scheme that relates vegetation and snowpack reflectance to top-of-canopy bi-directional reflectance. A logistic vegetation phenology model is used to parameterize temporal dynamics of canopy leaf area index. A per-pixel particle filter with a 30 day moving window is applied to assimilation observations corresponding to 1km resolution visible band directional reflectance and normalized difference vegetation index and 24km CMC daily snow depth and monthly snow density fields. The assimilation is forced using daily air temperature and precipitation fields. Validation of the datasets has been performed by comparison to MODIS snow cover maps and in-situ snow depth stations across Canada. Validation suggests similar accuracy to MODIS snow cover products over relatively flat terrain. Validation over mountainous regions is ongoing.

This data contains location information for 1 of Ontario’s snow monitoring networks: * Surface Water Monitoring Centre (SWMC) Snow course data is collected by: * conservation authorities * Ontario Power Generation * Ministry of Natural Resources (MNR) districts Data is collected twice a month from November 15 until May 15. The Surface Water Monitoring Centre uses this data to assess: * current snow cover * frozen ground conditions * snowpack * potential snowmelt * contributions to streamflow The snow data is located in a corporate water and climate database. This data helps MNR and conservation authorities assess the potential for flood at the local and provincial scale.

The daily climate records database, also known as Long Term Climate Extremes (LTCE), was developed to address the fragmentation of climate information due to station changes (opening, closing, relocation, etc.) over time. For approximately 750 locations in Canada, "virtual" climate stations have been developed by joining (threading) climate data for an urban location, from nearby stations to make long-term records. Each long-term record consists of the extremes (record values) of daily maximum/minimum temperatures, total precipitation and snowfall for each day of the year. Many of the longest data sets of extremes date as far back as the 1800s. This data provides the daily extremes of record for Snowfall for each day of the year. Daily elements include: Greatest Snowfall.

This mosaic is calculated over the North American domain with a horizontal spatial resolution of 1 km. This mosaic therefore includes all the Canadian and American radars available in the network and which can reach a maximum of 180 contributing radars. To better represent precipitation over the different seasons, this mosaic renders in mm/h to represent rain and in cm/h to represent snow. For the two precipitation types (rain and snow), we use two different mathematical relationships to convert the reflectivity by rainfall rates (mm/h rain cm/h for snow). This is a hybrid mosaic from DPQPE (Dual-Pol Quantitative Precipitation Estimation) for S-Band radars. For the US Nexrad radars, ECCC uses the most similar product from the US Meteorological Service (NOAA). This product displays radar reflectivity converted into precipitation rates, using the same formulas as the Canadian radars.

Fire weather refers to weather conditions that are conducive to fire. These conditions determine the fire season, which is the period(s) of the year during which fires are likely to start, spread and do sufficient damage to warrant organized fire suppression.The length of fire season is the difference between the start- and end-of-fire-season dates. These are defined by the Canadian Forest Fire Weather Index (FWI; http://cwfis.cfs.nrcan.gc.ca/) start-up and end dates. Start-up occurs when the station has been snow-free for 3 consecutive days, with noon temperatures of at least 12°C. For stations that do not report significant snow cover during the winter (i.e., less than 10 cm or snow-free for 75% of the days in January and February), start-up occurs when the mean daily temperature has been 6°C or higher for 3 consecutive days. The fire season ends with the onset of winter, generally following 7 consecutive days of snow cover. If there are no snow data, shutdown occurs following 7 consecutive days with noon temperatures lower than or equal to 5°C.Historical climate conditions were derived from the 1981–2010 Canadian Climate Normals. Future projections were computed using two different Representative Concentration Pathways (RCP). RCPs are different greenhouse gas concentration trajectories adopted by the Intergovernmental Panel on Climate Change (IPCC) for its fifth Assessment Report. RCP 2.6 (referred to as rapid emissions reductions) assumes that greenhouse gas concentrations peak between 2010-2020, with emissions declining thereafter. In the RCP 8.5 scenario (referred to as continued emissions increases) greenhouse gas concentrations continue to rise throughout the 21st century.Provided layer: the fire season length across Canada for a reference period (1981-2010).

Fire weather refers to weather conditions that are conducive to fire. These conditions determine the fire season, which is the period(s) of the year during which fires are likely to start, spread and do sufficient damage to warrant organized fire suppression.The length of fire season is the difference between the start- and end-of-fire-season dates. These are defined by the Canadian Forest Fire Weather Index (FWI; http://cwfis.cfs.nrcan.gc.ca/) start-up and end dates. Start-up occurs when the station has been snow-free for 3 consecutive days, with noon temperatures of at least 12°C. For stations that do not report significant snow cover during the winter (i.e., less than 10 cm or snow-free for 75% of the days in January and February), start-up occurs when the mean daily temperature has been 6°C or higher for 3 consecutive days. The fire season ends with the onset of winter, generally following 7 consecutive days of snow cover. If there are no snow data, shutdown occurs following 7 consecutive days with noon temperatures lower than or equal to 5°C.Historical climate conditions were derived from the 1981–2010 Canadian Climate Normals. Future projections were computed using two different Representative Concentration Pathways (RCP). RCPs are different greenhouse gas concentration trajectories adopted by the Intergovernmental Panel on Climate Change (IPCC) for its fifth Assessment Report. RCP 2.6 (referred to as rapid emissions reductions) assumes that greenhouse gas concentrations peak between 2010-2020, with emissions declining thereafter. In the RCP 8.5 scenario (referred to as continued emissions increases) greenhouse gas concentrations continue to rise throughout the 21st century.Provided layer: difference in projected fire season length for the short-term (2011-2040) under the RCP 8.5 (continued emissions increases) compared to reference period across Canada.

Fire weather refers to weather conditions that are conducive to fire. These conditions determine the fire season, which is the period(s) of the year during which fires are likely to start, spread and do sufficient damage to warrant organized fire suppression.The length of fire season is the difference between the start- and end-of-fire-season dates. These are defined by the Canadian Forest Fire Weather Index (FWI; http://cwfis.cfs.nrcan.gc.ca/) start-up and end dates. Start-up occurs when the station has been snow-free for 3 consecutive days, with noon temperatures of at least 12°C. For stations that do not report significant snow cover during the winter (i.e., less than 10 cm or snow-free for 75% of the days in January and February), start-up occurs when the mean daily temperature has been 6°C or higher for 3 consecutive days. The fire season ends with the onset of winter, generally following 7 consecutive days of snow cover. If there are no snow data, shutdown occurs following 7 consecutive days with noon temperatures lower than or equal to 5°C.Historical climate conditions were derived from the 1981–2010 Canadian Climate Normals. Future projections were computed using two different Representative Concentration Pathways (RCP). RCPs are different greenhouse gas concentration trajectories adopted by the Intergovernmental Panel on Climate Change (IPCC) for its fifth Assessment Report. RCP 2.6 (referred to as rapid emissions reductions) assumes that greenhouse gas concentrations peak between 2010-2020, with emissions declining thereafter. In the RCP 8.5 scenario (referred to as continued emissions increases) greenhouse gas concentrations continue to rise throughout the 21st century.Provided layer: difference in projected fire season length for the long-term (2071-2100) under the RCP 8.5 (continued emissions increases) compared to reference period across Canada.

Seasonal and annual multi-model ensembles of projected change (also known as anomalies) in snow depth based on an ensemble of twenty-eight Coupled Model Intercomparison Project Phase 5 (CMIP5) global climate models are available for 1900-2100. Projected change in snow depth is with respect to the reference period of 1986-2005 and expressed as a percentage (%). The 5th, 25th, 50th, 75th and 95th percentiles of the ensemble of snow depth change are available for the historical time period, 1900-2005, and for emission scenarios, RCP2.6, RCP4.5 and RCP8.5, for 2006-2100. Twenty-year average changes in snow depth (%) for four time periods (2021-2040; 2041-2060; 2061-2080; 2081-2100), with respect to the reference period of 1986-2005, for RCP2.6, RCP4.5 and RCP8.5 are also available in a range of formats. The median projected change across the ensemble of CMIP5 climate models is provided. Note: Projections among climate models can vary because of differences in their underlying representation of earth system processes. Thus, the use of a multi-model ensemble approach has been demonstrated in recent scientific literature to likely provide better projected climate change information.

Snow and ice are important hydrological resources. Their minimum spatial extent here referred to as annual minimum snow/ice (MSI) cover, plays a very important role as an indicator of long-term changes and baseline capacity for surface water storage. The MSI probability is derived from sequence of seventeen 10-day clear-sky composites corresponding to April, 1 to September, 20 warm period for each year since 2000. Data from Moderate Resolution Imaging Spectroradiometer (MODIS) on Terra satellite for the period since 2000 have been processed with the special technology developed at the Canada Centre for Remote Sensing (CCRS) as described in Trishchenko, 2016; Trishchenko et al., 2016; 2009, 2006, Trishchenko and Ungureanu, 2021, Khlopenkov and Trishchenko, 2008, Luo et al., 2008. The presence of snow or ice is determined for each pixel of the image based on snow/ice scene identification procedure and the probability if computed for the entire warm season as a ratio of number of snow/ice flags to the total number of pixels available (less or equal to 17). The minimum snow and ice extent can be derived from the probability map by applying a certain threshold. New data version V5.0 replaces previous version V4.0 for all data available since 2000. All MSI files were reprocessed for all MODIS input data based on collection 6.1. The output format has not changed since previous version. It is described in Trishchenko (2024). The impact of input data change is small and can be detected only for time interval 2000-2015. Data starting 2016 has been already derived using MODIS collection 6.1 input.The differences between the MSI data based on MODIS Collection 5 (i.e. MSI V4) versus MODIS Collection 6.1 (i.e. MSI V5), on average, are quite small. The region-wide relative difference in the MSI extent varies from -3.97% to +1.75%. The mean value is -0.14%, the median value is 0.18% and standard deviation is 1.83%. As such, we do not expect any sizeable impact of the version change on our previous conclusions regarding trends and climate variations, except for refining the relative values of statistical parameters within the range of a few percents. References:TRISHCHENKO, A.P., 2024: Probability maps of the annual minimum snow and ice (MSI) presence over April,1 to September, 20 period since 2000 derived from MODIS 250m imagery over Canada and neighbouring regions. Data format description. CCRS, NRCan. 4pp.