iBet uBet web content aggregator. Adding the entire web to your favor.
iBet uBet web content aggregator. Adding the entire web to your favor.



Link to original content: https://doi.org/10.3390/rs14195006
Wind Field Retrieval with Rain Correction from Dual-Polarized Sentinel-1 SAR Imagery Collected during Tropical Cyclones
Next Article in Journal
Quantifying Plant Species α-Diversity Using Normalized Difference Vegetation Index and Climate Data in Alpine Grasslands
Next Article in Special Issue
Evaluation of Multi-Incidence Angle Polarimetric Gaofen-3 SAR Wave Mode Data for Significant Wave Height Retrieval
Previous Article in Journal
Estimating the Legacy Effect of Post-Cutting Shelterbelt on Crop Yield Using Google Earth and Sentinel-2 Data
Previous Article in Special Issue
Improving the Modeling of Sea Surface Currents in the Persian Gulf and the Oman Sea Using Data Assimilation of Satellite Altimetry and Hydrographic Observations
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Wind Field Retrieval with Rain Correction from Dual-Polarized Sentinel-1 SAR Imagery Collected during Tropical Cyclones

1
College of Marine Sciences, Shanghai Ocean University, Shanghai 201306, China
2
Key Laboratory of Marine Ecological Monitoring and Restoration Technologies, Ministry of Natural Resources, Shanghai 201206, China
3
Dipartimento di Ingegneria, Università degli Studi di Napoli Parthenope, 80143 Napoli, Italy
4
Southern Marine Science and Engineering Guangdong Laboratory (Guangzhou), Guangzhou 511458, China
5
National Satellite Ocean Application Service, Ministry of Natural Resources of the People’s Republic of China, Beijing 100081, China
*
Author to whom correspondence should be addressed.
Remote Sens. 2022, 14(19), 5006; https://doi.org/10.3390/rs14195006
Submission received: 1 September 2022 / Revised: 26 September 2022 / Accepted: 30 September 2022 / Published: 8 October 2022

Abstract

:
The purpose of this study is to include rain effects in wind field retrieval from C-band synthetic aperture radar (SAR) imagery collected under tropical cyclone conditions. An effective and operationally attractive approach to detect rain cells in SAR imagery is proposed and verified using four Sentinel-1 (S-1) SAR images collected in dual-polarized (vertical-vertical (VV) and vertical-horizontal (VH)) interferometric-wide swath imaging mode during the Satellite Hurricane Observation Campaign. SAR images were collocated with ancillary observations that include sea surface wind and rain rate from the Stepped-Frequency Microwave Radiometer (SFMR) on board of the National Oceanic and Atmospheric Administration aircraft. The winds are inverted from VV- and VH-polarized S-1 image using the CMOD5.N and S1IW.NR geophysical model functions (GMFs), respectively. Location and radius of cyclone’s eye, together with the TC central pressure, are calculated from the VV-polarized SAR-derived wind and a parametric model. A cost function is proposed that consists of the difference between the measured VV-polarized SAR normalized radar cross section (NRCS) and the NRCS predicted using CMOD5.N forced with the wind speed retrieved by the VH-polarized SAR images using S1IW.NR GMF and the wind direction retrieved from the patterns visible in the SAR image. This cost function is related to the SFMR rain rate. Experimental results show that the difference between measured and predicted NRCS values range from 0.5 dB to 5 dB within a distance of 100 km from the cyclone’s eye, while the difference increases spanning from 3 dB to 6 dB for distances larger than 100 km. Following this rationale, first the rain bands are extracted from SAR imagery and, then, the composite wind fields are reconstructed by replacing: (1) dual-polarized SAR-derived winds over the rain-free regions; (2) winds simulated using the radial-vortex model over the rain-affected regions. The validation of the composite wind speed against SFMR winds yields a <2 m s−1 and >0.7 correlation (COR) at all flow directions up to retrieval speeds of 70 m s−1. This result outperforms the winds estimated using the VH-polarized S1IW.NR GMF, which call for high error accuracy, such as about 4 m s−1 with a 0.45 COR ranged from 330° to 360°.

1. Introduction

