RADARSAT Constellation Mission Research Articles

Generating large-scale sea ice motion from Sentinel-1 and the RADARSAT Constellation Mission using the Environment and Climate Change Canada automated sea ice tracking system

Abstract. As Arctic sea ice extent continues to decline, remote sensing observations are becoming even more vital for the monitoring and understanding of sea ice. Recently, the sea ice community has entered a new era of synthetic aperture radar (SAR) satellites operating at C-band with the launch of Sentinel-1A in 2014 and Sentinel-1B (S1) in 2016 and the RADARSAT Constellation Mission (RCM) in 2019. These missions represent five spaceborne SAR sensors that together routinely cover the pan-Arctic sea ice domain. Here, we describe, apply, and validate the Environment and Climate Change Canada automated sea ice tracking system (ECCC-ASITS) that routinely generates large-scale sea ice motion (SIM) over the pan-Arctic domain using SAR images from S1 and RCM. We applied the ECCC-ASITS to the incoming image streams of S1 and RCM from March 2020 to October 2021 using a total of 135 471 SAR images and generated new SIM datasets (7 d 25 km and 3 d 6.25 km) by combining the image stream outputs of S1 and RCM (S1 + RCM). Results indicate that S1 + RCM SIM provides more coverage in the Hudson Bay, Davis Strait, Beaufort Sea, Bering Sea, and directly over the North Pole compared to SIM from S1 alone. Based on the resolvable S1 + RCM SIM grid cells, the 7 d 25 km spatiotemporal scale is able to provide the most complete picture of SIM across the pan-Arctic from SAR imagery alone, but considerable spatiotemporal coverage is also available from 3 d 6.25 SIM products. S1 + RCM SIM is resolved within the narrow channels and inlets of the Canadian Arctic Archipelago, filling a major gap from coarser-resolution sensors. Validating the ECCC-ASITS using S1 and RCM imagery against buoys indicates a root-mean-square error (RMSE) of 2.78 km for dry ice conditions and 3.43 km for melt season conditions. Higher speeds are more apparent with S1 + RCM SIM as comparison with the National Snow and Ice Data Center (NSIDC) SIM product and the Ocean and Sea Ice Satellite Application Facility (OSI SAF) SIM product indicated an RMSE of u=4.6 km d−1 and v=4.7 km d−1 for the NSIDC and u=3.9 km d−1 and v=3.9 km d−1 for OSI SAF. Overall, our results demonstrate the robustness of the ECCC-ASITS for routinely generating large-scale SIM entirely from SAR imagery across the pan-Arctic domain.

Eastern Arctic Sea Ice Sensing: First Results from the RADARSAT Constellation Mission Data

Sea ice monitoring plays a vital role in secure navigation and offshore activities. Synthetic aperture radar (SAR) has been widely used as an effective tool for sea ice remote sensing (e.g., ice type classification, concentration and thickness retrieval) for decades because it can collect data by day and night and in almost all weather conditions. The RADARSAT Constellation Mission (RCM) is a new Canadian SAR mission providing several new services and data, with higher spatial coverage and temporal resolution than previous Radarsat missions. As a very deep convolutional neural network, Normalizer-Free ResNet (NFNet) was proposed by DeepMind in early 2021 and achieved a new state-of-the-art accuracy on the ImageNet dataset. In this paper, the RCM data are utilized for sea ice detection and classification using NFNet for the first time. HH, HV and the cross-polarization ratio are extracted from the dual-polarized RCM data with a medium resolution (50 m) for an NFNet-F0 model. Experimental results from Eastern Arctic show that destriping in the HV channel is necessary to improve the quality of sea ice classification. A two-level random forest (RF) classification model is also applied as a conventional technique for comparisons with NFNet. The sea ice concentration estimated based on the classification result from each region was validated with the corresponding polygon of the Canadian weekly regional ice chart. The overall classification accuracy confirms the superior capacity of the NFNet model over the RF model for sea ice monitoring and the sea ice sensing capacity of RCM.

