Measurements of Sea Surface Currents in the Baltic Sea Region Using Spaceborne Along-Track InSAR

The main challenging problems in ocean current retrieval from along-track interferometric (ATI)-synthetic aperture radar (SAR) are phase calibration and wave bias removal. In this paper, a method based on differential InSAR (DInSAR) technique for correcting the phase offset and its variation is proposed. The wave bias removal is assessed using two different Doppler models and two different wind sources. In addition to the wind provided by an atmospheric model, the wind speed used for wave correction in this work is extracted from the calibrated SAR backscatter. This demonstrates that current retrieval from ATI-SAR can be completed independently of atmospheric models. The retrieved currents, from four TanDEM-X (TDX) acquisitions over the 6resund channel in the Baltic Sea, are compared to a regional ocean circulation model. It is shown that by applying the proposed phase correction and wave bias removal, a good agreement in spatial variation and current direction is achieved. The residual bias, between the ocean model and the current retrievals, varies between 0.013 and 0.3 m/s depending on the Doppler model and wind source used for wave correction. This paper shows that using SAR as a source of wind speed reduces the bias and root-mean-squared-error (RMSE) of the retrieved currents by 20% and 15%, respectively. Finally, the sensitivity of the sea current retrieval to Doppler model and wind errors are discussed

Retrieval of Ocean Surface Currents and Winds Using Satellite SAR backscatter and Doppler frequency shift

Ocean surface winds and currents play an important role for weather, climate, marine life, ship navigation, oil spill drift and search and rescue. In-situ observations of the ocean are sparse and costly. Satellites provide a useful complement to these observations. Synthetic aperture radar (SAR) is particularly attractive due to its high spatial resolution and its capability to extract both sea surface winds and currents day and night and almost independent of weather.The work in this thesis involves processing of along-track interferometric SAR (ATI-SAR) data, analysis of the backscatter and Doppler frequency shift, and development of wind and current retrieval algorithms. Analysis of the Doppler frequency shift showed a systematic bias. A calibration method was proposed and implemented to correct for this bias. Doppler analysis also showed that the wave contribution to the SAR Doppler centroid often dominates over the current contribution. This wave contribution is estimated using existing theoretical and empirical Doppler models. For wind and current retrieval, two methods were developed and implemented.The first method, called the direct method, consists of retrieval of the wind speed from SAR backscatter using an empirical backscatter model. In order to retrieve the radial current, the retrieved wind speed is used to correct for the wave contribution. The current retrieval was assessed using two different (theoretical and empirical) Doppler models and wind inputs (model and SAR-derived). It was found that the results obtained by combining the Doppler empirical model with the SAR-derived wind speed were more consistent with ocean models.The second method, called Bayesian method, consists of blending the SAR observables (backscatter and Doppler shift) with an atmospheric and an oceanic model to retrieve the total wind and current vector fields. It was shown that this method yields more accurate estimates, i.e. reduces the models biases against in-situ measurements. Moreover, the method introduces small scale features, e.g. fronts and meandering, which are weakly resolved by the models.The correlation between the surface wind vectors and the SAR Doppler shift was demonstrated empirically using the Doppler shift estimated from over 300 TanDEM-X interferograms and ECMWF reanalysis wind vectors. Analysis of polarimetric data showed that theoretical models such as Bragg and composite surface models over-estimate the backscatter polarization ratio and Doppler shift polarization difference. A combination of a theoretical Doppler model and an empirical modulation transfer function was proposed. It was found that this model is more consistent with the analyzed data than the pure theoretical models.The results of this thesis will be useful for integrating SAR retrievals in ocean current products and assimilating SAR observables in the atmospheric, oceanic or coupled models. The results are also relevant for preparation studies of future satellite missions

Sea Surface Current Measurements Using Along-Track Interferometric SAR

