The Canadian long term satellite data record (LTDR) derived from 1-km resolution Advanced Very High Resolution Radiometer (AVHRR) data was produced by the Canada Center for Remote Sensing (CCRS). Processing included: geolocation, calibration, and compositing using Earth Observation Data Manager (Latifovic et al. 2005), cloud screening (Khlopenkov and Trishchenko, 2006), BRDF correction (Latifovic et. al., 2003), atmosphere and other corrections as described in Cihlar et. al. (2004). For temporal analysis of vegetation cross-sensor correction of Latifovic et al. (2012) is advised. Data collected by the AVHRR instrument on board the National Oceanic and Atmospheric Administration (NOAA) 9,11,14,16,17,18 and 19 satellites were used to generate Canada-wide 1-km 10-day AVHRR composites. Data are available starting in 1985. It is important to note that there are three types of AVHRR sensors: (i) AVHRR-1 flown onboard TIROS-N, NOAA-6, NOAA-8, and NOAA-10; (ii) AVHRR-2 flown onboard NOAA-7, NOAA-9, NOAA-11, NOAA-12, and NOAA-14; and (iii) AVHRR-3 currently operational onboard NOAA-15, NOAA-16, NOAA-17, NOAA-18 and NOAA-19. The AVHRR-1 has four channels, AVHRR-2 has five channels and the AVHRR-3 has six channels, although only five channels of AVHRR-3 can be operational at any one time. As such, channels 3A (1.6 m) and 3B (3.7 m) work interchangeably. The processing procedure was designed to minimize artefacts in AVHRR composite images. There are thirty six 10-day image composites per year. The following three processing levels are provided: P1) top of atmosphere reflectance and brightness temperature, P2) reflectance at surface and surface temperature and P3) reflectance at surface normalized to a common viewing geometry (BRDF normalization). The processing level P1 and P2 are provided for all 36 composites while level P3 is provided for 21 composites from April – October.

Land cover classification image for the Aspen Parkland ecoregion of Saskatchewan with a spatial resolution of 10m. The goal of this land cover classification was to distinguish native from tame grasslands. The classification was based on Sentinel-1 and Sentinel-2 imagery using machine learning analysis in the Google Earth Engine platform. The classification was conducted on imagery acquired in 2022, and the classification model was built with field data collected in 2020 - 2022. There are eight classes in total: native grassland, tame grassland, mixed/altered grassland, cropland, shrubs, trees, water, and urban area. Download: here The Prairie Landscape Inventory (PLI) aims to develop improved methods of assessing land cover and land use for conservation. Native grassland has historically been one of the hardest to map at-risk ecosystems because of the challenges in distinguishing native grassland from tame grassland land cover using remotely sensed imagery. This classification distinguishes native grassland from tame grassland and will provide valuable information for habitat suitability for native grassland species, identifying high biodiversity potential and invasion risk potential. The classification map has eight (8) classes:  1. Cropland This class represents all cultivated areas with crop commodities, including corn, pulse, soybeans, canola, grains, and summer-fallow. 2. Native grassland This class represents the native grassland areas that are composed of at least 75% native grass, sedge and forb species, such as the needle grasses and wheatgrasses along with June grass and blue grama grass. Unbroken grassland that is invaded by species like Kentucky bluegrass, crested wheatgrass or smooth brome, such that native cover is less than 75%, is not considered native for the purpose of this project.  3. Mixed/altered grassland This class represents a grassland with a mix of less than 75% native grass, sedge and forb species or less than 75% tame species. These are grassland areas that do not fit into either of the native or tame grassland definitions.  4. Tame grassland This class represents the tame grassland areas that are composed of at least 75% seeded or planted species with introduced grasses and forb species such as crested wheatgrass, smooth brome, Kentucky bluegrass, alfalfa, and sweet clover. 5. Water This class represents permanent water locations such as lakes and rivers. 6. Shrubs This class represents the sites dominated by woody vegetation of relatively low height (generally +/-2 meters) with shrub canopy typically >20% of total vegetation cover. 7. Trees This class represents the coniferous/deciduous trees, mixed-wood area, and other trees >2 meters height with tree canopy typically >20% of total vegetation cover. 9. Urban area This class represents both urban municipalities and buffered roads. Urban municipalities was used to mask the urban/developed area class of the Annual Crop Inventory 2021 (Agriculture Agri-Food Canada).   Colour Classes: Value Label Red Green Blue 1 Cropland 255 255 190 2 Native grassland 168 168 0 3 Mixed/altered grassland 199 215 158 4 Tame grassland 245 202 122 5 Water 190 232 255 6 Shrubs 205 102 153 7 Trees 66 128 53 9 Urban area 128 128 128 Accuracy metrics This model has an overall accuracy of 73 per cent. The table below summarizes the user’s accuracy, producer’s accuracy, and F1-score of the model on the validation dataset. Class User’s accuracy (%) Producer’s accuracy (%) F1-score Cropland 91.2 94.5 0.93 Native grassland 74.8 73.1 0.74 Mixed grassland 44.7 44.1 0.44 Tame grassland 67.9 72.8 0.70 Water 94.8 91.3 0.93 Shrubs 61.2 51.1 0.56 Trees 89.7 94.6 0.92

