Evaluation of HY-2A Scatterometer Ocean Surface Wind Data during 2012–2018

This study focuses on the evaluation of global Haiyang-2A satellite scatterometer (HSCAT) operational wind products from 2012 to 2018. In order to evaluate HSCAT winds, HSCAT operational wind products were collocated with buoy measurements and rainfall data. Error varieties under different atmospheric stratification and rainfall conditions were taken into consideration. After data quality control, the average bias and root mean square error (RMSE) between buoys and HSCAT data were 0.1 m/s and 1.3 m/s for wind speed, and 1° and 27° for wind direction, respectively. Especially, the varieties of the wind direction difference change a lot under non-neutral atmospheric conditions. HSCAT wind speeds are overestimated with an increasing rainfall rate while wind directions tend to be perpendicular to buoys’. In brief, the HSCAT wind product qualities are not stable during 2012 to 2018, especially for the data in 2015 and 2016. Atmospheric stratification and rain effects should be considered in wind retrieval and marine application.


Introduction
The sea surface wind field is one of the direct factors affecting small-scale and large-scale seawater movement and plays an important role in marine environmental monitoring, weather forecasting, as well as other marine and atmospheric sciences studies. Owing to the spatial and temporal constraints, traditional wind measurement methods, such as ships, buoys, weather stations, and other in situ observations, are not satisfied with the requirements of high-quality and continuous wind observations for ocean research. Since the 1970s, satellite observations have become an indispensable approach to collect sea surface wind field data [1].
There are four main satellite microwave sensors for observing the global sea surface winds: Microwave scatterometer, microwave radiometer, altimeter, and synthetic aperture radar [1]. The microwave scatterometer is the main sensor for providing wind data with its all-weather and all-time capacity [2]. The earliest satellite scatterometer is the Ku-band fan beam Seasat-A Satellite Scatterometer (SASS), launched by National Aeronautics and Space Administration (NASA) in 1978. SASS opened a new era of satellite-observing sea surface wind fields [3]. Afterwards, series C-band and Ku-band scatterometers were successfully launched [4][5][6][7][8][9]. The Ku-band scatterometer of the first Chinese microwave ocean environment satellite Haiyang-2A (HSCAT), launched on 16 August 2011, is a pencil beam scatterometer with a 1700-km swath and 25-km spatial resolution [10].
Evaluating the quality of scatterometer data is indispensable for sea surface wind field study. There are abundant evaluation investigations of different sensors studied by worldwide researchers. Masuko et al. (2000) evaluated NSCAT wind vectors with Japan Meteorological Agency (JMA) buoys, products with Level 2B (L2B) winds, numerical weather prediction winds, and buoy winds. The results showed that the new retrieval algorithm can improve the spatial consistency of UHR winds [25].
At present, studies of satellite scatterometer wind field evaluation are mostly regional and in a short-time period. Our research focused on a global and long-term evaluation to study the stability of the HSCAT-operated wind product during 2012-2018. Additionally, atmospheric stratification and rain effects on scatterometer wind retrieval were taken into account.
In this paper, we used HSCAT operational Level 2B wind products from January 2012 to December 2018 and buoy data provided by the European Centre for Medium-Range Weather Forecasts (ECMWF) MARS archive to evaluate and analyze HSCAT wind products. In addition, we also combined HSCAT data with the Special Sensor Microwave Imager Sounder (SSMIS) rainfall data to evaluate the rain effects on HSCAT winds. Scatterometer winds, buoy winds, and SSMIS rainfall data as well as their collocation methods are introduced in Section 2; the main product analysis, atmospheric stratification, and rain effect results are discussed in Section 3; and summaries of this study are provided in Section 4.  Table 1 [10]. The Ku-band scatterometer on HY-2A (HSCAT) operates at 13.3 GHz with the pencil beam conical scanning mode. The inner antenna works in horizontal polarization (HH) at a 41 • incidence angel over a 1350-km swath, and the outer antenna works in vertical polarization (VV) at a 48 • incidence angel over a 1700-km swath. HSCAT Level 2B (L2B) wind vectors are retrieved based on NSCAT-2 Geophysical Model Function (GMF) by maximum likelihood estimation (MLE) and the circle median filter is used to remove wind direction ambiguities [26]. The L2B products provide daily sea surface information in orbits with a 25-km resolution, and every single data file consists of 1624 rows and 76 columns. The wind direction in L2B products is defined as 0 degrees when wind is blowing toward the north and increases in the clockwise [27]. The HSCAT L2B wind data from 2012 to 2018 were used in this paper.

ECMWF MARS Buoy Data
In situ wind observations are used to evaluate HSCAT winds. Buoy wind vectors are distributed by the Global Telecommunication System (GTS), which have been retrieved and quality controlled by the ECMWF MARS archive. The temporal resolution of buoy winds is one hour by averaging the wind data over 10 min. The wind direction is defined as 0 degrees when wind is blowing from north and increases in the clockwise direction [23]. The data of 248 buoy stations consisting of 30 nearshore buoy Remote Sens. 2019, 11, 2968 4 of 22 stations (within 50 km from the coast) and 218 offshore stations (more than 50 km away from the coast) were collected in this study and their locations are shown in Figure 1.
Remote Sens. 2019, 11, x FOR PEER REVIEW 4 of 22 and increases in the clockwise direction [23]. The data of 248 buoy stations consisting of 30 nearshore buoy stations (within 50 km from the coast) and 218 offshore stations (more than 50 km away from the coast) were collected in this study and their locations are shown in Figure 1. As scatterometer wind data are regarded as 10-m neutrally stable wind, all buoy winds were converted to 10-m neutral winds (marked as ⃗ ) by using the LKB model proposed by Liu, Katsaros and Businger in 1979 [28,29]. However, the scatterometer measures the sea surface normalized radar cross-section (NRCS), which is affected by wind friction directly. Therefore, we converted ⃗ to stress-equivalent wind (marked as ⃗ ). The relation between the stress-equivalent wind and neutral winds is [30]: where is the air density (kg/m 3 ); and 〈 〉 is the global average air density and equals to 1.225 kg/m 3 here.

