[go: up one dir, main page]

Next Article in Journal / Special Issue
Unambiguous Imaging of Static Scenes and Moving Targets with the First Chinese Dual-Channel Spaceborne SAR Sensor
Previous Article in Journal
A Brief Review of OPT101 Sensor Application in Near-Infrared Spectroscopy Instrumentation for Intensive Care Unit Clinics
Previous Article in Special Issue
Fast Vessel Detection in Gaofen-3 SAR Images with Ultrafine Strip-Map Mode
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Preliminary Assessment of Wind and Wave Retrieval from Chinese Gaofen-3 SAR Imagery

1
Marine Science and Technology College, Zhejiang Ocean University, Zhoushan 316000, China
2
Key Laboratory for Earth Observation of Hainan Province, Hainan 572029, China
3
Physical Oceanography Laboratory/CIMST, Ocean University of China and Qingdao National Laboratory for Marine Science and Technology, Qingdao 266100, China
*
Author to whom correspondence should be addressed.
Sensors 2017, 17(8), 1705; https://doi.org/10.3390/s17081705
Submission received: 30 May 2017 / Revised: 13 July 2017 / Accepted: 18 July 2017 / Published: 25 July 2017
(This article belongs to the Special Issue First Experiences with Chinese Gaofen-3 SAR Sensor)
Figure 1
<p>The quick-look image of calibrated Gaofen-3 (GF-)3 image in vertical–vertical (VV)-polarization around the Hawaiian islands acquired in Quad-Polarization Stripmap (QPS) mode at 16:22 UTC on 20 December 2016.</p> ">
Figure 2
<p>(<b>a</b>) European Centre for Medium-Range Weather Forecasts (ECMWF) wind field at 18:00 UTC and; (<b>b</b>) Significant wave height (SWH) field from WaveWatch-III at 18:00 UTC on 20 December 2016. The black rectangle represents the coverage of GF-3 synthetic aperture radar (SAR) image in <a href="#sensors-17-01705-f001" class="html-fig">Figure 1</a>.</p> ">
Figure 3
<p>SAR-derived wind field from VV-polarization GF-3 SAR image in Stander Stripmap (SS) mode at 02:17 UTC on 29 September 2016, in which the white circle represents the location of National Data Buoy Center (NDBC) in situ buoy (ID: 46013).</p> ">
Figure 4
<p>The comparison between SAR-derived wind speed U<sub>10</sub> and measurements from NDBC buoys. (<b>a</b>) Comparison for 16 VV-polarization GF-3 SAR images; and (<b>b</b>) Comparison for 42 horizontal–horizontal (HH)-polarization GF-3 SAR images.</p> ">
Figure 5
<p>(<b>a</b>) The two-dimensional SAR spectrum of sub-scene covering the buoy (ID:46013) in <a href="#sensors-17-01705-f003" class="html-fig">Figure 3</a>; (<b>b</b>) The corresponding one-dimensional spectrum with the fitted result by using Gaussian fit function.</p> ">
Figure 6
<p>The comparison between SAR-derived H<sub>s</sub> and the measurements from NDBC buoys. (<b>a</b>) Comparison for 16 VV-polarization GF-3 SAR images; and (<b>b</b>) Comparison for 42 HH-polarization GF-3 SAR images.</p> ">
Figure 7
<p>The comparison between SAR-derived wind speed U<sub>10</sub> and ECMWF re-analysis gridded winds at 1 m/s bins. (<b>a</b>) Comparison for 96 VV-polarization GF-3 SAR images; and (<b>b</b>) Comparison for 70 HH-polarization GF-3 SAR images.</p> ">
Figure 8
<p>The comparison between SAR-derived SWH H<sub>s</sub> and WaveWatch-III data at 0.5 m/s bins. (<b>a</b>) Comparison of 96 VV-polarization GF-3 SAR images; and (<b>b</b>) Comparison of 70 HH-polarization GF-3 SAR images.</p> ">
Versions Notes

Abstract

:
The Chinese Gaofen-3 (GF-3) synthetic aperture radar (SAR) launched by the China Academy of Space Technology (CAST) has operated at C-band since September 2016. To date, we have collected 16/42 images in vertical-vertical (VV)/horizontal-horizontal (HH) polarization, covering the National Data Buoy Center (NDBC) buoy measurements of the National Oceanic and Atmospheric Administration (NOAA) around U.S. western coastal waters. Wind speeds from NDBC in situ buoys are up to 15 m/s and buoy-measured significant wave height (SWH) has ranged from 0.5 m to 3 m. In this study, winds were retrieved using the geophysical model function (GMF) together with the polarization ratio (PR) model and waves were retrieved using a new empirical algorithm based on SAR cutoff wavelength in satellite flight direction, herein called CSAR_WAVE. Validation against buoy measurements shows a 1.4/1.9 m/s root mean square error (RMSE) of wind speed and a 24/23% scatter index (SI) of SWH for VV/HH polarization. In addition, wind and wave retrieval results from 166 GF-3 images were compared with the European Centre for Medium-Range Weather Forecasts (ECMWF) re-analysis winds, as well as the SWH from the WaveWatch-III model, respectively. Comparisons show a 2.0 m/s RMSE for wind speed with a 36% SI of SWH for VV-polarization and a 2.2 m/s RMSE for wind speed with a 37% SI of SWH for HH-polarization. Our work gives a preliminary assessment of the wind and wave retrieval results from GF-3 SAR images for the first time and will provide guidance for marine applications of GF-3 SAR.

1. Introduction

