Investigation of the microwave signatures of the Baltic Sea ice

It is essential for winter shipping in the Baltic Sea to get reliable and up-to-date information of its rapidly changing ice conditions. Spaceborne synthetic aperture radar (SAR) images are the only way to produce this information operationally in fine scale independent of daylight and nearly independent of weather conditions. Currently, classification algorithms for the RADARSAT-1 and ENVISAT SAR images utilize mainly the image structure and only limited information on sea ice geophysics and empirical statistics of backscattering signatures of various ice types are utilized. Therefore, interpretation of the classification results is often difficult. Both classification results and their interpretation should very likely improve with the addition of this information. Spaceborne microwave radiometer data are not suitable for the operational Baltic Sea ice monitoring aiding ship navigation due to their coarse spatial resolution, but they can provide an independent data source on the sea ice conditions for validation of the SAR classification algorithms. Both SAR and radiometer data based sea ice products can also be utilized in the geophysical studies of the Baltic Sea ice. In order to support development of operational classification algorithms for SAR and radiometer data, basic research on the microwave remote sensing of the Baltic Sea ice has been conducted in this work. The research work included the following topics: (1) statistics of C- and X-band backscattering signatures of various ice types, (2) statistics of L- and C-band polarimetric discriminants of various ice types, (3) radar incidence angle dependence of backscattering coefficient (σ°) in RADARSAT-1 SAR images, (4) dependence between standard deviation and measurement length for σ° signatures and its usability in sea ice classification, (5) comparison between SAR σ° time series and results from a thermodynamic snow/ice model, and (6) statistics of passive microwave signatures of various ice types. Additionally, a comprehensive literature review of the previous work on the microwave remote sensing of the Baltic Sea ice is presented. The main results of this work include the following. It is not possible to discriminate open water and various ice types using the level of σ°, co- or cross-polarization ratio, or standard deviation of σ°. C-band VH-polarized σ° at high incidence angle provides slightly better ice type discrimination accuracy than any other combination of C- and X-band radar parameters. VH-polarization is more suitable for estimating the degree of ice deformation than co-polarizations. Snow wetness has a large effect on the σ° statistics. Notably, when snow cover is wet then the σ° contrasts between various ice types are smaller than in the dry snow case. Incidence angle dependence of the C-band HH-polarized σ° was derived for level ice and deformed ice. It is utilized in the operational SAR classification algorithms developed by Finnish Institute of Marine Research. The method for deriving the σ° incidence angle dependence is applicable for any SAR sensor. There is a large variation of level ice σ° with changing weather conditions. A 1-D high-resolution thermodynamic snow/ice model generally helps to interpret changes in the σ° time series. The modeled snow and ice surface temperature, cases of snow melting, and evolution of snow and ice thickness are related to the changes in σ°. It was found out that the standard deviation of σ° for various ice types depends on the length of measurement. This may be utilized in the SAR image classification. It is not possible to resolve concentrations of thin new ice and all other ice types combined in the Baltic Sea using radiometer data as has been done for the Arctic seasonal ice zones.Talvimerenkulku Itämerellä tarvitsee luotettavaa ja ajantasaista informaatiota Itämeren nopeasti muuttuvista jääoloista. Synteettisen apertuurin tutkan (SAR) kuvat ovat ainoa tapa tuottaa operatiivisesti tarvittavaa jääinformaatiota riippumatta päivänvalon määrästä ja lähes riippumatta sääolosuhteista. RADARSAT-1 ja ENVISAT SAR-tutkakuvien luokittelualgoritmit perustuvat tällä hetkellä lähinnä kuvien rakenteeseen, eikä merijään geofysiikkaa ja empiiristä tilastotietoa eri jäätyyppien sirontavasteista hyödynnetä kuin rajallisesti. SAR-kuvien luokittelutulosten tulkitseminen on siten usein vaikeaa. Sekä itse luokittelutulokset, että niiden tulkinta parantuisivat, jos luokittelualgorimit hyödyntäisivät edellä mainittua tietoa. Satelliittiradiometrien kuvat eivät sovellu Itämeren jään operatiiviseen monitorointiin niiden karkean spatiaalisen resoluution vuoksi. Niillä kuitenkin voitaisiin validoida SAR-kuvien luokittelualgoritmeja, koska ne ovat SAR-kuvista riippumaton datalähde Itämeren jääoloista. Tässä työssä on suoritettu seuraavaa perustutkimusta Itämeren jään mikroaaltokaukokartoituksessa, minkä tarkoituksena on tukea SAR- ja radiometrikuvien operatiivisten luokittelualgoritmien kehitystyötä: (1) eri jäätyyppien C- ja X-kanavien sirontakertoimien (σ°) statistiikka, (2) eri jäätyyppien L- ja C-kanavien polarimetristen diskriminanttien statistiikka, (3) σ°:n mittauskulmariippuvuus RADARSAT-1 SAR-kuvissa, (4) σ°:n keskihajonnan ja mittausmatkan välinen riippuvuus ja hyödyntäminen jäätyyppiluokittelussa, (5) SAR-kuvien sirontakerroinaikasarjojen vertailu merijään termodynamiikkamalliin, ja (6) eri jäätyyppien kirkkauslämpötilojen statistiikka. Työssä saavutettiin seuraavia merkittäviä tuloksia. Eri jäätyyppien ja avoveden luokittelu ei ole mahdollista käyttäen sirontakerrointa, yhdensuuntais- ja ristipolarisaatiosuhdetta tai σ° keskihajontaa. C-kanavan VH-polarisaation σ° suurella mittauskulmalla luokittelee eri jäätyypit hieman paremmin kuin mikään muu C- ja X-kanavan tutkaparametrikombinaatio. Merijään deformoitumisasteen estimointiin sopii paremmin VH-polarisaation σ° kuin yhdensuuntaispolarisaation. Lumipeitteen kosteudella on suuri vaikutus sirontakerroinstatistiikkaan; erityisesti, kun lumipeite on märkä on sirontakerroinkontrasti eri jäätyyppien välillä pienempi kun lumipeite on kuiva. C-kanavan HH-polarisaation σ°:n mittauskulmariippuvuus määritettiin tasaiselle ja deformoituneelle jäälle. Mittauskulmariippuvuuden laskentamenetelmää voidaan käyttää mille tahansa SAR-tutkakuvalle. Muuttuvat sääolosuhteet aiheuttavat suuria muutoksia tasaisen jään σ°:ssa. Merijään termodynamiikkamalli yleensä auttaa selittämään muutoksia σ°:n aikasarjassa. σ°:n muutokset ovat yhteydessä termodynamiikkamallilla laskettuihin lumen ja jään parametreihin. σ°:n keskihajonnan havaittiin riippuvan etäisyydestä. Tätä riippuvuutta voitaneen hyödyntään SAR-kuvien luokittelussa. Itämerellä satelliittiradiometridatalla pystytään määrittämään vain merijään kokonaiskonsetraatio, toisin kuin arktisten merien kausiluontoisilla merijääalueilla, missä myös eri jäätyyppien konsentraatioiden määrittäminen on mahdollista.reviewe