SSMIS Data
To study rain effects on HSCAT wind retrieval, we used the Special Sensor Microwave Imager Sounder (SSMIS) rainfall data to extract HSCAT wind data contaminated by rain. The Special Sensor Microwave Imager (SSM/I) and SSMIS are series microwave radiometers payloaded on the nearpolar orbiting Defense Meteorological Satellite Program (DMSP) satellites since 1987. SSMIS works at 19.35, 22.235, 37, and 85.5 GHz. The rainfall data are provided with a gridded dataset of a 0.25° spatial resolution generated by Remote Sensing System (REMSS) [31]. To match with the HSCAT descending time, we chose the DMSP F17 platform, of which the descending time (6:37 a.m.) is the closest to HSCAT's (6:00 a.m.) and quite stable during the chosen time period. The ascending and descending equatorial crossing time for different instruments can be inquired from REMSS [32].

Data Collocation Method
Before collocation, data selection and quality control need to be done firstly. Because of the dependences between the azimuth dispersion of rotating scan scatterometer and the accuracies of wind vector retrieval, small azimuth dispersion causes a lower quality in the winds [33][34][35]. Crosstrack wind cell can be divided as the outer swath (OS, or outer region (OR)) and inner swath (IS) and As scatterometer wind data are regarded as 10-m neutrally stable wind, all buoy winds were converted to 10-m neutral winds (marked as − −− → V 10n ) by using the LKB model proposed by Liu, Katsaros and Businger in 1979 [28,29]. However, the scatterometer measures the sea surface normalized radar cross-section (NRCS), which is affected by wind friction directly. Therefore, we converted − −− → V 10n to stress-equivalent wind (marked as − −− → V 10s ). The relation between the stress-equivalent wind and neutral winds is [30]: where ρ is the air density (kg/m 3 ); and ρ is the global average air density and equals to 1.225 kg/m 3 here.

SSMIS Data
To study rain effects on HSCAT wind retrieval, we used the Special Sensor Microwave Imager Sounder (SSMIS) rainfall data to extract HSCAT wind data contaminated by rain. The Special Sensor Microwave Imager (SSM/I) and SSMIS are series microwave radiometers payloaded on the near-polar orbiting Defense Meteorological Satellite Program (DMSP) satellites since 1987. SSMIS works at 19.35, 22.235, 37, and 85.5 GHz. The rainfall data are provided with a gridded dataset of a 0.25 • spatial resolution generated by Remote Sensing System (REMSS) [31]. To match with the HSCAT descending time, we chose the DMSP F17 platform, of which the descending time (6:37 a.m.) is the closest to HSCAT's (6:00 a.m.) and quite stable during the chosen time period. The ascending and descending equatorial crossing time for different instruments can be inquired from REMSS [32].

Data Collocation Method
Before collocation, data selection and quality control need to be done firstly. Because of the dependences between the azimuth dispersion of rotating scan scatterometer and the accuracies of wind vector retrieval, small azimuth dispersion causes a lower quality in the winds [33][34][35]. Cross-track wind cell can be divided as the outer swath (OS, or outer region (OR)) and inner swath (IS) and IS consists of a sweet region (SR) and nadir region (NR), as shown in Figure 2. Since only vertical polarization exists in the outside edge region of the swath, HSCAT OS data were not included in the collocation dataset. For buoy data, some of the buoys changed their moored locations during measurement. We selected the buoy data measured at the same time but with different locations, and averaged the buoy data as one measurement. As for SSMIS rainfall data, we removed all the invalid data marked by REMSS. IS consists of a sweet region (SR) and nadir region (NR), as shown in Figure 2. Since only vertical polarization exists in the outside edge region of the swath, HSCAT OS data were not included in the collocation dataset. For buoy data, some of the buoys changed their moored locations during measurement. We selected the buoy data measured at the same time but with different locations, and averaged the buoy data as one measurement. As for SSMIS rainfall data, we removed all the invalid data marked by REMSS. After data quality control, we made the buoy and scatterometer collocation datasets first. HSCAT wind field data and buoy wind data were collocated within 25 √2 ⁄ km (about 17.7 km) and 30 min. In this study, wind speeds were used within 0.1 to 25 m/s and HSCAT wind directions were turned 180° to keep with the buoy wind directions before analysis. As for rainfall data, we used the location of HSCAT wind data to find F17 SSMIS rain data within 1 h.
Due to the sea surface status being affected by the terrain in shallow water and the contamination of the island, NRCS is greatly affected near the coast. We classified buoy data as the offshore dataset and nearshore dataset according to the distance between buoy stations and coastline. There are 30 nearshore buoy stations within 50 km from the coastline and others are offshore buoy stations, as shown in Figure 1. Besides, it was found that the amount of nearshore data is about 4% of the offshore data. RMSEs of the nearshore data were about 1.58 m/s and 35.52°while the results of the offshore data were 1.30 m/s and 27.40°during the whole study period. The deviation of the nearshore data is large and therefore only the offshore dataset was used in this study.

Statistical Analysis of HSCAT Wind Variations from 2012 to 2018 Year by Year
To evaluate the quality of the HSCAT wind data, HSCAT winds were collocated with buoy data year by year from 2012 to 2018. The spatial-temporal collocation windows are 17.7 km and 30 min as introduced in Section 2.

The Distribution of Collocated Data and Statistical Parameters
The whole-year HSCAT wind data were compared with the contemporaneous buoy wind data and the wind speed and wind direction collocation contours of histograms are shown in Figure 3 from 2012 to 2018. Figure 3a  After data quality control, we made the buoy and scatterometer collocation datasets first. HSCAT wind field data and buoy wind data were collocated within 25/

√
2 km (about 17.7 km) and 30 min. In this study, wind speeds were used within 0.1 to 25 m/s and HSCAT wind directions were turned 180 • to keep with the buoy wind directions before analysis. As for rainfall data, we used the location of HSCAT wind data to find F17 SSMIS rain data within 1 h.
Due to the sea surface status being affected by the terrain in shallow water and the contamination of the island, NRCS is greatly affected near the coast. We classified buoy data as the offshore dataset and nearshore dataset according to the distance between buoy stations and coastline. There are 30 nearshore buoy stations within 50 km from the coastline and others are offshore buoy stations, as shown in Figure 1. Besides, it was found that the amount of nearshore data is about 4% of the offshore data. RMSEs of the nearshore data were about 1.58 m/s and 35.52 • while the results of the offshore data were 1.30 m/s and 27.40 • during the whole study period. The deviation of the nearshore data is large and therefore only the offshore dataset was used in this study.