The RADARSAT Constellation Mission Core Applications: First Results

The Canadian RADARSAT Constellation Mission (RCM) has passed its early operation phase with the performance evaluation being currently active. This evaluation aims to confirm that the innovative design of the mission’s synthetic aperture radar (SAR) meets the expectations of intended users. In this study, we provide an overview of initial results obtained for three high-priority applications; flood mapping, sea ice analysis, and wetland classification. In our study, the focus is on results obtained using not only linear polarization, but also the adopted Compact Polarimetric (CP) architecture in RCM. Our study shows a promising level of agreement between RCM and RADARSAT-2 performance in flood mapping using dual-polarized HH-HV SAR data over Red River, Manitoba, suggesting smooth continuity between the two satellite missions for operational flood mapping. Visual analysis of coincident RCM CP and RADARSAT-2 dual-polarized HH-HV SAR imagery over the Resolute Passage, Canadian Central Arctic, highlighted an improved contrast between sea ice classes in dry ice winter conditions. A statistical analysis using selected sea ice samples confirmed the increased contrast between thin and both rough and deformed ice in CP SAR. This finding is expected to enhance Canadian Ice Service’s (CIS) operational visual analysis of sea ice in RCM SAR imagery for ice chart production. Object-oriented classification of a wetland area in Newfoundland and Labrador by fusion of RCM dual-polarized VV-VH data and Sentinel-2 optical imagery revealed promising classification results, with an overall accuracy of 91.1% and a kappa coefficient of 0.87. Marsh presented the highest user’s and producer’s accuracies (87.77% and 82.08%, respectively) compared to fog, fen, and swamp.

Ocean Surface Wind Speed Retrieval From the RADARSAT Constellation Mission Co- and Cross-Polarization Images Without Wind Direction Input

New techniques for automated retrieval of ocean surface wind speed from the RADARSAT Constellation mission (RCM) without input wind direction have been developed and tested. We collected a large number of collocated and coincided ocean buoy wind measurements and RCM co- and cross-polarization observations acquired over the West and East coasts of Canada from February 1, 2020 to July 25, 2021: 1190 data points for RCM ScanSAR 50 m (SC50M) VV–VH, 3819 points for SC50M HH–HV, and 397 points for ScanSAR Low Noise (SCLN) HH–HV images. The observations captured a wide range of wind conditions with the maximum wind speeds of 22.3, 23.8, and 19.7 m/s for SC50M VV–VH, SC50M HH–HV, and SCLN HH–HV, respectively. For these three types of RCM data, the ocean surface wind speed was modeled as a function of the co-polarization (VV or HH) and cross-polarization backscatter (VH or HV) as well as the noise floor and the incidence angle. The models were tested on completely independent subsets of data, and their performance was compared against CMOD5.N (for VV) and CMODH (for HH) models which require input wind direction. The root-mean square errors (RMSEs) for the proposed models in the testing subsets are 1.20, 1.52, and 1.26 m/s, while CMOD models showed 1.37, 1.91, and 3.06 m/s for SC50M VV–VH, SC50M HH–HV, and SCLN HH–HV, respectively. The proposed models are being integrated in various Environment and Climate Change Canada applications including ice and wind data assimilation systems and the National SAR Winds program.

Using the Pauli Scattering Mechanisms for Analysis of Polarimetric SAR Images from RADARSAT Constellation and TanDEM-X Missions

