Regional Ocean Modeling System

Last updated
ROMS logo.png

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.

Contents

ROMS is used to model how a given region of the ocean responds to physical forcings such as heating or wind. It can also be used to model how a given ocean system responds to inputs like sediment, freshwater, ice, or nutrients, requiring coupled models nested within the ROMS framework.

Framework

ROMS is a 4D modeling system. It is a 3-dimensional model (a 2D horizontal grid and a vertical grid) that can be run over a given amount of time, time being the 4th dimension. It is gridded into vertical levels that make up the water column and horizontal cells that make up the coordinates of the 2D cartesian plane of the model region.

Kernel

Central to the ROMS framework are four models that form what is called the dynamical/numerical core or kernel:

  1. Non-Linear Model kernel (NLM): NLROMS [1] [2]
  2. Perturbation Tangent Linear Model kernel (TLM): TLROMS
  3. Finite-amplitude tangent linear Representer Model kernel (RPM): RPROMS
  4. Adjoint Model kernel (ADM): ADROMS [3]

Vertical grid

The vertical grid is a hybrid stretched grid. It is hybrid in that its stretching intervals fall somewhere between the two extremes of 1) the evenly-spaced sigma grid used by the Princeton Ocean Model and 2) a true z-grid with a static depth interval . The vertical grid can be squeezed or stretched to increase or decrease the resolution for an area of interest, such as a thermocline or bottom boundary layer. Grid stretching in the vertical direction follows bottom topography, allowing for the idealized flow of water over features such as seamounts. [4] The numbering of the vertical grid goes from the bottom waters upward to the air-water interface: the bottom water level is level 1 and the topmost surface water level is the highest number (such as level 20). With a coupled sediment module, the numbering of the sediment seabed levels goes from the sediment-water interface downward: the topmost seabed level is level 1 and the deepest seabed level is the highest number.

Horizontal grid

The horizontal grid is a structured grid, meaning that it has a rectangular 4-sided grid cell structure. The horizontal grid is also an orthogonal curvilinear grid, meaning that it maximizes ocean grid cells of interest and minimizes extra land grid cells. The horizontal grid is also a staggered grid or Arakawa-C grid, where the velocities in the north-south and east-west directions are calculated at the edges of each grid cell, while the values for scalar variables such as density are calculated at the center of each grid cell, known as "rho-points."

Physics

In both the vertical and horizontal directions, the default equations use centered, second-order finite difference schemes. Higher order schemes are available if desired, for example using parabolic spline reconstruction. [2]

In general, the physical schemes used by ROMS are based on three governing equations:

  1. Continuity
  2. Conservation of momentum (Navier-Stokes)
  3. Transport equations of tracer variables (such as salinity and temperature)

Equations are coupled to solve for five unknowns at each location in the model grid using numerical solutions:

Source code

ROMS uses an open-access source code that can be downloaded by filling out an online request form. It runs on C-processing and was developed for shared computing uses. To download the source code a user must create an account and file a request with the developers on the ROMS website.

Input and output

Input

Boundaries such as coastlines can be specified for a given region using land- and sea-masking. The top vertical boundary, the air-sea interface, uses an interaction scheme developed by Fairall et al. (1996). [5] The bottom vertical boundary, the sediment-water interface, uses a bottom stress or bottom-boundary-layer scheme developed by Styles and Glenn (2000). [6]

Inputs that are needed for an implementer to run ROMS for a specific ocean region include:

The programming framework of ROMS is split into three parts: Initialize, Run, and Finalize, which is standard for the Earth System Modeling Framework (ESMF). "Run" is the largest of these three parts, where the user chooses which options they want to use and assimilates data if desired. [7] The model run must be initialized or compiled before it is run.

Output

The output format of model run files is netCDF. Model output is often visualized using independent secondary programming software such as MATLAB or Python. Simple visualization software such as NASA's Panoply Data Viewer can also be used to visualize model output for teaching or demonstration purposes.

User options

The general approach of ROMS gives model implementers a high level of freedom and responsibility. One approach cannot meet the needs of all the diverse applications the model is currently used for. Therefore, it is up to each model implementer (either an individual or a research group) to choose how they want to use each of the available options. Options include choices such as:

When using ROMS, if an implementer runs into a problem or bug, they can report it to the ROMS forum.

Applications

