ICONE.GIF (1472 octets) .../...

4.2 Comparison with buoy data

The aim of this section is to estimate the accuracy of the weekly and monthly wind speed and direction in comparison with buoy wind data. This is achieved by using : the National Data Buoy Center (NDBC), the Tropical Atmosphere Ocean (TAO), and the European Buoys (ODAS) buoy networks (Figure 9). More than 90 buoys covering Atlantic and Pacific ocean areas betwenn 10°S and 57°N.

424.GIF (14540 octets)

Figure 9 : Buoy network location

For the validation of the scatterometer average wind field, the buoy wind data are referenced to 10m height, assuming a logarithm wind profile, Von Karman's constant of 0.4, neutral stratification and, a wind speed dependent drag coefficient (Ezraty 1987).

For each week and each month, mean values of buoy wind speed, zonal and meridional components are computed arithmetically. Weekly and monthly means are computed for all ERS-1, ERS-2 and NSCAT periods for which at least 3.5 days and 15 days buoy measurements are collected, respectively. For each averaging period, the closest scatterometer grid point (1° ´ 1° ) to each buoy location is selected. Therefore, a collocated data sets between scatterometer gridded wind fields (averaging objective method) and buoy averaged winds are performed for NDBC, TAO and ODAS buoy networks. Results are then compared using the following standard statistic data analysis :

The wind speed, zonal component and merdional component are assumed as a random variables wich could be characterized by their moments. For this purpose, the four conventional (C moments) and linear moments (L moments) of each variable are estimated.

Let is W a wind variable (wind speed, zonal component, merdional component or wind difference). The corresponding four C moments are determined as :

                   F406.GIF (2618 octets)            (4)

SYM46.GIF (875 octets), SYM19.GIF (60 octets)W, SW and KW are the W mean (bias), standart deviation, skweness and kurtosis, respectively. Variance and rms values are dereived from SYM46.GIF (875 octets) and SYM19.GIF (60 octets)W estimates.

The L moments (Hosking, 1990} are defined by :

   F407.GIF (1245 octets)            (5)

SYM12.GIF (64 octets) n is the nth linear moment of W

F408.gif (907 octets) is the shifted nth Legendre polynomial. It is related to Legendre polynomial F413.GIF (896 octets) by :

F409.GIF (1146 octets)                      (6)

F is the probability function of wind variable W

Q(F), called quantile function, is provided by the following equation :

F410.GIF (1393 octets)       (7)

The meaning of C moments and L moment are similar as can be shown through the equations. The main advantage of L moments is their relative small sensitivity to data errors generaly producing outliers in data series.

The statistical significance of the fisrt and second moment is evaluated by Student test (T-test) and Fisher test (F-Test), respectivelly. Troughout this paper, the significance is estimated for 95% confidence.

Moreever, the linear regression parameters are estimated to assess the comparisons between satellite gridded wind fields and buoy averaged winds. In this paper we provide the following parameters :

F411.GIF (2973 octets)                  (8)

Where  F412.GIF (1166 octets)

x and y denote the buoy and scatterometer wind estimates, respectively. b is the slope and a is the intercept on the y axis : y = bx + a. bs is the slope of symmetric regression line. SYM18.GIF (62 octets) is the correlation coefficient. Its calculation involves the residual, SYM05.GIF (58 octets), between y and linear regression model. SYM19.GIF (60 octets) p1, and SYM19.GIF (60 octets) p2 are the rms deviations of the first and second principal component of x and y distribution. They provide a measurement of the major and minor axis of the elliptical x and y distribution.

 

4.3 Global comparisons

