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.

<span class="mw-page-title-main">Loop Current</span> Ocean current between Cuba and Yucatán Peninsula

A parent to the Florida Current, the Loop Current is a warm ocean current that flows northward between Cuba and the Yucatán Peninsula, moves north into the Gulf of Mexico, loops east and south before exiting to the east through the Florida Straits and joining the Gulf Stream. The Loop Current is an extension of the western boundary current of the North Atlantic subtropical gyre. Serving as the dominant circulation feature in the Eastern Gulf of Mexico, the Loop Currents transports between 23 and 27 sverdrups and reaches maximum flow speeds of from 1.5 to 1.8 meters/second.

<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. The different types of models run are thermotropic, barotropic, hydrostatic, and nonhydrostatic. Some of the model types make assumptions about the atmosphere which lengthens the time steps used and increases computational speed.

<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 several different initial chaotic conditions evolve in phase space in a way that never repeats, so all chaos is unpredictable. This underscores that chaotic systems can be completely deterministic and yet still be inherently unpredictable over long periods of time. Because chaos continually increases in systems, it is impossible to predict the future of systems well. For instance, even the small flap of a butterfly's wings could set the world on a vastly different trajectory, such as by causing a hurricane. 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.

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.

<span class="mw-page-title-main">HyCOM</span>

The Hybrid Coordinate Ocean Model (HyCOM) is an open-source ocean general circulation modeling system. HyCOM is a primitive equation type of ocean general circulation model. The vertical levels of this modeling system are slightly different than other models, because the vertical coordinates remain isopycnic in the open stratified ocean, smoothly transitioning to z-level coordinates in the weakly stratified upper-ocean mixed layer, to terrain-following sigma coordinates in shallow water regions, and back to z-level coordinates in very shallow water. Therefore, the setup is a “hybrid” between z-level and terrain-following vertical levels. HyCOM outputs are provided online for the global ocean at a spatial resolution of 0.08 degrees (approximately 9 km) from 2003 to present. HyCOM uses netCDF data format for model outputs.

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.

<span class="mw-page-title-main">Ronald J. Stouffer</span> American climate scientist

Ronald J. Stouffer is a meteorologist and adjunct professor at the University of Arizona, formerly Senior Research Climatologist and head of the Climate and Ecosystems Group at the Geophysical Fluid Dynamics Laboratory (GFDL), part of NOAA. He has also served on the faculty of Princeton University.

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. http://www.ccpo.odu.edu/~tezer/POMDB/references.php
  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.