This data publication contains two collections of raster maps of forest attributes across Canada, the first collection for year 2001, and the second for year 2011. The 2001 collection is actually an improved version of an earlier set of maps produced also for year 2001 (Beaudoin et al 2014, DOI: https://doi.org/10.1139/cjfr-2013-0401) that is itself available through the web site “http://nfi-nfis.org”. Each collection contains 93 maps of forest attributes: four land cover classes, 11 continuous stand-level structure variables such as age, volume, biomass and height, and 78 continuous values of percent composition for tree species or genus. The mapping was done at a spatial resolution of 250m along the MODIS grid. Briefly the method uses forest polygon information from the first version of photoplots database from Canada’s National Forest Inventory as reference data, and the non-parametric k-nearest neighbors procedure (kNN) to create the raster maps of forest attributes. The approach uses a set of 20 predictive variables that include MODIS spectral reflectance data, as well as topographic and climate data. Estimates are carried out on target pixels across all Canada treed landmass that are stratified as either forest or non-forest with 25% forest cover used as a threshold. Forest cover information was extracted from the global forest cover product of Hansen et al (2013) (DOI: https://doi.org/10.1126/science.1244693). The mapping methodology and resultant datasets were intended to address the discontinuities across provincial borders created by their large differences in forest inventory standards. Analysis of residuals has failed to reveal residual discontinuities across provincial boundaries in the current raster dataset, meaning that our goal of providing discontinuity-free maps has been reached. The dataset was developed specifically to address strategic issues related to phenomena that span multiple provinces such as fire risk, insect spread and drought. In addition, the use of the kNN approach results in the maintenance of a realistic covariance structure among the different variable maps, an important property when the data are extracted to be used in models of ecosystem processes. For example, within each pixel, the composition values of all tree species add to 100%. * Details on the product development and validation can be found in the following publication: Beaudoin, A., Bernier, P.Y., Villemaire, P., Guindon, L., Guo, X.-J. 2017. Tracking forest attributes across Canada between 2001 and 2011 using a kNN mapping approach applied to MODIS imagery, Canadian Journal of Forest Research 48: 85–93. DOI: https://doi.org/10.1139/cjfr-2017-0184 * Please cite this dataset as: Beaudoin A, Bernier PY, Villemaire P, Guindon L, Guo XJ. 2017. Species composition, forest properties and land cover types across Canada’s forests at 250m resolution for 2001 and 2011. Natural Resources Canada, Canadian Forest Service, Laurentian Forestry Centre, Quebec, Canada. https://doi.org/10.23687/ec9e2659-1c29-4ddb-87a2-6aced147a990 * This dataset contains these NFI forest attributes: ## LAND COVER : landbase vegetated, landbase non-vegetated, landcover treed, landcover non-treed ## TREE STRUCTURE : total above ground biomass, tree branches biomass, tree foliage biomass, stem bark biomass, stem wood biomass, total dead trees biomass, stand age, crown closure, tree stand heigth, merchantable volume, total volume ## TREE SPECIES : abies amabilis (amabilis fir), abies balsamea (balsam fir), abies lasiocarpa (subalpine fir), abies spp. (unidentified fir), acer macrophyllum (bigleaf maple), acer negundo (manitoba maple, box-elder), acer pensylvanicum (striped maple), acer rubrum (red maple), acer saccharinum (silver maple), acer saccharum (sugar maple), acer spicatum (mountain maple), acer spp. (unidentified maple), alnus rubra (red alder), alnus spp. (unidentified alder), arbutus menziesii (arbutus), betula alleghaniensis (yellow birch), betula papyrifera (white birch), betula populifolia (gray birch), betula spp. (unidentified birch), carpinus caroliniana (blue-beech), carya cordiformis (bitternut hickory), chamaecyparis nootkatensis (yellow-cedar), fagus grandifolia (american beech), fraxinus americana (white ash), fraxinus nigra (black ash), fraxinus pennsylvanica (red ash), juglans cinerea (butternut), juglans nigra (black walnut), juniperus virginiana (eastern redcedar), larix laricina (tamarack), larix lyallii (subalpine larch), larix occidentalis (western larch), larix spp. (unidentified larch), malus spp. (unidentified apple), ostrya virginiana (ironwood, hop-hornbeam), picea abies (norway spruce), picea engelmannii (engelmann spruce), picea glauca (white spruce), picea mariana (black spruce), picea rubens (red spruce), picea sitchensis (sitka spruce), picea spp. (unidentified spruce), pinus albicaulis (whitebark pine), pinus banksiana (jack pine), pinus contorta (lodgepole pine), pinus monticola (western white pine), pinus ponderosa (ponderosa pine), pinus resinosa (red pine), pinus spp. (unidentified pine), pinus strobus (eastern white pine), pinus sylvestris (scots pine), populus balsamifera (balsam poplar), populus grandidentata (largetooth aspen), populus spp. (unidentified poplar), populus tremuloides (trembling aspen), populus trichocarpa (black cottonwood), prunus pensylvanica (pin cherry), prunus serotina (black cherry), pseudotsuga menziesii (douglas-fir), quercus alba (white oak), quercus macrocarpa (bur oak), quercus rubra (red oak), quercus spp. (unidentified oak), salix spp. (unidentified willow), sorbus americana (american mountain-ash), thuja occidentalis (eastern white-cedar), thuja plicata (western redcedar), tilia americana (basswood), tsuga canadensis (eastern hemlock), tsuga heterophylla (western hemlock), tsuga mertensiana (mountain hemlock), tsuga spp. (unidentified hemlock), ulmus americana (white elm), unidentified needleaf, unidentified broadleaf, broadleaf species, needleaf species, unknown species

Landsat Thematic Mapper (TM) and Enhanced Thematic Mapper (ETM+) sensors were used to generate the circa 2010 Mosaic of Canada at 30 m spatial resolution. All scenes were processed to Standard Terrain Correction Level 1T by the United States Geological Survey (USGS). Further processing performed by the Canada Centre for Remote Sensing included conversion of sensor measurements to top of atmosphere reflectance, cloud and cloud shadow detection, re-projection, selection of best measurements, mosaic generation ,noise removal and quality control. To provide a clear sky measurement for each location in Canada, data from the years 2009, 2010, and 2011 were used, but 2010 was preferentially selected. Bands 3 (0.63-0.69 µm), 4 (0.76-0.90 µm), 5 (1.55-1.75 µm), and 7 (2.08-2.35 µm) are provided in this version as significant atmosphere effects strongly limit the quality of the blue (0.45-0.52 µm) and green (0.52-0.60 µm) bands. Multi-criteria compositing was used for the selection of the most representative pixel. For ETM+ onboard Landsat 7 a scan line malfunction caused missing lines of data in all scenes collected after May 2003. Atmosphere and target variability between scenes cause these lines to have significant radiometric differences in some cases. A Fourier transformation approach was applied to correct this occurrence. This mosaic was developed for land cover and biophysical mapping applications across Canada. Other applications of these data are also possible, but should consider the temporal and spectral limitations of the product. Research to enhance the spatial, spectral and temporal aspects are in development for future versions of moderate resolution products from historical Landsat sensors, Landsat 8, and Sentinel 2 data.

Note: To visualize the data in the viewer, zoom into the area of interest. The National Air Photo Library (NAPL) of Natural Resources Canada archives over 6 million aerial photographs covering all of Canada, some of which date back to the 1920s. This collection includes Time Series of aerial orthophoto mosaics over a selection of major cities or targeted areas that allow the observation of various changes that occur over time in those selected regions. These mosaics are disseminated through the Data Cube Platform implemented by NRCan using geospatial big data management technologies. These technologies enable the rapid and efficient visualization of high-resolution geospatial data and allow for the rapid generation of dynamically derived products. The data is available as Cloud Optimized GeoTIFF (COG) for direct access and as Web Map Services (WMS) or Web Coverage Services (WCS) with a temporal dimension for consumption in Web or GIS applications. The NAPL mosaics are made from the best spatial resolution available for each time period, which means that the orthophotos composing a NAPL Time Series are not necessarily coregistrated. For this dataset, the spatial resolutions are: 25 cm for the year 1950, 75 cm for the year 1960, 50 cm for the year 1964, 75 cm for the year 1973, 75 cm for the year 1994 and 50 cm for the year 2001. The NAPL indexes and stores federal aerial photography for Canada, and maintains a comprehensive historical archive and public reference centre. The Earth Observation Data Management System (EODMS) online application allows clients to search and retrieve metadata for over 3 million out of 6 million air photos. The EODMS online application enables public and government users to search and order raw Government of Canada Earth Observation images and archived products managed by NRCan such as aerial photos and satellite imagery. To access air photos, you can visit the EODMS web site: https://eodms-sgdot.nrcan-rncan.gc.ca/index-en.html

The dataset includes two data products derived from the Visible Infrared Imaging Radiometer Suite (VIIRS) imager operated by the US National Oceanic and Atmospheric Administration (NOAA) onboard Suomi National Polar-Orbiting Partnership (SNPP) satellite: 1) Normalized Difference Vegetation Index (NDVI) 2) Snow Mask (Snow) with supplementary information about data quality and scene identification Each product, NDVI and Snow, has been derived at two spatial resolutions: 1) I-band resolution for 250-m spatial grid (VIIRS image bands I1 and I2) 2) M-band resolution for 500-m spatial grid (VIIRS moderate resolution bands M5 and M7) Datasets are produced with a daily temporal frequency, i.e. one file per day. The study area with the size of 5,700 km × 4,800 km covers Canada and neighboring regions (Trishchenko, 2019). The VIIRS time series are produced from VIIRS /SNPP imagery at CCRS from January 1, 2017.

