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

High-resolution polar low winds obtained from unsupervised sar wind retrieval

High-resolution sea surface observations by spaceborne synthetic aperture radar (SAR) instruments are sorely neglected resources for meteorological applications in polar regions. Such radar observations provide information about wind speed and direction based on wind-induced roughness of the sea surface. The increasing coverage of SAR observations in polar regions calls for the development of SAR-specific applications that make use of the full information content of this valuable resource. Here we provide examples of the potential of SAR observations to provide details of the complex, mesoscale wind structure during polar low events, and examine the performance of two current wind retrieval methods. Furthermore, we suggest a new approach towards accurate wind vector retrieval of complex wind fields from SAR observations that does not require a priori wind direction input that the most common retrieval methods are dependent on. This approach has the potential to be particularly beneficial for numerical forecasting of weather systems with strong wind gradients, such as polar lows

Ocean surface currents derived from Sentinel-1 SAR Doppler shift measurements

Reliable information about ocean surface currents is crucial for operational oceanography, regulating weather development, and climate research (e.g., UN SDG 13). Upper-ocean currents are also key for monitoring life below water, including conservation of marine biodiversity at every trophic level (e.g., UN SDG 14). Locating upper ocean currents “with the right strength at the right place and time” is moreover critically needed to support the maritime transport sector, renewable marine energy, and maritime safety operations as well as for monitoring and tracking of marine pollution. In spite of this, upper ocean currents and their variability are mostly indirectly estimated and often without quantitative knowledge of uncertainties. In this thesis, Sentinel-1 Synthetic Aperture Radar (SAR) based Doppler frequency shift observations are examined for the retrievals of ocean surface current velocity in the radar line-of-sight direction. In the first study (Paper 1), Sentinel-1 A/B Interferometric Wide (IW) data acquired along the northern part of the Norwegian coastal zone from October-November 2017 at a spatial resolution of 1.5 km are compared with independent in-situ data, ocean model fields, and coastal High-Frequency Radar observations. Although only a limited dataset was available, the findings and results reveal that the strength of the meandering Norwegian Coastal Current derived from the SAR Doppler frequency shift observations are consistent with observations. However, limitations are encountered due to insufficient calibration and lack of ability to properly partition the geophysical signals into wave and current contributions. A novel approach for calibration of the attitude contribution to the Sentinel-1B Wave Mode (WV) Doppler frequency shift emerged for a test period in December 2017 - January 2018. Building on this calibrated dataset, an empirical model function (CDOP3S) for prediction of the sea state-induced contribution to the Doppler shift observations is developed for the global open ocean in Paper 2. The assessment against collocated surface drifter data are promising and suggest that the Sentinel-1B WV acquisitions can be used to study the equatorial ocean surface currents at a monthly timescale with a 20 km spatial resolution. The calibrated dataset combined with the new geophysical model function developed in Paper 2 also allowed for the study (Paper 3) of ocean surface current retrievals from the high-resolution Sentinel-1B IW swath data acquired along the coastal zone on northern Norway. In this case, the geophysical model function had to be trained and adjusted for fetch limited coastal sea state conditions. The results demonstrate that the Sentinel-1B SAR-derived ocean surface currents significantly improved, compared to the findings reported in Paper 1. Although the thesis builds on a limited period of observations, constrained by the availability of experimental attitude calibration, the results are all in all promising. Reprocessing of the full Sentinel-1 A/B SAR Doppler shift dataset using the novel attitude bias correction is therefore strongly recommended for further improvement of the empirical model function. Regular use of the Sentinel-1 A/B SAR for ocean surface current monitoring would thus be feasible, leading to advances in studies of upper ocean dynamics in support to the Copernicus Marine Environment Monitoring Service (CMEMS) program and the United Nations (UN) Decade of Ocean Sciences.Doktorgradsavhandlin

Remote Sensing Retrieval Study of the Surface Kinetic Parameters in the Yangtze Estuary and Its Adjacent Waters

