Next Article in Journal
Tree-Structured Parzan Estimator–Machine Learning–Ordinary Kriging: An Integration Method for Soil Ammonia Spatial Prediction in the Typical Cropland of Chinese Yellow River Delta with Sentinel-2 Remote Sensing Image and Air Quality Data
Next Article in Special Issue
Study on the Polarization Pattern Induced by Wavy Water Surfaces
Previous Article in Journal
High Spatial Resolution Fractional Vegetation Coverage Inversion Based on UAV and Sentinel-2 Data: A Case Study of Alpine Grassland
Previous Article in Special Issue
Toward a High-Resolution Wave Forecasting System for the Changjiang River Estuary
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A New Approach for Ocean Surface Wind Speed Retrieval Using Sentinel-1 Dual-Polarized Imagery

1
Faculty of Information Science and Engineering, Ocean University of China, Qingdao 266100, China
2
Hainan Institute, Zhejiang University, Sanya 572025, China
*
Author to whom correspondence should be addressed.
Remote Sens. 2023, 15(17), 4267; https://doi.org/10.3390/rs15174267
Submission received: 11 July 2023 / Revised: 25 August 2023 / Accepted: 29 August 2023 / Published: 30 August 2023

Abstract

:
A synthetic aperture radar (SAR) has the capability to observe ocean surface winds with a high spatial resolution, even under extreme conditions. The purpose of this work was to develop a new method for wind speed retrieval with the combination of SAR dual-polarized signals. In this study, we collected 28 tropical cyclone imageries observed using the Sentinel-1 dual-polarization mode. These imageries were collocated with radiometer wind speed measurements and reanalysis of wind vector products. In the new method, the wind speed was set as the output. VV-polarized (vertical transmitting–vertical receiving polarized) normalized radar cross section (NRCS), incident angle, VH-polarized (vertical transmitting–horizontal receiving polarized) NRCS, and wind direction were set as the inputs. Based on different output combinations, wind retrieval models were developed with multiple linear regression (MLR). According to the validation and comparison, the proposed models performed better than the traditional piecewise VH-polarization geophysical model functions (GMFs). The impact of thermal noise on the retrieval of low wind speeds (<10 m/s) could be partially reduced. The input of wind direction is unnecessary if the combination of VV- and VH-polarized imageries has been utilized. These results suggest that the use of MLR and the dual-polarization combination can improve SAR wind retrieval accuracy. Compared with SMAP measurements, our SAR retrievals can provide fine structures of TC wind fields.

1. Introduction

The spaceborne synthetic aperture radar (SAR) has the capability to acquire backscatter from the ocean surface in all-weather conditions [1,2]. This imagery has been widely used to detect ocean dynamic phenomena, such as wind vectors [3,4], surface waves [5,6,7,8], surface currents [9,10], and internal waves [11,12]. In general, ocean surface roughness is positively dependent upon wind speed. SAR sensors can receive more Bragg or non-Bragg backscatter when surface wind speed increases. Based on this theory, wind speeds could be retrieved from SAR imagery [13]. With the advantages of high resolution, large coverage, and high accuracy, SAR wind retrievals have contributed to oceanic and atmospheric investigations, for example, the analysis of sea surface wind distribution [14,15], the estimation of tropical cyclone (TC) intensity, size and fullness [16,17], and the initialization of numerical models [18,19].
C-band SAR instruments in orbit, such as RADARSAT-2, Sentinel-1, and Gaofen-3, generally have a dual-polarization sensing mode. This model has the advantage of receiving signals in both co- (i.e., VV- and HH-polarization) and cross-polarization (i.e., VH- and HV-polarization), which contain abundant backscattering information to observe the ocean surface more effectively. The geophysical model function (GMF), which relates the microwave normalized radar cross section (NRCS) to the wind speed, has played a key role in the wind retrieval from co-polarization and cross-polarization channels of SAR. To date, many empirical functions have been proposed. For co-polarized data, there is the CMOD family, including CMOD4 [20], CMOD5.N [21], CMOD7 [22], etc., originally developed for scatterometers, together with the C_SARMOD [23] and C_SARMOD2 [24] models re-tuned for SAR. For cross-polarized SAR imagery, there are C-2PO [2], C-2POD [25], C-3PO [26], and MS1A [27], and a large number of modified models, such as S1IW.NR [28], S1ER.NR [29] and MMS1A [17]. For most GMFs, the sigma nought ( σ 0 ) is on the left side of the equation. The right side is a function of wind speed, wind direction, and the incident angle. In wind retrieval, people usually have to pre-set a group of wind speed values and then input them into a GMF to obtain a first-guess σ 0 , and finally match them with SAR-measured NRCS to determine wind speeds. The difficulty of inversing GMFs makes wind retrieval less direct and leads to greater computational complexity.
According to the relationships between NRCS and wind speed, the most significant difference between co- and cross-polarization GMFs is the wind retrieval range. Since co-polarized signals saturate under high wind conditions (>25 m/s), the imagery cannot recognize the high wind area; therefore, they are often used to retrieve low-to-moderate wind speeds (<25 m/s) [29]. Meanwhile, cross-polarized acquisitions keep high sensitivity with high winds but are not good at low-wind retrieval due to the strong impact of thermal noise compared with conventional co-polarized acquisitions [27,29,30]. In addition, the wind direction parameter is another difference between co- and cross-polarization GMFs. Many studies have indicated that cross-polarized NRCS is almost insensitive to wind direction and, therefore, most cross-polarization GMFs do not require a wind direction input, making retrieval easier [31,32]. However, co-polarized NRCS is sensitive to wind direction since wind direction can change the local incident angle to influence Bragg contribution; thus, the GMFs for co-polarized images generally need a wind direction input [20,22,33].
With dual-polarized imagery, the differing sensitivity between contemporaneous co- and cross-polarized SAR signals can be advantageously combined to retrieve winds ranging from low to high as well as wind directions. Zhang et al. employed CMOD5.N and C-2POD to construct a cost function to merge information from both VV and VH channels [25]. The minimization of the cost function allows optimum estimates for wind vectors. This method performs well for the RADARSAT-2 imagery of TCs according to the comparison of retrievals with collocated QuikSCAT scatterometer winds. Mouche et al. also developed a cost function for the determination of the maximum probability to obtain a wind vector based on CMOD5.N, MS1A, and the priori information of wind vectors from ECMWF (spatial resolution is 0.125° with a time step of 3 h) [27]. The combination model performed better than the single GMF in two aspects. On the one hand, compared with the single use of the CMOD5.N model, the combination method improved the accuracy of SAR retrieval for winds larger than 25–30 m/s using the cross-polarization channel of the Sentinel-1 C-SAR sensor. On the other hand, despite noise equivalent sigma zero (NESZ) correction, the impact of NESZ remains significant at sub-swath boundaries for low wind regions. The combination model can solve this issue with the contribution of CMOD5.N and ECMWF references. Overall, excellent consistency was found between the Soil Moisture Active Passive (SMAP) winds and the results of their method, indicating that the joint use of co- and cross-polarization can indeed improve wind retrieval under TC conditions.
In this paper, a new approach was explored to retrieve ocean surface wind speed with the combination of co- and cross-polarized imageries, inspired by the well-designed empirical wave retrieval model family CWAVE [34,35,36], instead of a cost function. We focused on the extra wide swath (EW) mode and the interferometric wide swath (IW) mode of the Sentinel-1 A/B sensors and collected 28 TC cases with multiple wind collocations for model development and validation. Different schemes were investigated for inputs. Multiple linear regression (MLR) was utilized in modeling. A comparison was made between the proposed models and some traditional models. In addition, wind retrieval bias was assessed against the rain rate. The remaining sections of this paper are organized as follows. Section 2 introduces the Sentinel-1 imagery, wind references, and data collocation. Section 3 presents a series of functions for the new model. In Section 4, the new models are validated with TC wind measurements and analyses. Discussion and conclusions are made in Section 5 and Section 6, respectively.

