User:Liekegeessinck/sandbox

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 (observer moves with particle) and the geophysical fluid dynamics are simulated with Eulerian modelling (observer remains stationary) 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. [1][2][3]

Concept

edit

OceanParcels takes a flow field given through Eulerian methods or experimental data (e.g. provided through land-based measurements, satellite imagery or Doppler currents) as input. This data can be given on various grid structures. To obtain data at an arbitrary point, interpolation is used. Lagrangian modelling uses the data given by the flow field and interpolation to simulate the dynamics of an object. [1]

Grid Structure

edit
 
Figure 1: Grid structures in horizontal and vertical direction. [2]
 
Figure 2: left: A grid, middle: B grid, right: C grid[4]

In the horizontal plane, i.e. in the x-y-plane, OceanParcels allows for rectilinear and curvilinear data grids (see Fig. 1 a and b respectively). In the vertical plane, i.e. the x-z-plane, z-levels and s-levels are supported (see Fig. 1 c and d respectively). All four combinations of these grids are supported in three-dimensional space.[2]

OceanParcels also supports so-called staggered A, B and C grids, which are common in ocean modelling. These take into account that the different variables can be measured at different grid points. The relevant variables in Ocean analysis are the meriodial velocity (v), the zonal velocity (u) and the Tracers (q). In an A grid all of these are evaluated at the same grid points. In a B grid, u and v are evaluated at the grid point in the center between the grid point where T is being evaluated, that is there are two overlying grids. In a C grid all variables are evaluated on their own grid, that is there are three overlying grids. See Figure 2 for a visualization.[2]

Interpolation Methods

edit

To obtain field data at a point that is not on a grid node, interpolation methods are used. OceanParcels supports several different such methods:

Lagrangian Modelling

edit

The Lagrangian trajectory of a particle can be calculated by numerically integrating over its position   as a function of time, given by  . Here   is the starting time and   is the time at which the position of the particle is evaluated.   is the velocity of the particle at time  .   captures the movement of the particle caused by its specific characteristics (e.g. an iceberg moves differently in the ocean than a non-frozen water mass). Unless otherwise specified, OceanParcels uses a fourth-order Runge-Kutta scheme for the integration of this function. As alternatives, OceanParcels supports Euler-forward integration or adaptive Runge-Kutta-Fehlberg integration. [1]

Kernels

edit

OceanParcels provides the option to model individual behavior of a particle. Individual behavior is behavior that is defined by characteristics of a particle, for example fish has some control over its movement as long it is alive while plastic has not, leading to different dynamics. OceanParcels implements this using short pieces of code that are executed whenever a ParticleSet (class defining the particles used) is executed. Those pieces of code are called kernels. Some commonly used kernels are predefined in OceanParcels, such as kernels modeling Brownian motion. However, the user may also define own problem specific kernels. This allows the user to model various objects, examples of which are given in the section "Applications". For implementing customized kernels the user can choose between using Python with a possibility to automatically translate to C or using a C-library. [1][2][3]

Applications

edit

Lagrangian analysis can be used to simulate the pathways of virtual particles, which can represent water masses, temperature tracers, salinity tracers, nutrient tracers, kelp, coral larvae, plastics, fish, icebergs, plankton, etc. The pathways of the different particles can be modelled using specific kernels, for example that fish can swim and that icebergs can melt. [1]

Water masses

edit

The movement of water masses can be calculated numerically using velocity fields, which are simulated by a primitive equation model (fine resolution Antarctic model). Using this model, it can be estimated how many times the water mass has traveled around Antarctica and how much time it has spend in the Southern Ocean. It turned out that a water mass travels around Antarctica six times, on average, before it is ventilated. A particle is ventilated when it reaches the surfaces at least once. [5]

Surface temperatures can be analyzed using proxies, which can reconstruct environmental changes over the last 600 years. These reconstructions are based on alkenones, Membrane lipids and stable isotope ratios of globigerinoides ruber. [6] It is argued that ocean currents are transporting these proxies far away from there origin, creating a bias towards a too high or low temperature. Studies have shown that these currents have a negligible effect on the recorded temperature. [7]

The proxies used for salinity are the ratios between elements and calcium and stable oxygen isotopes ratios of foraminifera. These isotope ratios can be calculated by the following formula:  , where y is the common isotope and x the more rare isotope. Both ratios have a positive correlation in the Mediterranean Sea, which has a strong west to east salinity gradient. Salinity tracers are important for reconstructing the ocean circulation, because together with temperature, the density can be determined, which can reconstruct large-scale circulation patterns, including meridional overturning circulation.[8]

The influence of the Equatorial Undercurrent (EUC), New Guinea Coastal Undercurrent (NGCU) and the New Ireland Coastal Undercurrent were examined using models to elucidate the high iron concentrations in the Pacific ocean. During El Niño, the latter two currents strengthen, which enhances the iron concentrations in the Pacific. [9]

Plastics are harmful for marine life, fisheries and tourism. Lagrangian particle models have shown that about three quarters of buoyant marine plastics released from land are ending up in coastal waters, with the highest concentration in Southeast Asia. However, these simulations are overestimated with respect to field measurements. [10]

The transport of fish larvae by ocean currents is an important dispersal mechanism, because the timing has a lot of influence on the location where the larvae are settled. Pomatomus saltatrix is a fish species that is present globally, and its dispersal is examined by using particle tracking simulations. Especially the East Australian Current is important for maximum dispersal of the larvae. [11]