Statistical Analysis of HSCAT Wind Variations from 2012 to 2018 Year by Year
To evaluate the quality of the HSCAT wind data, HSCAT winds were collocated with buoy data year by year from 2012 to 2018. The spatial-temporal collocation windows are 17.7 km and 30 min as introduced in Section 2.

The Distribution of Collocated Data and Statistical Parameters
The whole-year HSCAT wind data were compared with the contemporaneous buoy wind data and the wind speed and wind direction collocation contours of histograms are shown in Figure 3      The quantitative statistics of the wind vector bias in the collocated winds are given in Table 2. Invalid data were removed, the speeds of which are twice as high as the speed's standard deviation(σ) by using the 2σ quality control method proposed by Lin et al. (2013) [26]. The statistical   The quantitative statistics of the wind vector bias in the collocated winds are given in Table 2. Invalid data were removed, the speeds of which are twice as high as the speed's standard deviation(σ) by using the 2σ quality control method proposed by Lin et al. (2013) [26]. The statistical parameters were calculated using the following equation (here, '< >' represents the statistical average) [9,36]: Wind direction statistical parameters were calculated by using Equations (4)-(7) to ensure the continuity at 0 • and 360 • [37]: The statistical parameters in Table 2 explain the quantification of the variation trends shown in All in all, the HSCAT wind speed accuracy can satisfy the scatterometer mission requirements (2 m/s) but the wind direction accuracy is a little worse. Besides, the wind data in 2015 and 2016 have the greatest offsets. A strict quality control is needed in the application of these datasets.

The Probability Distribution Functions (PDFs) of Wind Speeds and Directions
In order to analyze the entire distributions of the collocated data between HSCAT and buoys, we plotted the probability distribution functions (PDFs) of wind speeds and directions for each year. The PDFs are shown in Figure 4 in which a-g represent the PDF curves for each year from 2012 to 2018; number 1 and 2 represent the wind speed and direction PDFs; and the dotted line is the PDF for the buoy data and the solid line for HSCAT PDF.  For wind speed PDFs, the peak values of the buoy data are about 6 to 7 m/s and HSCAT data are distributed at 6 to 8 m/s except in 2016, which is less than 5 m/s. Within the wind speed ranges from 0 to 4 m/s, the numbers of the HSCAT wind speed are slightly higher except in 2015 and 2016; within 4 to 8 m/s, the numbers of the HSCAT wind speed are about 2% lower than the buoys; within 8 to 13 m/s, the numbers of the HSCAT wind speed are generally higher; and when speeds are higher than 13 m/s, the HSCAT speed PDFs are basically consistent with the wind speed PDFs of the buoy data.
For wind direction PDFs, both the HSCAT and buoy peak values are distributed stably about 90 to 110° among seven years. The numbers of the HSCAT wind direction are a little higher when the wind directions are less than 90°; in the range 90-100°, the numbers of the HSCAT direction at the peak values are less than the buoys; the numbers of the HSCAT direction are lower in the range 120-220°, and then become higher in 220 to 330° but become lower again when directions are larger than 330°. From 2014, the HSCAT wind direction PDFs show two little peaks (70-80° and 110-120°) on the main central peak region but disappear in 2018.
In a word, for the annual wind vector PDF distributions, the numbers of HSCAT and buoy wind data are basically the same, but the number distributions vary in different wind speed and wind direction bins.

The Trends of Residual Variations of Wind Speed and Wind Direction
The dependencies of the wind speed and direction residual ( − ) from 2012 to 2018 are shown in Figure 5. In order to study wind speed and direction bias variations in different speed bins, the residuals were calculated in the speed bins of 1 m/s. To avoid a binning effect, we used the average wind speed, ( + ) 2 ⁄ , as the binning wind speed [38]. All collocated data were quality controlled by the 2σ method . Figure 5a and b represent the wind speed and direction residual distributions; the error bars represent the wind data biases and RMSEs, respectively. For wind speed PDFs, the peak values of the buoy data are about 6 to 7 m/s and HSCAT data are distributed at 6 to 8 m/s except in 2016, which is less than 5 m/s. Within the wind speed ranges from 0 to 4 m/s, the numbers of the HSCAT wind speed are slightly higher except in 2015 and 2016; within 4 to 8 m/s, the numbers of the HSCAT wind speed are about 2% lower than the buoys; within 8 to 13 m/s, the numbers of the HSCAT wind speed are generally higher; and when speeds are higher than 13 m/s, the HSCAT speed PDFs are basically consistent with the wind speed PDFs of the buoy data.
For wind direction PDFs, both the HSCAT and buoy peak values are distributed stably about 90 to 110 • among seven years. The numbers of the HSCAT wind direction are a little higher when the wind directions are less than 90 • ; in the range 90-100 • , the numbers of the HSCAT direction at the peak values are less than the buoys; the numbers of the HSCAT direction are lower in the range 120-220 • , and then become higher in 220 to 330 • but become lower again when directions are larger than 330 • . From 2014, the HSCAT wind direction PDFs show two little peaks (70-80 • and 110-120 • ) on the main central peak region but disappear in 2018.
In a word, for the annual wind vector PDF distributions, the numbers of HSCAT and buoy wind data are basically the same, but the number distributions vary in different wind speed and wind direction bins.

