High resolution forest change for Canada (Binary Change/No-change) 1985-2011The forest change data included in this product is national in scope (entire forested ecosystem) and represents the first wall-to-wall characterization of wildfire and harvest in Canada at a spatial resolution commensurate with human impacts. The information outcomes represent 27 years of stand replacing change in Canada’s forests, derived from a single, consistent spatially-explicit data source, derived in a fully automated manner. This demonstrated capacity to characterize forests at a resolution that captures human impacts is key to establishing a baseline for detailed monitoring of forested ecosystems from management and science perspectives. Time series of Landsat data were used to characterize national trends in stand replacing forest disturbances caused by wildfire and harvest for the period 1985–2011 for Canada's 650 million hectare forested ecosystems (https://authors.elsevier.com/sd/article/S0034425717301360 ). Landsat data has a 30m spatial resolution, so the change information is highly detailed and is commensurate with that of human impacts. These data represent annual stand replacing forest changes. The stand replacing disturbances types labeled are wildfire and harvest, with lower confidence wildfire and harvest, also shared. The distinction and sharing of lower class membership likelihoods is to indicate to users that some change events were more difficult to allocate to a change type, but are generally found to be in the correct category. For an overview on the data, image processing, and time series change detection methods applied, as well as information on independent accuracy assessment of the data, see Hermosilla et al. (2016; http://www.tandfonline.com/doi/full/10.1080/17538947.2016.1187673).The data available is, 1. a binary change/no-change; 2. Change year; and, 3. Change type. When using this data, please cite as: White, J.C., M.A. Wulder, T. Hermosilla, N.C. Coops, and G. Hobart. (2017). A nationwide annual characterization of 25 years of forest disturbance and recovery for Canada using Landsat time series. Remote Sensing of Environment. 192: 303-321. DOI: 10.1016/j.rse.2017.03.035. https://authors.elsevier.com/sd/article/S0034425717301360Geographic extent: Canada's forested ecosystems (~ 650 Mha)Time period: 1985–2011

High resolution forest change for Canada (Change Type) 1985-2011The forest change data included in this product is national in scope (entire forested ecosystem) and represents the first wall-to-wall characterization of wildfire and harvest in Canada at a spatial resolution commensurate with human impacts. The information outcomes represent 27 years of stand replacing change in Canada’s forests, derived from a single, consistent spatially-explicit data source, derived in a fully automated manner. This demonstrated capacity to characterize forests at a resolution that captures human impacts is key to establishing a baseline for detailed monitoring of forested ecosystems from management and science perspectives. Time series of Landsat data were used to characterize national trends in stand replacing forest disturbances caused by wildfire and harvest for the period 1985–2011 for Canada's 650 million hectare forested ecosystems (https://authors.elsevier.com/sd/article/S0034425717301360 ). Landsat data has a 30m spatial resolution, so the change information is highly detailed and is commensurate with that of human impacts. These data represent annual stand replacing forest changes. The stand replacing disturbances types labeled are wildfire and harvest, with lower confidence wildfire and harvest, also shared. The distinction and sharing of lower class membership likelihoods is to indicate to users that some change events were more difficult to allocate to a change type, but are generally found to be in the correct category. For an overview on the data, image processing, and time series change detection methods applied, as well as information on independent accuracy assessment of the data, see Hermosilla et al. (2016; http://www.tandfonline.com/doi/full/10.1080/17538947.2016.1187673).The data available is, 1. a binary change/no-change; 2. Change year; and, 3. Change type. When using this data, please cite as: White, J.C., M.A. Wulder, T. Hermosilla, N.C. Coops, and G. Hobart. (2017). A nationwide annual characterization of 25 years of forest disturbance and recovery for Canada using Landsat time series. Remote Sensing of Environment. 192: 303-321. DOI: 10.1016/j.rse.2017.03.035. https://authors.elsevier.com/sd/article/S0034425717301360Geographic extent: Canada's forested ecosystems (~ 650 Mha)Time period: 1985–2011