TanDEM-X multiparametric data features in sea ice classification over the Baltic sea

In this study, we assess the potential of X-band Interferometric Synthetic Aperture Radar imagery for automated classification of sea ice over the Baltic Sea. A bistatic SAR scene acquired by the TanDEM-X mission over the Bothnian Bay in March of 2012 was used in the analysis. Backscatter intensity, interferometric coherence magnitude, and interferometric phase have been used as informative features in several classification experiments. Various combinations of classification features were evaluated using Maximum likelihood (ML), Random Forests (RF) and Support Vector Machine (SVM) classifiers to achieve the best possible discrimination between open water and several sea ice types (undeformed ice, ridged ice, moderately deformed ice, brash ice, thick level ice, and new ice). Adding interferometric phase and coherence-magnitude to backscatter-intensity resulted in improved overall classification performance compared to using only backscatter-intensity. The RF algorithm appeared to be slightly superior to SVM and ML due to higher overall accuracies, however, at the expense of somewhat longer processing time. The best overall accuracy (OA) for three methodologies were achieved using combination of all tested features were 71.56, 72.93, and 72.91% for ML, RF and SVM classifiers, respectively. Compared to OAs of 62.28, 66.51, and 63.05% using only backscatter intensity, this indicates strong benefit of SAR interferometry in discriminating different types of sea ice. In contrast to several earlier studies, we were particularly able to successfully discriminate open water and new ice classes.Peer reviewe