The Trends of Residual Variations of Wind Speed and Wind Direction
The dependencies of the wind speed and direction residual (SCAT − BUOY) from 2012 to 2018 are shown in Figure 5. In order to study wind speed and direction bias variations in different speed bins, the residuals were calculated in the speed bins of 1 m/s. To avoid a binning effect, we used the average wind speed, (spd SCAT + spd BUOY )/2, as the binning wind speed [38]. All collocated data were quality controlled by the 2σ method . Figure 5a The collocated wind speed data are mostly distributed in the range 3-8 m/s. Average wind speed biases generally increase as wind speeds increase. When wind speeds are in the range 3-11 m/s, wind speed biases are approximately equal to 0. The wind speed RMSEs are less than 2 m/s and have no obvious changes with the increasing speeds. When wind speeds are higher than 10 m/s, wind speed biases are positive; that is, HSCAT wind speeds are higher than buoy wind speeds. Wind speed biases and RMSEs increase with the increasing wind speeds and the average bias and RMSE are about 2 and 3 m/s when the wind speed is 20 m/s. However, the increasing trend of the speed biases and RMSEs can be related with the decreasing collocation numbers. Generally, less data can cause a large offset.
In terms of wind direction, the whole distributions of wind direction bias remain nearly unchanged around zero as wind speeds increase except the bias at 23 m/s (about 120°), which may be caused by the small collocation numbers at the high wind speed region. When wind speeds are less than 5 m/s, the wind direction RMSEs decrease with the increasing speeds and drop from 60° to 20° rapidly. When speeds are higher than 5 m/s, the wind direction RMSEs remain nearly unchanged around 18°.
Both the wind speed and wind direction average biases are near 0 and the RMSEs are less than 2 m/s and 20° at moderate speeds, which means HSCAT wind products basically meet the scatterometer mission requirement. However, larger offsets are found in high speeds, which may be caused by the small numbers of collocated wind data.

HSCAT Wind Field Data Monthly Variations during 2012 to 2018
To further study HSCAT wind field data variations over time, wind speed and direction residual monthly variations are shown in Figure 6. Figure 6a   The collocated wind speed data are mostly distributed in the range 3-8 m/s. Average wind speed biases generally increase as wind speeds increase. When wind speeds are in the range 3-11 m/s, wind speed biases are approximately equal to 0. The wind speed RMSEs are less than 2 m/s and have no obvious changes with the increasing speeds. When wind speeds are higher than 10 m/s, wind speed biases are positive; that is, HSCAT wind speeds are higher than buoy wind speeds. Wind speed biases and RMSEs increase with the increasing wind speeds and the average bias and RMSE are about 2 and 3 m/s when the wind speed is 20 m/s. However, the increasing trend of the speed biases and RMSEs can be related with the decreasing collocation numbers. Generally, less data can cause a large offset.
In terms of wind direction, the whole distributions of wind direction bias remain nearly unchanged around zero as wind speeds increase except the bias at 23 m/s (about 120 • ), which may be caused by the small collocation numbers at the high wind speed region. When wind speeds are less than 5 m/s, the wind direction RMSEs decrease with the increasing speeds and drop from 60 • to 20 • rapidly. When speeds are higher than 5 m/s, the wind direction RMSEs remain nearly unchanged around 18 • .
Both the wind speed and wind direction average biases are near 0 and the RMSEs are less than 2 m/s and 20 • at moderate speeds, which means HSCAT wind products basically meet the scatterometer mission requirement. However, larger offsets are found in high speeds, which may be caused by the small numbers of collocated wind data.

HSCAT Wind Field Data Monthly Variations during 2012 to 2018
To further study HSCAT wind field data variations over time, wind speed and direction residual monthly variations are shown in Figure 6. Figure 6a The collocated wind speed data are mostly distributed in the range 3-8 m/s. Average wind speed biases generally increase as wind speeds increase. When wind speeds are in the range 3-11 m/s, wind speed biases are approximately equal to 0. The wind speed RMSEs are less than 2 m/s and have no obvious changes with the increasing speeds. When wind speeds are higher than 10 m/s, wind speed biases are positive; that is, HSCAT wind speeds are higher than buoy wind speeds. Wind speed biases and RMSEs increase with the increasing wind speeds and the average bias and RMSE are about 2 and 3 m/s when the wind speed is 20 m/s. However, the increasing trend of the speed biases and RMSEs can be related with the decreasing collocation numbers. Generally, less data can cause a large offset.
In terms of wind direction, the whole distributions of wind direction bias remain nearly unchanged around zero as wind speeds increase except the bias at 23 m/s (about 120°), which may be caused by the small collocation numbers at the high wind speed region. When wind speeds are less than 5 m/s, the wind direction RMSEs decrease with the increasing speeds and drop from 60° to 20° rapidly. When speeds are higher than 5 m/s, the wind direction RMSEs remain nearly unchanged around 18°.
Both the wind speed and wind direction average biases are near 0 and the RMSEs are less than 2 m/s and 20° at moderate speeds, which means HSCAT wind products basically meet the scatterometer mission requirement. However, larger offsets are found in high speeds, which may be caused by the small numbers of collocated wind data.

HSCAT Wind Field Data Monthly Variations during 2012 to 2018
To further study HSCAT wind field data variations over time, wind speed and direction residual monthly variations are shown in Figure 6. Figure 6a

HSCAT Wind Field Data Variations in 24 Hours
The wind speed and direction residual variations during 24 h are shown in Figure 7. These statistics were gathered by the buoy time. Figure 7a Wind direction biases are distributed stably around zero and the RMSEs of the wind direction are about 25°. However, from the second half of 2015, the wind direction biases tend to a 2° positive deviation and RMSEs also increase and increase to 30° in the end of 2016, which are far away from the mission requirement (20°). In 2017, the direction biases and RMSEs both decrease to about 2° and 20° in 2018. The variations of the wind direction are also similar with the statistical parameters mentioned above.

HSCAT Wind Field Data Variations in 24 Hours
The wind speed and direction residual variations during 24 h are shown in Figure 7. These statistics were gathered by the buoy time. Figure 7a and b represent the wind speed and direction error bar graphs and collocated data histograms during 2012-2018, respectively. The error bars indicate the wind data biases and RMSEs.  Wind speed bias variations are stably distributed around 0 in 24 h except in the periods that have less collocation pairs and wind speed RMSEs are generally less than 2 m/s. Wind direction biases are also distributed around 0 and fluctuate within ±2°. However, the RMSEs of the direction are stably distributed around 20 to 30°, which does not satisfy the mission requirement. This suggests that the wind speeds of HSCAT can provide stable speed products in one day with accurate qualities. Unlike the stable daily wind speed residual variations, the wind direction products of HSCAT in one day are not stable or accurate enough to be applied in marine research and these are also similar to the statistical parameters mentioned above.  This suggests that the wind speeds of HSCAT can provide stable speed products in one day with accurate qualities. Unlike the stable daily wind speed residual variations, the wind direction products of HSCAT in one day are not stable or accurate enough to be applied in marine research and these are also similar to the statistical parameters mentioned above.

