Land cover information is necessary for a large range of environmental applications related to climate impacts and adaption, emergency response, wildlife habitat, etc. In Canada, a 2008 user survey indicated that the most practical land cover data is provided in a nationwide 30 m spatial resolution format, with an update frequency of five years. In response to this need, the Canada Centre for Remote Sensing (CCRS) has generated a 30 m land cover map of Canada for the base year 2010, as the first of a planned series of maps to be updated every five years, or more frequently. This land cover dataset is also the Canadian contribution to the 30 m spatial resolution 2010 Land Cover Map of North America, which is produced by Mexican, American and Canadian government institutions under a collaboration called the North American Land Change Monitoring System (NALCMS). This paper describes the mapping approach used for generating this land cover dataset for Canada from Thematic Mapper (TM) and Enhanced Thematic Mapper (ETM+) Landsat sensor observations. The innovative part of the mapping approach is the local optimization of the land cover classifier, which has resulted in increased spatial consistency and accuracy. Training and classifying with locally confined reference samples over a large number of partially overlapping areas (i.e., moving windows) ensures the optimization of the classifier to a local land cover distribution, and decreases the negative effect of signature extension. A weighted combination of labels, which is determined by the classifier in overlapping windows, defines the final label for each pixel. Since the approach requires extensive computation, it has been developed and deployed using the Government of Canada’s High-Performance Computing Center (HPC). An accuracy assessment based on 2811 randomly distributed samples shows that land cover data produced with this new approach has achieved 76.60% accuracy with no marked spatial disparities.

Ice-wedge networks underlie polygonal terrain and comprise the most widespread form of massive ground ice in continuous permafrost. Here, we show that climate-driven thaw of hilltop ice-wedge networks is rapidly transforming uplands across Banks Island in the Canadian Arctic Archipelago. Change detection using high-resolution WorldView images and historical air photos, coupled with 32-year Landsat reflectance trends, indicate broad-scale increases in ponding from ice-wedge thaw on hilltops, which has significantly affected at least 1500 km² of Banks Island and over 3.5% of the total upland area. Trajectories of change associated with this upland ice-wedge thermokarst include increased micro-relief, development of high-centred polygons, and, in areas of poor drainage, ponding and potential initiation of thaw lakes. Millennia of cooling climate have favoured ice-wedge growth, and an absence of ecosystem disturbance combined with surface denudation by solifluction has produced high Arctic uplands and slopes underlain by ice-wedge networks truncated at the permafrost table. The thin veneer of thermally-conductive mineral soils strongly links Arctic upland active-layer responses to summer warming. For these reasons, widespread and intense ice-wedge thermokarst on Arctic hilltops and slopes contrast more muted responses to warming reported in low and subarctic environments. Increasing field evidence of thermokarst highlights the inherent climate sensitivity of the Arctic permafrost terrain and the need for integrated approaches to monitor change and investigate the cascade of environmental consequences.

Extreme flood events in recent years in Canada have highlighted the need for historical information to better manage future flood risk. In this paper, a methodology to generate flood maps from Landsat to determine historical inundation frequency is presented for a region along the St-John River, New Brunswick, Canada that experiences annual springtime flooding from snowmelt and river ice. 1985–2016 Landsat data from the USGS archive were classified by combining See5 decision trees to map spectrally variable water due to spring ice and sediment, and image thresholding to map inundated floodplains. Multiple scenes representing each year were overlaid to produce seasonal time-series of spring (March–May) and summer (June–August) maximum annual water extents. Comparisons of annual surface water maps were conducted separately for each season against historical hydrometric water depth as a measure of relative springtime flood severity, and 1 m water masks from digital orthophotos were used to perform a formal accuracy assessment of summer water. Due to Landsat’s 16-day revisit time, peak flood depth was poorly related to flood extent; however, spring depth measured during Landsat acquisitions was significantly related to extent (tau = 0.6, p-value < 0.001). Further, summer maps validated against 30 m water fractions scaled from 1 m water masks were over 97% accurate. Limitations with respect to the assessment of flood extent from depth, timing differences between peak flood depth and extent due to Landsat revisit time and cloud cover, and suggestions to overcome limitations through multi-sensor integration including radar are discussed.