2. Dataset Description

2.1. Sentinel-1 SAR Imagery

The C-SAR sensor on board the European Space Agency (ESA) Sentinel-1 A/B satellites can be operated in four different modes: the Stripmap (SM) mode, the Interferometric Wide swath (IW) mode, the Extra Wide Swath (EW) mode and the Wave (WV) mode. For TC observations, the EW mode and the IW mode are generally used to acquire extensive wind field information with the largest and second-largest swaths in four modes. Some imagery parameters are shown in Table 1. Except for these parameters, EW and IW imagery have different NESZ distributions, leading to different NRCS–wind relationships when surface wind is weak, which should be particularly considered in model development. In this study, 16 and 12 TC cases were, respectively, observed in EW and IW modes from 2016 to 2021, with dual-polarization (VV + VH). The TCs were over the Atlantic Ocean (AL), the East Pacific Ocean (EP), the Southern Hemisphere Ocean (SH), and the West Pacific Ocean (WP). Sentinel Application Platform (SNAP) 7.0 software was utilized for GRD border noise removal, thermal noise removal, and imagery calibration.

2.2. SFMR Wind and Rainfall Measurements

The National Oceanic and Atmospheric Administration (NOAA) Stepped-Frequency Microwave Radiometer (SFMR) on board reconnaissance aircraft can measure the brightness temperature from the sea surface at six C-band frequencies along the flight track over a storm [37]. Ocean surface wind speeds are retrieved according to a function between the wind speed and surface emissivity. Rain rate retrieval is based on the relationship between microwave absorption and rain rate. The SFMR wind speed measurements have an RMSE of ~3.9 m/s compared with the dropsonde measurements [38]. SFMR wind speeds and rain rate with a spatial resolution of 0.01° have been used as collocations in the establishment and evaluation of the SAR wind and rainfall retrieval models [39,40]. Our study used data for model development.

2.3. SMAP Wind Measurements

The L-band radiometer on board the SMAP satellite can measure the bright temperature to retrieve the ocean surface wind speed without being affected by rainfall [41,42]. Its 1000 km swath provides it with a chance to make a complete observation of a storm. Remote Sensing Systems (www.remss.com/missions/smap/ (accessed on 21 March 2022)) provided the SMAP Final Wind Speed products and considered the equivalent neutral winds at a 10 m height above the ocean surface. Their grid spacing was 0.25° × 0.25°. For winds greater than 25 m/s, a consistency was found between SMAP and SFMR measurements, with a bias of 0.5 m/s and a standard deviation of ~3 m/s [38]. In addition, the evaluation of size values proved that the SMAP winds could be used for TC wind field monitoring and modeling [43]. In this paper, we used the SMAP Final Wind Speed products as an independent data source for validation.

2.4. Dropsonde Wind Vector Measurements

To validate the proposed wind retrieval models, we collected NOAA Hurricane Research Division (HRD) dropsonde data. Dropsondes were released from the aircraft along the flight track and drifted down on a parachute, measuring vertical profiles of pressure, temperature, humidity, and wind as they fell. At a 10 m height, the accuracy of the wind speed measurement was about 0.5–2.0 m/s [44]. In this work, only sea surface wind vector measurements were used whose pressure level of observation was 1070.0 hPa. They acquired over 8 TCs. There were 11 observations collocated with EW mode SAR images and 12 observations collocated with IW mode SAR images. The time difference between dropsonde observations and SAR images was less than 0.5 h.

2.5. ECMWF ERA5 Wind Vectors

ERA5 is the fifth-generation ECMWF atmospheric reanalysis of the global climate, covering the period from 1950 to the present. Based on data assimilation, ERA5 combines model data with observations from across the world into a global dataset with an hourly temporal resolution and a spatial resolution of 0.25° × 0.25° [45]. ECMWF wind vectors have been used as references in investigations of SAR wind retrieval with both signal and dual polarizations. In this study, all wind direction values were collected from ERA5, except a case study on Tropical Storm Karl (2016). In addition, most wind speed values lower than 15 m/s were from ERA5 in model development as a supplement to SFMR measurements.

2.6. H*Wind Data

H*Wind is an integrated system that provides a view of the strength and extent of hurricane wind fields according to the input of all available surface weather observations with respect to the storm and is used operationally to improve the damage assessment of hurricane intensity [46]. In general, the total uncertainty of the H*Wind products in hurricanes is approximately 6% near the storm center, increasing to nearly 13% near the tropical storm force wind radius [47]. In this paper, a case study of Tropical Storm Karl was carried out with the wind vectors provided by Risk Management Solutions H*Wind. The spatial resolution was 0.01° × 0.01°.

2.7. Collocation

To construct the dataset, the Sentinel-1 imagery of TCs was collocated with the wind data introduced above. There were 13 SFMR and 16 SMAP collocations. All SAR imageries were matched with ERA5. The time differences between SAR and SFMR were controlled within 0.5 h. The time difference between SAR and SMAP was just a few minutes. The time difference between SAR and ERA5 was within 1 h. According to the best tracks provided by the National Hurricane Center (NHC) and the Joint Typhoon Warning Center (JTWC), the TC intensities of different matches are shown with the Saffir–Simpson hurricane wind scale in Figure 1. For Tropical Storm Karl, the time difference between SAR and H*Wind was about 1 h and 37 min. According to the TC motion vector estimated from the best track and the time difference, the locations of H*Wind were shifted to match SAR imagery using the method presented in [28].
To note, a SAR image of Hurricane Michael observed at 23:44 UTC on 9 October 2018 could be matched both with SFMR and SMAP (Figure 2), where the storm eye was obvious. Figure 2a,b are VV- and VH-polarized images from the EW mode, respectively. The white and black lines crossing the storm eye are the tracks of flights belonging to the Air Force Reconnaissance (AFRC) and NOAA, respectively. Figure 2c is the collocated SMAP final wind with a maximum value of 53.69 m/s. Figure 3 shows an example of IW mode dual-polarized imagery acquired over Hurricane Irma at 9:07 UTC on 29 October 2017. The white line is the AFRC flight track.

3. Model Development