The Sea Surface Temperature and Air Temperature Impact on HSCAT Wind Products
Although scatterometers can measure sea surface parameters with all-weather, all-time capabilities, the actual measurements are the sea surface normalized radar cross-section (NRCS), which is affected by both sea and atmospheric conditions. To study the sea and atmospheric parameter impacts on HSCAT wind field retrieval, we classified the collocated data into two groups according to the difference between the sea surface temperature and air temperature measured by buoys: The neutral group (the absolute values of temperature difference are less than or equal to 1 • C) and non-neutral group (the absolute values of temperature difference are higher than 1 • C). After data quality control, we obtained 217,359 pairs of collocations from 2012 to 2018 consisting of 129,604 pairs of the neutral group and 87,755 pairs of the non-neutral group.

The Probability Distribution Functions (PDFs) of Wind Speeds and Directions
We plotted the PDFs of the wind speed and direction for different groups. The PDFs are shown in Figure 8 in which (a) and (b) represent wind speed and direction PDFs; number 1, 2, and 3 represent the neutral group, non-neutral group, and all collocated data PDFs; and the dotted line is the PDF for the buoy data and the solid line for HSCAT PDF, respectively. Although scatterometers can measure sea surface parameters with all-weather, all-time capabilities, the actual measurements are the sea surface normalized radar cross-section (NRCS), which is affected by both sea and atmospheric conditions. To study the sea and atmospheric parameter impacts on HSCAT wind field retrieval, we classified the collocated data into two groups according to the difference between the sea surface temperature and air temperature measured by buoys: The neutral group (the absolute values of temperature difference are less than or equal to 1 ℃) and non-neutral group (the absolute values of temperature difference are higher than 1 ℃). After data quality control, we obtained 217,359 pairs of collocations from 2012 to 2018 consisting of 129,604 pairs of the neutral group and 87,755 pairs of the non-neutral group.

The Probability Distribution Functions (PDFs) of Wind Speeds and Directions
We plotted the PDFs of the wind speed and direction for different groups. The PDFs are shown in Figure 8 in which (a) and (b) represent wind speed and direction PDFs; number 1, 2, and 3 represent the neutral group, non-neutral group, and all collocated data PDFs; and the dotted line is the PDF for the buoy data and the solid line for HSCAT PDF, respectively. As for wind speed, HSCAT PDFs are consistent with the buoy PDFs except the numbers of HSCAT winds are less at the peak values (6-8 m/s). The HSCAT speed PDF is not in good agreement with the buoy PDF at 6 to 8 m/s for the neutral group while the non-neutral group is more consistent. In the neutral group, HSCAT has more collocation numbers when speeds are less than 5 m/s but less numbers than the buoys at 6 to 8 m/s. However, in the non-neutral group, HSCAT has less numbers  . (a,b) represent the wind speed and direction PDFs; the number 1, 2, and 3 represent the neutral group, non-neutral group, and all collocation data PDFs; the dotted line is the PDF for the buoy data and the solid line is for HSCAT PDF, respectively.
As for wind speed, HSCAT PDFs are consistent with the buoy PDFs except the numbers of HSCAT winds are less at the peak values (6-8 m/s). The HSCAT speed PDF is not in good agreement with the buoy PDF at 6 to 8 m/s for the neutral group while the non-neutral group is more consistent. In the neutral group, HSCAT has more collocation numbers when speeds are less than 5 m/s but less numbers than the buoys at 6 to 8 m/s. However, in the non-neutral group, HSCAT has less numbers when speeds are less than 5 m/s. When speeds are higher than 10 m/s, both the neutral and non-neutral group have consistent collocation results.
As for wind direction, the HSCAT and buoy wind directions are mainly distributed in the range 60-120 • . Besides, HSCAT has another peak value distributed in the range 240-300 • . In the non-neutral group, there is a double peak characteristic of both the HSCAT and buoy PDFs, and the HSCAT PDFs have larger peak values. In the neutral group, HSCAT has more collocation numbers when the directions are less than 70 • but less numbers at 90 • , where the buoy direction collocation numbers are the maximum.
Generally, the consistency of speed collocation in the neutral group is worse than the non-neutral group but the direction collocation consistency is better in the neutral group. Since the neutral group is the main part of all the collocated data, the overall PDF distribution during 2012 to 2018 is similar to the neutral group. Figure 9 shows the variations of the wind speed and direction residual on the average wind speed from 2012 to 2018 under different stratification conditions. We also used the average wind speed as the binning wind speed with 1-m/s intervals to avoid the binning effect. In Figure 9, a and b represent the wind speed and direction residual distributions; 1, 2, and 3 represent the neutral group, non-neutral group, and all collocation data residual distributions; and the error bars represent the wind data biases and RMSEs, respectively. when speeds are less than 5 m/s. When speeds are higher than 10 m/s, both the neutral and nonneutral group have consistent collocation results. As for wind direction, the HSCAT and buoy wind directions are mainly distributed in the range 60-120°. Besides, HSCAT has another peak value distributed in the range 240-300°. In the non-neutral group, there is a double peak characteristic of both the HSCAT and buoy PDFs, and the HSCAT PDFs have larger peak values. In the neutral group, HSCAT has more collocation numbers when the directions are less than 70° but less numbers at 90°, where the buoy direction collocation numbers are the maximum.