The Chinese Gaofen-3 (GF-3) satellite was launched on 10 August 2016 by the China Academy of Space Technology (CAST), and carries a C-band (~5.3 GHz) synthetic aperture radar (SAR) sensor with different polarizations. It has 12 imaging modes with a spatial resolution of image ranging from 1 m to 500 m and a swath coverage ranging from 10 to 650 km. GF-3 SAR operates in different polarizations, including single-, dual- and quad-polarization. Through a cooperation project between the National Satellite Ocean Application Service (NSOAS) and our institutes, a number of images have been recorded of whole open seas, particularly covering the National Data Buoy Center (NDBC) in situ buoys of the National Oceanic and Atmospheric Administration (NOAA) around U.S. western coastal areas. Wind and wave monitoring are the two main marine applications of SAR in all-weather conditions, especially in tropical cyclones [1,2].
It is well known that Bragg waves, which have the wavelength with the order of centimeters, backscatter the microwaves of SAR. Bragg waves are produced by sea surface winds. Through studying co-polarization (vertical-vertical (VV) and horizontal-horizontal (HH)) Spaceborne Imaging Radar (SIR) data on Seasat launched in 1978 [3,4], the geophysical model function (GMF) has been exploited. GMF describes an empirical relationship between normalized radar cross-section (NRCS) in VV-polarization and a wind vector. C-band GMFs, e.g., CMOD4 [5], CMOD-IFR developed at Institut Francais de Recherche pour Exploitation de la MER (IFREMER) [6], CMOD5 [7] and CMOD5N [8]. These GMFs provide a convenient application for wind retrieval from SAR [9,10,11] and have been successfully implemented for various C-band SAR data over the last few decades, e.g., ERS-1/2 [12], Envisat-ASAR [13], Radarsat-1/2 [14] and Sentinel-1A/1B [15] within about 2 m/s error of wind speed. CMOD4 and CMOD-IFR work at wind speeds smaller than 20 m/s, due to no higher wind source being available in the tuning process. The formulation of an improved C-band GMF CMOD5 was essentially redesigned with a number of ERS-2 images taken in tropical cyclones and corresponding European Centre for Medium-Range Weather Forecasts (ECMWF) re-analysis winds. In particular, some case studies have shown that CMOD5 can be used in hurricanes [16,17,18] to some extent. Later, CMOD5 was retuned for neutral winds, denoted as CMOD5N, which takes a correction for CMOD5 can be more stable applied. As GMF is applied for HH-polarization SAR, NRCS in HH-polarization has to be converted into NRCS in VV-polarization by using a polarization ratio (PR) model [19,20,21,22,23]. The latest achievement of the PR model was proposed by the authors of [23] and the improvement is that the dependence of sea surface wind speed on PR is included in the model. More recently, using collocated scatterometer measurements on aboard Metop-A/B and ECMWF winds with co-polarization Envisat-ASAR and Sentinel-1A/1B SAR data, a new C-band GMF for wind retrieval from co-polarization SAR was developed in [24], denoted as C-SARMOD. Although the accuracy of wind speeds retrieved from Sentinel-1A/1B images in VV-polarization by using C-SARMOD showed a 1.6 m/s STD of wind speed against moored buoy measurements [25], the validation of its application for C-band SAR in HH-polarization has not yet been systematically investigated.
Traditionally, the methodology of wave retrieval from SAR is based on the SAR mapping mechanism [26]. The major modulations of waves on SAR have been well studied over many decades, e.g., tilt modulation, hydrodynamic modulation and velocity bunching [27,28]. Both tilt and hydrodynamic modulations are linear mapping mechanisms and the two modulation transfer functions (MTF) were proposed by the authors of [29,30]. However, velocity bunching is caused by the relative motion between satellite platform and sea surface, resulting in a Doppler frequency shift in the azimuth direction (the direction of satellite flight is defined as the azimuth direction and radar look direction is the range direction). This velocity bunching causes wavelengths smaller than a specific value (or cutoff wavelength) in the azimuth direction, which is not detectable and the peak of the SAR spectrum rotates toward the range direction. These two effects create the difficulty of directly inverting SAR intensity spectra to wave spectra. The first solution was established by the authors of [28], with an algorithm known as the Max–Planck Institute (MPI) algorithm. The basic scheme of the MPI is described as follows: (1) a first-guess wave spectrum is produced from numeric wave modes, e.g., Simulating WAves Nearshore (SWAN) and WaveWatch-III; (2) a simulated SAR spectrum is obtained by mapping the first-guess wave spectrum; (3) the reality of the wave spectrum is inverted by minimizing the simulated and real SAR spectrum through a cost function. The Semi Parametric Retrieval Algorithm (SPRA) [31] is more applicable than the MPI, because the SPRA employs wind measurements from a scatterometer in order to produce the first-guess spectrum by using the Jonswap parametric model. Later the Parameterized First-guess Spectrum Method (PFSM) algorithm was proposed by the authors of [32,33] and recent studies show that the PFSM works for X-band TerraSAR-X [34] and C-band Sentinel-1 SAR data [25]. The improvement of the PFSM involves the wind-wave and swell information on SAR being separated by calculating the wave number threshold of the SAR intensity spectra. A 0.54 m standard deviation (STD) of SWH is exhibited in study [25] by using the algorithm PFSM as the SAR-derived SWH from Sentinel-1 SAR data validated against the buoys. In addition, there are other wave algorithms, e.g., the Partition Rescaling and Shift Algorithm (PARSA) [35] for SAR complex data and unconstrained algorithms [36,37], in which the MTF of velocity bunching is bypassed and it only works for long wave dominated regions.
Empirical wave retrieval algorithms, without calculating complex MTF, were developed by the SAR group at the German Aerospace Center (DLR), and include CWAVE_ERS for ERS-2 SAR [38], CWAVE_ENV for Envisat-ASAR [39] and XWAVE for X-band TerraSAR-X SAR [40]. CWAVE allows direct retrieval of significant wave height (SWH) from SAR wave mode data at a fixed incidence angle of about 23°. However, these are not conducive to the operational application of wave retrieval from various SAR data. Interestingly, the cutoff wavelength in azimuth direction caused by velocity bunching is theoretically related with SWH [28]. Several recent studies have made great efforts to retrieve wave parameters through the cutoff wavelength [41,42,43,44,45]. In our previous study [46], the four existing algorithms, including PFSM and three other empirical algorithms, have been implemented for HH-polarization Sentinel-1 SAR images. It was found that the empirical algorithm, herein called CSAR_WAVE, has a good performance as the retrieval results compared with moored buoys.
In this study, we give a preliminary assessment of wind and wave retrieval from the new Chinese C-band GF-3 SAR for the first time. After employing ECMWF wind direction, wind speed is retrieved from a VV-polarization GF-3 SAR image using CMOD5N. The PR model proposed in [23] is used together with CMOD5N for wind speed retrieval from an HH-polarization GF-3 SAR image. SWH is retrieved using CSAR_WAVE for both VV and HH polarization GF-3 SAR without any prior knowledge.
The remainder of this paper is organized as follows. The description of collected C-band co-polarization GF-3 SAR images and the validation sources, including NOAA buoy measurements, available ECMWF re-analysis winds and wave computations from the WaveWatch-III model provided by the IFREMER group, are briefly introduced in Section 2. Section 3 shows the methodology of wind and wave retrieval algorithms used in this study. Then retrieval results and discussions are presented in Section 4 and Section 5, respectively. Conclusions are summarized in Section 6.

