15.4 Coastal Models
The great economic importance of the coastal zone has led to the development of many different numerical models for describing coastal currents, tides, and storm surges. The models extend from the beach to the continental slope, and they can include a free surface, realistic coasts and bottom features, river runoff, and atmospheric forcing. Because the models don't extend very far into deep water, they need additional information about deep-water currents or conditions at the shelf break.
The many different coastal models have many different goals, and many different implementations. Several of the models described above, including MOM and MICOM, have been used to model coastal processes. But many other specialized models have also been developed. Heaps (1987), Lynch et al., (1996), and Haidvogel (1998) provide good overviews of the subject. Rather than look at a menu of models, let's look at two typical models.
Princeton Ocean Model POM
where z = h (x, y, t) is the sea surface, and z = -H (x, y) is the bottom.
Sub-grid turbulence is parameterized using a closure scheme proposed by Mellor and Yamada (1982) whereby eddy diffusion coefficients vary with the size of the eddies producing the mixing and the shear of the flow.
The model is driven by wind stress and heat and water fluxes from meteorological models. The model uses known geostrophic, tidal, and Ekman currents at the outer boundary.
The model has been used to calculate the three-dimensional distribution of velocity, salinity, sea level, temperature, and turbulence for up to 30 days over a region roughly 100-1000 km on a side with grid spacing of 1-50 km.
Dartmouth Gulf of Maine Model
The model uses roughly 13,000 triangles to cover the Gulf of Maine and nearby waters of the North Atlantic (Figures 15.3). Minimum size of the elements is roughly one kilometer. The model has 10 to 40 horizontal layers. The vertical spacing of the layers is not uniform. Layers are closer together near the top and bottom and they are more widely spaced in the interior. Minimum spacing is roughly one meter in the bottom boundary layer.
The model integrates the three-dimensional, primitive equations, in shallow water form. The model has a simplified equation of state and a depth-averaged continuity equation, and it uses the hydrostatic and Boussinesq assumptions. Sub-grid mixing of momentum, heat and mass is parameterized using the Mellor and Yamada (1982) turbulence-closure scheme which gives vertical-mixing coefficients that vary with stratification and velocity shear. Horizontal mixing coefficients were calculated from Smagorinski (1963). A carefully chosen, turbulent, eddy viscosity is used in the bottom boundary layer. The model is forced by wind, heating, and tidal forcing from the deep ocean.
The model is spun up from rest for a few days using a specified density field at all grid points, usually from a combination of CTD data plus historical data. This gives a velocity field consistent with the density field. The model is then forced with local winds and heat fluxes to calculate the evolution of the density and velocity fields.
Comments on Coastal Models
Hackett et al., (1995) compared the ability of two of the five models to describe observed flow on the Norwegian shelf. They conclude that
Calculating storm surges is not easy. Here are some reasons, in a rough order of importance.
To reduce errors, models are tuned to give results that match conditions seen in past storms. Unfortunately, those past conditions are not well known. Changes in sea level and wind speed are rarely recorded accurately in storms except at a few, widely paced locations. Yet storm-surge heights can change by more than a meter over distances of tens of kilometers.
Despite these problems, models give very useful results. Let's look at one, commonly-used model.
Sea, Lake, and Overland Surges Model SLOSH
The model is the result of a lifetime of work by Chester Jelesnianski. In developing the model, Jelesnianski paid careful attention to the relative importance of errors in the model. He worked to reduce the largest errors, and ignored the smaller ones. For example, the distribution of winds in a hurricane is not well known, so it makes little sense to use a spatially varying drag coefficient for the wind. Thus, Jelesnianski used a constant drag coefficient in the air, and a constant eddy stress coefficient in the water.
SLOSH calculates water level from depth-integrated, quasi-linear, shallow water equations. Thus it ignores stratification. It also ignores river in flow, rain, and tides. The latter may seem strange, but the model is designed for forecasting. The time of landfall cannot be forecast accurately, and hence the height of the tides is mostly unknown. Tides can be added to the calculated surge, but the nonlinear interaction of tides and surge is ignored.
The model is forced by idealized hurricane winds. It needs only atmospheric pressure at the center of the storm, the distance from the center to the area of maximum winds, the forecast storm track and speed along the track.
In preparation for hurricanes coming ashore near populated areas, the model has been adapted for 27 basins from Boston Harbor Massachusetts to Laguna Madre Texas. The model uses a fixed polar mesh. Mesh spacing begins with a fine mesh near the pole, which is located near the coastal city for which the model is adapted. The grid stretches continuously to a coarse mesh at distant boundaries of a large basin. Such a mesh gives high resolution in bays and near the coast where resolution is most needed. Using measured depths at sea and elevations on land, the model allows flooding of land, overtopping of levees and dunes, and sub-grid flow through channels between offshore islands.
Sea level calculated from the model has been compared with heights measured by tide gauges for 13 storms, including Betsy (1965), Camile (1969), Donna (1960), and Carla (1961). The overall accuracy is ± 20%.
Advanced Circulation Model ADCIRC
The model is forced by:
The model successfully forecast the Hurricane Katrina storm surge, giving values in excess of 6.1 m near New Orleans.
|Department of Oceanography, Texas A&M University
Robert H. Stewart, email@example.com
All contents copyright © 2005 Robert H. Stewart,
All rights reserved
Updated on November 5, 2008