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.

This collection is a legacy product that is no longer supported. It may not meet current government standards. Land Cover information is the result of vectorization of raster thematic data originating from classified Landsat 5 and Landsat 7 ortho-images, for agricultural and forest areas of Canada, and for Northern Territories. The forest cover was produced by the Earth Observation for Sustainable Development (EOSD) project, an initiative of the Canadian Forest Service (CFS) with the collaboration of the Canadian Space Agency (CSA) and in partnership with the provincial and territorial governments. The agricultural coverage is produced by the National Land and Water Information Service (NLWIS) of Agriculture and Agri-Food Canada (AAFC). Northern Territories land cover was realized by the Canadian Centre of Remote Sensing (CCRS). Land Cover data are classified according to a harmonized legend build from the partner's legends. This legend is principally based on the legend described in following publication: EOSD publication: EOSD Land Cover Classification Legend Report, on which CFS and AAFC collaborated. Some classes related to Northern environments where added in order to meet the interpretation of the Northern land cover classification experts. Initially, Land Cover vector data are closest as possible to the source (original raster data). Slight differences can occur because the raster data goes through a data portrayal before being vectorized in order to enhance visual representation such as minimum size, smoothness of polygons and geometry.

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.

Each pixel value corresponds to the mean historical “Best-quality” Max-NDVI value for a given week, as calculated from the previous 20 years in the MODIS historical record (i.e. does not include data from the current year). These data are also often referred to as “weekly baselines” or “weekly normals”.

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

Each pixel value corresponds to the quality control, cloud cover and snow fraction value for each pixel in the Best-Quality Max-NDVI product.

Each pixel value corresponds to the day-of-week (1-7) from which the Weekly Best-Quality NDVI retrieval is obtained (1 = Monday, 7 = Sunday).

Each pixel value corresponds to the difference (anomaly) between the mean “Best-Quality” Max-NDVI of the week specified (e.g. Week 18, 2000-2014) and the “Best-Quality” Max-NDVI of the same week in a specific year (e.g. Week 18, 2015). Max-NDVI anomalies < 0 indicate where weekly Max-NDVI is lower than normal. Anomalies > 0 indicate where weekly Max-NDVI is higher than normal. Anomalies close to 0 indicate where weekly Max-NDVI is similar to normal.

Land cover classification image for the Cypress Upland 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 2023, and the classification model was built with field data collected in 2023. There are seven classes in total: native grassland, tame 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 seven (7) classes. The mixed grassland class included in the PLI land cover classification for other prairie ecoregions was not modelled in the Cypress Upland. 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.  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 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 92 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 96 96 0.96 Native grassland 90 93 0.92 Tame grassland 93 71 0.82 Water 100 100 1.00 Shrubs 77 88 0.83 Trees 96 996 0.96

Each pixel value corresponds to the best quality maximum NDVI recorded within that pixel over the week specified. Poor quality pixel observations are removed from this product. Observations whose quality is degraded by snow cover, shadow, cloud, aerosols, and/or low sensor zenith angles are removed (and are assigned a value of “missing data”). In addition, negative Max-NDVI values, occurring where R reflectance > NIR reflectance, are considered non-vegetated and assigned a value of 0. This results in a Max-NDVI product that should (mostly) contain vegetation-covered pixels. Max-NDVI values are considered high quality and span a biomass gradient ranging from 0 (no/low biomass) to 1 (high biomass).