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

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

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

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

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

In-depth verification of Sentinel-1 and TerraSAR-X geolocation accuracy using the Australian Corner Reflector Array

This article shows how the array of corner reflectors (CRs) in Queensland, Australia, together with highly accurate geodetic synthetic aperture radar (SAR) techniques—also called imaging geodesy—can be used to measure the absolute and relative geometric fidelity of SAR missions. We describe, in detail, the end-to-end methodology and apply it to TerraSAR-X Stripmap (SM) and ScanSAR (SC) data and to Sentinel-1interferometric wide swath (IW) data. Geometric distortions within images that are caused by commonly used SAR processor approximations are explained, and we show how to correct them during postprocessing. Our results, supported by the analysis of 140 images across the different SAR modes and using the 40 reflectors of the array, confirm our methodology and achieve the limits predicted by theory for both Sentinel-1 and TerraSAR-X. After our corrections, the Sentinel-1 residual errors are 6 cm in range and 26 cm in azimuth, including all error sources. The findings are confirmed by the mutual independent processing carried out at University of Zurich (UZH) and German Aerospace Center (DLR). This represents an improve�ment of the geolocation accuracy by approximately a factor of four in range and a factor of two in azimuth compared with the standard Sentinel-1 products. The TerraSAR-X results are even better. The achieved geolocation accuracy now approaches that of the global navigation satellite system (GNSS)-based survey of the CRs positions, which highlights the potential of the end-to-end SAR methodology for imaging geodesy

High-accuracy digital elevation model generation and ship monitoring from synthetic aperture radar images: innovative techniques and experimental results.

In this Thesis several state-of-the-art and innovative techniques for Digital Elevation Model (DEM) generation from Synthetic Aperture Radar (SAR) images are deeply analyzed, with a special focus on the methods which allow the improvement of the accuracy of the DEM product, which is directly related to the geolocation accuracy of geocoded images and is considered as an enabling factor for a large series of civilian and Defence applications. Furthermore, some of the proposed techniques, which are based both on phase and amplitude information, are experimented on real data, i.e. COSMO-SkyMed (CSK) data, assessing the achievable performances compared with the state-of-the-art, and pointing out and quantitatively highlighting the acquisition and processing strategies which would allow to maximize the quality of the results. Moreover, a critical analysis is performed about the main errors affecting the applied techniques, as well as the limitations of the orbital configurations, identifying several complementary techniques which would allow to overcome or mitigate the observed drawbacks. An innovative procedure for on-demand DEM production from CSK SAR data is elaborated and proposed, as well as an auto-validation technique which would enable the validation of the produced DEM also where vertical ground truths are not available. Based on the obtained results and on the consequent critical analysis, several interferometric specifications for new generation SAR satellites are identified. Finally, a literature review is proposed about the main state-of-the-art ship monitoring techniques, considered as one of the main fields of application which takes benefit from SAR data, based on single/multi-platform multi-channel SAR data, with a focus on TanDEM-X (TDX). In particular, in Chapter 1 the main concepts concerning SAR operating principles are introduced and the main characteristics and performances of CSK and TDX satellite systems are described; in Chapter 2 a review is proposed about the state-of-the-art SAR interferometric techniques for DEM generation, analyzing all the relevant processing steps and deepening the study of the main solutions recently proposed in the literature to increase the accuracy of the interferometric processing; in Chapter 3 complementary and innovative techniques respect to the interferometric processing are analyzed to mitigate disadvantages and to improve performances; in Chapter 4 experimental results are presented, obtained in the generation of high accuracy DEM by applying to a dataset of CSK images properly selected state-of-the-art interferometric techniques and innovative methods to improve DEM accuracy, exploring relevant limitations, and pointing out innovative acquisition and processing strategies. In Chapter 5, the basic principles of Ground Moving Target Indication (GMTI) are described, focusing on Displaced Phase Center Antenna (DPCA) and Along-Track Interferometry (ATI) techniques

Utilization of bistatic TanDEM-X data to derive land cover information

Forests have significance as carbon sink in climate change. Therefore, it is of high importance to track land use changes as well as to estimate the state as carbon sink. This is useful for sustainable forest management, land use planning, carbon modelling, and support to implement international initiatives like REDD+ (Reducing Emissions from Deforestation and Degradation). A combination of field measurements and remote sensing seems most suitable to monitor forests. Radar sensors are considered as high potential due to the weather and daytime independence. TanDEM-X is a interferometric SAR (synthetic aperture radar) mission in space and can be used for land use monitoring as well as estimation of biophysical parameters. TanDEM-X is a X-band system resulting in low penetration depth into the forest canopy. Interferometric information can be useful, whereas the low penetration can be considered as an advantage. The interferometric height is assumable as canopy height, which is correlated with forest biomass. Furthermore, the interferometric coherence is mainly governed by volume decorrelation, whereas temporal decorrelation is minimized. This information can be valuable for quantitative estimations and land use monitoring. The interferometric coherence improved results in comparison to land use classifications without coherence of about 10% (75% vs. 85%). Especially the differentiation between forest classes profited from coherence. The coherence correlated with aboveground biomass in a R² of about 0.5 and resulted in a root mean square error (RSME) of 14%. The interferometric height achieved an even higher correlation with the biomass (R²=0.68) resulting in cross-validated RMSE of 7.5%. These results indicated that TanDEM-X can be considered as valuable and consistent data source for forest monitoring. Especially interferometric information seemed suitable for biomass estimation

Parameters affecting interferometric coherence and implications for long-term operational monitoring of mining-induced surface deformation

Includes abstract.Includes bibliographical references.Surface deformation due to underground mining poses risks to health and safety as well as infrastructure and the environment. Consequently, the need for long-term operational monitoring systems exists. Traditional field-based measurements are point-based meaning that the full extent of deforming areas is poorly understood. Field-based techniques are also labour intensive if large areas are to be monitored on a regular basis. To overcome these limitations, this investigation considered traditional and advanced differential radar interferometry techniques for their ability to monitor large areas over time, remotely. An area known to be experiencing mining induced surface deformation was used as test case. The agricultural nature of the area implied that signal decorrelation effects were expected. Consequently, four sources of data, captured at three wavelengths by earth-orbiting satellites were obtained. This provided the opportunity to investigate different phase decorrelation effects on data from standard imaging platforms using real-world deformation phenomenon as test-case. The data were processed using standard dInSAR and polInSAR techniques. The deformation measurement results together with an analysis of parameters most detrimental to long-term monitoring were presented. The results revealed that, contrary to the hypothesis, polInSAR techniques did not provide an enhanced ability to monitor surface deformation compared to dInSAR techniques. Although significant improvements in coherence values were obtained, the spatial heterogeneity of phase measurements could not be improved. Consequently, polInSAR could not overcome ecorrelation associated with vegetation cover and evolving land surfaces. However, polarimetric information could be used to assess the scattering behaviour of the surface, thereby guiding the definition of optimal sensor configuration for long-term monitoring. Despite temporal and geometric decorrelation, the results presented demonstrated that mining-induced deformation could be measured and monitored using dInSAR techniques. Large areas could be monitored remotely and the areal extent of deforming areas could be assessed, effectively overcoming the limitations of field-based techniques. Consequently, guidelines for the optimal sensor configuration and image acquisition strategy for long-term operational monitoring of mining-induced surface deformation were provided

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