Video in French

Ocean Wave Model - ECWAM

Purpose

The ECMWF Ocean Wave Model (ECWAM) describes the development and evolution of wind generated surface waves and their height, direction and period.  Its domain extends across the full globe. 

ECWAM is coupled to:

• the atmospheric forecast model at all configurations (medium range ensemble, extended range ensemble, Seasonal),
• the dynamic ocean model NEMO.   (Note: the IFS analysis and the short range forecasts that are part of the IFS analysis cycle are not coupled with ECWAM).
  

ECWAM is solely concerned with ocean wave forecasting and does not model the ocean itself: dynamical modelling of the ocean is done by NEMO.

Structure

 ECWAM evaluates the 2-dimensional surface wave spectrum, in both oceanic and coastal (but not inshore) waters.  This describes how much wave energy is present for given sea wave frequencies and associated propagation directions.  The part of the spectrum under the direct influence of the local wind is called “wind-wave" or "wind-sea”; the remaining part is usually referred to as “swell”.

Changes in the wave spectrum are derived from the processes of:

• wave advection,
• wave refraction,
• wind-wave generation,
• wave dissipation due to white capping and bottom friction,
• non-linear wave interactions.

Interaction with atmospheric and ocean models

ECWAM has two-way interaction with the Atmospheric models:

• ECWAM supplies surface roughness (according to the forecast sea state) which is used in the evaluation of boundary layer winds.
• Atmospheric models supply surface wind conditions, and these dominate sea-surface wave development.

ECWAM has two-way interaction with NEMO and its sub-program LIM2:

• ECWAM supplies surface stress, Stokes drift, and turbulent energy flux to the ocean surface.
• NEMO and the LIM2 sub-program supply information on the concentration and thickness of ice.

Note: ECMWF uses LIM2 which is an earlier version of the Louvain-la-Neuve sea ice model currently available (Version 3.6)

Wave Data Assimilation

• ECWAM assimilates space-borne altimeter wave height data.
• ECWAM does not assimilate buoy wave data; instead, these data serve as an independent check on the quality of modelled wave parameters.

Output from ECWAM

ECWAM is run as:

• HRES WAM twice daily giving forecasts to Day10 associated with HRES forecasts based on 00UTC, 06UTC, 12UTC and 18UTC data times.  Global coverage 0.25^o x 0.25^o latitude-longitude grid (Ocean wave model ensemble 15-day forecast (Set IV ENS-WAM).
• ENS-WAM twice daily giving forecasts to Day15 associated with ensemble forecasts  based on 00UTC and 12UTC data times.  Global coverage 0.25^o x 0.25^o latitude-longitude grid (Ocean wave model ensemble 15-day forecast (Set IV ENS-WAM).
• Extended ENS-WAM daily giving extended range forecasts from Day16 to Day46 associated with the extended range forecasts based on 00UTC data times.  Global coverage 0.5^o x 0.5^o latitude-longitude grid (Atmospheric Model Ensemble extended forecast (Set VI - ENS extended), section VI-v-c).
• SEAS-WAM monthly for forecasts to 7 months ahead associated with the seasonal forecast model (System 5).  Global coverage 1.0^o x 1.0^o latitude-longitude grid (Atmospheric model Seasonal 7-month forecast (Set V - SEAS), section V-v-e).
• SEAS-WAM quarterly for forecasts to 1 year ahead associated with the seasonal forecast model (System 5).  Global coverage 1.0^o x 1.0^o latitude-longitude grid (Atmospheric model Seasonal 7-month forecast (Set V - SEAS), section V-v-e).

       Output is in the form of wave and swell height, direction and period.  Also available are wave energy flux, mean direction and magnitude (important for assessment of the impact of the waves on coastlines and offshore structures). 

Graphical and chart presentation of wave forecasts

Fig2.2-1: ECMWF forecast entry page.

Wave output on ecCharts.

Fig2.2-2: Procedure to load wave parameter charts on ecCharts.  Click on "Show Layers List" icon (1); Select "Add Layers" option (2); Input Wave into the "Layer select" box (3); Select desired chart by clicking on the icon (4).

Fig2.2-3: Wave parameter charts available on ecCharts (see Fig2.2-1and Fig2.2-2 above) and may be displayed by clicking on the desired icon.

Wavegram output on ecCharts

Wavegrams are also available to show a time series of significant wave height, mean wave direction, and mean wave period for any sea location.

Fig2.2-4: Wavegrams for any oceanic location are available on ecCharts. Choose location using the Probe icon (1); Click on "Views" (2); Select "Meteograms" option on the dropdown menu (3); Select "More" on the option page that appears(4); input Wave into the "Meteogram select" box (5); select desired chart(s) by clicking on the icon (6).

Wave output on Opencharts

Open Charts.

Fig2.2-5: Menu to select wave parameter charts from Open Charts (See Fig2.2-1 above).  Select "Range" (here medium and extended ranges); Select "Ocean Waves".

Fig2.2-6: Sample wave chart. Click on any point to produce a meteogram and/or wavegram.

Fig2.2-7:  Wavegram on Opencharts associated with Fig2.2-6 at location near 52N 27W.

Fig2.2-8: Wave parameter charts available on Open Charts.  These may be displayed by clicking on the desired icon (the above ensemble products are available at post-processing steps, 12-hourly from T+0h to T+168h).

Fig2.2-9: Seasonal forecast charts for Tropical Storm, Hurricane, Typhoons frequencies are available on Open Charts by selecting Long option in the menu (See Fig2.2-5)  and then may be displayed by clicking on the desired icon.

Fig2.2-10: Sequence of ocean wave forecasts.  Significant wave height forecast (colours) and 10m wind (arrows) from data time 00UTC 4 January 2014, step 12 hours.  Wave heights at 1.25m intervals as scale.

Swell propagates outwards well away from the source.  Increased swell (e.g. reaching a coast) can give forewarning of a storm system well before any indication in the atmosphere.  Fishermen have long used the arrival of long-period swell as an indication of an approaching storm even if the sky is clear.  Surfers often benefit from significantly large swell in calm conditions well away from the swell source region.

Fig2.2-11:  180h forecast for significant wave height (contours) for all waves with periods between 21 and 25 seconds (shading), initialised at 00UTC 2 December 2016.  The highest significant wave heights (contours) are still confined to the storm location in the Atlantic south of Iceland, while long waves from that storm are already affecting coastlines from Iberia to South Greenland (coloured). 

Waves with different periods propagate with different speeds - longer periods travel fastest.  These can be tracked through the forecast period and areas where different wave trains potentially interact can be identified.

Fig2.2-12:  Chart showing forecast significant wave heights for several ranges of wave periods (Blue,10-12s; Green, 12-14s; Yellow,14-17s; Red,17-21s).  Forecast data based on data time 00UTC 25 October 2017.  The faster southward propagation of the long period waves over the shorter period waves from their source off NW Africa is clear.

The Extreme Forecast Index (EFI) can be used to indicate the significance of forecast significant wave heights when compared with the range of wave heights that might usually be expected as defined by the M-climate. 

Fig2.2-13: In this example the colours west of Ireland denote a low-point in wave heights, or potentially a form of 'weather window' for certain types of marine/shipping operations.  Equally this EFI can signify periods with anomalously big waves (yellow to red shading).

Considerations when using output from ECWAM

Interaction of wind-sea and swell

Use of the mean wave height and direction is the simplest method of describing the forecast wave regime in a given area and it is easy to be beguiled into just using this output for forecasts to customers.  However, the mean wave direction and height is made up of contributions from wind-sea and swell with different wave periods and they interact in a complex manner.  It is important to investigate the forecast wind-sea and swell separately to give an understanding of likely sea conditions in an area (e.g. for a ship requiring a particularly smooth passage) or at a location (e.g. an oil rig).  

When wind-sea and swell move in similar directions the wave heights can give information on the likely sea state as one is superimposed on the other, particularly where both have a significant and comparable wave height.  On occasion the swell and wind-sea may be moving in opposite directions (an opposing sea) and wave heights give information on the likely rougher sea state to be expected.  Often the wind-sea and swell are at right-angles (a cross sea).  Where the wind-sea and swell heights are similar the sea can be very disturbed and difficult for shipping.   An illustration is given in Figs2.2-14(a to e) 

Fig2.2-14(a): ecChart of mean wave direction (wave height is indicated by the length of the arrow).  On ecCharts, wave height may be shown by use of the probe tool or more graphically by superimposing mean wave heights).  This chart gives an overview of wave conditions.  Northwesterly waves (i.e. moving towards the northwest) are indicated near point A.  Easterly waves i.e. moving towards the east) are indicated near point C.   However, it is important to investigate the contributions to the mean wave directions and heights from inspection of the wind-sea and swell at this time.