2. Description of Datasets

In total, we collected 224 GF-3 SAR images in co-polarization (VV and HH polarization) through September 2016 to March 2017 at seas. These single look complex (SLC) images were acquired in Stander Stripmap (SS) or Quad-Polarization Stripmap (QPS) mode. In the 224 images, there are 166 images, which are a matchup with 0.125 × 0.125° grids ECMWF re-analysis winds at intervals of six hours and 0.5 × 0.5° grids waves from the WaveWatch-III model at intervals of 3 h. The time difference between the imaging time of those images and ECMWF/WaveWatch-III data was within two hours. They were used to investigate the accuracy of wind and wave retrieval results in our study. The following equation is used for calculating the NRCS of co-polarization GF-3 SAR intensity image.
σ 0 = DN 2 ( M 32767 ) 2 N
wherein σ° is the NRCS united in dB, DN is the intensity derived from GF-3 SAR Level-1A data, M is the external calibration factor and N is the offset constant stored in the annotation file. As an example, a quick-look image of a calibrated GF-3 SAR image in VV-polarization around the Hawaiian islands acquired in QPS mode at 16:22 UTC on 20 December 2016 is shown in Figure 1. The ECMWF wind field at 18:00 UTC is shown in Figure 2a and SWH from the WaveWatch-III model at 18:00 UTC are shown in Figure 2b. The black rectangle represents the coverage of the case image. In the matchups, we only used homogenous sub-scenes derived from images for wave validation, that is, the good-quality SAR intensity spectra can be obtained by using the two-dimensional Fast Fourier Transform (FFT) method.
Out of the 166 GF-3 SAR image matchups, 58 images cover the NDBC in situ buoys of NOAA around U.S. western coastal areas. It is necessary to calculate that the wind speeds measured by the NDBC in situ buoys are at a height of 5 m above the sea surface, so we use Equation (2) to convert buoy-measured wind speeds to values at 10 m height as neutral winds,
U 10 U 5   =   ln ( 10 / z 0 ) ln ( 5 / z 0 )
wherein U10 is the wind speed at 10 m height, U5 is the wind speed measured by NDBC in situ buoy and z0 is the roughness length taken as a constant 1.52 × 10−4 [25,46].

3. Wind and Wave Retrieval Algorithms for C-Band SAR

In this section, the methodology of wind and wave retrieval from SAR is presented, including the C-band GMF with PR model for wind retrieval and a cutoff-wavelength-based empirical algorithm for wave retrieval.

3.1. Wind Retreival Algorithm

Until now, the CMOD family has successsfully been applied for wind retreival from various C-band SAR data [9,10,11]. The CMOD family takes the general formulation:
σ 0 = B 0 ( 1   +   B 1 cos ϕ   +   B 2 cos 2 ϕ ) p
wherein σ° is the linear SAR-measured NRCS, ϕ is the wind direction related to range direction, p is a parametric constant, B0 to B2 are the functions of wind speed at 10 m above sea surface U10 and radar incidence angle is θ. CMOD5N is the latest version and the accuracy of wind speeds retrieved from a number of Sentinel-1 SAR images using CMOD5N have been investigated in [15], showing a 1.35 m/s STD of wind speed validated against scatterometer measurements on aboard Metop-A/B. Because there are two unknown variables, e.g., wind speed and wind direction, it is impossible to solve Equation (3) to invert the wind vector. It has been found that homogenous wind streaks at the kilometer scale, which are parallel to the wind direction, can be retrieved from some SAR images [47]. However, a SAR-derived wind direction has a 180° ambiguity and external information is required so as to remove that ambiguity. Moreover, wind streaks do not appear in SAR images. Therefore, we directly employed ECMWF wind directions in the wind retrieval process.
PR models, which are described as the ratio between NRCS in VV and HH polarization, are usually used together with CMODs as wind retrieval for C-band HH-polarization SAR [19,20,21,22,23]. Recent research has shown that C-band SAR PR has a linear relation with wind speed [23], in addition to the incidence angle. It is not surprising that the PR model, which includes dependence of sea surface wind speed, performs better than other PR models involving only the dependence of incidence angle, which is stated as follows:
PR   =   σ VV 0 σ HH 0   =   P ( θ ) U 10 Q ( θ )
where
P ( θ )   =   P 1 θ 2   +   P 2 θ   +   P 3
and
Q ( θ )   =   Q 1 θ   +   Q 2
wherein σ VV 0 and σ HH 0 are the SAR-measured linear NRCS in VV- and HH-polarization respectively, coefficients P1 to P3 and Q1 to Q2 are tuned by an amount of quad-polarization Radarsat-2 images and the collocated buoy measurements. The dependence of wind speed on X-band PR has been investigated through dual-polarization TerraSAR-X images and ECMWF re-analysis winds [48] following the PR model proposed by the authors of [23], which was adopted for X-band SAR.

3.2. Wave Retreival Algorithm

