Wind wave model


In fluid dynamics, wind wave modeling describes the effort to depict the sea state and predict the evolution of the energy of wind waves using numerical techniques. These simulations consider atmospheric wind forcing, nonlinear wave interactions, and frictional dissipation, and they output statistics describing wave heights, periods, and propagation directions for regional seas or global oceans. Such wave hindcasts and wave forecasts are extremely important for commercial interests on the high seas. For example, the shipping industry requires guidance for operational planning and tactical seakeeping purposes.
For the specific case of predicting wind wave statistics on the ocean, the term ocean surface wave model is used.
Other applications, in particular coastal engineering, have led to the developments of wind wave models specifically designed for coastal applications.

Historical overview

Early forecasts of the sea state were created manually based upon empirical relationships between the present state of the sea, the expected wind conditions, the fetch/duration, and the direction of the wave propagation. Alternatively, the swell part of the state has been forecasted as early as 1920 using remote observations.
During the 1950s and 1960s, much of the theoretical groundwork necessary for numerical descriptions of wave evolution was laid. For forecasting purposes, it was realized that the random nature of the sea state was best described by a spectral decomposition in which the energy of the waves was attributed to as many wave trains as necessary, each with a specific direction and period. This approach allowed to make combined forecasts of wind seas and swells. The first numerical model based on the spectral decomposition of the sea state was operated in 1956 by the French Weather Service, and focused on the North Atlantic. The 1970s saw the first operational, hemispheric wave model: the spectral wave ocean model at the Fleet Numerical Oceanography Center.
First generation wave models did not consider nonlinear wave interactions. Second generation models, available by the early 1980s, parameterized these interactions. They included the "coupled hybrid" and "coupled discrete" formulations. Third generation models explicitly represent all the physics relevant for the development of the sea state in two dimensions. The wave modeling project, an international effort, led to the refinement of modern wave modeling techniques during the decade 1984-1994.
Improvements included two-way coupling between wind and waves, assimilation of satellite wave data, and medium-range operational forecasting.
Wind wave models are used in the context of a forecasting or hindcasting system. Differences in model results arise from: differences in wind and sea ice forcing, differences in parameterizations of physical processes, the use of data assimilation and associated methods, and the numerical techniques used to solve the wave energy evolution equation.
In the aftermath of World War II, the study of wave growth garnered significant attention. The global nature of the war, encompassing battles in the Pacific, Atlantic, and Mediterranean seas, necessitated the execution of landing operations on enemy-held coasts. Safe landing was paramount, given that choppy waters posed the danger of capsizing landing craft. Consequently, the precise forecasting of weather and wave conditions became essential, prompting the recruitment of meteorologists and oceanographers by the warring nations.
During this period, both Japan and the United States embarked on wave prediction research. In the U.S., comprehensive studies were carried out at the Scripps Institution of Oceanography affiliated with the University of California. Under the guidance of Harald Svedrup, Walter Munk devised an avant-garde wave calculation methodology for the United States Navy and later refined this approach for the Office of Naval Research.
This pioneering effort led to the creation of the significant wave method, which underwent subsequent refinements and data integrations. The method, in due course, came to be popularly referred to as the SMB method, an acronym derived from its founders Sverdrup, Munk, and Charles L. Bretschneider.
Between 1950 and 1980, various formulae were proposed. Given that two-dimensional field models had not been formulated during that time, studies were initiated in the Netherlands by Rijkswaterstaat and the Technische Adviescommissie voor de Waterkeringen to discern the most appropriate formula to compute wave height at the base of a dike. This work concluded that the 1973 Bretschneider formula was the most suitable. However, subsequent studies by Young and Verhagen in 1997 suggested that adjusting certain coefficients enhanced the formula's efficacy in shallow water regions.

General strategy

Input

