Hurricane Wind Speed Estimation Using WindSat 6 and 10 GHz Brightness Temperatures

The realistic and accurate estimation of hurricane intensity is highly desired in many scientific and operational applications. With the advance of passive microwave polarimetry, an alternative opportunity for retrieving wind speed in hurricanes has become available. A wind speed retrieval algorithm for wind speeds above 20 m/s in hurricanes has been developed by using the 6.8 and 10.7 GHz vertically and horizontally polarized brightness temperatures of WindSat. The WindSat measurements for 15 category 4 and category 5 hurricanes from 2003 to 2010 and the corresponding H*wind analysis data are used to develop and validate the retrieval model. In addition, the retrieved wind speeds are also compared to the Remote Sensing Systems (RSS) global all-weather product and stepped-frequency microwave radiometer (SFMR) measurements. The statistical results show that the mean bias and the overall root-mean-square (RMS) difference of the retrieved wind speeds with respect to the H*wind analysis data are 0.04 and 2.75 m/s, respectively, which provides an encouraging result for retrieving hurricane wind speeds over the ocean surface. The retrieved wind speeds show good agreement with the SFMR measurements. Two case studies demonstrate that the mean bias and RMS difference are 0.79 m/s and 1.79 m/s for hurricane Rita-1 and 0.63 m/s and 2.38 m/s for hurricane Rita-2, respectively. In general, the wind speed retrieval accuracy of the new model in hurricanes ranges from 2.0 m/s in light rain to 3.9 m/s in heavy rain.


