Canada's National Forest Inventory (NFI) sampling program is designed to support reporting on forests at the national scale. On the other hand, continuous maps of forest attributes are required to support strategic analyses of regional policy and management issues. We have therefore produced maps covering 4.03 × 106 km2 of inventoried forest area for the 2001 base year using standardised observations from the NFI photo plots (PP) as reference data. We used the k nearest neighbours (kNN) method with 26 geospatial data layers including MODIS spectral data and climatic and topographic variables to produce maps of 127 forest attributes at a 250 × 250 m resolution. The stand-level attributes include land cover, structure, and tree species relative abundance. In this article, we report only on total live aboveground tree biomass, with all other attributes covered in the supplementary data (http://nrcresearchpress.com/doi/suppl/10.1139/cjfr-2013-0401). In general, deviations in predicted pixel-level values from those in a PP validation set are greater in mountainous regions and in areas with either low biomass or sparse PP sampling. Predicted pixel-level values are overestimated at small observed values and underestimated at large ones. Accuracy measures are improved through the spatial aggregation of pixels to 1 km2 and beyond. Overall, these new products provide unique baseline information for strategic-level analyses of forests (https://nfi.nfis.org) Collection: - **[Canada's National Forest Inventory (NFI) 2006](https://open.canada.ca/data/en/dataset/e2fadaeb-3106-4111-9d1c-f9791d83fbf4)**

14 Class - Canadian Ecological Domain Classification from Satellite Data. Satellite derived data including 1) topography, 2) landscape productivity based on photosynthetic activity, and 3) land cover were used as inputs to create an environmental regionalization of the over 10 million km2 of Canada’s terrestrial land base. The outcomes of this clustering consists of three main outputs. An initial clustering of 100 classes was generated using a two-stage multivariate classification process. Next, an agglomerative hierarchy using a log-likelihood distance measure was applied to create a 40 and then a 14 class regionalization, aimed to meaningfully group ecologically similar components of Canada's terrestrial landscape. For more information (including a graphical illustration of the cluster hierarchy) and to cite this data please use: Coops, N.C., Wulder, M.A., Iwanicka, D. 2009. An environmental domain classification of Canada using earth observation data for biodiversity assessment. Ecological Informatics, Vol. 4, No. 1, Pp. 8-22, DOI: https://doi.org/10.1016/j.ecoinf.2008.09.005. ( Coops et al. 2009).

This dataset provides a Canada-wide map of vegetation height and the delineation of the northern forest limit. It is developed within the framework of Canada’s National Terrestrial Ecosystem Monitoring System (NTEMS). Vegetation height estimates were derived from ICESat-2 LiDAR observations, integrated with Landsat time series and topographic variables to model spatial patterns. The northern forest limit represents the transition between boreal forest and tundra, an ecologically significant zone for monitoring climate change impacts and biodiversity. Vegetation height was modeled for six time-periods including 1985-1995, 1990-2000, 1995-2005, 2000-2010, 2005-2015 and 2010-2021. Predictions for each time period represent the median conditions for that period. Predictions of height and the probability of canopy presence were generated using Random Forests models trained on spaceborne-lidar data collected by ICESat-2 from 2019-2021 and overlapping Landsat satellite imagery from 2010-2021. These Random Forests models were then applied to the entire archive of Landsat imagery, representing a period of ~35 years. This dataset provides spatially explicit prediction of vegetation height (m) along the Canadian northern forest limit at 30 m spatial resolution. Pixels with a low (< 50 %) probability of containing a vegetation canopy have been assigned a height of 0 m. The science and methods for this dataset were the result of a collaboration between the Canadian Forest Service of Natural Resources Canada, partnered with the Integrated Remote Sensing Studio (IRSS) in the Faculty of Forestry at the University of British Columbia. When using this data, please cite: Travers-Smith, H., Coops, N. C., Mulverhill, C., Wulder, M. A., Ignace, D., Lantz, T. C. (2024). Mapping vegetation height and identifying the northern forest limit across Canada using ICESat-2, Landsat time series and topographic data. Remote Sensing of Environment, 305, 114097. https://doi.org/10.1016/j.rse.2024.114097 (Travers-Smith et al. 2024). Additional details outlining application of the model to the time-series of Landsat data can be found here: Travers-Smith, H., Coops, N., Mulverhill, C., Wulder, M. A., Lantz, T. C., Ignace, D. (2025). Satellite observations reveal stable forest limits and shrub expansion across the Canadian forest-tundra ecotone. Environmental Research Letters, 20(10). https://doi.org/10.1088/1748-9326/adfc7f (Travers-Smith et al. 2025).

