Global Wave Climatology - Ifremer

18 downloads 0 Views 458KB Size Report
measurement artifacts of different nature when very high waves are in question. ... last century and to further attribute them to the wind storms, potentially freak ...
Extreme waves in visual wave observations by VOS Vika Grigorieva and Sergey K. Gulev P.P. Shirshov Institute of Oceanology, Russian Academy of Sciences, Moscow, Russia [email protected] http://www.sail.msk.ru Being a disastrous phenomenon, extreme waves (including those classified as rouge waves) are critically important for marine structures. Although there are many works on the theory of rogue waves [e.g. 1, 2, 3, 10, 11, 12, 13], their physics of is still poorly understood and they are hardly detectable in different wind wave data. Furthermore, we have very tentative estimates of rogue waves and in most wave data it is hard to discriminate between the observational truth and the measurement artifacts of different nature when very high waves are in question. In this paper we will analyse the collection of visual observations worldwide in order to quantify the most extreme wave cases and to identify their nature. Besides model hindcasts, satellite and buoy measurements Voluntary Observing Ship (VOS) data represent now an important source of global wave information, including that of extreme wave events. There is a general concern that VOS observations are quite inaccurate and, therefore, are hardly applicable for quantitative wind wave analysis. However this concern results primarily from a poor knowledge of the actual errors and uncertainties inherent in the VOS wave observations. This accuracy was more or less quantified during the last years and many errors and uncertainties were either minimized or removed from the VOS collection [5]. Visual wave observations were extensively used for the description of climatological characteristics of wind waves [5, 6, 7] and for the assessment of the long term tendencies in wave parameters [8, 9]. However, they were not yet employed for the analysis of freak waves. We used these data to quantify the highest observed wave cases in the World Ocean during the last century and to further attribute them to the wind storms, potentially freak waves, or observational artifacts. Voluntary observing ships provide massive visual wave information worldwide covering the period from 1885 onwards [14, 15]. Before 1950 marine officers reported the highest wave component (sea or swell) and its upper value was limited by 16 m according to the existing coding system. From 1950s observers started to report separately wind sea height and swell height estimates. Furthermore, the coding system was changed and allowed for the detection of the maximum wave height up to 25 meters for both sea and swell. From our collection of more than several million reports we selected all cases reporting wave height exceeding 15 meters. Altogether, after 1950 there has been found more than 1000 of such cases for seas and about 3000 cases for swells (Fig.1). Importantly, we considered separately the subsets of reports providing information about all wave parameters (wind sea and swell heights, their periods and directions) and of incomplete reports which provide typically only one wave parameter (e.g. wind sea height). In the first instance we consider wave statistics based on all selected reported without any census. -100

-80

-60

-40

-20

0

20

40

60

80

100

120

140

160

-100

180 -160 200 -140 220 240 -120 -100 260 -80 280

80

80

80

60

60

60

40

40

wind sea january

40 20 0 -100

-80

-60

-40

-20

0

20

40

60

80

100

120

140

160

180

-160 220 -140 200

16 to 1820to 20 to 220 to -120 260 -100 280 -80 240

18 2020 22 0 25-100

-80

-60

-40

-20

0

20

40

60

80

100

120

140

160

180 -160 200 -140 220 240 -120 -100 260 -80 280 80 60

swell january

40 20

-80

-60

-40

-20

0

20

40

60

80

100

120

140

160

180

-160 220 -140 240 -120 260 -100 280 -80 200

0

Fig.1 Extreme wind sea (left panel) and extreme swell (right panel), m, Northern Hemisphere, January, 1950-2006

Figure 2 shows the occurrence histograms of the extreme wind wave and swell heights, their periods, wind speed for all selected extreme cases (i.e. reported wave height is higher than 15 m). Additionally we also estimated the occurrences of wave periods, steepness and wave lengths. Furthermore, wave parameters were analyzed within the radius of 100 km around the location of extreme wave report.

0.5

0.3

0.5

0.5

wind sea height

0.4

corresponding wind

0.4

0.5

0.5

max swell height

wave length

0.4

corresponding wind

0.4

wave length

0.4

0.2

0.3

0.3

0.3

0.2

0.2

0.2

0.1

0.1

0.1

0

0 15 16

17 18 19 20 21 22 23

24 (m)

15

20

25 30

35

40

0

0.4

