This data publication contains a set of files in which areas affected by fire or by harvest from 1984 to 2015 are identified at the level of individual 30m pixels on the Landsat grid. Details of the product development can be found in Guindon et al (2018). The change detection is based on reflectance-corrected yearly summer (July and August) Landsat mosaics from 1984 to 2015 created from individual scenes developed from USGS reflectance products (Masek et al, 2006; Vermote et al, 2006). Briefly, the change detection method uses a six-year temporal signature centered on the disturbance year to identify fire, harvest and no change. The signatures were derived from visually-interpreted disturbance or no-change polygons that were used to fit a decision tree model. The method detects about 91% of the areas harvested and 85% of the areas burned across Canada’s forests over the study period, but overestimates areas disturbed in the two initial and mostly in the two final years of the 1985 to 2015 time series. This is caused by the absence of appropriate pre-disturbance and post-disturbance data for the model-based detection and attribution. Disturbance coverage in those four years should therefore be used with caution. As in Guindon et al (2014), the method was designed to minimize commission errors and has a disturbance class attribution success rate of about 98%. The attribution success rate of disturbance year for fire is of about 69% for the exact year and of about 99% when attribution to the following year is also considered as a success. This common one-year lag is mostly due to the use of mid-summer Landsat mosaics for the analysis that will cause spring and fall events of the same year to be attributed to successive years. For example, a fire that occurred in the fall of 2004 (after July and August), will be detected and attributed to 2005, while for a fire that occurred in the spring of 2004 will be detected and attributed to 2004. The presence of clouds and shadows or image availability causes 10% of missing data annually and therefore can too delay the capture of events. The data provides uniform spatial and temporal information on fire and harvest across all provinces and territories of Canada and is intended for strategic-level analysis. Since no attention was given to other minor disturbances such as mining, road or flooding, the product should not be used for their identification. Finally, calibration datasets were developed for only three major forest pests (mountain pine beetle, eastern spruce budworm and forest tent caterpillar), but were folded within the “no-change” class in order to minimize commission errors for fire and harvest . Less common pests for which validation datasets are hard to develop were not considered and therefore could in rare circumstances generate false fire events. Considering that area having two (3.3%) to three disturbances (less than 1%) events are not common, only the most recent disturbance is provided, overlapping older disturbances in these rare case. ## Please cite this dataset as: Guindon, L., P. Villemaire, R. St-Amant, P.Y. Bernier, A. Beaudoin, F. Caron, M. Bonucelli and H. Dorion. 2017. Canada Landsat Disturbance (CanLaD): a Canada-wide Landsat-based 30-m resolution product of fire and harvest detection and attribution since 1984. https://doi.org/10.23687/add1346b-f632-4eb9-a83d-a662b38655ad ## Scientific article citation: The creation, validation and limitations of the CanLaD product are described in the Supplementary Material file associated with the following article: Guindon, L.; Bernier, P.Y.; Gauthier, S.; Stinson, G.; Villemaire, P.; Beaudoin, A. 2018. Missing forest cover gains in boreal forests explained. Ecosphere, 9 (1) Article e02094. doi:10.1002/ecs2.2094. ## Cited references: 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. 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.

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)

In 2013, the Earth Observation Team of the Science and Technology Branch (STB) at Agriculture and Agri-Food Canada (AAFC) repeated the process of generating annual crop inventory digital maps using satellite imagery to for all of Canada, in support of a national crop inventory. A Decision Tree (DT) based methodology was applied using optical (Landsat-8) and radar (RADARSAT-2) based satellite images, and having a final spatial resolution of 30m. In conjunction with satellite acquisitions, ground-truth information was provided by provincial crop insurance companies and point observations from the BC Ministry of Agriculture and our regional AAFC colleagues.

In 2015, the Earth Observation Team of the Science and Technology Branch (STB) at Agriculture and Agri-Food Canada (AAFC) repeated the process of generating annual crop inventory digital maps using satellite imagery to for all of Canada, in support of a national crop inventory. A Decision Tree (DT) based methodology was applied using optical (Landsat-8) and radar (RADARSAT-2) based satellite images, and having a final spatial resolution of 30m. In conjunction with satellite acquisitions, ground-truth information was provided by provincial crop insurance companies and point observations from the BC Ministry of Agriculture and our regional AAFC colleagues.