Tropical cyclones (TCs) always result in strong winds, extreme waves and heavy rainfall, which are significant disasters on near-shore waters [1]. The met-ocean conditions associated to TCs make real-time observations from ships and moored buoys unfeasible. Instead, remote-sensed technology is useful for rainfall monitoring [2]. The Tropical Rainfall Measuring Mission (TRMM)—a well-known international cooperation project designed for global precipitation observations—operating at microwaves has released three-dimensional data since November 1997 [3,4]. However, remote-sensed TRMM rainfall data are not suitable for the meso- and small-scale studies concerning the air–sea interactions due to the relative coarse spatial resolution of approximately 0.25° grid (about 25 km).
At present, synthetic aperture radar (SAR) images is an active and fine spatial resolution (down to 1 m) sensor with an all-day and almost all-weather capability of sea surface monitoring in polar region [5] and TCs, e.g., the morphology of cyclones [6], strong winds [7,8] and extreme waves [9]. Rain footprints over sea surface were also qualitative observed in SAR imagery [10]. Traditionally, in SAR images, three main aspects rule the rain-induced effects on sea surface roughness [11,12]: disturbances on wind-generated short waves winds [13,14]; splash products, i.e., volume scattering and attenuation due to raindrops at the air-sea layer [15]; large-scale downdraft roughness modulated by rain-induced wind flow [16]. Accordingly, rain cells affect the homogeneity of the SAR image distribution and the normalized radar cross-sections (NRCS) measured by the SAR. Even though the rain-induced characteristics on co-polarized SAR backscattering at C- and X-band was investigated [17,18], the rain attenuation of the SAR wind signals at co- and cross-polarized channel is complicated. Moreover, the correction of rain cells is needed for SAR-derived wind retrieval under extreme weather using a practical scheme.
The wind retrieval methodology from SAR image is well developed in last decades [8,19]. SAR-based wind retrieval relies on the relationship between sea surface wind and co- and cross-polarized SAR measurements [20], through the development of geophysical model functions (GMFs) [21,22,23]. The validation of wind retrieval against the measurements from scatterometer [24] and moored buoys yields a 2 m s−1 accuracy of wind speed using co-polarized C-band [25,26] and X-band SAR [27] GMF. However, it is revealed that co-polarized GMFs suffer the saturation at high winds, i.e., wind speed larger than 25 m s−1 [28]. Utilizing the RADARSAT-2 (R-2) fine quad mode data collocated with operational weather buoys [29], it is found that the cross-polarized SAR backscattering signal is linearly related with wind speed up to 55 m s−1 [30,31] and is independent of incidence angle and wind direction. Although cross-polarized GMFs work well at very high winds [7,32], the accuracy of wind field retrieval from cross-polarized NRCS SAR measurements at low-to-moderate winds is strictly related to the signal-to-noise ratio [33], which is particularly important for low-backscatter cross-pol measurements at small incidence angles.
In this study, four dual-polarized Sentinel-1 (S-1) C-band SAR images are considered, which were taken by the Satellite Hurricane Observation Campaign (SHOC), providing high-quality data for wind [18] and wave research [34] in TC. These images were collocated with the wind field and rain rate observations provided by the Stepped-Frequency Microwave Radiometer (SFMR) on-board of the National Oceanic and Atmospheric Administration (NOAA) aircraft. The TC parameters, including location of cyclone’s eye, central pressure, maximum wind speed and maximum wind radius, are calculated based on the VV-polarized SAR wind retrieval algorithm and a parametric model developed by Holland (1980). A practical approach for the identification of rain cells is developed when analyzing the cost function of the SAR backscattering related to the SFMR rain rate. The information on TCs is used to produce a wind field using a radial-vortex model [35] and the regions identified as raining is replaced by the simulated wind field. Finally, a composite wind field is derived from S-1 SAR-derived wind field at VV- and VH- polarization channel after taking the corrections from the radial-vortex modeling wind field.
The remainder of this paper is organized as follows: the dataset is briefly described in Section 2, including the four dual-polarized S-1 IW SAR images collocate with the SFMR measurements in TC and the SAR-derived wind fields from VV-polarized and VH-polarized images. Besides, the methodology of the rain cell identification from SAR data, estimation of TC parameters and the parametric wind profile model is introduced in Section 2. Section 3 presents the retrieved wind speeds after rain corrections. The conclusions are summarized in Section 4.

2. Materials and Methods

2.1. Materials