Canada's National Forest Inventory (NFI) sampling program is designed to support reporting on forests at the national scale. On the other hand, continuous maps of forest attributes are required to support strategic analyses of regional policy and management issues. We have therefore produced maps covering 4.03 × 106 km2 of inventoried forest area for the 2001 base year using standardised observations from the NFI photo plots (PP) as reference data. We used the k nearest neighbours (kNN) method with 26 geospatial data layers including MODIS spectral data and climatic and topographic variables to produce maps of 127 forest attributes at a 250 × 250 m resolution. The stand-level attributes include land cover, structure, and tree species relative abundance. In this article, we report only on total live aboveground tree biomass, with all other attributes covered in the supplementary data (http://nrcresearchpress.com/doi/suppl/10.1139/cjfr-2013-0401). In general, deviations in predicted pixel-level values from those in a PP validation set are greater in mountainous regions and in areas with either low biomass or sparse PP sampling. Predicted pixel-level values are overestimated at small observed values and underestimated at large ones. Accuracy measures are improved through the spatial aggregation of pixels to 1 km2 and beyond. Overall, these new products provide unique baseline information for strategic-level analyses of forests (https://nfi.nfis.org) Collection: - **[Canada's National Forest Inventory (NFI) 2006](https://open.canada.ca/data/en/dataset/e2fadaeb-3106-4111-9d1c-f9791d83fbf4)**

The probability (likelihood) of ice freeze days, the number of days in the forecast period with a minimum temperature below the frost temperature, -5°C for herbaceous crops over the non-growing season (ifd_herb_nogrow_prob). Week 1 and week 2 forecasted probability is available daily from November 1 to March 31. Week 3 and week 4 forecasted probability is available weekly (Thursday) from November 1 to March 31. Over-wintering crops are biennial and perennial field crops such as herbaceous plants (strawberry, alfalfa, timothy, and many other forage crops) and woody fruit trees (apple, pear, peach, cherry, plum, apricot, chestnut, pecan, grape, etc.). These crops normally grow and develop in the growing season and become dormant in the non-growing season. However, extreme weather and climate events such as cold waves in the growing season and ice freezing events during the winter are a major constraint for their success of production and survival in Canada. The winter survival of these plants depends largely on agrometeorological conditions from late autumn to early spring, especially ice-freezing damage during the winter season. The optimum temperature for such crops is 25°C. Agriculture and Agri-Food Canada (AAFC) and Environment and Climate Change Canada (ECCC) have together developed a suite of extreme agrometeorological indices based on four main categories of weather factors: temperature, precipitation, heat, and wind. The extreme weather indices are intended as short-term prediction tools and generated using ECCC’s medium range forecasts to create a weekly index product on a daily basis.

Frost free days are the number of days in the forecast period with a minimum temperature above the frost temperature; the temperature at which frost damage occurs. This temperature is 0°C for warm season crops (ffd_warm). Week 1 and week 2 forecasted index is available daily from April 1 to October 31. Week 3 and week 4 forecasted index is available weekly (Thursday) from April 1 to October 31. Warm season crops require a relatively warm temperature condition. Typical examples include bean, soybean, corn and sweet potato. They normally grow during the summer season and early fall, then ripen in late fall in southern Canada only. Other agricultural regions in Canada do not always experience sufficiently long growing seasons for these plants to achieve maturity. The optimum temperature for such crops is 30°C. Agriculture and Agri-Food Canada (AAFC) and Environment and Climate Change Canada (ECCC) have together developed a suite of extreme agrometeorological indices based on four main categories of weather factors: temperature, precipitation, heat, and wind. The extreme weather indices are intended as short-term prediction tools and generated using ECCC’s medium range forecasts to create a weekly index product on a daily and weekly basis.

Heat Wave Days are the number of days in the forecast period with a maximum temperature above the cardinal maximum temperature, the temperature at which crop growth ceases. This temperature is 35°C for warm season crops (dhw_warm). Week 1 and week 2 forecasted index is available daily from April 1 to October 31. Week 3 and week 4 forecasted index is available weekly (Thursday) from April 1 to October 31. Warm season crops require a relatively warm temperature condition. Typical examples include bean, soybean, corn and sweet potato. They normally grow during the summer season and early fall, then ripen in late fall in southern Canada only. Other agricultural regions in Canada do not always experience sufficiently long growing seasons for these plants to achieve maturity. The optimum temperature for such crops is 30°C. Agriculture and Agri-Food Canada (AAFC) and Environment and Climate Change Canada (ECCC) have together developed a suite of extreme agrometeorological indices based on four main categories of weather factors: temperature, precipitation, heat, and wind. The extreme weather indices are intended as short-term prediction tools and generated using ECCC’s medium range forecasts to create a weekly index product on a daily and weekly basis.

The Probability (likelihood) of cool wave days for cool season/overwintering crops occurring Cool Wave Days are the number of days in the forecast period with a minimum temperature below the cardinal minimum temperature, the lowest temperature at which crop growth will begin (dcw_cool_prob). This temperature is 5°C for cool season crops. Week 1 and week 2 forecasted probability is available daily from April 1 to October 31. Week 3 and week 4 forecasted probability is available weekly (Thursday) from April 1 to October 31. Cool season crops require a relatively low temperature condition. Typical examples include wheat, barley, canola, oat, rye, pea, and potato. They normally grow in late spring and summer, and mature between the end of summer and early fall in the southern agricultural areas of Canada. The optimum temperature for such crops is 25°C. Agriculture and Agri-Food Canada (AAFC) and Environment and Climate Change Canada (ECCC) have together developed a suite of extreme agrometeorological indices based on four main categories of weather factors: temperature, precipitation, heat, and wind. The extreme weather indices are intended as short-term prediction tools and generated using ECCC’s medium range forecasts to create a weekly index product on a daily and weekly basis.

Dominant Species Map 2015 The data represent dominant tree species for British Columbia forests in 2015, are based upon Landsat data and modeling, with results mapped at 30 m spatial resolution. It is developed within the framework of Canada’s National Terrestrial Ecosystem Monitoring System (NTEMS). The map was generated with the Random Forests classifier that used predictor variables derived from Landsat time series including surface reflectance, land cover, forest disturbance, and forest structure, and ancillary variables describing the topography and position. Training and validation samples were derived from the Vegetation Resources Inventory (VRI), from a pool of polygons with homogeneous internal conditions and with low discrepancies with the remotely sensed predictions. Local models were applied over 100x100 km tiles that considered training samples from the 5x5 neighbouring tiles to avoid edge effects. An overall accuracy of 72% was found for the species which occupy 80% of the forested areas. Satellite data and modeling have demonstrated the capacity for up-to-date, wall-to-wall, forest attribute maps at sub-stand level for British Columbia, Canada. BC Species Likelihood 2015 The tree species class membership likelihood distribution data included in this product focused on the province of British Columbia, based upon Landsat data and modeling, with results mapped at 30 m spatial resolution. The data represent tree species class membership likelihood in 2015. The map was generated with the Random Forests classifier that used predictor variables derived from Landsat time series including surface reflectance, land cover, forest disturbance, and forest structure, and ancillary variables describing the topography and position. Training and validation samples were derived from the Vegetation Resources Inventory (VRI) selecting from a stratified pool of polygons with homogeneous internal conditions and with low discrepancies when related to remotely sensed information. Local models were applied over 100x100 km tiles that, to avoid edge effects, considered training samples from the 5x5 neighbouring tiles. An overall accuracy of 72% was found for the species which occupy 80% of the forested areas. As an element of the mapping process, we also obtain the votes received for each class by the Random Forest models. The votes can be understood as analogous to class membership likelihoods, providing enriched information on land cover class uncertainty for use in modeling. Tree species class membership likelihoods lower than 5% have been masked and converted to zero. When using this data, please cite as: Shang, C., Coops, N.C., Wulder, M.A., White, J.C., Hermosilla, T., 2020. Update and spatial extension of strategic forest inventories using time series remote sensing and modeling. International Journal of Applied Earth Observation and Geoinformation 84, 101956. DOI: 10.1016/j.jag.2019.101956 ( Shang et al. 2020).