Detection and classification of sea ice from spaceborne multi-frequency synthetic aperture radar imagery and radar altimetry

The sea ice cover in the Arctic is undergoing drastic changes. Since the start of satellite observations by microwave remote sensing in the late 1970\u27s, the maximum summer sea ice extent has been decreasing and thereby causing a generally thinner and younger sea ice cover. Spaceborne radar remote sensing facilitates the determination of sea ice properties in a changing climate with the high spatio-temporal resolution necessary for a better understanding of the ongoing processes as well as safe navigation and operation in ice infested waters.The work presented in this thesis focuses on the one hand on synergies of multi-frequency spaceborne synthetic aperture radar (SAR) imagery for sea ice classification. On the other hand, the fusion of radar altimetry observations with near-coincidental SAR imagery is investigated for its potential to improve 3-dimensional sea ice information retrieval.Investigations of ice/water classification of C- and L-band SAR imagery with a feed-forward neural network demonstrated the capabilities of both frequencies to outline the sea ice edge with good accuracy. Classification results also indicate that a combination of both frequencies can improve the identification of thin ice areas within the ice pack compared to C-band alone. Incidence angle normalisation has proven to increase class separability of different ice types. Analysis of incidence angle dependence between 19-47\ub0 at co- and cross-polarisation from Sentinel-1 C-band images closed a gap in existing slope estimates at cross-polarisation for multiyear sea ice and confirms values obtained in other regions of the Arctic or with different sensors. Furthermore, it demonstrated that insufficient noise correction of the first subswath at cross-polarisation increased the slope estimates by 0.01 dB/1\ub0 for multiyear ice. The incidence angle dependence of the Sentinel-1 noise floor affected smoother first-year sea ice and made the first subswath unusable for reliable incidence angle estimates in those cases.Radar altimetry can complete the 2-dimensional sea ice picture with thickness information. By comparison of SAR imagery with altimeter waveforms from CryoSat-2, it is demonstrated that waveforms respond well to changes of the sea ice surface in the order of a few hundred metres to a few kilometres. Freeboard estimates do however not always correspond to these changes especially when mixtures of different ice types are found within the footprint. Homogeneous ice floes of about 10 km are necessary for robust averaged freeboard estimates. The results demonstrate that multi-frequency and multi-sensor approaches open up for future improvements of sea ice retrievals from radar remote sensing techniques, but access to in-situ data for training and validation will be critical

Compaction of C-band synthetic aperture radar based sea ice information for navigation in the Baltic Sea