Introduction
Passive microwave remote sensing is an important tool for studying atmospheric and oceanographic processes and can provide both daytime and nighttime observations of geophysical parameters [1], such as atmospheric water vapor, cloud liquid water, rain rate, sea surface temperature, sea surface wind, and sea ice.This tool can provide scientists and forecasters with vital information for understanding and studying global weather and climate change.The first Special Sensor Microwave Imager (SSM/I) aboard the Defense Meteorological Satellite Program (DMSP) was launched in July 1987 [2], making it possible to obtain the global ocean surface wind from a spaceborne passive microwave instrument [3].Since then, a successive dataset of wind speed has been obtained from the passive microwave radiometry [3][4][5].
Algorithms [3,[6][7][8][9] have been developed for passive microwave radiometry that are able to retrieve the wind speeds of low to moderate winds in no-rain areas with a good degree of accuracy.
However, these algorithms seem to break down completely when the scenes are under severe weather conditions, including both high winds and precipitation.It is difficult to measure wind speeds in these conditions because the physics of ocean surface remote sensing is still poorly understood under these weather conditions.Intense rain not only influences the atmospheric attenuation but also changes the sea surface dielectric properties in a complicated manner [10,11].Thus, it is very difficult to accurately model brightness temperatures in these situations.Fortunately, the brightness temperature (TB) acquired at the lower frequencies, such as the L-band and C-band, are far from saturation, which is the physical reason that researchers can "see" the ocean surface and derive its properties [11,12].A new capability for ocean surface wind speed retrieval is performed using the L-band brightness temperatures provided by the Soil Moisture and Ocean Salinity (SMOS) mission [11].The study [11] shows that the L-band ocean emissivity dependence with wind speed appears to be less sensitive to roughness and foam changes than at the higher C-band microwave frequencies.A physical method is developed to retrieve sea surface temperature and wind speed for hurricane predictions in [13].In [13], the effects of precipitation emission and scattering on the measurements are both taken into account.However, the used ocean microwave emissivity model is performed in clear sky conditions and at low wind speeds (smaller than 25 m/s).Hence, the performance of the model needs to be further tested at high wind speeds and intensely rainy conditions.In addition, the C-band H-polarization TBs provided by AMSR-E and AMSR2 are used to retrieve wind speed in hurricanes using the small slope approximation method and an ocean surface roughness spectrum [14].In [14], the retrieved wind speeds are vulnerable to rain because the model does not account for the effects of absorption and scattering caused by rain.
It is very difficult to accurately model brightness temperatures in rain because of the high variability in rainy atmospheres.Thus, the statistical algorithm seems to be a good candidate for retrieval of wind speeds in this situation.For accurate radiometer retrievals of wind speeds under rainy conditions, channel combination methods at different frequencies have been used to retrieve wind speeds [10,12,15,16].The core of these algorithms makes it possible to reduce the signal coming from rain without considerably reducing the wind speed signal.Such a technique has been employed successfully for wind speed retrieval with the airborne microwave polarimetry.In addition, an artificial neural network technique has been used to derive the relationship between the AMSR2 TBs and sea surface wind speeds in hurricanes [17][18][19].However, the retrieved results are mainly compared with the in situ surface wind speed observations below 30 m/s in [17][18][19].They are not validated using the high wind measurements.
In this study, a new model for retrieving high wind speeds in hurricanes using WindSat C-band and X-band vertically and horizontally polarized brightness temperatures is proposed.The core of the new model is also based on the channel combination at different frequencies.Referring to [15], we define two new parameters, W6H and W6V, which represent an increment at 6H and 6V due to ocean wind, respectively.Finally, a wind speed retrieval algorithm for wind speeds above 20 m/s in hurricanes has been developed by using the two new parameters.The WindSat measurements for 15 category 4 and category 5 hurricanes from 2003 to 2010 and the corresponding H*wind analysis data are used to develop and validate the retrieval model.This paper is organized as follows.In Section 2, the datasets, including WindSat TBs, H*wind analysis data, SFMR measurements, RSS global all-weather product and rain rate, are introduced.The methodology is presented in Section 3. Section 4 gives detailed results.The last section presents the conclusions.

Datasets
In this paper, WindSat TB measurements over hurricanes, and H*wind analysis data are used to develop and test a new wind speed retrieval algorithm for inside hurricanes.The SFMR along-track wind speed measurements and the RSS global all-weather wind speed product are compared with the retrieved wind speeds.In addition, the rain effect is also evaluated by using the RSS rain rate product.The selected datasets are described in the following sections.

Wind Sat Brightness Temperature
WindSat is a conically scanning passive microwave radiometer aboard the U.S. Department of Defense Coriolis satellite, which orbits the Earth in an 840-km circular sun-synchronous orbit.Its primary objective is to demonstrate the capability of a polarimetric microwave radiometer to measure the ocean surface wind vector from space.WindSat operates in 5 channels: 6.8, 10.7, 18.7, 23.8, and 37 GHz.All are fully polarimetric except for the 6.8 and 23.8 GHz channels, which feature dual polarization (vertical and horizontal).The integrated field of view (IFOV) is 40 × 60 km, 25 × 38 km, 16 × 27 km, 12 × 20 km, and 8 × 13 km for the 6.8, 10.7, 18.7, 23.8, and 37 GHz channels, respectively.WindSat takes TB observations during both the forward-looking and the aft-looking scans.The forward-looking swath is approximately 1000 km and the aft-looking swath is approximately 400 km.During each scan, every channel is of a different location and resolution.Hence, the TB measurements must be resampled and averaged to a common size for wind vector retrieval.The selected TB measurements have been averaged to the 6.8 GHz footprint size with the WindSat ground data processing version 2.0 algorithm.A detailed description of the WindSat instrument can be found in [20].In this paper, we only use the vertically and horizontally polarized measurements for wind speed retrieval.

H*wind Products
Validation of hurricane surface wind fields is inherently challenging because there is no realizable "surface truth" data.Fortunately, a relatively reliable and reasonable estimation of tropical cyclone (TC) wind fields has been developed through the H*wind analysis system.The wind fields produced by this system have been extensively used to study the wind speed retrieval algorithms under severe weather conditions [10,11,15,[21][22][23][24].H*wind is a software tool that provides an objective analysis of the 1-min sustained-wind field at a height of 10 m above the sea surface and was developed by the Hurricane Research Division (HRD) of the Atlantic Oceanographic and Meteorological Laboratory [25].The H*wind analysis system assimilates all available surface observations, such as buoys, ships, C-MAN platforms, and surface airways, and all available aircraft and remote sensing data sources, which typically includes the SFMR, QuikSCAT, SSM/I, Geostationary Operational Environmental Satellite (GOES) low-level cloud drift winds, and Global Positioning System (GPS) dropwindsondes from aircraft.Typical wind speed errors in an H*wind analysis are estimated to be 10% to 20% [26].The wind speed errors in this analysis will vary depending on the quantity and quality of data that are available, as well as on the degree of quality control employed by the analyst.Although each wind speed analysis produced by this system may be inaccurate, the ensemble average over a large number of collocations can minimize the error in the overall analysis.
A hurricane is defined as storm with wind speeds exceeding 65 knots (33.5 m/s) [15].1.These data can be obtained from the Automated Tropical Cyclone Forecast system [27], with products available at the National Oceanic Atmospheric Administration (NOAA) National Hurricane Center (NHC).

Airborne Stepped-Frequency Microwave Radiometer (SFMR) Measurements
An SFMR is an airborne remote sensing instrument designed for retrieving surface wind speed and rainfall in TCs [28].The first operational SFMR was flown on the WP-3D aircraft in 2004 [29].An SFMR observes the nadir brightness temperatures at 6 C-band frequencies, i.e., 4.55, 5.06, 5.64, 6.34, 6.96, and 7.22 GHz.The TB measurements are processed with a 10-second running mean to retrieve ocean surface wind speed using a geophysical model function that relates the surface emissivity to wind speed.Global Positioning System dropwindsonde measurements have been used to validate the SFMR wind speed estimates, and the RMS error is approximately 4 m/s for measured wind speeds ranging from 10 to 70 m/s [29].A previous study [30] shows that the SFMR wind speed in heavy precipitation has been further improved by using the new relationship between microwave absorption and rain rate.The advantage of SFMR is that it can potentially retrieve ocean surface wind speed at a relatively high spatial resolution (1.5 km).However, it can only provide point measurements along the aircraft flight track.

Other Sources
Using WindSat TB measurements, a global all-weather wind speed retrieval algorithm has been established by RSS [10].It is a semi-statistical algorithm that can be applied in all rainy atmospheres over the whole Earth.The estimated RMS accuracy of this algorithm ranges from approximately 2.0 m/s in light rain to approximately 5.0 m/s in heavy rain, if C-band channels and higher frequencies

Airborne Stepped-Frequency Microwave Radiometer (SFMR) Measurements
An SFMR is an airborne remote sensing instrument designed for retrieving surface wind speed and rainfall in TCs [28].The first operational SFMR was flown on the WP-3D aircraft in 2004 [29].An SFMR observes the nadir brightness temperatures at 6 C-band frequencies, i.e., 4.55, 5.06, 5.64, 6.34, 6.96, and 7.22 GHz.The TB measurements are processed with a 10-second running mean to retrieve ocean surface wind speed using a geophysical model function that relates the surface emissivity to wind speed.Global Positioning System dropwindsonde measurements have been used to validate the SFMR wind speed estimates, and the RMS error is approximately 4 m/s for measured wind speeds ranging from 10 to 70 m/s [29].A previous study [30] shows that the SFMR wind speed in heavy precipitation has been further improved by using the new relationship between microwave absorption and rain rate.The advantage of SFMR is that it can potentially retrieve ocean surface wind speed at a relatively high spatial resolution (1.5 km).However, it can only provide point measurements along the aircraft flight track.

Other Sources
Using WindSat TB measurements, a global all-weather wind speed retrieval algorithm has been established by RSS [10].It is a semi-statistical algorithm that can be applied in all rainy atmospheres over the whole Earth.The estimated RMS accuracy of this algorithm ranges from approximately 2.0 m/s in light rain to approximately 5.0 m/s in heavy rain, if C-band channels and higher frequencies are available [10].It is a good candidate for comparison with our retrieved wind speed results using WindSat TB measurements because the sufficient matchups can be obtained from them.We obtain the latest RSS v7.0.1 all-weather wind speed products from the RSS website (http://www.remss.com/).
To analyze the effect of rain, the RSS WindSat v7.0.1 rain rate product is used to directly study the wind speed retrieval errors.The other independent rain rate products are not selected because the time difference may result in a spatial mismatch [31].For example, a TC with a forward motion of 10 m/s could travel 24 km in 40 min, comparable to the mean spatial resolution of the selected TB.In addition, the distribution of rainfall inside TCs is complex and changeable.To some extent, the mismatch could introduce some errors.We obtain the rain rate product from the RSS website (http://www.remss.com/),and the product was developed using the Unified Microwave Ocean Retrieval Algorithm (UMORA) [32].

Matchup
For training and testing the wind speed retrieval algorithm in TCs, we collected TB measurements over the H*wind analysis fields.For WindSat, we collected data from 15 hurricanes during the period from 2003 to 2010.The time of the H*wind analysis data and the average time of the WindSat overpass are within 3 h.No temporal interpolation is performed as the time intervals are too large in most cases.Within the 3-h time window between the products, the hurricane could have moved over a relatively large distance in some cases, which could introduce significant error.Thus, we shift the H*wind analysis data for each TC based on the "best track" fixes from NHC so that the eye of the H*wind analysis data coincides with the eye of the WindSat measurements.In addition, the H*wind analysis fields are performed at a high resolution, which is approximately 5 km [25].
We resample the H*wind analysis data to make them comparable to the low-resolution (6.9 GHz) WindSat TB measurements.Using a nearest neighbor interpolation, all H*wind data points within a radius of 30 km (around the size of the WindSat 6.8 GHz IFOV) of the WindSat TB locations are weighted, inversely proportional to distance, and then averaged: The function of weight is where the "dis" is the distance (in kilometers) from an H*wind point to a WindSat grid point.The form of weighting function is selected because it is based on the fact that the power that is received by the radiometer falls off roughly exponentially with distance and that wind speed is roughly proportional to the brightness temperature (power).
A previous study has shown that the radiometer TB is directly sensitive to surface roughness rather than wind speed [10].The H*wind analysis reports peak wind speeds for 1-min sustained winds.Thus, we need to convert from 1-min sustained H*wind data to 10 min sustained satellite wind speeds because the surface winds for producing the surface roughness of the satellite measurements need to be sustained a longer time, which is at least 10 min.A derived scale factor of 0.88 is used here and by the U.S. Navy [33].
Table 1 lists the statistics of the selected TCs in the datasets, which includes the TC name, year, maximum wind speed, the latitudinal and longitudinal shifts of the TC eye, and the population of the selected wind speed above 20 m/s.Finally, the collocations contain 26,242 matchups for the WindSat.

Polarimetric Microwave Radiation
The Stokes parameters are a set of values that characterize the polarization state of the electromagnetic radiation.Four parameters, I, Q, U, and V, were defined by George Gabriel Stokes in 1852.Since passive microwave radiometers usually provide vertically and horizontally polarized brightness temperatures (T V and T H ) to observe the Earth, an alternate Stokes vector can be represented as T V , T H , T 3 , and T 4 : where T V and T H describe the brightness temperatures of vertical and horizontal polarizations, respectively.T 3 and T 4 are the cross-correlation terms between these two orthogonal polarizations.In addition, c is a constant relating the brightness temperature to the electric energy density [34,35].
In recent years, various geophysical parameters, such as the wind vector, sea surface temperature, atmospheric water vapor, cloud liquid water, rain rate, and soil moisture, have been estimated by using the polarimetric brightness temperatures [8,9,36,37].The brightness temperatures measured by the spaceborne polarimetric radiometer are the sum of the upwelling atmospheric radiation (T BU ), the reflected cosmic background radiation (T BC ) and the downwelling atmospheric radiation (T BD ), and the emission of the ocean surface.The atmosphere attenuates the emission of the ocean surface and the reflected downwelling radiation.Thus, the four Stokes parameters received by a spaceborne polarimetric radiometer at each frequency can be expressed as follows [34]: where e p refers to the ocean surface emissivity at polarization p (p = v, h), and the corresponding ocean surface reflectivity is r p = 1 − e p .T S is the sea surface temperature, and τ is the transmissivity of the atmosphere.The Ω p (p = v, h) term is a correction factor to account for nonspecular reflection of the T BD from the rough ocean surface [38].Figure 2 describes the brightness temperatures received by a radiometer.
Remote Sens. 2016, 8, 721 7 of 20 ( ) ( ) where ep refers to the ocean surface emissivity at polarization p (p = v, h), and the corresponding ocean surface reflectivity is rp = 1 − ep.TS is the sea surface temperature, and is the transmissivity of the atmosphere.The (p = v, h) term is a correction factor to account for nonspecular reflection of the TBD from the rough ocean surface [38].Figure 2 describes the brightness temperatures received by a radiometer.

Channel Selection
Under clear weather conditions, it is common to select higher frequencies such as 36 GHz rather than lower frequencies for retrieving the low to moderate ocean surface winds because the sensitivity to ocean surface wind at higher frequencies is better than that at lower frequencies [15,39].However, the higher frequencies are almost invalid for retrieval of high winds in the severe weather conditions because the atmospheric opacity for higher frequencies exceeds that for lower frequencies.Figure 3 shows the WindSat brightness temperatures for hurricane Isabel on 10 September 2003, along with the best track (piecewise red curve) analysis from the NHC.The top three panels of Figure 3 show the results of horizontally polarized brightness temperature for the 10.7, 18.7, and 37 GHz channels, respectively.Note that we use a different color scale for the 10.7, 18.7, and 37 GHz images to highlight similar spatial features with high radiances, which are apparent signatures of rain.A rain band extending from the center to the northeast can be clearly found in these panels.The 18.7 and 37 GHz channels are nearly saturated close to the eye of the hurricane because of the influence of precipitation.However, the maximum value of the 10.7 GHz channel near the eye of the hurricane is approximately 210 K, which is fairly far from saturation.Similar characteristics can be found at lower frequency, such as 6.8 GHz.In addition, similar results can also be found in the vertically polarized brightness temperatures in the bottom panels of Figure 3. Therefore, we can use the lower frequencies, such as the 6.8 and 10.7 GHz channels, to develop an algorithm for retrieval of high winds under severe weather conditions.

Channel Selection
Under clear weather conditions, it is common to select higher frequencies such as 36 GHz rather than lower frequencies for retrieving the low to moderate ocean surface winds because the sensitivity to ocean surface wind at higher frequencies is better than that at lower frequencies [15,39].However, the higher frequencies are almost invalid for retrieval of high winds in the severe weather conditions because the atmospheric opacity for higher frequencies exceeds that for lower frequencies.Figure 3 shows the WindSat brightness temperatures for hurricane Isabel on 10 September 2003, along with the best track (piecewise red curve) analysis from the NHC.The top three panels of Figure 3 show the results of horizontally polarized brightness temperature for the 10.7, 18.7, and 37 GHz channels, respectively.Note that we use a different color scale for the 10.7, 18.7, and 37 GHz images to highlight similar spatial features with high radiances, which are apparent signatures of rain.A rain band extending from the center to the northeast can be clearly found in these panels.The 18.7 and 37 GHz channels are nearly saturated close to the eye of the hurricane because of the influence of precipitation.However, the maximum value of the 10.7 GHz channel near the eye of the hurricane is approximately 210 K, which is fairly far from saturation.Similar characteristics can be found at lower frequency, such as 6.8 GHz.In addition, similar results can also be found in the vertically polarized brightness temperatures in the bottom panels of Figure 3. Therefore, we can use the lower frequencies, such as the 6.8 and 10.7 GHz channels, to develop an algorithm for retrieval of high winds under severe weather conditions.

Wind Speed Retrieval Model
Previous studies have shown that the sensitivity to the ocean wind at horizontal polarization is larger than that at vertical polarization [7,40].The study of [15] only used 6.9 and 10.7 GHz horizontal polarization (hereafter, 6H and 10H) data from AMSR-E to retrieve wind speed inside hurricanes.However, we retain the 6.9 and 10.7 GHz vertical polarization (hereafter, 6 V and 10 V) data in addition to the horizontal polarization data in this study.We find that the vertical polarization data also shows the same characteristics as the horizontal polarization data described in [15].In this study, we provide an improved algorithm for retrieving wind speed inside TCs by using 6 V(H) and 10V(H) brightness temperatures.Based on the study [15], we also assess the wind-induced emissivity changes at both 6 V(H) and 10 V(H) because we retrieve wind speed by combining data from these channels.
Firstly, the brightness temperatures are simulated using Equations ( 4) and (5).Generally, three processes should be taken into account: absorption, emission, and scattering.The atmospheric absorption includes three components: water vapor, oxygen, and liquid water in the form of clouds or rain [7].The above components also emit microwave energy according to the ambient temperature.Meanwhile, scattering must be taken into account unless the microwave wavelength is much larger than the size of the raindrops.Previous studies have shown that the raindrop-induced scattering can be neglected for frequencies below 10 GHz when the rain rate is less than 12 mm/h [15,41].Scattering has been neglected in our simulation because the selected profiles are in almost clear weather conditions.
The 60-level ECMWF model outputs are used to simulate the brightness temperatures [42] and were generated by the ERA-40 assimilation system, which relies on the 3-dimensional variational scheme described in [43].The sampled database gathers profiles corresponding to various ocean conditions, and the highest surface pressure is 1049 hPa.The wind speed ranges from 0.15 to 21.43 m/s, and the sea surface temperature ranges from 278.34 to 305.02 K.In addition, the water vapor ranges from 4 to 55.64 mm, and the liquid water ranges from 0 to 0.19 mm.The atmospheric transmittance at different channels is computed from the standard radiative transfer equations [39,44].Likewise, the downwelling and upwelling atmospheric brightness temperatures are also calculated from the standard radiative transfer equations [39,44].Note that the oxygen, water vapor, and cloud liquid water absorption coefficient use O2-MPM93 expression [45], H2O-PWR98

Wind Speed Retrieval Model
Previous studies have shown that the sensitivity to the ocean wind at horizontal polarization is larger than that at vertical polarization [7,40].The study of [15] only used 6.9 and 10.7 GHz horizontal polarization (hereafter, 6H and 10H) data from AMSR-E to retrieve wind speed inside hurricanes.However, we retain the 6.9 and 10.7 GHz vertical polarization (hereafter, 6V and 10V) data in addition to the horizontal polarization data in this study.We find that the vertical polarization data also shows the same characteristics as the horizontal polarization data described in [15].In this study, we provide an improved algorithm for retrieving wind speed inside TCs by using 6V(H) and 10V(H) brightness temperatures.Based on the study [15], we also assess the wind-induced emissivity changes at both 6V(H) and 10V(H) because we retrieve wind speed by combining data from these channels.
Firstly, the brightness temperatures are simulated using Equations ( 4) and (5).Generally, three processes should be taken into account: absorption, emission, and scattering.The atmospheric absorption includes three components: water vapor, oxygen, and liquid water in the form of clouds or rain [7].The above components also emit microwave energy according to the ambient temperature.Meanwhile, scattering must be taken into account unless the microwave wavelength is much larger than the size of the raindrops.Previous studies have shown that the raindrop-induced scattering can be neglected for frequencies below 10 GHz when the rain rate is less than 12 mm/h [15,41].Scattering has been neglected in our simulation because the selected profiles are in almost clear weather conditions.
The 60-level ECMWF model outputs are used to simulate the brightness temperatures [42] and were generated by the ERA-40 assimilation system, which relies on the 3-dimensional variational scheme described in [43].The sampled database gathers profiles corresponding to various ocean conditions, and the highest surface pressure is 1049 hPa.The wind speed ranges from 0.15 to 21.43 m/s, and the sea surface temperature ranges from 278.34 to 305.02 K.In addition, the water vapor ranges from 4 to 55.64 mm, and the liquid water ranges from 0 to 0.19 mm.The atmospheric transmittance at different channels is computed from the standard radiative transfer equations [39,44].Likewise, the downwelling and upwelling atmospheric brightness temperatures are also calculated from the standard radiative transfer equations [39,44].Note that the oxygen, water vapor, and cloud liquid water absorption coefficient use O2-MPM93 expression [45], H2O-PWR98 expression [46], and clw-MPM93 expression [45], respectively.In addition, the ocean surface emissivity and the reflectivity are calculated with the RTTOV radiative transfer model [47], and the correction factors are derived from RTTOV v11.3.The cosmic background radiation is approximately 2.7 K.
To account for wind speed in ocean emission, one parameter is defined as 6(10)P* (P = V, H), which represents the difference between the WindSat 6(10)P and calm ocean emission, corrected for atmospheric effects.This parameter represents the quantity of ocean microwave emission changed by ocean wind at 6(10)P, which is defined by Equation (8). 6(10) P * = WindSat_6 (10) P − Atmos_6 (10) P − CalmOcean_6 (10) P (8) CalmOcean_6 (10 where the WindSat_6(10)P is the brightness temperature measurement for the 6( 10) GHz channel at P-polarization (P = V, H) and CalmOcean_6(10)P is the ocean emission for 6(10) GHz channel under calm ocean conditions.The SST used here is from ECMWF.In addition, Atmos_6(10)P is calculated by referring to [48]. Figure 4a,b illustrate the relationships between 6(10)V* and ECMWF wind speed.The 6(10)V* decreases initially because the shielding effect has been found in the calculation of ocean emissivity at V-polarization and then increases with increasing wind speed.The sensitivity of 6(10)V* to wind speed is approximately 0.37(0.5)K/(m/s) when the wind speed is from 8 to 20 m/s. Figure 4d,e show the relationships between 6(10)H* and ECMWF wind speed.The 6(10)H* increases with increasing wind speed, which is different from 6(10)V*.The sensitivity of 6(10)H* to wind speed is better than that of 6(10)V*, and it is approximately 0.85(0.93)K/(m/s).The previous study [15] showed that the sensitivity of 6H* to wind speed is approximately 1 K/(m/s), which is higher than our calculated result.The diverse data resources may result in this difference.On the contrary, our calculated result of 6H* is higher than the experimental result [49], which indicates that the sensitivity at approximately 6-7H is approximately 0.5 K/(m/s), even though the experiments are limited to ground-based observations.Figure 4c illustrate the relationships between the 6V* and 10V*.When the wind speed is from 0 to 20 m/s, the ratio of 6V* to 10V* is a constant which is approximately 0.82.The ratio of 6H* and 10H* has similar characteristics (Figure 4f), and the value is approximately 0.92.These results are in good agreement with the previous results [15,28].
expression [46], and clw-MPM93 expression [45], respectively.In addition, the ocean surface emissivity and the reflectivity are calculated with the RTTOV radiative transfer model [47], and the correction factors are derived from RTTOV v11.3.The cosmic background radiation is approximately 2.7 K.
To account for wind speed in ocean emission, one parameter is defined as 6(10)P* (P = V, H), which represents the difference between the WindSat 6(10)P and calm ocean emission, corrected for atmospheric effects.This parameter represents the quantity of ocean microwave emission changed by ocean wind at 6(10)P, which is defined by Equation ( 8).where the WindSat_6(10)P is the brightness temperature measurement for the 6( 10) GHz channel at P-polarization (P = V, H) and CalmOcean_6(10)P is the ocean emission for 6(10) GHz channel under calm ocean conditions.The SST used here is from ECMWF.In addition, Atmos_6(10)P is calculated by referring to [48]. Figure 4a,b illustrate the relationships between 6(10)V* and ECMWF wind speed.The 6(10)V* decreases initially because the shielding effect has been found in the calculation of ocean emissivity at V-polarization and then increases with increasing wind speed.The sensitivity of 6(10)V* to wind speed is approximately 0.37(0.5)K/(m/s) when the wind speed is from 8 to 20 m/s.Figures 4d and 4e show the relationships between 6(10)H* and ECMWF wind speed.The 6(10)H* increases with increasing wind speed, which is different from 6(10)V*.The sensitivity of 6(10)H* to wind speed is better than that of 6(10)V*, and it is approximately 0.85(0.93)K/(m/s).The previous study [15] showed that the sensitivity of 6H* to wind speed is approximately 1 K/(m/s), which is higher than our calculated result.The diverse data resources may result in this difference.On the contrary, our calculated result of 6H* is higher than the experimental result [49], which indicates that the sensitivity at approximately 6-7H is approximately 0.5 K/(m/s), even though the experiments are limited to ground-based observations.Figure 4c illustrate the relationships between the 6V* and 10V*.When the wind speed is from 0 to 20 m/s, the ratio of 6V* to 10V* is a constant which is approximately 0.82.The ratio of 6H* and 10H* has similar characteristics (Figure 4f), and the value is approximately 0.92.These results are in good agreement with the previous results [15,28].Making use of the fact that the ratio between 6H* and 10H* is a constant, Shibata [15] established a wind speed retrieval algorithm by assuming that this ratio still holds at higher wind speeds and under rainy conditions.Based on the Shibata method and the above discussions, we develop an improved algorithm with two new parameters known as W6V and W6H.The definitions of W6V and W6H are similar to the previous study [15]: 6 (10) where the WindSat_6(10)P and CalmOcean_6(10)P are the same as Equations ( 8) and ( 9).The definition of other parameters in Equations ( 10) to ( 13) can be found in Figure 5.The brightness temperatures of group A (Figure 5a) are the results under calm ocean conditions.The brightness temperatures of group B (Figure 5a) are the results under rough ocean conditions.In Figure 5a, the relationship between 6H and 10H is slightly nonlinear, but we have approximated it as linear for simplicity.By calculating 6(10) H¯, we can obtain Figure 5b, in which line A is approximated as a linear line passing through the point O(a 1 ,b 1 ) with the slope c 1 .Point F in Figure 5b corresponds to an arbitrary point on line B in Figure 5a.Then, we can obtain an intersecting point E made by two lines OP and EF in Figure 5b.The slope of the line EF is sl 1 .Finally, we can obtain the parameter "W6H" (unit: Kelvin) based on the length of EF* divided by the atmospheric effect (=fac 1 ).A detailed description can be found in [15].The parameters for V-polarization in Equations ( 14) to ( 17) are similar to those of H-polarization in Equations ( 10) to (13).
Finally, we propose a new algorithm that relates W6V and W6H to wind speed (m/s): The coefficients a 1 ~f1 and a 2 ~f2 in Equations ( 10) to ( 17) are listed in Table 2.

Brightness Temperature 6(10)P¯ Comparisons
Before calculating the fitted coefficients of the model, we must estimate the 6(10) P¯ (P = V, H) with Equation (13) and Equation (17).We first need to calculate the ocean TBs under the calm sea conditions, which are a function of frequency, polarization, incidence angle, SST and salinity.The SST derived from the Reynolds weekly analysis was used to estimate the ocean microwave emission at 6(10) GHz under the calm sea conditions with the Klein and Swift model [50].For the salinity effect, the previous study [15] demonstrated that the sensitivity of TBs at 6.8 GHz H-polarization is 0.003 K/PSU at 0 °C SST and −0.039 K/PSU at 30 °C SST for a salinity range of 30 to 35 PSU.The values for 10.7 GHz H-polarization are 0.006 K/PSU and −0.012 K/PSU, respectively.Thus, the salinity effect is also negligible in our study.
Figure 6 shows the brightness temperatures of 6H¯ (V¯) and 10H¯ (V¯) inside and around hurricanes.The cases in which the center of the hurricane contains land areas were removed because 6H¯ (V¯) and 10H¯ (V¯) cannot be evaluated over land.Figure 6a shows the WindSat 6H¯ and 10H¯ from hurricane Fabian on 3 September 2003, and the maximum wind speed is 47 m/s. Figure 6b displays the WindSat 6H¯ and 10H¯ on 31 August 2004, with a maximum wind speed of 51 m/s for hurricane Frances. Figure 6c presents the WindSat 6H¯ and 10H¯ on 9 July 2005, with a maximum wind speed of 51 m/s for hurricane Dennis; Figure 6d-f illustrate the results of the vertically polarized channels.In Figure 6a, the data marked by arrow A correspond to the outer area around Fabian, and the data marked by B correspond to the area inside Fabian.The situation is similar for the other cases.In addition, the arrow B moves upward as the wind speed increases, and the above results are in good agreement with the previous study [15].

Brightness Temperature 6(10)P¯Comparisons
Before calculating the fitted coefficients of the model, we must estimate the 6(10) P¯(P = V, H) with Equation (13) and Equation (17).We first need to calculate the ocean TBs under the calm sea conditions, which are a function of frequency, polarization, incidence angle, SST and salinity.The SST derived from the Reynolds weekly analysis was used to estimate the ocean microwave emission at 6(10) GHz under the calm sea conditions with the Klein and Swift model [50].For the salinity effect, the previous study [15] demonstrated that the sensitivity of TBs at 6.8 GHz H-polarization is 0.003 K/PSU at 0 • C SST and −0.039 K/PSU at 30 • C SST for a salinity range of 30 to 35 PSU.The values for 10.7 GHz H-polarization are 0.006 K/PSU and −0.012 K/PSU, respectively.Thus, the salinity effect is also negligible in our study.
Figure 6 shows the brightness temperatures of 6H¯(V¯) and 10H¯(V¯) inside and around hurricanes.The cases in which the center of the hurricane contains land areas were removed because 6H¯(V¯) and 10H¯(V¯) cannot be evaluated over land.Figure 6a shows the WindSat 6H¯and 10Hf rom hurricane Fabian on 3 September 2003, and the maximum wind speed is 47 m/s. Figure 6b displays the WindSat 6H¯and 10H¯on 31 August 2004, with a maximum wind speed of 51 m/s for hurricane Frances. Figure 6c presents the WindSat 6H¯and 10H¯on 9 July 2005, with a maximum wind speed of 51 m/s for hurricane Dennis; Figure 6d-f illustrate the results of the vertically polarized channels.In Figure 6a, the data marked by arrow A correspond to the outer area around Fabian, and the data marked by B correspond to the area inside Fabian.The situation is similar for the other cases.In addition, the arrow B moves upward as the wind speed increases, and the above results are in good agreement with the previous study [15].

Wind Speed Retrieval and Validation
The new model is established by using 4862 matchups between WindSat TBs and H*wind analysis data for Fabian (2003), Frances (2004), and Dennis (2005).We have checked that other groups of hurricanes give similar fitted coefficients for wind speed retrieval.The results show that the mean bias and RMS difference of the new model (Figure 7a) are 0.07 m/s and 2.63 m/s, respectively.Furthermore, we also compare the new model to the model of Shibata (2006) model with the old coefficients [15] and the new coefficients derived from our matchups, respectively.The performance of the new model is improved.The mean bias and RMS difference of the Shibata method with the new coefficients (Figure 7b) are −0.22 m/s and 2.81 m/s, respectively, which are better than those retrieved with the Shibata model using the old coefficients that were derived from the simulations (with mean bias and RMS difference of 1.23 m/s and 5.13 m/s, respectively).In addition, the error increases with increasing wind speed for both the new model and the Shibata model.Through the above discussion, we can conclude the added V-polarization channels contribute to improving the retrieval performance.

Wind Speed Retrieval and Validation
The new model is established by using 4862 matchups between WindSat TBs and H*wind analysis data for Fabian (2003), Frances (2004), and Dennis (2005).We have checked that other groups of hurricanes give similar fitted coefficients for wind speed retrieval.The results show that the mean bias and RMS difference of the new model (Figure 7a) are 0.07 m/s and 2.63 m/s, respectively.Furthermore, we also compare the new model to the model of Shibata (2006) model with the old coefficients [15] and the new coefficients derived from our matchups, respectively.The performance of the new model is improved.The mean bias and RMS difference of the Shibata method with the new coefficients (Figure 7b) are −0.22 m/s and 2.81 m/s, respectively, which are better than those retrieved with the Shibata model using the old coefficients that were derived from the simulations (with mean bias and RMS difference of 1.23 m/s and 5.13 m/s, respectively).In addition, the error increases with increasing wind speed for both the new model and the Shibata model.Through the above discussion, we can conclude the added V-polarization channels contribute to improving the retrieval performance.

Wind Speed Retrieval and Validation
The new model is established by using 4862 matchups between WindSat TBs and H*wind analysis data for Fabian (2003), Frances (2004), and Dennis (2005).We have checked that other groups of hurricanes give similar fitted coefficients for wind speed retrieval.The results show that the mean bias and RMS difference of the new model (Figure 7a) are 0.07 m/s and 2.63 m/s, respectively.Furthermore, we also compare the new model to the model of Shibata (2006) model with the old coefficients [15] and the new coefficients derived from our matchups, respectively.The performance of the new model is improved.The mean bias and RMS difference of the Shibata method with the new coefficients (Figure 7b) are −0.22 m/s and 2.81 m/s, respectively, which are better than those retrieved with the Shibata model using the old coefficients that were derived from the simulations (with mean bias and RMS difference of 1.23 m/s and 5.13 m/s, respectively).In addition, the error increases with increasing wind speed for both the new model and the Shibata model.Through the above discussion, we can conclude the added V-polarization channels contribute to improving the retrieval performance.For hurricane Bill, in particular, the new model largely overestimates wind speeds, with a mean bias of 2.89 m/s.On the contrary, the largest negative bias (−1.45 m/s) is found in hurricane Bertha.On the whole, the retrieved results have a good agreement with the H*wind analysis wind speed.
It is challenging to retrieve wind speeds under severe weather conditions, especially in typhoons and hurricanes.Rain not only increases the atmospheric attenuation but also changes the sea surface roughness in a complicated manner.Thus, it is difficult to accurately model brightness temperatures under rainy conditions.Hence, the statistical algorithm is a relatively good candidate that can be used to retrieve wind speed under severe weather conditions.In general, the new wind speed retrieval algorithm under rain gives an encouraging accuracy compared to the H*wind wind speed.
hurricane Bill, in particular, the new model largely overestimates wind speeds, with a mean bias of 2.89 m/s.On the contrary, the largest negative bias (−1.45 m/s) is found in hurricane Bertha.On the whole, the retrieved results have a good agreement with the H*wind analysis wind speed.
It is challenging to retrieve wind speeds under severe weather conditions, especially in typhoons and hurricanes.Rain not only increases the atmospheric attenuation but also changes the sea surface roughness in a complicated manner.Thus, it is difficult to accurately model brightness temperatures under rainy conditions.Hence, the statistical algorithm is a relatively good candidate that can be used to retrieve wind speed under severe weather conditions.In general, the new wind speed retrieval algorithm under rain gives an encouraging accuracy compared to the H*wind wind speed.

Comparison with Remote Sensing Systems (RSS) All-Weather Data
A global all-weather wind speed product has been developed using the WindSat TBs [10], which is a good candidate for comparison to our retrieved wind speed results based on WindSat TBs inside hurricanes.Table 3 illustrates the retrieved wind speed using our new algorithm and the RSS global algorithm, respectively.The overall RMS difference of wind speed between our new model and H*wind analysis data is lower than that between the RSS wind speed and H*wind analysis data.In addition, results are similar for each individual case expect for hurricane Bill.All biases are negative between the RSS wind speed and H*wind analysis data.That means the RSS wind speed underestimates wind speed on the whole, especially for high wind speeds.The mean bias is −2.53 m/s between them.Compared to the H*wind analysis data, the retrieved wind speed using the new method overestimates wind speed slightly, and the mean bias is 0.04 m/s for the selected 15 hurricanes.The SFMR measurements are compared with our retrieved wind speed results, although these data can only provide point observations along the aircraft flight track.The 30-min window is utilized to collocate the matchups between the WindSat data and SFMR measurements.Since the SFMR measurements are at a high spatial resolution of 1.5 km, a method similar to that described in section 2.5 is also performed.Finally, we obtained 938 wind speed matchups for the selected 15 hurricanes.Figure 9 suggests that the derived wind speed has a relatively good agreement with the SFMR measurements, with a mean bias of −0.68 m/s and an RMS difference of 3.14 m/s.The general trend is that the retrieved wind speed errors increase with increasing SFMR rain rate.

Case Studies
As a case study (denoted as Rita-1), two different models are used to retrieve wind speed using the WindSat TBs for hurricane Rita on 22 September 2005, at 1330 UTC. Figure 10a presents the H*wind analysis wind speed, and Figure 10b shows the retrieved wind speed using our new model.Figure 10c shows the retrieved wind speed using the RSS global algorithm, and Figure 10d illustrates the corresponding RSS rain rate.Figure 10e shows the difference between the retrieved wind speed and the H*wind analysis data.The difference between the RSS wind speed and the H*wind analysis data is presented in Figure 10f.In addition, the corresponding W6H and W6V are presented in Figure 10g,h, respectively.In this case, the H*wind analysis wind speed is obtained for 22 September 2005, at 1330 UTC.At that time, the center of typhoon's eye was located at 25.197°N and −88.531°E.However, the center of typhoon's eye was located at 25.141°N and −88.237°E at the time of the WindSat observation, which is obtained from the NHC.Therefore, we shifted the H*wind analysis wind speed so that the eye of the H*wind analysis coincides with eye of the WindSat measurement.Comparing W6H with W6V, we find that W6H seems to have a better correlation than W6V.The distribution of W6H is similar to that of wind speed.The performance of the new model is satisfied, and the distribution of retrieved wind speed is very close to the H*wind analysis wind speed, with the mean bias of 0.79 m/s and RMS difference of 1.79 m/s.As for RSS wind speeds, the influence of rain is obvious.On the contrary, our new model can reduce the influence of rain to some extent, and the retrieved results (Figure 10e) are in good agreement with the H*wind analysis wind speed.In general, the RSS global wind speed underestimates the wind speed in most regions, and the retrieved results are vulnerable to rain (Figure 10f).In addition, it is seen that the retrieved wind speed overestimates wind speed in the center of the eye for our retrieved results and RSS wind speeds.This may be because there are still strong waves that increase the emission from ocean surface, even though the wind speed is relatively low in the eye of the typhoon.
In a second case study (denoted as Rita-2), the two different models are also utilized to retrieve wind speed with the WindSat TBs for hurricane Rita on 22 September 2005, at 0130 UTC. Figure 11b illustrates the retrieved wind speed using our new model.Figure 11c shows the retrieved wind speed using the RSS global algorithm, and Figure 11d shows the corresponding RSS rain rate.Figure 11e-f present the difference between the retrieved wind speeds and the H*wind analysis data.In addition, Figures 11g,h show the corresponding W6H and W6V results, respectively.In this case, the center of the typhoon's eye is located at 24.536°N and −87.238°E for the H*wind analysis data; however, the center of typhoon's eye is located at 24.535°N and −86.808°E at the time of the WindSat observation.Hence, we shifted the H*wind analysis wind speed so that the eye of the H*wind analysis wind speed coincides with eye of the WindSat measurement.The results show that the performance of the new model is satisfied, and the mean bias and RMS difference of the new model are 0.63 m/s and 2.38 m/s, respectively.Near the south of the typhoon's eye, our new model seems to underestimate wind speed when the rain rate is relatively low, and similar results are also found in the RSS global algorithm.

Case Studies
As a case study (denoted as Rita-1), two different models are used to retrieve wind speed using the WindSat TBs for hurricane Rita on 22 September 2005, at 1330 UTC. Figure 10a presents the H*wind analysis wind speed, and Figure 10b shows the retrieved wind speed using our new model.Figure 10c shows the retrieved wind speed using the RSS global algorithm, and Figure 10d illustrates the corresponding RSS rain rate.Figure 10e shows the difference between the retrieved wind speed and the H*wind analysis data.The difference between the RSS wind speed and the H*wind analysis data is presented in Figure 10f.In addition, the corresponding W6H and W6V are presented in Figure 10g,h, respectively.In this case, the H*wind analysis wind speed is obtained for 22 September 2005, at 1330 UTC.At that time, the center of typhoon's eye was located at 25.197 • N and −88.531 • E. However, the center of typhoon's eye was located at 25.141 • N and −88.237 • E at the time of the WindSat observation, which is obtained from the NHC.Therefore, we shifted the H*wind analysis wind speed so that the eye of the H*wind analysis coincides with eye of the WindSat measurement.Comparing W6H with W6V, we find that W6H seems to have a better correlation than W6V.The distribution of W6H is similar to that of wind speed.The performance of the new model is satisfied, and the distribution of retrieved wind speed is very close to the H*wind analysis wind speed, with the mean bias of 0.79 m/s and RMS difference of 1.79 m/s.As for RSS wind speeds, the influence of rain is obvious.On the contrary, our new model can reduce the influence of rain to some extent, and the retrieved results (Figure 10e) are in good agreement with the H*wind analysis wind speed.In general, the RSS global wind speed underestimates the wind speed in most regions, and the retrieved results are vulnerable to rain (Figure 10f).In addition, it is seen that the retrieved wind speed overestimates wind speed in the center of the eye for our retrieved results and RSS wind speeds.This may be because there are still strong waves that increase the emission from ocean surface, even though the wind speed is relatively low in the eye of the typhoon.
In a second case study (denoted as Rita-2), the two different models are also utilized to retrieve wind speed with the WindSat TBs for hurricane Rita on 22 September 2005, at 0130 UTC. Figure 11b illustrates the retrieved wind speed using our new model.Figure 11c shows the retrieved wind speed using the RSS global algorithm, and Figure 11d shows the corresponding RSS rain rate.Figure 11e-f present the difference between the retrieved wind speeds and the H*wind analysis data.In addition, Figure 11g,h show the corresponding W6H and W6V results, respectively.In this case, the center of the typhoon's eye is located at 24.536 • N and −87.238 • E for the H*wind analysis data; however, the center of typhoon's eye is located at 24.535 • N and −86.808 • E at the time of the WindSat observation.Hence, we shifted the H*wind analysis wind speed so that the eye of the H*wind analysis wind speed coincides with eye of the WindSat measurement.The results show that the performance of the new model is satisfied, and the mean bias and RMS difference of the new model are 0.63 m/s and 2.38 m/s, respectively.Near the south of the typhoon's eye, our new model seems to underestimate wind speed when the rain rate is relatively low, and similar results are also found in the RSS global algorithm.

Rain Effects on Wind Retrieval
The passive radiometer measurement is less sensitive to the ocean surface wind speed under rainy conditions.Rain increases the atmospheric attenuation, which is obvious at higher frequencies (Figure 3).In addition, it is very difficult to accurately simulate the brightness temperatures in rain, especially for the high variability of rainy atmospheres.The brightness temperatures depend on cloud type and the distribution of rain within the footprint [32,51].Then, the full Mie absorption theory needs to be applied to calculate the atmospheric radiative transfer equation, which requires additional input parameters, such as size and form of the rain drops; however, these parameters are not readily available.Therefore, previous studies have mainly focused on the statistical algorithm for retrieving wind speed under rainy conditions [10,15,17].

Rain Effects on Wind Retrieval
The passive radiometer measurement is less sensitive to the ocean surface wind speed under rainy conditions.Rain increases the atmospheric attenuation, which is obvious at higher frequencies (Figure 3).In addition, it is very difficult to accurately simulate the brightness temperatures in rain, especially for the high variability of rainy atmospheres.The brightness temperatures depend on cloud type and the distribution of rain within the footprint [32,51].Then, the full Mie absorption theory needs to be applied to calculate the atmospheric radiative transfer equation, which requires additional input parameters, such as size and form of the rain drops; however, these parameters are not readily available.Therefore, previous studies have mainly focused on the statistical algorithm for retrieving wind speed under rainy conditions [10,15,17].

Rain Effects on Wind Retrieval
The passive radiometer measurement is less sensitive to the ocean surface wind speed under rainy conditions.Rain increases the atmospheric attenuation, which is obvious at higher frequencies (Figure 3).In addition, it is very difficult to accurately simulate the brightness temperatures in rain, especially for the high variability of rainy atmospheres.The brightness temperatures depend on cloud type and the distribution of rain within the footprint [32,51].Then, the full Mie absorption theory needs to be applied to calculate the atmospheric radiative transfer equation, which requires additional input parameters, such as size and form of the rain drops; however, these parameters are not readily available.Therefore, previous studies have mainly focused on the statistical algorithm for retrieving wind speed under rainy conditions [10,15,17].
Using the RSS rain rate products, the rain effect is evaluated in the new model.Table 4 presents the biases and RMS errors of the difference between the wind speeds from the retrievals and the corresponding H*wind analysis wind speeds as a function of the WindSat rain rate.In general, the biases are negative when the rain rate is less than 4 mm/h or more than 14 mm/h.Therefore, the retrieved wind speed underestimates wind speed in this situation.The biases become positive when the rain rate ranges from 4 to 14 mm/h.In general, the wind speed retrieval errors increase with increasing rain rate.The wind speed retrieval accuracy of our new model inside hurricane ranges from 2.0 m/s in light rain to 3.9 m/s in heavy rain.

Conclusions
In this study, we assess the wind-induced emissivity changes in C-band and X-band channels using simulated brightness temperatures.Four parameters (6V*, 6H*, 10V*, and 10H*) are defined, and they present the quantity of ocean microwave emission changed by ocean wind.The results show that the sensitivity of 6(10)V* to wind speed is approximately 0.37(0.5)K/(m/s) and that the sensitivity of 6(10)H* to wind speed is approximately 0.85(0.93)K/(m/s).Furthermore, we find that the ratio of 6V(H)* to 10V(H)* is a constant, which is approximately 0.82(0.92).These results are in good agreement with previous studies [15,28].Based on the work in [15], we define two new parameters, W6H and W6V, which represent increments in 6H and 6V due to ocean wind, respectively.Finally, a wind speed retrieval algorithm for wind speeds above 20 m/s in hurricanes has been developed using the two new parameters.The C-band and X-band vertically and horizontally polarized brightness temperatures are used to retrieve high wind speeds because the brightness temperature acquired at higher frequencies become saturated inside hurricanes.
The WindSat measurements over 15 hurricanes from 2003 to 2010 and the corresponding H*wind analysis data are used to develop and validate the retrieval model.The performance is encouraging, with a mean bias of 0.04 m/s and an overall RMS difference of 2.75 m/s.In general, the retrieved wind speeds have a good agreement with the H*wind analysis wind speeds.Two case studies are performed, resulting in a mean bias and RMS difference of 0.79 m/s and 1.79 m/s for hurricane Rita-1 and 0.63 m/s and 2.38 m/s for hurricane Rita-2.Compared to the RSS global all-weather wind speeds, our retrieved wind speeds is closer to the H*wind analysis data, and the RSS global all-weather wind speeds underestimates the wind speed.In addition, the retrieved wind speeds are compared to the SFMR wind speeds.Good agreement is found between them.
It is a challenging job to retrieve wind speed under rainy conditions, especially for typhoons and hurricanes.The core of the new algorithm for determining the wind speed under rainy conditions involves using different channel combinations that reduce the signal coming from rain without greatly reducing the wind speed signal.With respect to the effect of rain, the performance of the new model inside hurricanes ranges from 2.0 m/s in light rain to 3.9 m/s in heavy rain.The performance of retrieving wind speed under rainy conditions is satisfactory.
This new algorithm can be used to retrieve wind speeds in hurricanes under rainy conditions with an encouraging degree of accuracy, although we cannot expect the algorithm to have the same accuracy as in rain-free cases.It is noted that the retrieved wind speed always overestimates wind speed in the eye of the hurricane.That may be because there are still strong waves that increase the emission from the ocean surface, even if the wind speed is relatively low in the eye of hurricane.For practical applications of the new wind speed algorithm, we should update the current method for hurricanes in which the areas of high wind and high rain are not fully observed.In addition, the relationship between wind speed and wind-induced emissivity under different sea states, wave ages, and air-sea stability in hurricanes should be studied.Furthermore, the wind direction retrieval will be studied using the polarimetric brightness temperatures in the future.

Figure 2 .
Figure 2. The brightness temperatures received by a radiometer.

Figure 2 .
Figure 2. The brightness temperatures received by a radiometer.

Figure 3 .
Figure 3. Brightness temperatures measured by WindSat orbit #3510 on 10 September 2003, over hurricane Isabel (a) at 10.7 GHz, horizontal polarization; (b) at 18.7 GHz, horizontal polarization; (c) at 37.0 GHz, horizontal polarization; (d) at 10.7 GHz, vertical polarization; (e) at 18.7 GHz, vertical polarization and (f) at 37.0 GHz, vertical polarization.The red curve denotes the best track, which was obtained from the NHC.

Figure 3 .
Figure 3. Brightness temperatures measured by WindSat orbit #3510 on 10 September 2003, over hurricane Isabel (a) at 10.7 GHz, horizontal polarization; (b) at 18.7 GHz, horizontal polarization; (c) at 37.0 GHz, horizontal polarization; (d) at 10.7 GHz, vertical polarization; (e) at 18.7 GHz, vertical polarization and (f) at 37.0 GHz, vertical polarization.The red curve denotes the best track, which was obtained from the NHC.

Figure 5 .
Figure 5. (a) Brightness temperatures of 6H and 10H: group A for calm ocean conditions and group B for rough ocean conditions; (b) Thematic figure drawn from (a).The two figures can be found in [15].

Figure 5 .
Figure 5. (a) Brightness temperatures of 6H and 10H: group A for calm ocean conditions and group B for rough ocean conditions; (b) Thematic figure drawn from (a).The two figures can be found in [15].

Figure 7 .
Figure 7. Retrieved wind speed of three hurricanes using our new method and the Shibata method with the new coefficients versus H*wind analysis data.(a) Our new method and (b) the method of Shibata (2006) with the new coefficients.

4. 2 . 1 .
Figure 8 presents the retrieved wind speed using the new model versus the corresponding H*wind analysis data for each selected hurricane.The mean bias and RMS error are given in the top of each single panel.For wind speeds above 20 m/s, the mean bias and overall RMS difference of the new model are 0.04 m/s and 2.75 m/s, respectively, for 15 hurricanes from 2003 to 2010.Positive biases are found in some hurricanes, including Dennis, Katrina, Rita, Wilma, Helene, Ike, Bill, and Ida.For

Figure 7 .
Figure 7. Retrieved wind speed of three hurricanes using our new method and the Shibata method with the new coefficients versus H*wind analysis data.(a) Our new method and (b) the method of Shibata (2006) with the new coefficients.

4. 2 . 1 .
Figure 8 presents the retrieved wind speed using the new model versus the corresponding H*wind analysis data for each selected hurricane.The mean bias and RMS error are given in the top of each single panel.For wind speeds above 20 m/s, the mean bias and overall RMS difference of the new model are 0.04 m/s and 2.75 m/s, respectively, for 15 hurricanes from 2003 to 2010.Positive biases are found in some hurricanes, including Dennis, Katrina, Rita, Wilma, Helene, Ike, Bill, and Ida.For

Figure 7 .
Figure 7. Retrieved wind speed of three hurricanes using our new method and the Shibata method with the new coefficients versus H*wind analysis data.(a) Our new method and (b) the method of Shibata (2006) with the new coefficients.

4. 2 . 1 .
Figure 8 presents the retrieved wind speed using the new model versus the corresponding H*wind analysis data for each selected hurricane.The mean bias and RMS error are given in the top of each single panel.For wind speeds above 20 m/s, the mean bias and overall RMS difference of the new model are 0.04 m/s and 2.75 m/s, respectively, for 15 hurricanes from 2003 to 2010.Positive biases are found in some hurricanes, including Dennis, Katrina, Rita, Wilma, Helene, Ike, Bill, and Ida.

Figure 9 .
Figure 9. (a) The comparison between SFMR wind speeds and the retrieved wind speeds with our new method; (b) Errors of the retrieved wind speeds versus SFMR rain rate.

Figure 9 .
Figure 9. (a) The comparison between SFMR wind speeds and the retrieved wind speeds with our new method; (b) Errors of the retrieved wind speeds versus SFMR rain rate.

Figure 10 .
Figure 10.(a) H*wind wind speed of Rita on 22 September 2005, at 1330 UTC; (b) Retrieved wind speed of Rita with our new model on22 September 2005, at 1155 UTC; (c) RSS all-weather wind speed of Rita on 22 September 2005, at 1155 UTC; (d) RSS rain rate of Rita on 22 September 2005, at 1155 UTC; (e) The bias between H*wind and Retrieved wind speed with our new model; (f) The bias between H*wind and RSS all-weather wind speed; (g) W6H over Rita on 22 September 2005, at 1155 UTC; (h) W6V over Rita on 22 September 2005, at 1155 UTC.

Figure 11 .
Figure 11.(a) H*wind wind speed for Rita on 22 September 2005, at 0130 UTC; (b) Retrieved wind speed based on our new model for Rita on 21 September 2005, at 2328 UTC; (c) RSS all-weather wind speed for Rita on 21 September 2005, at 2328 UTC; (d) RSS rain rate for Rita on 21 September 2005, at 2328 UTC; (e) The bias between H*wind and the wind speed retrieved with our new model; (f) The bias between H*wind and the RSS all-weather wind speed; (g) W6H over Rita on 21 September 2005, at 2328 UTC; (h) W6V over Rita on 21 September 2005, at 2328 UTC.

Figure 10 . 20 Figure 10 .
Figure 10.(a) H*wind wind speed of Rita on 22 September 2005, at 1330 UTC; (b) Retrieved wind speed of Rita with our new model on22 September 2005, at 1155 UTC; (c) RSS all-weather wind speed of Rita on 22 September 2005, at 1155 UTC; (d) RSS rain rate of Rita on 22 September 2005, at 1155 UTC; (e) The bias between H*wind and Retrieved wind speed with our new model; (f) The bias between H*wind and RSS all-weather wind speed; (g) W6H over Rita on 22 September 2005, at 1155 UTC; (h) W6V over Rita on 22 September 2005, at 1155 UTC.

Figure 11 .
Figure 11.(a) H*wind wind speed for Rita on 22 September 2005, at 0130 UTC; (b) Retrieved wind speed based on our new model for Rita on 21 September 2005, at 2328 UTC; (c) RSS all-weather wind speed for Rita on 21 September 2005, at 2328 UTC; (d) RSS rain rate for Rita on 21 September 2005, at 2328 UTC; (e) The bias between H*wind and the wind speed retrieved with our new model; (f) The bias between H*wind and the RSS all-weather wind speed; (g) W6H over Rita on 21 September 2005, at 2328 UTC; (h) W6V over Rita on 21 September 2005, at 2328 UTC.

Figure 11 .
Figure 11.(a) H*wind wind speed for Rita on 22 September 2005, at 0130 UTC; (b) Retrieved wind speed based on our new model for Rita on 21 September 2005, at 2328 UTC; (c) RSS all-weather wind speed for Rita on 21 September 2005, at 2328 UTC; (d) RSS rain rate for Rita on 21 September 2005, at 2328 UTC; (e) The bias between H*wind and the wind speed retrieved with our new model; (f) The bias between H*wind and the RSS all-weather wind speed; (g) W6H over Rita on 21 September 2005, at 2328 UTC; (h) W6V over Rita on 21 September 2005, at 2328 UTC.
The hurricanes selected for this paper are Fabian and Isabel in 2003; Frances and Ivan in 2004; Dennis, Katrina, Rita, and Wilma in 2005; Helene in 2006; Felix in 2007; Bertha and Ike in 2008; Bill and Ida in 2009; and Igor in 2010.The hurricane best track analysis and the corresponding intensity are depicted in Figure

Table 1 .
Statistics of matchups between WindSat observations and H*wind analysis data for each selected TC.

Table 3 .
Statistics of errors between the new retrieved wind speed and H*wind wind speed and between RSS all-weather wind speed and H*wind wind speed for H*wind wind speeds above 20 m/s for the selected hurricanes from 2003 to 2010.

Table 4 .
The statistics for different rain intervals between H*wind wind speed and the wind speed retrieved with our new method.