After data collocation, the distribution of NRCS was compared with the incident angle and wind speed (Figure 4). The winds of the unfilled points were measured using SFMR. Since under TC conditions, most of these winds were larger than 15 m/s, the low winds of the filled points were collected from ERA5 reanalysis as a supplement for data fitting in an entire wind speed range. Figure 4a–d shows the VV- and VH-polarization for the EW and IW modes, respectively. Low winds have a strong linear correlation with VV NRCS. The range of VV NRCS is wide enough to distinguish specific wind speeds. However, for the VH channel, the signals of low winds are seriously affected by NESZ (shown with black curves in Figure 4b,d), making the retrieval difficult. These are the main reasons why VV-polarized imagery performed well in the retrieval of low wind speed, but VH-polarized imagery did not. Interestingly, the opposite happened in moderate and high winds. The linear correlation between NRCS and wind speed was stronger in VH imagery than in VV imagery and weakly impacted by the incident angle and NESZ. So, if signals from VV- and VH-polarization channels are combined for model development, the new model could take advantage of both channels.
Inspired by the CWAVE family of the SAR empirical wave retrieval models [34,35,36], we present a new approach, i.e., Equation (1), to establish the model. Sentinel-1 dual-polarized signals were combined as fitted data. The output is the wind speed at 10 m height ( U 10 ). There are four inputs: VH NRCS ( σ 0 VH ), incident angle ( θ ), VV NRCS ( σ 0 VV ), and wind direction ( φ ). To note, the wind direction is the angle between the real and azimuth directions. We used MLR to automatically decide which parameter was important for wind retrieval and expressed the relationship by determining coefficients. Three different schemes were tested for model inputs. We named them Model 1 (inputs: σ 0 VH , θ ), Model 2 (inputs: σ 0 VH , θ , σ 0 VV ) and Model 3 (inputs: σ 0 VH , θ , σ 0 VV , φ ).
U 10 = i , j = 1 n A i , j X i X j + i = 1 n A i X i + A 0 ,
where A 0 , A i and A i , j are coefficients. For their values, please refer to Table A1. X i and X j are inputs. ji. X 1 = σ 0 VH . X 2 = θ . X 3 = σ 0 VV . X 4 = φ . n = 2, 3, 4 for Model 1, 2, and 3, respectively. Wind retrieval models were separately developed for the EW mode and the IW mode. Then, for evaluation, the new models were used to simulate wind speeds from the fitted data. Figure 5 and Figure 6 show the comparison of the simulations and wind references for two modes.
After MLR fitting, an additional linear regression was utilized to correct the bias, especially for the retrieval of high winds. The final retrieval was U 10 final = a U 10 b . a and b are coefficients. The values of a and b are listed in Table A2. Finally, the new models were evaluated again through an error calculation based on the fitted data, including the bias (retrievals minus fitted data), Pearson correlation coefficient (Cor), and root mean squared error (RMSE). The calculation methods of Cor and RMSE are given in Equations (2) and (3). The results are illustrated in Table 2. In general, the error reduced if a proposed model had more inputs, indicating that the combination of VV- and VH-polarized signals, together with wind direction, have the potential to improve retrieval accuracy.
Cor = cov U 10 final ,       U 10 reference std U 10 final std U 10 reference ,
RMSE = i = 1 N U 10 , i final U 10 , i reference 2 N
where U 10 final is the SAR-simulated wind speed and U 10 reference is the reference wind speed. cov U 10 final ,   U 10 reference stands for the covariance of them. std U 10 final and std U 10 reference are the standard deviations of U 10 final and U 10 reference .

4. Validation and Comparison

4.1. Validation with SMAP Winds

To evaluate the accuracy of the proposed models, validation was carried out in this section. Wind speeds were retrieved via different models from 16 TC imageries observed using the Sentinel-1 EW and IW modes. Then, retrievals were compared with the collocated SMAP wind speeds, ensuring the objectivity of validation. ERA5 wind directions were used as inputs in EW Model 3 and IW Model 3. The sample number distributions of SMAP wind speeds and ERA5 wind directions are shown in Figure 7. The maximum wind was up to about 70 m/s. Wind directions were distributed from 0 to 360 degrees.
Models were validated separately for the Sentienl-1 EW and IW modes. For the EW mode, wind speeds were retrieved by EW Model 1, 2, 3, MMS1A, and S1EW.NR. MMS1A and S1EW.NR are piecewise functions. The MMS1A model was modified from MS1A through a correction of SFMR data [17]. The S1EW.NR model was fitted from ERA5 wind speeds [29]. For wind retrieval, the two models just required VH NRCS and incident angle inputs.
Comparisons are shown in Figure 8. According to the scatters, error statistics were carried out for winds lower than 10 m/s, higher than 10 m/s, and for all winds. The results are illustrated in Table 3. Overall, EW Model 2 performed best. Better consistencies could be found between the results of the proposed models and SMAP winds. The main reason was the modeling method: MLR was able to develop a single function that could make retrievals change more smoothly in adjacent sub-swaths of SAR imagery. However, this is difficult to realize with a piecewise function. In addition, S1EW.NR has a more serious underestimation. Compared with the S1EW.NR model, the MS1A model, was closer to the new models. It showed another factor that affects model performance: the fitted data. To note, all the tested models underestimated winds higher than 30 m/s. This issue was induced by the high wind difference between SFMR, ERA5, and SMAP.
For the Sentienl-1 IW mode, wind speeds were retrieved via IW Model 1, 2, 3, MMS1A, and S1IW.NR. The S1IW.NR model is linear and has a power piecewise function, which is established with SFMR and SMAP data [28]. Figure 9 shows the comparisons of the retrievals and SMAP winds. According to the scatters, the performance metrics are illustrated in Table 4. Just like the EW situation, when compared with the two traditional models, the proposed models had smaller RMSE values. The retrievals of high winds were smaller than SMAP data. Their reasons are the same as the EW situation. In addition, since MMS1A was mainly designed for the EW mode, it had the largest error for winds lower than 10 m/s, indicating that one GMF is unable to fit low winds for both EW and IW modes. This is because the low VH NRCS is seriously affected by different thermal noises for different modes.

4.2. Validation with Dropsonde Winds

An experiment was carried out to validate the proposed models with the winds measured using Dropsonde. The wind directions used in Model 3 were acquired from Dropsonde wind vectors. Comparisons are shown in Figure 10. The EW and IW Models 1, 2, and 3 are the proposed models. The MMS1A, S1EW.NR and S1IW.NR are traditional models. Table 5 and Table 6 show the validation results for the EW and IW modes, respectively. On average, the new models performed better than the traditional models. However, due to a small number of samples, this improvement was not obvious from Models 1 to 3.

4.3. A Case Study of Hurricane Michael

In this study, Hurricane Michael (2018) is very special. This is because it is the only sample whose SAR image was collocated with both SFMR and SMAP measurements, making it possible to compare our retrievals with winds and different grid spacings. Based on EW Model 3, the surface winds of Hurricane Michael were retrieved, as shown in Figure 11a. It had a grid spacing of 1 km. The storm eye, eyewall, and rain bands could be identified easily. This is the advantage of our model. It can provide wind fields with small structures. Figure 11b is the collocated SMAP wind map. Due to a 0.25°-grid-spacing, it looks much coarser. There are no storm eyes and other fine structures.
In Figure 11a, the white line crossing the storm eye represents the AFRC flight track, on which the winds were measured by SFMR. The time difference between the SFMR measurements and the SAR image was less than 0.5 h. To further show the difference between the three wind products, winds along the track were compared in Figure 11b. SFMR and SAR retrievals had the same grid spacing; however, SFMR wind values were slightly lower than SAR. They could observe very low wind speeds in the storm center and radial wind profile. In the region from the eye wall to the outer region, although the distribution of SMAP data was similar to SAR and SFMR, it could not show low winds in the center and the location of maximum wind speed when facing a 50 km-diameter storm eye.

4.4. A Case Study of Tropical Storm Karl

