Princeton Ocean Model

Last updated

The Princeton Ocean Model (POM) is a community general numerical model for ocean circulation that can be used to simulate and predict oceanic currents, temperatures, salinities and other water properties. POM-WEB [1] and POMusers.org [2]

Contents

Development

The model code [3] was originally developed at Princeton University (G. Mellor and Alan Blumberg) in collaboration with Dynalysis of Princeton (H. James Herring, Richard C. Patchen). The model incorporates the Mellor–Yamada turbulence scheme developed in the early 1970s by George Mellor and Ted Yamada; this turbulence sub-model is widely used by oceanic and atmospheric models. At the time, early computer ocean models such as the Bryan–Cox model (developed in the late 1960s at the Geophysical Fluid Dynamics Laboratory, GFDL, and later became the Modular Ocean Model, MOM)), were aimed mostly at coarse-resolution simulations of the large-scale ocean circulation, so there was a need for a numerical model that can handle high-resolution coastal ocean processes. The Blumberg–Mellor [4] model (which later became POM) thus included new features such as free surface to handle tides, sigma vertical coordinates (i.e., terrain-following) to handle complex topographies and shallow regions, a curvilinear grid to better handle coastlines, and a turbulence scheme to handle vertical mixing. At the early 1980s the model was used primarily to simulate estuaries such as the Hudson–Raritan Estuary (by Leo Oey) and the Delaware Bay (Boris Galperin), but also first attempts to use a sigma coordinate model for basin-scale problems have started with the coarse resolution model of the Gulf of Mexico (Blumberg and Mellor) and models of the Arctic Ocean (with the inclusion of ice-ocean coupling by Lakshmi Kantha and Sirpa Hakkinen).

In the early 1990s when the web and browsers started to be developed, POM became one of the first ocean model codes that were provided free of charge to users through the web. [5] The establishment of the POM users group and its web support (by Tal Ezer) resulted in a continuous increase in the number of POM users which grew from about a dozen U.S. users in the 1980s to over 1000 users in 2000 and over 4000 users by 2009; there are users from over 70 different countries. In the 1990s the usage of POM expands to simulations of the Mediterranean Sea (Zavatarelli [6] ) and the first simulations with a sigma coordinate model of the entire Atlantic Ocean [7] for climate research (Ezer). The development of the Mellor–Ezer optimal interpolation data assimilation scheme [8] that projects surface satellite data into deep layers allows the construction of the first ocean forecast systems for the Gulf Stream [9] and the U.S. east coast running operationally at the NOAA's National Weather Service (Frank Aikman and others [10] ). Operational forecast system for other regions such as the Great Lakes, the Gulf of Mexico (Oey), the Gulf of Maine (Huijie Xue) and the Hudson River (Blumberg) followed. For more information on applications of the model, see the searchable database of over 1800 POM-related publications. [11]

Derivatives and other models

In the late 1990s and the 2000s many other terrain-following community ocean models have been developed; some of their features can be traced back to features included in the original POM, other features are additional numerical and parameterization improvements. Several ocean models are direct descendants of POM such as the commercial version of POM known as the estuarine and coastal ocean model (ECOM), the navy coastal ocean model (NCOM) and the finite-volume coastal ocean model (FVCOM). Recent developments in POM include a generalized coordinate system that combines sigma and z-level grids (Mellor and Ezer), inundation features that allow simulations of wetting and drying (e.g., flood of land area) (Oey), and coupling ocean currents with surface waves (Mellor). Efforts to improve turbulent mixing also continue (Galperin, Kantha, Mellor and others).

Users' meetings

POM users' meetings were held every few years, and in recent years the meetings were extended to include other models and renamed the International Workshop on Modeling the Ocean (IWMO). Meeting Pages: [12] List of meetings:

Reviewed papers from the IWMO meetings are published by Ocean Dynamics [13] in special issues (IWMO-2009 Part-I, [14] IWMO-2009 Part-II, [15] IWMO-2010, [16] IWMO-2011, [17] IWMO-2012, [18] IWMO-2013, [19] IWMO-2014 [20] ).

Related Research Articles

<span class="mw-page-title-main">Geophysical Fluid Dynamics Laboratory</span> American physics research center

The Geophysical Fluid Dynamics Laboratory (GFDL) is a laboratory in the National Oceanic and Atmospheric Administration (NOAA) Office of Oceanic and Atmospheric Research (OAR). The current director is Venkatachalam Ramaswamy. It is one of seven Research Laboratories within NOAA's OAR.

A stochastic differential equation (SDE) is a differential equation in which one or more of the terms is a stochastic process, resulting in a solution which is also a stochastic process. SDEs have many applications throughout pure mathematics and are used to model various behaviours of stochastic models such as stock prices, random growth models or physical systems that are subjected to thermal fluctuations.

