This dataset provides treed area dynamics across Canada's 650 Mha forested ecosystems from 1984 to 2022, derived from Landsat-based annual land cover layers at a 30-m spatial resolution. It is developed within the framework of Canada’s National Terrestrial Ecosystem Monitoring System (NTEMS). This dataset identifies areas that remained treed, transitioned to treed (newly treed), or transitioned to other cover that is not treed vegetation (was-treed). The data enable national and regional assessments of long-term changes in treed area, capturing trends in treed area, post-disturbance recovery, and shifts in forest extent. When using this data, please cite as: Hermosilla, T., Wulder, M.A., White, J.C., Bater, C.W., Baral, S.K., Leach, J.A., 2025. Expansion of treed area over Canada’s forested ecosystems: spatial and temporal trends. Forestry: An International Journal of Forest Research 98(5) 786-799. https://doi.org/10.1093/forestry/cpaf015. (Hermosilla et al. 2025)

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)

Harvest changes occurred from 1985 to 2022 displaying the year of greatest harvest disturbance. It is developed within the framework of Canada’s National Terrestrial Ecosystem Monitoring System (NTEMS). The information outcomes represent 38 years of harvest activity in Canada's forests, derived from a single, consistent, spatially explicit data source in a fully automated manner. Time series of Landsat data with 30 m spatial resolution were used to characterize national trends in stand replacing forest disturbances caused by harvest for the period 1985-2022 for Canada's 650-million-hectare forested ecosystems. When using this 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. https://doi.org/10.1080/17538947.2016.1187673 ( Hermosilla et al. 2016). See references below for an overview on the data processing, metric calculation, change attribution, and time series change detection methods applied, as well as information on independent accuracy assessment of the data. Hermosilla, T., Wulder, M. A., White, J. C., Coops, N.C., Hobart, G.W., (2015). An integrated Landsat time series protocol for change detection and generation of annual gap-free surface reflectance composites. Remote Sensing of Environment 158, 220-234. ( Hermosilla et al. 2015a). Hermosilla, T., Wulder, M.A., White, J.C., Coops, N.C., Hobart, G.W., (2015). Regional detection, characterization, and attribution of annual forest change from 1984 to 2012 using Landsat-derived time-series metrics. Remote Sensing of Environment 170, 121-132. ( Hermosilla et al. 2015b). Hermosilla, T., M.A. Wulder, J.C. White, N.C. Coops, G. W. Hobart, (2017). Updating Landsat time series of surface-reflectance composites and forest change products with new observations. International Journal of Applied Earth Observation and Geoinformation. 63,104-111. https://doi.org/10.1016/j.jag.2017.07.013 (Hermosilla et al. 2017)

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)

Wildfire change year 1985-2022. Wildfire changes occurred from 1985 to 2022 displaying the year of greatest wildfire disturbance. It is developed within the framework of Canada’s National Terrestrial Ecosystem Monitoring System (NTEMS). The information outcomes represent 38 years of wildfires in Canada's forests, derived from a single, consistent, spatially explicit data source in a fully automated manner. Time series of Landsat data with 30 m spatial resolution were used to characterize national trends in stand replacing forest disturbances caused by wildfire for the period 1985-2022 for Canada's 650-million-hectare forested ecosystems. When using this 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. https://doi.org/10.1080/17538947.2016.1187673 (Hermosilla et al. 2016). See references below for an overview on the data processing, metric calculation, change attribution, and time series change detection methods applied, as well as information on independent accuracy assessment of the data.. Hermosilla, T., Wulder, M. A., White, J. C., Coops, N.C., Hobart, G.W., (2015). An integrated Landsat time series protocol for change detection and generation of annual gap-free surface reflectance composites. Remote Sensing of Environment 158, 220-234. https://doi.org/10.1016/j.rse.2014.11.005 (Hermosilla et al. 2015a). Hermosilla, T., Wulder, M.A., White, J.C., Coops, N.C., Hobart, G.W., (2015). Regional detection, characterization, and attribution of annual forest change from 1984 to 2012 using Landsat-derived time-series metrics. Remote Sensing of Environment 170, 121-132. https://doi.org/10.1016/j.rse.2015.09.004 (Hermosilla et al. 2015b). Hermosilla, T., M.A. Wulder, J.C. White, N.C. Coops, G. W. Hobart, (2017). Updating Landsat time series of surface-reflectance composites and forest change products with new observations. International Journal of Applied Earth Observation and Geoinformation. 63,104-111. https://doi.org/10.1016/j.jag.2017.07.013 (Hermosilla et al. 2017).

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. 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. This product was developed within the framework of Canada’s National Terrestrial Ecosystem Monitoring System (NTEMS). 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).

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

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)

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

The wetland year count data included in this product is national in scope (entire forested ecosystem) and represents a wall to wall wetland characterization for 1984-2016 (Wulder et al. 2018). It is developed within the framework of Canada’s National Terrestrial Ecosystem Monitoring System (NTEMS). This product was generated using both annual gap free composite reflectance images and annual forest change maps following the Virtual Land Cover Engine (VLCE) process (see Hermosilla et al. 2018), over the 650 million ha forested ecosystems of Canada. Elements of the VLCE classification approach are inclusion of disturbance information in the processes as well as ensuring class transitions over time are logical. Further, a Hidden Markov Model is implemented to assess individual year class likelihoods to reduce variability and possible noise in year-on-year class assignments (for instances when class likelihoods are similar). The values can range from 0 to 33 denoting the number of years between 1984 and 2016 that a pixel was classified as wetland or wetland-treed in the VLCE data cube. For an overview on the data, image processing, and time series change detection methods applied, as well as information on independent accuracy assessment of the data, see Hermosilla et al. (2016; http://www.tandfonline.com/doi/full/10.1080/17538947.2016.1187673). A detailed description of the VLCE process and the subsequently generated land cover product, including an accuracy assessment, please see Hermosilla et al. (2018). The focused wetland analyses can be found described in Wulder et al (2018). Geographic extent: Canada's forested ecosystems (~ 650 Mha) Time period: 1985–2011