Ocean currents affect the weather, the climate and the marine ecosystem. Observing ocean currents is important for understanding the upper-ocean layer dynamics and its interaction with the other components of the climate system. In-situ measurements are sparse and their deployment and maintenance is costly. Satellite remote sensing with large spatial coverage offers a good complement to the in-situ observations. In this work we have studied the spaceborne Along-Track Interferometric SAR (ATI-SAR) for measuring sea surface currents. The measurement principle is based on the fact that the phase difference between two SAR acquisitions is directly related to radial (line-of-sight) velocity of the illuminated surface. Previous studies based on similar systems were carried out in areas with well defined and strong tidal currents ( ~1 - 3 m/s). In this work we demonstrate thecapability of ATI-SAR, through several study cases, in areas with weak currents ( <0.5 m /s). This is challenging for the satellite measurements of surface currents because it requires very accurate processing and retrieval algorithms. In addition, it has been found that wave motion contribution, systematically dominates the measured ATI-SAR radial velocity in these weak current areas. Estimation of the wave motion contribution relies on high-resolution and accurate wind data. Thus, a wind speed retrieval algorithm from SAR is needed to support the ATI-SAR current retrieval. We have shown that with an appropriate processing of the ATI-SAR phase and with applying the necessary corrections to the measured velocity a good agreement with ocean circulation models is achieved (rmse =0.1 m /s). These corrections include phase calibration and wind compensation to correct for instrument and geophysical systematic errors, respectively. Finally, a novel method for removing the wind direction ambiguity, based on the ATI-SAR phase, is presented. In previous methods, the wind ambiguity removal was based on external information, e.g. an atmospheric model or on visual observation of wind shadows

Ocean Surface Observations Using the TanDEM-X Satellite Formation

The TanDEM-X SAR satellite formation permits improved ocean surface observations by means of bistatic along-track interferometry (ATI) when compared to single-satellite systems. The flexible imaging geometry of its two cooperating SAR sensors forms an interferometer that can achieve very high sensitivity to motions of objects on ground. This way, radar imaging of surface currents with unprecedented accuracy, high spatial resolution and wide coverage at the same time becomes possible. We demonstrate the capabilities of the sensor in the contexts of tidal current mapping, measurement of thermohaline and wind-driven ocean currents as well as detection of areas with surface films. We have developed a dedicated postprocessing system for TanDEM-X image products that allows extracting surface current information from the data. By this paper, we address bistatic data acquisition and processing aspects for sea surface imaging with TanDEM-X like interferometric baseline geometry, temporal decorrelation, and phase calibration. We present a variety of examples of data evaluation that clearly demonstrate the application potential of the methodology

Empirical Relationship Between the Doppler Centroid Derived From X-Band Spaceborne InSAR Data and Wind Vectors

One of the challenges in ocean surface current retrieval from synthetic aperture radar (SAR) data is the estimation and removal of the wave-induced Doppler centroid (DC). This article demonstrates empirically the relationship between the dc derived from spaceborne X-band InSAR data and the ocean surface wind and waves. In this study, we analyzed over 300 TanDEM-X image pairs. It is found that the general characteristics of the estimated dc follow the theoretically expected variation with incidence angle, wind speed, and wind direction. An empirical geophysical model function (GMF) is fit to the estimated dc and compared to existing models and previous experiments. Our GMF is in good agreement (within 0.2 m/s) with other models and data sets. It is found that the wind-induced Doppler velocity contributes to the total Doppler velocity with about 15% of the radial wind speed. This is much larger than the sum of the contributions from the Bragg waves (~0.2 m/s) and the wind-induced drift current (~3% of wind speed). This indicates a significant (dominant) contribution of the long wind waves to the SAR dc. Moreover, analysis of dual-polarized data shows that the backscatter polarization ratio (PR=σ⁰VV/σ⁰HH) and the dc polarization difference (PD=|dcVV|-|dcHH|) are systematically larger than 1 and smaller than 0 Hz, respectively, and both increase in magnitude with incidence angle. The estimated PR and PD are compared to other theoretical and empirical models. The Bragg scattering theory-based (pure Bragg and composite surface) models overestimate both PR and PD, suggesting that other scattering mechanisms, e.g., wave breaking, are involved. In general, it is found that empirical models are more consistent with both backscatter and Doppler data than theory-based models. This motivates a further improvement of SAR dc GMFs