A daily, global Sea Surface Temperature (SST) data set from December 2013 produced at 1-km resolution (also known as ultra-high resolution) by the JPL ROMS group. SST 20131220 blended Global.png
A daily, global Sea Surface Temperature (SST) data set from December 2013 produced at 1-km resolution (also known as ultra-high resolution) by the JPL ROMS group.

The versatility of ROMS has been proven in its diverse applications to different systems and regions. It is best applied to mesoscale systems, [9] or those systems that can be mapped at high resolution, such as 1-km to 100-km grid spacing.

Coupled model applications

Biogeochemical, bio-optical, sea ice, sediment, and other models can be embedded within the ROMS framework to study specific processes. These are usually developed for specific regions of the world's oceans but can be applied elsewhere. For example, the sea ice application of ROMS was originally developed for the Barents Sea Region. [10]

ROMS modeling efforts are increasingly being coupled with observational platforms, such as buoys, satellites, and ship-mounted underway sampling systems, to provide more accurate forecasting of ocean conditions.

Regional applications

There is an ever-growing number of applications of ROMS to particular regions of the world's oceans. These integrated ocean modeling systems use ROMS for the circulation component, and add other variables and processes of interest. A few examples are:

See also

Related Research Articles

Climate model Quantitative methods used to simulate climate

Numerical climate models use quantitative methods to simulate the interactions of the important drivers of climate, including atmosphere, oceans, land surface and ice. They are used for a variety of purposes from study of the dynamics of the climate system to projections of future climate. Climate models may also be qualitative models and also narratives, largely descriptive, of possible futures.

General circulation model Type of climate model

A general circulation model (GCM) is a type of climate model. It employs a mathematical model of the general circulation of a planetary atmosphere or ocean. It uses the Navier–Stokes equations on a rotating sphere with thermodynamic terms for various energy sources. These equations are the basis for computer programs used to simulate the Earth's atmosphere or oceans. Atmospheric and oceanic GCMs are key components along with sea ice and land-surface components.

The oxygen minimum zone (OMZ), sometimes referred to as the shadow zone, is the zone in which oxygen saturation in seawater in the ocean is at its lowest. This zone occurs at depths of about 200 to 1,500 m (660–4,920 ft), depending on local circumstances. OMZs are found worldwide, typically along the western coast of continents, in areas where an interplay of physical and biological processes concurrently lower the oxygen concentration and restrict the water from mixing with surrounding waters, creating a "pool" of water where oxygen concentrations fall from the normal range of 4–6 mg/l to below 2 mg/l.

Numerical weather prediction Weather prediction using mathematical models of the atmosphere and oceans

Numerical weather prediction (NWP) uses mathematical models of the atmosphere and oceans to predict the weather based on current weather conditions. Though first attempted in the 1920s, it was not until the advent of computer simulation in the 1950s that numerical weather predictions produced realistic results. A number of global and regional forecast models are run in different countries worldwide, using current weather observations relayed from radiosondes, weather satellites and other observing systems as inputs.

Sea spray Sea water particles that are formed directly from the ocean

Sea spray are aerosol particles formed from the ocean, mostly by ejection into Earth's atmosphere by bursting bubbles at the air-sea interface. Sea spray contains both organic matter and inorganic salts that form sea salt aerosol (SSA). SSA has the ability to form cloud condensation nuclei (CCN) and remove anthropogenic aerosol pollutants from the atmosphere. Also, coarse sea spray inhibits lightning.

Eddy (fluid dynamics) Swirling of a fluid and the reverse current created when the fluid is in a turbulent flow regime

In fluid dynamics, an eddy is the swirling of a fluid and the reverse current created when the fluid is in a turbulent flow regime. The moving fluid creates a space devoid of downstream-flowing fluid on the downstream side of the object. Fluid behind the obstacle flows into the void creating a swirl of fluid on each edge of the obstacle, followed by a short reverse flow of fluid behind the obstacle flowing upstream, toward the back of the obstacle. This phenomenon is naturally observed behind large emergent rocks in swift-flowing rivers.

Atmospheric model

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.

Thin layers (oceanography) Congregations of plankton