As introduced in Section 1, the many SAR mapping mechanism-based algorithms [28,30,31,32,33,34,35,36] need a prior ‘first-guess’ wave spectrum provided from the numeric wave mode or produced by the parametric wave function using wind speed. Therefore, none of them can be operationally applied for wave retrieval from SAR. Although the two empirical algorithms, e.g., CWAVE_ERS [39] and CWAVE_ENV [40], work well for SAR data at a specific incidence angle of about 23°, they have yet to be validated for SAR data in imaging mode at various incidence angles. Based on the theoretical relationship between cutoff wavelength in azimuth direction and SWH [28,30], a new empirical algorithm, denoted as CSAR_WAVE, was exploited through 93 C-band Stripmap mode Sentinel-1 SAR images in VV-polarization and collocated wave measurements from NDBC in situ buoys [45]. The validation against buoy measurements showed a 18.6% SI of SWH in our previous study. Recently, CSAR_WAVE has been adopted for HH-polarization Sentinel-1 SAR in [46]. Moreover, CSAR_WAVE provides a convenient empirical method to retrieve SWH for C-band SAR, including, but not limited to, StripMap mode data. The functions of CSAR_WAVE are designed as follows,
H s   =   ( λ c β ) ( A 1   +   A 2 sin θ   +   A 3 cos 2 φ )   +   A 4
β   =   R V
wherein Hs is the SWH, λc is the cutoff wavelength in azimuth direction, β is the satellite range-to-velocity parameter, R is the slant range, V is the satellite flight velocity, θ is the radar incidence angle, φ is the wave propagation angle relative to range direction ranging from 0 to 90° and the coefficients A1 to A4 are determined from the C-band Sentinel-1 SAR image and collocated NDBC in situ buoys [45] and ECMWF re-analysis wave data at 0.125° [46]. It is necessary to figure out that CSAR_WAVE can operate without any prior knowledge. The advantage of CSAR_WAVE is its application can be implemented without using SAR-derived wind speeds. Therefore, wind and wave are simultaneously measured from co-polarization GF-3 SAR data in our work.

4. Method and Results

In this study, wind and wave retrieval results from 58 GF-3 images in co-polarization were validated against NDBC in situ buoy measurements around the U.S. western coastal area, including wind speed and SWH. In order to perform the matchup, each GF-3 imagery was divided into several sub-scenes with a spatial coverage of about 3 × 3 km in azimuth and range direction, respectively. We chose the sub-scenes, covering the moored buoys for studying the accuracy of winds and SWH retrieval results from co-polarization GF-3 images.

4.1. Validation of Wind Retreival Results

After employing ECMWF wind directions, wind speeds were retrieved from VV-polarization GF-3 SAR images by using CMOD5N. The PR model involving Equations (4)–(6) is used together with CMOD5N for wind speeds retrieval from HH-polarization GF-3 SAR images. As an example, the information retrieved from VV-polarization GF-3 SAR image in SS mode at 02:17 UTC on 29 September 2016 is shown in Figure 3. The SAR averaged wind speed located at the sub-scene of about 3 × 3 km coverage covering the NDBC in situ buoy (ID: 46013) is 11.6 m/s and the buoy-measured wind speed is 10.8 m/s. The difference between the SAR-derived wind speed and the buoy measurement is only 0.8 m/s.
We collected a number of sub-scenes from 16/42 GF-3 images in VV/HH polarization matchup with NDBC in situ buoys. The comparison of the matchups is shown in Figure 4. Validation shows a 1.4 m/s/1.9 m/s root mean square error (RMSE) of wind speed with a bias of −0.4 m/s for VV/HH polarization. ECMWF re-analysis winds are popularly used for tuning and validating the wind retrieval algorithms for SAR [7,8]. The accuracy of retrieved wind speeds from co-polarization GF-3 images has a good performance with around 2 m/s RMSE of wind speed. Therefore, it is found that winds retrieval from GF-3 SAR have a similar accuracy to that of the other C-band SARs, which have a up to 1.78 m/s stander deviation (STD) for wind speeds as validated against buoys, scatterometer and numeric models [10,13,14].

4.2. Validation of Wave Retreival Results

In this study, we used CSAR_WAVE to retrieve SWH from co-polarization GF-3 SAR images. Two variables, including wave propagation angle relative to range direction φ and cutoff wavelength in azimuth direction λc, are derived from the SAR spectrum. The two-dimensional SAR spectrum is calculated from a SAR intensity image by using the two-dimensional FFT method. φ is directly obtained from a two-dimensional SAR spectrum, which ranges from 0 to 90° in the CSAR_WAVE model. Then, we employ a Gaussian fit function to integrate the two-dimensional SAR spectrum in the range direction. The Gaussian fit function is stated as,
exp { π k × λ c 2 π }
in which, kx is the wavenumber in the azimuth direction.
The two-dimensional SAR spectrum of sub-scene covering the NDBC in situ buoy (ID: 46013) in Figure 3 and the corresponding one-dimensional spectrum with the fitted result are shown in Figure 5a,b, respectively. The retrieved SWH is 2.2 m and the buoy-measured SWH is 2.1 m. Figure 6a,b show the further comparisons between SAR-derived SWHs and measurements from NDBC in situ buoys, showing a 24%/23% scatter index (SI) of SWH with a 0.58 m/0.57 m RMSE of SWH for VV/HH polarization.
The retrieval SWH from co-polarization GF-3 SAR using CSAR_WAVE has a close accuracy to that from other C-band SAR by using theoretical-based algorithms and empirical CWAVE models, e.g., about 20% SI of SWH validated against moored buoys or measurements from altimeter [25,38,39] and around 38% SI of SWH as compared to the SWH retrieval results with simulation results from the WAM model [49].

5. Discussion

