This data product aimed to extend the existing pre-1985 disturbance history record by mapping wildfire, harvest, and insect outbreaks in Canadian forests between 1965 and 1984. Our geospatial data processing methodology relied on multi-layer perceptrons (MLP) trained on spectral recovery signatures to map and age these disturbances. Specific years were not assigned to insect outbreaks due to the lack of dependable training and validation data. In order to provide a more accurate data product that is compatible with existing datasets (e.g. provincial forest inventories), we used these reliable, but incomplete datasets to correct our predictions of disturbance type and year whenever they were available. Coupled with the updated Canada Landsat Disturbance (CanLaD) data product (Guindon et al. 2017), we are thus able to obtain a pan-Canadian 30m resolution disturbance history record from 1965 until 2023. The full description of the methodology and the exhaustive validation analyses are described in detail in Correia et al. (2024). The following limitations should be taken into account when using this dataset: • It is recommended to group disturbance age predictions into age classes, as this should reduce the noise present in the disturbance age estimation models. • Fire-harvest misclassification seems to be particularly common in transition zones like the southern edge of the Boreal Shield, where fires and harvest are both relatively common. • There seems to be an overestimation of 1965 fires due to a misclassification of burnt areas older than 1965 in northern, less productive areas as belonging to the beginning of our time series. • We likely detected mostly high-severity burnt areas that depict complete mortality, since the faster recovery of low-severity burns makes them more challenging to detect. • Insect outbreak detections were mostly associated with the historic eastern spruce budworm outbreak of the 1970s. Even though pixel-level insect disturbance year was not predicted, realistic estimates can be obtained by cross-checking our data product with historic reports. The following raster layers are available: • canlad_1965_1984_disturbanceType: Estimated disturbance type o 2 = Fire o 3 = Harvest o 4 = Insect • canlad_1965_1984_disturbanceYear: Estimated disturbance year o Numeric value from 1965 to 1984 • canlad_1965_1984_correctionMask: Raster indicating which predictions have been corrected with external datasets o 0 = Unconfirmed disturbance o 1 = Confirmed fire o 2 = Confirmed harvest Please cite this data product as: Correia, D. L. P., L. Guindon, and M. A. Parisien. 2024. Canada-wide Landsat-based 30-m resolution product of disturbance detection prior to 1984. https://doi.org/10.23687/660b7c6a-cdec-4c02-90c7-d63e91825c42 References: Correia, D. L. P., L. Guindon, and M. A. Parisien. 2024. Extending Canadian forest disturbance history maps prior to 1985. Ecosphere [in press]. 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

The greatest 10-day precipitation total expected over the next 2 weeks (p10d). Week 1 and week 2 forecasted index is available daily from September 1 to August 31. Week 3 and week 4 forecasted index is available weekly (Thursday) from September 1 to August 31. Units: cm/10 days Precipitation (moisture availability) establishes the economic yield potential and product quality of field crops. Both dry and wet precipitation extremes have the ability to inhibit proper crop growth. The greatest daily precipitation index covers the risk of excessive precipitation in the short term, while the other indices pertain to longer term moisture availability. Agriculture and Agri-Food Canada (AAFC) and Environment and Climate Change Canada (ECCC) have together developed a suite of extreme agrometeorological indices based on four main categories of weather factors: temperature, precipitation, heat, and wind. The extreme weather indices are intended as short-term prediction tools and generated using ECCC’s medium range forecasts to create a weekly index product on a daily and weekly basis.

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)

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.

Probability of total precipitation above 25mm over the forecast period (pweek25_prob). Week 1 and week 2 forecasted probability is available daily from September 1 to August 31. Week 3 and week 4 forecasted probability is available weekly (Thursday) from September 1 to August 31. Precipitation (moisture availability) establishes the economic yield potential and product quality of field crops. Both dry and wet precipitation extremes have the ability to inhibit proper crop growth. The greatest daily precipitation index covers the risk of excessive precipitation in the short term, while the other indices pertain to longer term moisture availability. Agriculture and Agri-Food Canada (AAFC) and Environment and Climate Change Canada (ECCC) have together developed a suite of extreme agrometeorological indices based on four main categories of weather factors: temperature, precipitation, heat, and wind. The extreme weather indices are intended as short-term prediction tools and generated using ECCC’s medium range forecasts to create a weekly index product on a daily and weekly basis.

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

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