Scatterometer hurricane wind speed retrievals using cross polarization

Introduction Conclusions References


Introduction
Wind speeds and/or directions over oceans are routinely measured with high spatial and temporal resolution by active radar instruments such as scatterometers and Synthetic Aperture Radars (SAR).The winds can be inferred from the measured Normalized Radar Cross Section (NRCS) signal strength (sigma0 or σ 0 ) as these are linked to the ocean roughness and its azimuthal variation, even though the exact mechanisms behind this are still not fully understood theoretically (e.g.Romeiser et al., 1997).
In the case of scatterometers, the wind speed and direction can be retrieved using three collocated measurements from different azimuth viewing angles combined with the geophysical C-band model (CMOD) function that empirically links the radar backscatter to the local wind vector and incidence angle (e.g.Stoffelen and Portabella, 2006).In the case of SAR instruments, additional wind direction (or speed) information is required to determine the absolute wind speed using co-polarization (vertically transmitted and vertically polarized, VV or its horizontal equivalent, HH) signals.
Near real-time global wind estimates from various scatterometer instruments provide spatial resolutions on the order of tens of kilometers and are proven to be very accurate (Vogelzang et al., 2011).Scatterometers are used in diverse applications including marine meteorology, Numerical Weather Prediction (NWP), oceanography and climate (Bourassa et al., 2009;Vogelzang and Stoffelen, 2012).The high spatial resolution (tens of meters) provided by SAR instruments, can be used to determine local wind speed phenomena in bays, fjords and along the coast line provided the prevailing wind direction is uniform (Portabella et al., 2002).Besides their limited coverage, another main limitation of SAR systems is their calibration stability.Scatterometers obtain a precision of < 0. and a rather simple and useful relationship to the severe wind speed has emerged.
The operational use of cross-polarization signals for wind retrievals over the ocean will only be possible in future meteorological satellite systems.The Meteorological Operational satellite programme Second Generation (MetOp-SG) satellites will replace the current MetOp system in the 2020+ time frame.One instrument to be carried on-board the MetOP-SG is a C-band scatterometer (SCA), similar to ASCAT on MetOp, but with higher spatial resolution, increased coverage and stability.The prime objective of SCA, herein referred to as ASCAT-SG, will be to measure wind speed and direction over the oceans, following its predecessor ASCAT on MetOp.Another innovation in the design of ASCAT-SG is the inclusion of a single cross-polarization beam, which can be used for the retrieval of severe wind speeds and therefore be a prime improvement in the nowcasting of hurricane intensities, for the protection of coastal residents and infrastructure.
In preparation for ASCAT-SG, the cross polarization (VH) data from the RADARSAT-2 satellite are compared to wind speeds (Sect.2) measured by NOAA's Hurricane Hunter flights (Sect.2.2) and to wind speeds provided by ECMWF forecasts (Sect.2.4).The results indicate that VH backscatter is capable of retrieving strong to extreme/severe wind speeds without any dependence on wind direction.As such, VH backscatter could be used to complement the standard VV products by extending wind speed retrieval beyond 25 m s −1 .In Sect.3, a new VH Geophysical Model Function is described and used to retrieve the wind speeds in hurricane Earl (2010).Sect. 4 describes the use of RADARSAT-2 data to nowcast the Hurricane best track estimate, as defined by NOAA, thereby showing how cross polarization measurements can be exploited to warn people about severe wind events.Finally Sect. 5 provides the main conclusions of this work.Figures