Wind and current are significant parameters in the hydrodynamic processes, making a significant effect on the expansion of the Yangtze (Changjiang River) Diluted Water, sediment transport, resuspension and shelf circulation in the Yangtze Estuary. They are indispensable as input parameters in the numerical simulation of these phenomena. Synthetic aperture radar (SAR) can acquire data with different resolutions (down to 1 m) and coverage (up to 400 km) over a site during day or night time under all weather conditions, being capable of providing ocean surface kinetic parameters with high resolution. SAR images were collected to verify and improve the validity of wind direction retrieval by 2D fast Fourier transformation (FFT) method, wind speed by CMOD4 model and current by Doppler frequency method. These SAR-retrieved wind and current results were analyzed and assessed against in situ data and corresponding numerically simulated surface wind and current fields. Comparisons to the in situ and simulations show that 1) SAR can measure sea surface wind fields with a high resolution at sub-km scales and provide a powerful complement to conventional wind measurement techniques. 2) The Doppler shift anomaly measurements from SAR images are able to capture quantitative surface currents, thus are helpful to reveal the multi-scale upper layer dynamics around the East China Sea

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

Wind-Wave induced velocity in ATI SAR Ocean Surface Currents: First experimental evidence from an airborne campaign

Conventional and along-track interferometric (ATI) Synthetic Aperture Radar (SAR) sense the motion of the ocean surface by measuring the Doppler shift of reflected signals. Measurements are affected by a Wind-wave induced Artefact Surface Velocity (WASV) which was modelled theoretically in past studies and has been estimated empirically only once before with Envisat ASAR by Mouche et al., (2012). An airborne campaign in the tidally dominated Irish Sea served to evaluate this effect and the current retrieval capabilities of a dual-beam SAR interferometer known as Wavemill. A comprehensive collection of Wavemill airborne data acquired in a star pattern over a well-instrumented validation site made it possible for the first time to estimate the magnitude of the WASV, and its dependence on azimuth and incidence angle from data alone. In light wind (5.5 m/s) and moderate current (0.7 m/s) conditions, the wind-wave induced contribution to the measured ocean surface motion reaches up to 1.6 m/s upwind, with a well-defined 2nd order harmonic dependence on direction to the wind. The magnitude of the WASV is found to be larger at lower incidence angles. The airborne WASV results show excellent consistency with the empirical WASV estimated from Envisat ASAR. These results confirm that SAR and ATI surface velocity estimates are strongly affected by WASV and that the WASV can be well characterized with knowledge of the wind knowledge and of the geometry. These airborne results provide the first independent validation of Mouche et al., 2012, and confirm that the empirical model they propose provides the means to correct airborne and spaceborne SAR and ATI SAR data for WASV to obtain accurate ocean surface current measurements. After removing the WASV, the airborne Wavemill retrieved currents show very good agreement against ADCP measurements with a root mean square error (RMSE) typically around 0.1 m/s in velocity and 10° in direction

Towards Retrieving Reliable Ocean Surface Currents in the Coastal Zone From the Sentinel-1 Doppler Shift Observations

