Statistical comparison of SAR backscatter from icebergs embedded in sea ice and in open water using RADARSAT-2 images of in Newfoundland waters and the Davis Strait

Icebergs are considered a threat to marine operations. Satellite monitoring of icebergs is one option to aid in the development of iceberg hazard maps. Satellite synthetic aperture radar (SAR) is an obvious choice because of its relative weather independence, day and night operation. Nonetheless, the detection of icebergs in SAR can be a challenge, particularly with high iceberg areal density, heterogeneous background clutter and the presence of sea ice. This thesis investigates and compares polarimetric signatures of icebergs embedded in sea ice and icebergs in open water. In this thesis, RADARSAT-2 images have been used for analysis, which was acquired over locations near the coastline (approximately 3-35 km) of the islands of Newfoundland and Greenland. All icebergs considered here are in the lower incident angle range (below 30 degrees) of the SAR acquisition geometry. For analysis, polarimetry parameters such as co- (HH) and cross- (HV) polarization and several decomposition techniques, specifically Pauli, Freeman-Durden, Yamaguchi, Cloud-Pottier and van Zyl classification, have been used to determine the polarimetric signatures of icebergs and sea ice. Statistical hypothesis tests were used to determine the differences among backscatters from different icebergs. Statistical results tend to show a dominant surface scattering mechanism for icebergs. Moreover, icebergs in open water produce larger volume scatter than icebergs in sea ice, while icebergs in sea ice produce larger surface scatter than icebergs in open water. In addition, there appear to be minor observable differences between icebergs in Greenland and icebergs in Newfoundland

Electromagnetic backscatter modelling of icebergs at c-band in an ocean environment

This thesis outlines the development of an electromagnetic (EM) backscatter model of icebergs. It is a necessary first step for the generation of in-house synthetic aperture radar (SAR) data of icebergs to support optimum iceberg/ship classifier design. The EM modelling was developed in three stages. At first, an EM backscatter model was developed to generate simulated SAR data chips of iceberg targets at small incidence angles. The model parameters were set to mimic a dual polarized dataset collected at C-Band with the Sentinel-1A satellite. The simulated SAR data chips were compared with signatures and radiometric properties of the satellite data, including total radar cross section (TRCS). A second EM model was developed to mimic the parameters of a second SAR data collection with RADARSAT-2; this second data collection was at larger incidence angles and was fully polarimetric (four channels and interchannel phase). The full polarimetric SAR data allowed for a comparison of modelled TRCS and polarimetric decompositions. Finally, the EM backscatter models were tested in the context of iceberg/ship classification by comparing the performance of various computer vision classifiers using both simulated and real SAR image data of iceberg and vessel targets. This step is critical to check the compatibility of simulated data with the real data, and the ability to mix real and simulated SAR imagery for the generation of skilled classifiers. An EM backscatter modelling tool called GRECOSAR was used for the modelling work. GRECOSAR includes the ability to generate small scenes of the ocean using Pierson-Moskowitz spectral parameters. It also allows the placement of a 3D target shape into that ocean scene. Therefore, GRECOSAR is very useful for simulating SAR targets, however it can only model single layer scattering from the targets. This was found to be limiting in that EM scattering throughout volume of the iceberg could not be generated. This resulted in EM models that included only surface scattering of the iceberg. In order to generate realistic SAR scenes of icebergs on the ocean, 3D models of icebergs were captured in a series of field programs off the coast of Newfoundland and Labrador, Canada. The 3D models of the icebergs were obtained using a light detection and ranging (LiDAR) and multi-beam sonar data from a specially equipped vessel by a team of C-CORE. While profiling the iceberg targets, SAR images from satellites were captured for comparison with the simulated SAR images. The analysis of the real and simulated SAR imagery included comparisons of TRCS, SAR signature morphology and polarimetric decompositions of the targets. In general, these comparisons showed a good consistency between the simulated and real SAR scene. Simulations were also performed with varying target orientation and sea conditions (i.e., wind speed and direction). A wide variability of the TRCS and SAR signature morphology was observed with varying scene parameters. Icebergs were modelled using a high dielectric constant to mimic melting iceberg surfaces as seen during field work. Given that GRECOSAR could only generate surface backscatter, a mathematical model was developed to quantify the effect of melt water on the amount of surface and volume backscatter that might be expected from the icebergs. It was found that the icebergs in a high state of melt should produce predominantly surface scatter, thus validating the use of GRECOSAR for icebergs in this condition. Once the simulated SAR targets were validated against the real SAR data collections, a large dataset of simulated SAR chips of ships and icebergs were created specifically for the purpose of target classification. SAR chips were generated at varying imaging parameters and target sizes and passed on to an iceberg/ship classifier. Real and simulated SAR chips were combined in varying quantities (or targets) resulting in a series of different classifiers of varying skill. A good agreement between the classifier’s performance was found. This indicates the compatibility of the simulated SAR imagery with this application and provides an indication that the simulated data set captures all the necessary physical properties of icebergs for ship and iceberg classification