The increasing frequency of flooding worldwide has driven research to improve near real-time flood mapping from remote-sensing data. Improved automation and processing speed to map both open water and vegetated area flooding have resulted from these research efforts. Despite these achievements, flood mapping in urban areas where a significant number of overall impacts are felt remains a challenge. Near real-time data availability, shadowing caused by manmade infrastructure, spatial resolution, and cloud cover inhibiting optical transmission, are all factors that complicate detailed urban flood mapping needed to inform response efforts. This paper uses numerous data sources collected during two major flood events that impacted the same region of Eastern Canada in 2017 and 2019 to test different urban flood mapping approaches presented as case studies in three separate urban boroughs. Cloud-free high-resolution 3 m PlanetLab optical data acquired near peak-flood in 2019 were used to generate a maximum flood extent product for that year. Approaches using new Lidar Digital Elevation Models (DEM)s and water height estimated from nineteen RADARSAT-2 flood maps, point-based flood perimeter observations from citizen geographic information, and simulated traffic camera or other urban sensor network data were tested and verified using independent data. Coherent change detection (CCD) using multi-temporal Interferometric Wide (IW) Sentinel-1 data was also tested. Results indicate that while clear-sky high-resolution optical imagery represents the current gold standard, its availability is not guaranteed due to timely coverage and cloud cover. Water height estimated from 8 to 12.5 m resolution RADARSAT-2 flood perimeters were not sufficiently accurate to flood adjacent urban areas using a Lidar DEM in near real-time, but all nineteen scenes combined captured boroughs that flooded at least once in both flood years. CCD identified flooded boroughs and roughly captured their flood extents, but lacked timeliness and sufficient detail to inform street-level decision-making in near real-time. Point-based flood perimeter observation, whether from in-situ sensors or high-resolution optical satellites combined with Lidar DEMs, can generate accurate full flood extents under certain conditions. Observed point-based flood perimeters on manmade features with low topographic variation produced the most accurate flood extents due to reliable water height estimation from these points.

Mapping the historical occurrence of flood water in time and space provides information that can be used to help mitigate damage from future flood events. In Canada, flood mapping has been performed mainly from RADARSAT imagery in near real-time to enhance situational awareness during an emergency, and more recently from Landsat to examine historical surface water dynamics from the mid-1980s to present. Here, we seek to integrate the two data sources for both operational and historical flood mapping. A main challenge of a multi-sensor approach is ensuring consistency between surface water mapped from sensors that fundamentally interact with the target differently, particularly in areas of flooded vegetation. In addition, automation of workflows that previously relied on manual interpretation is increasingly needed due to large data volumes contained within satellite image archives. Despite differences between data received from both sensors, common approaches to surface water and flooded vegetation mapping including multi-channel classification and region growing can be applied with sensor-specific adaptations for each. Historical open water maps from 202 Landsat scenes spanning the years 1985–2016 generated previously were enhanced to improve flooded vegetation mapping along the Saint John River in New Brunswick, Canada. Open water and flooded vegetation maps were created over the same region from 181 RADARSAT 1 and 2 scenes acquired between 2003–2016. Comparisons of maps from different sensors and hydrometric data were performed to examine consistency and robustness of products derived from different sensors. Simulations reveal that the methodology used to map open water from dual-pol RADARSAT 2 is insensitive to up to about 20% training error. Landsat depicts open water inundation well, while flooded vegetation can be reliably mapped in leaf-off conditions. RADARSAT mapped approximately 8% less open water area than Landsat and 0.5% more flooded vegetation, while the combined area of open water and flooded vegetation agreed to within 0.2% between sensors. Derived historical products depicting inundation frequency and trends were also generated from each sensor’s time-series of surface water maps and compared.

When severe flooding occurs in Canada, the Emergency Geomatics Service (EGS) is tasked with creating and disseminating maps that depict flood extents in near real time. EGS flood mapping methods were created with efficiency and robustness in mind, to allow maps to be published quickly, and therefore have the potential to generate high-repeat water products that can enhance frequent wetland monitoring. The predominant imagery currently used is synthetic aperture radar (SAR) from RADARSAT-2 (R2). With the commissioning phase of the RADARSAT Constellation Mission (RCM) complete, the EGS is adapting its methods for use with this new source of SAR data. The introduction of RCM’s circular-transmit linear-receive (CTLR) beam mode provides the option to exploit compact polarimetric (CP) information not previously available with R2. The aim of this study was to determine the most effective CP parameters for use in mapping open water and flooded vegetation, using current EGS methodologies, and compare these products to those created by using R2 data. Nineteen quad-polarization R2 scenes selected from three regions containing wetlands prone to springtime flooding were used to create reference flood maps, using existing EGS tools. These scenes were then used to simulate 22 RCM CP parameters at different noise floors and spatial resolutions representative of the three RCM beam modes. Using multiple criteria, CP parameters were ranked in order of importance and entered into a stepwise classification procedure, for evaluation against reference R2 products. The top four CP parameters —m-chi-volume or m-delta-volume, RR intensity, Shannon Entropy intensity (SEi), and RV intensity—achieved a maximum agreement with baseline R2 products of upward of 98% across all 19 scenes and three beam modes. Separability analyses between flooded vegetation and other land-cover classes identified four candidate CP parameters—RH intensity, RR intensity, SEi, and the first Stokes parameter (SV0)—suitable for flooded-vegetation-region growing. Flooded-vegetation-region-growing CP thresholds were found to be dependent on incidence angle for each of these four parameters. After region growing using each of the four candidate CP parameters, RH intensity was deemed best to map flooded vegetation, based on our evaluations. The results of the study suggest a set of suitable CP parameters to generate flood maps from RCM data, using current EGS methodologies that must be validated further as real RCM data become available.