Recent developments on calibration and partitioning of the signal between the wave and current contributions significantly improve the accuracy of geophysical retrievals from Sentinel-1 Synthetic Aperture Radar-based Doppler shift measurements in the open ocean. In this study, we revise the Sentinel-1B Interferometric Wide products acquired from December 2017 to January 2018 along the coastal zone of northern Norway. We find that the satellite attitude is responsible for 30% of the variation in the Doppler shift observations, while the antenna pattern can describe an additional 15%. The residual variation after recalibration is about 3.8 Hz, corresponding to 0.21–0.15 m/s radial velocity (RVL) depending on the incidence angle. Using recalibrated Sentinel-1 observations, collocated with near-surface wind from MetCoOp-Ensemble Prediction System and sea state from MyWaveWAM, we develop an empirical function (CDOP3SiX) for estimating the sea-state-induced Doppler shift. CDOP3SiX improves the accuracy of sea state contribution estimates under mixed wind fetch conditions and demonstrates that the Norwegian Coastal Current can be detected in the Sentinel-1 derived ocean surface current RVL maps. Moreover, two anticyclonic mesoscale eddies with radial velocities of about 0.5 m/s are detected. The surface current patterns are consistent with the collocated sea surface temperature observations. The Doppler shift observations from Sentinel-1 can therefore be used to study ocean surface currents in the coastal zone with a 1.5 km spatial resolution. Key Points The Sentinel-1 Doppler shift observations are used to retrieve information about the ocean surface currents in the coastal zone Mesoscale eddies are detected in the Synthetic Aperture Radar-derived ocean surface current radial velocity fields Combination of the wind and wave information from collocated models improves the accuracy of the wave-induced contribution estimates Plain Language Summary Knowledge of ocean surface currents is crucial for studies of volume, heat and salt transport, tracking pollutants, and fisheries. The Doppler shift from Sentinel-1 Synthetic Aperture Radar (SAR) observations can be used to retrieve information about ocean surface currents. Challenging calibration and lack of algorithms for separating the wave and current contributions have limited the application of this observation-based method. Recent developments on calibration showed promising improvements in the accuracy of the signal. In this study, we apply this recent calibration method to Sentinel-1B scenes and develop an algorithm applicable for the challenging conditions in the coastal zone. We found that the signal from the Norwegian Coastal Current can be detected in the Sentinel-1 derived ocean surface current radial velocity fields. Also, we demonstrated the potential of SAR data for observing eddies with diameter of about 40–70 km. The Sentinel-1 derived surface currents express meandering structures and boundaries in consistence with the satellite-based sea surface temperature field. Comparison with the ocean model also reveals reasonable agreement, especially for the major surface current features. Therefore, given accurate calibration and new algorithm for removal of the wind and wave contribution, the Sentinel-1 observations can be used for monitoring ocean surface currents in the coastal zone with high spatial resolution.publishedVersio

GNSS transpolar earth reflectometry exploriNg system (G-TERN): mission concept

The global navigation satellite system (GNSS) Transpolar Earth Reflectometry exploriNg system (G-TERN) was proposed in response to ESA's Earth Explorer 9 revised call by a team of 33 multi-disciplinary scientists. The primary objective of the mission is to quantify at high spatio-temporal resolution crucial characteristics, processes and interactions between sea ice, and other Earth system components in order to advance the understanding and prediction of climate change and its impacts on the environment and society. The objective is articulated through three key questions. 1) In a rapidly changing Arctic regime and under the resilient Antarctic sea ice trend, how will highly dynamic forcings and couplings between the various components of the ocean, atmosphere, and cryosphere modify or influence the processes governing the characteristics of the sea ice cover (ice production, growth, deformation, and melt)? 2) What are the impacts of extreme events and feedback mechanisms on sea ice evolution? 3) What are the effects of the cryosphere behaviors, either rapidly changing or resiliently stable, on the global oceanic and atmospheric circulation and mid-latitude extreme events? To contribute answering these questions, G-TERN will measure key parameters of the sea ice, the oceans, and the atmosphere with frequent and dense coverage over polar areas, becoming a “dynamic mapper”of the ice conditions, the ice production, and the loss in multiple time and space scales, and surrounding environment. Over polar areas, the G-TERN will measure sea ice surface elevation (<;10 cm precision), roughness, and polarimetry aspects at 30-km resolution and 3-days full coverage. G-TERN will implement the interferometric GNSS reflectometry concept, from a single satellite in near-polar orbit with capability for 12 simultaneous observations. Unlike currently orbiting GNSS reflectometry missions, the G-TERN uses the full GNSS available bandwidth to improve its ranging measurements. The lifetime would be 2025-2030 or optimally 2025-2035, covering key stages of the transition toward a nearly ice-free Arctic Ocean in summer. This paper describes the mission objectives, it reviews its measurement techniques, summarizes the suggested implementation, and finally, it estimates the expected performance.Peer ReviewedPostprint (published version