In this work operational sea ice synthetic aperture radar (SAR) data products were improved and developed. A SAR instrument is transmitting electromagnetic radiation at certain wavelengths and measures the radiation which is scattered back towards the instrument from the target, in our case sea and sea ice. The measured backscattering is converted to an image describing the target area through complex signal processing. The images, however, differ from optical images, i.e. photographs, and their visual interpretation is not straightforward. The main idea in this work has been to deliver the essential SAR-based sea ice information to end-users (typically on ships) in a compact and user-friendly format. The operational systems at Finnish Institute of Marine Research (FIMR) are currently based on the data received from a Canadian SAR-satellite, Radarsat-1. The operational sea ice classification, developed by the author with colleagues, has been further developed. One problem with the SAR data is typically that the backscattering varies depending on the incidence angle. The incidence angle is the angle in which the transmitted electromagnetic wave meets the target surface and it varies within each SAR image and between different SAR images depending on the measuring geometry. To improve this situation, an incidence angle correction algorithm to normalize the backscattering over the SAR incidence angle range for Baltic Sea ice has been developed as part of this work. The algorithm is based on SAR backscattering statistics over the Baltic Sea. To locate different sea ice areas in SAR images, a SAR segmentation algorithm based on pulse-coupled neural networks has been developed and tested. The parameters have been tuned suitable for the operational data in use at FIMR. The sea ice classification is based on this segmentation and the classification is segment-wise rather than pixel-wise. To improve SAR-based distinguishing between sea ice and open water an open water detection algorithm based on segmentation and local autocorrelation has been developed. Also ice type classification based on higher-order statistics and independent component analysis have been studied to get an improved SAR-based ice type classification. A compression algorithm for compressing sea ice SAR data for visual use has been developed. This algorithm is based on the wavelet decomposition, zero-tree structure and arithmetic coding. Also some properties of the human visual system were utilized. This algorithm was developed to produce smaller compressed SAR images, with a reasonable visual quality. The transmission of the compressed images to ships with low-speed data connections in reasonable time is then possible. One of the navigationally most important sea ice parameters is the ice thickness. SAR-based ice thickness estimation has been developed and evaluated as part of this work. This ice thickness estimation method uses the ice thickness history derived from digitized ice charts, made daily at the Finnish Ice Service, as its input, and updates this chart based on the novel SAR data. The result is an ice thickness chart representing the ice situation at the SAR acquisition time in higher resolution than in the manually made ice thickness charts. For the evaluation of the results a helicopter-borne ice thickness measuring instrument, based on electromagnetic induction and laser altimeter, was used.reviewe

Sea ice local surface topography from single-pass satellite InSAR measurements: a feasibility study

Quantitative parameters characterizing the sea ice surface topography are needed in geophysical investigations such as studies on atmosphere–ice interactions or sea ice mechanics.Recently, the use of space-borne single-pass interferometric synthetic aperture radar (InSAR) for retrieving the ice surface topography has attracted notice among geophysicists. In this paper the potential of InSAR measurements is examined for several satellite configurations and radar frequencies, considering statistics of heights and widths of ice ridges as well as possible magnitudes of ice drift. It is shown that, theoretically, surface height variations can be retrievedwith relative errors < 0.5 m. In practice, however, the sea ice drift and open water leads may contribute significantly to the measured interferometric phase. Another essential factor is the dependence of the achievable interferometric baseline on the satellite orbit configurations. Possibilities to assess the influence of different factors on the measurement accuracy are demonstrated: signal-to-noise ratio, presence of a snow layer, and the penetration depth into the ice. Practical examples of sea surface height retrievals from bistatic SAR images collectedduring the TanDEM-X Science Phase are presented

Sea ice classification of TerraSAR-X ScanSAR images for the MOSAiC expedition incorporating per-class incidence angle dependency of image texture

We provide sea ice classification maps of a subweekly time series of single (horizontal–horizontal, HH) polarization X-band TerraSAR-X scanning synthetic aperture radar (TSX SC) images from November 2019 to March 2020, covering the Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAiC) expedition. This classified time series benefits from the wide spatial coverage and relatively high spatial resolution of TSX SC data and is a useful basic dataset for future MOSAiC studies on physical sea ice processes and ocean and climate modeling. Sea ice is classified into leads, young ice with different backscatter intensities, and first-year ice (FYI) or multiyear ice (MYI) with different degrees of deformation. We establish the per-class incidence angle (IA) dependencies of TSX SC intensities and gray-level co-occurrence matrix (GLCM) textures and use a classifier that corrects for the class-specific decreasing backscatter with increasing IAs, with both HH intensities and textures as input features. Optimal parameters for texture calculation are derived to achieve good class separation while maintaining maximum spatial detail and minimizing textural collinearity. Class probabilities yielded by the classifier are adjusted by Markov random field contextual smoothing to produce classification results. The texture-based classification process yields an average overall accuracy of 83.70 % and good correspondence to geometric ice surface roughness derived from in situ ice thickness measurements (correspondence consistently close to or higher than 80 %). A positive logarithmic relationship is found between geometric ice surface roughness and TSX SC HH backscatter intensity, similar to previous C- and L-band studies. Areal fractions of classes representing ice openings (leads and young ice) show prominent increases in middle to late November 2019 and March 2020, corresponding well to ice-opening time series derived from in situ data in this study and those derived from satellite synthetic aperture radar (SAR) and optical data in other MOSAiC studies