Back Close
Full 2 Cross-polarization wind speed relationships 2.1 Low-to-strong wind speed measurements using buoys and RADARSAT-2 In the last few years, a number of papers (e.g.Zhang et al., 2011;Zhang and Perrie, 2011;Hwang et al., 2010;Vachon and Wolfe, 2011) on using the cross-polarized (VH or HV) synthetic aperture radar (SAR) imagery for the retrieval of wind speeds have been published.These papers are based on data measured by the RADARSAT-2 satellite, which was launched in December 2007 for the Canadian Space Agency.RADARSAT-2 deploys a C-band SAR which can be configured in multiple polarization modes, including a fully polarimetric mode in which HH, HV, VV and VH polarized data and phase information can be acquired (Slade, 2011).This fine-quad polarization mode has an extremely low noise floor and cross-talk correction between the different channels.The RADARSAT-2 mode with the widest swath (ScanSAR Wide A [SCWA]) has a nominal width of 500 km, an imaging resolution of 100 m, operates in dual-polarization receive mode (VV and VH, or HH and HV), and shows a significantly higher noise floor (> −30 dB) when compared to the fine-quad polarization mode.
In the low-to-strong wind speed regime (≤ 20 m s −1 ) the cross-polarized backscatter by RADARSAT-2 fine quad-polarization (HH, HV, VH, and VV) mode (see above mentioned papers) has been compared to wind vectors measured at buoys.It was found that the cross-polarized backscatter in dB was insensitive to incidence angle and wind direction resulting in a simple linear geophysical relationships that may be directly used to derive absolute wind speeds.In contrast to this, the VV backscatter shows a strong incidence angle and azimuthal wind direction dependence.This VV dependence complicates the retrieval of wind speeds considerably but does enable the retrieval of wind direction (Wackerman et al., 1996;Portabella et al., 2002;Horstmann and Koch, 2005;Stoffelen and Portabella, 2006), which is not possible for the VH channel.The absolute wind speed relationships in the low-to-strong wind speed regime based on cross polarization data resulted in very similar results by both Vachon and Wolfe (2011)  (2) Their results showed the potential of a robust wind speed retrieval method based on the cross-polarization signals alone for Hurricane wind speeds.
The results in both Zhang et al. (2011) and Hwang et al. (2010), where RADARSAT-2 (SCWA) images of hurricanes were used in combination with VV-wind retrievals and NWP model winds, showed that the VH [dB]-wind relationship for strong-to-severe wind speeds follow a different slope compared to that found in Eqs. ( 1) and ( 2).These evaluations were performed however with a very limited hurricane dataset.The derivation of a statistical VH relationship for strong-to-severe (> 20 m s −1 ) wind speeds based on a large number of Hurricane data is one of the main objectives of this work.
In this work, all presented wind speed fits are related to VH measurements in dB, resulting in the wind speed relationships being described as linear with respect to VH [dB].Note that this means that the VH measurements themselves increase exponentially with respect to wind speed.As the SFMR measures the wind speed below the aircraft, the data provides transects/tracks of the wind speeds through the individual hurricanes and not the full wind speed pattern throughout the hurricane system.The NOAA Hurricane Hunters flew a total of 18 collocated flight tracks through 9 of the 19 available hurricanes with the SFMR instrument.In Table 1, the hurricane images used in this study are indicated with the dates and times of RADARSAT-2 overpasses including the number of collocated flight tracks performed within the RADARSAT-2 images which were close enough in space and time to be representative (see below). Figure 1 shows a ScanSAR Wide-A (SCWA) RADARSAT-2 image (converted to σ 0 in dB) of hurricane Earl as observed on 2 September 2010 with the SAR configured in dual-pol mode with VV (bottom panel)

Strong
and VH (top panel).The VV σ 0 is observed (and expected) to decrease with increasing incidence angle (i.e. from left-to-right in the image); whereas, the VH σ 0 does not appear to have any dependence with increasing incidence angle.This is also apparent in the well defined eye of the hurricane in the VH image in comparison to the one in the VV image.At low incidence angles the VV-signal saturates close to 0 dB (left side of the image).The VH signal strength is lower compared to the co-polarized VV signals.
The full image is created by combining 4 individual beams (W1, W2, W3 and S7) to generate SCWA image mode.The four beams are essentially stitched together with a blend operation performed at the beam seams to attempt to mitigate undesirable beam seams.Beam seams can be evident in low signal-to-noise ratio regions as is observed in the VH image as linear features that run over the full along-track extent on either side of the hurricane eye.The beam seams are not as evident in the VV image due to higher signal-to-noise values.The data provided by the RADARSAT-2 team is a summation of the measured geophysical signal (σ 0 ) and the Noise Equivalent Sigma Zero (NESZ).The NESZ is the noise (bias) using the known and modeled properties of the instrument and has a value around −30 dB, which may vary by 1-2 dB.The local value of the NESZ, which is incidence-angle dependent, is assumed to be correct within 1 dB.larger than NESZ + 1 dB are taken into account within this work, before the NESZ is subtracted from the measured signals to obtain the noise-corrected σ 0 .

Collocating RADARSAT-2 with NOAA Hurricane-Hunter flights
The construction of wind speed versus cross-polarization relationships require a validated wind product.For a number of hurricanes there is the availability of hurricane surface winds from the NOAA WP-3D research aircraft (Uhlhorn et al., 2007).A typical NOAA flight takes around 6 h and within this time a number of flight legs are flown through the eye of the hurricane.Compared to the near instantaneous measurements of RADARSAT-2, both data streams will never be fully collocated in time and space, especially when taking into account the Hurricanes intrinsic movement during the flight.
Therefore, each of the individual legs (lasting about 20 min) has to be collocated separately.The number of individual flight legs for each hurricane, in which a NOAA flight was close enough in space and time is presented in Table 1.To determine if a flight leg is a candidate for inclusion in the dataset, the time difference between the RADARSAT-2 and NOAA flight acquisitions is first calculated.The median results for the RADARSAT-2 collocated ECMWF-fields 1 , the HWIND model 2 and the retrieved wind from the VV channel using the CMOD5.ngeophysical model function (Verhoef, 2008) are overplotted on the distribution.The ECMWF wind fields are obtained from short-range forecasts with three-hour intermittency.HWIND is a model field which provides an integrated view of the extent and strength of the wind field and is used operationally to improve the assessment of hurricane intensity (Powell et al., 1998).The computation of the VV-GMF wind speed requires knowledge of the VV(σ 0 ), incidence angle and wind direction.Since RADARSAT-2 detects only one azimuth viewing direction the wind retrievals can not be determined independently.Within this study, two different wind direction descriptions, required for the VV-GMF calculations are used.In the first the forecasted ECMWF wind direction of the 10 m wind is assumed, secondly the assumption that the wind rotates perfectly circular around Introduction