17 18 19 20 21 22 23

0.3

0.2

0.1 0 5 10

15

20

25 30

35

40

0.5

0.3

average swell

swell periods

steepness*100

0.4

0.2

0.2

0.3

0.2

40 80 120 160 200 240 280 320 400 (m)

50 (m/s)

Max swell 1950-2006 january Northern Hemisphere 2324 reports

average wind sea

wind sea periods 0.3

24 (m)

0.3

0.4

steepness*100

0.2

0 15 16

Max wind sea 1950-2006 january Northern Hemisphere 667 reports 0.3

0.2 0.1

40 80 120 160 200 240 280 320 400 (m)

50 (m/s)

0.3

0.1

0 5 10

0.3

0.2 0.2 0.1

0.1 0.1

0

4

8

12

16 20 24 28 32

36

40

0.1

0.1

0

0 2

4

6

8

10 12 14

16

18 (s)

0.1 0 1

2

3

4

5

6

7

8

10 (m)

4

8

12

16 20 24 28 32

36

40

0

0 2

4

6

8

10 12 14

16

18 (s)

1

2

3

4

5

6

7

8

10 (m)

Fig.2 Occurrence histograms for the wind sea (left panel) and swell (right panel), characteristics over the Northern Hemisphere for January (1950-2006).

Preliminary results derived from the subset of all selected extreme cases show quite reliable probability distributions for the whole range of wave parameters implying that the distribution of extreme wave characteristics should now represent the tails of these distributions. Maximum wind speed changes from 5 to 10 m/s, with the wave lengths being about 40 meters for the wind sea and 80-160 m for swell (70% of total subset). Corresponding wave periods change in a wide range from 1 to 25 seconds with the peak of probability at 5 seconds for seas and at 8-10 seconds for swell. We have to note, however, that wave periods are considered to be the most uncertain parameters in VOS and their errors can reach of 30 % from average values. Two-dimensional occurrence distributions of extreme wave cases and corresponding wind speed are presented in Figure 3. They demonstrate the bimodal structure implying 2 maxima of the occurrences at wave heights of 15-16 m and 20-21 m. Multiply analysis of coding errors nevertheless confirmed the existence of this structure. Very tentatively we can hypothesize that the first maximum is associated with the high waves originating from wind storms and the second one results from the observational artifacts or the freak waves. Analysis of the wave steepness derived from the estimates of wave heights and periods allows for the discrimination of some observational artifacts. Usually wave steepness is defined as δ = Hw /λ = 2*π* Hw/(g*Pw2), where Hw is wind sea height, λ is wave length, Pw is wind sea period, g is the gravitational acceleration. According to the wave theory δ > 1/7 results in increasing wave steepness and the loss of stability, forcing wave breaking. In the case of multi-system sea, however, wave height and wave length are statistical parameters computed from the wave spectra. This implies some uncertainty in this definition. For instance for some extreme cases identified in the VOS data steepness can reach unrealistic values, particularly due to uncertainties in the estimation of wave periods. For example, for Pw= 1 sec δ ≈ 0.64Hw. Therefore, we consider this criterion to be effective enough and adopted it for the further analysis.

50

50

50

50

wind sea Jan

swell Jan 40

40

40

40

5

W I N D (m/S)

30 2.5

1.5

20

0.5

10

10

15

20

10

10

25

15

H E I G H T S (m)

20

25

H E I G H T S (m)

50

50

50

50

wind sea Jul

swell Jul

40

40

40

40 9

W I N D (m/S)

5

30

2.5

20

1.5

0.5

10

10

15

20

10

25

10

15

20

25

Fig.3 Two-dimensional probability distributions of the wind sea heights and wind speed for January and July over the Northern Hemisphere 0.4

Hmax/Hmean 0.3

0.2

0.1

0 1

2

3

4

5

6

7

8

9 10

Fig. 4 Histogram of the ratio η for the Northern Hemisphere, January 1950-2006 One of the key parameters used for the detection and quantification of the rogue waves is the ratio between maximum and the averaged wave height in the record. We used this criterion, taking the average value from all observations within the radius of 100 kilometers from location of maximum. In this case we consider the ratio η=Hmax/Hmean where Hmax is the reported extreme wave height and Hmean is the average over all simultaneous reports within 100 km radius. Occurrence histograms for these mean estimates are presented in Figure 4. Usually the threshold of η=2 is considered for the detection of extreme cases, the values of η>4-5 are

