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.

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

Annual mapping of national level forest harvesting for Canada detected inclusive of 1985 to 2015 from Landsat satellite imagery.This dataset is composed of two layers: (1) binary harvest mask, and (2) year of harvest disturbance detection. The information outcomes represent 31 years of harvesting activity in Canada’s forests, derived from a single, consistent, spatially-explicit data source in an automated manner. Time series of Landsat data with 30-m spatial resolution were used to characterize national trends in stand replacing forest disturbances, including those attributed to harvest for the period 1985–2015 for Canada's 650 million hectare forested ecosystems (Hermosilla et al. 2016). See references below for an overview regarding the data, image processing, and time-series change detection methods applied, as well as information on independent accuracy assessment of the data. 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) For additional resources on the data used and methods applied, please see: 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., Wulder, M.A., White, J.C., Coops, N.C., Hobart, G.W., 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.( Hermosilla et al. 2017)

Polygons containing the date of capture of the Landsat images used to create the first version of the Baseline Thematic Mapping v1 (BTM1). This spatial view is only meaningful in conjunction with the satellite images or the BTM data derived from the satellite images. The images were captured from 1990 to 1997