Table 2, 3, and 4 provide the main statistical parameters characterizing wind speed comparisons. The wind speed correlation coefficients ranging from 0.85 to 0.89 indicate a good consistency between satellite and buoy averaged winds. The rms values of the differences buoy-satellite wind speeds do not exceed 1.16m/s over NDBC and TAO networks. Results derived from ODAS/satellite comparisons show higher rms values : 1.48m/s for NSCAT, and 1.66m/s for ERS-2. The latter are mainly due to a poor number of comparison data points, and to the high wind variability in ODAS area (Figure 9). Furthermore, the statistics calculated by several meteorological centers (ECMWF, CMM, UKMet) indicate that ODAS buoy wind speed tend to be underestimated according to meteorological wind analysis (see ftp://ftp.shom.fr/meteo/qc-stats, site maintained by P. Blouch).

The results of the regression analyses carried out on collocated data, show that the slopes calculated over each buoy network and against buoy wind estimates, are quite similar for the three averaged scatterometer wind speeds. In NDBC area (Table 2), buoy and scatterometer wind speeds agree quite closely, which is expressed by slopes of about 1 and intercepts of about zero. Comparisons between buoy and scatterometer winds in Pacific tropical ocean give regression line slopes of about 0.80, suggesting an overestimation of low wind speed and underestimation of high wind speed by scatterometer wind fields compared to TAO winds. In north Atlantic area, the slopes are very close to 1, whereas the intercepts are of about 0.50, indicating that the scatterometer wind fields are consistently high compared to ODAS week-averaged wind speeds. The calculation of the statistical parameters according to the buoy wind speed ranges, show that their values are made variable by the outlying points at low and high wind speeds.

For the wind direction, no systematic bias is found, and the overall bias and standard deviation about the mean angular difference are less 8° and 38° , respectively. These results are consistent with the calibration/validation of the scatterometers against buoy (Graber et al, 1996 and 1997; Caruso et al, 1999). For instance, in Pacific tropical area, where the wind direction is quite steady, the standard deviation calculated for buoy wind speed higher than 5m/s, does not exceed 17° .

Table 2 : Comparison of averaged weekly wind speed and direction estimated from NDBC buoy measurements and from ERS-1, ERS-2 and NSCAT scatterometer observations.

Data SET

BuoyWind Speed Range (m/s)

Length

Wind Speed (m/s)

Wind Direction

Bias
(m/s)

Rms
(m/s)

r

b

a

bs

s p1

s p2

Bias
(deg)

Std
(deg)

NDBC
/ ERS-1

0-24

3281

0.02

1.16

0.88

0.99

0.00

1.16

2.87

0.78

3

35

 

0-5

320

-0.14

1.03

0.74

0.87

0.68

2.12

1.14

0.47

5

47

 

5-10

2603

0.05

1.16

0.83

1.01

-0.14

1.35

2.04

0.72

3

34

 

> 10

358

-0.0

1.31

0.76

0.97

0.32

1.80

1.64

0.69

3

30

NDBC
/
ERS-2

0-24

1921

0.35

1.15

0.89

0.96

-0.07

1.12

2.76

0.75

6

33

 

0-5

142

0.06

0.82

0.75

0.87

0.50

1.85

0.96

0.42

0

47

 

5-10

1581

0.37

1.16

0.83

0.98

-0.23

1.30

1.97

0.71

6

33

 

> 10

198

0.40

1.26

0.77

0.82

1.61

1.42

1.60

0.75

6

25

NDBC
/NSCAT

0-24

522

-0.38

1.02

0.90

0.96

0.68

1.09

2.58

0.65

8

25

 

0-5

28

-0.54

0.94

0.76

1.08

0.17

1.95

0.95

0.37

3

29

 

5-10

444

-0.37

1.01

0.85

0.96

0.69

1.21

1.87

0.62

8

26

 

> 10

50

-0.32

1.15

0.79

0.78

2.68

1.24

1.62

0.74

7

15

 

Table 3 : Comparison of averaged weekly wind speed and direction estimated from TAO buoy measurements and from ERS-1, ERS-2 and NSCAT scatterometer observations.

Data SET

BuoyWind Speed Range
(m/s)

Length

Wind Speed (m/s)

Wind Direction

Bias
(m/s)

Rms
(m/s)

r

b

a

bs

s p1

s p2

Bias
(deg)

Std
(deg)

TAO /
ERS-1

0-24

10047

0.29

0.89

0.89

0.80

0.85

0.94

2.15

0.59

3

31

 

0-5

3262

-0.14

0.85

0.76

0.70

1.31

1.32

1.06

0.54

1

51

 

5-10

6693

0.47

0.91

0.84

0.86

0.42

1.12

1.51

0.54

5

17

 

> 10

92

0.24

0.92

0.70

0.24

7.66

2.69

0.86

0.31

8

9

TAO /
ERS-2

0-24

6737

0.56

1.03

0.89

0.80

0.63

0.93

2.26

0.60

3

27

 

0-5

1925

0.06

0.84

0.75

0.67

1.22

1.35

1.01

0.54

4

45

 

5-10

4736

0.75

1.10

0.85

0.87

0.12

1.11

1.60

0.55

5

16

 

> 10

76

0.76

1.14

0.78

1.74

-8.48

2.94

1.02

0.26

7

10

TAO /
NSCAT

0-24

1780

-0.26

0.92

0.88

0.80

1.47

0.94

2.20

0.62

5

20

 

0-5

515

-0.70

1.18

0.74

0.71

1.81

1.61

1.10

0.55

2

33

 

5-10

1246

-0.08

0.79

0.83

0.80

1.39

1.07

1.47

0.55

7

11

 

> 10

19

0.03

0.82

0.78

1.85

-8.92

3.09

0.99

0.24

10

5

 

Table 4 : Comparison of averaged weekly wind speed and direction estimated from ODAS buoy measurements and from ERS-2 and NSCAT scatterometer observations

Data SET

BuoyWind Speed Range
(m/s)

Length

Wind Speed (m/s)

Wind Direction

Bias
(m/s)

Rms
(m/s)

r

b

a

bs

s p1

s p2

Bias
(deg)

Std
(deg)

ODAS
/ ERS-2

0-24

222

-0.70

1.66

0.88

1.02

0.51

1.20

3.58

0.99

1

38

 

0-5

10

-1.26

2.01

0.72

0.63

2.65

1.96

1.65

0.76

31

75

 

5-10

155

-0.61

1.68

0.80

1.18

-0.76

1.73

2.31

0.82

3

39

 

> 10

57

-0.83

1.50

0.80

0.78

3.31

1.20

1.91

0.85

4

22

ODAS
/ NSCAT

0-24

194

-0.63

1.48

0.91

1.00

0.55

1.13

3.82

0.91

2

30

 

0-5

6

-1.29

2.07

0.72

0.51

3.26

1.96

1.65

0.79

14

76

 

5-10

118

-0.62

1.44

0.81

1.12

-0.37

1.58

2.05

0.73

1

30

 

> 10

70

-0.57

1.47

0.86

1.11

-0.84

1.37

2.73

0.82

9

22

 

A413.GIF (12194 octets)

B414.GIF (17765 octets)

C415.GIF (17042 octets)

Figure 10 :

a/ Scatter plot and frequency of wind speeds, x-axis shows averaged TAO wind speed,
    y-axis shows the gridded scatterometer wind speed
b/ as Figure 12a but for zonal component
c/ as Figure 12a but for meridional component

The geographical features of the difference between scatterometer and buoy weekly wind estimates have been investigated. At each TAO buoy location, statistical parameters of the difference series are computed. For instance the rms values of the wind speed and direction differences are shown in Figure 11. The main result is that the rms values are generally higher within the 150E and 170E Pacific band. A study of the correlation between in-situ and sensor matchups indicates that there is a significant difference between the correlation coefficients computed in the eastern Pacific and those computed in the western Pacific with 95 % confidence. Using TAO buoy measurements, it was determined that in the western Pacific there is 6 times more energy compared to the eastern Pacific (Mangum, 1992) . This high variability of the wind is mainly explained by the high surface temperature and the convective activity of this zone. Averaging procedures give less representative results in these wind conditions.

a)