identified very seldom. The occurrence histogram of η (Fig.4) show that the ratio varies from 1 to 25 with the peak values of 8-9. Thus, this ratio is not effective enough for the detection of the observational artifacts and we need another criterion to check our data quality. In the next step we performed an accurate quality control of the selected extreme wave cases. First, we selected only the reports providing information about all wave parameters (wind sea and swell heights and periods) and excluded incomplete reports which provide frequently only one wave parameter. This census resulted in the reduction of the number of reports by about 30%. Then we excluded the reports without wind speed and reports which give suspicious characteristics of the wave steepness. We used the threshold of 0.2 for wind sea steepness and of 0.15 for swell steepness. This limitation broke up to 80% reports. Finally 129 extreme wind sea cases and about 800 extreme swell cases were left for the further analysis. The frequency distributions derived from this subset show already quite significant differences from those derived from the whole census (Fig.5). For the wind waves they show one major maximum at 15-16 m with the smaller secondary one at more than 20 meters. Remarkably the wind speed distribution is represented by the wide plato with the modal values varying within 535 m/sec. However, when we consider the extreme swells, the two maxima in the distributions remain, being still centered at the same values (15-16 and 20-21 meters). 0.6

0.3

0.3

max wind sea

wind sea steepness*100

corresponding wind

0.4

0.2

0.2

0.2

0.1

0.1

0

0 15 16

17 18 19 20 21 22 23

24 (m)

0.5

0 5 10

15

20

25 30

35

40

4

6

8

10

12 14 16 18 20

0.3

0.3

swell steepness*100

corresponding wind

max swell height

0.4

2

45 (m/s)

0.2

0.2

0.1

0.1

0.3 0.2 0.1 0

0 15 16

17 18 19 20 21 22 23

24 (m)

0 5 10

15

20

25 30

35

40

45 (m/s)

2

4

6

8

10

12 14 16 18 20

Fig.5 Occurrence histograms after the quality control for the wind sea (upper pictures) and swell (down pictures) parameters, Northern Hemisphere, January 1950-2006 In the next instance we analyzed individual cases in order to attribute extremes to the storm conditions, artifacts or freak waves. For the seas we requested that storm wind speed for all cases exceeds at least 10 m/sec and that the sea level pressure should be deeper than 1002 hPa. These requirements eliminated about 98% of all reports which can be with a high probability associated with severe storms. After that the rest consisted of few cases (e.g. 19 cases for January). They are located in the Northern hemisphere mid latitudes and few in the Atlantic tropics where then are likely associated with the tropical cyclone development. All these cases report wind waves higher than 16 meters with the highest waves (of 22 meters and more) identified in the eastern mid-latitudinal Atlantic and northeast Pacific. The data from these reports were further analyzed with respect to the wave age [4]: a = Cp / Vef, where Cp is deep water phase speed at spectral peak, which is derived from the peak wave period Pw as Cp = (g/2π) Pw .

Here g is the gravitational acceleration, Vef = V10*Cos θ is the component of the wind in the wave direction, θ is the angle between wave and wind directions (in 95% cases those angle are equal), V10 is the wind at 10 m anemometer height and neutral stability. For all extreme seas wave age a > 1.2. So, those reports are either observational artifacts or swell, which was reported as wind sea. Furthermore, we performed a detailed synoptic analysis of all cases with extreme wind waves using SLP and wind speed data from the NCEP/NCAR re-analysis. Reanalysed winds and SLP also do not hint on the extremes in the VOS data. Thus, extreme heights which passed the quality control of wind waves can not be attributed to the freak waves. Similar analysis of extreme swells is much more complicated, since swell does not depend on the local wind and, thus, local wind cannot be used as additional information for the analysis. We first limited our analysis to the cases when extreme swell is associated with wind sea height smaller than 5 meter. This is reasonable, because high wind seas of 10-15 meters are unlikely to be associated with simultaneously high swells. In other words, sailors cannot accurately identify these cases and discriminate high wind seas and high swells. Secondly, we eliminated from the analysis all cases of extreme swells identified in the marginal seas or in the regions where the fetch is short enough to provide swell growth up to extreme values. Next, for every report of extreme swell we analyzed visual swell reports in the vicinity of about 200 km from this point. We required that the neighboring swell reports should imply swells of at least higher than monthly mean values. All these procedures filtered out substantial amount of reports and left for the further analysis about 1% cases of extreme swells. In the Northern Hemisphere about 50% of extreme swell cases were identified in the North Atlantic and North Pacific mid latitudes and subtropics where they can be attributed to the propagation of the swell parcels from storms. Tropical cases can be associated with the propagating swells form the tropical cyclone systems. For the Southern Hemisphere the number of extreme wind sea and swell cases is generally smaller than for the Northern Hemisphere (Fig.6). This can be reasonably explained by generally smaller number of observations here compared to the Northern Hemisphere. Altogether there have been identified 86 cases of extreme wind seas of which the quality control procedures left 9 cases to be further considered with confidence. For swell data we first selected 227 extreme cases of which only 7 cases were left after the critical analysis. -100 0