X-Band Interferometric SAR Observations of Baltic Fast Ice

Detailed mapping of fast-ice deformation can be used to characterize the rheological behavior of fast ice and subsequently improve sea ice modeling. This study uses interferometric synthetic aperture radar to map fast-ice deformation with unprecedented spatial resolution (meter range) and sensitivity (cm-mm range). Two interferometric acquisitions, each with a temporal baseline of 24 h, were performed by the X-band SAR satellite constellation Cosmo-SkyMed over the northeast Bay of Bothnia in the middle of the 2012 ice season. The first interferogram shows deformation of the fast ice due to force from impinging drift ice, and the normal strain within the fast ice is measured. Complementary intensity correlation measurements reveal a slow movement of the drift ice toward the fast ice. The second interferogram exhibits a low fringe rate over the fast ice with fringes being aligned along the coastline. Deformation appears to be stronger around leads, skerries, and grounded ice ridges. It is also observed that the coherence images provide information that is complementary to the information in the backscatter images

Investigation of Sea Ice Using Multiple Synthetic Aperture Radar Acquisitions

The papers of this thesis are not available in Munin. Paper I: Yitayew, T. G., Ferro-Famil, L., Eltoft, T. & Tebaldini, S. (2017). Tomographic imaging of fjord ice using a very high resolution ground-based SAR system. Available in IEEE Transactions on Geoscience and Remote Sensing, 55 (2):698-714. Paper II: Yitayew, T. G., Ferro-Famil, L., Eltoft, T. & Tebaldini, S. (2017). Lake and fjord ice imaging using a multifrequency ground-based tomographic SAR system. Available in IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing, 10(10):4457-4468. Paper III: Yitayew, T. G., Divine, D. V., Dierking, W., Eltoft, T., Ferro-Famil, L., Rosel, A. & Negrel, J. Validation of Sea ice Topographic Heights Derived from TanDEMX Interferometric SAR Data with Results from Laser Profiler and Photogrammetry. (Manuscript).The thesis investigates imaging in the vertical direction of different types of ice in the arctic using synthetic aperture radar (SAR) tomography and SAR interferometry. In the first part, the magnitude and the positions of the dominant scattering contributions within snow covered fjord and lake ice layers are effectively identified by using a very high resolution ground-based tomographic SAR system. Datasets collected at multiple frequencies and polarizations over two test sites in Tromsø area, northern Norway, are used for characterizing the three-dimensional response of snow and ice. The presented experimental results helped to improve our understanding of the interaction between radar waves and snow and ice layers. The reconstructed radar responses are also used for estimating the refractive indices and the vertical positions of the different sub-layers of snow and ice. The second part of the thesis deals with the retrieval of the surface topography of multi-year sea ice using SAR interferometry. Satellite acquisitions from TanDEM-X over the Svalbard area are used for analysis. The retrieved surface height is validated by using overlapping helicopter-based stereo camera and laser profiler measurements, and a very good agreement has been found. The work contributes to an improved understanding regarding the potential of SAR tomography for imaging the vertical scattering distribution of snow and ice layers, and for studying the influence of both sensor parameters such as its frequency and polarization and scene properties such as layer stratification, air bubbles and small-scale roughness of the interfaces on snow and ice backscattered signal. Moreover, the presented results reveal the potential of SAR interferometry for retrieving the surface topography of sea ice

Direct comparison of sea surface velocity estimated from Sentinel-1 and TanDEM-X SAR data