High resolution forest change for Canada (Change Year) 1985-2011The forest change data included in this product is national in scope (entire forested ecosystem) and represents the first wall-to-wall characterization of wildfire and harvest in Canada at a spatial resolution commensurate with human impacts. The information outcomes represent 27 years of stand replacing change in Canada’s forests, derived from a single, consistent spatially-explicit data source, derived in a fully automated manner. This demonstrated capacity to characterize forests at a resolution that captures human impacts is key to establishing a baseline for detailed monitoring of forested ecosystems from management and science perspectives. Time series of Landsat data were used to characterize national trends in stand replacing forest disturbances caused by wildfire and harvest for the period 1985–2011 for Canada's 650 million hectare forested ecosystems (https://authors.elsevier.com/sd/article/S0034425717301360 ). Landsat data has a 30m spatial resolution, so the change information is highly detailed and is commensurate with that of human impacts. These data represent annual stand replacing forest changes. The stand replacing disturbances types labeled are wildfire and harvest, with lower confidence wildfire and harvest, also shared. The distinction and sharing of lower class membership likelihoods is to indicate to users that some change events were more difficult to allocate to a change type, but are generally found to be in the correct category. For an overview on the data, image processing, and time series change detection methods applied, as well as information on independent accuracy assessment of the data, see Hermosilla et al. (2016; http://www.tandfonline.com/doi/full/10.1080/17538947.2016.1187673).The data available is, 1. a binary change/no-change; 2. Change year; and, 3. Change type. When using this data, please cite as: White, J.C., M.A. Wulder, T. Hermosilla, N.C. Coops, and G. Hobart. (2017). A nationwide annual characterization of 25 years of forest disturbance and recovery for Canada using Landsat time series. Remote Sensing of Environment. 192: 303-321. DOI: 10.1016/j.rse.2017.03.035. https://authors.elsevier.com/sd/article/S0034425717301360Geographic extent: Canada's forested ecosystems (~ 650 Mha)Time period: 1985–2011

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.

Wildfire Year/dNBR/Mask 1985-2015Wildfire change magnitude 85-15. Spectral change magnitude for wildfires that occurred from 1985 and 2015. 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

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)

Deforestation in Canada is estimated with the National Deforestation Monitoring System (NDMS). Details describing this process are published here: https://cfs.nrcan.gc.ca/publications?id=36042. Deforestation is the direct human-induced conversion of forested land to non-forested land use. Canada’s National Deforestation Monitoring System (NDMS) was designed and implemented to provide information needed by Canada to meet its obligation under the United Nations Framework Convention on Climate Change (UNFCCC) to report the areas affected annually by deforestation. It also provides important information for the public, government policy makers, and scientists. To provide information about the amount of deforestation and why, where, and when it occurred in Canada, the NDMS uses deforestation mapped on a system of sample areas. The mapping is based on visual interpretation of satellite imagery supported by available ancillary information, such as high resolution imagery, forest inventory, and industrial databases, and informed by records-based information and expert knowledge. Accurate detection and mapping of deforestation events involves manual interpretation of satellite remote sensing imagery by specialized analysts. A key factor in the mapping is to distinguish deforestation from other forest cover losses that occur in Canada. The NDMS was designed to make use of all available lines of evidence and be flexible to accommodate variable resourcing levels. This system has been producing national deforestation monitoring results annually since 2006. The flexibility of the NDMS’s design makes it possible to adapt to future changes in data and resource availability, and positions the program well for sustained operational delivery into the future.

Monthly, seasonal and annual trends of total precipitation change (millimetres) based on adjusted station data (AHCCD) are available. Trends are calculated using the Theil-Sen method using the station’s full period of available data. The availability of precipitation trends will vary by station; if more than 5 consecutive years are missing data or more than 10% of the data within the time series is missing, a trend was not calculated.

Monthly, seasonal and annual trends of daily minimum, mean and maximum surface air temperature change (degrees Celsius) based on homogenized station data (AHCCD) are available. Trends are calculated using the Theil-Sen method using the station’s full period of available data. The availability of temperature trends will vary by station; if more than 5 consecutive years are missing data or more than 10% of the data within the time series is missing, a trend was not calculated.

Data include a collection of annual land cover maps derived from MODIS 250 m spatial resolution remotely sensed imagery for the period 2000 to 2011. Processing of the time series was designed to reduce the occurrence of false change between maps. The method was based on change updating as described in Pouliot et al. (2011, 2013). Change detection accounted for both abrupt changes such as forest harvesting and more gradual changes such as recurrent insect defoliation. To determine the new label for a pixel identified as change, an evidential reasoning approach was used to combine spectral and contextual information. The 2005 MODIS land cover of Canada at 250 m spatial resolution described in Latifovic et al. (2012) was used as the base map. It contains 39 land cover classes, which for time series development was considered too detailed and was reduced to 25 and 19 class versions. The 19 class version corresponds to the North America Land Change Monitoring System (NALCMS) Level 2 legend as described in Latifovic et al. (2012). Accuracy assessment of time series is difficult due to the need to assess many maps. For areas of change in the time series accuracy was found to be 70% based on the 19 class thematic legend.This time series captures the spatial distribution of dominant land cover transitions. It is intended for use in modeling, development of remote sensing products such as leaf area index or land cover based albedo retrievals, and other exploratory analysis. It is not appropriate for use in any rigorous reporting or inventory assessments due to the accuracy of the land cover classification and uncertainty as to the capture of all relevant changes for an application.