-80

-60

-40

-20

0

20

40

60

80

100

120

140

160

-140 240 -120 260 -100 280 -80 180 -160 200 220

0

-100 0

-80

-60

-40

-20

0

20

40

60

80

100

120

140

160

180 -160 200 220 -140 240 -120 260 -100 280 -80

0

-20

-20

-20

-20

-40

-40

-40

-40

-60 -80 -100

max seas 1950-2006, Aug,m -80

-60

-40

-20

0

20

40

60

80

100

120

140

160

180

15 to -60 18 to 20-80to 22 to

-160 220 -140 240 -120 260 -100 280 -80 200

18 -60 20 22 -80 25

-100

-60 -80

max swell 1950-2006, Aug ,m -80

-60

-40

-20

0

20

40

60

80

100

120

140

160

180

-160 220 -140 240 -120 260 -100 280 -80 200

Fig.6 Extreme wind sea (left panel) and extreme swell (right panel), m, Southern Hemisphere, January, 1950-2006 The further analysis of wave age and synoptic situations revealed that all extreme wind seas in the South Ocean can be attributed to either severe storm conditions or observational errors. None of the remained 7 extreme swell reports was found to be suspicious. This work was focused on the explanation of extreme waves in the VOS-collection. These waves do not influence the global climatology, but may play an important role in the extreme wave statistics. We tried to find signatures of the known freak events in VOS data and to show that not all 25 meter heights are the observational artifacts. The VOS collection does not include any reports to be verified with such events as the New Year Wave and the extreme wave in the south-east Africa region. Nevertheless, one report was found near the location of the catastrophe

on 15 February 1982 near the Grand Banks of Newfoundland (Lawton G.) This report fits to all our requirements. Remarkably, it reports swell with the height of 15 meters, period of 15 sec and the length of 351 m implying steepness of 0.04. The corresponding wind sea height was 5 m with the period being 6 sec. We also analyzed wave reports over all well-known banks with specific topography and wave regime. Only the data from the Dogger Bank and Chatham rise (Fig. 7) have been accepted for analysis after the quality control. We have to stress, however, that the period reported (about 17 seconds), should be considered with caution. Although the analysis of synoptic situation shows strong gradients of sea level pressure in this location, we do believe that wind sea of 15 meters height and 9 sec period can be considered as a reliable indication of the storm in this area. 160 -30

170

180

190 -30

Dogger bank

H=16m P=9s V=9.8 m/s

H=20m P=10s V=3 m/s

24.02.1974 -40

-40

H=16m P=10s V=13m/s

-5015 to 18 to 160

19.08.1976

18 20 20 to 22 22 to 25 160

-30

170

180

-50 190

170

180

190 -30

P=10s V=3 m/s -4016.05.1977 06.02.1971

170

H=18.5m P=10s V=10m/s

H=15m P=17s V=7m/s

a>3 !

09.03.1975 -40 01.11.1978

-50 160

H=20.5m P=10s

180

-50 190

a>5 ! H=18.5m P=17s V=3 m/s

H=15m P=8s V=9.9 m/s

Fig. 7 Left panel: Extreme cases on the Chatam rise. Wind sea data are shown on the top and swell data are in the bottom panel. Right panel: extreme cases on the Dogger-bank. Blue color indicates wind sea, red color stands for the swell data.

