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

Temporal analysis of changes in Winnipeg, Manitoba, based on GeoAI features automatically extracted from satellite images acquired in 2013 and 2023. Simple geospatial analysis enables the detection of features present in 2023 that were not already there in 2013. The addition of new buildings is a good indicator of urban development and/or sprawl. Complementarily, an analysis of changes in the forest coverage from the GeoAI datasets is done. This analysis reflects the gains and losses between both dates. ‌ GeoAI enables temporal coverage of various areas in Canada, thus providing a useful tool for change detection and trend analysis at high resolution. While the series is still fairly new, and such examples are limited for the time being, NRCan strives to gradually increase its GeoAI data offering for both spatial and temporal coverage. For more information about the GeoAI - GeoBase Series, please visit the following link: https://open.canada.ca/data/en/dataset/74738ff5-5367-5958-9aee-98fffdcd1876

The 1 cm resolution vegetation digital height model was extracted using a bare earth model and digital surface model (DSM) derived from unmanned aerial vehicle (UAV) imagery acquired from a single day survey on July 28th 2016, in Cambridge Bay, Nunavut. The mapping product covers 525m2 and 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

This layer represents Land use polygons as determined by a combination of analytic techniques, mostly using Landsat 5 image mosaics . BTM 1 was done on a federal satellite image base that was only accurate to about 250m. The images were geo-corrected, not ortho-corrected, so there is distortion in areas of high relief. This is not a multipart feature

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

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.

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

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

Organic soils in the boreal forest commonly store as much carbon as the vegetation above ground. While recent efforts through the National Forest Inventory has yielded new spatial datasets of forest structure across the vast area of Canada’s boreal forest, organic soils are poorly mapped. In this geospatial dataset, we produce a map primarily of forested and treed peatlands, those with more than 40 cm of peat accumulation and over 10% tree canopy cover. National Forest Inventory ground plots were used to identify the range of forest structure that corresponds to the presence of over 40 cm of peat soils. Areas containing that range of forest cover were identified using the National Forest Inventory k-NN forest structure maps and assigned a probability (0-100% as integer) of being a forested or treed peatland according to a statistical model. While this mapping product captures the distribution of forested and treed peatlands at a 250 m resolution, open, completely treeless peatlands are not fully captured by this mapping product as forest cover information was used to create the maps. The methodology used in the creation of this product is described in: Thompson DK, Simpson BN, Beaudoin A. 2016. Using forest structure to predict the distribution of treed boreal peatlands in Canada. Forest Ecology and Management, 372, 19-27. https://cfs.nrcan.gc.ca/publications?id=36751 This distribution uses an updated forest attribute layer current to 2011 from: 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 Additionally, this distribution varies slightly from the original published in 2016 in that here slope data is derived from the CDEM: https://open.canada.ca/data/en/dataset/7f245e4d-76c2-4caa-951a-45d1d2051333 The above peatland probability map was further processed to delineate bogs vs fens (based on mapped Larix content via the k-NN maps), as well as an approximation of the extent of open peatlands using EOSD data. The result is a 9-type peatland map with a more complete methodology as detailed in: Webster, K. L., Bhatti, J. S., Thompson, D. K., Nelson, S. A., Shaw, C. H., Bona, K. A., Hayne, S. L., & Kurz, W. A. (2018). Spatially-integrated estimates of net ecosystem exchange and methane fluxes from Canadian peatlands. Carbon Balance and Management, 13(1), 16. https://doi.org/10.1186/s13021-018-0105-5 In plain text, the legend for the 9-class map is as follows: value="0" label="not peat" alpha="0" value="1" label="Open Bog" alpha="255" color="#0a4b32" value="2" label="Open Poor Fen" alpha="255" color="#5c5430" value="3" label="Open Rich Fen" alpha="255" color="#792652" value="4" label="Treed Bog" alpha="255" color="#6a917b" value="5" label="Treed Poor Fen" alpha="255" color="#aba476" value="6" label="Treed Rich Fen" alpha="255" color="#af7a8f" value="7" label="Forested Bog" alpha="255" color="#aad7bf" value="8" label="Forested Poor Fen" alpha="255" color="#fbfabc" value="9" label="Forested Rich Fen" alpha="255" color="#ffb6db" This colour scale is given in qml/xml format in the resources below. The 9-type peatland map from Webster et al 2018 was further refined slightly following two simple conditions: (1) any 250-m raster cell with greater than 40% pine content is classified as upland (non-peat); (2) all 250-m raster cells classified as water or agriculture via the NRCan North American Land Cover Monitoring System (https://doi.org/10.3390/rs9111098) is also classified as non-peatland (value of zero in the 9-class map. This mapping scheme was used at a regional scale in the following paper: Thompson, D. K., Simpson, B. N., Whitman, E., Barber, Q. E., & Parisien, M.-A. (2019). Peatland Hydrological Dynamics as A Driver of Landscape Connectivity and Fire Activity in the Boreal Plain of Canada. Forests, 10(7), 534. https://doi.org/10.3390/f10070534 And is reproduced here at a national scale. Note that this mapping product does not fully capture all permafrost peatland features covered by open canopy spruce woodland with lichen ground cover. Nor are treeless peatlands near the northern treeline captured in the training data, resulting in unknown mapping quality in those regions.