Fig2.2-14(b): The forecast wind-sea has developed in response to the forecast winds around a depression in mid-Atlantic.  Waves move northwestwards near point A and southeastwards near point B.  The length of the arrows near points A & B suggest wave heights are around 3m (wave heights are also available as charts, not shown here).  Near point C wind-sea waves are relatively small and move towards the east-northeast.

Fig2.2-14(c): The forecast total swell has been developed in response to earlier weather systems elsewhere and has propagated across the Atlantic.  Swell is moving northwards near point A and southwestwards near point C with arrow length suggesting wave heights of around 2m (wave heights are also available as charts, not shown here).  Near point B swell waves are relatively small and move towards the southeast.

Fig2.2-14(d): The forecast wind-sea (blue) and swell (black) shown on a single chart.  To the north of point B the wind-sea and swell waves have a similar direction of travel; to the east of point B wind-sea dominates with only weak swell contribution but almost at a right-angle.  Near points A and C the wind-sea and swell waves differ widely in direction but with similar heights (a cross sea). 

Fig2.2-14(e):The forecast mean wave directions derived from the wind-sea and mean swell (as shown in Fig2.2-14(a)) superimposed on the previous chart (Fig2.2-14(d)).  This illustrates the important additional information that is gained from consideration of the wind-sea and mean swell forecasts together.  The mean wave directions (Fig2.2-14(a)) give no indication of that a sea passage to the west of Portugal is likely to be through confused rough seas.

Waves and swell with a long period

Large swell waves with a long period breaking on a beach slope tend to have a large swash with water washed well up the beach.  This is often unexpected and takes people by surprise and can cause damage and casualties.  Users should note when large waves with a long period are forecast to run onto an exposed coast.  Extreme forecast index products for waves can alert users to the potential problem. 

Fig2.2-15: Wave energy forecast VT 09UTC 16 Sep 2023, DT 12UTC 15 Sep 2023.   The chart shows extreme wave energy flux (~1300KW/m) being driven towards the exposed southern coast of South Africa.

Fig2.2-16: Wave energy forecast VT 09UTC 16 Sep 2023, DT 12UTC 15 Sep 2023.  The chart shows the extreme forecast index for significant wave height at 0.9 to 1.0 alerting to the significant nature of the ocean swell.

Fig2.2-17: Wave ENSgram DT DT 00UTC 16 Sep 2023. Significant wave heights were forecast to exceed 9 m with wave period above 15 s.  These values were also observed at a nearby buoy.  

Very large waves and swell were induced by a deep depression in southern Atlantic associated with very strong winds. The waves became larger as they approached the coasts and coincided with spring tides. Significant wave heights exceeded 8 m in many places suggesting maximum waves were significantly higher.  There was considerable coastal damage and some loss of life.

Sea-surface currents

The interaction of waves with sea-surface currents is not yet included in the operational version of the model.  In particular areas, (e.g. Gulf Stream or Agulhas current), the current effect may give rise to localised changes of up to a metre in the wave height.

Users should also note that whilst ECMWF does provide some ocean current output, from its ocean model (as "sea water velocity fields"), the current 1/4 degree resolution of that model is insufficient to allow strong gradients in western boundary currents to be captured. This means that stronger currents that are observed around the world tend to be underestimated in this output, sometimes substantially so.

Shoaling

Shoaling is the deformation of waves as they move from the ocean into shallow waters causing the waves to become steeper, increase in height, and have shorter wavelength .  The basic equations in ECWAM do represent the effect when the waves propagate from deep to shallow water, but the effect is not dramatic over most coastal waters.   Waves inshore and at the beach, where shoaling is very strong, are not represented since the resolution of ECWAM (~ 10km) cannot represent the actual beach slope.  Wave products near coasts, and, to a lesser extent, within small and enclosed basins (e.g. Baltic Sea) may be of lower quality than for the open ocean.  This may be due to incomplete resolution the detail of the coast by the land-sea mask.  Small islands too may not be identified and hence allow waves to propagate unhindered.   Note, however, that the wave model has a scheme that attempts to represent the impact of unresolved islands on the global propagation of waves. 

