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

Since 1988, the governments of Canada and Quebec have been working together to conserve, restore, protect and develop the St. Lawrence River under the St. Lawrence Action Plan (SLAP). One of the projects identified under the theme of biodiversity conservation is the development of an integrated plan for the conservation of the natural environments and biodiversity of the St. Lawrence River. The identification of priority sites for conservation has been the first step of this planning exercise. Conservation planning of natural environments requires a reliable, accurate and up-to-date image of the spatial distribution of ecosystems in the study area. In order to produce an Atlas of Priority Sites for Conservation in the St. Lawrence Lowlands, an updated cartography of the land cover of this vast territory was undertaken. This project required obtaining reliable information on the natural environments of the St. Lawrence Lowlands. Although several land cover mapping projects have been conducted for specific types of habitats, it was particularly important to obtain a homogeneous product that would cover the entire territory and that would provide the most detailed information on its various thematic components: agricultural, aquatic, human-modified and forest environments, wetlands as well as old fields and bare ground. The methodology used to produce the land cover mapping of the St. Lawrence Lowlands thus relied mainly on combining and enhancing the best existing products for each theme. This project was made in collaboration with MDDELCC as part of the St. Lawrence Action Plan (SLAP).

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

The 0.34 cm resolution orthomosaic was created from unmanned aerial vehicle (UAV) imagery acquired from a single day survey, July 28th 2016, in Cambridge Bay, Nunavut. Five control points taken from a Global Differential Positioning System were positioned in the corners and the center of the vegetation survey. The orthomosaic covering 525m2 was produced by Canada Centre for Remote Sensing /Canada Centre for Mapping and Earth Observation. The UAV survey was completed in collaboration with the Canadian High Arctic Research Station (CHARS) for northern vegetation monitoring research. For more information, refer to our current Arctic vegetation research: Fraser et al; "UAV photogrammetry for mapping vegetation in the low-Arctic" Arctic Science, 2016, 2(3): 79-102. http://www.nrcresearchpress.com/doi/abs/10.1139/AS-2016-0008.

The MODIS surface albedo dataset was produced by the Canada Center for Remote Sensing (CCRS), Natural Resources Canada. The dataset represents the solar shortwave broadband surface albedo and it is at a 10-day interval covering the entire Canadian landmass as well as northern USA, Alaska, and the Greenland. The dataset was derived from the Moderate Resolution Imaging Spectroradiometer (MODIS) onboard the TERRA satellite which provides a global coverage every 1-2 days in 36 spectral bands ranging from visible to infrared and to thermal wavelengths between 405 and 14,385 nm, and was available since 2000. For the estimation of surface albedo, the first seven spectral bands of B1 to B7 ranging from 459 nm to 2155 nm were used. B1 and B2 have a 250 meter resolution and B3 to B7 have a 500 meter resolution. A downscaling method using a regression and normalization scheme was employed to downscale the bands B3 to B7 to 250 meter resolution while preserving radiometric properties of the original data. To obtain clear-sky observations from MODIS, composite images for a 10 day period were generated by using a series of advanced algorithms (Luo et al., 2008). The 10-day composites of B1-B7 reflectance were then used to retrieve spatially continuous spectral albedo by using a combined land/snow BRDF (Bi-directional Reflectance Distribution Function) model. In that method, the modified RossThick-LiSparse BRDF model (Maignan et al., 2004) for land and Kokhanovsky and Zege’s model (2004) for snow are linearly combined for mixed surface conditions. They are weighted by snow fraction (0.0 ~ 1.0). The seven spectral albedo were then converted into the shortwave broadband surface albedo using the empirical MODIS polynomial conversion equation of Liang et al. (1999). The data product is in LCC (Lambert Conformal Conic) projection with a 250m pixel resolution. There are 36 albedo images per year. A dataset representing the pixel state (e.g. cloud/shadow, snow/ice, water, land, et al.) was also generated for each 10-day corresponding to the surface albedo product. References: Kokhanovsky, A. A. and Zege, E. P., 2004, Scattering Optics of Snow, Applied Optics, 43, 1589-1602, doi:10.1364/AO.43.001589, 20. Liang, S., Strahler, A.H., Walthall, C., 1999. Retrieval of land surface albedo from satellite observations: a simulation study. J. Appl. Meteorol. 38, 712–725. Luo, Y., Trishchenko, A.P., Khlopenkov, K.V., 2008. Developing clear-sky, cloud and cloud shadow mask for producing clear-sky composites at 250-meter spatial resolution for the seven MODIS land bands over Canada and North America. Remote Sens. Environ. 112, 4167–4185. Maignan, F., F.M. Bréon and R. Lacaze, 2004, Bidirectional reflectance of Earth targets : evaluation of analytical models using a large set of spaceborne measurements with emphasis with the hot spot, Remote Sens. Environ., 90, 210-220.