Conclusions References
Tables Figures

Back Close
Full the center of a hurricane.Comparisons showed that there is no significant difference between the U 10 (VV, φ ECMWF ) and U 10 (VV, φ rotate ) versus VH [dB] distributions.A third method to determine wind direction, by using surface features (Lehner et al., 1998;Horstmann and Koch, 2005), has not been adopted in this work.
The ECMWF and HWIND models, VV-GMF and the median SFMR measurements show a similar change in slope for severe wind speeds (> 20 m s −1 ).The SFMR measurements are between 1 and 3 dB lower than the VV-GMF as the wind speed increases from 20 to 40 m s −1 respectively.Furthermore, below 20 m s −1 , the slope for both the VV-GMF and the median SFMR measurements are similar to the Vachon and Wolfe relationship.The VV-GMF median fit depicts the highest slope, which indicates its lack of sensitivity for extreme wind retrievals.The collocated ECMWF forecast data has a similar slope up to 20 m s −1 but shows a lower sensitivity in the extreme wind regime.This is most likely due to the lower spatial resolution in the model, thereby smoothing the largest gradients close to the outer eye wall (no ECMWF winds within the eye wall were used).The HWIND model shows a larger wind speed at all collocated be a shift in parameter space between small incidence angle (top-left of the distribution) and large incidence angle (bottom right); however, not enough independent data is available from the SFMR data to determine any underlying incidence angle dependence directly.In the next section, the incidence angle dependence is analyzed in more detail by comparing the VH data to the ECMWF wind velocity forecasts.

Incidence and azimuthal angle dependence using ECMWF forecasts
An independent method to assess the incidence angle and wind direction dependencies is performed by comparing the VH measurements to the NWP model calculations or wind retrievals from the VV-band.By incorporating all individual VH pixels for the 19 hurricane images enough data is acquired to retrieve these dependencies statistically.
Both comparisons have their pros and cons.The VV wind speed retrievals are perfectly collocated to the VH observations but have sensitivity issues when retrieving extreme wind speeds beyond 25 m s −1 .On the other hand, one could use forecasts from an NWP model at all speeds.However, a global NWP model, in this case the ECMWF model, has a lower resolution and is not collocated in space and time.The lower resolution will result in a smoothed wind field, potentially missing small regions with high wind speeds or high wind shear.The ECMWF model is known to provide good forecasts in the case of tropical hurricanes; therefore, it is expected that the model will provide a good statistical description of the general wind field except close to the hurricane eye where the spatial variability is largest.
In this work, the ECMWF model comparison was preferred to the VV retrieved wind speeds because the VV retrieved wind speed depends on both the incidence and azimuth angles.If these dependencies are not perfectly defined in the geophysical model function (VV-GMF), then any offset will feed into a newly developed VH-GMF parameterization.In contrast, the ECMWF model is independent of the RADARSAT-2 geometry and therefore, distributions of collocated data should not produce a bias in any potential azimuthal or incidence angle dependence.The ECMWF wind speeds may however, under-estimate the maximum sustained winds as measured by RADARSAT-2.

Back Close
Full For all the Hurricanes, the ECMWF forecast was retrieved on a 0.25 × 0.25 degrees grid.Since the ECMWF provides three hour outputs, the closest time slot to the RADARSAT-2 image overpass was used.This results in a small offset in collocation which has to be corrected for.The correction offset method was preferable to interpolating the ECMWF fields to the correct time which induced stretching of the Hurricane shape.For each of the ECMWF fields, the center of the hurricane was determined by looking for the hurricane eye in the modeled velocity, after which, the ECMWF field was translated to match the VH hurricane center.The ECMWF field was subsequently bi-linearly interpolated (up-scaled) to the RADARSAT-2 resolution.Since the ECMWF resolution is lower, compared to the RADARSAT-2 data, it will miss the high variations in wind speed close to the hurricane eye.To reduce the amount of potential missing features, the radius of maximum wind speed (R ecmwf ) in the ECMWF wind field is first calculated.All wind-velocities within 1.1 times R ecmwf are removed from the comparison to reduce the regimes for which the ECMWF model will clearly be biased low.This precaution filters out the highest wind velocities from the ECMWF model fields.
In Fig. 4, the number density histogram of the combined hurricane distribution is shown.The individual hurricanes overlap in the same part of parameter space resulting in a tight distribution with a correlation of 0.85 in U 10 vs. VH [dB].In the bottom panel, the median incidence angle, binned in wind speed and VH strength, is plotted showing a very similar incidence dependence as was suggested by the SFMR data (Fig. 3).This could indicate a VH GMF incidence-angle dependency and/or sampling or calibration problem in the RADARSAT-2 data, either, of which would reduce the correlation in the plots.

