In this dataset, we share maps of annual dominant tree species (also known as leading tree species) from 1984-2022 covering the entirety of Canada's 650 Mha forested ecosystems using Landsat time-series imagery at a 30-m spatial resolution. It is developed within the framework of Canada’s National Terrestrial Ecosystem Monitoring System (NTEMS). Classifications are based on regionally representative Random Forests model using local training samples from Canada's National Forest Inventory (Hermosilla et al., 2024). Descriptive metrics provide information on spectral, geographic, climatic, and topographic characteristics. Initial annual tree species classifications were subjected to a time series post-classification process using the forward-backward Hidden Markov Model to improve the temporal consistency of tree species transitions within the time series. Assessment of the annual species maps using independent validation data resulted in an overall accuracy of 86.1% ± 0.14% (95%-confidence interval). These data allow consistent comparison of trends and rates of change in tree species composition nationally and across regions using a common time frame, spatial resolution, and analytical approach. Hermosilla, T., Wulder, M.A., White, J.C., Coops, N.C., Bater, C.W., Hobart, G.W., 2024. Characterizing long-term tree species dynamics in Canada's forested ecosystems using annual time series remote sensing data. Forest Ecology and Management, 122313. https://doi.org/10.1016/j.foreco.2024.122313 (Hermosilla et al. 2024)

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)

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. It is 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).

High-resolution false-color Landsat image composite of Canada's forested ecosystems (2022). This national image product represents the Composite to Change (C2C) proxy composite image derived from thousands of Landsat images acquired between July 1 and August 30, 2022. It is developed within the framework of Canada’s National Terrestrial Ecosystem Monitoring System (NTEMS). The overall process followed is described in (Hermosilla et al. 2016 ) with details on the generation of gap-free surface reflectance composites in ( Hermosilla et al. 2015). Following the motivation and rationale presented in White et al. (White et al. 2014), Landsat imagery is subjected to a series of processing steps to remove clouds and shadows as well as haze and other unwanted atmospheric effects. Year-on-year time series of Landsat imagery are interrogated to avoid missing values, and to ensure exhaustive spatial coverage of the national surface reflectance composites. False-colour 3-channel image (bands: shortwave infrared, SWIR1; near infrared; red) When using these data, please cite as: Hermosilla, T., M.A. Wulder, J.C. White, N.C. Coops, G.W. Hobart, L.B. Campbell, 2016. Mass data processing of time series Landsat imagery: pixels to data products for forest monitoring. International Journal of Digital Earth 9(11), 1035-1054 (Hermosilla et al. 2016 ).

Fish Habitat Assessment Output: 4 of 16 High Water Level (75.4m ASL) - Nursery Habitat - Low Vegetation Association Species (All Temperature Windows) Habitat suitability was assessed for the Bay of Quinte Area of Concern, at a 3 m grid resolution, using the Habitat Ecosystem Assessment Tool (HEAT), temperature algorithms, vegetation models, and water level input. Habitat classifications were based on three variables: depth (elevation), vegetation, and substrate; and modified by temperature suitabilities. The final suitability maps were based on documented habitat and temperature associations for the fish in the area. Different life stages (spawning requirements, nursery habitat, adult habitat) were modeled for the years of 1972-2011. Suitability values were scaled from 0 (not suitable) to 1 (highly suitable) and converted to suitability classes of very low, low, medium, and high. The final maps for each guild – life stage combination are maximum suitability values from the 39-year period modelled.

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) Related Products (16): - **[Poplars, Aspens and Cottonwoods (Genus Populus) in Canada 2006](https://ouvert.canada.ca/data/en/dataset/08620b3f-0bda-46f2-968d-e47d5a6032de)** - **[Birches (Genus Betula) in Canada 2006](https://ouvert.canada.ca/data/en/dataset/1410d784-ffde-43c8-a816-14cdfa0aa9c4)** - **[Treed land in Canada 2006](https://ouvert.canada.ca/data/en/dataset/1f1806b9-3927-496c-8c91-8789809f4472)** - **[Merchantable forest volume in Canada 2006](https://ouvert.canada.ca/data/en/dataset/2b3569c6-ff95-40a5-a958-dc68e3aa558b)** - **[Needle-leaved species in Canada 2006](https://ouvert.canada.ca/data/en/dataset/39ffee48-f89b-4b65-af03-58a706bac7a1)** - **[Forest height in Canada 2006](https://ouvert.canada.ca/data/en/dataset/3b860e37-32e6-4f47-a423-a7519ffa4429)** - **[Total live above-ground biomass in Canada 2006](https://ouvert.canada.ca/data/en/dataset/53af4b0e-015b-405e-8de8-e7eb4498eda1)** - **[Total forest volume in Canada 2006](https://ouvert.canada.ca/data/en/dataset/5b6b60d5-8299-45d4-8bd2-c274e75bc115)** - **[Spruces (Genus Picea) in Canada 2006](https://ouvert.canada.ca/data/en/dataset/72af5640-bb51-4c7a-8f41-6b71227a598f)** - **[Tree Crown Closure in Canada 2006](https://ouvert.canada.ca/data/en/dataset/a1510fe3-8ef0-4130-9ee4-8a7ea1f9a22d)** - **[Forest Composition across Canada 2006](https://ouvert.canada.ca/data/en/dataset/a42bd5d6-83a7-4fb7-a257-389dcf7ea48d)** - **[True Firs (Genus Abies) in Canada 2006](https://ouvert.canada.ca/data/en/dataset/d845f357-e2b5-494c-821f-064dba664427)** - **[Hemlocks (Genus Tsuga) in Canada 2006](https://ouvert.canada.ca/data/en/dataset/e676e5ab-b709-46ba-b471-4e982dce0c07)** - **[Broad-leaved species in Canada 2006](https://ouvert.canada.ca/data/en/dataset/e7b9b34a-70f8-47c6-9498-94361b9febbf)** - **[Maples (Genus Acer) in Canada 2006](https://ouvert.canada.ca/data/en/dataset/ed296d14-b222-45e7-9dcb-0ca8015207ad)** - **[Cedars (Genus Thuja) in Canada 2006](https://ouvert.canada.ca/data/en/dataset/fc8bb212-9ffe-447f-9152-e26baff7a735)**

Fish Habitat Assessment Output: 11 of 16 Average Water Level (75.0m ASL) - Nursery Habitat - High Vegetation Association Species (All Temperature Windows) Habitat suitability was assessed for the Bay of Quinte Area of Concern, at a 3 m grid resolution, using the Habitat Ecosystem Assessment Tool (HEAT), temperature algorithms, vegetation models, and water level input. Habitat classifications were based on three variables: depth (elevation), vegetation, and substrate; and modified by temperature suitabilities. The final suitability maps were based on documented habitat and temperature associations for the fish in the area. Different life stages (spawning requirements, nursery habitat, adult habitat) were modeled for the years of 1972-2011. Suitability values were scaled from 0 (not suitable) to 1 (highly suitable) and converted to suitability classes of very low, low, medium, and high. The final maps for each guild – life stage combination are maximum suitability values from the 39-year period modelled.

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 Elevation(Ht) Stddev 2015 Standard deviation of height of lidar first returns (m). Represents the variability in canopy heights. 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) Geographic extent: Canada's forested ecosystems (~ 650 Mha) Time period: 1985–2011