9 N       1.00                  
8 N   1.04           1.13          
5 N     0.87 0.81 0.92 0.89   1.09 0.95   1.07    
2 N 0.92 1.09 1.03 0.85 0.86 0.85   1.19         1.16
0 N   0.86 0.75 0.87 0.80 0.76 0.73 0.88 0.92 0.84   1.02  
2 S   0.97 0.83 0.86 0.77 0.89   1.16 1.06        
5 S   0.77 0.68 0.80 0.89                
8 S         0.89     1.03          
  95W 110W 125W 140W 155W 170W 170E 165E 156E 154E 147E 143E 137E

b)

9 N       36                  
8 N   33           40          
5 N     24 17 26 31   54 39   54    
2 N 34 33 19 14 22 28   62         76
0 N   26 28 22 19 33 52 62 63 73   63  
2 S   27 21 23 18 25   62 55        
5 S   20 19 22 28                
8 S         23     42          
  95W 110W 125W 140W 155W 170W 170E 165E 156E 154E 147E 143E 137E
Figure 11 :

a/ Geographical behavior of the rms difference between gridded scatterometer and averaged TAO wind speeds

b/ Geographical behavior of the rms difference between gridded scatterometer and averaged TAO wind directions

[R]


.../... ICONE.GIF (1472 octets)