Fig2.2-18: An example chart of wind-wave and swell. Some shoaling is possible towards the French and British coasts as the sea becomes less deep but forecast values cannot be absolutely relied upon. See Fig2.2-19 for detail around the Azores.  No parameters are shown on coasts nor where ice cover >30% (i.e. where some of the grid points used in interpolation of wave data for display are on land or ice. Users should identify whether the ice areas or coastal zones are the cause. See section regarding wave parameters near sea ice) 

 Fig2.2-19: The same example chart of wind-wave and swell as in Fig2.2-18, magnified near the Azores.  There are some areas around the islands where wave parameters are not forecast (i.e. where some of the grid points used in interpolation of wave data are on land) but the detail of coast may not be fully resolved.  ECWAM shows re-build of wind waves to the lee of the islands as the wind fetch increases and also the penetration of larger waves through the inter-island straits.

Near sea ice

Sea ice is not static but forms or extends with low air temperature or sea-surface temperatures, and can move with winds and sea current.  NEMO passes information to ECWAM regarding the extent and movement of the sea ice field forecast by LIM2, allowing a more realistic definition of what is open sea throughout the forecast period.  In the current operational version of the wave model, the interaction between waves and sea-ice is not actually represented.  Where sea ice cover >30% all wave parameters are set to missing (i.e. no valid values).  Wave products near ice-edges may be of lower quality than for the open ocean.  This may be due to uncertainty in sea-ice cover, or in the detail of an ice edge and consequently also in the boundary of the water area.  Spurious areas of ice or incorrect extent of ice will act as if a coastline or island and stop waves from propagating correctly, possibly decaying the waves completely and incorrectly sheltering an otherwise exposed location.

Fig2.2-20: Illustration of the importance of distinguishing between ice cover and shallow water when an area of wave parameters is missing.

Fig2.2-21: Significant height of combined wind waves and swell (Hs).  The coloured areas show the difference between the heights derived from the 2d spectra (used where >30% sea ice cover is forecast) and from the wave model (as if open sea).  Some large differences are evident, illustrating the need to treat the values of wave height with caution where sea ice is present.

Waves near tropical storms

In the IFS, there is an active two-way coupling between the atmosphere and ocean waves - surface wind stress generates the waves and in turn the waves modulate the wind stress.  The ECWAM generally forecasts realistic wave parameters (wave height, period etc). 

With sufficiently strong winds, the drag on the low-level air flow is modelled using the Charnock parameter which is used to specify an aerodynamical roughness length scale.  In the IFS, this roughness length scale tends to increase with 10m wind speed.  However, observational evidence suggests that for exceptionally strong winds the coupling between the ocean surface and the wind becomes less efficient at transferring momentum.  So, for exceptionally strong (mean) 10m winds (33m/s) the roughness length scale should be capped.  This limitation avoids the effect of too much drag on the lower atmosphere and enables more realistic (stronger) winds to be forecast in the vicinity of relatively intense tropical cyclones.  It should also be noted that tropical cyclone development with strengthened winds has only a limited effect upon the size and character of waves developed by the Wave model.  The model change to limit the roughness length scale at very high wind speeds was introduced in cycle 47r1 in June 2020.  It is important not to confuse to the aerodynamic roughness length scale with the actual roughness of the sea.  The aerodynamic roughness length scale is a parameter concerning the amount of momentum that is exchanged between the atmosphere and the ocean surface via the surface stress.   The actual roughness of the sea is concerned with wave heights, wave steepness and how much wave breaking is happening.

When considering forecast wave parameters in the vicinity of typhoons, hurricanes etc., it should be remembered that IFS still has difficulties in producing some intense tropical cyclones and their subsequent motion.

Additional Sources of Information

(Note: In older material there may be references to issues that have subsequently been addressed)

• Definitive information on the ECMWF WAVE MODEL products and availability can be found:
  • for HRES WAM in the Catalogue of ECMWF real-time products (Set II).
  • for ENS WAM in the Catalogue of ECMWF real-time products (Set IV)
• Read an in depth view of the structure of ECWAM which gives a description of the theory behind the ECWAM model.
• Read about techniques regarding satellite measurement of wave height.
• Read more on ocean wave modelling and model output parameters.
• Fact sheet - Ocean wave forecasting.
• ECMWF - Lead centre for wave forecast verification and information on national models and wave verification results. 
• Watch a comprehensive lecture on ocean wave forecasting at ECMWF or air-sea Interaction and earth system modelling.