Kara and Barents sea ice thickness estimation based on CryoSat-2 radar altimeter and Sentinel-1 dual-polarized synthetic aperture radar

We present a method to combine CryoSat-2 (CS2) radar altimeter and Sentinel-1 synthetic aperture radar (SAR) data to obtain sea ice thickness (SIT) estimates for the Barents and Kara seas. From the viewpoint of tactical navigation, along-track altimeter SIT estimates are sparse, and the goal of our study is to develop a method to interpolate altimeter SIT measurements between CS2 ground tracks. The SIT estimation method developed here is based on the interpolation of CS2 SIT utilizing SAR segmentation and segmentwise SAR texture features. The SIT results are compared to SIT data derived from the AARI ice charts; to ORAS5, PIOMAS and TOPAZ4 ocean-sea ice data assimilation system reanalyses; to combined CS2 and Soil Moisture and Ocean Salinity (SMOS) radiometer weekly SIT (CS2SMOS SIT) charts; and to the daily MODIS (Moderate Resolution Imaging Spectro-radiometer) SIT chart. We studied two approaches: CS2 directly interpolated to SAR segments and CS2 SIT interpolated to SAR segments with mapping of the CS2 SIT distributions to correspond to SIT distribution of the PIOMAS ice model. Our approaches yield larger spatial coverage and better accuracy compared to SIT estimates based on either CS2 or SAR data alone. The agreement with modelled SIT is better than with the CS2SMOS SIT. The average differences when compared to ice models and the AARI ice chart SIT were typically tens of centimetres, and there was a significant positive bias when compared to the AARI SIT (on average 27 cm) and a similar bias (24 cm) when compared to the CS2SMOS SIT. Our results are directly applicable to the future CRISTAL mission and Copernicus programme SAR missions.Peer reviewe

Contrasting seasonality in optical-biogeochemical properties of the Baltic Sea

Optical-biogeochemical relationships of particulate and dissolved organic matter are presented in support of remote sensing of the Baltic Sea pelagic. This system exhibits strong seasonality in phytoplankton community composition and wide gradients of chromophoric dissolved organic matter (CDOM), properties which are poorly handled by existing remote sensing algorithms. Absorption and scattering properties of particulate matter reflected the seasonality in biological (phytoplankton succession) and physical (thermal stratification) processes. Inherent optical properties showed much wider variability when normalized to the chlorophyll-a concentration compared to normalization to either total suspended matter dry weight or particulate organic carbon. The particle population had the largest optical variability in summer and was dominated by organic matter in both seasons. The geographic variability of CDOM and relationships with dissolved organic carbon (DOC) are also presented. CDOM dominated light absorption at blue wavelengths, contributing 81% (median) of the absorption by all water constituents at 400 nm and 63% at 442 nm. Consequentially, 90% of water-leaving radiance at 412 nm originated from a layer (z90) no deeper than approximately 1.0 m. With water increasingly attenuating light at longer wavelengths, a green peak in light penetration and reflectance is always present in these waters, with z90 up to 3.0–3.5 m depth, whereas z90 only exceeds 5 m at biomass < 5 mg Chla m-3. High absorption combined with a weakly scattering particle population (despite median phytoplankton biomass of 14.1 and 4.3 mg Chla m-3 in spring and summer samples, respectively), characterize this sea as a dark water body for which dedicated or exceptionally robust remote sensing techniques are required. Seasonal and regional optical-biogeochemical models, data distributions, and an extensive set of simulated remote-sensing reflectance spectra for testing of remote sensing algorithms are provided as supplementary data