Leaf area index (LAI) quantified the density of vegetation irrespective of land cover. LAI quantifies the total foliage surface area per groud surface area. LAI has been identified by the Global Climate Observing System as an essential climate variable required for ecosystem,weather and climate modelling and monitoring. This product consists of annual maps of the maximum LAI during a grownig season (June-July-August) at 100m resolution covering Canada's land mass.

Index Grid for NTS 1:250,000 scale maps

Note: To visualize the data in the viewer, zoom into the area of interest. The National Air Photo Library (NAPL) of Natural Resources Canada archives over 6 million aerial photographs covering all of Canada, some of which date back to the 1920s. This collection includes Time Series of aerial orthophoto mosaics over a selection of major cities or targeted areas that allow the observation of various changes that occur over time in those selected regions. These mosaics are disseminated through the Data Cube Platform implemented by NRCan using geospatial big data management technologies. These technologies enable the rapid and efficient visualization of high-resolution geospatial data and allow for the rapid generation of dynamically derived products. The data is available as Cloud Optimized GeoTIFF (COG) for direct access and as Web Map Services (WMS) or Web Coverage Services (WCS) with a temporal dimension for consumption in Web or GIS applications. The NAPL mosaics are made from the best spatial resolution available for each time period, which means that the orthophotos composing a NAPL Time Series are not necessarily coregistrated. For this dataset, the spatial resolutions are: 100 cm for the year 1932 and 50 cm for the year 1950. The NAPL indexes and stores federal aerial photography for Canada, and maintains a comprehensive historical archive and public reference centre. The Earth Observation Data Management System (EODMS) online application allows clients to search and retrieve metadata for over 3 million out of 6 million air photos. The EODMS online application enables public and government users to search and order raw Government of Canada Earth Observation images and archived products managed by NRCan such as aerial photos and satellite imagery. To access air photos, you can visit the EODMS web site: https://eodms-sgdot.nrcan-rncan.gc.ca/index-en.html

The 2005 AAFC Land Use is a culmination and curated metaanalysis of several high-quality spatial datasets produced between 1990 and 2021 using a variety of methods by teams of researchers as techniques and capabilities have evolved. The information from the input datasets was consolidated and embedded within each 30m x 30m pixel to create consolidated pixel histories, resulting in thousands of unique combinations of evidence ready for careful consideration. Informed by many sources of high-quality evidence and visual observation of imagery in Google Earth, we apply an incremental strategy to develop a coherent best current understanding of what has happened in each pixel through the time series.