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)**

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)**

Forest Lorey's Height 2015 Lorey's mean height. It is developed within the framework of Canada’s National Terrestrial Ecosystem Monitoring System (NTEMS). Average height of trees weighted by their basal area (m). Products relating the structure of Canada's forested ecosystems have been generated and made openly accessible. The shared products are based upon peer-reviewed science and relate aspects of forest structure including: (i) metrics calculated directly from the lidar point cloud with heights normalized to heights above the ground surface (e.g., canopy cover, height), and (ii) modelled inventory attributes, derived using an area-based approach generated by using co-located ground plot and ALS data (e.g., volume, biomass). Forest structure estimates were generated by combining information from lidar plots (Wulder et al. 2012) with Landsat pixel-based composites (White et al. 2014; Hermosilla et al. 2016) using a nearest neighbour imputation approach with a Random Forests-based distance metric. These products were generated for strategic-level forest monitoring information needs and are not intended to support operational-level forest management. All products have a spatial resolution of 30 m. For a detailed description of the data, methods applied, and accuracy assessment results see Matasci et al. (2018). When using this data, please cite as follows: Matasci, G., Hermosilla, T., Wulder, M.A., White, J.C., Coops, N.C., Hobart, G.W., Bolton, D.K., Tompalski, P., Bater, C.W., 2018b. Three decades of forest structural dynamics over Canada's forested ecosystems using Landsat time-series and lidar plots. Remote Sensing of Environment 216, 697-714. Matasci et al. 2018) Geographic extent: Canada's forested ecosystems (~ 650 Mha) Time period: 1985–2011

This dataset provides wall-to-wall maps of forest structure across Canada's 650 million hectare forested ecosystems for the year 2022, generated at a spatial resolution of 30 m. It is developed within the framework of Canada’s National Terrestrial Ecosystem Monitoring System (NTEMS). Structure estimates include key attributes such as canopy height, canopy cover, and aboveground biomass, derived using a combination of airborne lidar and Landsat-based spectral composites. Structure models were trained using the - lidar-plot framework - (Wulder et al. 2012), which integrates co-located airborne lidar data and ground plot measurements with Landsat time-series composites (Hermosilla et al. 2016). A Nearest Neighbour imputation approach was applied to estimate structural attributes across the full extent of Canada's forested area. These nationally consistent products are intended to support strategic-level forest monitoring and assessment and are not designed for operational forest management. For further details on the methods, accuracy assessment, and source data, see Matasci et al. (2018). Matasci, G., Hermosilla, T., Wulder, M.A., White, J.C., Coops, N.C., Hobart, G.W., Bolton, D.K., Tompalski, P., Bater, C.W., 2018. Three decades of forest structural dynamics over Canada's forested ecosystems using Landsat time-series and lidar plots. Remote Sensing of Environment, 216, 697-714. https://doi.org/10.1016/j.rse.2018.07.024 (Matasci et al. 2018)

Forest Basal Area 2015 Cross-sectional area of tree stems at breast height. It is developed within the framework of Canada’s National Terrestrial Ecosystem Monitoring System (NTEMS). The sum of the cross-sectional area (i.e. basal area) of each tree in square metres in a plot, divided by the area of the plot (units = m2ha). Products relating the structure of Canada's forested ecosystems have been generated and made openly accessible. The shared products are based upon peer-reviewed science and relate aspects of forest structure including: (i) metrics calculated directly from the lidar point cloud with heights normalized to heights above the ground surface (e.g., canopy cover, height), and (ii) modelled inventory attributes, derived using an area-based approach generated by using co-located ground plot and ALS data (e.g., volume, biomass). Forest structure estimates were generated by combining information from lidar plots (Wulder et al. 2012) with Landsat pixel-based composites (White et al. 2014; Hermosilla et al. 2016) using a nearest neighbour imputation approach with a Random Forests-based distance metric. These products were generated for strategic-level forest monitoring information needs and are not intended to support operational-level forest management. All products have a spatial resolution of 30 m. For a detailed description of the data, methods applied, and accuracy assessment results see Matasci et al. (2018). When using this data, please cite as follows: Matasci, G., Hermosilla, T., Wulder, M.A., White, J.C., Coops, N.C., Hobart, G.W., Bolton, D.K., Tompalski, P., Bater, C.W., 2018b. Three decades of forest structural dynamics over Canada's forested ecosystems using Landsat time-series and lidar plots. Remote Sensing of Environment 216, 697-714. Matasci et al. 2018) Geographic extent: Canada's forested ecosystems (~ 650 Mha) Time period: 1985–2011