In order to further investigate the accuracy of winds and waves retrieved from GF-3 SAR images. We also compared SAR-derived wind speeds with 0.125 × 0.125° grids ECMWF re-analysis wind speeds and SAR-derived SWH with those from the 0.5 × 0.5° grids WaveWatch-III model through more GF-3 co-polarization images at seas.
Out of 58 GF-3 SAR images matchup with buoys, there are additional 166 images in our collections, including 96 images in VV-polarization and 70 images in HH-polarization. The retrieval wind speeds from these images were compared with ECMWF re-analysis gridded winds at 1 m/s bins, showing a 2.0 and 2.2 m/s RMSE of wind speed for VV-polarization and HH-polarization respectively, as exhibited in Figure 7.
The computations run by the WaveWatch-III model at intervals of three hours. In our previous study, that data has been used for the validation of wave retrieval results from X-band SAR by using algorithm PFSM [34]. Although the open-accessed simulations from the WaveWatch-III model have a lower spatial resolution of 0.5° grid, the model results can be used for independent validation purposes. We applied the CSAR_WAVE model to a total of 96/70 GF-3 SAR images in VV/HH polarization and then compared the retrieval results with those from the WaveWatch-III model at 0.5 m bins of SWH. Figure 8 shows the SI of SWH is 36% with a 0.74 m RMSE for VV-polarization and the SI of SWH is 37% with a 0.74 m RMSE for HH-polarization.
As a result, it is found that SAR-derived winds have a good agreement with ECMWF re-analysis winds for wind speeds between 5 and 10 m/s. However, SAR-derived winds are larger than ECMWF re-analysis winds for wind speeds smaller than 5 m/s and wind speeds greater than 10 m/s. We think it is probably due to the change of atmospheric stability at low winds and high wind conditions. Compared with SWH from WaveWatch-III model, the SWH retrieval results are generally larger. It is necessary to figure out that herein statistical results have a larger error than the result as exhibited in Figure 6, due to the coarse spatial resolution of wave data from WaveWatch-III model. To give a better performance, the bias can somehow be improved with the simple expediency of subtracting 0.6 for VV-polarization or 0.7 for HH-polarization.

6. Conclusions

Wind speed and SWH are two of most the important parameters in oceanography research, and SAR has the capability to observe wind and wave in wide spatial coverage. The Chinese GF-3 satellite carries a C-band SAR sensor, and has been operating in 12 imaging modes with a fine spatial resolution of up to 1 m since September 2016. Recently, the validation of SAR-derived winds from GF-3 images has been presented in [50]. In our work, the accuracy of wind speed and SWH retrieval results from GF-3 SAR have been simultaneously investigated.
We employed wind directions directly from ECMWF re-analysis data. Then CMOD5N, together with the PR model, involving the dependence on wind speed and incidence angle, were used for retrieving wind speeds. Validations against NDBC in situ buoys showed a 1.4 and 1.9 m/s RMSE of wind speed through 16 VV-polarization and 42 HH-polarization GF-3 SAR images. SAR-derived wind speeds from an additional 96/70 GF-3 images in VV/HH polarization were compared with ECMWF re-analysis winds, showing a 2.0/2.2 m/s RMSE of wind speed, respectively.
Theoretical-based algorithms and empirical CWAVE models have been developed in recent years. However, all of these rely on either wind speed or computations from numeric wave models. The new empirical CSAR_WAVE model is designed based on the relation between SAR cutoff wavelength in azimuth direction and SWH, and can be applicable without any prior information. The comparisons between SWH retrieval results and measurements from NDBC in situ buoys show a 24 and 23% SI of SWH with a 0.58 and 0.57 m RMSE of SWH for VV-polarization and HH polarization, respectively. We also compared the SWH retrieval results with the computations from the WaveWatch-III model, showing a 36/37% SI of SWH through 96/70 GF-3 images in VV/HH polarization. Therefore, CSAR_WAVE is suitable for wave retrieval from Chinese C-band GF-3 SAR.
In summary, the proposed methods are operationally applicable for wind and wave retrieval from C-band GF-3 SAR images in co-polarization. Moreover, the independent extraction of wind and wave from co-polarization GF-3 SAR images due to SAR-derived wind speed is no longer required in the process of wave retrieval when using the empirical CSAR_WAVE model. It is concluded that the performance achieves the design requirements of GF-3 SAR referred to the preliminary assessment of winds and waves presented in this study.

Acknowledgments

We appreciate the cooperation project between NSOAS and our institutes. And the GF-3 data are downloaded via http://dds.nsoas.org.cn as authorized account issued by NSOAS under the contract of Specific Project of Chinese High Resolution Earth Observation System (No. 41-Y20A14-9001-15/16). Xinzhe Yuan and Limin Cui provide great help in order to deliver and process the GF-3 data. NOAA buoy measurements were accessed via http://www.ndbc.noaa.gov/. ECMWF re-analysis wind data were openly accessed via http://www.ecmwf.int. Wave data from the WaveWatch-III model were kindly provided by the professional IFREMER group and the data were collected through the server: ftp.ifremer.fr. This research is partly supported by Hainan Provincial Department of Science and Technology under grant No. ZDKJ2016021, the National Natural Science Foundation of China under grant No. 41376010.

Author Contributions