This work describes techniques for extraction of information in polarimetric images from RADARSAT constellation mission (RCM) and TanDEM-X (TDX) and TerraSAR-X (TSX) satellites. Using the boundary conditions of Maxwell's equations it is shown how the single bounce and double bounce scattering can be extracted from quad-pol and compact polarimetric (CP) SAR. We develop new equations to calculate the Pauli decomposition from simulated CP data using RCM quad-pol data. The RCM CP option is circular transmit, linear receive (CTLR). We find that using quad-pol data performs better than using CP data for Pauli decomposition, especially extraction of the HV component. Estimation of single bounce and double bounce scattering from RCM CP data is feasible and yields more information about scattering properties than the CP components alone. We show how more information can be extracted from bistatic TSX/TDX quad-pol data by estimation of the coherence of the different Pauli components and optimization of the interferometric coherence using ships as detection objects. This is a unique opportunity to use optimized coherence on moving ships for characterizing their structure.

Sea-Ice Mapping of RADARSAT-2 Imagery by Integrating Spatial Contexture With Textural Features

Mapping different types of sea ice that form, grow, and melt in polar oceans is essential for shipping navigation, climate change modeling, and local community safety. Currently, ice charts are manually generated by analysts at the Canadian Ice Service (CIS) based on dual-polarized RADARSAT-2/RADARSAT Constellation Mission (RCM) imagery on a daily basis. Inspired by the demand for a computer-based mapping system, we have developed an automatic sea ice classification method that integrates spatial contexture (unsupervised segmentation) with textural features (supervised pixel-level labeling). First, the full-scene image is oversegmented, and the segments are merged into homogeneous regions across the entire scene. Second, pixel-based classifiers (support vector machine, random forest) are compared for their ability to label the generated homogeneous regions. Finally, the segmentation and labeling are combined using a proposed energy function. The proposed method was tested on 18 dual-polarization RADARSAT-2 scenes acquired over the Beaufort Sea. This dataset contains water, young ice, first-year ice, and multi-year ice covering melt, summer, and freeze-up seasons. The proposed method obtains an average classification accuracy of 86.33% based on the leave-one-out validation. The experimental results show that the proposed method achieves promising classification results in both quantity and quality measurements compared to benchmark methods. The robustness against incidence angle variance indicates that the proposed method is well-qualified for operational sea ice mapping.

A multi-year cross-validation experiment for estimating rice plant area index (PAI) over the JECAM-India test site from simulated RADARSAT constellation mission (RCM) compact polarimetric SAR data

ABSTRACT Using the cross-validation approach, strategies for estimating biophysical parameters are still pre-operational with synthetic aperture radar (SAR) data. In this regard, the Joint Experiment for Crop Assessment and Monitoring (JECAM) SAR inter-comparison experiments provide an opportunity for the potential implementation of cross-validation strategies for biophysical parameters retrieval utilizing the next-generation compact polarimetric (CP) modes available from the RADARSAT Constellation Mission (RCM). This work first uses the conventional semi-empirical Water Cloud Model (WCM) modified by exploiting the scattering power decompositions of CP measurements to estimate the Plant Area Index (PAI) for rice. The modified WCM (MWCM) is then inverted using the scattering power components from the decomposition. We compare the PAI estimates using MWCM- between the estimates obtained from (1) the conventional WCM using the RH and RV backscatter intensities and (2) MWCM- decomposition scattering powers. We exploit a time series of simulated compact-pol SAR data over the JECAM test site in Vijayawada, India, throughout 2018 and 2019. We use the C-band RADARSAT-2 full-pol data to simulate the RADARSAT Constellation Mission (RCM) compact-pol mode data. Utilizing the advantage of systematically collected multi-year SAR data and in-situ measurements, the present research also assesses the calibrated model transferability performances to another data set and cross-validation of a model in a multi-year experiment setting. The comparative analysis indicates potential improvements in PAI estimation with MWCM- scattering powers. A high range of correlation coefficient () between the estimated and observed PAI is observed with good Root Mean Square Error (RMSE) of m2 m−2, and Mean Absolute Error (MAE) of m2 m−2.

Sentinel-1 and RADARSAT Constellation Mission InSAR Assessment of Slope Movements in the Southern Interior of British Columbia, Canada