<span class="mw-page-title-main">Lattice Boltzmann methods</span> Class of computational fluid dynamics methods

The lattice Boltzmann methods (LBM), originated from the lattice gas automata (LGA) method (Hardy-Pomeau-Pazzis and Frisch-Hasslacher-Pomeau models), is a class of computational fluid dynamics (CFD) methods for fluid simulation. Instead of solving the Navier–Stokes equations directly, a fluid density on a lattice is simulated with streaming and collision (relaxation) processes. The method is versatile as the model fluid can straightforwardly be made to mimic common fluid behaviour like vapour/liquid coexistence, and so fluid systems such as liquid droplets can be simulated. Also, fluids in complex environments such as porous media can be straightforwardly simulated, whereas with complex boundaries other CFD methods can be hard to work with.

<span class="mw-page-title-main">Atmospheric model</span> Mathematical model of atmospheric motions

In atmospheric science, an atmospheric model is a mathematical model constructed around the full set of primitive, dynamical equations which govern atmospheric motions. It can supplement these equations with parameterizations for turbulent diffusion, radiation, moist processes, heat exchange, soil, vegetation, surface water, the kinematic effects of terrain, and convection. Most atmospheric models are numerical, i.e. they discretize equations of motion. They can predict microscale phenomena such as tornadoes and boundary layer eddies, sub-microscale turbulent flow over buildings, as well as synoptic and global flows. The horizontal domain of a model is either global, covering the entire Earth, or regional (limited-area), covering only part of the Earth. Atmospheric models also differ in how they compute vertical fluid motions; some types of models are thermotropic, barotropic, hydrostatic, and non-hydrostatic. These model types are differentiated by their assumptions about the atmosphere, which must balance computational speed with the model's fidelity to the atmosphere it is simulating.

<span class="mw-page-title-main">Lorenz system</span> System of ordinary differential equations with chaotic solutions

The Lorenz system is a system of ordinary differential equations first studied by mathematician and meteorologist Edward Lorenz. It is notable for having chaotic solutions for certain parameter values and initial conditions. In particular, the Lorenz attractor is a set of chaotic solutions of the Lorenz system. The term "butterfly effect" in popular media may stem from the real-world implications of the Lorenz attractor, namely that tiny changes in initial conditions evolve to completely different trajectories. This underscores that chaotic systems can be completely deterministic and yet still be inherently impractical or even impossible to predict over longer periods of time. For example, even the small flap of a butterfly's wings could set the earth's atmosphere on a vastly different trajectory, in which for example a hurricane occurs where it otherwise would have not. The shape of the Lorenz attractor itself, when plotted in phase space, may also be seen to resemble a butterfly.

The material point method (MPM) is a numerical technique used to simulate the behavior of solids, liquids, gases, and any other continuum material. Especially, it is a robust spatial discretization method for simulating multi-phase (solid-fluid-gas) interactions. In the MPM, a continuum body is described by a number of small Lagrangian elements referred to as 'material points'. These material points are surrounded by a background mesh/grid that is used to calculate terms such as the deformation gradient. Unlike other mesh-based methods like the finite element method, finite volume method or finite difference method, the MPM is not a mesh based method and is instead categorized as a meshless/meshfree or continuum-based particle method, examples of which are smoothed particle hydrodynamics and peridynamics. Despite the presence of a background mesh, the MPM does not encounter the drawbacks of mesh-based methods which makes it a promising and powerful tool in computational mechanics.

<span class="mw-page-title-main">Wind wave model</span> Numerical modelling of the sea state

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.

The moving particle semi-implicit (MPS) method is a computational method for the simulation of incompressible free surface flows. It is a macroscopic, deterministic particle method developed by Koshizuka and Oka (1996).

<span class="mw-page-title-main">Sigma coordinate system</span> Coordinate system used in fluid dynamics

The sigma coordinate system is a common coordinate system used in computational models for oceanography, meteorology and other fields where fluid dynamics are relevant. This coordinate system receives its name from the independent variable used to represent a scaled pressure level.

<span class="mw-page-title-main">Gerris (software)</span> Computer Software

Gerris is computer software in the field of computational fluid dynamics (CFD). Gerris was released as free and open-source software, subject to the requirements of the GNU General Public License (GPL), version 2 or any later.

<span class="mw-page-title-main">Regional Ocean Modeling System</span> Free-surface, terrain-following, primitive equations ocean model

Regional Ocean Modeling System (ROMS) is a free-surface, terrain-following, primitive equations ocean model widely used by the scientific community for a diverse range of applications. The model is developed and supported by researchers at the Rutgers University, University of California Los Angeles and contributors worldwide.

Ocean general circulation models (OGCMs) are a particular kind of general circulation model to describe physical and thermodynamical processes in oceans. The oceanic general circulation is defined as the horizontal space scale and time scale larger than mesoscale. They depict oceans using a three-dimensional grid that include active thermodynamics and hence are most directly applicable to climate studies. They are the most advanced tools currently available for simulating the response of the global ocean system to increasing greenhouse gas concentrations. A hierarchy of OGCMs have been developed that include varying degrees of spatial coverage, resolution, geographical realism, process detail, etc.

In mathematics, error analysis is the study of kind and quantity of error, or uncertainty, that may be present in the solution to a problem. This issue is particularly prominent in applied areas such as numerical analysis and statistics.

<span class="mw-page-title-main">River plume</span> Mix of fresh river water and seawater

A river plume is a freshened water mass that is formed in the sea as a result of mixing of river discharge and saline seawater. River plumes are formed in coastal sea areas at many regions in the World. River plumes generally occupy wide-but-shallow sea surface layers bounded by sharp density gradients. The area of a river plume is 3-5 orders of magnitude greater than its depth; therefore, even small rivers with discharge rates ~1–10 m/s form river plumes with horizontal spatial extents ~10–100 m. Areas of river plumes formed by the largest rivers are ~100–1000 km2. Despite the relatively small volume of total freshwater runoff to the World Ocean, river plumes occupy up to 21% of shelf areas of the ocean, i.e., several million square kilometers.

<span class="mw-page-title-main">ADCIRC</span> High-performance, cross-platform numerical ocean circulation model

The ADCIRC model is a high-performance, cross-platform numerical ocean circulation model popular in simulating storm surge, tides, and coastal circulation problems. Originally developed by Drs. Rick Luettich and Joannes Westerink, the model is developed and maintained by a combination of academic, governmental, and corporate partners, including the University of North Carolina at Chapel Hill, the University of Notre Dame, and the US Army Corps of Engineers. The ADCIRC system includes an independent multi-algorithmic wind forecast model and also has advanced coupling capabilities, allowing it to integrate effects from sediment transport, ice, waves, surface runoff, and baroclinicity.

In cosmology, Gurzadyan-Savvidy (GS) relaxation is a theory developed by Vahe Gurzadyan and George Savvidy to explain the relaxation over time of the dynamics of N-body gravitating systems such as star clusters and galaxies. Stellar systems observed in the Universe – globular clusters and elliptical galaxies – reveal their relaxed state reflected in the high degree of regularity of some of their physical characteristics such as surface luminosity, velocity dispersion, geometric shapes, etc. The basic mechanism of relaxation of stellar systems has been considered the 2-body encounters, to lead to the observed fine-grained equilibrium. The coarse-grained phase of evolution of gravitating systems is described by violent relaxation developed by Donald Lynden-Bell. The 2-body mechanism of relaxation is known in plasma physics. The difficulties with description of collective effects in N-body gravitating systems arise due to the long-range character of gravitational interaction, as distinct of plasma where due to two different signs of charges the Debye screening takes place. The 2-body relaxation mechanism e.g. for elliptical galaxies predicts around years i.e. time scales exceeding the age of the Universe. The problem of relaxation and evolution of stellar systems and the role of collective effects are studied by various techniques, see. Among the efficient methods of study of N-body gravitating systems are the numerical simulations, particularly, Sverre Aarseth's N-body codes are widely used.

Tides in marginal seas are tides affected by their location in semi-enclosed areas along the margins of continents and differ from tides in the open oceans. Tides are water level variations caused by the gravitational interaction between the Moon, the Sun and the Earth. The resulting tidal force is a secondary effect of gravity: it is the difference between the actual gravitational force and the centrifugal force. While the centrifugal force is constant across the Earth, the gravitational force is dependent on the distance between the two bodies and is therefore not constant across the Earth. The tidal force is thus the difference between these two forces on each location on the Earth.

<span class="mw-page-title-main">Kuroshio Current Intrusion</span> Movement of water from the Pacific to the West Philippine/South China Sea

The Kuroshio Current is a northward flowing Western Boundary Current (WBC) in the Pacific Ocean. It is a bifurcation arm of the North Equatorial Current and consists of northwestern Pacific Ocean water. The other arm is the southward flowing Mindanao Current. The Kuroshio Current flows along the eastern Philippine coast, up to 13.7 Sv... of it leaking into the Luzon Strait - the gap between the Philippines and Taiwan - before continuing along the Japanese coast. Some of the leaked water manages to intrude into the South China Sea (SCS). This affects the heat and salt budgets and circulation and eddy generation mechanisms in the SCS. There are various theories about possible intrusion paths and what mechanisms initiate them.

Lagrangian ocean analysis is a way of analysing ocean dynamics by computing the trajectories of virtual fluid particles, following the Lagrangian perspective of fluid flow, from a specified velocity field. Often, the Eulerian velocity field used as an input for Lagrangian ocean analysis has been computed using an ocean general circulation model (OGCM). Lagrangian techniques can be employed on a range of scales, from modelling the dispersal of biological matter within the Great Barrier Reef to global scales. Lagrangian ocean analysis has numerous applications, from modelling the diffusion of tracers, through the dispersal of aircraft debris and plastics, to determining the biological connectivity of ocean regions.