Weizeng Shao and Jian Sun came up the original idea and designed the experiments; Weizeng Shao and Yexin Shen analyzed the data; Yexin Shen collected the data; all authors contributed to the writing and revising of the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Li, X.F.; Pichel, W.G.; He, M.X.; Wu, S.; Friedman, K.; Clemente-Colon, P.; Zhao, C.F. Observation of hurricane-generated ocean swell refraction at the gulf stream north wall with the RADARSAT-1 synthetic aperture radar. IEEE Trans. Geosci. Remote Sens. 2002, 40, 2131–2142. [Google Scholar]
  2. Li, X.F.; Zhang, J.A.; Yang, X.F.; Pichel, W.G.; DeMaria, M.; Long, D.; Li, Z.W. Tropical cyclone morphology from spaceborne synthetic aperture radar. Bull. Am. Meteorol. Soc. 2013, 94, 215–230. [Google Scholar] [CrossRef]
  3. Gerling, T. Structure of the surface wind field from the Seasat SAR. J. Geophys. Res. 1986, 91, 2308–2320. [Google Scholar] [CrossRef]
  4. Masuko, H.; Okamoto, K.; Shimada, M.; Niwa, S. Measurement of microwave backscattering signatures of the ocean surface using X band and Ka band airborne scatterometers. J. Geophys. Res. 1986, 91, 13065–13083. [Google Scholar] [CrossRef]
  5. Stoffelen, A.; Anderson, D. Scatterometer data interpretation: Derivation of the transfer function CMOD4. J. Geophys. Res. 1997, 102, 5767–5780. [Google Scholar] [CrossRef]
  6. Quilfen, Y.; Bentamy, A.; Elfouhaily, T.M.; Katsaros, K.; Tournadre, J. Observation of tropical cyclones by high-resolution scatterometry. J. Geophys. Res. 1998, 103, 7767–7786. [Google Scholar] [CrossRef]
  7. Hersbach, H.; Stoffelen, A.; Haan, S.D. An improved C-band scatterometer ocean geophysical model function: CMOD5. J. Geophys. Res. 2007, 112. [Google Scholar] [CrossRef]
  8. Hersbach, H. Comparison of C-Band scatterometer CMOD5. N equivalent neutral winds with ECMWF. J. Atmos. Ocean. Technol. 2010, 27, 721–736. [Google Scholar] [CrossRef]
  9. Lehner, S.; Horstmann, J. Mesoscale wind measurements using recalibrated ERS SAR images. J. Geophys. Res. 1998, 103, 7847–7856. [Google Scholar] [CrossRef]
  10. Monaldo, F.M.; Thompson, D.R.; Pichel, W.G.; Clemente-Colon, P. A systematic comparison of QuikSCAT and SAR ocean surface wind speeds. IEEE Trans. Geosci. Remote Sens. 2004, 42, 283–291. [Google Scholar] [CrossRef]
  11. Shao, W.Z.; Sun, J.; Guan, C.L.; Sun, Z.F. A method for sea surface wind field retrieval from SAR image mode data. J. Ocean Univ. China 2014, 13, 198–204. [Google Scholar] [CrossRef]
  12. Lehner, S.; Schulz-Stellenfleth, J.; Schattler, B.; Breit, H.; Horstmann, J. Wind and wave measurements using complex ERS-2 SAR wave mode data. IEEE Trans. Geosci. Remote Sens. 2000, 38, 2246–2257. [Google Scholar] [CrossRef]
  13. Yang, X.F.; Li, X.F.; Pichel, W.; Li, Z.W. Comparison of ocean surface winds from ENVISAT ASAR, Metop ASCAT scatterometer, buoy Measurements, and NOGAPS model. IEEE Trans. Geosci. Remote Sens. 2011, 49, 4743–4750. [Google Scholar] [CrossRef]
  14. Yang, X.F.; Li, X.F.; Zheng, Q.A.; Gu, X.; Pichel, W.; Li, Z.W. Comparison of ocean-surface winds retrieved from QuikSCAT scatterometer and Radarsat-1 SAR in offshore waters of the U.S. west coast. IEEE Geosci. Remote Sens. Lett. 2011, 8, 163–167. [Google Scholar] [CrossRef]
  15. Monaldo, F.M.; Jackson, C.; Li, X.F.; Pichel, W.G. Preliminary evaluation of Sentinel-1A wind speed retrievals. IEEE J. Sel. Top. Appl. Earth Observ. Remote Sens. 2016, 9, 2638–2642. [Google Scholar] [CrossRef]
  16. Horstmann, J.; Thompson, D.R.; Monaldo, F.; Iris, S.; Graber, H.C. Can synthetic aperture radars be used to estimate hurricane force winds? Geophys. Res. Lett. 2005, 32, 1–5. [Google Scholar] [CrossRef]
  17. Shen, H.; Perrie, W.G.; He, Y.J. A new hurricane wind retrieval algorithm for SAR images. Geophys. Res. Lett. 2006, 33, 221–256. [Google Scholar] [CrossRef]
  18. Li, X.F. The first Sentinel-1 SAR image of a typhoon. Acta Oceanol. Sin. 2015, 34, 1–2. [Google Scholar] [CrossRef]
  19. Thompson, D.R.; Elfouhaily, T.M.; Chapron, B. Polarization ratio for microwave backscattering from the ocean surface at low to moderate incidence angles. In Proceedings of the IEEE International Geoscience and Remote Sensing Symposium, Seattle, WA, USA, 6–10 July 1998; pp. 1671–1673. [Google Scholar]
  20. Vachon, P.W.; Dobson, F.W. Wind retrieval from RADARSAT SAR images selection of a suitable C-band HH polarization wind retrieval model. Can. J. Remote Sens. 2000, 26, 2122–2131. [Google Scholar] [CrossRef]
  21. Horstmann, J.; Koch, W.; Lehner, S.; Tonboe, R. Wind retrieval over the ocean using synthetic aperture radar with C-band HH polarization. IEEE Trans. Geosci. Remote Sens. 2000, 38, 306–313. [Google Scholar] [CrossRef]
  22. Mouche, A.A.; Hauser, D.; Daloze, J.F.; Guerin, C. Dual polarization measurements at C-band over the ocean: Results from airborne radar observations and comparison with ENVISAT ASAR data. IEEE Trans. Geosci. Remote Sens. 2005, 43, 753–769. [Google Scholar] [CrossRef]
  23. Zhang, B.; Perrie, W.G.; He, Y.J. Wind speed retrieval from RADARSAT-2 quad-polarization images using a new polarization ratio model. J. Geophys. Res. 2011, 116, 1318–1323. [Google Scholar] [CrossRef]
  24. Mouche, A.A.; Chapron, B. Global C-Band Envisat, RADARSAT-2 and Sentinel-1 SAR measurements in co-polarization and cross-polarization. J. Geophys. Res. 2015, 88, 3–9. [Google Scholar]
  25. Lin, B.; Shao, W.Z.; Li, X.F.; Li, H.; Du, X.Q.; Ji, Q.Y.; Cai, L.N. Development and validation of an ocean wave retrieval algorithm for VV-polarization Sentinel-1 SAR data. Acta Oceanol. Sin. 2017, 36, 95–101. [Google Scholar] [CrossRef]
  26. Chapron, B.; Johnsen, H.; Garello, R. Wave and wind retrieval from SAR images of the ocean. In Annales des Telecommunications; Springer: Berlin, Germany, 2001; pp. 682–699. [Google Scholar]
  27. Alpers, W.; Bruning, C. On the relative importance of motion-related contributions to SAR imaging mechanism of ocean surface waves. IEEE Trans. Geosci. Remote Sens. 1986, 24, 873–885. [Google Scholar] [CrossRef]
  28. Hasselmann, K.; Hasselmann, S. On the nonlinear mapping of an ocean wave spectrum into a synthetic aperture radar image spectrum. J. Geophys. Res. 1991, 96, 10713–10729. [Google Scholar] [CrossRef]
  29. Alpers, W.; Ross, D.B.; Rufenach, C.L. On the detectability of ocean surface waves by real and synthetic radar. J. Geophys. Res. 1981, 86, 10529–10546. [Google Scholar] [CrossRef]
  30. Hasselmann, K.; Raney, R.K.; Plant, W.J.; Alpers, W.; Shuchman, R.A.; Lyzenga, D.R. Theory of synthetic aperture radar ocean imaging: A MARSEN view. J. Geophys. Res. 1985, 90, 4659–4686. [Google Scholar] [CrossRef]
  31. Mastenbroek, C.; de Valk, C.F. A semi-parametric algorithm to retrieve ocean wave spectra from synthetic aperture radar. J. Geophys. Res. 2000, 105, 3497–3516. [Google Scholar] [CrossRef]
  32. Sun, J.; Guan, C.L. Parameterized first-guess spectrum method for retrieving directional spectrum of swell-dominated waves and huge waves from SAR images. Chin. J. Oceanol. Limnol. 2006, 24, 12–20. [Google Scholar]
  33. Sun, J.; Kawamura, H. Retrieval of surface wave parameters from SAR images and their validation in the coastal seas around Japan. J. Oceanogr. 2009, 65, 567–577. [Google Scholar] [CrossRef]
  34. Shao, W.Z.; Li, X.F.; Sun, J. Ocean wave parameters retrieval from TerraSAR-X images validated against buoy measurements and model results. Remote Sens. 2015, 7, 12815–12828. [Google Scholar] [CrossRef]
  35. Schulz-Stellenfleth, J.; Lehner, S.; Hoja, D. A parametric scheme for the retrieval of two-dimensional ocean wave spectra from synthetic aperture radar look cross spectra. J. Geophys. Res. 2005, 110, 297–314. [Google Scholar] [CrossRef]
  36. Johnsen, H.; Engen, G.; Hogda, K.A.; Chapron, B.; Desnos, Y.L. Validation of ENVISAT-ASAR wave mode Level 1 and 2 products using ERS SAR data. In Proceedings of the IEEE 2000 International Geoscience and Remote Sensing Symposium, Honolulu, HI, USA, 24–28 July 2000; pp. 1498–1500. [Google Scholar]
  37. Lyzenga, D.R. Unconstrained inversion of wave height spectra from SAR images. IEEE Trans. Geosci. Remote Sens. 2002, 40, 261–270. [Google Scholar] [CrossRef]
  38. Schulz-Stellenfleth, J.; Konig, T.; Lehner, S. An empirical approach for the retrieval of integral ocean wave parameters from synthetic aperture radar data. J. Geophys. Res. 2007, 112, 1–14. [Google Scholar] [CrossRef]
  39. Li, X.M.; Lehner, S.; Bruns, T. Ocean wave integral parameter measurements using Envisat ASAR wave mode data. IEEE Trans. Geosci. Remote Sens. 2011, 49, 155–174. [Google Scholar] [CrossRef] [Green Version]
  40. Pleskachevsky, A.L.; Rosenthal, W.; Lehner, S. Meteo-Marine parameters for highly variable environment in coastal regions from satellite radar images. ISPRS J. Photogramm. Remote Sens. 2016, 119, 464–484. [Google Scholar] [CrossRef]
  41. Ren, L.; Yang, J.S.; Zheng, G.; Wang, J. Significant wave height estimation using azimuth cutoff of C-band RADARSAT-2 single-polarization SAR images. Acta Oceanol. Sin. 2015, 12, 1–9. [Google Scholar] [CrossRef]
  42. Stopa, J.E.; Ardhuin, F.; Collard, F.; Chapron, B. Estimating wave orbital velocities through the azimuth cut-off from space borne satellites. J. Geophys. Res. 2016, 120, 7616–7634. [Google Scholar] [CrossRef]
  43. Shao, W.Z.; Zhang, Z.; Li, X.F.; Li, H. Ocean wave parameters retrieval from Sentinel-1 SAR imagery. Remote Sens. 2016, 8, 707. [Google Scholar] [CrossRef]
  44. Grieco, G.; Lin, W.; Migliaccio, M.; Nirchio, F.; Portabella, M. Dependency of the Sentinel-1 azimuth wavelength cut-off on significant wave height and wind speed. Int. J. Remote Sens. 2016, 37, 5086–5104. [Google Scholar] [CrossRef]
  45. Stopa, J.E.; Mouche, A.A. Significant wave heights from Sentinel-1 SAR: Validation and applications. J. Geophys. Res. 2017, 122, 1827–1848. [Google Scholar] [CrossRef]
  46. Shao, W.Z.; Lin, B.; Li, H.; Ji, Q.Y. Development and validation of wave retrieval algorithms for Sentinel-1 synthetic aperture radar in HH-polarization. J. App. Remote Sens. 2017. under review. [Google Scholar]
  47. Alpers, W.; Brummer, B. Atmospheric boundary layer rolls observed by the synthetic aperture radar aboard the ERS-1 satellite. J. Geophys. Res. 1994, 99, 12613–12621. [Google Scholar] [CrossRef]
  48. Shao, W.Z.; Zhang, Z.; Li, X.M.; Wang, W.L. Sea surface wind speed retrieval from TerraSAR-X HH-polarization data using an improved polarization ratio model. IEEE J. Sel. Top. Appl. Earth Observ. Remote Sens. 2016, 9, 4991–4997. [Google Scholar] [CrossRef]
  49. Monaldo, F.M.; Beal, R.C. Comparison of SIR-C SAR wavenumber spectra with WAM model predictions. J. Geophys. Res. 1998, 103, 18815–18825. [Google Scholar] [CrossRef]
  50. Wang, H.; Yang, J.S.; Mouche, A.A.; Shao, W.Z.; Zhu, J.H.; Ren, L.; Xie, C.H. GF-3 SAR ocean wind retrieval: The first view and preliminary assessment. Remote Sens. 2017, 9, 694. [Google Scholar] [CrossRef]
