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

High-resolution map of leading tree species distribution for Canada’s forested ecosystems (2019). Leading tree species map produced from a 2019 Landsat image composite, geographic and climate data, elevation derivatives, and remote sensing derived phenology following the framework described in Hermosilla et al. (2022). Regional classification models were generated based on Canada’s National Forest Inventory using a 150x150 km tiling system. The leading tree species are defined by representing the most voted tree species from the Random Forests classification models (i.e. the class with the highest class membership probability).The data represents leading tree species of Canada's forested ecosystems in 2019. An image compositing window of August 1 ± 30 days was used to generate the best-available-pixel (BAP) image composites utilized as source data for the classification. The science and methods developed to generate the information outcomes shown here, that track and characterize the history of Canada’s forests, were led by Canadian Forest Service of Natural Resources Canada, developed within the framework of Canada’s National Terrestrial Ecosystem Monitoring System (NTEMS), partnered with the University of British Columbia, augmented by processing capacity from Digital Research Alliance of Canada. For an overview on the data, image processing, and methods applied, as well as information on independent accuracy assessment of the data, see Hermosilla et al. (2022) https://doi.org/10.1016/j.rse.2022.113276 When using this data, please cite as: Hermosilla, T., Bastyr, A., Coops, N.C., White, J.C., Wulder, M.A., 2022. Mapping the presence and distribution of tree species in Canada’s forested ecosystems. Remote Sensing of Environment 282, 113276.

Fish Habitat Assessment Output: 1 of 16 High Water Level (75.4m ASL) - Spawning 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.

Wildfire change magnitude 1985-2022. Spectral change magnitude for wildfires that occurred from 1985 to 2022. It is developed within the framework of Canada’s National Terrestrial Ecosystem Monitoring System (NTEMS). The wildfire change magnitude included in this product is expressed via differenced Normalized Burn Ratio (dNBR), computed as the variation between the spectral values before and after a given change event. Higher dNBR values are related to higher burn severity. 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. (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).