In this study, four S-1 SAR images acquired in dual-polarized (VV+VH) interferometric-wide (IW) swath imaging mode during TCs are considered with the clear cyclone eye, see Table 1. They are collocated with sea surface wind and rain rate observations obtained by SFMR, which benefits the comprehensive TC studies using SAR data [36]. The SFMR provide near real-time observations of the wind speed and rain rate within and around tropical cyclone [37]. The SAR wind retrieval is essential for the identification of rain cell. The well-developed GMF at C-band, denoted as CMOD4 [38], has been used for SAR wind retrieval since 1997. The performance of GMF is improved by the updated CMOD family [39] using the following formulation:
σ 0 = B 0 ( U 10 , θ ) ( 1 + B 1 ( U 10 , θ ) cos ϕ + B 2 ( U 10 , θ ) cos 2 ϕ )
where σ 0 is the linear NRCS measured by SAR, U10 is the wind speed at 10-m height, d φ represents the angle between the wind direction and the radar look direction and matrix B (B0, B1, and B2) are functions of the incidence angle θ and the wind speed at 10 m above the sea surface. Because SAR operates using a single look radar beam, a wind vector is impossible to be simultaneously retrieved due to two unknown variables (wind speed and direction) in the GMF. In practice, the wind direction is extracted from the two-dimensional SAR intensity spectrum at wavelengths between 800 m and 3000 m [40] or directly measured following the wind streaks [41]. In this study, the wind directions are estimated from wind-induced low-frequency SAR image features and the 180° ambiguity is removed using external European Centre for Medium-Range Weather Forecasts (ECMWF)-interim (ERA-5) data with a spatial resolution of approximately 0.25° grid. The complete field of wind directions that will be considered in the SAR wind speed retrieval process is obtained using the interpolation method proposed in [42] and we do not repeat it here.
As a showcase, the quick-looks of the VV-polarized S-1 SAR images belonging to the dataset are shown in Figure 1, where subfigures refer to TCs: (a) Irma at 10:30 UTC on 7 September 2017, (b) Dorian at 22:46 UTC on 30 August 2019, (c) Lsaias at 23:19 UTC on 2 August 2020 and (d) Delta at 00:07 UTC on 8 October 2020. The red lines represent the SFMR tracks. Recently, aimed to wind retrieval from VH-polarized S-1 SAR images after thermal noise removal, an algorithm termed as S-1 IW Mode Wind Speed Retrieval Model after Noise Removal (S1IW.NR) was developed [43,44] and the specific formulation is expressed as follows, considering the different performance at various incidence angle θ, e.g., IW1: [31.0°–35.9°], IW2: [35.9°–41.3°] and IW3: [41.3°–46.0°]:
σ 0 V H [ d B ] = { 0.22 U 10 29.68 I W 1 4.67 U 10 0.39 41.02 I W 2 56.67 U 10 0.26 I W 3
where σ 0 V H is the VH-polarized NRCS united in dB and U10 is the wind speed at 10-m height. The cross-polarized GMF is conveniently applied for SAR strong wind retrieval without be prior information on wind direction. Although maximum wind speeds that can be accurately retrieved from VH-polarized SAR images could be up to 50 m s−1 [7,8], the rain effect is not accounted for in the inversion algorithm and the retrieval accuracy is reduced compared to that from VV-polarized SAR images at low-to-moderate winds due to the low signal-to-noise ratio. The inverted four wind maps obtained from VV-polarized S-1 SAR images using the CMOD5.N GMF and VH-polarized S-1 SAR images using the S1IW.NR are shown in Figure 2 and Figure 3, respectively. It can be clearly observed that the VH-polarized SAR-derived wind speeds are significantly greater than retrieval results from VV-polarized SAR images. The statistical analysis of Taylor diagram between VH-polarized SAR-derived wind speeds and the available SFMR measurements in term of flow directions of TC eyes is presented in Figure 4. The flow direction discrete at interval of a 30°. It is found that maximum root mean square error (RMSE) of wind speed is about 4 m s−1 with a 0.45 correlation (COR) ranged from 330° to 360°; however, the RMSE of wind speed is less than 3 m s−1 at other flow directions. Although this result suggests that, in this study, the VH-polarized SAR-derived wind speeds could be applied up to 70 m s−1 due to the retrieval floor, the accuracy is still further improved at special flow directions.

2.2. Method

In this section, the methodology proposed to identify rain cells and to develop TC wind maps from dual-polarized S-1 SAR imagery is described. Since the SAR-derived wind fields with rainfall parameter are not reliable due to the attenuations in dual-polarized SAR imagery, we used the Holland model to construct the TC wind profile.

2.2.1. Rain Cells Identification

The rainfall affects significantly the microwave signal interacting with rain cells. Moreover, the attenuation due to the rainfall is polarization selective, i.e., the co-polarized backscatter exhibits a sensitivity to the rainfall larger than the cross-polarized one. A recent study proposed by [17] showed that the radar backscattering increases with rain rate at low-to-moderate winds. The basic principle is based on the complementary sensitivity of the SAR signal induced by rain cell. Following this rationale, the scheme of detecting rain in SAR data is designed by a complicated cost function relying on four parameters derived from the local gradient of SAR intensity image [45]. Besides, a simple approach is proposed in [35], defined as the difference between the measured VV-polarized NRCS σ V V 0 , and the predicted one σ S 0 ; however, the influence of radar incidence angle is excluded in this method. This latter one is predicted using the CMOD5.N that is forced with a wind speed retrieved from the VH-polarized SAR scene using the S1IW.NR and a wind direction derived from the joint analysis of the pattern visible in the SAR image and using ERA-5 winds. Once the cost function is defined, a threshold value T is set as follows:
T = | σ V V 0 σ S 0 |   [ dB ]
To summarize, the following three steps are needed: (1) wind speed retrieval from VH-polarized SAR image using GMF S1IW.NR; (2) prediction of the GMF CMOD5.N-simulated NRCS using SAR-derived wind direction and VH-polarized wind speed retrieved as in (1); (3) applying Equation (3) in term of radar incidence angle from SAR image to calculate the T. It is noted that the simple identification method has worse performance than that using local gradient method based on backscattering difference between co-polarized and cross-polarized S-1 image [46], because the information on weak rainfall at regular sea state is missing using simple identification method. However, the efficiency of the approach herein has been confirmed using R-2 images in hurricanes, where rainfall is strong at such a condition.

2.2.2. Estimation of TC Parameters

Although the rain flag is a valuable assessment for rain monitoring, quality-flagged wind fields are usually required for TC analyses, i.e., assimilating SAR-derived winds into comprehensive numerical prediction models and the distribution of extreme waves produced by high winds. Within this framework, rainfall corrections on the rain-flagged winds are needed to obtain quick estimations of the complete TC wind field. The well-known parametric Holland wind model [47] is an analytical model for the radial profile of wind speed under TC conditions:
V g   = [ A B ( p n p c ) e A / r B ρ r 2 + r 2 f 2 4 ] 1 / 2   r f 2 + V T C  
where Vg is the wind speed at radius r, pc and pn are the central and ambient pressure, respectively, f is the Coriolis force parameter, ρ is the air density constant, VTC is the movement velocity of TC and the parameters A and B represent the location relative to the origin and the shape related to the maximum wind radius rmax and maximum wind radius vmax:
A = r max B
B = ρ e v m a x 2 p n   p c
In this study, the method proposed in [48] to estimate TC parameters is employed. It is a three-step procedure that consists of: (1) an edge detection technique based on the “Daubechies D4” wavelet applied on the VV-polarized NRCS SAR image to determine the TC’s eye, because the eye of the TC is detected as a low-frequency feature at this scale; (2) select two cuts in the range and azimuth directions through TC eyes and the average distance of the two peak values in each cut, that is maximum wind radius vmax; (3) the wind speed retrieval using the CMOD5.N GMF and the selection of regions where wind speed is lower than 20 m s−1 which are less affected by rain; (4) the minimization of the cost function:
J ( p c ,   v m a x ) = i j | σ V V , i j 0 σ Holland , ij 0 ( p c ,   v m a x ) |
where σ Holland , ij 0 ( p c ,   v m a x ) is the simulated NRCS by CMOD5.N GMF using Holland winds and σ VV , ij 0 is the VV-polarized NRCS measured by the SAR. Results relevant to the comparison between SAR-derived TC intensity, e.g., central pressure pc, maximum wind speed vmax and maximum wind radius rmax, and the SFMR observations are presented in Figure 5. It is found that the RMSEs of central pressure and maximum wind speed are 15.6 hPa and 10.8 m s−1, respectively, comparing the retrieval results with the best-track information of NOAA.

2.2.3. TC Wind Radial Profile Model

A TC is mesoscale cyclonic system with a low-pressure center and strong wind speeds around an eyewall, which maintains a unique wind profile along the radial direction. According to the rational physics, the TC wind profile has strict relation with the strength of the TC [47]. This kind of radial profile model has been widely used to reconstruct cyclostrophic wind field as knowing TC best tracks. The revised radial profile model was recently developed [49] consisting of an outer wind and surface pressure, sea surface temperature, radius of maximum winds, and central pressure. Less additional information is required in the SAR wind retrieval process; thus, the following equation is employed for reconstructing radial wind filed:
V r   = { v max ( r r max )   r < r max v max ( r r max ) 0.5 r r max  
where Vr is the radial wind speed at the distance r from TC eye. As mentioned in [35], the parameters in Equation (8) were calculated by fitting the function to the VH-polarized SAR wind data along each radial direction by a least square method depending on the option of radial direction. In this study, the maximum wind speed vmax and maximum wind radius rmax are estimated using the method in the above section.
Collectively, the composite wind field after rain correction is reconstructed by VV-polarized SAR winds retrieved using CMOD5.N, VH-polarized SAR winds retrieved using S1IW.NR and the vortex-modeling TC wind field. The flowchart of the methodology is depicted in Figure 6.

3. Results

In this section, the identification results of rain cell for four S-1 SAR image in TCs are first presented. Then, a composited wind fields from dual-polarized images with rain corrections are validated against the measurements from SFMR.
The SAR-derived wind speeds from VH-polarized images and the SAR-derived wind directions from VV-polarized images are used as inputs in the GMF CMOD5.N. Then, the difference between SAR-measured and CMOD5.N-simulated NRCSs are computed. These samples are collocated with the measurements from SFMR in term of radar incidence angle. The relationship between the difference in NRCS and SFMR rain rate is shown in Figure 7, where the color bar represents the distance between the location of SFMR measurements and the TC’s eye. Traditionally, the rainfall dominates at the inner circle of a TC (radius from TC eye less than 100 km) and the rainfall is relative weak at the outer circle of a TC (radius from TC eye greater than 100 km). The difference is divided into two regions: samples at the inner and outer circle. It is found that the difference in NRCS ranges from 0.5 dB to 5 dB (3 dB to 6 dB) when the distance from TC’s eye is lower (larger) than 100 km. This is likely due to the fact that the wind speed is relatively smaller at the outer part of TCs, resulting in the rain cells highly affecting the sea surface roughness. According to these results, the identification of rain cells is implemented for the S-1 SAR dataset, see Figure 8, where they are marked as grey spots. Note that the sharp pattern of rain cells at the edge of S-1 swaths is due to the SAR winds derived from VH-polarized NRCS; see Figure 8d for the Hurricane Delta.
The TC wind fields are composited by VH-polarized SAR-derived winds, VV-polarized SAR-derived winds and radial-vortex modeling winds. Each image is classified in a region distorted by the rain cells, where model-simulated winds are selected, and a rain-free region, where SAR-derived winds are selected. In particular, the VH-polarized SAR-derived winds at low-to-moderate condition (wind speed lower than 25 m s−1) is replaced by the VV-polarized SAR-derived winds. Figure 9 shows wind speed in two radial profiles on the images of TC Irma overlaid by radial-vortex modeling winds at two flow directions: (a) 45° and (b) 135°, in which the red curve is the VH-polarized SAR-derived without rain, the blue curve is with rain and the black line is the model-simulated result. It is found that the pattern of the radial-vortex modeling wind profile performs well. Besides, the sharp gradient of VH-polarized SAR-derived wind in the radial profiles at the flow direction of 135° (see Figure 9b) is clearly observed, which is caused by the noise floor at the edge of swaths, e.g., 100 km radial distance away from TC eye. The composite winds speed maps obtained for the S-1 SAR dataset (see Figure 1) are shown in Figure 10. It can be noted that the pattern of extreme winds around the TC’s eye (e.g., see Figure 10b relevant to the TC Dorian) is clearer than that in the VV-polarized and VH-polarized wind maps (e.g., see Figure 2b and Figure 3b, respectively). In addition, rain bands far away from the TC’s eye (e.g., see Figure 10d relevant to the TC Delta) are also observed. The statistical analysis of composite winds validated against with SFMR measurements in term of flow directions of TC eyes is shown in Figure 11. In general, the accuracy of composited wind speed is improved after rain correction, e.g., the COR is greater than 0.7 at all flow directions. In particular, the RMSE of wind speed is reduced to be 1.8 m s−1 with a 0.72 COR at the flow directions of 0°–30° and 330°–360°. This behavior is probably caused by the heavy rain in front of TC eyes. In this sense, we can conclude that the TC wind fields retrieved from dual-polarized SAR image are more reliable after rain correction.

4. Conclusions

Heavy rainfall, strong winds and extreme ocean waves are effects associated to TCs. Remote sensing tools as the microwave imager of the TRMM are promising for rainfall monitoring of the oceans. S-1 SAR mission operating at C-band has great capabilities providing the sea surface observations over large swath coverage and fine spatial resolution at both co- and cross-polarization. At present, sea winds associated to TCs could be inverted from VV- and VH-polarized SAR imagery by conventional GMFs that describe the semi-empirical relationships between NRCS and the wind field. Nonetheless, since rain cells directly modulate the ocean roughness resulting in changes of the measured NRCS, the wind retrieval accuracy reduces when heavy rainfall associated to TCs is in place.
In this study, a novel approach to wind retrieval from dual-polarized SAR imagery after rain correction is proposed. Four S-1 C-band SAR collected during TCs are considered, which are collocated with the wind speed and rain rate SFMR measurements. The SAR winds are retrieved from VH- and VV-polarized SAR images by S1IW.NR and CMOD5.N GMFs, respectively. The comparison between VH-polarized SAR-derived winds and the measurements of SFMR yields that maximum RMSE of wind speed is about 4 m s−1 with a 0.45 COR ranged from 330° to 360°, whereas the RMSE of wind speed is less than 3 m s−1 at other flow directions. The TC parameters are estimated by minimizing of the cost function, which relates the CMOD5.N-simulated NRCS taking VH-polarized SAR-derived wind speed <20 m s−1 with the SAR-measured NRCS. The retrieval TC intensity are compared with the best-track information of NOAA, indicating that the RMSE of central pressure and maximum wind speed is 15.6 hPa and 10.8 m s−1, respectively. In the meantime, the difference between VV-polarized NRCS measured by the SAR and the NRCS simulated using the CMOD5.N GMF taking VH-polarized SAR winds as input is related with the SFMR rain rate. Rain cells are identified on the basis that the NRCS difference ranges from 0.5 dB to 5 dB (3 dB to 6 dB) when the distance from TC’s eye is lower (larger) than 100 km. Hence, a composite wind field is reconstructed by using dual-polarized SAR-derived winds and winds simulated by the parametric Holland model for rain-free and rain-affected regions, respectively. The validation results of composite winds against with SFMR show the RMSE < 2 m s−1 and COR > 0.7 at all flow direction range, indicating that rain correction provides better performance in TC conditions.
As mentioned in [46], the automatic detection of rain cell necessitates complicated treatment especially for weak rain identification at regular sea state. Future work will include the improvement of rain identification and development of a new rain rate retrieval algorithm by means of a larger SAR dataset collected during TCs.

Author Contributions

Conceptualization, W.S. and Z.L.; methodology, W.S. and Z.L.; validation, W.S. and Z.L.; formal analysis, W.S., X.J. and J.Z.; investigation, W.S. and Z.L.; resources, W.S.; writing—original draft preparation, W.S., Z.L., F.N. and A.B.; writing—review and editing, W.S. and Z.L.; visualization, Z.L.; funding acquisition, W.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China, grant number 42076238, 42176012 and 42130402, and the Shanghai Frontiers Research Center of the Hadal Biosphere.

Data Availability Statement

Not applicable.

Acknowledgments

The S-1 SAR images were provided by ESA via https://scihub.copernicus.eu (accessed on 1 September 2022). The ERA-5 winds were accessed via http://www.ecmwf.int (accessed on 1 September 2022). We also appreciate the provision of NOAA for SFMR data and the information on hurricane best-tracks.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

References

  1. Young, I.R. Observations of the Spectra of Hurricane Generated Waves. Ocean Eng. 1998, 25, 261–276. [Google Scholar] [CrossRef]
  2. Chen, T.C.; Wu, C.C. The Remote Effect of Typhoon Megi (2010) on the Heavy Rainfall over Northeastern Taiwan. Mon. Weather Rev. 2016, 144, 3109–3131. [Google Scholar] [CrossRef]
  3. Koo, M.S.; Hong, S.Y.; Kim, J. Evaluation of the Tropical Rainfall Measuring Mission (TRMM) Multi-satellite Precipitation Analysis (TMPA) Data over South Korea. Asia-Pac. J. Atmos. Sci. 2009, 45, 265–282. [Google Scholar]
  4. Kummerow, C.; Hong, Y.; Olson, W.S.; Yang, S.; Adler, R.F.; McCollum, J.; Ferraro, R.; Petty, G.; Shin, D.B.; Wilheit, T.T. The Evolution of the Goddard Profiling Algorithm (GPROF) for Rainfall Estimation from Passive Microwave Sensors. J. Appl. Meteorol. 2001, 40, 1801–1820. [Google Scholar] [CrossRef]
  5. Shao, W.Z.; Zhao, C.; Jiang, X.W.; Wang, W.L.; Shen, W.; Zuo, J.C. Preliminary Analysis of Wave Retrieval from Chinese Gaofen-3 SAR Imagery in the Arctic Ocean. Eur. J. Remote Sens. 2022, 55, 456–468. [Google Scholar] [CrossRef]
  6. 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] [Green Version]
  7. Zhang, B.; Perrie, W. Cross-polarized Synthetic Aperture Radar: A New Potential Measurement Technique for Hurricanes. Bull. Am. Meteorol. Soc. 2012, 93, 531–541. [Google Scholar] [CrossRef] [Green Version]
  8. 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]
  9. Shao, W.Z.; Ding, Y.Y.; Li, J.C.; Gou, S.P.; Nunziata, F.; Yuan, X.Z.; Zhao, L.B. Wave Retrieval under Typhoon Conditions Using a Machine Learning Method Applied to Gaofen-3 SAR Imagery. Can. J. Remote Sens. 2019, 45, 723–732. [Google Scholar] [CrossRef]
  10. Atlas, D.; Black, P.G. The Evolution of Convective Storms from Their Footprints on the Sea as Viewed by Synthetic Aperture Radar from Space. Bull. Am. Meteorol. Soc. 1994, 75, 1183–1190. [Google Scholar] [CrossRef]
  11. Melsheimer, C.; Alpers, W.; Gade, M. Simultaneous Observations of Rain Cells Over the Ocean by the Synthetic Aperture Radar aboard the ERS Satellites and by Surface-based Weather Radars. J. Geophys. Res.-Oceans 2001, 106, 4665–4677. [Google Scholar] [CrossRef]
  12. Shi, J.; Hu, J.C.; Shao, W.Z.; Wang, X.Q.; Yuan, X.Z.; Zhao, L.B.; Li, X.F. The Impact of Rain to Observed Signal from Chinese Gaofen-3 Synthetic Aperture Radar in Typhoons. Acta Oceanol. Sin. 2019, 38, 121–133. [Google Scholar] [CrossRef]
  13. Fritz, J.P.; Chandrasekar, V. A Fully Polarimetric Characterization of the Impact of Precipitation on Short Wavelength Synthetic Aperture Radar. IEEE Trans. Geosci. Remote Sens. 2012, 50, 2037–2048. [Google Scholar] [CrossRef]
  14. Zhao, X.B.; Shao, W.Z.; Zhao, L.B.; Gao, Y.; Hu, Y.Y.; Yuan, X.Z. Impact of Rain on Wave Retrieval from Sentinel-1 Synthetic Aperture Radar Images in Tropical Cyclones. Adv. Space Res. 2021, 67, 3072–3086. [Google Scholar] [CrossRef]
  15. Nie, C.; Long, D.G. A C-band scatterometer Simultaneous Wind/rain Retrieval Method. IEEE Trans. Geosci. Remote Sens. 2008, 46, 3618–3631. [Google Scholar] [CrossRef]
  16. Xu, F.; Li, X.F.; Wang, P.; Yang, J.S.; Pichel, W.G.; Jin, Y.Q. A Backscattering Model of Rainfall over Rough Sea Surface for Synthetic Aperture Radar. IEEE Trans. Geosci. Remote Sens. 2015, 53, 3042–3054. [Google Scholar] [CrossRef]
  17. Yuan, X.Z.; Shao, W.Z.; Han, B.; Wang, X.C.; Wang, X.Q.; Gao, Y. Rain-Induced Characteristics of the Observed Signals from Synthetic Aperture Radar at the C and X-Bands in Tropical Cyclones. Remote Sens. Lett. 2021, 12, 573–584. [Google Scholar] [CrossRef]
  18. 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] [Green Version]
  19. Yao, R.; Shao, W.Z.; Jiang, X.W.; Yu, T. Wind Speed Retrieval from the Chinese Gaofen-3 Synthetic Aperture Radar Using an Analytical Approach in the Nearshore Waters of China’s seas. Int. J. Remote Sens. 2022, 43, 3028–3048. [Google Scholar] [CrossRef]
  20. 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.-Oceans 1986, 91, 13065–13083. [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.A.; Vogelzang, J.; Verhoef, A. The CMOD7 Geophysical Model Function for ASCAT and ERS Wind Retrievals. IEEE J. Sel. Top. Appl. Earth Observ. Remote Sens. 2017, 10, 2123–2134. [Google Scholar] [CrossRef]
  23. Shao, W.Z.; Nunziata, F.; Zhang, Y.G.; Corcione, V.; Migliaccio, M. Wind Speed Retrieval from the Gaofen-3 Synthetic Aperture Radar for VV- and HH-polarization Using a Re-tuned Algorithm. Eur. J. Remote Sens. 2021, 54, 318–337. [Google Scholar] [CrossRef]
  24. Yang, X.F.; Li, X.F.; Zheng, Q.A.; Gu, X.F.; Pichel, W.G.; 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]
  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. Wang, H.; Li, H.M.; Lin, M.S.; Zhu, J.H.; Wang, J.; Li, W.W.; Cui, L.M. Calibration of the Copolarized Backscattering Measurements from Gaofen-3 Synthetic Aperture Radar Wave Mode Imagery. IEEE J. Sel. Top. Appl. Earth Observ. Remote Sens. 2019, 12, 1748–1762. [Google Scholar] [CrossRef]
  27. 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]
  28. Hwang, P.; Plant, W.J. An Analysis of the Effects of Swell and Surface Roughness Spectra on Microwave Backscatter from the Ocean. J. Geophys. Res.-Oceans 2010, 115, C04014. [Google Scholar] [CrossRef] [Green Version]
  29. Vachon, P.W.; Wolfe, J. C-band Cross-polarization Wind Speed Retrieval. IEEE Geosci. Remote Sens. Lett. 2011, 8, 456–459. [Google Scholar] [CrossRef]
  30. Hwang, P.; Fois, F. Surface Roughness and Breaking Wave Properties Retrieved from Polarimetric Microwave Radar Backscattering. J. Geophys. Res.-Oceans 2015, 120, 3640–3657. [Google Scholar] [CrossRef]
  31. Shao, W.Z.; Li, X.F.; Hwang, P.; Zhang, B.; Yang, X.F. Bridging the Gap Between Cyclone Wind and Wave by C-Band SAR Measurements. J. Geophys. Res.-Oceans 2017, 122, 6714–6724. [Google Scholar] [CrossRef]
  32. Zhang, G.S.; Li, X.F.; Perrie, W.; Hwang, P.A.; Zhang, B.; Yang, X.F. 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]
  33. Zhu, S.; Shao, W.Z.; Marino, A.; Sun, J.; Yuan, X.Z. Semi-empirical Algorithm for Wind Speed Retrieval from Gaofen-3 Quad-polarization Strip Mode SAR Data. J. Ocean Univ. 2020, 19, 23–35. [Google Scholar] [CrossRef]
  34. Shao, W.Z.; Hu, Y.Y.; Yang, J.S.; Nunziata, F.; Sun, J.; Li, H.; Zuo, J.C. 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] [Green Version]
  35. Shen, H.; Seitz, C.; Perrie, W.; He, Y.J.; Powell, M. Developing a Quality Index Associated with Rain for Hurricane Winds from SAR. Remote Sens. 2018, 10, 1783. [Google Scholar] [CrossRef] [Green Version]
  36. Combot, C.; Mouche, A.; Knaff, J.; Zhao, Y.L.; Vinour, L.; Quilfen, Y.; Chapron, B. Extensive high-resolution Synthetic Aperture Radar (SAR) data analysis of Tropical Cyclones: Comparisons with SFMR flights and Best-Track. Mon. Weather Rev. 2020, 148, 4545–4563. [Google Scholar] [CrossRef]
  37. Sapp, J.W.; Alsweiss, S.O.; Jelenak, Z.; Chang, P.S.; Carswell, J. Stepped Frequency Microwave Radiometer Wind-Speed Retrieval Improvements. Remote Sens. 2019, 11, 214. [Google Scholar] [CrossRef] [Green Version]
  38. Stoffelen, A.; Anderson, D. Scatterometer Data Interpretation: Estimation and Validation of the Transfer Function CMOD4. J. Geophys. Res.-Oceans 1997, 102, 5767–5780. [Google Scholar] [CrossRef]
  39. Lu, Y.; Zhang, B.; Perrie, W.; Mouche, A.A.; Li, X.F.; Wang, H. A C-band Geophysical Model Function for Determining Coastal Wind Speed Using Synthetic Aperture Radar. IEEE J. Sel. Top. Appl. Earth Observ. Remote Sens. 2018, 11, 2417–2428. [Google Scholar] [CrossRef] [Green Version]
  40. Alpers, W.; Brümmer, B. Atmospheric Boundary Layer Rolls Observed by the Synthetic Aperture Radar aboard the ERS-1 Satellite. J. Geophys. Res.-Atmos. 1994, 99, 12613–12621. [Google Scholar] [CrossRef]
  41. Lehner, S.; Horstmann, J.; Koch, W.; Rosenthal, W. Mesoscale Wind Measurements Using Recalibrated ERS SAR Images. J. Geophys. Res.-Oceans 1998, 103, 7847–7856. [Google Scholar] [CrossRef]
  42. Shao, W.Z.; Li, X.M.; Lehner, S.; Guan, C.L. Development of Polarization Ratio Model for Sea Surface Wind Field Retrieval from TerraSAR-X HH Polarization Data. Int. J. Remote Sens. 2014, 35, 4046–4063. [Google Scholar] [CrossRef]
  43. Gao, Y.; Guan, C.L.; 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]
  44. Gao, Y.; Sun, J.; Zhang, J.; Guan, C.L. Extreme Wind Speeds Retrieval Using Sentinel-1 IW Mode SAR Data. Remote Sens. 2021, 13, 1867. [Google Scholar] [CrossRef]
  45. Koch, W. Directional Analysis of SAR Images Aiming at Wind Direction. IEEE Trans. Geosci. Remote Sens. 2004, 42, 702–710. [Google Scholar] [CrossRef]
  46. Zhao, Y.; Longepe, N.; Mouche, A.; Husson, R. Automated Rain Detection by Dual-polarization Sentinel-1 Data. Remote Sens. 2021, 13, 3155. [Google Scholar] [CrossRef]
  47. Holland, G.J. An Analytical Model of the Wind and Pressure Profiles in Hurricanes. Mon. Weather Rev. 1980, 108, 1212–1218. [Google Scholar] [CrossRef]
  48. Reppucci, A.; Lehner, S.; Schulz-Stellenfleth, J.; Brusch, S. Tropical Cyclone Intensity Estimated from Wide-swath SAR images. IEEE Trans. Geosci. Remote Sens. 2010, 48, 1639–1649. [Google Scholar] [CrossRef]
  49. Holland, G.J.; Belanger, J.I.; Fritz, A. A Revised Model for Radial Profiles of Hurricane Winds. Mon. Weather Rev. 2010, 138, 4393–4401. [Google Scholar] [CrossRef]