Incidence angle dependence
In order to resolve the potential incidence angle dependence, a number of linear regression lines (Eq. were performed in between 7 and 20 m s −1 and not the entire wind speed regime.There is a clear regime change at around 20 m s −1 from the low-to-strong to strong-to-severe wind speeds, a single linear relationship would not suffice for the entire distribution.It is expected that, based on the lower panel of Fig. 4, a similar type of incidence angle dependence is required for nominal and severe wind speeds.There are insufficient points in the U 10 > 20 m s −1 part of the distribution available (especially in the low incidence angle part) to retrieve an incidence-angle relationship.Since there is a change in slope, and therefore the VH dependence on the ocean surface, one cannot exclude that the incidence-angle dependence is different for the two regimes.Additional high wind hurricane images are required to probe the strong-to-severe wind speed regime more thoroughly in a future effort.
The resulting A and B coefficients in the 7 < U 10 < 20 m s −1 regime are plotted and a second order polynomial, e.g.A(θ) = a 0 + a 1 θ + a 2 θ 2 , was fitted through the points as depicted in Fig. 5 (top left and right panels respectively).The fitted parameters are provided in Table 2. Based on this parameterization, the measured VH signals can be globally made incidence angle independent by shifting all measurements to a single incidence angle, in this case 35 degrees, using the following equations: where ∆dB is calculated using Eq. ( 4), the fitted second order polynomial and θ 0 = 35.Note that by describing the correction in this manner, it does not depend on the a 0 and b 0 values but only on the θ dependent coefficients.This reduces any bias between the Introduction

Conclusions References
Tables Figures

Back Close
Full ECMWF forecast wind-speeds and the VH wind speeds for each incidence angle.As such, the resulting correlation between the ECMWF velocity and VH signal improves from 0.85 to 0.89, which, by itself, is not significant.Also, the distance between the 0.2 and 0.8 VH percentiles decrease to a minimum of 60 % at 10 m s −1 up to 80 % at 35 m s −1 compared to the original VH-U 10 (ECMWF) distribution.The decrease in distribution width is readily visible in the bottom panel of Fig. 5.
Both the SFMR and ECMWF comparisons show an incidence-angle dependence in contrast to what was found using the buoy data for wind speeds below 20 m s −1 (Vachon and Wolfe, 2011).Correcting the data for this dependence improves the description of the VH incidence-angle distribution and tightens the U 10 (ECMWF)-VH distribution.The joint wind speed and VH distribution is obtained through a combination of the matched sampled position and strengths of the wind speed.In this study, the wind speeds are not random with respect to their incidence angle, since hurricanes were targeted, which is the case for the comparison against buoys.As the data consists of a limited number (19) of hurricane images, it is conceivable that the coincident hurricane positions and strengths (distribution) (Fig. 6) within the RADARSAT-2 swath, as conceived by the different measurement systems introduces an apparent incidence angle dependency.
There is a relatively large number of hurricanes in the high incidence angle regime (35-40 • ), where also the incidence angle correction peaks (Fig. 5, upper panels).On the other hand, the incidence angle dependence in the U 10 -VH distribution (Fig. 4) lower panel shows a very similar behavior throughout the velocity range.In an additional test, performed to check for a potentially induced incidence angle dependence, the VH-U 10 (ECMWF) distribution was tested using a bootstrap technique.In this test, ten thousand sample distributions were created between 7 and 20 m s −1 , each consisting of 19 randomly drawn hurricane images using simple random sampling with replacement, i.e. an individual sample can consist of a single hurricane more than once.The incidence angle relationship (a i and b i ; Eq. 4) coefficients where determined for each of the samples and the moments for each of the coefficients were calculated.For all coefficients, the median and mean of the distribution where within 7 % of the values derived Introduction

Conclusions References
Tables Figures

Back Close
Full in Table 2.The individual parameters show a normal distribution (small skewness and kurtosis) with no additional secondary peaks, a potential indicator of one or more statistical outliers or a sample induced parameterization within the 19 hurricanes.The test indicates that the retrieved incidence-angle dependence is a well-defined feature of the hurricane dataset; it does not necessarily mean the dependency is actually caused by the hurricane dataset.As such additional data is required to ensure that the retrieved incidence-angle dependence is due to a physical process or an offset introduced by the sampling of hurricanes within the study.
The retrieved incidence angle dependence relationship is used only when using ECMWF forecast winds below 20 m s −1 because it was observed as a well defined feature of the ECMWF-VH distribution.Any incidence angle dependence in the wind speed above 20 m s −1 , or when utilizing literature relationships (e.g.Vachon and Wolfe, 2011), is ignored and left for future work.

Azimuthal angle dependence
The VH data can subsequently be checked for a wind direction dependence.In Fig. 7, the data is binned with respect to the wind direction and wind speed.only.It does require additional VV or HH measurements to retrieve both direction and wind speed.For comparison, the center panel shows the median VV signal collocated with the VH measurements in the top panel.The combined distribution based on the 19 hurricanes shows a chaotic image due to the VV dependence.Finally, one can check the VV-GMF wind velocity/direction distribution.In the bottom panel, the median VH signals are shown within each wind direction and wind speed bin, where the wind speed is retrieved using the GMF CMOD5.n(Hersbach, 2009 based on the collocated VV measurements and the ECMWF wind direction).Note that the maximum wind speed is lower in this latter case due to the lower sensitivity in the VV band.

VH-GMF for strong-to-severe wind speeds
Based on the various data sources it is evident that there are two wind regimes that exhibit a change in slope around 20 m s −1 .As such, it is possible to construct GMF models/relationships from both the ECMWF winds and the SFMR winds: 1. ECMWF winds collocated with RADARSAT-2 VH data results in a well-defined distribution between 7 and 37 m s −1 .In Fig. 8, the ECMWF winds vs. VH [dB] data distribution (represented by the greyscale contours) reveals a distinct trend/dependence in the low-to-strong regime (dashed-blue line) where we observe a much steeper slope when compared to the strong-to-severe (solid blue line) wind speed regime.The two linear relationships (U 10 -VH [dB]) cross at approximately 20 m s −1 .For the low-to-strong wind speed regime an incidence angle dependence has been retrieved.since it has been calibrated using drop sondes during the flights at strong and severe winds.In the low-to-strong wind speed regime additional data is available by assuming the validity of the wind speed relationship from Vachon and Wolfe (2011) or Zhang et al. (2011) which is based on collocated VH data with in-situ buoy measurements.
The four relevant linear relationships are provided in Table 3 and Fig. 8 where the lowto-strong wind speed relationship based on the ECMWF data has been corrected for incidence angle.The main difference in the two strong-to-severe regime parameterizations is an offset of 1 dB (see Fig. 8 and Table 3).This results in an approximate difference of 3.7 to 4.8 m s −1 , between VH: −21 and −30 dB respectively, for which the retrieved wind speeds based on SFMR data exceeds the ECMWF based retrieval.The apparent transition between the low-to-strong wind speeds and the strong-to-severe wind speeds has to be taken into account in the VH-GMF.Two potential methods are considered here to ensure this: -Calculate the wind speed using the relationship in the low-to-strong wind speeds and strong-to-severe wind speeds and calculate the local maximum of the two functions.This ensures that there is an easy transition from one to the other relationship, but the procedure will result in a sharp edge at the transition point.
- that this relationship does not do a good job of defining the high wind speeds found in Hurricanes.The wind speed in the top right panel shows the retrieval based on the ECMWF fits.The bottom left panel uses the SFMR data for the strong-to-severe winds and the Vachon and Wolfe relationship for the low-to-strong wind speeds.In the bottom right panel the SFMR velocities along a flight leg are plotted, the transect itself is shown as the black line in the top left panel.Overlaid are the collocated retrieved wind speeds from the three retrievals and the VV-GMF.The higher wind speeds in the bottom left panel vs. the top right panel in Fig. 9 show the influence of the offset in the VH-GMF when using the SFMR + Vachon and Wolfe fits versus the ECMWF fits.
It also provides a good indication of the uncertainties involved when using the cross polarization signals to retrieve the wind speed.Since the SFMR instrument is specifically calibrated for Hurricane winds, the Vachon and Wolfe parameterization is based on in-situ measurements and the maximum velocity on which the ECMWF fit is based is only 37 m s −1 , the proposed VH-GMF is based on the combination of the Vachon and Wolfe and SFMR fit (Table 3).The strong sensitivity of the retrieved wind speed to the exact linear fits and the need to assess an incidence angle dependence for the strong to severe winds requires additional data not available at this point.Future work should take into account a large number of additional hurricane (and typhoon) images with collocated SFMR measurements, both in time and space.The addition of new cross polarization channels to the NOAA aircraft instrument should, together with satellite data, assist in a more accurate definition of a future VH-GMF.Nonetheless, with the currently available data, one can conclude that cross-polarization signals from satellites will enable the retrieval of strong-to-severe wind speeds where VV(/HH) data has saturated.

Potential of near real-time assessment of the NOAA best track estimates
In the previous sections, a VH-GMF was presented showing a linear relationship of VH [dB] with respect to the local wind speed.Since the amount of data available is Introduction

Conclusions References
Tables Figures

Back Close
Full not yet sufficient to retrieve a unique relationship (with low error estimates) for the entire range of potential wind speeds, it is important to validate the retrieved VH-GMF.
The correspondence between the SFMR data and VH-GMF retrievals (Fig. 9) is, in part, due to the input data, as the Hurricane image was used in both the ECMWF and SFMR fitting.Even though the bootstrap tests show that the resulting fits do not depend or change due to a single Hurricane image, the comparison is not independent.Furthermore, the simple relationship found in this and previous studies suggest the potential use of the VH measurements for a near real time wind speed assessment in the case of Hurricanes and Typhoons, thereby enabling hurricane prediction centers to improve their forecasts and potential damage assessments.In this section, the VH measurements are compared to the Tropical Cyclone Best Track Information product from the NHC GIS Archive 3 .This Best Track is a subjectively-smoothed representation of the tropical cyclone's location and intensity.The intensity is described as its maximum 1 min sustained surface wind speed and minimum sea-level pressure at 6 h intervals over its lifetime.The track and intensity values are determined as a post storm assessment of all available data.The track data itself is intended for regional-scale analysis and will not show the Hurricane's erratic motions over its lifetime (see website for more detailed information on the product).The VH images and their retrieved velocities are independent of the NOAA best track assessment data, but also describe the hurricane in a very different way.The VH-GMF results are based on a high spatial resolution snapshot; whereas, the NOAA best track is a single intensity value of maximum sustained wind depending on time, and when comparing them, the variability has to be treated in a different way.For the VH-GMF, only the tail of the image distribution is important since it describes the highest wind velocities.In the case of the one-minute maximum sustained winds the most important issue is its stability in time, since this describes how stable the wind speed estimate is at the time of the RADARSAT-2 overpass.For each Hurricane image, the NOAA best track intensity is estimated by interpolating the maximum velocity to the satellite overpass time and the associated error is defined as Introduction

Conclusions References
Tables Figures

Back Close
Full the standard deviation of the 5 closest 6 h points of the NOAA track.This means that the more variable the Hurricane is in its lifetime, e.g. it is increasing or decreasing in strength within a 24 h window, the larger the local error bar.For the retrieved VH-GMF velocities, a land-sea mask is applied to disregard high VH values due to land return.If one would subsequently take the highest wind speed in the image, the value could depend on high returns from non-sea-surface targets, e.g.coastlines or small islands that weren't removed by the land-sea mask operation, ships or an occasional single high noise pixel.Instead of comparing the maximum retrieved wind speed, the 0.995 and 0.9995 percentile wind speed values are retrieved with the center defined as its mean velocity value to ensure a robust retrieved value.These values are specific for the used swath of the RADARSAT-2 instrument and could change for other cross polarization instruments.The top panel of Fig. 10 has a plot that indicates which hurricanes have SFMR data (blue dots) and those that do not have SFMR data (red dots).The linear fit (solid line) was estimated taking into account both the variability (error) estimates in x and y and has a correlation of 0.83 and a slope of 1.39.The difference in slope between the VH-GMF and the 1 : 1 line could be caused by a number of issues including the lack of an incidence angle correction, the difference in the "meaning" of the velocity value (high spatial resolution versus an averaged time signal) and finally the use of percentile values lower than 1.0.Since the VH-GMF at these wind speeds is described using a linear relationship in the strong-to-severe wind regime, one can also The width of the distribution can be decreased when more hurricane SFMR data is available in order to define the incidence-angle dependence in the strong-to-severe wind speed regime and by adding more points to this plot, i.e. images of typhoons or hurricanes.Note that this is not the same relationship as was defined earlier (Eq.3), which was defined using all available data instead of being based on statistics of the extreme values.The comparisons between the ECMWF data, SFMR and best-track wind speeds vs. the measured VH data sets are independent and bring forward two important conclusions: 1.The cross polarization radar return in C-band increases exponentially with respect to wind velocity (linear U 10 -VH [dB]) as has been found by three independent wind speed data records, but now up to 60 m s −1 .
2. A near real-time assessment of maximum sustained wind speed is possible using cross polarization measurements.This will be extremely useful for Hurricane forecast warnings and hurricane model run initialization.

Conclusions
In The upper limit of 45 m s −1 is set by the currently available collocated data.The validity of the retrieved linear relationship (U 10 -VH [dB]) has been evaluated by comparing the cross-polarized data to two additional independent wind speed datasets, i.e. the short-range ECMWF model forecast and the NOAA best estimate one-minute maximum sustained winds from their GIS database.Comparing these data sets results in the following conclusions: -The VH backscatter relationship to wind speed shows a different behavior in the low-to-strong wind speed regime (U 10 < 20 m s −1 ) compared to the strong-tosevere wind speed regime (U 10 > 20 m s −1 ).The low-to-strong wind speed regime is described by the relationship found by Vachon and Wolfe (2011) and the strongto-severe wind regime by the relationship retrieved here (Eq.3).
-The VH backscatter shows a modest incidence angle dependence.A correction for this dependence is only feasible for wind speeds up to 20 m s −1 due to the low amount of severe wind speed data.Additional data is needed to determine the incidence-angle dependence in the strong-to-severe wind regime and determine its origin (geophysical or calibration).
-There is no wind direction dependence in the VH backscatter signals.This means that ocean wind speeds can be directly retrieved from a VH channel in the C-band.
For an operational wind product, three additional VV or HH beams are required to retrieve the wind speed and direction up to 24 m s −1 .Combining the information from both the VV and VH signals should enable the calculation of wind direction for larger wind speeds.In the case of hurricane force winds, the knowledge of wind speeds is more important than the wind direction, as the former denotes the class and devastation of the hurricane.In those cases the wind direction can be defined by the general flow pattern (e.g. in the case of hurricanes and typhoons the wind rotates around the eye).Introduction

Conclusions References
Tables Figures

Back Close
Full -There is a strong correlation between the highest VH measurements within a hurricane SAR image and the one-minute maximum sustained wind speeds from the NOAA best track product and can be described by a linear relationship (U 10 -VH [dB]).Since this relationship is only based on a combination of the most extreme VH measurements, it is different from the more general VH-GMF which is found by comparing all available SFMR collocated measurements.
-An operational satellite with a VH channel will be able to retrieve, in near real-time, the one-minute maximum sustained wind speed from a hurricane or typhoon, as well as be able to make an assessment of the spatial wind speed distribution.Up to the launch of the MetOP-SG satellites this may be performed by the RADARSAT-2 satellite, as long as this is possible within its measurement schedule and science mandate.Introduction

Conclusions References
Tables Figures

Conclusions References
Tables Figures

Back Close
Full  Full 2 dB; whereas, SAR systems typically achieve < 2 dB.A main limitation of the ASCAT fan-beam scatterometer is the low sensitivity of co-polar signal strength with respect to severe wind speeds (> 25 m s −1 ).This reduces the ability to derive wind speeds during severe wind events like Hurricanes and Typhoons.In recent years, RADARSAT-2 cross-polarization signals (VH) have been studied in more detail, Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | VH [dB] = 0.592 • U 10 [m s −1 ] − 35.6 (Vachon and Wolfe) (1) VH [dB] = 0.580 • U 10 [m s −1 ] − 35.652 (Zhang).
-to-severe wind speed measurements over hurricanes using NOAA Hurricane-Hunter winds and RADARSAT-2 σ 0 A total of 19 RADARSAT-2 dual polarized SAR images (VV and VH) were acquired, through the Canadian Space Agency's Hurricane Watch program, during strong-tosevere wind events between 2008 and 2011.Hurricane transects and hurricane surface winds were also obtained from the NOAA WP-3D research aircraft carrying the airborne Stepped-Frequency Microwave Radiometer (SFMR)(Uhlhorn et al., 2007) instrument.The retrieval of wind velocity is based on nadir microwave emissions from the sea surface in 6 C-band frequencies along the flight track.The retrieval algorithm is validated using global positioning system (GPS) dropwindsonde measurements and has been shown to retrieve wind speeds within 4 m s −1 RMS error from 10 to about For both the VV and VH signals only those combined signals which were Discussion Paper | Discussion Paper | Discussion Paper | If the time difference is < 24 h the flight is separated in to different flight legs by looking for large changes in flight direction (> 90 • ).Overall most flight legs occur within 3 h of the RADARSAT-2 acquisition time.Only the Hurricane Bill image acquired on 22 August exceeded the 24 h maximum (by 3 h).Each flight leg is checked for a passage through the center of the hurricane by matching the measured velocity patterns through the hurricane eye with the position of the eye as determined in the RADARSAT-2 image.For those flight legs which contain a hurricane eye, the entire SFMR leg is shifted in longitude and latitude to match the centers, after which the VH measurements are collocated onto the measured points of the SFMR.Finally, the chosen flight legs with the collocated VH data are checked for large inconsistencies (i.e. did the aircraft clearly pass through a different part of the hurricane).Since the rotation of the hurricane is not taken into account in this procedure, a comparison between the two data sources can only be made statistically.Discussion Paper | Discussion Paper | Discussion Paper | combining a large number of flight legs with collocated VH data, the collocation errors should increase the width of the distribution without inducing a bias.In Fig. 2 a comparison of a single flight leg through hurricane Earl is shown.Plotted are the individual retrieved SFMR wind speeds (black line) measurements along the flight track and the collocated RADARSAT-2 measured VH signal (blue line).A clear correlation is visible between the two signals.By combining all the available legs from the 9 measured hurricanes the 2-D joint distribution shown in Fig. 3 is created.Due to the high instrument noise floor in the SCWA mode the lower part of the distribution below −35 dB is missing compared to the data used in the Vachon and Wolfe (2011) and Zhang and Perrie (2011) papers, which used the more sensitive (but smaller spatial extent) fine quad-polarization RADARSAT-2 mode.This affects the velocity distribution below 20 m s −1 and only a hint of the fit by Vachon and Wolfe (2011) (Eq.1; dasheddotted line) in this velocity regime is visible.Above 20 m s −1 a different slope compared to the low wind speed fit is apparent.
Discussion Paper | Discussion Paper | Discussion Paper | VH measurement from 10 to 45 m s −1 .The exact transition from strong to severe wind speeds, i.e. the transition from the Vachon and Wolfe relationship to the severe wind speed regime, can not be assessed at this point based on these data only.The SFMR distribution is well defined up to 40 m s −1 with a correlation of 0.70 and can be related to VH in between 21 and 40 m s −1 as VH(21 < U 10 < 40 m s −1 ) = 0.218U 10 − 29.07.(3) Theoretical derivations by Hwang et al. (2010) and Voronovich and Zavorotny (2011), showed that the cross polarization signals should depict an incidence angle dependence.This was also indicated by measurements using a limited set of the RADARSAT-Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | 4) were fitted to the VH [dB] measurements versus ECMWF velocity for different incidence angle bins.The linear fits, Figures Back Close Full Screen / Esc Printer-friendly Version Interactive Discussion Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | The color scale shows the median value of the measured VH (top panel) within each bin.It is observed in the plot that the VH versus wind speed shows no dependence to the wind direction angles, as observed byHwang et al. (2010),Zhang and Perrie (2011)  andVachon and Wolfe (2011) in the < 20 m s −1 regime.The same observation is found/confirmed in the strong-to-severe (> 20 ms) wind speed regime.In contrast to the earlier cited papers and this work, cross polarization data from L-band radar systems (1-2 GHz) do show a wind direction dependence(Yueh et al., 2010).The solid contour lines shown in Fig. 7 (top panel) depicts the VH cross polarization signal variation derived in the L-band, which is much larger than any variations depicted by the colored contours.It is unknown why L-band exhibits a dependence on wind direction and C-band does not for crosspolarized data.The lack of dependence seen in the C-band cross polarization makes it a more robust wavelength for absolute wind speed retrievals using a single observation 7959 Discussion Paper | Discussion Paper | Discussion Paper |