This paper presents the first direct comparison of the sea surface radial velocity (RVL) derived from the two satellite SAR systems Sentinel-1 and TanDEM-X, operating at different frequencies and imaging modes. The RVL is derived from the Doppler centroid (Dc) provided in the Sentinel-1 OCN product and from the along-track interferometric phase of the TanDEM-X. The comparison is carried out using opportunistic acquisitions, collocated in space and time, over three different sites. First, it is observed that the RVL derived from both satellites is biased, thus calibration is applied using the land as a reference. The comparison shows that the correlation and the mean RVL bias between the two datasets depend on the differences in acquisition time, incidence angle and azimuth angle, and on wind and current speed and direction. It is found that, given a time difference of < 20 min, the spatial correlation coefficient is relatively high (between 0.7 and 0.93), which indicates that the two SAR systems observe similar sea surface current fields. The spatial correlation degrades primarily due to increasing time difference and decreasing current magnitudes. The mean RVL bias increases primarily with the radial wind speed, which suggests that the RVL bias is mainly due to the wave-induced Doppler shift. This study shows that under certain conditions, i.e. similar acquisition geometry and short time delay, a good agreement between the two independently derived RVL is achieved. This encourages a synergistic use of the sea surface velocity estimated from different C- and X-band SAR systems

Remote sensing of sea ice properties and dynamics using SAR interferometry