Landslides are the most common natural hazard in British Columbia. The province has recorded the largest number of historical landslide fatalities in Canada, and damage to infrastructure comes at a great cost. In order to understand the potential impacts of landslides, radar remote sensing has become a cost-effective method for detecting downslope movements. This study investigates downslope movements in the Southern Interior of British Columbia, Canada, with Sentinel-1 and RADARSAT Constellation Mission (RCM) interferometric synthetic aperture radar (InSAR) data. The 2-dimensional time-series analysis with Sentinel-1 ascending and descending InSAR pairs from October 2017 to June 2021 observed distinct earthflow movements of up to ~15 cm/year in the east–west direction. The Grinder Creek, Red Mountain, Yalakom River, and Retaskit Creek earthflows previously documented are still active, with east–west movements of ~30 cm over the past four years. New RCM data acquired from June 2020 to September 2020 with a 4-day revisit capability were compared to 12-day Sentinel-1 InSAR pairs. The 4-day RCM InSAR pairs at higher spatial resolution showed better performance by detecting relatively small-sized slope movements within a few hundred meters, which were not clearly observed by Sentinel-1. The temporal variabilities observed from the RCM InSAR showed great potential for observing detailed slope movements within a narrower time window.

Benchmarked RADARSAT-2, SENTINEL-1 and RADARSAT Constellation Mission Change-Detection Monitoring at North Slide, Thompson River Valley, British Columbia: ensuring a Landslide-Resilient National Railway Network

In this research note, we demonstrate the applicability of interferometric analyses (InSAR) of RADARSAT 2 (RS2), SENTINEL 1 (S1) and RADARSAT Constellation Mission (RCM) datasets to characterize and monitor landslides along a high-risk section of the national railway transportation corridor traversing the Thompson River valley, British Columbia. As a geomorphically active landform, the North Slide is an ideal case study for field-testing and evaluating slope change-detection monitoring incorporating satellite, aerial and ground-based geospatial technologies. RS2, S1 and RCM InSAR datasets provide valuable baseline spatial and temporal information on movement of the landslide near critical railway infrastructure when benchmarked with real-time kinematic (RTK) global navigation satellite system (GNSS) measurements, uninhabited aerial vehicle (UAV) photogrammetry, bathymetric soundings, and ground observations. We demonstrate that monitoring unstable slopes and infrastructure at risk with multiple high spatial- and temporal-resolution satellite SAR platforms is a cost-effective natural hazard management practice that also provides important geoscience information to help develop appropriate mitigation and climate adaptation measures.

Detection of Actionable Oil using RADARSAT-2 and Simulated RADARSAT Constellation Mission Data

Abstract 688970 For oil spill response, one of the key parameters is detection of actionable oil. Actionable oil (AO), which tends to be thick, emulsified oil, refers to oil that can be cleaned-up versus non-actionable oil (sheen) which cannot be readily cleaned-up. Previous studies by MDA of a controlled oil spills in the North Sea have shown the capability of RADARSAT-2 quad polarized data to detect AO using the entropy parameter (H) derived from the Cloude-Pottier decomposition. H → 0 for oil-free water and H → 1 in the presence of an oil slick. To further test the detection of AO, RADARSAT-2, ASTER, WorldView imagery and in situ measurement were acquired April 25, 2017 at the MC20 site off Louisiana. Five oil thickness classes ranging from 1 micron to 100 microns were derived from a maximum likelihood classifier based on the ASTER and WorldView images and in situ samples. For oil-free water, the average H was 0.13 and 0.62 for the slick. There was good correlation between the variability of H and the oil thickness classes. Specifically, larger H was correlated with oil thickness in the 50 – 200 micron range and smaller H was correlated with oil in the 1 micro range. Not surprising there was overlap in H for area were the oil thickness was ~ 1 micron and the area deemed to be slick-free. The results indicate that actionable oil can be discriminated from non-actionable oil based on the relative difference of H. Although these results are encouraging, one of the operational limitations is based on the use of the relatively small swath-width (25 km – 50 km) of the RADARSAT-2 quad polarized mode. The impact of the swath-width can be mitigated with data from the RADARSAT Constellation Mission (RCM). The RCM has a compact polarimetry (CP) mode that provides more polarimetric information than dual polarized modes, less than quad polarized modes, but is available for swath widths up to 500 km. Analysis of simulated SC50 RCM data (50 m resolution, 350 km swath width) derived from the aforementioned RADARSAT-2 image shows similar oil-slick variability that was observed in the RADARSAT-2 image and hence the capability to detect AO.