Sea-ice information and forecast needs for industry maritime stakeholders

Profound changes in Arctic sea-ice, a growing desire to utilize the Arctic’s abundant natural resources, and the potential competitiveness of Arctic shipping routes, all provide for increased industry marine activity throughout the Arctic Ocean. This is anticipated to result in further challenges for maritime safety. Those operating in ice-infested waters require various types of information for sea-ice and iceberg hazards. Ice information requirements depend on regional needs and whether the stakeholder wants to avoid ice all together, operate near or in the Marginal Ice Zone, or areas within the ice pack. An insight into user needs demonstrates how multiple spatial and temporal resolutions for sea-ice information and forecasts are necessary to provide information to the marine operating community for safety, planning, and situational awareness. Although ship-operators depend on sea-ice information for tactical navigation, stakeholders working in route and capacity planning can benefit from climatological and long-range forecast information at lower spatial and temporal resolutions where the interest is focused on open-water season. The advent of the Polar Code has brought with it additional information requirements, and exposed gaps in capacity and knowledge. Thus, future satellite data sources should be at resolutions that support both tactical and planning activities

Sea-ice information and forecast needs for industry maritime stakeholders

Profound changes in Arctic sea-ice, a growing desire to utilize the Arctic’s abundant natural resources, and the potential competitiveness of Arctic shipping routes, all provide for increased industry marine activity throughout the Arctic Ocean. This is anticipated to result in further challenges for maritime safety. Those operating in ice-infested waters require various types of information for sea-ice and iceberg hazards. Ice information requirements depend on regional needs and whether the stakeholder wants to avoid ice all together, operate near or in the Marginal Ice Zone, or areas within the ice pack. An insight into user needs demonstrates how multiple spatial and temporal resolutions for sea-ice information and forecasts are necessary to provide information to the marine operating community for safety, planning, and situational awareness. Although ship-operators depend on sea-ice information for tactical navigation, stakeholders working in route and capacity planning can benefit from climatological and long-range forecast information at lower spatial and temporal resolutions where the interest is focused on open-water season. The advent of the Polar Code has brought with it additional information requirements, and exposed gaps in capacity and knowledge. Thus, future satellite data sources should be at resolutions that support both tactical and planning activities

Iceberg topography and volume classification using TanDEM-X interferometry

Icebergs in polar regions affect water salinity, alter marine habitats, and impose serious hazards on maritime operations and navigation. These impacts mainly depend on the iceberg volume, which remains an elusive parameter to measure. We investigate the capability of TanDEM-X bistatic single-pass synthetic aperture radar interferometry (InSAR) to derive iceberg subaerial morphology and infer total volume. We cross-verify InSAR results with Operation IceBridge (OIB) data acquired near Wordie Bay, Antarctica, as part of the OIB/TanDEM-X Antarctic Science Campaign (OTASC). While icebergs are typically classified according to size based on length or maximum height, we develop a new volumetric classification approach for applications where iceberg volume is relevant. For icebergs with heights exceeding 5 m, we find iceberg volumes derived from TanDEM-X and OIB data match within 7 %. We also derive a range of possible iceberg keel depths relevant to grounding and potential impacts on subsea installations. These results suggest that TanDEM-X could pave the way for future single-pass interferometric systems for scientific and operational iceberg mapping and classification based on iceberg volume and keel depth

Использование спутниковой информации для обнаружения айсбергов и оценки айсберговой угрозы