Figure 1. The quick-look image of calibrated Gaofen-3 (GF-)3 image in vertical–vertical (VV)-polarization around the Hawaiian islands acquired in Quad-Polarization Stripmap (QPS) mode at 16:22 UTC on 20 December 2016.
Figure 1. The quick-look image of calibrated Gaofen-3 (GF-)3 image in vertical–vertical (VV)-polarization around the Hawaiian islands acquired in Quad-Polarization Stripmap (QPS) mode at 16:22 UTC on 20 December 2016.
Sensors 17 01705 g001
Figure 2. (a) European Centre for Medium-Range Weather Forecasts (ECMWF) wind field at 18:00 UTC and; (b) Significant wave height (SWH) field from WaveWatch-III at 18:00 UTC on 20 December 2016. The black rectangle represents the coverage of GF-3 synthetic aperture radar (SAR) image in Figure 1.
Figure 2. (a) European Centre for Medium-Range Weather Forecasts (ECMWF) wind field at 18:00 UTC and; (b) Significant wave height (SWH) field from WaveWatch-III at 18:00 UTC on 20 December 2016. The black rectangle represents the coverage of GF-3 synthetic aperture radar (SAR) image in Figure 1.
Sensors 17 01705 g002
Figure 3. SAR-derived wind field from VV-polarization GF-3 SAR image in Stander Stripmap (SS) mode at 02:17 UTC on 29 September 2016, in which the white circle represents the location of National Data Buoy Center (NDBC) in situ buoy (ID: 46013).
Figure 3. SAR-derived wind field from VV-polarization GF-3 SAR image in Stander Stripmap (SS) mode at 02:17 UTC on 29 September 2016, in which the white circle represents the location of National Data Buoy Center (NDBC) in situ buoy (ID: 46013).
Sensors 17 01705 g003
Figure 4. The comparison between SAR-derived wind speed U10 and measurements from NDBC buoys. (a) Comparison for 16 VV-polarization GF-3 SAR images; and (b) Comparison for 42 horizontal–horizontal (HH)-polarization GF-3 SAR images.
Figure 4. The comparison between SAR-derived wind speed U10 and measurements from NDBC buoys. (a) Comparison for 16 VV-polarization GF-3 SAR images; and (b) Comparison for 42 horizontal–horizontal (HH)-polarization GF-3 SAR images.
Sensors 17 01705 g004
Figure 5. (a) The two-dimensional SAR spectrum of sub-scene covering the buoy (ID:46013) in Figure 3; (b) The corresponding one-dimensional spectrum with the fitted result by using Gaussian fit function.
Figure 5. (a) The two-dimensional SAR spectrum of sub-scene covering the buoy (ID:46013) in Figure 3; (b) The corresponding one-dimensional spectrum with the fitted result by using Gaussian fit function.
Sensors 17 01705 g005
Figure 6. The comparison between SAR-derived Hs and the measurements from NDBC buoys. (a) Comparison for 16 VV-polarization GF-3 SAR images; and (b) Comparison for 42 HH-polarization GF-3 SAR images.
Figure 6. The comparison between SAR-derived Hs and the measurements from NDBC buoys. (a) Comparison for 16 VV-polarization GF-3 SAR images; and (b) Comparison for 42 HH-polarization GF-3 SAR images.
Sensors 17 01705 g006
Figure 7. The comparison between SAR-derived wind speed U10 and ECMWF re-analysis gridded winds at 1 m/s bins. (a) Comparison for 96 VV-polarization GF-3 SAR images; and (b) Comparison for 70 HH-polarization GF-3 SAR images.
Figure 7. The comparison between SAR-derived wind speed U10 and ECMWF re-analysis gridded winds at 1 m/s bins. (a) Comparison for 96 VV-polarization GF-3 SAR images; and (b) Comparison for 70 HH-polarization GF-3 SAR images.
Sensors 17 01705 g007
Figure 8. The comparison between SAR-derived SWH Hs and WaveWatch-III data at 0.5 m/s bins. (a) Comparison of 96 VV-polarization GF-3 SAR images; and (b) Comparison of 70 HH-polarization GF-3 SAR images.
Figure 8. The comparison between SAR-derived SWH Hs and WaveWatch-III data at 0.5 m/s bins. (a) Comparison of 96 VV-polarization GF-3 SAR images; and (b) Comparison of 70 HH-polarization GF-3 SAR images.
Sensors 17 01705 g008

Share and Cite

MDPI and ACS Style

Shao, W.; Sheng, Y.; Sun, J. Preliminary Assessment of Wind and Wave Retrieval from Chinese Gaofen-3 SAR Imagery. Sensors 2017, 17, 1705. https://doi.org/10.3390/s17081705

AMA Style

Shao W, Sheng Y, Sun J. Preliminary Assessment of Wind and Wave Retrieval from Chinese Gaofen-3 SAR Imagery. Sensors. 2017; 17(8):1705. https://doi.org/10.3390/s17081705

Chicago/Turabian Style

Shao, Weizeng, Yexin Sheng, and Jian Sun. 2017. "Preliminary Assessment of Wind and Wave Retrieval from Chinese Gaofen-3 SAR Imagery" Sensors 17, no. 8: 1705. https://doi.org/10.3390/s17081705

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop