15.2 Numerical Models in Oceanography Numerical models are used for many purposes in oceanography. For our purpose we can divide models into two classes: Mechanistic models are simplified models used for studying processes. Because the models are simplified, the output is easier to interpret than output from more complex models. Many different types of simplified models have been developed, including models for describing planetary waves, the interaction of the flow with seafloor features, or the response of the upper ocean to the wind. These are perhaps the most useful of all models because they provide insight into the physical mechanisms influencing the ocean. The development and use of mechanistic models is, unfortunately, beyond the scope of this book. Simulation models are used for calculating realistic circulation of oceanic regions. The models are often very complex because all important processes are included, and the output is difficult to interpret. The first simulation models were developed by Kirk Bryan and Michael Cox at the Geophysical Fluid Dynamics laboratory in Princeton. Their model (Bryan, 1969) calculated the 3dimensional flow in the ocean using the continuity and momentum equation with the hydrostatic and Boussinesq approximations and a simplified equation of state. Such models are called primitive equation models because they use the most basic, or primitive form of the equations of motion. The equation of state allows the model to calculate changes in density due to fluxes of heat and water through the sea surface, so the model includes thermodynamic processes. The BryanCox model used large horizontal and vertical viscosity and diffusion to eliminate turbulent eddies having diameters smaller about 500 km, which is a few grid points in the model. It had complex coastlines, smoothed seafloor features, and a rigid lid. The rigid lid was needed to eliminate oceansurface waves, such as tides and tsunamis, that move far too fast for the coarse time steps used by all simulation models. The rigid lid had, however, disadvantages. Islands substantially slowed the computation, and the seafloor features were smoothed to eliminate steep gradients. The first simulation model was regional. It was quickly followed by a global model (Cox, 1975) with a horizontal resolution of 2° and with 12 levels in the vertical. The model ran far too slowly even on the fastest computers of the day, but it laid the foundation for more recent models. The coarse spatial resolution required that the model have large values for viscosity, and even regional models were too viscous to have realistic western boundary currents or mesoscale eddies. Since those times, the goal has been to produce models with ever finer resolution, more realistic modeling of physical processes, and better numerical schemes. Computer technology is changing rapidly, and models are evolving rapidly. The output from the most recent models of the North Atlantic, which have resolution of 0.03° look very much like the real ocean. Models of other areas show previously unknown currents near Australia and in the South Atlantic. Ocean and Atmosphere Models use very different spacing of grid points. As a result, ocean modeling lags about a decade behind atmosphere modeling. Dominant ocean eddies are 1/30 the size of dominant atmosphere eddies (storms). But, ocean features evolve at a rate that is 1/30 the rate in the atmosphere. Thus ocean models running for say one year have 30 × more horizontal grid points than the atmosphere, but they have 1/30 the number of time steps. Both have about the same number of grid points in the vertical. As a result, ocean models run 30 times slower than atmosphere models of the same complexity.


Department of Oceanography, Texas A&M University Robert H. Stewart, stewart@ocean.tamu.edu All contents copyright © 2005 Robert H. Stewart, All rights reserved Updated on November 3, 2008 