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)

The data shared are spatially explicit projections of wildfire burn probability across Canada’s forested ecozones under multiple future climate scenarios at a 30-m spatial resolution. It is developed within the framework of Canada’s National Terrestrial Ecosystem Monitoring System (NTEMS). Four future climate scenarios were used to examine the spatiotemporal distribution of burn probability in the 21st century based on climate, vegetation, and topographic conditions ( Mulverhill et al. 2024). Projected burn probability is provided for four Shared Socioeconomic Pathways (SSP1-2.6, SSP2-4.5, SSP3-7.0, and SSP5-8.5) and four future time periods, including 2021-2040, 2041-2060, 2061-2080, and 2081-2100, along with a baseline period representing average climate conditions and burn probability between 1991 and 2020. Outputs represent the probability that the conditions (climate, vegetation, topography) of a given pixel resemble those of historically burned areas. All non-climate variables were held static; therefore, projections represent burn probability under future climate scenarios given contemporary (2020) forest conditions. When using this dataset, please cite Mulverhill et al. (2025), as below.Mulverhill, C., Coops, N. C., Wulder, M. A., Hermosilla, T., White, J. C., & Bater, C. W. (2025). Projected Future Changes in Burn Probability in Canada’s Forests and Communities Under Different Climate Change Scenarios. Canadian Journal of Remote Sensing, 51(1). https://doi.org/10.1080/07038992.2025.2560347(Mulverhill et al. 2025).For a detailed description of the source data and methods applied to the baseline period to enable the Mulverhill et al. (2025) projections, see:Mulverhill, C., Coops, N.C., Wulder, M.A., White, J.C., Hermosilla, T., and Bater, C.W. 2024. “Multidecadal mapping of status and trends in annual burn probability over Canada’s forested ecosystems.” ISPRS Journal of Photogrammetry and Remote Sensing, Vol. 209 pp. 279–295. https://doi.org/10.1016/j.isprsjprs.2024.02.006(Mulverhill et al. 2024).

The third generation of high resolution 10-m wetland inventory map of Canada, covering an approximate area of one billion hectares, was generated using multi-year (2016-2020), multi-source imagery (Sentinel-1, Sentinel-2, ALOS PALSAR-2, and SRTM) Earth Observation (EO) data as well as environmental features. Over 8800 wetland polygons were processed within an object-based random forest classification scheme on the Google Earth Engine cloud computing platform. The average overall accuracy of 90.5% is an increase of 4.7% over CWIM2.CWIM Versions:The Canadian Wetland Inventory Map (CWIM) is an extension of work started at Memorial University to produce a Newfoundland and Labrador wetland inventory during 2015-2018 which was significantly funded by Environment and Climate Change Canada. The first national CWIM was produced 2018-2019 as a collaboration between Memorial University, C-CORE, and Natural Resources Canada. Dr. Brian Brisco was instrumental in connecting ground truth from multiple sources to the project and providing guidance. Version 2 was produced in 2020 which included more training data and processing by Canada’s ecozones rather than provinces to take advantage of the commonality of landscape ecological features within ecozones to improve the accuracy. Version 3 produced in 2021 continued adding more data sources to further improve accuracy specifically an overestimation of wetland area as well as introducing a confidence map. Version 3A completed in 2022 updates only the arctic ecozones due to their relatively lower accuracy and added hydro-physiographic data layers. Currently work is underway to create a northern circumpolar wetland inventory map to be published in 2025.Paper on Newfoundland and Labrador Wetland Inventory:Mahdianpari, M.; Salehi, B.; Mohammadimanesh, F.; Homayouni, S.; Gill, E. The First Wetland Inventory Map of Newfoundland at a Spatial Resolution of 10 m Using Sentinel-1 and Sentinel-2 Data on the Google Earth Engine Cloud Computing Platform. Remote Sens. 2019, 11, 43. https://doi.org/10.3390/rs11010043Paper on CWIM1:Mahdianpari, M., Salehi, B., Mohammadimanesh, F., Brisco, B., Homayouni, S., Gill, E., … Bourgeau-Chavez, L. (2020). Big Data for a Big Country: The First Generation of Canadian Wetland Inventory Map at a Spatial Resolution of 10-m Using Sentinel-1 and Sentinel-2 Data on the Google Earth Engine Cloud Computing Platform. Canadian Journal of Remote Sensing, 46(1), 15–33. https://doi.org/10.1080/07038992.2019.1711366Paper on CWIM2:Mahdianpari, M., Brisco, B., Granger, J. E., Mohammadimanesh, F., Salehi, B., Banks, S., … Weng, Q. (2020). The Second Generation Canadian Wetland Inventory Map at 10 Meters Resolution Using Google Earth Engine. Canadian Journal of Remote Sensing, 46(3), 360–375. https://doi.org/10.1080/07038992.2020.1802584Paper on CWIM3:M. Mahdianpari et al., "The Third Generation of Pan-Canadian Wetland Map at 10 m Resolution Using Multisource Earth Observation Data on Cloud Computing Platform," in IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing, vol. 14, pp. 8789-8803, 2021, doi: 10.1109/JSTARS.2021.3105645.Paper on Arctic ecoregion enhancement for CWIM3A:Michael Merchant, et al., ”Leveraging google earth engine cloud computing for large-scale arctic wetland mapping,” in International Journal of Applied Earth Observation and Geoinformation, vol. 125, 2023, https://doi.org/10.1016/j.jag.2023.103589.

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

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