Considering Chatham rise, we can conclude that extreme wind sea estimates are quite reliable for this region, although extreme swells are likely overestimated. However, we can argue that South Ocean provides conditions for long enough fetch and, therefore, hypothesize that extreme swell can also contribute to the freak events. Our analysis of extreme wave heights in the VOS collection over the last several decades shows that VOS data when properly quality controlled and pre-processed can be capable of capturing even single extreme wave events. Among all extreme wind waves in the VOS collection about 20% can be attributed to wind wave storms and more than 80% are considered as observational artifacts. Frequency distributions of extreme wave heights clearly show maxima peaked at 15-16 meters and 20-21 meters and these maxima hold even after very strict quality control, indirectly hinting on the existence of freak events in the VOS records. At the same time we have to stress that all data reporting extreme swells are highly doubtful. However, few cases (about 1-2% of identified events) can really be considered as confirmation of extremely high swells, potentially of the freak nature.

References

1. Annenkov S., Badulin S. Multi-wave resonance and Formation of High-Amplitude Waves in the Ocean. Proceedings of a Workshop “Rogue waves”. Brest, France, 2000, p.205-214. 2. Boukhanovsky A.V., Lopatoukhin L.I., Rozhkov V.A. Approaches, experience and some results of wind wave climate investigations. III. Exterme and freak waves. (in Russian) Saint-Petersburg University, 2006, volume 7,issue 1, pp. 58-70 3. Davidan I.N., Lopatoukhin L.J., Wind waves in the World ocean. Saint Petersburg (Russia): Hydromet. Publishing; 1985. 278. [In Russian]. 4. Dobson, F.W., S.D.Smith, and R.J.Anderson, Measuring the relationship between wind stress and sea state in the open ocean in the Presence of swell. Atmosphere-Ocean, 32, 61-81, 1994. 5. Gulev, S.K., V.Grigorieva, A.Sterl, and D. Woolf, 2003: Assessment fo the reliability of wave observations from voluntary observing ships: insights from the validation of a global wind wave climatology based on voluntary observing ship data. J. Geophys. Res. – Oceans, 108(C7), 3236, doi:10,1029/2002JC001437 6. V.Grigorieva, S.Gulev and A.Sterl: Global climatology of wind waves based on VOS observation. Symposium Internatiol “In route to GODAE”,205-206, 2002,Biarritz 7. Gulev, S.K., V. Grigorieva, K.M.Selemenov, and O. Zolina, 2003: Evaluation of surface winds and waves from volunatary observing ship data. Advances in Applications of Marine Climatology, WMO, Geneva, 53-67 8. Gulev, S.K. and V. Grigorieva, 2004: Last century changes in ocean wind wave height from global visual wave data. Geophys. Res. Lett., 31, L24302, doi: 10.1029/2004GL021040 9. Gulev, S.K. and V. Grigorieva, 2006: Variability of the winter wind waves and swell in the North Atlantic and North Pacific as revealed by the Voluntary Observing Ship data, J. Climate, 19, p. 5667-5685 10. Lopatoukhin L.J., Rozhkov V.A., Ryabinin V.E., Swail V.R, Boukhanovsky A.V., Degtyarev A.B. Estimation of extreme wind wave heights.// World Meteorological Organisation (WMO). WMO/TD-No. 1041, 2000, JCOMM Technical Report, 71p. 11. Lopatoukhin, L.J., Boukhanovsky, A.V., 2004, "Freak Wave Generation and their Probability", International Shipbuilding Progress, vol. 51, no. 2/3, pp. 157-171. 12. Pelinofsky E., Kharif T., Talipova T., Damiani T. Nonlinear Wave Focusing as a Mechanism of the Freak Wave Generation in the Ocean. Proceedings of a Workshop “Rogue waves”. Brest, France, 2000, p.193-204. 13. Rogue waves 2000. Proceedings of a Workshop organized by Ifremer. Brest, France. 2930 November. 2000. 396p. 14. Woodruff, S. D., H. F. Diaz, J. D. Elms, and S. J. Worley, COADS release 2. Data and metadata enhancements for improvements of marine surface fluxfields, Phys. Chem. Earth, 23, 517– 526, 1998. 15. Woodruff, S. D., S. J. Worley, J. A. Arnott, H. F. Diaz, J. D. Elms, M. Jackson, S. J. Lubker, and D. E. Parker, COADS Updates and the Blend With the UK Meteorological Office Marine Data Bank, in Achievements in Marine Climatology, edited by V. Swail, pp. 23– 31