Figure 1. VV-polarized normalized radar cross section (NRCS) images relevant to the Sentinel-1 (S-1) synthetic aperture radar (SAR) dataset collected over the tropical cyclones (TCs): (a) Irma; (b) Dorian; (c) Lsaias and (d) Delta. The red lines represent the he Stepped-Frequency Microwave Radiometer (SFMR) tracks provided by the National Oceanic and Atmospheric Administration aircraft.
Figure 1. VV-polarized normalized radar cross section (NRCS) images relevant to the Sentinel-1 (S-1) synthetic aperture radar (SAR) dataset collected over the tropical cyclones (TCs): (a) Irma; (b) Dorian; (c) Lsaias and (d) Delta. The red lines represent the he Stepped-Frequency Microwave Radiometer (SFMR) tracks provided by the National Oceanic and Atmospheric Administration aircraft.
Remotesensing 14 05006 g001aRemotesensing 14 05006 g001b
Figure 2. Inverted wind maps derived from of VV-polarized S-1 SAR images using the CMOD5.N geophysical model function (GMF). (ad) The SAR scenes shown in Figure 1a–d.
Figure 2. Inverted wind maps derived from of VV-polarized S-1 SAR images using the CMOD5.N geophysical model function (GMF). (ad) The SAR scenes shown in Figure 1a–d.
Remotesensing 14 05006 g002aRemotesensing 14 05006 g002b
Figure 3. Inverted wind maps derived from of VH-polarized S-1 SAR images using the S1IW.NR GMF. (ad) The SAR scenes shown in Figure 1a–d.
Figure 3. Inverted wind maps derived from of VH-polarized S-1 SAR images using the S1IW.NR GMF. (ad) The SAR scenes shown in Figure 1a–d.
Remotesensing 14 05006 g003aRemotesensing 14 05006 g003b
Figure 4. Taylor diagram of the matchups between VH-polarized SAR-derived wind speeds and SFMR observation.
Figure 4. Taylor diagram of the matchups between VH-polarized SAR-derived wind speeds and SFMR observation.
Remotesensing 14 05006 g004
Figure 5. Comparison between SAR-derived TC parameters and SFMR observations. (a) Central pressure and (b) maximum wind speed.
Figure 5. Comparison between SAR-derived TC parameters and SFMR observations. (a) Central pressure and (b) maximum wind speed.
Remotesensing 14 05006 g005
Figure 6. Flowchart of the proposed methodology.
Figure 6. Flowchart of the proposed methodology.
Remotesensing 14 05006 g006
Figure 7. Relationship between the VV-polarized NRCS difference (SAR measurements versus CMOD5.N simulations) and SFMR rain rate. The color bar represents the distance between the location of SFMR measurements and the TC’s center.
Figure 7. Relationship between the VV-polarized NRCS difference (SAR measurements versus CMOD5.N simulations) and SFMR rain rate. The color bar represents the distance between the location of SFMR measurements and the TC’s center.
Remotesensing 14 05006 g007
Figure 8. Maps of rain cell identification: (ad) refer to the SAR scenes shown in Figure 1a–d.
Figure 8. Maps of rain cell identification: (ad) refer to the SAR scenes shown in Figure 1a–d.
Remotesensing 14 05006 g008aRemotesensing 14 05006 g008b
Figure 9. Wind speed in two radial profiles on the images of TC Irma overlaid by radial-vortex modeling winds at two flow directions: (a) 45° and (b) 135°. The red curve is the VH-polarized SAR-derived without rain, the blue curve is with rain and the black line is the model-simulated result.
Figure 9. Wind speed in two radial profiles on the images of TC Irma overlaid by radial-vortex modeling winds at two flow directions: (a) 45° and (b) 135°. The red curve is the VH-polarized SAR-derived without rain, the blue curve is with rain and the black line is the model-simulated result.
Remotesensing 14 05006 g009
Figure 10. Maps of composite SAR-derived winds: (ad) refer to the SAR scenes shown in Figure 1a–d.
Figure 10. Maps of composite SAR-derived winds: (ad) refer to the SAR scenes shown in Figure 1a–d.
Remotesensing 14 05006 g010
Figure 11. Taylor diagram of SAR-derived wind speeds versus SFMR observations.
Figure 11. Taylor diagram of SAR-derived wind speeds versus SFMR observations.
Remotesensing 14 05006 g011
Table 1. Main SAR and ancillary information of the dataset.
Table 1. Main SAR and ancillary information of the dataset.
TCCategoryDate and UTC Time of the SAR
Acquisitions
Pixel Size,
Range × Azimuth (m)
Maximum Central Pressure (hPa)Maximum Wind Speed Radius (km)
Irma407-09-2017
10:30 UTC
10 × 1092135
Dorian530-08-2019
22:46 UTC
10 × 1094954
Lsaias102-08-2020
23:19 UTC
10 × 10995111
Delta408-10-2020
00:07 UTC
10 × 1097373
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

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. https://doi.org/10.3390/rs14195006

AMA Style

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 Sensing. 2022; 14(19):5006. https://doi.org/10.3390/rs14195006

Chicago/Turabian Style

Shao, Weizeng, Zhengzhong Lai, Ferdinando Nunziata, Andrea Buono, Xingwei Jiang, and Juncheng Zuo. 2022. "Wind Field Retrieval with Rain Correction from Dual-Polarized Sentinel-1 SAR Imagery Collected during Tropical Cyclones" Remote Sensing 14, no. 19: 5006. https://doi.org/10.3390/rs14195006

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