The Trends of Residual Variations of Wind Speed and Wind Direction
Generally, the consistency of speed collocation in the neutral group is worse than the nonneutral group but the direction collocation consistency is better in the neutral group. Since the neutral group is the main part of all the collocated data, the overall PDF distribution during 2012 to 2018 is similar to the neutral group. Figure 9 shows the variations of the wind speed and direction residual on the average wind speed from 2012 to 2018 under different stratification conditions. We also used the average wind speed as the binning wind speed with 1-m/s intervals to avoid the binning effect. In Figure 9, a and b represent the wind speed and direction residual distributions; 1, 2, and 3 represent the neutral group, non-neutral group, and all collocation data residual distributions; and the error bars represent the wind data biases and RMSEs, respectively.  For both the neutral and non-neutral group, the average speed biases are distributed around 0 in 0 to 12 m/s and the biases increase with increasing speeds. The speed RMSEs of the neutral group are about 1 m/s, but the RMSEs of the non-neutral group are larger (about 2 m/s). When the wind speed is higher than 12 m/s, the average speed biases start increasing in both groups. The speed biases of the neutral group are about 1.5 m/s and the RMSEs increase to about 1.5 m/s in the range 16-20 m/s. Additionally, in this range, the non-neutral group's average wind speed biases and RMSEs are larger. When the speed equals 20 m/s, the bias and RMSE increase to about 2 m/s and 2.5 to 3 m/s, respectively.

The Trends of Residual Variations of Wind Speed and Wind Direction
For wind direction, the biases remain nearly unchanged around zero. In 0 to 5 m/s, both RMSEs of the neutral and non-neutral groups decrease from 45 • to 20 • . When the wind speed is higher than 5 m/s, the RMSEs of the neutral group are about 25 • while the RMSEs of the non-neutral group are about 30 • . Both the wind speed and direction biases of the non-neutral group have large changes at 23 m/s, which may be connected with the small collocation numbers at 23 m/s.

The Statistical Parameters of Collocated Data
For further study of the winds' residual variations under neutral and non-neutral conditions, the quantitative statistics of the wind vector difference between HSCAT and buoys are given in Table 3. Table 3. The statistical parameters of the wind vector bias between HSCAT and buoys under different stratification conditions. After data quality control, generally, the number of the neutral group is 1.3 to 1.8 times the non-neutral group number. Those statistical parameters also validate the results in Sections 3.3.1 and 3.3.2 that the winds' average biases and RMSEs of the neutral group are less than the non-neutral group; both groups of the HSCAT wind speed can meet the mission requirement but the RMSEs of the wind direction are much higher than 20 • , especially in the non-neutral group, which can reach 30 • .

BMSE-Non-
For wind speed, the biases of the non-neutral group are about 0.2 m/s higher than the neutral group except in 2016, which is about 0.4 m/s higher, but they all meet the mission requirement. In 2013, 2015, and 2016, both the neutral and non-neutral groups' offsets are much larger than other years. Except in 2015, the average wind speed bias of the neutral group is similar to other years but much larger under non-neutral conditions. Those large offsets decrease in 2017 and 2018.
For wind direction, the average biases are quite small from 2012 to 2015. However, after 2015, the average biases are more than 1 • and the RMSEs are quite large whether under neutral conditions or non-neutral conditions, although the RMSEs of the neutral group are less 4 to 6 • compared to the non-neutral group. From 2012 to 2014, the direction RMSEs of the neutral group are less than 20 • and start increasing from 2015. In 2016, both the neutral and non-neutral groups' RMSEs increase to 30 • and then decrease in 2017, but wind direction RMSEs under non-neutral conditions are still about 7 • higher than the mission requirement.
HSCAT wind products can meet the mission requirement under neutral conditions. Under non-neutral conditions, wind speeds can also satisfy the requirement while wind direction offsets are large. In terms of the variations of statistics, although wind speed biases in 2013, 2015, and 2016 are larger than the other years, the data can still meet the requirement. Wind direction biases and RMSEs increase from 2015 and RMSEs of the non-neutral group in 2016 can get to 30 • . Under non-neutral conditions, variations of the wind direction are larger than the wind speed, which shows that HSCAT wind direction retrieval is more susceptible than the wind speed to sea and atmospheric parameters. Besides, wind speed and direction offsets are large whether for the neutral group or non-neutral group in 2016 and the offset differences are less than other years, which illustrates the worst quality for the 2016 product in an indirect way. After 2015, HSCAT wind products need strict quality control, especially the wind direction products when HSCAT wind data are used in marine studies.

The Rain Impact on HSCAT Wind Products
As mentioned above, the sea surface wind field is retrieved from scatterometer NRCS, which is affected by sea surface roughness. At different rainfall rates, the changes of the sea surface roughness are complicated. Besides, the rain drop size and density also have an impact on the return signal scattering. Therefore, dealing with rain contaminated data is one of the difficulties in scatterometer wind retrieval. It is especially important for scatterometer data usage under rainy conditions to study HSCAT wind field variations at different rainfall rates [20].
We collocated the HSCAT and buoy collocation dataset with SSMIS data within 1 h. We divided rain collocation data into four groups according to different rainfall rates: Rain free (rainfall rate equals to 0 mm/h), light rain (rainfall rate is from 0 to 3 mm/h), moderate rain (rainfall rate is from 3 to 8 mm/h), and heavy rain (rainfall rate is higher than 8 mm/h). To avoid the extra offsets caused by small numbers of rain sample data, we studied all the rain collocated data from 2012 to 2018.