Thin layers are concentrated aggregations of phytoplankton and zooplankton in coastal and offshore waters that are vertically compressed to thicknesses ranging from several centimeters up to a few meters and are horizontally extensive, sometimes for kilometers. Generally, thin layers have three basic criteria: 1) they must be horizontally and temporally persistent; 2) they must not exceed a critical threshold of vertical thickness; and 3) they must exceed a critical threshold of maximum concentration. The precise values for critical thresholds of thin layers has been debated for a long time due to the vast diversity of plankton, instrumentation, and environmental conditions. Thin layers have distinct biological, chemical, optical, and acoustical signatures which are difficult to measure with traditional sampling techniques such as nets and bottles. However, there has been a surge in studies of thin layers within the past two decades due to major advances in technology and instrumentation. Phytoplankton are often measured by optical instruments that can detect fluorescence such as LIDAR, and zooplankton are often measured by acoustic instruments that can detect acoustic backscattering such as ABS. These extraordinary concentrations of plankton have important implications for many aspects of marine ecology, as well as for ocean optics and acoustics. Zooplankton thin layers are often found slightly under phytoplankton layers because many feed on them. Thin layers occur in a wide variety of ocean environments, including estuaries, coastal shelves, fjords, bays, and the open ocean, and they are often associated with some form of vertical structure in the water column, such as pycnoclines, and in zones of reduced flow.

The MEMO model is a Eulerian non-hydrostatic prognostic mesoscale model for wind-flow simulation. It was developed by the Aristotle University of Thessaloniki in collaboration with the Universität Karlsruhe. The MEMO Model together with the photochemical dispersion model MARS are the two core models of the European zooming model (EZM). This model belongs to the family of models designed for describing atmospheric transport phenomena in the local-to-regional scale, frequently referred to as mesoscale air pollution models.

James C. McWilliams is a professor at the UCLA Institute of Geophysics and Planetary Physics and Department of Atmospheric and Oceanic Sciences.

In oceanography, a front is a boundary between two distinct water masses. The formation of fronts depends on multiple physical processes and small differences in these lead to a wide range of front types. They can be as narrow as a few hundreds of metres and as wide as several tens of kilometres. While most fronts form and dissipate relatively quickly, some can persist for long periods of time.

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.

River plume

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 layer bounded by sharp density gradient. 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 World rivers are ~100–1000 km2. Despite relatively small volume of total freshwater runoff to the World Ocean, river plumes occupy up to 21% of shelf areas of the World Ocean, i.e., several million square kilometers.

ADCIRC 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.

CICE is a computer model that simulates the growth, melt and movement of sea ice. It has been integrated into many coupled climate system models as well as global ocean and weather forecasting models and is often used as a tool in Arctic and Southern Ocean research. CICE development began in the mid-1990s by the United States Department of Energy (DOE), and it is currently maintained and developed by a group of institutions in North America and Europe known as the CICE Consortium. Its widespread use in earth system science in part owes to the importance of sea ice in determining Earth's planetary albedo, the strength of the global thermohaline circulation in the world's oceans, and in providing surface boundary conditions for atmospheric circulation models, since sea ice occupies a significant proportion (4-6%) of earth's surface. CICE is a type of cryospheric model.