In 2014, the Earth Observation Team of the Science and Technology Branch (STB) at Agriculture and Agri-Food Canada (AAFC) repeated the process of generating annual crop inventory digital maps using satellite imagery to for all of Canada, in support of a national crop inventory. A Decision Tree (DT) based methodology was applied using optical (Landsat-8) and radar (RADARSAT-2) based satellite images, and having a final spatial resolution of 30m. In conjunction with satellite acquisitions, ground-truth information was provided by provincial crop insurance companies and point observations from the BC Ministry of Agriculture and our regional AAFC colleagues.

In 2017, the Earth Observation Team of the Science and Technology Branch (STB) at Agriculture and Agri-Food Canada (AAFC) repeated the process of generating annual crop inventory digital maps using satellite imagery to for all of Canada, in support of a national crop inventory. A Decision Tree (DT) based methodology was applied using optical (Landsat-8, Sentinel-2, Gaofen-1) and radar (RADARSAT-2) based satellite images, and having a final spatial resolution of 30m. In conjunction with satellite acquisitions, ground-truth information was provided by: provincial crop insurance companies in Alberta, Saskatchewan, Manitoba, & Quebec; point observations from the BC Ministry of Agriculture, & the Ontario Ministry of Agriculture, Food and Rural Affairs; and data collection supported by our regional AAFC Research and Development Centres in St. John’s, Kentville, Charlottetown, Fredericton, Guelph, and Summerland

In 2016, the Earth Observation Team of the Science and Technology Branch (STB) at Agriculture and Agri-Food Canada (AAFC) repeated the process of generating annual crop inventory digital maps using satellite imagery to for all of Canada, in support of a national crop inventory. A Decision Tree (DT) based methodology was applied using optical (Landsat-8, Sentinel-2, Gaofen-1) and radar (RADARSAT-2) based satellite images, and having a final spatial resolution of 30m. In conjunction with satellite acquisitions, ground-truth information was provided by: provincial crop insurance companies in Alberta, Saskatchewan, Manitoba, & Quebec; point observations from the BC Ministry of Agriculture, & the Ontario Ministry of Agriculture, Food and Rural Affairs; and data collection supported by our regional AAFC Research and Development Centres in St. John’s, Kentville, Charlottetown, Fredericton, Guelph, and Summerland.

In 2010 the Earth Observation Team of the Science and Technology Branch (STB) at Agriculture and Agri-Food Canada (AAFC) continued the process of generating annual crop inventory digital maps using satellite imagery. Focusing on the Prairie Provinces, a Decision Tree (DT) based methodology was applied using both optical (AWiFS, Landsat-5, DMC) and radar (RADARSAT-2) based satellite imagery, and having a final spatial resolution of 56m. Methods were also developed to enhance the optical classification with RADARSAT-2 imagery, addressing issues associated with cloud cover. In conjunction with satellite acquisitions, ground-truth information was provided by provincial crop insurance companies and point observations from our regional AAFC colleagues. The overall process for Crop Inventory Map includes: satellite data acquisition; field data acquisition for classification training and accuracy assessment; and, operational implementation of the classification methodology.

In 2009 the Earth Observation Team of the Science and Technology Branch (STB) at Agriculture and Agri-Food Canada (AAFC) began the process of generating annual crop inventory digital maps using satellite imagery. Focusing on the Prairie Provinces, a Decision Tree (DT) based methodology was applied using both optical (AWiFS, Landsat-5) and radar (RADARSAT-2) based satellite imagery, and having a final spatial resolution of 56m. Methods were also developed to enhance the optical classification with RADARSAT-2 imagery, addressing issues associated with cloud cover. In conjunction with satellite acquisitions, ground-truth information was provided by provincial crop insurance companies and point observations from our regional AAFC colleagues. The overall process for Crop Inventory Map includes: satellite data acquisition; field data acquisition for classification training and accuracy assessment; and, operational implementation of the classification methodology.The initial methodology was developed in partnership with AAFC Research Branch, and supported in part by the Canadian Space Agency. The long-term objective of this endeavour is to expand from the Prairies and produce an annual crop inventory of the entire agricultural extent of Canada.

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