The Distribution of Collocated Data and Statistical Parameters
We compared HSCAT wind data with buoy wind data, and the wind speed and wind direction collocation contours of histograms are shown in Figure 10. Figure 10a,b represent the wind speed and direction collocation contours of histograms; number 1 to 4 represent collocation contours of histograms of the rain free group, light rain group, moderate rain group, and heavy rain group, respectively.
It is evident that the collocation numbers decrease with an increasing rainfall rate and both the wind speed and direction collocation contours appear as losing trends. As for wind speed, the speed contours move to the top of the figure with the increasing rainfall rate, which means HSCAT-retrieved wind speeds are increasingly overestimated. HSCAT wind speeds are almost twice the buoy speeds and have large deviation in heavy rain. As for wind direction, HSCAT wind direction numbers increase in the range 60-120 • and 240-300 • , which means that the HSCAT wind directions tend to be perpendicular to the buoy directions. The two-peak characteristic is evidently shown in Figure 10, which increases the wind direction estimation errors.
We also calculated the statistics of the wind vector biases in different rainfall rates, which are given in Table 4. The statistical parameters also illustrate the decreasing collocation numbers and increase of the wind speed and direction offsets with an increasing rainfall rate. The average biases of the wind speed increased and RMSEs also multiplied with the increasing rainfall rate, i.e., the higher the rainfall rate, the worse the HSCAT-retrieved wind speed product. Besides, the positive biases of the wind speed also increase, which reflects that the overestimation of the retrieved speeds is also increasing with the rainfall rate. In the rain-free group, the wind speed and direction can satisfy the mission requirement.
In the light rain group, the wind speed RMSEs are over 2 m/s, which reflects that the HSCAT-retrieved speed data need quality control even under light rain conditions. Similar to the wind speed, the wind direction data quality greatly decreases with the increasing rainfall rate. The wind direction RMSEs are almost 40 • under heavy rain conditions, i.e., twice the mission requirement. The wind direction is also a little larger than the requirement even in the rain-free group.
affected by sea surface roughness. At different rainfall rates, the changes of the sea surface roughness are complicated. Besides, the rain drop size and density also have an impact on the return signal scattering. Therefore, dealing with rain contaminated data is one of the difficulties in scatterometer wind retrieval. It is especially important for scatterometer data usage under rainy conditions to study HSCAT wind field variations at different rainfall rates [20].
We collocated the HSCAT and buoy collocation dataset with SSMIS data within 1 h. We divided rain collocation data into four groups according to different rainfall rates: Rain free (rainfall rate equals to 0 mm/h), light rain (rainfall rate is from 0 to 3 mm/h), moderate rain (rainfall rate is from 3 to 8 mm/h), and heavy rain (rainfall rate is higher than 8 mm/h). To avoid the extra offsets caused by small numbers of rain sample data, we studied all the rain collocated data from 2012 to 2018.  It is evident that the collocation numbers decrease with an increasing rainfall rate and both the wind speed and direction collocation contours appear as losing trends. As for wind speed, the speed contours move to the top of the figure with the increasing rainfall rate, which means HSCATretrieved wind speeds are increasingly overestimated. HSCAT wind speeds are almost twice the buoy speeds and have large deviation in heavy rain. As for wind direction, HSCAT wind direction Figure 10. Contours of histograms between HSCAT and buoy wind data in different rainfall rate. (a,b) represent the wind speed and direction collocation contours of histograms; number 1 to 4 represent the collocation contours of histograms of the rain free group, light rain group, moderate rain group, and heavy rain group, respectively. Therefore, HSCAT can basically satisfy the marine study requirement under rain-free conditions, although HSCAT wind directions have larger offsets. In rainfall conditions, both the speed and direction offsets are larger than the requirement. When the rainfall rate is higher than 3 mm/h, we do not recommend the use of HSCAT wind products due to their large deviations. It is also a key problem to be resolved in order to improve the accuracy of HSCAT wind retrieval in rainfall conditions.

The Probability Distribution Functions (PDFs) of Wind Speeds and Directions
We plotted the PDFs of wind speeds and directions in different rainfall rates from 2012 to 2018 as shown in Figure 11 with a speed of 1 m/s and a 10 • direction interval. (a) and (b) represent the wind speed and direction PDFs; number 1 to 4 represent the rain-free group, light rain group, moderate rain group, and heavy rain group PDFs; and the dotted line is the PDF for the buoy data and the solid line is for HSCAT PDF, respectively. rainfall rate, the worse the HSCAT-retrieved wind speed product. Besides, the positive biases of the wind speed also increase, which reflects that the overestimation of the retrieved speeds is also increasing with the rainfall rate. In the rain-free group, the wind speed and direction can satisfy the mission requirement. In the light rain group, the wind speed RMSEs are over 2 m/s, which reflects that the HSCAT-retrieved speed data need quality control even under light rain conditions. Similar to the wind speed, the wind direction data quality greatly decreases with the increasing rainfall rate. The wind direction RMSEs are almost 40° under heavy rain conditions, i.e., twice the mission requirement. The wind direction is also a little larger than the requirement even in the rain-free group.
Therefore, HSCAT can basically satisfy the marine study requirement under rain-free conditions, although HSCAT wind directions have larger offsets. In rainfall conditions, both the speed and direction offsets are larger than the requirement. When the rainfall rate is higher than 3 mm/h, we do not recommend the use of HSCAT wind products due to their large deviations. It is also a key problem to be resolved in order to improve the accuracy of HSCAT wind retrieval in rainfall conditions.