A wave model requires as initial conditions information describing the state of the sea. An analysis of the sea or ocean can be created through data assimilation, where observations such as buoy or satellite altimeter measurements are combined with a background guess from a previous forecast or climatology to create the best estimate of the ongoing conditions. In practice, many forecasting system rely only on the previous forecast, without any assimilation of observations.
A more critical input is the "forcing" by wind fields: a time-varying map of wind speed and directions. The most common sources of errors in wave model results are the errors in the wind field. Ocean currents can also be important, in particular in western boundary currents such as the Gulf Stream, Kuroshio or Agulhas current, or in coastal areas where tidal currents are strong. Waves are also affected by sea ice and icebergs, and all operational global wave models take at least the sea ice into account.
Image:Wave model current induced refraction.png|300px|thumb|This figures show an example of the effects of currents on the wave heights. This example is adapted from scientific paper published in the Journal of Physical Oceanography. The top panels show the tidal currents at 3 AM and 11 AM on 28 October 2008, off the West coast of France, around the island of Ouessant, which lies 20 km from the mainland. The bottom panel show the heights and directions of waves, computed with the numerical model WAVEWATCH III, using a triangular mesh with variable resolution. The strong currents south of Ouessant deflect the waves away from the measuring buoy at low tide.

Representation

The sea state is described as a spectrum; the sea surface can be decomposed into waves of varying frequencies using the principle of superposition. The waves are also separated by their direction of propagation. The model domain size can range from regional to the global ocean. Smaller domains can be nested within a global domain to provide higher resolution in a region of interest. The sea state evolves according to physical equations – based on a spectral representation of the conservation of wave action – which include: wave propagation / advection, refraction, shoaling, and a source function which allows for wave energy to be augmented or diminished. The source function has at least three terms: wind forcing, nonlinear transfer, and dissipation by whitecapping. Wind data are typically provided from a separate atmospheric model from an operational weather forecasting center.
For intermediate water depths the effect of bottom friction should also be added. At ocean scales, the dissipation of swells - without breaking - is a very important term.

Output

The output of a wind wave model is a description of the wave spectra, with amplitudes associated with each frequency and propagation direction. Results are typically summarized by the significant wave height, which is the average height of the one-third largest waves, and the period and propagation direction of the dominant wave.

Coupled models

Wind waves also act to modify atmospheric properties through frictional drag of near-surface winds and heat fluxes. Two-way coupled models allow the wave activity to feed back upon the atmosphere. The European Centre for Medium-Range Weather Forecasts coupled atmosphere-wave forecast system described below facilitates this through exchange of the Charnock parameter which controls the sea surface roughness. This allows the atmosphere to respond to changes in the surface roughness as the wind sea builds up or decays.

Examples

WAVEWATCH

The operational wave forecasting systems at NOAA are based on the WAVEWATCH III model. This system has a global domain of approximately 50 km resolution, with nested regional domains for the northern hemisphere oceanic basins at approximately 18 km and approximately 7 km resolution. Physics includes wave field refraction, nonlinear resonant interactions, sub-grid representations of unresolved islands, and dynamically updated ice coverage. Wind data is provided from the GDAS data assimilation system for the GFS weather model. Up to 2008, the model was limited to regions outside the surf zone where the waves are not strongly impacted by shallow depths.
The model can incorporate the effects of currents on waves from its early design by Hendrik Tolman in the 1990s, and is now extended for near shore applications.

WAM

The wave model WAM was the first so-called third generation prognostic wave model where the two-dimensional wave spectrum was allowed to evolve freely with no constraints on the spectral shape. The model underwent a series of software updates from its inception in the late 1980s. The last official release is Cycle 4.5, maintained by the German Helmholtz Zentrum, Geesthacht.
ECMWF has incorporated WAM into its deterministic and ensemble forecasting system., known as the Integrated Forecast System. The model currently comprises 36 frequency bins and 36 propagation directions at an average spatial resolution of 25 km. The model has been coupled to the atmospheric component of IFS since 1998.