Forest Elevation(Ht) Mean 2015 Mean height of lidar first returns (m). Represents the mean canopy height. It is developed within the framework of Canada’s National Terrestrial Ecosystem Monitoring System (NTEMS). Products relating the structure of Canada's forested ecosystems have been generated and made openly accessible. The shared products are based upon peer-reviewed science and relate aspects of forest structure including: (i) metrics calculated directly from the lidar point cloud with heights normalized to heights above the ground surface (e.g., canopy cover, height), and (ii) modelled inventory attributes, derived using an area-based approach generated by using co-located ground plot and ALS data (e.g., volume, biomass). Forest structure estimates were generated by combining information from lidar plots (Wulder et al. 2012) with Landsat pixel-based composites (White et al. 2014; Hermosilla et al. 2016) using a nearest neighbour imputation approach with a Random Forests-based distance metric. These products were generated for strategic-level forest monitoring information needs and are not intended to support operational-level forest management. All products have a spatial resolution of 30 m. For a detailed description of the data, methods applied, and accuracy assessment results see Matasci et al. (2018). When using this data, please cite as follows: Matasci, G., Hermosilla, T., Wulder, M.A., White, J.C., Coops, N.C., Hobart, G.W., Bolton, D.K., Tompalski, P., Bater, C.W., 2018b. Three decades of forest structural dynamics over Canada's forested ecosystems using Landsat time-series and lidar plots. Remote Sensing of Environment 216, 697-714. Matasci et al. 2018) Wulder et al. 2018) Geographic extent: Canada's forested ecosystems (~ 650 Mha) Time period: 1985–2011

The probability (likelihood) of ice freeze days, the number of days in the forecast period with a minimum temperature below the frost temperature, -10°C for woody crops over the non-growing season (ifd_wood_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. 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.

Landsat-derived forest age for Canada 2022 Satellite-based forest age map for 2022 across Canada's forested ecozones at a 30-m spatial resolution, developed within the framework of Canada’s National Terrestrial Ecosystem Monitoring System (NTEMS). Remotely sensed data from Landsat (disturbances, surface reflectance composites, forest structure) and MODIS (Gross Primary Production) are utilized to determine age. Age can be determined where disturbance can be identified directly (disturbance approach) or inferred using spectral information (recovery approach) or using inverted allometric equations to model age where there is no evidence of disturbance (allometric approach). The disturbance approach is based upon satellite data and mapped changes and is the most accurate. The recovery approach also avails upon satellite data plus logic regarding forest succession, with an accuracy that is greater than pure modeling. Given the lack of widespread recent disturbance over Canada's forests, the allometric approach is required over the greatest area (86.6%). Using information regarding realized heights and growth and yield modeling, ages are estimated where none are otherwise possible. Trees of all ages are mapped, with trees >150 years old combined in an - old tree - category. See Maltman et al. (2023) for an overview of the methods, data, image processing, as well as information on agreement assessment using Canada's National Inventory (NFI). Maltman, J.C., Hermosilla, T., Wulder, M.A., Coops, N.C., White, J.C., 2023. Estimating and mapping forest age across Canada's forested ecosystems. Remote Sensing of Environment 290, 113529. ( Maltman et al. 2023).

Wall-to-wall map of water bodies across Canada's forested ecosystems for the year 2022, derived from the "water" class of the annual Virtual Land Cover of Engine (VLCE) product. It is developed within the framework of Canada’s National Terrestrial Ecosystem Monitoring System (NTEMS). The VLCE maps are based on Landsat image time-series composites and represent annual land cover classifications from 1984 to 2022 at a spatial resolution of 30 m. The classification process integrates forest change information and ancillary topographic and hydrologic variables, applying a regional modeling framework based on a 150x150 km tiling system ( Hermosilla et al., 2022). Training data are drawn from multiple land cover sources and selected proportionally to land cover distributions using a distance-weighted approach. Classifications are refined over time using a Hidden Markov Model to ensure consistency and reduce classification noise between years. Hermosilla, T., Wulder, M.A., White, J.C., Coops, N.C. 2022. Land cover classification in an era of big and open data: Optimizing localized implementation and training data selection to improve mapping outcomes. Remote Sensing of Environment. 268, 112780. https://doi.org/10.1016/j.rse.2021.112780. ( Hermosilla et al., 2022) Hermosilla, T., Wulder, M.A., White, J.C., Coops, N.C., Hobart, G.W. 2018. Disturbance-Informed Annual Land Cover Classification Maps of Canada's Forested Ecosystems for a 29-Year Landsat Time Series. Canadian Journal of Remote Sensing. 44(1) 67-87. https://doi.org/10.1080/07038992.2018.1437719.( Hermosilla et al., 2018)