Te methods of satellite monitoring of dangerous ice formations, namely icebergs in the Arctic seas, representing a threat to the safety of navigation and economic activity on the Arctic shelf are considered. Te main objective of the research is to develop methods for detecting icebergs using satellite radar data and high space resolution images in the visible spectral range. Te developed method of iceberg detection is based on statistical criteria for fnding gradient zones in the analysis of two-dimensional felds of satellite images. Te algorithms of the iceberg detection, the procedure of the false target identifcation, and determination the horizontal dimensions of the icebergs and their location are described. Examples of iceberg detection using satellite information with high space resolution obtained from Sentinel-1 and Landsat-8 satellites are given. To assess the iceberg threat, we propose to use a model of their drif, one of the input parameters of which is the size of the detected objects. Tree possible situations of observation of icebergs are identifed, namely, the «status» state of objects: icebergs on open water; icebergs in drifing ice; and icebergs in the fast ice. At the same time, in each of these situations, the iceberg can be grounded, that prevents its moving. Specifc features of the iceberg monitoring at various «status» states of them are considered. Te «status» state of the iceberg is also taken into account when assessing the degree of danger of the detected object. Te use of iceberg detection techniques based on satellite radar data and visible range images is illustrated by results of monitoring the coastal areas of the Severnaya Zemlya archipelago. Te approaches proposed to detect icebergs from satellite data allow improving the quality and efciency of service for a wide number of users with ensuring the efciency and safety of Arctic navigation and activities on the Arctic shelf.Рассмотрены методы обработки оперативной спутниковой информации об айсбергах. Применённая методика обнаружения айсбергов основана на статистических критериях поиска градиентных зон при анализе двухмерных полей. Использование методик обнаружения айсбергов по спутниковым радиолокационным данным и снимкам видимого диапазона показано на материалах мониторинга прибрежных районов Северной Земли

Arctic Sea Ice Characterization using Spaceborne Fully Polarimetric L-, C- and X-Band SAR with Validation by Airborne Measurements

Accepted manuscript version. Published version available at https://doi.org/10.1109/TGRS.2018.2809504.In recent years, spaceborne synthetic aperture radar (SAR) polarimetry has become a valuable tool for sea ice analysis. Here, we employ an automatic sea ice classification algorithm on two sets of spatially and temporally near coincident fully polarimetric acquisitions from the ALOS-2, Radarsat-2, and TerraSAR-X/TanDEM-X satellites. Overlapping coincident sea ice freeboard measurements from airborne laser scanner data are used to validate the classification results. The automated sea ice classification algorithm consists of two steps. In the first step, we perform a polarimetric feature extraction procedure. Next, the resulting feature vectors are ingested into a trained neural network classifier to arrive at a pixelwise supervised classification. Coherency matrix-based features that require an eigendecomposition are found to be either of low relevance or redundant to other covariance matrix-based features, which makes coherency matrix-based features dispensable for the purpose of sea ice classification. Among the most useful features for classification are matrix invariant-based features (geometric intensity, scattering diversity, and surface scattering fraction). Classification results show that 100% of the open water is separated from the surrounding sea ice and that the sea ice classes have at least 96.9% accuracy. This analysis reveals analogous results for both X-band and C-band frequencies and slightly different for the L-band. The subsequent classification produces similarly promising results for all four acquisitions. In particular, the overlapping image portions exhibit a reasonable congruence of detected sea ice when compared with high-resolution airborne measurements

Quad polarimetric synthetic aperture radar analysis of icebergs in Greenland and Svalbard

