Evaluation of wave retrieval for Chinese Gaofen-3 synthetic aperture radar

The goal of this study was to investigate the performance of a spectral-transformation wave retrieval algorithm and confirm the accuracy of wave retrieval from C-band Chinese Gaofen-3 (GF-3) Synthetic Aperture Radar (SAR) images. More than 200 GF-3 SAR images of the coastal China Sea and the Japan Sea for dates from January to July 2020 were acquired in the Quad-Polarization Strip (QPS) mode. The images had a swath of 30 km and a spatial resolution of 8 m pixel size. They were processed to retrieve Significant Wave Height (SWH), which is simulated from a numerical wave model called Simulating WAves Nearshore (SWAN). The first-guess spectrum is essential to the accuracy of Synthetic Aperture Radar (SAR) wave spectrum retrieval. Therefore, we proposed a wave retrieval scheme combining the theocratic-based Max Planck Institute Algorithm (MPI), a Semi-Parametric Retrieval Algorithm (SPRA), and the Parameterized First-guess Spectrum Method (PFSM), in which a full wave-number spectrum and a non-empirical ocean spectrum proposed by Elfouhaily are applied. The PFSM can be driven using the wind speed without calculating the dominant wave phase speed. Wind speeds were retrieved using a Vertical-Vertical (VV) polarized geophysical model function C-SARMOD2. The proposed algorithm was implemented for all collected SAR images. A comparison of SAR-derived wind speeds with European Center for Medium-Range Weather Forecasts (ECMWF) ERA-5 data showed a 1.95 m/s Root-Mean-Squared Error (RMSE). The comparison of retrieved SWH with SWAN-simulated results demonstrated a 0.47 m RMSE, which is less than the 0.68 m RMSE of SWH when using the PFSM algorithm.Output Status: Forthcoming/Available Onlin

Retrieval of the ocean wave spectrum in open and thin ice covered ocean waters from ERS Synthetic Aperture Radar images

This paper concerns with the task of retrieving ocean wave spectra from imagery provided by space-borne SAR systems such as that on board ERS satellite. SAR imagery of surface wave fields traveling into open ocean and into thin sea ice covers composed of frazil and pancake icefields is considered. The major purpose is to gain insight on how the spectral changes can be related to sea ice properties of geophysical interest such as the thickness. Starting from SAR image cross spectra computed from Single Look Complex (SLC) SAR images, the ocean wave spectrum is retrieved using an inversion procedure based on the gradient descent algorithm. The capability of this method when applied to satellite SAR sensors is investigated. Interest in the SAR image cross spectrum exploitation is twofold: first, the directional properties of the ocean wave spectra are retained; second, external wave information needed to initialize the inversion procedure may be greatly reduced using only information included in the SAR image cross spectrum itself. The main drawback is that the wind waves spectrum could be partly lost and its spectral peak wave number underestimated. An ERS-SAR SLC image acquired on April 10, 1993 over the Greenland Sea was selected as test image. A pair of windows that include open-sea only and sea ice cover, respectively, were selected. The inversions were carried out using different guess wave spectra taken from SAR image cross spectra. Moreover, care was taken to properly handle negative values eventually occurring during the inversion runs. This results in a modification of the gradient descending technique that is required if a non-negative solution of the wave spectrum is searched for. Results are discussed in view of the possibility of SAR data to detect ocean wave dispersion as a means for the retrieval of ice thickness

On the remote sensing of oceanic and atmospheric convection in the Greenland Sea by synthetic aperture radar

In this paper we discuss characteristic properties of radar signatures of oceanic and atmospheric convection features in the Greenland Sea. If the water surface is clean (no surface films or ice coverage), oceanic and atmospheric features can become visible in radar images via a modulation of the surface roughness, and their radar signatures can be very similar. For an unambiguous interpretation and for the retrieval of quantitative information on current and wind variations from radar imagery with such signatures, theoretical models of current and wind phenomena and their radar imaging mechanisms must be utilized. We demonstrate this approach with the analysis of some synthetic aperture radar (SAR) images acquired by the satellites ERS-2 and RADARSAT-1. In once case, an ERS-2 SAR image an a RADARSAT-1 ScanSAR image exhibit pronounced cell-like signatures with length scales on the order of 10-20 km and modulation depths of about 5-6 dB and 9-10 dB, respectively. Simulations with a numerical SAR imagaing model and various input current and wind fields reveal that the signatures in both images can be expained consistently by wind variations on the order of±2.5 ms, but not by surface current variations on realistic orders of magnitude. Accordingly, the observed features must be atmospheric convection cells. This is confirmed by visible typical cloud patterns in a NOAA AVHRR image of the test scenario. In another case, the presence of an oceanic convective chimney is obvious from in situ data, but no signatures of it are visible in an ERS-2 SAR image. We show by numerical simulations with an oceanic convection model and our SAR imaging model that this is consistent with theoretical predictions, since the current gradients associated with the observed chimney are not sufficiently strong to give rise to significant signatures in an ERS-2 SAR image under the given conditions. Further model results indicate that it should be generally difficult to observe oceanic convection features in the Greenland Sea with ERS-2 or RADARSAT-1 SAR, since their signatures resulting from pure wave-current interaction will be too weak to become visible in the noisy SAR images in most cases. This situation will improve with the availability of future high-resolution SARs such as RADARSAT-2 SAR in fine resolution mode (2004) and TerraSAR-X (2005) which will offer significantly reduced speckle noise fluctuations at comparable spatial resolutions and thus a much better visibility of small image variations on spatial scales on the order of a few hundred meters