Compact Polarimetry Synthetic Aperture Radar Ocean Wind Retrieval: Model Development and Validation

AbstractWe have developed C-band compact polarimetry geophysical model functions for RADARSAT Constellation Mission ocean surface wind speed retrieval. A total of 1594 RADARSAT-2 images acquired in quad-polarization SAR imaging mode were collocated with in situ buoy observations. This dataset is first used to simulate compact polarimetric data and to examine their dependencies on radar incidence angle and wind vectors. We find that right circular transmit, right circular receive (RR-pol) radar backscatters are less sensitive to incidence angles and wind directions but are more dependent on wind speeds, compared to right circular transmit, horizontal receive (RH-pol), right circular transmit, vertical receive (RV-pol), and right circular transmit, left circular receive (RL-pol). Subsequently, the matchup data pairs are used to derive the coefficients of the transfer functions for the proposed compact polarimetric geophysical model (CMOD) functions, and to validate the associated wind speed retrieval accuracy. Statistical comparisons show that the retrieved wind speeds from CMODRH, CMODRV, CMODRL, and CMODRR are in good agreement with buoy measurements, with root-mean-square errors of 1.38, 1.51, 1.47, and 1.25 m s−1, respectively. The results suggest that compact polarimetry is a good alternative to linear polarization for wind speed retrieval. CMODRR is more appropriate to retrieve high wind speeds than CMODRH, CMODRV or CMODRL.

Multi-Source EO for Dynamic Wetland Mapping and Monitoring in the Great Lakes Basin

Wetland managers, citizens and government leaders are observing rapid changes in coastal wetlands and associated habitats around the Great Lakes Basin due to human activity and climate variability. SAR and optical satellite sensors offer cost effective management tools that can be used to monitor wetlands over time, covering large areas like the Great Lakes and providing information to those making management and policy decisions. In this paper we describe ongoing efforts to monitor dynamic changes in wetland vegetation, surface water extent, and water level change. Included are assessments of simulated Radarsat Constellation Mission data to determine feasibility of continued monitoring into the future. Results show that integration of data from multiple sensors is most effective for monitoring coastal wetlands in the Great Lakes region. While products developed using methods described in this article provide valuable management tools, more effort is needed to reach the goal of establishing a dynamic, near-real-time, remote sensing-based monitoring program for the basin.

CP-IRGS: A Region-Based Segmentation of Multilook Complex Compact Polarimetric SAR Data

The Canadian RADARSAT constellation mission (RCM) is represented by three synthetic aperture radar (SAR) satellites, each of which includes a compact polarimetry (CP) mode. CP is advantageous because it provides increased backscatter information relative to single and conventional dual-polarized modes and has larger swath widths relative to a quad polarization mode. CP captures single-look complex data which can be used to derive the multilook complex (MLC) coherence matrix, or, equivalently, the Stokes vector data of the backscattered field. The challenge is to develop computer vision algorithms that can be used to effectively segment the scene using this new data source. An unsupervised region-based segmentation approach has been designed and implemented that utilizes the complex Wishart distribution characteristic of the MLC CP data. The segmentation method is based on the iterative region growing with semantics algorithm originally designed for single and dual pol intensity SAR data. The algorithm has been tested using both simulated CP SAR images and a pair of available quad polarization SAR images. The results demonstrate that the CP-IRGS algorithm provides more accurate segmentation images than those using only the RH and RV channel intensity images.