Plankton are important archives to reconstruct past conditions of the ocean surfaces. While the plankton sinks to the sediments, it can be advected by turbulent ocean currents. Ocean simulations can be used to define ocean bottom provinces based on the particle surface origin locations. [12]

The influence of icebergs on sea ice can be modeled and compared to observations. When icebergs melt, a large fresh water flux takes place. Around Antarctica, a negative fresh water flux takes place due to sea ice freezing.[13]

Features

edit

Development

edit

OceanParcels was first introduced by Erik van Sebille and Michael Lange in 2017. The latest version can accessed via GitHub. [3] It is constantly being improved, updated and expanded. The code design ensures that this can be done without effecting the user interface. [1] Below you see a list of some possible improvements:

  • Implementing support for using Arakawa D and E grids in the interpolation[2]
  • Implementing more advanced interpolation methods[1]
  • Implementing support for using unstructured grids in the interpolation process[2]
  • Implementing new more efficient programming methods[1]
  • Implementing further numerical integration methods for the Lagrangian modelling[2]

Funding

edit

OceanParcels is being funded by the following means:

Alternatives to OceanParcels

edit

The following list contains some code packages with (partially) similar functions as OceanParcels.

References

edit
  1. ^ a b c d e f g h i Parcels v0.9: prototyping a Lagrangian Ocean Analysis framework for the petascale age Lange, M and E van Sebille (2017), Geoscientific Model Development, 10, 4175-4186
  2. ^ a b c d e f g h i Delandmeter, Philippe; van Sebille, Erik (2019) [19.08.2019]. "The Parcels v2.0 Lagrangian framework: new field interpolation schemes". Geoscientific Model Development. 12 (8): 3571–3584.
  3. ^ a b c d "OceanParcels". Retrieved 6 April 2022.
  4. ^ "Ocean general circulation model", Wikipedia, 2021-12-14, retrieved 2022-04-05
  5. ^ Döös, K.: Interocean Exchange of Water Masses, J. Geophys. Res.- Oceans, 100, 13499–13514, 1995.
  6. ^ Grauel, A. L., Leider, A., Goudeau, M. L. S., Müller, I. A., Bernasconi, S. M., Hinrichs, K. U., et al. (2013). What do SST proxies really tell us? A high-resolution multiproxy (UK'37, TEXH86 and foraminifera δ18O) study in the Gulf of Taranto, central Mediterranean Sea. Quaternary Science Reviews, 73, 115–131. https://doi.org/10.1016/j.quascirev.2013.05.007
  7. ^ Limited lateral transport bias during export of sea surface temperature proxy carriers in the Mediterranean Sea, Rice, A, PB Nooteboom, E van Sebille, F Peterse, M Ziegler, A Sluijs (2022), Geophysical Research Letters, 49, e2021GL096859.
  8. ^ Dämmer, L. K., de Nooijer, L., van Sebille, E., Haak, J. G., and Reichart, G.-J.: Evaluation of oxygen isotopes and trace elements in planktonic foraminifera from the Mediterranean Sea as recorders of seawater oxygen isotopes and salinity, Clim. Past, 16, 2401–2414, https://doi.org/10.5194/cp-16-2401-2020, 2020.
  9. ^ Qin, X., Menviel, L., Sen Gupta, A., and van Sebille, E.: Iron sources and pathways into the Pacific Equatorial Undercurrent, Geophys. Res. Lett., 43, 9843–9851, https://doi.org/10.1002/2016GL070501, 2016.
  10. ^ Global simulations of marine plastic transport show plastic trapping in coastal zones Onink, V, C Jongedijk, MJ Hoffman, E van Sebille, C Laufkötter (2021), Environmental Research Letters, 16, 064053
  11. ^ Multiple spawning events promote increased larval dispersal of a predatory fish in a western boundary current, Schilling, HT, JD Everett, JA Smith, J Stewart, JM Hughes, M Roughan, C Kerry, and IM Suthers (2020), Fisheries Oceanography, 29, 309-323
  12. ^ Nooteboom, P. D., Bijl, P. K., Kehl, C., van Sebille, E., Ziegler, M., von der Heydt, A. S., and Dijkstra, H. A.: Sedimentary microplankton distributions are shaped by oceanographically connected areas, Earth Syst. Dynam., 13, 357–371, https://doi.org/10.5194/esd-13-357-2022, 2022.
  13. ^ Marsh, R., Ivchenko, V. O., Skliris, N., Alderson, S., Bigg, G. R., Madec, G., Blaker, A. T., Aksenov, Y., Sinha, B., Coward, A. C., Le Sommer, J., Merino, N., and Zalesny, V. B.: NEMO-ICB (v1.0): interactive icebergs in the NEMO ocean model globally configured at eddy-permitting resolution, Geosci. Model Dev., 8, 1547–1562, https://doi.org/10.5194/gmd-8-1547-2015, 2015.
  14. ^ Dagestad, Knut-Frode; Röhrs, Johannes; Breivik, Oyvind; Adlandsvik, Bjorn (13 April 2018). "OpenDrift v1.0: a generic framework for trajectory modelling" (PDF). Geoscientific Model Development. 11 (4): 1405–1420.
  15. ^ "TRACMASS". Retrieved 6 April 2022.