OceanParcels, “Probably A Really Computationally Efficient Lagrangian Simulator”, is a set of Python classes and methods that is used to track particles like water, plankton and plastics. It uses the output of ocean general circulation model (OGCMs). OceanParcels main goal is to process the increasingly large amounts of data that is governed by OGCM's. The flow dynamics are simulated using Lagrangian modelling and the geophysical fluid dynamics are simulated with Eulerian modelling or provided through experimental data. OceanParcels is dependent on two principles, namely the ability to read external data sets from different formats and customizable kernels to define particle dynamics.

References

  1. "T h e P r i n c e t o n O c e a n M o d e L".
  2. https://www.pomusers.org/
  3. see POM web
  4. Blumberg, A. F. and G. L. Mellor, A description of a three-dimensional coastal ocean circulation model. Three-Dimensional Coastal ocean Models, edited by N. Heaps, 208 pp., American Geophysical Union., 1987
  5. The Princeton Ocean Model, The Program in Atmospheric and Oceanic Sciences (AOS), Princeton University, archived from the original on July 2, 2010, retrieved November 13, 2010
  6. Zavatarelli, M. and G. L. Mellor, A numerical study of the Mediterranean Sea circulation. J. Phys. Oceanogr., Vol. 25, No. 6, Part II, 1384–1414, 1995
  7. Ezer, T. and G. L. Mellor, Simulations of the Atlantic Ocean with a free surface sigma coordinate ocean model. J. Geophys. Res., 102(C7), 15,647–15,657, 1997.
  8. Mellor, G. L. and T. Ezer, A Gulf Stream model and an altimetry assimilation scheme, J. Geophys. Res., 96, 8779–8795, 1991.
  9. Ezer, T. and G. L. Mellor, A numerical study of the variability and the separation of the Gulf Stream, induced by surface atmospheric forcing and lateral boundary flows. J. Phys. Oceanogr., 22, 660–682, 1992
  10. Aikman, F., G. L. Mellor, T. Ezer, D. Shenin, P. Chen, L. Breaker, and D. B. Rao, Toward an operational nowcast/forecast system for the U.S. east coast, In: Modern Approaches to Data Assimilation in Ocean Modeling, P. Malanotte-Rizzoli Ed., Elsevier Oceanog. Ser., 61, 347–376, 1996.
  11. "T h e P r i n c e t o n O c e a n M o d e L".
  12. "IWMO".
  13. Ocean Dynamics
  14. Oey, L.-Y., Ezer, T., C.-R. Wu and Y. Miyazawa, Editorial- International Workshop on Modeling the Ocean (IWMO) special issue in Ocean Dynamics, Ocean Dynamics, doi:10.1007/s10236-010-0281-7, 60(2), 299–300, 2010.
  15. Oey, L.-Y., Ezer, T., C.-R. Wu and Y. Miyazawa, Editorial- International Workshop on Modeling the Ocean (IWMO) special issue part 2 in Ocean Dynamics, Ocean Dynamics, doi:10.1007/s10236-010-0338-7, 60(5), 1271–1272, 2010.
  16. Ezer, T, L.-Y. Oey, H. Xue and X. H. Wang, Editorial- The 2nd International Workshop on Modeling the Ocean (IWMO-2010), Ocean Dyn., doi:10.1007/s10236-011-0470-z, 61(9),1287–1289, 2011.
  17. Oey, L.-Y., Ezer, T, B. Qiu, J. Berntsen and R. He, Editorial- The 3rd International Workshop on Modeling the Ocean (IWMO-2011), Ocean Dyn., doi:10.1007/s10236-013-0595-3, 2013.
  18. Oey, L.-Y., Y. Miyazawa, N. Aiki, Y. Masumoto, Ezer, T and T. Waseda (2013), Editorial- The 4th International Workshop on Modeling the Ocean (IWMO-2012), Ocean Dyn., 63(11–12), 1345–1347, doi:10.1007/s10236-013-0658-5.
  19. Berntsen, J., L.-Y. Oey, T. Ezer, R. Greatbatch, H. Xue, and Y. Miyazawa (2014), Editorial- The 5th International Workshop on Modeling the Ocean (IWMO-2013), Topical Collection, Ocean Dynamics, doi:10.1007/s10236-014-0764-z
  20. Oey, L.-Y., T. Ezer, J. Sheng, F. Chai, J. Gan, K. Lamb and Y. Miyazawa (2016), Editorial – The 6th International Workshop on Modeling the Ocean (IWMO 2014), Ocean Dynamics, doi:10.1007/s10236-016-1028-x.