Wetland Hydroperiod Change Along the Upper Columbia River Floodplain, Canada, 1984 to 2019

Increasing air temperatures and changing hydrological conditions in the mountainous Kootenay Region of British Columbia, Canada are expected to affect floodplain wetland extent and function along the Columbia River. The objective of this study was to determine the seasonally inundated hydroperiod for a floodplain section (28.66 km2) of the Upper Columbia River wetlands complex using time series satellite image observations and binary open water mask extraction. A mid pixel resolution (30 m) optical satellite image time series of 61 clear sky scenes from the Landsat Thematic Mapper (TM) and Operational Land Imager (OLI) sensors were used to map temporal variations in floodplain open water wetland extent during the April to October hydrologically active season from 1984 to 2019 (35 years). The hydroperiod from the first 31 scenes (T1: 18 years) was compared to the second 30 (T2: 16 years) to identify changes in the permanent and seasonal open water bodies. The seasonal variation in open water extent and duration was similar across the two time periods but the permanent water body extent diminished by ~16% (or ~3.5% of the floodplain). A simple linear model (r2 = 0.87) was established to predict floodplain open water extent as a function of river discharge downstream of the case study area. Four years of Landsat Multi-Spectral Scanner (MSS) data from 1992 to 1995 (12 scenes) were examined to evaluate the feasibility of extending the hydroperiod record back to 1972 using lower resolution (60 m) archive data. While the MSS hydroperiod produced a similar pattern of open water area to duration to the TM/OLI hydroperiod, small open water features were omitted or expanded due to the lower resolution. While MSS could potentially extend the TM/OLI hydroperiod record, this was not performed as the loss of features like the river channel diminished its value for change detection purposes. Radarsat 2 scenes from 2015 to 2019 were examined to evaluate the feasibility of continued mountain valley hydroperiod monitoring using higher spatial and temporal resolution sensors like the Radarsat Constellation Mission (RCM). From the available horizontal transmit/receive (HH) single polarization sample set (8 scenes), the hydroperiod pattern of open water extent to duration was similar to the longer Landsat time series and possessed greater feature detail, but it was significantly reduced in seasonal inundation area due to the systematic omission of open water areas containing emergent vegetation. However, accepting that differences exist in sensor-based hydroperiod attributes, the higher temporal resolution of RCM will be suited to mountain floodplain inundation monitoring and open water hydroperiod analysis.

C-band synthetic aperture radar (SAR) imagery for the classification of diverse cropping systems

ABSTRACT Cloudy conditions reduce the utility of optical imagery for crop monitoring. New constellations of satellites – including the RADARSAT Constellation Mission (RCM) and Sentinel-1A/B, both available under free and open data policies – can be used to create stacks of dense seasonal C-band Synthetic Aperture Radar (SAR) data. Yet to date, the contribution of SAR imagery to operational crop mapping is often limited to that of a gap-filler, compensating for optical data obscured by clouds. The Joint Experiment for Crop Assessment and Monitoring (JECAM) SAR Inter-Comparison Experiment is a multi-year, multi-partner project focused on evaluating methods for SAR-based crop classification. Stacks of dense time-series SAR imagery, from RADARSAT-2 and Sentinel-1 satellites, were acquired for 10 sites located in six countries. Decision Tree (DT) and Random Forest (RF) classification methodologies were applied to these SAR data-stacks, as well as to data-stacks of optical only, and optimized SAR/optical data combinations. For the dense time-series SAR stacks, overall classification accuracies above 85% and 80% were obtained for 6 of 10 and 8 of 10 sites, respectively. For maize, the SAR-only data delivered user’s and producer’s accuracies greater than 90% for half the sites. For soya bean, accuracies greater than 80% were reported for 5 of 9 sites and classification accuracies were greater than 80% for wheat on half the sites. Classification results were influenced by the mix and number of agriculture classes present at each site, the available SAR imagery, as well as the training and validation data sets for individual crop types. These results have important operational implications for regions of the world dominated by cloudy conditions and the lack of adequate amounts of optical imagery to support satellite-based crop monitoring.