This publication contains vector data (shapefile) of the post-harvest forest residues in Canada for the bioenergy/bioproducts sector in oven-dry tonnes per year (ODT/yr) over the next 20 years. The maps were produced using different remote sensing products. We used forest attribute data at 250 m MODIS for the years 2001 and 2011 (Beaudoin et al. 2014 and 2018) combined with forest cover changes for the years 1985 to 2015 contained in the CanLaD dataset at 30 m Landsat(Guindon et al. 2017 and 2018). Results of available biomass (in the form of harvest residues) were reported at the 10 km x 10 km scale, while the map of mature forests in Canada was prepared at the forest management unit (FMU) level. Briefly, our methodology consisted of three steps: 1- create a map of mature forests for the year 2011, based on 2001-2010 average cut volumes within FMUs; 2- develop an annual cut rate from the area harvested within FMUs from 1985 to 2015 and; 3- define the amount of biomass in the form of forest residues available for the bioenergy sector. The biomass of branches and leaves of forest attribute data was used as a proxy to define the biomass of forest residues available. Nationally, the average biomass of forest residues available after harvest is 26 ± 16 ODT/ha, while the total annual availability for all managed forests in Canada was 21 million ODT/yr. A scientific article gives additional details on the methodology: Barrette J, Paré D, Manka F, Guindon L, Bernier P, Titus B. 2018. Forecasting the spatial distribution of logging residues across the Canadian managed forest. Can. J. For. Res. 48: http://www.nrcresearchpress.com/doi/10.1139/cjfr-2018-0080 Reference for this dataset: Barrette J, Paré D, Manka F, Guindon L, Bernier P, Titus B. 2018. Maps forecasting the availability of logging residues in Canada. Natural Resources Canada, Canadian Forest Service, Laurentian Forestry Centre, Quebec, Canada. https://doi.org/10.23687/5072c495-240c-42a3-ad55-c942ab37c32a