Landfast ice is attached to the coastline and islands and stays immobile over most of the ice season. It is an important element of polar ecosystems and plays a vital role as a marine habitat and in life of local people and economy through offshore technology. Landfast ice is routinely used for on-ice traffic, tourism, and industry, and it protects coasts from storms in winter from erosion. However, landfast ice can break or experience deformation in order of centimeters to meters, which can be dangerous for the coastline and man-made structures, beacons, on-ice traffic, and represents a safety risk for working on the ice and local people. Therefore, landfast ice deformation and stability are important topics in coastal engineering and sea ice modeling. In the framework of this dissertation, InSAR (SAR Interferometry) technology has been applied for deriving landfast ice displacements (publication I), and mapping sea ice morphology, topography and its temporal change (publication III). Also, advantages of InSAR remote sensing in sea ice classification compared to backscatter intensity were demonstrated (publications II and IV). In publication I, for the first time, Sentinel-1 repeat-pass InSAR data acquired over the landfast ice areas were used to study the landfast ice displacements in the Gulf of Bothnia. An InSAR pair with a temporal baseline of 12 days acquired in February 2015 was used. In the study, the surface of landfast ice was stable enough to preserve coherence over the 12-day period, enabling analysis of the interferogram. The advantage of this long temporal baseline is in separating the landfast ice from drift ice and detecting long-term trends in deformation maps. The interferogram showed displacements of landfast ice on the order of 40 cm. The main factor seemed to be compression by drift ice, which was driven against the landfast ice boundary by strong winds from southwest. Landfast ice ridges can hinder ship navigation, but grounded ridges help to stabilize the ice cover. In publication III, ridge formation and displacements in the landfast ice near Utqiaġvik, Alaska were examined. The phase signatures of two single-pass bistatic X-band SAR (Synthetic Aperture Radar) image pairs acquired by TanDEM-X (TerraSAR-X add-on for Digital Elevation Measurements) satellite on 13 and 24 January 2012 were analyzed. Altogether six cases were identified with ridge displacement in four and formation in two cases under onshore compression. The ridges moved approximately 0.6 and 3.7 km over the study area and ridge formation reached up to 1 meter in upward. The results well corresponded with the locations identified as convergence zones retrieved from the drift algorithm generated by a SAR-based sea ice-tracking algorithm, backscatter intensity images and coastal radar imagery. This method could potentially be used in future to evaluate sea ice stability and ridge formation. A bistatic InSAR pair acquired by the TanDEM-X mission in March 2012 over the Bothnian Bay was used in two further studies (publications II and IV). The potential of X-band InSAR imagery for automated sea ice classification was evaluated. The first results were presented in publication II and the data were further elaborated in publication IV. The backscatter intensity, coherence magnitude and InSAR-phase features, as well as their different combinations, were used as the informative features in classification experiments. In publication II, the purpose was to assess ice properties on the scale used in ice charting, with ice types based on ice concentration and sea ice morphology, while in publication IV, a detailed small-scale analysis was performed. In addition, the sampling design was different in these publications. In publication II, to achieve the best discrimination between open water and several sea-ice types, RF (Random Forests) and ML (Maximum likelihood) classifiers were employed. The best overall accuracies were achieved by combining backscatter intensity & InSAR-phase using RF approach and backscatter intensity & coherence-magnitude using ML approach. The results showed the advantage of adding InSAR features to backscatter intensity for sea ice classification. In the further study (publication IV), a set of state-of-the-art classification approaches including ML, RF and SVM (Support Vector Machine) classifiers were used to achieve the best discrimination between open water and several sea-ice types. Adding InSAR-phase and coherence magnitude to backscatter intensity improved the OA (Overall Accuracy) compared to using only backscatter intensity. The RF and SVM algorithms gave somewhat larger OA compared to ML at the expense of a somewhat longer processing time. Results of publications II and IV demonstrate InSAR features have potential to improve sea ice classification. InSAR could be used by operational ice services to improve mapping accuracy of automated sea ice charting with statistical and machine learning classification approaches.Viime vuosikymmeninä satelliittivälitteisestä SAR-tutkasta on tullut erittäin tärkeä työkalu merijään kaukokartoituksessa. Tämän tutka perustuu sähkömagneettisten aaltojen sirontaan kiinnostavasta kohteesta takaisin tutkaan, mitä seuraa signaalin voimakkuuden mittaaminen. SAR-tutkat käyttävät synteettistä antennia, joka perustuu satelliitin liikkeeseen, mikä mahdollistaa tarkkojen, korkean erotuskyvyn kuvien tuottamisen. SAR-anturit mittaavat myös signaalin vaihetta, jota käytetään interferometria tekniikassa pinnan topografian ja siirtymien laskemiseen eri sovelluksissa, kuten maan muodonmuutoksissa, tarkassa kartoituksessa, maanjäristyksen arvioinnissa ja tulivuorenpurkauksien tarkkailussa. Interferometri tekniikkaa käytettiin tässä opinnäytetyössä pienten jäänsiirtymien analysointiin kiintojäävyöhykkeellä, joka on kiinni rantaviivassa ja saarissa eikä juuri liiku tuulien tai virtausten mukana. Kiintojääalueilla on pohjaan tarttuneita jäävalleja, jotka edistävät kiintojääpeitteen vakautumista. Kiintojäällä on tärkeä rooli merellisenä elinympäristönä, maankäytön kysymyksissä sekä paikallisten ihmisten elämässä ja meritekniikassa. Kiintojää voi murtua liikahdella useita metrejä, mikä voi olla vaarallista rakenteille, majakoille ja jäällä liikkujille. Tässä väitöskirjassa Sentinel-1A ja TanDEM-X satelliitteja ja interferometri tekniikkaa on käytetty arktisilla alueilla ja Itämerellä mittaamaan kiintojään muodonmuutoksia ja siirtymiä sekä niihin liittyviä mekanismeja. Lisäksi on tutkittu automaattista merijääluokitusta interferometrian apuohjelmiston avulla, mikä laajentaa operatiivisten merijääpalvelujen tutkahavaintojen käyttöä. Sentinel-1A:n avulla voitiin tarkastella 12 päivän pituisia muutoksia Pohjanlahden kiintojäävyöhykkeellä, kun interferometria tekniikka mittasi voimakkaan tuulen aiheuttaman 40 cm:n siirtymiä. Pohjoisella jäämerellä voitiin tunnistaa jäävallien siirtymiä ja muodostumia. Vallit siirtyivät noin 0,6 ja 3,7 km matkoja ja muodostuessaan ne kasvoivat metrin korkeuteen. Interferometri tekniikan lisääminen tutkakuvauksen analyysiin osoitti potentiaalin parantaa automaattisen merijääkartoituksen kartoituksen tarkkuutta tilastollisilla ja koneoppimiseen perustuvan luokittelun menetelmillä. Tulevaisuuden työnä merijään luokituksessa ja vallitutkimuksissa olisi suositeltavaa käyttää erilaisia ja useampia tutkakuvauksen geometrioita sekä erilaisia jääolosuhteita eri sääolosuhteiden vallitessa

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