NEW VEGETATION INDICES FOR FULL AND COMPACT POLARIMETRIC SAR DATA: IN PREPARATION FOR THE RADARSAT CONSTELLATION MISSION (RCM)

Abstract. In this paper, we present two radar vegetation indices for full-pol and compact-pol SAR data, respectively. Both are derived using the notion of a geodesic distance between observation and well-known scattering models available in the literature. While the full-pol version depends on a generalized volume scattering model, the compact-pol version uses the ideal depolariser to model the randomness in the vegetation. We have utilized the RADARSAT Constellation Mission (RCM) time-series data from the SAMPVEX16-MB campaign in the Manitoba region of Canada for comparing and assessing the indices in terms of the change in the biophysical parameters as well. The compact-pol data for comparison is simulated from the full-pol RCM time series. Both the indices show better performance at correlating with biophysical parameters such as Plant Area Index (PAI) and Volumetric Water Content (VWC) for wheat and soybean crops, in comparison to the state-of-art Radar Vegetation Index (RVI) of Kim and van Zyl. These indices are timely for the upcoming release of the data from the RCM, which will provide data in both full and compact-pol modes, aimed at better crop monitoring from space.

Hybrid Compact Polarimetric SAR for Environmental Monitoring with the RADARSAT Constellation Mission

Canada’s successful space-based earth-observation (EO) radar program has earned widespread and expanding user acceptance following the launch of RADARSAT-1 in 1995. RADARSAT-2, launched in 2007, while providing data continuity for its predecessor’s imaging capabilities, added new polarimetric modes. Canada’s follow-up program, the RADARSAT Constellation Mission (RCM), launched in 2019, while providing continuity for its two predecessors, includes an innovative suite of polarimetric modes. In an effort to make polarimetry accessible to a wide range of operational users, RCM uses a new method called hybrid compact polarization (HCP). There are two essential elements to this approach: (1) transmit only one polarization, circular; and (2) receive two orthogonal polarizations, for which RCM uses H and V. This configuration overcomes the conventional dual and full polarimetric system limitations, which are lacking enough polarimetric information and having a small swath width, respectively. Thus, HCP data can be considered as dual-pol data, while the resulting polarimetric classifications of features in an observed scene are of comparable accuracy as those derived from the traditional fully polarimetric (FP) approach. At the same time, RCM’s HCP methodology is applicable to all imaging modes, including wide swath and ScanSAR, thus overcoming critical limitations of traditional imaging radar polarimetry for operational use. The primary image data products from an HCP radar are different from those of a traditional polarimetric radar. Because the HCP modes transmit circularly polarized signals, the data processing to extract polarimetric information requires different approaches than those used for conventional linearly polarized polarimetric data. Operational users, as well as researchers and students, are most likely to achieve disappointing results if they work with traditional polarimetric processing tools. New tools are required. Existing tutorials, older seminar notes, and reference papers are not sufficient, and if left unrevised, could succeed in discouraging further use of RCM polarimetric data. This paper is designed to provide an initial response to that need. A systematic review of studies that used HCP SAR data for environmental monitoring is also provided. Based on this review, HCP SAR data have been employed in oil spill monitoring, target detection, sea ice monitoring, agriculture, wetland classification, and other land cover applications.

Compact Polarimetry Response to Modeled Fast Sea Ice Thickness