Landcover dataset created for the northern part of Saskatchewan based on a combination of Landsat 5 Thematic Mapper (TM) and Landsat 7 Enhanced Thematic Mapper (ETM+) data representing circa 2000 conditions. Download: here It is a priority of the Saskatchewan and Canadian government to assess and monitor the health and sustainability of Canada's Forest. The North Digital Land Cover Classification (NDLC) will provide Saskatchewan's contribution to Canada's Earth Observation for Sustainable Development of Forests (EOSD) initiative, helping Canada fulfill it's obligation to the Kyoto Protocol. The NDLC supports the mission and directives of the Saskatchewan provincial government by providing an essential dataset which will enable researchers, natural resource managers and government to assess the health and sustainability of our forests, perform research in the area of climate change, manage natural resources and create policy. The NDLC will be based on a combination of Landsat 5 Thematic Mapper (TM) and Landsat 7 Enhanced Thematic Mapper (ETM+) data representing circa 2000 conditions. The NDLC is being produced through a collaboration of federal, provincial, and territorial governments, agencies and industry. Classification Value Background 0 Agriculture 1 Not Assigned 2 Pasture Upland Herbaceous Graminoid 3 Not Assigned 4 Not Assigned 5 Hardwood Open Canopy 6 Hardwood Closed Canopy 7 Jack Pine Closed Canopy 8 Jack Pine Open Canopy 9 Spruce Closed Canopy 10 Spruce Open Canopy 11 Mixed Hardwoods/Softwoods, Softwood/Hardwood Open and Closed Canopy 12 Treed Rock 13 Recent Burn 14 Revegetating/Regenerating Burn 15 Cutovers 16 Water 17 Marsh 18 Herbaceous Fen 19 Mud Sand Saline 20 Shrub Fen 21 Treed Bog 22 Open Bog 23 Not Assigned 24 Settlements/Roads 25 Barren Land 26 Mixed Softwoods Open and Closed 27 Cloud/Shadow/Haze 28 Unclassified 29 0. Background: Where pixels values are equal to 0 in all channels of satellite image data.  1. Agriculture: Cropland and agricultural clearing areas  2. Not Assigned: Empty Class  3. Pasture Upland Herbaceous Graminoid: Lands containing known pastures, tame or native grasses and herbaceous vegetation. May contain low-lying shrubs with less then 10% tree cover. 4. Not Assigned: Empty Class  5. Not Assigned: Empty Class  6. Hardwood Open Canopy: Trembling Aspen, White Birch, Balsam Poplar composes greater than 75% of species by area, Crown Closure: greater than 10% and less than or equal to 55% (SE crown closure classes A and B).  7. Hardwood Closed Canopy: Trembling Aspen, White Birch, Balsam Poplar composes greater than 75% of species by area, Crown Closure: greater than 55% (SE crown closure classes C and D).  8. Jack Pine Closed Canopy: Jack Pine composes greater than 75% of species by area, Crown Closure: greater than 55% (SE crown closure classes C and D).  9. Jack Pine Open Canopy: Jack Pine composes greater than 75% of species by area, Crown Closure: greater than 10 and less than or equal to 55% (SE crown closure classes C and D).  10. Spruce Closed Canopy: White Spruce, Black Spruce composes greater than 75% of species by area, Crown Closure: greater than 55% (SE crown closure classes C and D).  11. Spruce Open Canopy: White Spruce, Black Spruce composes greater than 75% of species by area, Crown Closure: greater than 10 and less than or equal to 55% (SE crown closure classes C and D).  12. Mixed Hardwoods/Softwoods, Softwood/Hardwood Open and Closed Canopy: An area of hardwood and softwood combinations in which neither hardwood nor softwood account for greater than 75% of species by area and where the crown closure is greater than 10%.  13. Treed Rock: Forest vegetation less than 10%.  14. Recent Burn: An area showing evidence of recent burning natural or prescribed and there is little to no regeneration or revegetation visible.  15. Revegetating/Regenerating Burn: An area showing evidence of natural or prescribed burning and where regeneration or revegetation is visible.  16. Cutovers: An area of deforestation, vegetated and non-vegetated. Ancillary data required to correctly classify due to the anthropogenic land cover/land use class.  17. Water: These areas include lakes, rivers, streams, and reservoirs.  18. Marsh: A periodically wet or continually flooded but non peat-forming area supporting grasses, sedges and reeds.  19. Herbaceous Fen: A wetland area consisting of nutrient rich water and decomposing peat supporting vascular and nonvascular plants grasses, sedges, and reeds.  20. Mud Sand Saline: Water saturated soil, salt water and sand containing no vegetation.  21. Shrub Fen: A wetland area consisting of nutrient rich water and decomposing peat supporting low shrubs, forbs, grass, moss, and a sparse tree cover.  22. Treed Bog: A wetland area consisting of decomposing peat moss, lichen and shrubs with 10% to 25% tree cover of stunted black spruce and tamarack.  23. Open Bog: A wetland area consisting of low nutrient water and decomposing peat moss, lichen, and sparse tree cover.  24. Not Assigned: Empty Class  25. Settlements/Roads: Anthropogenic land cover consisting of urban, commercial, industrial, major roads, highways, surface mines, gravel pits and spoil piles.  26. Barren Land: With the exception of the settlements and Roads class, any area of exposed rock, soil or non-vegetated land.  27. Mixed Softwoods Open and Closed: Jack Pine/Spruce, Spruce/Jack Pine Open and Closed, an area of softwood combinations in which neither Jack Pine or Spruce account for greater than 75% of species by area and where crown closure is greater than 10%.  28. Cloud/Shadow/Haze: An area of cloud, shadow, haze.  29. Unclassified: An area of unidentifiable land cover, indicates no work done/not classified, wrong information, missing data and possible new class greater than 3 pixels. 

Topographic maps produced by Natural Resources Canada conform to the National Topographic System (NTS) of Canada. Indexes are available in three standard scales: 1:1,000,000, 1:250,000 and 1:50,000. The area covered by a given mapsheet is determined by its latitude and longitude. 1:1,000,000 mapsheets are identified by a combination of three numbers (e.g. 098). 1:250,000 mapsheets are identified by a combination of numbers, and letters ranging from A through P (e.g. 098C). Sixteen smaller segments (1 to 16) form blocks used for 1:50,000 mapping (e.g. 098C03).

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

AAFC’s Canadian Ag-Land Monitoring System (CALMS), operational since 2009, was developed by AAFC’s Earth Observation Service (EOS) to deliver weekly NDVI-based maps of crop condition in near-real-time. The CALMS uses data collected by the Moderate Resolution Imaging Spectro-radiometer (MODIS), a sensor mounted onboard NASA’s Terra satellite that has been acquiring data since February 2000. The state-of-the-art radiometric, spectral and spatial resolutions of MODIS Terra make it particularly well-suited for large-scale vegetation mapping and assessment. Crop condition (NDVI) maps are generated weekly by AAFC throughout Canada’s growing season, the period defined as the six-month period stretching from the start of Julian week 12 (end of March) to the end of Julian week 44 (late October). Weeks of the year are defined according to the ISO 8601 week-numbering standard, where weeks start on a Monday and end the following Sunday. CALMS products are generated in the MODIS native Integrated Sinusoidal (ISIN) projection for the region covering the twelve MODIS tiles h09v03 to h14v03 and h09v04 to h14v04.