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

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

The AAFC Infrastructure Flood Mapping in Saskatchewan 20 centimeter colour orthophotos is a collection of georeferenced color digital orthophotos with 20 cm pixel size. The imagery was delivered in GeoTIF and ECW formats. The TIF and ECW mosaics were delivered in the same 1 km x 1 km tiles as the LiDAR data, and complete mosaics for each area in MrSID format were also provided. The digital photos were orthorectified using the ground model created from the DTM Key Points. With orthorectification, only features on the surface of the ground are correctly positioned in the orthophotos. Objects above the surface of the ground, such as building rooftops and trees, may contain horizontal displacement due to image parallax experienced when the photos were captured. This is sometimes apparent along the cut lines between photos. For positioning of above-ground structures it is recommended to use the LiDAR point clouds for accurate horizontal placement.

This land cover data set was derived from the Advanced Very High Resolution Radiometer (AVHRR) sensor operating on board the United States National Oceanic and Atmospheric Administration (NOAA) satellites. Information on the NOAA series of satellites can be found at www.noaa.gov/satellites.html The vegetation and land cover information set has been classified into twelve categories. Information on the classification of the vegetation and land cover, raster to vector conversion, generalization for cartographic presentations is included in the paper "The Canada Vegetation and Land Cover: A Raster and Vector Data Set for GIS Applications - Uses in Agriculture" (https://geogratis.cgdi.gc.ca/download/landcover/scale/gis95ppr.pdf). A soil quality evaluation was obtained by cross-referencing the AVHRR information with Census of Agriculture records and biophysical (Soil Landscapes of Canada) data and is also included in the above paper. AVHRR Land Cover Data approximates a 1:2M scale and was done originally for Agriculture Canada. The projection used is Lambert Conformal Conic (LCC) 49/77 with origin at 49N 95W.

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

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

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.

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.

The 1 cm resolution digital surface model (DSM) 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 DSM 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