Compact Polarimetric (CP) Synthetic Aperture Radar (SAR) is expected to gain more and more ground for Earth observation applications in the coming years. This comes in light of the recently launched RADARSAT Constellation Mission (RCM), which uniquely provides CP SAR imagery in operational mode. In this study, we present observations about the sensitivity of CP SAR imagery to thickness of thermodynamically-grown fast sea ice during early ice growth (September–December 2017) in the Resolute Bay area, Canadian Central Arctic. Fast ice is most suitable to use for this preliminary study since it exhibits only thermodynamic growth in absence of ice mobility and deformation. Results reveal that ice thickness up to 30 cm can be retrieved using several CP parameters from the tested set. This ice thickness corresponds to the thickness of young ice. We found the surface scattering mechanism to be dominant during the early ice growth, exposing an increasing tendency up to 30 cm thickness with a correlation coefficient with the thickness equal to 0.86. The degree of polarization was found to be the parameter with the highest correlation up to 0.95. While thickness retrieval within the same range is also possible using parameters from Full Polarimetric (FP) SAR parameters as shown in previous studies, the advantage of using CP SAR mode is the much larger swath coverage, which is an operational requirement.

Automated surface water detection from space: a Canada-wide, open-source, automated, near-real time solution

The goal of this research was to develop a fully automated method to map open water extent that is operationally practical on a national scale. Such a system needs to produce acceptable results in all regions of the country and particularly in the Prairie Potholes Region where understanding water surface dynamics is important for predicting flooding, agriculture/water availability and for evaporation calculations in weather models. A system was developed to automate, ingest, and process Radarsat-2 (RS2) imagery, from which mapping open water body extents in near-real time was carried out using a machine learning classification technique. A Random Forest classification algorithm was trained using the data extracted from the Global Surface Water (GSW) occurrence dataset. The GSW occurrence thresholds used to extract the training data were examined and there was little influence of uncertainty on the classification. The quality of classifications generated from RS2 Fine Wide mode imagery improved with increasing incidence angle. All Fine Quad incident angles produced acceptable results, but Standard and Wide mode imagery produced results below the accuracy thresholds deemed acceptable for this operational solution. Validation was carried out by comparing mapped water extents to temporally coincident high resolution multi-spectral imagery and to the USGS Global Land Cover Characteristics dataset, that is currently used as a land-water mask by ECCC in weather numerical weather modelling. The system that has been developed will allow new image datasets (e.g. Radarsat Constellation Mission) or training data that becomes available to be used to improve models. The open source code will be made available on Github.

A Radar Vegetation Index for Crop Monitoring Using Compact Polarimetric SAR Data

Crop growth monitoring using compact-pol synthetic aperture radar (CP-SAR) data is gaining attention with the rapid advancements toward operational applications. In this article, we propose a vegetation index for compact polarimetric (CP) SAR data [compact-pol radar vegetation index (CpRVI)]. The CpRVI is derived using the concept of a geodesic distance between the Kennaugh matrices projected on a unit sphere. This distance is utilized to compute a similarity measure between the observed Kennaugh matrix and the Kennaugh matrix of an ideal depolarizer (a realization of vegetation canopy). The similarity measure is then modulated with a scaled quantity derived from the scattering power ratio of the same and opposite sense polarization with respect to the transmitted circular polarization. In this article, we utilize time-series-simulated RADARSAT Constellation Mission (RCM) compact-pol SAR data (RH-RV) obtained from the full-pol RADARSAT-2 observations during the soil moisture active passive (SMAP) validation experiment 2016 (SMAPVEX16-MB) campaign in Manitoba, Canada, to assess the proposed vegetation index. Among the various crops grown in this region, in particular, we analyze the growth stages of wheat and soybean due to their different canopy structures. A temporal analysis of the proposed CpRVI with crop biophysical parameters [the plant area index (PAI) and vegetation water content (VWC)] at different phenological stages confirms the trend of CpRVI with the plant growth. Nevertheless, variations of CpRVI values are apparent with different plant densities for both the crop types. Also, the linear regression analysis confirms that the CpRVI values significantly correlate with PAI (r = 0.72 and 0.85) and VWC (r = 0.62 and 0.75) for both wheat and soybean. We observed good retrieval of PAI and VWC for both wheat and soybean.