1. Introduction
The sea surface wind field is one of the important ocean wave dynamic processes, and has a very significant impact on the generation of ocean waves and the movement of ocean currents [
1], as well as a very significant physical significance in the field of marine meteorology. In addition, the observation of the sea surface wind field is of great significance to maritime navigation, military defense, and other aspects. Therefore, the observation of the sea surface wind field is extremely significant.
The observation of the sea surface wind field is originally conducted by using meteorological stations, anemometer towers, buoys, ships and other on-site observation methods. These methods can obtain accurate sea surface wind field information, but there are some problems: (1) The range of information obtained by meteorological stations, anemometer towers, buoys and other methods is very limited. Generally, only offshore wind field information can be obtained. These methods cannot obtain large-area wind field information; (2) The maintenance cost of equipment such as anemometer towers and buoys is high, and widespread use of equipment is difficult; (3) Ship observation takes a long time and the data obtained are only distributed on the navigation track. Therefore, traditional methods such as meteorological stations, anemometer towers, and buoys cannot meet the needs of a large-scale real-time observation of the sea surface wind field [
2].
With the development of satellite remote sensing observation technology, it is possible to observe the wind field in a large area. Satellite remote sensing technology mainly includes a microwave scatterometer, radiometer and synthetic aperture radar. The scatterometer, altimeter and radiometer are first applied to the observation of the sea surface wind field, which can obtain a large area of wind field data, though their resolution is not high. Synthetic aperture radar (SAR) has a high resolution, which is not affected by light and weather conditions, and can realize all day, all-weather observation. In June 1978, the United States launched the first integrated ocean power satellite Seasat-A, which was equipped with an altimeter, a scatterometer and a synthetic aperture radar, and acquired a lot of accurate and effective information. In view of the great success of the US SAR satellite launch, the European Space Agency (ESA), the German Aerospace Center, and the Canadian Space Agency also launched remote sensing satellites equipped with synthetic aperture radars ERS-1, ERS-2, Sentinel-1, TerraSAR-X, RADARSAT-1, and RADARSAT-2, etc. China also launched the C-band synthetic aperture radar satellite GF-3. The use of a synthetic aperture radar to observe the sea surface wind field has now become the most important method.
The sea surface wind field inversion using a synthetic aperture radar is generally achieved using geophysical model functions (GMF). The geophysical model function is generally used for the sea surface wind field inversion of VV (vertical transmit and vertical receive) polarization SAR data. The most commonly used GMFs include: CMOD4 [
3], CMOD-IFR2 [
4], CMOD5 [
5], and CMOD5.N [
6]. In 1997, Stoffelen et al. [
3] developed the CMOD4 function based on the launched ERS satellite scatterometer data and the European Centre for Medium-range Weather Forecasts (ECMWF) numerical simulation data. In 1998, Quilfen et al. [
4] developed the CMOD-IFR2 function based on the US National Oceanic and Atmospheric Administration (NOAA) buoy data, ERS satellite scatterometer data and ECMWF numerical simulation data. Addressing the problem of low wind speed in CMOD4 and high wind speeds (greater than 25 m/s) in CMOD-IFR2 [
7,
8], Hersbach et al. (2007) [
5] improved the CMOD4 function and developed the CMOD5 function on the basis of this. In 2010, Hersbach [
6] adjusted the CMOD5 parameters based on ECMWF data under neutral atmospheric conditions and developed the CMOD5.N function.
For the wind field inversion of HH (horizontal transmit and horizontal receive) polarization SAR data, the polarization ratio model was usually used to convert the HH polarization backscatter coefficient to the VV polarization backscatter coefficient, and then the CMOD series of functions were used to perform the wind field inversion. In 1998, Thompson et al. [
9] developed a polarization ratio model based on the HH polarization data of RADARSAT-1, as follows:
where PR represents the polarization ratio;
represents the radar incident angle; and
is an empirical parameter, ranging from 0 to 2. When
= 0, it is the Bragg model; when
= 0.6, it is called the Thompson model [
9]; and when
= 1, it is called the Kirchhoff model. In 2005, Mouche et al. [
10] proposed a new polarization ratio model based on the C-band airborne real-aperture radar, which only depends on the radar incident angle, as shown below:
where PR is the polarization ratio;
is the radar incidence angle; and
,
and
are the calibration parameters of the model.
Some scholars have found that the co-polarization backscatter coefficient will show the saturation effect under high wind speed conditions, while the cross-polarization backscatter coefficient will not show the saturation effect under high wind speed conditions [
11]. In 2012, Zhang and Perrie [
12] developed a cross-polarization inversion model based on the RADARSAT-2 full-polarization data and named it C-2PO. In 2017, Zhang et al. [
13] found that the cross-polarization backscatter coefficient is related to the radar incidence angle, and proposed the C-3PO model on the basis of C-2PO.
The Sentinel-1 satellite launched by ESA and the GF-3 satellite launched by China are SAR imaging satellites operating in orbit and are widely used in marine research. The research aim of this paper is to use the SAR data of Sentinel-1 and GF-3 satellites to perform wind field inversion using CMOD5.N. The inverted wind speed is compared with the reanalysis data ERA5 and buoy data to compare the advantages and disadvantages of the inverted wind speed of the two satellites, and to compare the quality of Sentinel-1 SAR data and GF-3 SAR data.
The paper is organized as follows. We describe our materials and methods in
Section 2, including imaging area, SAR data reanalysis data, buoy data, spatio-temporal matching and wind field inversion method. In
Section 3, we present our results and discussion. The conclusions are provided in
Section 4.
3. Results and Discussion
3.1. Comparison of Inverted Wind Field and ERA5 Data
This study uses SAR data with imaging areas located in the Hormuz Strait, Malacca Strait, and the east and west coasts of the United States. Sentinel-1 SAR data have 95 scenes, and GF-3 SAR data have 57 scenes. The wind speed provided by ERA5 is divided into three intervals of 0~5 m/s, 5~10 m/s, and above 10 m/s for comparison. The root mean squared error (RMSE) and the scatter index (SI) are used for evaluation. The advantages and disadvantages of the wind field inverted by Sentinel-1 SAR data and GF-3 SAR data under different wind speeds are compared to evaluate the quality of the two satellites. The RMSE and the SI are defined as follows:
where
represents the wind speed inverted by SAR;
represents the reference wind speed, i.e., ERA5 wind speed; and
and
represent the average value of all
and
.
In this study, the backscatter coefficient of each scene is divided into 25 × 25 blocks, and each block represents a sub image. The backscatter coefficient of each sub image takes the average value as the backscattering coefficient of the sub image. The central latitude, longitude and incident angle of each sub image are used as the latitude, longitude and incident angle of the sub image. Each sub image is spatio-temporal matched with ERA5 data. The wind direction of ERA5 is used as input, and CMOD5. N is used for wind field inversion. The wind speed, incident angle and relative wind direction are input into CMOD5.N, the wind speed is traversed from 0 to 30 m/s at an interval of 1 m/s, and a set of backscatter coefficients are obtained. The wind speed corresponding to the backscatter coefficient closest to the backscatter coefficient of the sub image is the inverted wind speed. One scene of SAR data can be inverted to obtain 25 × 25 wind speeds, and the inverted wind speed is obtained after removing invalid data. The inverted wind speed is compared with the matched ERA5 wind speed, the RMSE and the SI are calculated, and the quality of Sentinel-1 SAR data and GF-3 SAR data applied to wind field inversion are compared.
3.1.1. Comparison of Inversion Results When Wind Speed of 0~5 m/s
When the wind speed of ERA5 is 0~5 m/s, the comparison result of the inversion wind speed and ECMWF data is shown in
Figure 5.
When the wind speed of ERA5 is 0~5 m/s, it can be seen from the above figure that compared with ECMWF data, the RMSE of Sentinel-1 SAR data is 1.66 m/s, and the SI is 0.29; the RMSE of GF-3 SAR data is 1.63 m/s, and the SI is 0.36.
When the wind speed of ERA5 is 0~5 m/s, the wind field inversion results of GF-3 SAR data are slightly better than the inversion results of Sentinel-1 SAR data, and the wind field inversion accuracy of Sentinel-1 SAR data and GF-3 SAR data meet the recognized index requirements, that is, the RMSE is less than 2 m/s. When the wind speed is in the range of 0~5 m/s, the RMSE of the two satellite inversions is large, and the inversion effect is poor. The possible reason for this is that when the wind speed is in the range of 0~5 m/s, the wind speed is low, the sea surface roughness is small, and the backscatter coefficient is not very accurate.
3.1.2. Comparison of Inversion Results When Wind Speed of 5~10 m/s
When the wind speed of ERA5 is 5~10 m/s, the comparison result of the inversion wind speed and ECMWF data is shown in
Figure 6.
When the wind speed of ERA5 is 5~10 m/s, it can be seen from the above figure that, compared with ECMWF data, the RMSE of Sentinel-1 SAR data is 1.37 m/s, and the SI is 0.19; the RMSE of GF-3 SAR data is 1.45 m/s, and the SI is 0.21.
It can be seen that when the wind speed of ERA5 is 5~10 m/s, the wind field inversion result of Sentinel-1 SAR data is slightly better than the inversion results of GF-3 SAR data. The wind field inversion accuracy of Sentinel-1 SAR data and GF-3 SAR data meets the recognized index requirements, that is, the RMSE is less than 2 m/s.
3.1.3. Comparison of Inversion Results When Wind Speed above 10 m/s
When the wind speed of ERA5 is above 10 m/s, the comparison result of the inversion wind speed and ECMWF data is shown in
Figure 7.
When the wind speed of ERA5 is above 10 m/s, it can be seen from the above figure that, compared with ECMWF data, the RMSE of Sentinel-1 SAR data is 1.49 m/s, and the SI is 0.12; the RMSE of GF-3 SAR data is 1.87 m/s, and the SI is 0.14.
When the wind speed of ERA5 is above 10 m/s, the wind field inversion results of Sentinel-1 SAR data are better than the inversion results of GF-3 SAR data, and the wind field inversion accuracy of Sentinel-1 SAR data and GF-3 SAR data meets the recognized index requirements, that is, the RMSE is less than 2 m/s. The RMSE of the inversion results of GF-3 SAR data is large, which may be due to the small number of GF-3 SAR data with high wind speed.
In summary, in the three wind speed ranges of 0~5 m/s, 5~10 m/s and above 10 m/s, the RMSE of Sentinel-1 SAR data is 1.66 m/s, 1.37 m/s, and 1.49 m/s, respectively, and the SI of Sentinel-1 is 0.29, 0.19 and 0.12 respectively; the RMSE of GF-3 SAR data is 1.63 m/s, 1.45 m/s, and 1.87 m/s, respectively, and the SI of GF-3 is 0.36, 0.21 and 0.14, respectively. The RMSE of the inversion results of Sentinel-1 SAR data and GF-3 SAR data in three wind speed ranges is less than 2 m/s, and the inversion accuracy is high. In general, under different wind speed conditions, the wind field inversion results of Sentinel-1 SAR data are slightly better than the wind field inversion results of GF-3 SAR data. When the two satellites are applied to wind field inversion, the SAR data quality of the Sentinel-1 satellite is slightly better than the SAR data quality of the GF-3 satellite.
3.2. Comparison of Inverted Wind Field and Buoy Data
Buoy observation is recognized as the most accurate observation method, and the buoy data are measured in reliable and accurate on-site data. In this study, we use ERA5 wind direction data matched with SAR data as the external input wind direction, and use CMOD5.N to invert the wind speed of Sentinel-1 and GF-3 SAR data, and compare the inverted wind speed with buoy data to further verify the quality of Sentinel-1 and GF-3 SAR data.
This study uses SAR data with imaging areas located near the east and west coasts of the United States, and uses RMSE for evaluation. A total of 20 scenes of Sentinel-1 SAR data and 10 scenes of GF-3 SAR data are used for wind speed inversion. The results of wind speed inversion and buoy data comparison are shown in
Figure 8.
As can be seen from
Figure 8, compared with the buoy data, the RMSE of Sentinel-1 SAR data is 1.20 m/s, and the RMSE of GF-3 SAR data is 1.48 m/s. The RMSE of wind field inversion results from Sentinel-1 and GF-3 SAR data is less than 2 m/s, and the inversion accuracy is high. The inversion results of Sentinel-1 SAR data are slightly better than the inversion results of GF-3 SAR data.
3.3. Quality Evaluation
The RMSE of Sentinel-1 SAR data in the three wind speed ranges of 0~5 m/s, 5~10 m/s, and above 10 m/s is less than 2 m/s. The RMSE of 0~5 m/s is the largest, at 1.66 m/s; the RMSE of 5~10 m/s is the smallest, at 1.37 m/s. The SI is greater than 0.2 in the range of 0~5 m/s, which is 0.29. The SI of the other two ranges is less than 0.2. The inversion results of Sentinel-1 SAR data are poor in low wind speeds, and the inversion results are better when the wind speed is greater than 5 m/s. The RMSE of the inversion results of GF-3 SAR data in the three wind speed ranges of 0~5 m/s, 5~10 m/s, and above 10 m/s is less than 2 m/s. The RMSE in the wind speed range of above 10 m/s is the largest, at 1.87 m/s. The RMSE in the range of 5~10 m/s is the greatest, at 1.45 m/s. The maximum SI is 0.36 in the range of 0~5 m/s, and the minimum is 0.14 in the range of above 10 m/s. In general, under different wind speeds, the wind field inverted from Sentinel-1 SAR data is slightly better than the wind field inverted from GF-3 SAR data.
Compared with the buoy data, the RMSE of Sentinel-1 SAR data is 1.20 m/s (similar to the accuracy in [
14]), and the RMSE of GF-3 SAR data is 1.48 m/s (as in [
17], RMSE is less than 2 m/s), and both are less than 2 m/s. The accuracy is high, but the inversion result of Sentinel-1 SAR data is slightly better than GF-3 SAR data.
The RMSE of Sentinel-1 SAR data and GF-3 SAR data under different wind speed conditions is less than 2 m/s, and both have a high inversion accuracy, but the inversion results of Sentinel-1 SAR data are slightly better than GF-3 SAR data. The comparison with the buoy data shows that the inversion accuracy of Sentinel-1 and GF-3 SAR data is higher, but the inversion results of Sentinel-1 SAR data are slightly better than GF-3 SAR data. In summary, for wind field inversion, the data quality of Sentinel-1 SAR and GF-3 SAR are both high, and the data quality of Sentinel-1 SAR is slightly better than the quality of GF-3 SAR.