Polarimetric synthetic aperture radar (PolSAR) has been widely used in ocean and cryospheric applications. This is because, PolSAR can be used in all-day operations and in areas of cloud cover, and therefore can provide valuable large-scale monitoring in polar regions, which is very helpful to shipping and offshore maritime operations. In the last decades, attention has turned to the potential of PolSAR to detect icebergs in the Arctic since they are a major hazard to vessels. However, there is a substantial lack of literature exploring the potentialities of PolSAR and the understanding of iceberg scattering mechanisms. Additionally, it is not known if high resolution PolSAR can be used to detect icebergs smaller than 120 metres. This thesis aims to improve the knowledge of the use of PolSAR scattering mechanisms of icebergs, and detection of small icebergs. First, an introduction to PolSAR is outlined in chapter two, and monitoring of icebergs is presented in chapter three. The first data chapter (Chapter 4) is focused on developing a multi-scale analysis of icebergs using parameters from the Cloude-Pottier and the Yamaguchi decompositions, the polarimetric span and the Pauli scattering vector. This method is carried out using ALOS-2 PALSAR quad polarimetric L-band SAR on icebergs in Greenland. This approach outlines the good potential for using PolSAR for future iceberg classification. One of the main important outcomes is that icebergs are composed by a combination of single targets, which therefore may require a more complex way of processing SAR data to properly extract physical information. In chapter five, the problem of detecting icebergs is addressed by introducing six state-of-the-art detectors previously applied to vessel monitoring. These detectors are the Dual Intensity Polarisation Ratio Anomaly Detector (iDPolRAD), Polarimetric Notch Filter (PNF), Polarimetric Matched Filter (PMF), reflection symmetry (sym), Optimal Polarimetric Detector (OPD) and the Polarimetric Whitening Filter (PWF). Cloude-Pottier entropy, and first and third eigenvalues (eig1 and eig3) of the coherency matrix are also utilised as parameters for comparison. This approach uses the same ALOS-2 dataset, but also evaluates detection performance in two scenarios: icebergs in open ocean, and in sea ice. Polarimetric modes (quad-pol, dual-pol, and single intensities) are also considered for comparison. Currently it is very difficult to detect icebergs less than 120 metres in length using this approach, due to the scattering mechanisms of icebergs and sea ice being very similar. However, it was possible to obtain detection performances of the OPD and PWF, which both showed a Probability of Detection (PF) of 0.99 when the Probability of False Alarms (PF) was set to 10-5 in open ocean. Similarly, in dual pol images, the PWF gave the best performance with a PD of 0.90. Results in sea ice found eig3 to be the best detector with a PD of 0.90 while in dual-pol mode, iDPolRAD gave a PD of 0.978. Single intensity detector performance found the HV channel gave the best detection with a PD of 0.99 in open ocean and 0.87 in sea ice. In the previous two approaches, only satellite data is used. However, in chapter six, data from a ground-based Ku-band Gamma Portable Radio Interferometer (GPRI) instrument is introduced, providing images that are synchronised with the satellite acquisitions. In this approach, the same six detectors are applied to three multitemporal RADARSAT-2 quad pol C-band SAR images on icebergs in Kongsfjorden, Svalbard to evaluate the detection performance within a changing fjord environment. As before, we also make use of Cloude-Pottier entropy, eig1 and eig3. Finally, we evaluate the target-to-clutter ratio (TCR) of the icebergs and check for correlation between the backscattering coefficients and the iceberg dimension. The results obtained from this thesis present original additions to the literature that contributes to the understanding of PolSAR in cryospheric applications. Although these methods are applied to PolSAR and ground-based radar on vessels, they have been applied for the first time on icebergs in this thesis. To summarise, the main findings are that icebergs cannot be represented as single or partial targets, but they do exhibit a collection of single targets clustered together. This result leads to the fact that entropy is not sufficient as a parameter to detect icebergs. Detection results show that the OPD and PWF detectors perform best in an open ocean setting and using quad-pol mode. These results are degraded in dual-pol mode, while single intensity detection is best in the HV cross polarisation channel. When these detectors are applied to the RADARSAT-2 in Svalbard, the OPD and PWF detectors also perform best with PD values ranging between 0.5-0.75 for a PF of 0.01-0.05. However, the sea ice present in the fjord degrades performance across all detectors. Correlation plots with iceberg size show that a regression is not straightforward and Computer Vision methodologies may work best for this

Ocean remote sensing techniques and applications: a review (Part II)

As discussed in the first part of this review paper, Remote Sensing (RS) systems are great tools to study various oceanographic parameters. Part I of this study described different passive and active RS systems and six applications of RS in ocean studies, including Ocean Surface Wind (OSW), Ocean Surface Current (OSC), Ocean Wave Height (OWH), Sea Level (SL), Ocean Tide (OT), and Ship Detection (SD). In Part II, the remaining nine important applications of RS systems for ocean environments, including Iceberg, Sea Ice (SI), Sea Surface temperature (SST), Ocean Surface Salinity (OSS), Ocean Color (OC), Ocean Chlorophyll (OCh), Ocean Oil Spill (OOS), Underwater Ocean, and Fishery are comprehensively reviewed and discussed. For each application, the applicable RS systems, their advantages and disadvantages, various RS and Machine Learning (ML) techniques, and several case studies are discussed.Peer ReviewedPostprint (published version