The maps show a multiyear ground deformation rate caused by small-scale deformation processes in Canada, measured in meters per year. Horizontal-east and vertical deformation components were computed from data acquired on ascending and descending orbits. This horizontal-east/vertical 2D decomposition is approximate and assumes constant viewing geometry and the absence of horizontal-north deformation. In the line-of-sight (LOS) map computed from ascending orbit data, a negative signal approximately corresponds to either subsidence or eastward motion, while a positive signal corresponds to uplift or westward motion. In the LOS map computed from descending orbit data, a negative signal approximately corresponds to either subsidence or westward motion, while a positive signal corresponds to uplift or eastward motion. In the horizontal-east map, a negative signal corresponds to westward motion, while a positive signal corresponds to eastward motion. In the vertical map, a negative signal indicates subsidence, while a positive signal indicates uplift. The maps were calculated from Sentinel-1 Synthetic Aperture Radar data collected between 2017 and 2024 during the snow-free season. Interferometric analysis of Sentinel-1 data was performed using GAMMA Software (https://www.gamma-rs.ch), and the long-term deformation rate was computed with the Multidimensional Small Baseline Subset (MSBAS) Software Version 10 (https://doi.org/10.1080/07038992.2024.2424753) at the Canada Centre for Mapping and Earth Observation, Natural Resources Canada. Long-wavelength signals caused by postglacial rebound and tectonic motion were filtered to enhance the visibility of small-scale deformation processes, such as those originating from landslides and mining. Field studies have confirmed only a few of these processes to date. The maps are expected to contain processing artifacts, which will be addressed in future work. References: Samsonov, S. V., & Feng, W. (2023). Deformation Retrievals for North America and Eurasia from Sentinel-1 DInSAR: Big Data Approach, Processing Methodology and Challenges. Canadian Journal of Remote Sensing, 49(1). https://doi.org/10.1080/07038992.2023.2247095 Samsonov, S. V. (2024). Multidimensional Small Baseline Subset (MSBAS) Software for Constrained and Unconstrained Deformation Analysis of Partially Coherent DInSAR and Speckle Offset Data. Canadian Journal of Remote Sensing, 50(1). https://doi.org/10.1080/07038992.2024.2424753 Limitation of Liability : The information contained on this website is provided on an “as is” basis and Natural Resources Canada makes no representations or warranties respecting the information, either expressed or implied, arising by law or otherwise, including but not limited to, effectiveness, completeness, accuracy or fitness for a particular purpose. Natural Resources Canada does not assume any liability in respect of any damage or loss based on the use of this website. In no event shall Natural Resources Canada be liable in any way for any direct, indirect, special, incidental, consequential, or other damages based on any use of this website or any other website to which this site is linked, including, without limitation, any lost profits or revenue or business interruption.

**Attention: there is a new version of this product (CanLaBS v2)** CanLaBS v2 can be found here: https:/doi.org/10.23687/2af751e7-79f9-4da8-9b45-14688818dca3 This data publication contains a set of files in which different variables related to fire burned severity (Canada Landsat Burned Severity, CanLaBS) were computed for all events in Canada between 1985 and 2015 as detected by the Canada Landsat Disturbance (CanLaD (Guindon et al. 2017 and 2018) product. Details on the creation of this product are available in Guindon et al. 2020 (https://doi.org/10.1139/cjfr-2020-0353) and in supplementary materials accompanying the publication. The current document is therefore a complement to the article and supplementary materials. The supplementary materials are referenced in the publication (cjfr-2020-0353suppla, cjfr-2020-0353supplb etc.). This is the first Canada-wide product that aims to promote nationwide research on fire severity by making available the data used in the article. The data is in the form of grids composed of pixels at a resolution of 30m. To simplify the distribution and manipulation of the data and considering that two or three fire occurrences within a given location is rare (respectively 2.3% and less than 0.01%), only the most recent fire data are considered in the final product. For these very rare cases, from 2015 to 1985, the most recent burned areas overlap the older data. Overlapping fire count can be found in layer “CanLaBS_Nbdisturb_v0”, multiple fire events in same areas have values equal to or greater than two. Landsat radiometric values for calculating the NBR index were derived from summer Landsat mosaics (July and August), for years 1984 to 2015 (Guindon et al. 2018). These mosaics were developed from individual USGS Landsat scenes with surface reflectance correction (Masek et al., 2006; Vermote et al., 2006). For each annual compound, the pixel with the less atmospheric opacity was selected. An algorithm was also developed to remove clouds that were not detected by the cloud masks provided with the USGS data. Here is a general description of the layers provided and a more technical description can be found in Table 1 (see "Ressources" section below): 1. NBR and dNBR. All these values are multiplied by 1000. The value of dNBR represents the value obtained for NBRpre - NBRpost. It is calculated for each pixel that was classified as a fire in CanLaD, according to the corrected year (see cjfr-2020-0353suppla). 2. Year of fire. The fire years detected in CanLaD (Guindon et al. 2018) was corrected using different fire databases, this layer contains the correct year. (see cjfr-2020-0353suppla) 3. Julian Days of the Fire, based on various high-resolution products. However, this variable is only available from 1989 onwards. 4. Presence of salvage logging one year after the fire. Classification of satellite images detecting scarified soils (see cjfr-2020-0353suppld). 5. Pre-fire forest attributes: Pre-fire forest attributes values were calculated for median mosaics, from 1985 to 2000. These attributes values were derived from NFI (national forest inventory) photo-plot attributes and were spatialized. Pre-fire attribute values were created to stratify the analyses (see cjfr-2020-0353supplc). The predicted variables are as follows: • Canopy density in percent. • Predicted living biomass in tonnes per hectare. • Percentage coniferous biomass proportion of total biomass. • Percentage hardwood biomass proportion of total biomass. • Percentage unknown species biomass proportion of total biomass. Note, as unknown species are found especially in northern areas, they are considered coniferous for the purpose of the article. 6. Missing remote sensing data, one year after the fire. The estimation of burned severity needs NBR data (NBRpost) in the next year after fire occurrences. NBRpost is available for 91% of the cases, but for the remaining 9%, no data were available due to the presence of clouds. For these cases, satellite data from the years following the fire were used with a regression radiometry correction. This gives values to missing data for year following the fire. This layer flags the areas that have derived data. The values of 1= one year after the fire (no regression), 2= two years after the fire (regression), 3= three years after the fire (regression) and 4= four years after the fire (no regression, set as missing data). (see cjfr-2020-0353supplb). 7. Areas with more than one fire disturbance between 1985 and 2015 (1=one single disturbance, 2=two or more, 3=three or more). ## Data citation: 1. Guindon, L., Villemaire P., Manka F., Dorion H. , Skakun R., St-Amant R., Gauthier S. : Canada Landsat Burned Severity (CanLaBS): a Canada-wide Landsat-based 30-m resolution product of burned severity since 1985 https://doi.org/10.23687/b1f61b7e-4ba6-4244-bc79-c1174f2f92cd 2. The creation, the validation and the limits of the CanLaBS product are describe in the text and supplementary material: Guindon, L., Gauthier, S., Manka, F., Parisien, MA, Whitman, E., Bernier, P., Beaudoin, A., Villemaire P., Skakun R. Trends in wildfire burn severity across Canada, 1985 to 2015 https://doi.org/10.1139/cjfr-2020-0353 ## References cited: 1. Guindon, L., Villemaire, P., St-Amant, R., Bernier, P.Y., Beaudoin, A., Caron, F., Bonucelli, M., and Dorion, H. 2017. Canada Landsat Disturbance (CanLaD): a Canada-wide Landsat-based 30m resolution product of fire and harvest detection and attribution since 1984. https://doi.org/10.23687/add1346b-f632-4eb9-a83d-a662b38655ad 2. Guindon, L., Bernier, P., Gauthier, S., Stinson, G., Villemaire, P., & Beaudoin, A. (2018). Missing forest cover gains in boreal forests explained. Ecosphere, 9(1), e02094. https://doi.org//10.1002/ecs2.2094 3. Masek, J.G., Vermote, E.F., Saleous N.E., Wolfe, R., Hall, F.G., Huemmrich, K.F., Gao, F., Kutler, J., and Lim, T-K. (2006). A Landsat surface reflectance dataset for North America, 1990–2000. IEEE Geoscience and Remote Sensing Letters 3(1):68-72. http://dx.doi.org/10.1109/LGRS.2005.857030. 4. Vermote, E., Justice, C., Claverie, M., & Franch, B. (2016). Preliminary analysis of the performance of the Landsat 8/OLI land surface reflectance product. Remote Sensing of Environment. http://dx.doi.org/10.1016/j.rse.2016.04.008.

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

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