The Probability Distribution Functions (PDFs) of Wind Speeds and Directions
We plotted the PDFs of wind speeds and directions in different rainfall rates from 2012 to 2018 as shown in Figure 11 with a speed of 1 m/s and a 10° direction interval. (a) and (b) represent the wind speed and direction PDFs; number 1 to 4 represent the rain-free group, light rain group, moderate rain group, and heavy rain group PDFs; and the dotted line is the PDF for the buoy data and the solid line is for HSCAT PDF, respectively. (a4) (b4) Figure 11. Probability distribution functions (PDFs) between HSCAT and buoys under different rainfall rates. (a) and (b) represent wind speed and direction PDFs; number 1 to 4 represent the rainfree group, light rain group, moderate rain group, and heavy rain group PDFs; dotted line is the PDF for buoy data and solid line for HSCAT PDF, respectively.
The wind speed and direction PDFs reconfirm the results in Section 3.4.1.
For wind speed, the distributions of the buoy speed are mainly distributed in 6 to 8 m/s and are nearly unchanged with the increasing rainfall rate. Therefore, the buoy speed PDFs' fluctuation can be regarded as changes of the ocean environment and the decreasing of collocations with the increasing rainfall rate. However, the entire distribution of the buoy speed is stable so that the buoy speed can be a reliable evaluation standard. However, the peaks of the HSCAT wind speed PDFs move from 5 to 6 m/s under the rain free condition to 16 to 18 m/s under the heavy rain condition. Figure 11. Probability distribution functions (PDFs) between HSCAT and buoys under different rainfall rates. (a,b) represent wind speed and direction PDFs; number 1 to 4 represent the rain-free group, light rain group, moderate rain group, and heavy rain group PDFs; dotted line is the PDF for buoy data and solid line for HSCAT PDF, respectively.
The wind speed and direction PDFs reconfirm the results in Section 3.4.1. For wind speed, the distributions of the buoy speed are mainly distributed in 6 to 8 m/s and are nearly unchanged with the increasing rainfall rate. Therefore, the buoy speed PDFs' fluctuation can be regarded as changes of the ocean environment and the decreasing of collocations with the increasing rainfall rate. However, the entire distribution of the buoy speed is stable so that the buoy speed can be a reliable evaluation standard. However, the peaks of the HSCAT wind speed PDFs move from 5 to 6 m/s under the rain free condition to 16 to 18 m/s under the heavy rain condition. The overestimation of the HSCAT speed increases as the rainfall rate increases, and the HSCAT speeds of the heavy rain group are almost twice the buoy speeds.
For wind direction, the rain free buoy wind directions are mainly distributed in the range 60-120 • and then slightly move to 180 • with an increasing rainfall rate. However, the movement is quite small so the buoy direction can be seen as a relatively stable evaluation standard. Similar with the results in Section 3.4.1, under the rain-free condition, the HSCAT and buoy direction PDFs are consistent well. However, with the rainfall rate increasing, the HSCAT wind directions have two peaks in the ranges 60-120 • and 240-300 • with growing peak values, and the extent of the two PDFs' separation increases. It illustrates that HSCAT-retrieved wind directions tend to be perpendicular to the true directions with an increasing rainfall rate, which causes large wind direction offsets under rainy conditions.

Conclusions
The stabilities of the HY-2A scatterometer in different times, atmospheric stratification conditions, and rainfall rates were analyzed by comparing HSCAT L2B wind products with buoy wind measurements from 2012 to 2018. In the study, HSCAT and buoy data were selected within the range 0.1-25 m/s and collocated within 25/

√
2 km (about 17.7 km) and 30 min. After data quality control, yearly matching data were 34,462,31,276,27,339,36,211,25,898,27,511, and 24,997 pairs from 2012 to 2018. According to the collocated data, wind speed biases and RMSEs between HSCAT and buoy measurements were about 0.1 to 0.2 m/s and 1.2 m/s while biases and RMSEs of the wind direction were 1 • and 27 • , respectively. These suggest that HSCAT wind direction accuracy is slightly worse than the scatterometer mission requirement (20 • ). The yearly PDF curves also supported these results. In 2015 and 2016, the HSCAT PDFs and buoy PDFs showed clear separation and then tended to agree with each other in 2017. From 2017, HSCAT winds and buoy winds had similar distributions and the offsets tended to decrease.
For wind speed and direction bias variations in different speed bins, the collocated data were mostly distributed in the speed range 6-8 m/s. The minimums of the wind speed bias and RMSE were distributed around 3 to 11 m/s while the minimums of the wind direction biases and RMSEs appeared when speeds are higher than 5 m/s. When speeds are higher than 10 m/s, the average speed biases and RMSEs were increased with increasing binning speeds. The wind direction offsets were barely changed, except the RMSEs of the wind direction increased slightly when speeds were higher than 16 m/s. According to the HSCAT wind field data variations over time, wind speed and direction average biases were distributed around 0. Wind speed RMSEs were about 2 m/s and wind direction RMSEs were about 20 to 30 • . HSCAT wind speed biases were a little larger in 2013, and HSCAT wind products were the most unstable in 2015 and 2016, which had the largest offsets and variations among the seven years.
To study the atmospheric stratification impacts on HSCAT wind products, we classified the collocated data into a neutral group (the absolute values of temperature difference are less than 1 • C) and non-neutral group (the absolute values of temperature difference are higher than 1 • C) according to the differences between the sea surface temperature and air temperature and then analyzed the wind field error features in the different groups. The wind PDFs and error characteristic variations with the wind speeds illustrated that the quality of HSCAT-retrieved wind products are better in neutral conditions. Wind direction has more of an influence from marine and atmosphere environmental parameters and owns large changes in the non-neutral condition. From the yearly statistical parameters of the collocated data, the HSCAT wind product qualities declined from 2015 but were improved from 2017. In 2016, the offset differences between the neutral and non-neutral group were less than other years, which reflects the worst quality in an indirect way and the data in 2016 need to be treated very carefully in applications.
We also studied the rain impacts on HSCAT wind products. We divided rain-collocated data into a rain-free group (0 mm/h), light rain group (0-3 mm/h), moderate rain group (3-8 mm/h), and heavy rain group (>8 mm/h) according to the rainfall rate. From the wind PDFs and collocated contours, we found that the accuracy of the HSCAT-retrieved wind field decreases as the rainfall rate increases: HSCAT-retrieved wind speeds and the RMSEs of collocated speeds increased with an increasing rainfall rate; HSCAT-retrieved wind directions tended to concentrate on the 60-120 • and 240-300 • ranges, which means HSCAT-retrieved wind directions tend to be perpendicular to the true directions. These results may provide references for the correction of scatterometer rain-contaminated data.
Overall, HSCAT L2B wind products are less accurate and less stable than the mission requirement especially for wind direction products during 2012 to 2018. The qualities of HSCAT operational wind products declined from 2012 to 2016 but were improved from 2017. HSCAT wind products should be treated carefully, especially in 2015 and 2016, or adopt other ways to improve data qualities when HSCAT wind data are used in marine studies.
Long-term evaluation for HSCAT wind products was analyzed in this study, and we also studied the error characteristics in different times, atmospheric stratification conditions, and rainfall rates primarily. However, only a little collocated data was distributed in the high wind speed region due to the lack of high speed data from buoys. The reasons for quality improvement from 2017 were not discussed here. Comparing HSCAT wind data with other satellite scatterometer data, studying the environmental parameter impacts on scatterometers working in different bands, finding improvement methods for HSCAT winds, and other related studies need to be investigated in the future.