A case study of Tropical Storm Karl was carried out to statistically and visually compare the retrieval differences between the proposed EW Model 1, 2, 3, MMS1A, and S1EW.NR and especially investigate the extent of thermal noise’s impact on retrieval. The Sentinel-1 VV- and VH-polarized imageries, as shown in Figure 12a,b, were acquired at 22:22 UTC 29 on September 2016. Thermal noise had no impact on the VV-polarized image. However, for the VH-polarized image, the impact was obvious in low wind regions and the leftmost sub-swath 1. The collocated H*Wind data are shown in Figure 12c. Its time was 24:00 UTC on 29 September 2016. The retrievals of EW Model 1, 2, 3, MMS1A, and S1EW.NR models are shown in Figure 12d–i.
In the area framed in red in Figure 12d, wind variations were averaged in the azimuth direction and then compared in Figure 12i. A smaller rise and fall in the winds retrieved by EW Model 2 and 3 demonstrate that the impact of thermal noise was partially weakened. Even though the EW Model 1 was developed with the same MLR method as Models 2 and 3, its retrievals had a larger variation. For all the points whose wind speeds were lower than 10 m/s, the RMSE values were 2.46, 2.06, 2.30, 3.52, and 2.52 m/s for EW Model 1, 2, 3, MMS1A, and S1EW.NR, respectively. These results show the benefit of the VV and VH combination.

5. Discussion

As shown in Section 4.1, on average, our SAR retrievals were lower than SMAP measurements if the wind exceeded about 30 m/s. When the wind increased, the difference became larger. To investigate this reason, we compared the relationships between SAR VH-polarized NRCS and the winds measured using SFMR and SMAP. As shown in Figure 4 and Table A1, among the parameters of NRCS, the incident angle, and wind direction, VH NRCS had the strongest impact on wind speed retrieval. Thus, we used it as a medium to compare two different wind measurements and showed the variation characteristics of retrieval in different wind regimes.
The results are shown in Figure 13. Black solid and dashed curves represent the mean variations of scatters with different colors. According to the two curves, the difference in SAR retrievals and SMAP winds had two main reasons. First, if the wind exceeded 30 m/s for the same VH NRCS, the SMAP measurement would be higher than SFMR on average. Second, the difference between SFMR and SMAP increased with increasing VH NRCS. As a result, when we used the model developed from SFMR data to retrieve wind, the results were close to SFMR and lower than SMAP data. This difference became larger with increasing wind.
In the proposed models, wind direction is one of the inputs in EW and IW Model 3. However, according to the MLR results, it may not be a key parameter in wind speed retrieval if VV- and VH-polarized imageries are combined. As shown in Table A1, most A 4 and A i , 4 are at least one or two magnitudes smaller than other coefficients. It means that wind speed and direction are generally not related or are weakly related to each other. This is true for randomly selected sea areas in nature. Wind direction just modulates the relationship between co-polarized backscatter and wind speed by influencing the local incident angle in a cosine way, which is hard to fit using MLR. As a result, although Model 3 has a wind direction coefficient, it is not much better than Model 2 (see Table 2), indicating that the influence process of wind direction on wind retrieval is not essential if VV- and VH-polarized imageries are combined, while Model 2 and 3 are an evident improvement with respect to Model 1, indicating that the dual-polarization combination is an effective method to improve retrieval accuracy.
Under TC conditions, extreme wind speeds are always accompanied by intense rainfall, which can generate an error for wind retrieval. Since the SFMR instrument can measure contemporaneous wind speed and rain rate along the flight track, its data are an excellent reference to evaluate the impact of rainfall on wind retrieval [39]. In this study, the wind retrieval bias, which was equal to the proposed model retrieval minus the SFMR wind speed, was compared with the SFMR rain rate and wind speed. These results were used to measure the impact of rainfall on wind retrieval. For three models of each mode, their bias distributions were similar. Thus, biases were averaged from the results of Models 1, 2, and 3. As shown in Figure 14, although high rain rates are generally accompanied by high winds, there is no evidence that the rainfall had a definite impact on wind retrieval. For all points, the bias was around zero on average. Due to this ambiguous pattern, reducing errors with rain rate or further retrieving rain rate is still difficult.

6. Conclusions

Inspired by the CWAVE wave retrieval models, this study proposes a new approach for SAR ocean surface wind speed retrieval in a grid spacing of 1 km with the MLR method. To develop and validate new retrieval models, Sentinel-1 dual-polarized (VV + VH) observations were acquired from 28 TC cases over different ocean basins. Collocated wind data were collected from multiple sources, such as SFMR, SMAP, dropsonde, ECMWF ERA5, and H*Wind.
According to the comparison between NRCS, incident angle, and wind speed, VV- and VH-polarized imagery can distinguish wind speeds well, respectively, within low wind regions and moderate-to-high wind regions. Based on these advantages, new models were designed with the combination of VV- and VH-polarized signals. In addition, the GMF form was changed from the traditional σ 0 output to the new U 10 output, making the retrieval straighter and faster. Three schemes were investigated for inputs, named after Model 1 (inputs: σ 0 VH , θ ), Model 2 (inputs: σ 0 VH , θ , σ 0 VV ) and Model 3 (inputs: σ 0 VH , θ , σ 0 VV , φ ). According to the MLR results, the model with more inputs could fit data better. However, the improvement from Model 2 to Model 3 was not as significant as from Model 1 to Model 2.
Validation shows better consistencies between the results of proposed models and SMAP and dropsonde winds compared with traditional VH-polarization GMFs in the piecewise function. The big difference between SAR retrievals and SMAP winds over 30 m/s was mainly caused by the observational difference in SFMR and SMAP. A retrieval case study on Hurricane Michael indicated that our models could provide the wind information of a fine structure, for example, TC eye, eyewall, and rainbands. A case study of Tropical Storm Karl suggested that the impact of thermal noise on wind retrieval could be partially reduced using Models 2 and 3. These results indicate that the use of MLR and the combination of VV- and VH-polarized imageries could improve wind retrieval accuracy. In addition, the input of wind direction is not essential if the combination has been utilized. Under TC conditions, although intense rainfall is often accompanied by high winds, our study showed no evidence that rainfall had a definite impact on wind retrieval. Even with rainfall, the retrieval bias was around zero on average. As a result, it is still difficult to reduce errors with rain rate or further retrieve rain rate.

Author Contributions

Initiation of the idea: Y.G.; Data processing and model proposing: Y.G.; Writing and editing: all authors contributed; Supervision: Y.W.; Funding acquisition: Y.G. All authors have read and agreed to the published version of the manuscript.

Funding

This study is supported by the Natural Science Foundation of Shandong (Grant ZR2022QF069).

Data Availability Statement

The data are available at https://pan.baidu.com/s/19A7SVd5fuQscLn2DiVruwg (accessed on 30 August 2023).

Acknowledgments