A baroclinic instability is a fluid dynamical instability of fundamental importance in the atmosphere and ocean. It can lead to the formation of transient mesoscale eddies, with a horizontal scale of 10-100 km. In contrast, flows on the largest scale in the ocean are described as ocean currents, the largest scale eddies are mostly created by shearing of two ocean currents and static mesoscale eddies are formed by the flow around an obstacle (as seen in the animation on eddy. Mesoscale eddies are circular currents with swirling motion and account for approximately 90% of the ocean's total kinetic energy. Therefore, they are key in mixing and transport of for example heat, salt and nutrients.

Lisa M. Beal is a professor at the University of Miami known for her work on the Agulhas Current. She is the editor-in-chief of the Journal of Geophysical Research: Oceans.

North Brazil Current Retroflection Ocean current retroflection phenomenon

The North Brazil Current (NBC) retroflects north-eastwards and merges into the North Equatorial Counter Current (NECC). The retroflection occurs in a seasonal pattern when there is strong retroflection from late summer to early winter. There is weakened or no retroflection during other times of the year. Just like in the Agulhas Current, the retroflection also sheds some eddies that make their way to the Caribbean Sea through the Lesser Antilles.

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 Models (OGCM's). 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. Shchepetkin, Alexander F. (2003). "A method for computing horizontal pressure-gradient force in an oceanic model with a nonaligned vertical coordinate". Journal of Geophysical Research. 108 (C3). doi: 10.1029/2001jc001047 . ISSN   0148-0227.
  2. 1 2 Shchepetkin, A.F.; McWilliams, J.C. (2005). The Regional Ocean Modeling System: A Split-Explicit, Free-Surface, Topography-Following-Coordinate Ocean Model, 2003. Los Angeles, California: University of California at Los Angeles: Institute of Geophysics and Planetary Physics.
  3. Moore, Andrew M.; Arango, Hernan G.; Di Lorenzo, Emanuele; Cornuelle, Bruce D.; Miller, Arthur J.; Neilson, Douglas J. (2004-01-01). "A comprehensive ocean prediction and analysis system based on the tangent linear and adjoint of a regional ocean model". Ocean Modelling. 7 (1–2): 227–258. doi:10.1016/j.ocemod.2003.11.001. ISSN   1463-5003.
  4. Song, Yuhe; Haidvogel, Dale (1994-11-01). "A Semi-implicit Ocean Circulation Model Using a Generalized Topography-Following Coordinate System". Journal of Computational Physics. 115 (1): 228–244. doi:10.1006/jcph.1994.1189. ISSN   0021-9991.
  5. Fairall, C. W.; Bradley, E. F.; Rogers, D. P.; Edson, J. B.; Young, G. S. (1996-02-15). "Bulk parameterization of air-sea fluxes for Tropical Ocean-Global Atmosphere Coupled-Ocean Atmosphere Response Experiment". Journal of Geophysical Research: Oceans. 101 (C2): 3747–3764. CiteSeerX   10.1.1.469.6689 . doi:10.1029/95jc03205. ISSN   0148-0227.
  6. Styles, Richard; Glenn, Scott M. (2000-10-15). "Modeling stratified wave and current bottom boundary layers on the continental shelf" (PDF). Journal of Geophysical Research: Oceans. 105 (C10): 24119–24139. doi: 10.1029/2000jc900115 . ISSN   0148-0227. Archived from the original (PDF) on 2019-12-30.
  7. "ROMS > start". www.myroms.org. Retrieved 2019-02-08.
  8. Hedstrom, Katherine S. (2016). "Technical Manual for a Coupled Sea-Ice/Ocean Circulation Model (Version 5)" (PDF). OCS Study BOEM 2016-037. Cooperative Agreement No. M15AC00011.
  9. "Met Office: Mesoscale modelling". 2010-12-29. Archived from the original on 2010-12-29. Retrieved 2018-04-26.
  10. Budgell, W. P. (2005-12-01). "Numerical simulation of ice-ocean variability in the Barents Sea region". Ocean Dynamics. 55 (3–4): 370–387. doi:10.1007/s10236-005-0008-3. ISSN   1616-7341.
  11. Warner, John C.; Armstrong, Brandy; He, Ruoying; Zambon, Joseph B. (2010-01-01). "Development of a Coupled Ocean–Atmosphere–Wave–Sediment Transport (COAWST) Modeling System" (PDF). Ocean Modelling. 35 (3): 230–244. doi:10.1016/j.ocemod.2010.07.010. hdl: 1912/4099 . ISSN   1463-5003.
  12. Feng, Yang; Friedrichs, Marjorie A. M.; Wilkin, John; Tian, Hanqin; Yang, Qichun; Hofmann, Eileen E.; Wiggert, Jerry D.; Hood, Raleigh R. (2015). "Chesapeake Bay nitrogen fluxes derived from a land-estuarine ocean biogeochemical modeling system: Model description, evaluation, and nitrogen budgets". Journal of Geophysical Research: Biogeosciences. 120 (8): 1666–1695. doi:10.1002/2015jg002931. PMC   5014239 . PMID   27668137.
  13. Combes, Vincent; Di Lorenzo, Emanuele (2007-10-01). "Intrinsic and forced interannual variability of the Gulf of Alaska mesoscale circulation". Progress in Oceanography. 75 (2): 266–286. doi:10.1016/j.pocean.2007.08.011. hdl: 1853/14532 . ISSN   0079-6611.
  14. "SOCIB System Description". www.socib.es. Retrieved 2022-08-14.