A Marine Radar Wind Sensor

A new method for retrieving the wind vector from radar-image sequences is presented. This method, called WiRAR, uses a marine X-band radar to analyze the backscatter of the ocean surface in space and time with respect to surface winds. Wind direction is found using wind-induced streaks, which are very well aligned with the mean surface wind direction and have a typical spacing above 50 m. Wind speeds are derived using a neural network by parameterizing the relationship between the wind vector and the normalized radar cross section (NRCS). To improve performance, it is also considered how the NRCS depends on sea state and atmospheric parameters such as air–sea temperature and humidity. Since the signal-to-noise ratio in the radar sequences is directly related to the significant wave height, this ratio is used to obtain sea state parameters. All radar datasets were acquired in the German Bight of the North Sea from the research platform FINO-I, which provides environmental data such as wind measurements at different heights, sea state, air–sea temperatures, humidity, and other meteorological and oceanographic parameters. The radar-image sequences were recorded by a marine X-band radar installed aboard FINO-I, which operates at grazing incidence and horizontal polarization in transmit and receive. For validation WiRAR is applied to the radar data and compared to the in situ wind measurements from FINO-I. The comparison of wind directions resulted in a correlation coefficient of 0.99 with a standard deviation of 12.8°, and that of wind speeds resulted in a correlation coefficient of 0.99 with a standard deviation of 0.41 m s^−1. In contrast to traditional offshore wind sensors, the retrieval of the wind vector from the NRCS of the ocean surface makes the system independent of the sensors’ motion and installation height as well as the effects due to platform-induced turbulence

Evaluation of Chinese Quad-polarization Gaofen-3 SAR Wave Mode Data for Significant Wave Height Retrieval

Our work describes the accuracy of Chinese quad-polarization Gaofen-3 (GF-3) synthetic aperture radar (SAR) wave mode data for wave retrieval and provides guidance for the operational applications of GF-3 SAR. In this study, we evaluated the accuracy of the SAR-derived significant wave height (SWH) from 10,514 GF-3 SAR images with visible wave streaks acquired in wave mode by using the existing wave retrieval algorithms, e.g., the theoretical-based algorithm parameterized first-guess spectrum method (PFSM), the empirical algorithm CSAR_WAVE2 for VV-polarization, and the algorithm for quad-polarization (Q-P). The retrieved SWHs were compared with the European Centre for Medium-Range Weather Forecasts (ECMWF) reanalysis field with 0.125° grids. The root mean square error (RMSE) of the SWH is 0.57 m, found using CSAR_WAVE2, and this RMSE value was less than the RMSE values for the analysis results achieved with the PFSM and Q-P algorithms. The statistical analysis also indicated that wind speed had little impact on the bias with increasing wind speed. However, the retrieval tended to overestimate when the SWH was smaller than 2.5 m and underestimate with an increasing SWH. This behavior provides a perspective of the improvement needed for the SWH retrieval algorithm using the GF-3 SAR acquired in wave mode

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

Development of a satellite SAR image spectra and altimeter wave height data assimilation system for ERS-1

The applicability of ERS-1 wind and wave data for wave models was studied using the WAM third generation wave model and SEASAT altimeter, scatterometer and SAR data. A series of global wave hindcasts is made for the surface stress and surface wind fields by assimilation of scatterometer data for the full 96-day SEASAT and also for two wind field analyses for shorter periods by assimilation with the higher resolution ECMWF T63 model and by subjective analysis methods. It is found that wave models respond very sensitively to inconsistencies in wind field analyses and therefore provide a valuable data validation tool. Comparisons between SEASAT SAR image spectra and theoretical SAR spectra derived from the hindcast wave spectra by Monte Carlo simulations yield good overall agreement for 32 cases representing a wide variety of wave conditions. It is concluded that SAR wave imaging is sufficiently well understood to apply SAR image spectra with confidence for wave studies if supported by realistic wave models and theoretical computations of the strongly nonlinear mapping of the wave spectrum into the SAR image spectrum. A closed nonlinear integral expression for this spectral mapping relation is derived which avoids the inherent statistical errors of Monte Carlo computations and may prove to be more efficient numerically

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