We thank the European Space Agency for making the Sentinel-1 data publicly available. The authors thank the NOAA/Atlantic Oceanographic and Meteorological Laboratory/Hurricane Research Division (https://www.aoml.noaa.gov/ (accessed on 20 February 2023)) for providing the SFMR and Dropsonde data and the Remote Sensing Systems (RSS, https://www.remss.com/ (accessed on 21 March 2022)) for providing SMAP data. We thank the European Centre for Medium-Range Weather Forecasts for ERA5 data and Risk Management Solutions for H*Wind data. We thank the National Hurricane Center and the Joint Typhoon Warning Centre for TC best tracks.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

In this section, the coefficients are listed in Table A1 for the proposed wind retrieval models, i.e., EW Model 1, 2, 3 and IW Model 1, 2, 3. The coefficients in their correction functions are listed in Table A2.
Table A1. Values of all coefficients in the proposed wind retrieval models.
Table A1. Values of all coefficients in the proposed wind retrieval models.
ModelCoefficient
EW Model 1 A 0 A 1 A 2 A 1 , 1 A 1 , 2
134.9485278.5359061.12939050.14220560.038811
A 2 , 2
0.003917
EW Model 2 A 0 A 1 A 2 A 3 A 1 , 1
143.81241311.0672082.355905−0.3078380.204342
A 1 , 2 A 1 , 3 A 2 , 2 A 2 , 3 A 3 , 3
0.036087−0.071111−0.023669−0.064649−0.035267
EW Model 3 A 0 A 1 A 2 A 3 A 4
147.34819811.3988982.377266−0.4406410.000234
A 1 , 1 A 1 , 2 A 1 , 3 A 1 , 4 A 2 , 2
0.2090360.035286−0.076520−0.000547−0.023973
A 2 , 3 A 2 , 4 A 3 , 3 A 3 , 4 A 4 , 4
−0.065019−0.000177−0.0339610.000105−0.000017
IW Model 1 A 0 A 1 A 2 A 1 , 1 A 1 , 2
185.59335712.4659331.3152790.141039−0.054268
A 2 , 2
−0.029085
IW Model 2 A 0 A 1 A 2 A 3 A 1 , 1
203.54922015.0886891.653653−0.7141530.249729
A 1 , 2 A 1 , 3 A 2 , 2 A 2 , 3 A 3 , 3
−0.015968−0.085755−0.027735−0.050190−0.034910
IW Model 3 A 0 A 1 A 2 A 3 A 4
217.78063616.3275312.159972−1.552834−0.163730
A 1 , 1 A 1 , 2 A 1 , 3 A 1 , 4 A 2 , 2
0.269266−0.016449−0.108816−0.003335−0.035309
A 2 , 3 A 2 , 4 A 3 , 3 A 3 , 4 A 4 , 4
−0.0411200.000859−0.0206040.0016880.000183
Table A2. Values of all coefficients in the correction functions of the proposed wind retrieval models.
Table A2. Values of all coefficients in the correction functions of the proposed wind retrieval models.
Correction CoefficientEWIW
Model 1Model 2Model 3Model 1Model 2Model 3
a0.730.740.740.700.720.74
b1.121.111.111.131.121.11

References

  1. Li, X. The first Sentinel-1 SAR image of a typhoon. Acta Oceanol. Sin. 2015, 34, 1–2. [Google Scholar] [CrossRef]
  2. Zhang, B.; Perrie, W. Cross-Polarized Synthetic Aperture Radar: A New Potential Measurement Technique for Hurricanes. Bull. Am. Meteor. Soc. 2012, 93, 531–541. [Google Scholar] [CrossRef]
  3. Horstmann, J.; Schiller, H.; Schulz-Stellenfleth, J.; Lehner, S. Global Wind Speed Retrieval From SAR. IEEE Trans. Geosci. Remote Sens. 2003, 41, 2277–2286. [Google Scholar] [CrossRef]
  4. Ye, X.; Lin, M.; Zheng, Q.; Yuan, X.; Liang, C.; Zhang, B.; Zhang, J. A Typhoon Wind-Field Retrieval Method for the Dual-Polarization SAR Imagery. IEEE Geosci. Remote. Sens. Lett. 2019, 16, 1511–1515. [Google Scholar] [CrossRef]
  5. Li, X.; Pichel, W.G.; He, M.; Wu, S.Y.; Zhao, C. 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]
  6. Shao, W.; Jiang, X.; Sun, Z.; Hu, Y.; Marino, A.; Zhang, Y. Evaluation of wave retrieval for Chinese Gaofen-3 synthetic aperture radar. Geo-Spat. Inf. Sci. 2022, 25, 229–243. [Google Scholar] [CrossRef]
  7. Shao, W.; Hu, Y.; Jiang, X.; Zhang, Y. Wave retrieval from quad-polarized Chinese Gaofen-3 SAR image using an improved tilt modulation transfer function. Geo-Spat. Inf. Sci. 2023, 1–19. [Google Scholar] [CrossRef]
  8. Jiang, T.; Shao, W.; Yuyi, H.; Gang, Z.; Wei, S. L-Band Analysis of the Effects of Oil Slicks on Sea Wave Characteristics. J. Ocean Univ. China 2023, 22, 9–20. [Google Scholar] [CrossRef]
  9. Kudryavtsev, V.; Kozlov, I.; Chapron, B.; Johannessen, J.A. Quad-polarization SAR features of ocean currents. J. Geophys. Res. Oceans 2014, 119, 6046–6065. [Google Scholar] [CrossRef]
  10. Liu, B.; He, Y.; Li, X. A New Concept of Full Ocean Current Vector Retrieval with Spaceborne SAR Based on Intrapulse Beam-Switching Technique. IEEE Trans. Geosci. Remote Sens. 2020, 58, 7682–7704. [Google Scholar] [CrossRef]
  11. Fan, K.; Fu, B.; Gu, Y.; Yu, X.; Liu, T.; Shi, A.; Xu, K.; Gan, X. Internal wave parameters retrieval from space-borne SAR image. Front. Earth Sci. 2015, 9, 700–708. [Google Scholar] [CrossRef]
  12. Magalhães, J.; Werner, A.; Santos-Ferreira, A.M.; da Silva, J. Surface Wave Breaking Caused by Internal Solitary Waves: Effects on Radar Backscattering Measured by SAR and Radar Altimeter. Oceanography 2021, 34, 166–176. [Google Scholar] [CrossRef]
  13. Zhang, B.; Perrie, W. Recent progress on high wind-speed retrieval from multi-polarization SAR imagery: A review. Int. J. Remote Sens. 2014, 35, 4031–4045. [Google Scholar] [CrossRef]
  14. Gao, Y.; Zhang, J.; Guan, C.; Sun, J. Analyzing Sea Surface Wind Distribution Characteristics of Tropical Cyclone Based on Sentinel-1 SAR Images. Remote Sens. 2021, 13, 4501. [Google Scholar] [CrossRef]
  15. Zhang, G.; Li, X.; Perrie, W.; Zhang, J. Tropical Cyclone Winds and Inflow Angle Asymmetry from SAR Imagery. Geophys. Res. Lett. 2021, 48, e2021GL095699. [Google Scholar] [CrossRef]
  16. Zhou, X.; Yang, X.F.; Li, Z.W.; Yu, Y.; Bi, H.B.; Ma, S.; Li, X.F. Estimation of tropical cyclone parameters and wind fields from SAR images. Sci. China Earth Sci. 2013, 56, 1977–1987. [Google Scholar] [CrossRef]
  17. Gao, Y.; Zhang, J.; Sun, J.; Guan, C. Application of SAR Data for Tropical Cyclone Intensity Parameters Retrieval and Symmetric Wind Field Model Development. Remote Sens. 2021, 13, 2902. [Google Scholar] [CrossRef]
  18. Lagasio, M.; Parodi, A.; Pulvirenti, L.; Meroni, A.; Boni, G.; Pierdicca, N.; Marzano, F.; Luini, L.; Venuti, G.; Realini, E.; et al. A Synergistic Use of a High-Resolution Numerical Weather Prediction Model and High-Resolution Earth Observation Products to Improve Precipitation Forecast. Remote Sens. 2019, 11, 2387. [Google Scholar] [CrossRef]
  19. Tiesi, A.; Pucillo, A.; Bonaldo, D.; Ricchi, A.; Carniel, S.; Miglietta, M. Initialization of WRF Model Simulations With Sentinel-1 Wind Speed for Severe Weather Events. Front. Mar. Sci. 2021, 8, 169. [Google Scholar] [CrossRef]
  20. Stoffelen, A.; Anderson, D. Scatterometer data interpretation: Estimation and validation of the transfer function CMOD4. J. Geophys. Res. 1997, 102, 5767–5780. [Google Scholar] [CrossRef]
  21. 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]
  22. Stoffelen, A.; Verspeek, J.; Vogelzang, J.; Verhoef, A. The CMOD7 Geophysical Model Function for ASCAT and ERS Wind Retrievals. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2017, 10, 2123–2134. [Google Scholar] [CrossRef]
  23. Mouche, A.; Chapron, B. Global C-Band Envisat, RADARSAT-2 and Sentinel-1 SAR measurements in copolarization and cross-polarization. J. Geophys. Res. Oceans 2015, 120, 7195–7207. [Google Scholar] [CrossRef]
  24. Lu, Y.; Zhang, B.; Perrie, W.; Li, X.; Wang, H. A C-Band Geophysical Model Function for Determining Coastal Wind Speed Using Synthetic Aperture Radar. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2018, 11, 2417–2428. [Google Scholar] [CrossRef]
  25. Zhang, B.; Perrie, W.; Zhang, J.; Uhlhorn, E.; He, Y. High-Resolution Hurricane Vector Winds from C-Band Dual-Polarization SAR Observations. J. Atmos. Ocean. Technol. 2014, 31, 272–286. [Google Scholar] [CrossRef]
  26. Zhang, G.; Li, X.; Perrie, W.; Hwang, P.A.; Zhang, B.; Yang, X. A Hurricane Wind Speed Retrieval Model for C-Band RADARSAT-2 Cross-Polarization ScanSAR Images. IEEE Trans. Geosci. Remote Sens. 2017, 55, 4766–4774. [Google Scholar] [CrossRef]
  27. Mouche, A.A.; Chapron, B.; Zhang, B.; Husson, R. Combined Co- and Cross-Polarized SAR Measurements under Extreme Wind Conditions. IEEE Trans. Geosci. Remote Sens. 2017, 55, 6746–6755. [Google Scholar] [CrossRef]
  28. Gao, Y.; Sun, J.; Zhang, J.; Guan, C. Extreme Wind Speeds Retrieval Using Sentinel-1 IW Mode SAR Data. Remote Sens. 2021, 13, 1867. [Google Scholar] [CrossRef]
  29. Gao, Y.; Guan, C.; Sun, J.; Xie, L. Tropical Cyclone Wind Speed Retrieval from Dual-polarization Sentinel-1 EW Mode Products. J. Atmos. Ocean. Technol. 2020, 37, 1713–1724. [Google Scholar] [CrossRef]
  30. Won Jeong, P.; Korosov, A.; Babiker, M.; Sandven, S.; Won, J.-S. Efficient Thermal Noise Removal for Sentinel-1 TOPSAR Cross-Polarization Channel. IEEE Trans. Geosci. Remote Sens. 2017, 56, 1555–1565. [Google Scholar] [CrossRef]
  31. Horstmann, J.; Falchetti, S.; Wackerman, C.; Maresca, S.; Caruso, M.J.; Graber, H.C. Tropical Cyclone Winds Retrieved from C-Band Cross-Polarized Synthetic Aperture Radar. IEEE Trans. Geosci. Remote. Sens. 2015, 53, 2887–2898. [Google Scholar] [CrossRef]
  32. Hwang, P.A.; Stoffelen, A.; Zadelhoff, G.V.; Perrie, W.; Zhang, B.; Li, H.; Shen, H. Cross-polarization geophysical model function for C-band radar backscattering from the ocean surface and wind speed retrieval. J. Geophys. Res. Oceans 2015, 120, 893–909. [Google Scholar] [CrossRef]
  33. Hersbach, H.; Stoffelen, A.; de Haan, S. An improved C-band scatterometer ocean geophysical model function: CMOD5. J. Geophys. Res. Oceans 2007, 112, C03006. [Google Scholar] [CrossRef]
  34. Li, X.; 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]
  35. Sheng, Y.; Shao, W.; Shuai, Z.; Sun, J.; Yuan, X.; Li, S.; Shi, J.; Juncheng, Z. Validation of significant wave height retrieval from co-polarization Chinese Gaofen-3 SAR imagery using an improved algorithm. Acta Oceanol. Sin. 2018, 37, 1–10. [Google Scholar] [CrossRef]
  36. Shao, W.; Hu, Y.; Yang, J.; Nunziata, F.; Sun, J.; Li, H.; Zuo, J. An Empirical Algorithm to Retrieve Significant Wave Height from Sentinel-1 Synthetic Aperture Radar Imagery Collected under Cyclonic Conditions. Remote Sens. 2018, 10, 1367. [Google Scholar] [CrossRef]
  37. Uhlhorn, E.W.; Black, P.G.; Franklin, J.L.; Goodberlet, M.; Carswell, J.; Goldstein, A.S. Hurricane Surface Wind Measurements from an Operational Stepped Frequency Microwave Radiometer. Mon. Weather Rev. 2007, 135, 3070–3085. [Google Scholar] [CrossRef]
  38. Meissner, T.; Ricciardulli, L.; Wentz, F.J. Capability of the SMAP Mission to Measure Ocean Surface Winds in Storms. Bull. Am. Meteor. Soc. 2017, 98, 1660–1677. [Google Scholar] [CrossRef]
  39. Mouche, A.; Chapron, B.; Knaff, J.; Zhao, Y.; Zhang, B.; Combot, C. Copolarized and Cross-Polarized SAR Measurements for High-Resolution Description of Major Hurricane Wind Structures: Application to Irma Category 5 Hurricane. J. Geophys. Res. Oceans 2019, 124, 3905–3922. [Google Scholar] [CrossRef]
  40. Shao, W.; Lai, Z.; Nunziata, F.; Buono, A.; Jiang, X.; Zuo, J. Wind Field Retrieval with Rain Correction from Dual-Polarized Sentinel-1 SAR Imagery Collected during Tropical Cyclones. Remote Sens. 2022, 14, 5006. [Google Scholar] [CrossRef]
  41. Reul, N.; Chapron, B.; Zabolotskikh, E.; Donlon, C.; Quilfen, Y.; Guimbard, S.; Piolle, J.F. A revised L-band radio-brightness sensitivity to extreme winds under Tropical Cyclones: The five year SMOS-storm database. Remote Sens. Environ. 2016, 180, 274–291. [Google Scholar] [CrossRef]
  42. Yueh, S.H.; Fore, A.G.; Tang, W.; Hayashi, A.; Stiles, B.; Reul, N.; Weng, Y.; Zhang, F. SMAP L-Band Passive Microwave Observations of Ocean Surface Wind during Severe Storms. IEEE Trans. Geosci. Remote Sens. 2016, 54, 7339–7350. [Google Scholar] [CrossRef]
  43. Fore, A.G.; Yueh, S.H.; Stiles, B.W.; Tang, W.; Hayashi, A.K. SMAP Radiometer-Only Tropical Cyclone Intensity and Size Validation. IEEE Geosci. Remote Sens. Lett. 2018, 15, 1480–1484. [Google Scholar] [CrossRef]
  44. Hock, T.F.; Franklin, J.L. The NCAR GPS dropwindsonde. Bull. Am. Meteor. Soc. 1999, 80, 407–420. [Google Scholar] [CrossRef]
  45. Hersbach, H.; Bell, B.; Berrisford, P.; Biavati, G.; Horányi, A.; Muñoz Sabater, J.; Nicolas, J.; Peubey, C.; Radu, R.; Rozum, I.; et al. ERA5 Hourly Data on Single Levels from 1940 to Present. Available online: https://cds.climate.copernicus.eu/cdsapp#!/dataset/reanalysis-era5-single-levels?tab=overview (accessed on 27 February 2023).
  46. Powell, M.; Houston, S.; Amat, L.; Morisseau-Leroy, N. The HRD real-time hurricane wind analysis system. J. Wind Eng. Ind. Aerodyn. 1998, 77–78, 53–64. [Google Scholar] [CrossRef]
  47. DiNapoli, S.; Bourassa, M.; Powell, M. Uncertainty and Intercalibration Analysis of H*Wind. J. Atmos. Ocean. Technol. 2012, 29, 822–833. [Google Scholar] [CrossRef]
