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

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

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

High-resolution annual forest land cover maps for Canada's forested ecosystems (1984-2022). The annual time series of forest land cover maps are national in scope (entire 650 million hectare forested ecosystem) and represent a wall-to-wall land cover characterization yearly from 1984 to 2022. These time-series land cover maps were produced from annual time-series of Landsat image composites, forest change information, and ancillary topographic and hydrologic data following the framework described in Hermosilla et al. (2022), which builds upon the approach introduced in Hermosilla et al. (2018). The methodological innovations included (i) a refined training pool derived from existing land cover products using airborne and spaceborne measures of forest structure; (ii) selection of training samples proportionally to the land cover distribution using a distance-weighted approach; and (iii) generation of regional classification models using a 150x150 km tiling system. Maps are post-processed using disturbance information to ensure logical class transitions over time using a Hidden Markov Model. Hidden Markov Models assess individual year class likelihoods to reduce variability and possible noise in year-on-year class assignments (for instances when class likelihoods are similar). 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. Vol. 268, No. 112780. https://doi.org/10.1016/j.rse.2021.112780. ( Hermosilla et al. 2022) Hermosilla, T., M.A. Wulder, J.C. White, N.C. Coops, G. W. Hobart, (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.DOI: 10.1080/07038992.2018.1437719 ( Hermosilla et al. 2018).

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)

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 data set provides forest age map products at 1-km resolution for Canada and the United States (U.S.A.). These continental forest age maps were compiled from forest inventory data, historical fire data, optical satellite data, and the images from the NASA Landsat Ecosystem Disturbance Adaptive Processing System (LEDAPS) project. These input data products have various sources and creation dates as described in the source paper by Pan et al. (2011). Canadian maps were produced with data available through 2004 and U.S.A. maps with data available through 2006. A supplementary map of the standard deviations for age estimates was developed for quantifying uncertainty.Note that the Pan et al. (2011) paper is included as a companion file with this data set and was the source of descriptions in the guide.Forest age, implicitly reflecting the past disturbance legacy, is a simple and direct surrogate for the time since disturbance and may be used in various forest carbon analyses that concern the impact of disturbances. By combining geographic information about forest age with estimated carbon dynamics by forest type, it is possible to conduct a simple but powerful analysis of the net CO2 uptake by forests, and the potential for increasing (or decreasing) this rate as a result of direct human intervention in the disturbance/age status.

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)

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

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)