Post-disturbance forest recovery data for Canada's forested ecosystems, representing a total area of ~650 million ha, captures the return of forests following wildfire and harvest that occurred between 1986 and 2012. It is developed within the framework of Canada’s National Terrestrial Ecosystem Monitoring System (NTEMS). These spatially-explicit outputs represent the rate of spectral recovery: the rate at which a pixel returns to 80% of its pre-disturbance value (White et al. 2017) within the observation period (1985-2017) using the Y2R or Years-to-Recovery metric derived from Landsat times series data. Baseline rates of spectral recovery (Y2R) were defined for each of Canada's 12 forested ecozones. These baselines were then used to identify spatial clusters of recovering pixels on the landscape where Y2R were either significantly faster or slower than their ecozonal baseline. Finally, areas that were disturbed by wildfire and harvest (1986-2012), but which had not recovered by the end of the observation period (2017) are also provided. Note that these areas are still recovering, but they had not yet recovered according to our metric of spectral recovery, by the end of the time series in 2017. For an overview of the methods, the validation of the Y2R metric, and interpretation of the derived trends, see White et al. (2022) and White et al. (2017).White, J.C., Hermosilla, T., Wulder, M.A., Coops, N.C., 2022. Mapping, validating, and interpreting spatio-temporal trends in post-disturbance forest recovery. Remote Sensing of Environment, 271, 112904. https://doi.org/10.1016/j.rse.2022.112904 ( White et al. 2022)White, J.C., Wulder, M.A., Hermosilla, T., Coops, N.C., Hobart, G.W. 2017. A nationwide annual characterization of 25 years of forest disturbance and recovery for Canada using Landsat time series. Remote Sensing of Environment, 194, pp. 303-321. DOI: https://doi.org/10.1016/j.rse.2017.03.035 .( White et al. 2017)

All available bathymetry and related information for Square Lake were collected and hard copy maps digitized where necessary. The data were validated against more recent data (Shuttle Radar Topography Mission 'SRTM' imagery and Indian Remote Sensing 'IRS' imagery) and corrected where necessary. The published data set contains the lake bathymetry formatted as an Arc ascii grid. Bathymetric contours and the boundary polygon are available as shapefiles.

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)

Wildfire Year/dNBR/Mask 1985-2015Wildfire change magnitude 85-15. Spectral change magnitude for wildfires that occurred from 1985 and 2015. 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 the change event. This dataset is composed of three layers: (1) binary wildfire mask, (2) year of greatest wildfire disturbance, and (3) differenced Normalized Burn Ratio (dNBR) transformed for data storage efficiency to the range 0-200. The actual dNBR value is derived as follows: dNBR = value / 100. Higher dNBR values are related to higher burn severity. The information outcomes represent 30 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-2015 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).Geographic extent: Canada's forested ecosystems (~ 650 Mha)Time period: 1985–2011

All available bathymetry and related information for Grassy Island Lake were collected and hard copy maps digitized where necessary. The data were validated against more recent data (Shuttle Radar Topography Mission 'SRTM' imagery and Indian Remote Sensing 'IRS' imagery) and corrected where necessary. The published data set contains the lake bathymetry formatted as an Arc ascii grid. Bathymetric contours and the boundary polygon are available as shapefiles.