Figure 1. Intensities and numbers of the TC cases observed by Sentinel-1 with two main wind matches. The intensities were represented by the Saffir–Simpson hurricane wind scale and estimated from the best tracks provided by NHC and JTWC.
Figure 1. Intensities and numbers of the TC cases observed by Sentinel-1 with two main wind matches. The intensities were represented by the Saffir–Simpson hurricane wind scale and estimated from the best tracks provided by NHC and JTWC.
Remotesensing 15 04267 g001
Figure 2. A collocation example of Hurricane Michael. (a) Sentinel-1 EW mode VV-polarized image acquired on 23:44 UTC 9 October 2018; (b) Contemporaneous VH-polarized image; and (c) Matched SMAP final winds measured on 23:36 UTC 9 October 2018. In (a,b), white and black lines crossing the storm eye are the tracks of flights from AFRC and NOAA, respectively. The SFMR observations along the tracks are within a 0.5 h time difference from Sentinel-1.
Figure 2. A collocation example of Hurricane Michael. (a) Sentinel-1 EW mode VV-polarized image acquired on 23:44 UTC 9 October 2018; (b) Contemporaneous VH-polarized image; and (c) Matched SMAP final winds measured on 23:36 UTC 9 October 2018. In (a,b), white and black lines crossing the storm eye are the tracks of flights from AFRC and NOAA, respectively. The SFMR observations along the tracks are within a 0.5 h time difference from Sentinel-1.
Remotesensing 15 04267 g002
Figure 3. A collocation example of Hurricane Irma. (a) Sentinel-1 IW mode VV-polarized image acquired on 9:07 UTC 29 October 2017; (b) Contemporaneous VH-polarized image. White lines show the AFRC flight track. The SFMR observations along the track are within a 0.5 h time difference from Sentinel-1.
Figure 3. A collocation example of Hurricane Irma. (a) Sentinel-1 IW mode VV-polarized image acquired on 9:07 UTC 29 October 2017; (b) Contemporaneous VH-polarized image. White lines show the AFRC flight track. The SFMR observations along the track are within a 0.5 h time difference from Sentinel-1.
Remotesensing 15 04267 g003
Figure 4. Comparison between NRCS, incident angle and wind speed for (a) EW mode VV-polarization; (b) EW mode VH-polarization; (c) IW mode VV-polarization; (d) IW mode VH-polarization. The winds of the unfilled points were measured using SFMR. The winds of the filled points are from the ERA5 reanalysis. Black curves stand for NESZ.
Figure 4. Comparison between NRCS, incident angle and wind speed for (a) EW mode VV-polarization; (b) EW mode VH-polarization; (c) IW mode VV-polarization; (d) IW mode VH-polarization. The winds of the unfilled points were measured using SFMR. The winds of the filled points are from the ERA5 reanalysis. Black curves stand for NESZ.
Remotesensing 15 04267 g004aRemotesensing 15 04267 g004b
Figure 5. Comparison of the model-simulated and fitted wind speeds for different models: (a) EW Model 1; (b) EW Model 2; (c) EW Model 3. Red curve stands for the additional linear regression.
Figure 5. Comparison of the model-simulated and fitted wind speeds for different models: (a) EW Model 1; (b) EW Model 2; (c) EW Model 3. Red curve stands for the additional linear regression.
Remotesensing 15 04267 g005
Figure 6. Comparison of the model-simulated and fitted wind speeds for different models: (a) IW Model 1; (b) IW Model 2; (c) IW Model 3. Red curve stands for the additional linear regression.
Figure 6. Comparison of the model-simulated and fitted wind speeds for different models: (a) IW Model 1; (b) IW Model 2; (c) IW Model 3. Red curve stands for the additional linear regression.
Remotesensing 15 04267 g006
Figure 7. Sample number distributions of the SMAP wind speeds used for (a) EW mode; (b) IW mode; the ERA5 wind directions used for (c) EW mode; (d) IW mode.
Figure 7. Sample number distributions of the SMAP wind speeds used for (a) EW mode; (b) IW mode; the ERA5 wind directions used for (c) EW mode; (d) IW mode.
Remotesensing 15 04267 g007aRemotesensing 15 04267 g007b
Figure 8. Comparison between SMAP wind speed measurements and SAR-retrieved wind speeds based on (a) EW Model 1; (b) EW Model 2; (c) EW Model 3; (d) MMS1A; (e) S1EW.NR. Red curve stands for the variation trend.
Figure 8. Comparison between SMAP wind speed measurements and SAR-retrieved wind speeds based on (a) EW Model 1; (b) EW Model 2; (c) EW Model 3; (d) MMS1A; (e) S1EW.NR. Red curve stands for the variation trend.
Remotesensing 15 04267 g008aRemotesensing 15 04267 g008b
Figure 9. Comparison between SMAP wind speed measurements and the SAR-retrieved wind speeds based on (a) IW Model 1; (b) IW Model 2; (c) IW Model 3; (d) MMS1A; (e) S1IW.NR. Red curve stands for the variation trend.
Figure 9. Comparison between SMAP wind speed measurements and the SAR-retrieved wind speeds based on (a) IW Model 1; (b) IW Model 2; (c) IW Model 3; (d) MMS1A; (e) S1IW.NR. Red curve stands for the variation trend.
Remotesensing 15 04267 g009
Figure 10. Comparison of Dropsonde measurements with wind speeds retrieved from Sentinel-1 (a) EW mode images; (b) IW mode images, respectively.
Figure 10. Comparison of Dropsonde measurements with wind speeds retrieved from Sentinel-1 (a) EW mode images; (b) IW mode images, respectively.
Remotesensing 15 04267 g010
Figure 11. Comparison of the SAR, SFMR and SMAP wind measurements in a case study on Hurricane Michael. (a) Surface wind speeds retrieved via EW Model 3 from Sentinel-1 image. (b) Collocated SMAP wind map. (c) Wind comparison along the SFMR observation track (i.e., the white line in (a)). The grid spacing of SAR and SFMR winds was about 1 km. The grid spacing of SMAP winds was 0.25°.
Figure 11. Comparison of the SAR, SFMR and SMAP wind measurements in a case study on Hurricane Michael. (a) Surface wind speeds retrieved via EW Model 3 from Sentinel-1 image. (b) Collocated SMAP wind map. (c) Wind comparison along the SFMR observation track (i.e., the white line in (a)). The grid spacing of SAR and SFMR winds was about 1 km. The grid spacing of SMAP winds was 0.25°.
Remotesensing 15 04267 g011
Figure 12. A case study of Tropical Storm Karl for model comparison. Sentinel-1 imageries in (a) VV-polarization; (b) VH-polarization. (c) The collocated H*Wind data. Wind retrievals of (d) EW Model 1; (e) EW Model 2; (f) EW Model 3; (g) MMS1A; (h) S1EW.NR. (i) Averaged retrieval variations in azimuth direction according to the area framed in red in (d).
Figure 12. A case study of Tropical Storm Karl for model comparison. Sentinel-1 imageries in (a) VV-polarization; (b) VH-polarization. (c) The collocated H*Wind data. Wind retrievals of (d) EW Model 1; (e) EW Model 2; (f) EW Model 3; (g) MMS1A; (h) S1EW.NR. (i) Averaged retrieval variations in azimuth direction according to the area framed in red in (d).
Remotesensing 15 04267 g012aRemotesensing 15 04267 g012b
Figure 13. Relationships between Sentinel-1 VH NRCS and wind speeds measured using SFMR and SMAP according to our samples. The black dashed curve stands for the mean variation in SAR and SFMR. The black solid curve stands for the mean variation in SAR and SMAP.
Figure 13. Relationships between Sentinel-1 VH NRCS and wind speeds measured using SFMR and SMAP according to our samples. The black dashed curve stands for the mean variation in SAR and SFMR. The black solid curve stands for the mean variation in SAR and SMAP.
Remotesensing 15 04267 g013
Figure 14. Comparisons between wind retrieval bias, SFMR rain rate and wind speed for the proposed models of (a) EW mode; (b) IW mode. For each mode, results are averaged from Model 1, 2, and 3.
Figure 14. Comparisons between wind retrieval bias, SFMR rain rate and wind speed for the proposed models of (a) EW mode; (b) IW mode. For each mode, results are averaged from Model 1, 2, and 3.
Remotesensing 15 04267 g014
Table 1. Some imagery parameters of the Sentinel-1 EW and IW modes.
Table 1. Some imagery parameters of the Sentinel-1 EW and IW modes.
Sensor ModeGround Swath (km)Spatial Resolution 1 (m)Incident Angle (°)Sub-Swaths Number
EW41030 × 4020–475
IW2505 × 2031–463
1 The spatial resolution is range × azimuth.
Table 2. Evaluation of the new models based on fitted data.
Table 2. Evaluation of the new models based on fitted data.
ModelBias (m/s)CorRMSE (m/s)
<10 m/s≥10 m/sAll<10 m/s≥10 m/sAll<10 m/s≥10 m/sAll
EW Model 10.75–0.320.130.040.920.932.664.093.57
EW Model 20.42–0.170.070.400.920.942.413.793.30
EW Model 30.38–0.130.080.400.920.942.423.773.28
IW Model 11.24–0.860.00450.270.900.912.775.184.36
IW Model 20.51–0.340.00970.670.900.932.594.303.69
IW Model 30.24–0.160.00490.640.910.942.344.053.45
Table 3. Validation results of the proposed EW Model 1, 2, 3, and traditional MMS1A and S1EW.NR models. Retrievals were compared with SMAP wind speeds.
Table 3. Validation results of the proposed EW Model 1, 2, 3, and traditional MMS1A and S1EW.NR models. Retrievals were compared with SMAP wind speeds.
ModelBias (m/s)CorRMSE (m/s)
<10 m/s≥10 m/sAll<10 m/s≥10 m/sAll<10 m/s≥10 m/sAll
EW Model1–0.70–1.47–1.360.260.940.951.763.603.37
EW Model20.10–1.53–1.270.390.950.961.743.303.10
EW Model3–0.16–1.64–1.410.410.950.961.793.343.15
MMS1A0.71–2.37–1.880.440.950.962.014.304.02
S1EW.NR–0.77–3.89–3.400.210.920.923.215.925.58
Table 4. Validation results of the proposed IW Model 1, 2, 3, and traditional MMS1A and S1IW.NR models. Retrievals were compared with SMAP wind speeds.
Table 4. Validation results of the proposed IW Model 1, 2, 3, and traditional MMS1A and S1IW.NR models. Retrievals were compared with SMAP wind speeds.
ModelBias (m/s)CorRMSE (m/s)
<10 m/s≥10 m/sAll<10 m/s≥10 m/sAll<10 m/s≥10 m/sAll
IW Model11.18–3.14–2.260.270.930.942.714.554.24
IW Model21.36–1.79–1.140.580.940.962.433.553.35
IW Model3–1.07–1.73–1.600.470.890.932.564.394.08
MMS1A6.45–4.04–1.88–0.010.930.936.966.036.22
S1IW.NR0.12–4.23–3.330.440.890.922.445.905.38
Table 5. Validation results of the proposed EW Model 1, 2, 3, and traditional MMS1A and S1EW.NR models. Retrievals were compared with Dropsonde wind speeds.
Table 5. Validation results of the proposed EW Model 1, 2, 3, and traditional MMS1A and S1EW.NR models. Retrievals were compared with Dropsonde wind speeds.
ModelBias (m/s)CorRMSE (m/s)
EW Model 10.050.906.69
EW Model 20.370.906.76
EW Model 30.300.906.71
MMS1A–1.520.907.26
S1EW.NR–3.600.889.37
Table 6. Validation results of the proposed IW Model 1, 2, 3, and traditional MMS1A and S1IW.NR models. Retrievals were compared with Dropsonde wind speeds.
Table 6. Validation results of the proposed IW Model 1, 2, 3, and traditional MMS1A and S1IW.NR models. Retrievals were compared with Dropsonde wind speeds.
ModelBias (m/s)CorRMSE (m/s)
IW Model 1–0.930.956.22
IW Model 2–1.620.946.75
IW Model 3–0.770.917.88
MMS1A–5.420.9610.60
S1IW.NR–4.850.949.03
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Gao, Y.; Wang, Y.; Wang, W. A New Approach for Ocean Surface Wind Speed Retrieval Using Sentinel-1 Dual-Polarized Imagery. Remote Sens. 2023, 15, 4267. https://doi.org/10.3390/rs15174267

AMA Style

Gao Y, Wang Y, Wang W. A New Approach for Ocean Surface Wind Speed Retrieval Using Sentinel-1 Dual-Polarized Imagery. Remote Sensing. 2023; 15(17):4267. https://doi.org/10.3390/rs15174267

Chicago/Turabian Style

Gao, Yuan, Yunhua Wang, and Weili Wang. 2023. "A New Approach for Ocean Surface Wind Speed Retrieval Using Sentinel-1 Dual-Polarized Imagery" Remote Sensing 15, no. 17: 4267. https://doi.org/10.3390/rs15174267

APA Style

Gao, Y., Wang, Y., & Wang, W. (2023). A New Approach for Ocean Surface Wind Speed Retrieval Using Sentinel-1 Dual-Polarized Imagery. Remote Sensing, 15(17), 4267. https://doi.org/10.3390/rs15174267

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