2.
In Fig. 8, the distribution of the SFMR wind speed vs. the VH data shows a very similar pattern (depicted by the black contour lines and the associated linear fit is the red solid line) compared to the ECMWF distribution in the strong-to-severe wind speed regime.The lower part of the distribution is incomplete in the U 10 < 20 m s −1 and VH < 30 dB region.The SFMR retrieved velocity is a superior product Discussion Paper | Discussion Paper | Discussion Paper | Create a continuous function by combining the two relationships in the following way U 10 = ((U 10 [LS]) p + (U 10 [SE]) p ) 1/p , where LS stands for low-to-strong relationship, SE for the strong-to-severe relationship and p is equal to 10.The high value of p ensures that the combined velocity is weighted towards the largest of the two function values.In the transition regime both functions have nonnegligible weights, resulting in a smooth transition with a continuous derivative.In Fig. 9, the different relationships of the wind speed are compared to the VH measurements for hurricane Earl.In the top left panel, the Vachon and Wolfe parameterization is adopted.This results in wind speeds lower than 30 m s −1 , which indeed indicates Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | directly plot the VH(0.995,0.9995) percentile values against the NOAA best estimate as is shown in the bottom panel of Fig. 10.As expected, the plot indicates a linear relationship with a correlation factor of 0.83, U 10 (NOAA) = 170.69+ 6.20 VH[0.995] + VH[0value in square brackets defines the cross polarization percentile value to be used from a single RADARSAT-2 image which has been land masked and incorporates the eye of the hurricane to ensure that the strongest wind speeds are within the image.Discussion Paper | Discussion Paper | Discussion Paper | this work, cross-polarized (VH) backscatter signals from RADARSAT-2 C-band SAR imagery acquired during severe wind speed (hurricane) events were compared to the collocated and time co-incident SFMR wind measurements acquired by NOAA's hurricane-hunter aircraft.From this data set, a relationship for wind speed retrieval using VH-polarization backscatter for strong-to-severe wind speeds (i.e.20 m s −1 < U 10 < 45 m s −1 ) was derived.Within this wind speed regime, the cross-polarized data showed no distinguishable loss of sensitivity and, as such, cross-polarized signals are a good candidate for the retrieval of severe wind speeds from satellite instruments.The VH signals increase exponentially with wind speed and are described as linear relationships of U 10 vs. VH [dB].Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Slade, B.: RadarSAT-2 product description, Tech.rep., MacDonald, Dettwiler and Associates Ltd., available at: gs.mdacorporation.com/SatelliteData/Radarsat2/Products.aspx(last access: August 2013), 2011.7949 Stoffelen, A. and Portabella, M.: On Bayesian scatterometer wind inversion, IEEE T. Geosci.

Fig. 2 .Fig. 3 .Fig. 4 .Fig. 6 .Fig. 8 .Fig. 9 .Fig. 10 .
Fig. 1.VH (top) and VV(bottom) σ 0 measurements for hurricane Earl (2 September 2010).The hurricane wind speed structure in the VV band is not clearly visible due to the incidence and azimuthal angle dependencies.The hurricane structure shows a well-defined eye and wind pattern in the VH band.Note the large difference in signal strength resulting in signals close to the noise level of the instrument.

Table 1 .
RADARSAT-2 hurricane images used in this study with acquisition dates and times.The NOAA aircraft made flights through a number of hurricanes for which SFMR wind velocities were measured.The table indicates the number of collocated flight tracks performed within the RADARSAT-2 images which were close enough in time to be representative.

Table 2 .
Fitted parameters for the retrieved functions A and B (Eq. 7).