Articles | Volume 13, issue 12
Geosci. Model Dev., 13, 5935–5957, 2020
https://doi.org/10.5194/gmd-13-5935-2020
© Author(s) 2020. This work is distributed under
the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
Geosci. Model Dev., 13, 5935–5957, 2020
https://doi.org/10.5194/gmd-13-5935-2020
© Author(s) 2020. This work is distributed under
the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
Development and technical paper 01 Dec 2020
Development and technical paper | 01 Dec 2020
Numerical integrators for Lagrangian oceanography
Tor Nordam and Rodrigo Duran
Related authors
Johannes Röhrs, Knut-Frode Dagestad, Helene Asbjørnsen, Tor Nordam, Jørgen Skancke, Cathleen E. Jones, and Camilla Brekke
Ocean Sci., 14, 1581–1601, https://doi.org/10.5194/os-14-1581-2018, https://doi.org/10.5194/os-14-1581-2018, 2018
Short summary
Short summary
Simulations of hypothetical oil spills are presented to investigate how the vertical mixing of oil affects transport towards various directions. It is shown that the horizontal transport of oil greatly varies for different oil types and weather conditions. These differences are a consequence of the entrainment of oil from the surface into the ocean. While oil spills often get entrained into the water by waves, we show that submerged oil typically resurfaces after a few hours or days.
Johannes Röhrs, Knut-Frode Dagestad, Helene Asbjørnsen, Tor Nordam, Jørgen Skancke, Cathleen E. Jones, and Camilla Brekke
Ocean Sci., 14, 1581–1601, https://doi.org/10.5194/os-14-1581-2018, https://doi.org/10.5194/os-14-1581-2018, 2018
Short summary
Short summary
Simulations of hypothetical oil spills are presented to investigate how the vertical mixing of oil affects transport towards various directions. It is shown that the horizontal transport of oil greatly varies for different oil types and weather conditions. These differences are a consequence of the entrainment of oil from the surface into the ocean. While oil spills often get entrained into the water by waves, we show that submerged oil typically resurfaces after a few hours or days.
Related subject area
Oceanography
The Regional Ice Ocean Prediction System v2: a pan-Canadian ocean analysis system using an online tidal harmonic analysis
Global storm tide modeling with ADCIRC v55: unstructured mesh design and performance
Development of a MetUM (v 11.1) and NEMO (v 3.6) coupled operational forecast model for the Maritime Continent – Part 1: Evaluation of ocean forecasts
Advanced parallel implementation of the coupled ocean–ice model FEMAO (version 2.0) with load balancing
The Meridionally Averaged Model of Eastern Boundary Upwelling Systems (MAMEBUSv1.0)
Model-driven optimization of coastal sea observatories through data assimilation in a finite element hydrodynamic model (SHYFEM v. 7_5_65)
A simplified atmospheric boundary layer model for an improved representation of air–sea interactions in eddying oceanic models: implementation and first evaluation in NEMO (4.0)
Performance of offline passive tracer advection in the Regional Ocean Modeling System (ROMS; v3.6, revision 904)
Implementation and assessment of a carbonate system model (Eco3M-CarbOx v1.1) in a highly dynamic Mediterranean coastal site (Bay of Marseille, France)
Towards Multiscale Modeling of Ocean Surface Turbulent Mixing Using Coupled MPAS-Ocean v6.3 and PALM v5.0
Multi-grid algorithm for passive tracer transport in the NEMO ocean circulation model: a case study with the NEMO OGCM (version 3.6)
HIDRA 1.0: Deep-Learning-Based Ensemble Sea Level Forecasting in the Northern Adriatic
Improved representation of river runoff in Estimating the Circulation and Climate of the Ocean Version 4 (ECCOv4) simulations: implementation, evaluation and impacts to coastal plume regions
Introducing LAB60: A 1∕60° NEMO 3.6 numerical simulation of the Labrador Sea
Impact of horizontal resolution on global ocean–sea ice model simulations based on the experimental protocols of the Ocean Model Intercomparison Project phase 2 (OMIP-2)
CSIRO Environmental Modelling Suite (EMS): scientific description of the optical and biogeochemical models (vB3p0)
Constraining the response of phytoplankton to zooplankton grazing and photo-acclimation in a temperate shelf sea with a 1-D model – towards S2P3 v8.0
Comparison of ocean vertical mixing schemes in the Max Plank Institute Earth System Model (MPI-ESM1.2)
Doppio – a ROMS (v3.6)-based circulation model for the Mid-Atlantic Bight and Gulf of Maine: configuration and comparison to integrated coastal observing network observations
Evaluation of global ocean–sea-ice model simulations based on the experimental protocols of the Ocean Model Intercomparison Project phase 2 (OMIP-2)
Representation of the Denmark Strait overflow in a z-coordinate eddying configuration of the NEMO (v3.6) ocean model: resolution and parameter impacts
A global eddying hindcast ocean simulation with OFES2
BFM17 v1.0: Reduced-Order Biogeochemical Flux Model for Upper Ocean Biophysical Simulations
Tracking water masses using passive-tracer transport in NEMO v3.4 with NEMOTAM: application to North Atlantic Deep Water and North Atlantic Subtropical Mode Water
Sensitivity study on the main tidal constituents of the Gulf of Tonkin by using the frequency-domain tidal solver in T-UGOm
DINCAE 1.0: a convolutional neural network with error estimates to reconstruct sea surface temperature satellite observations
Mitigation of model bias influences on wave data assimilation with multiple assimilation systems using WaveWatch III v5.16 and SWAN v41.20
MOMSO 1.0 – an eddying Southern Ocean model configuration with fairly equilibrated natural carbon
Simulating barrier island response to sea level rise with the barrier island and inlet environment (BRIE) model v1.0
Dealing with discontinuous meteorological forcing in operational ocean modelling: a case study using ECMWF-IFS and GETM (v2.5)
The Parcels v2.0 Lagrangian framework: new field interpolation schemes
The INALT family – a set of high-resolution nests for the Agulhas Current system within global NEMO ocean/sea-ice configurations
Sensitivity of deep ocean biases to horizontal resolution in prototype CMIP6 simulations with AWI-CM1.0
OceanMesh2D 1.0: MATLAB-based software for two-dimensional unstructured mesh generation in coastal ocean modeling
Ocean carbon and nitrogen isotopes in CSIRO Mk3L-COAL version 1.0: a tool for palaeoceanographic research
A high-resolution biogeochemical model (ROMS 3.4 + bio_Fennel) of the East Australian Current system
Nemo-Nordic 1.0: a NEMO-based ocean model for the Baltic and North seas – research and operational applications
Ecological ReGional Ocean Model with vertically resolved sediments (ERGOM SED 1.0): coupling benthic and pelagic biogeochemistry of the south-western Baltic Sea
Reanalysis of the PacIOOS Hawaiian Island Ocean Forecast System, an implementation of the Regional Ocean Modeling System v3.6
Thetis coastal ocean model: discontinuous Galerkin discretization for the three-dimensional hydrostatic equations
Data assimilation cycle length and observation impact in mesoscale ocean forecasting
Verification of the mixed layer depth in the OceanMAPS operational forecast model for Austral autumn
A global scavenging and circulation ocean model of thorium-230 and protactinium-231 with improved particle dynamics (NEMO–ProThorP 0.1)
Veros v0.1 – a fast and versatile ocean simulator in pure Python
Cohesive and mixed sediment in the Regional Ocean Modeling System (ROMS v3.6) implemented in the Coupled Ocean–Atmosphere–Wave–Sediment Transport Modeling System (COAWST r1234)
OpenDrift v1.0: a generic framework for trajectory modelling
A 4.5 km resolution Arctic Ocean simulation with the global multi-resolution model FESOM 1.4
AMM15: a new high-resolution NEMO configuration for operational simulation of the European north-west shelf
SedFoam-2.0: a 3-D two-phase flow numerical model for sediment transport
Parcels v0.9: prototyping a Lagrangian ocean analysis framework for the petascale age
Gregory C. Smith, Yimin Liu, Mounir Benkiran, Kamel Chikhar, Dorina Surcel Colan, Audrey-Anne Gauthier, Charles-Emmanuel Testut, Frederic Dupont, Ji Lei, François Roy, Jean-François Lemieux, and Fraser Davidson
Geosci. Model Dev., 14, 1445–1467, https://doi.org/10.5194/gmd-14-1445-2021, https://doi.org/10.5194/gmd-14-1445-2021, 2021
Short summary
Short summary
Canada's coastlines include diverse ocean environments. In response to the strong need to support marine activities and security, we present the first pan-Canadian operational regional ocean analysis system. A novel online tidal harmonic analysis method is introduced that uses a sliding-window approach. Innovations are compared to those from the Canadian global analysis system. Particular improvements are found near the Gulf Stream due to the higher model grid resolution.
William J. Pringle, Damrongsak Wirasaet, Keith J. Roberts, and Joannes J. Westerink
Geosci. Model Dev., 14, 1125–1145, https://doi.org/10.5194/gmd-14-1125-2021, https://doi.org/10.5194/gmd-14-1125-2021, 2021
Short summary
Short summary
We improve and test a computer model that simulates tides and storm surge over all of Earth's oceans and seas. The model varies mesh resolution (triangular element sizes) freely so that coastal areas, especially storm landfall locations, are well-described. We develop systematic tests of the resolution in order to suggest good mesh design criteria that balance computational efficiency with accuracy for both global astronomical tides and coastal storm tides under extreme weather forcing.
Bijoy Thompson, Claudio Sanchez, Boon Chong Peter Heng, Rajesh Kumar, Jianyu Liu, Xiang-Yu Huang, and Pavel Tkalich
Geosci. Model Dev., 14, 1081–1100, https://doi.org/10.5194/gmd-14-1081-2021, https://doi.org/10.5194/gmd-14-1081-2021, 2021
Short summary
Short summary
This article describes the development and ocean forecast evaluation of an atmosphere–ocean coupled prediction system for the Maritime Continent domain, which includes the eastern Indian and western Pacific oceans. The coupled system comprises regional configurations of the atmospheric model MetUM and ocean model NEMO, coupled using the OASIS3-MCT libraries. The model forecast deviation of selected fields relative to observations is within acceptable error limits of operational forecast models.
Pavel Perezhogin, Ilya Chernov, and Nikolay Iakovlev
Geosci. Model Dev., 14, 843–857, https://doi.org/10.5194/gmd-14-843-2021, https://doi.org/10.5194/gmd-14-843-2021, 2021
Short summary
Short summary
We describe the parallel implementation of the FEMAO model for an ice-covered sea with 2D Hilbert-curve domain decomposition. Load balancing is crucial because performance depends on the local depth. We propose, compare, and discuss four approaches to load balancing. The parallel library allowed us to modify the original sequential algorithm as little as possible. The performance increases almost linearly (tested with up to 996 CPU cores) for the model of the shallow White Sea.
Jordyn E. Moscoso, Andrew L. Stewart, Daniele Bianchi, and James C. McWilliams
Geosci. Model Dev., 14, 763–794, https://doi.org/10.5194/gmd-14-763-2021, https://doi.org/10.5194/gmd-14-763-2021, 2021
Short summary
Short summary
This project was created to understand the across-shore distribution of plankton in the California Current System. To complete this study, we used a quasi-2-D dynamical model coupled to an ecosystem model. This paper is a preliminary study to test and validate the model against data collected by the California Cooperative Oceanic Fisheries Investigations (CalCOFI). We show the solution of our model solution compares well to the data and discuss our model as a tool for further model development.
Christian Ferrarin, Marco Bajo, and Georg Umgiesser
Geosci. Model Dev., 14, 645–659, https://doi.org/10.5194/gmd-14-645-2021, https://doi.org/10.5194/gmd-14-645-2021, 2021
Short summary
Short summary
The problem of the optimization of ocean monitoring networks is tackled through the implementation of data assimilation techniques in a numerical model. The methodology has been applied to the tide gauge network in the Lagoon of Venice (Italy). The data assimilation methods allow identifying the minimum number of stations and their distribution that correctly represent the lagoon's dynamics. The methodology is easily exportable to other environments and can be extended to other variables.
Florian Lemarié, Guillaume Samson, Jean-Luc Redelsperger, Hervé Giordani, Théo Brivoal, and Gurvan Madec
Geosci. Model Dev., 14, 543–572, https://doi.org/10.5194/gmd-14-543-2021, https://doi.org/10.5194/gmd-14-543-2021, 2021
Short summary
Short summary
A simplified model of the atmospheric boundary layer (ABL) of intermediate complexity between a bulk parameterization and a full three-dimensional atmospheric model has been developed and integrated to the NEMO ocean model.
An objective in the derivation of such a simplified model is to reach an apt representation of ocean-only numerical simulations of some of the key processes associated with air–sea interactions at the characteristic scales of the oceanic mesoscale.
Kristen M. Thyng, Daijiro Kobashi, Veronica Ruiz-Xomchuk, Lixin Qu, Xu Chen, and Robert D. Hetland
Geosci. Model Dev., 14, 391–407, https://doi.org/10.5194/gmd-14-391-2021, https://doi.org/10.5194/gmd-14-391-2021, 2021
Short summary
Short summary
We modified the ROMS model to run in offline mode so that previously run fields of sea surface height and velocity fields are input to calculate tracer advection without running the full model with a larger time step; thus, it is faster. The code was tested with two advection schemes, and both are robust with over 99 % accuracy of the offline to online run after 14 simulation days. This allows for ROMS users to maximize use of new or existing output to quickly run additional tracer simulations.
Katixa Lajaunie-Salla, Frédéric Diaz, Cathy Wimart-Rousseau, Thibaut Wagener, Dominique Lefèvre, Christophe Yohia, Irène Xueref-Remy, Brian Nathan, Alexandre Armengaud, and Christel Pinazo
Geosci. Model Dev., 14, 295–321, https://doi.org/10.5194/gmd-14-295-2021, https://doi.org/10.5194/gmd-14-295-2021, 2021
Short summary
Short summary
A biogeochemical model of planktonic food webs including a carbonate balance module is applied in the Bay of Marseille (France) to represent the carbon marine cycle expected to change in the future owing to significant increases in anthropogenic emissions of CO2. The model correctly simulates the ranges and seasonal dynamics of most variables of the carbonate system (pH). This study shows that external physical forcings have an important impact on the carbonate equilibrium in this coastal area.
Qing Li and Luke Van Roekel
Geosci. Model Dev. Discuss., https://doi.org/10.5194/gmd-2020-262, https://doi.org/10.5194/gmd-2020-262, 2020
Revised manuscript accepted for GMD
Short summary
Short summary
Physical processes in the ocean span multiple spatial and temporal scales. Simultaneously resolving all these in a simulation is computationally challenging. Here we develop a more efficient technique to better study the interactions across scales, particularly focusing on the ocean surface turbulent mixing, by coupling a global ocean circulation model MPAS-Ocean and a large eddy simulation model PALM. The latter is customized and ported on GPU to further accelerate the computation.
Clément Bricaud, Julien Le Sommer, Gurvan Madec, Christophe Calone, Julie Deshayes, Christian Ethe, Jérôme Chanut, and Marina Levy
Geosci. Model Dev., 13, 5465–5483, https://doi.org/10.5194/gmd-13-5465-2020, https://doi.org/10.5194/gmd-13-5465-2020, 2020
Short summary
Short summary
In order to reduce the cost of ocean biogeochemical models, a multi-grid approach where ocean dynamics and tracer transport are computed with different spatial resolution has been developed in the NEMO v3.6 OGCM. Different experiments confirm that the spatial resolution of hydrodynamical fields can be coarsened without significantly affecting the resolved passive tracer fields. This approach leads to a factor of 7 reduction of the overhead associated with running a full biogeochemical model.
Lojze Žust, Anja Fettich, Matej Kristan, and Matjaž Ličer
Geosci. Model Dev. Discuss., https://doi.org/10.5194/gmd-2020-355, https://doi.org/10.5194/gmd-2020-355, 2020
Revised manuscript accepted for GMD
Short summary
Short summary
Adriatic basin sea level modeling is a challenging problem due to the interplay between terrain, weather, tides and seiches. Current state-of-the-art numerical models (e.g. NEMO) require large computational resources to produce reliable forecasts. In this study we propose HIDRA, a novel deep learning approach for sea level modeling, which drastically reduces the numerical cost while demonstrating predictive capabilities comparable to that of the NEMO model, outperforming it in many instances.
Yang Feng, Dimitris Menemenlis, Huijie Xue, Hong Zhang, Dustin Carroll, Yan Du, and Hui Wu
Geosci. Model Dev. Discuss., https://doi.org/10.5194/gmd-2020-321, https://doi.org/10.5194/gmd-2020-321, 2020
Revised manuscript accepted for GMD
Short summary
Short summary
Simulating coastal plume regions has been improved in global ECCOv4 with a series sensitivity tests. We find modeled SSS is closer to SMAP when using daily point-source runoff as well as increase the resolution from coarse to intermediate. Besides, the plume characteristics, freshwater transport, and critical water properties are modified greatly. But these may not happen when further increase to high resolution. The study will advance modelling land-ocean-atmosphere feedback seamlessly in ESMs.
Clark Pennelly and Paul G. Myers
Geosci. Model Dev., 13, 4959–4975, https://doi.org/10.5194/gmd-13-4959-2020, https://doi.org/10.5194/gmd-13-4959-2020, 2020
Short summary
Short summary
A high-resolution ocean simulation was carried out within the Labrador Sea, a region that low-resolution climate simulations may misrepresent. We show that small-scale eddies and their associated transport are better resolved at higher resolution than at lower resolution. These eddies transport important properties to the interior of the Labrador Sea, impacting the stratification and reducing the convection extent so that it is far more accurate when compared to what observations suggest.
Eric P. Chassignet, Stephen G. Yeager, Baylor Fox-Kemper, Alexandra Bozec, Frederic Castruccio, Gokhan Danabasoglu, Christopher Horvat, Who M. Kim, Nikolay Koldunov, Yiwen Li, Pengfei Lin, Hailong Liu, Dmitry V. Sein, Dmitry Sidorenko, Qiang Wang, and Xiaobiao Xu
Geosci. Model Dev., 13, 4595–4637, https://doi.org/10.5194/gmd-13-4595-2020, https://doi.org/10.5194/gmd-13-4595-2020, 2020
Short summary
Short summary
This paper presents global comparisons of fundamental global climate variables from a suite of four pairs of matched low- and high-resolution ocean and sea ice simulations to assess the robustness of climate-relevant improvements in ocean simulations associated with moving from coarse (∼1°) to eddy-resolving (∼0.1°) horizontal resolutions. Despite significant improvements, greatly enhanced horizontal resolution does not deliver unambiguous bias reduction in all regions for all models.
Mark E. Baird, Karen A. Wild-Allen, John Parslow, Mathieu Mongin, Barbara Robson, Jennifer Skerratt, Farhan Rizwi, Monika Soja-Woźniak, Emlyn Jones, Mike Herzfeld, Nugzar Margvelashvili, John Andrewartha, Clothilde Langlais, Matthew P. Adams, Nagur Cherukuru, Malin Gustafsson, Scott Hadley, Peter J. Ralph, Uwe Rosebrock, Thomas Schroeder, Leonardo Laiolo, Daniel Harrison, and Andrew D. L. Steven
Geosci. Model Dev., 13, 4503–4553, https://doi.org/10.5194/gmd-13-4503-2020, https://doi.org/10.5194/gmd-13-4503-2020, 2020
Short summary
Short summary
For 20+ years, the Commonwealth Science Industry and Research Organisation (CSIRO) has been developing a biogeochemical (BGC) model for coupling with a hydrodynamic and sediment model for application in estuaries, coastal waters and shelf seas. This paper provides a full mathematical description (equations, parameters), model evaluation and access to the numerical code. The model is particularly suited to applications in shallow waters where benthic processes are critical to ecosystem function.
Angela A. Bahamondes Dominguez, Anna E. Hickman, Robert Marsh, and C. Mark Moore
Geosci. Model Dev., 13, 4019–4040, https://doi.org/10.5194/gmd-13-4019-2020, https://doi.org/10.5194/gmd-13-4019-2020, 2020
Short summary
Short summary
The central Celtic Sea has previously been studied with a 1-D model called S2P3, showing discrepancies between observations and the model results due to poor representation of some processes. Therefore, the S2P3 model was developed to include zooplankton and phytoplankton cells' adaptation to changes in irradiance. Results demonstrate that better agreement with biological observations can be achieved when the model includes these processes and is adequately constrained.
Oliver Gutjahr, Nils Brüggemann, Helmuth Haak, Johann H. Jungclaus, Dian A. Putrasahan, Katja Lohmann, and Jin-Song von Storch
Geosci. Model Dev. Discuss., https://doi.org/10.5194/gmd-2020-202, https://doi.org/10.5194/gmd-2020-202, 2020
Revised manuscript accepted for GMD
Short summary
Short summary
We compare four ocean vertical mixing schemes in 100-year long coupled simulations with the Max Planck Institute Earth System Model (MPI-ESM1.2) and analyse their model biases. Overall, the mixing schemes modify biases in the ocean interior that vary with region and variable, but produce a similar global bias pattern. We therefore cannot classify any scheme as superior, but conclude that the chosen mixing scheme may be important for regional biases.
Alexander G. López, John L. Wilkin, and Julia C. Levin
Geosci. Model Dev., 13, 3709–3729, https://doi.org/10.5194/gmd-13-3709-2020, https://doi.org/10.5194/gmd-13-3709-2020, 2020
Short summary
Short summary
This article describes a regional circulation model, Doppio, for the Mid-Atlantic Bight and the Gulf of Maine. The model demonstrates useful skill in comparison to a comprehensive suite of observations. Development focused on achieving a model configuration that allows decadal-scale simulations of physical ocean circulation that can underpin studies of ecosystems and biogeochemistry. Doppio captures the temperature and salinity stratification well, along with the large-scale mean circulation.
Hiroyuki Tsujino, L. Shogo Urakawa, Stephen M. Griffies, Gokhan Danabasoglu, Alistair J. Adcroft, Arthur E. Amaral, Thomas Arsouze, Mats Bentsen, Raffaele Bernardello, Claus W. Böning, Alexandra Bozec, Eric P. Chassignet, Sergey Danilov, Raphael Dussin, Eleftheria Exarchou, Pier Giuseppe Fogli, Baylor Fox-Kemper, Chuncheng Guo, Mehmet Ilicak, Doroteaciro Iovino, Who M. Kim, Nikolay Koldunov, Vladimir Lapin, Yiwen Li, Pengfei Lin, Keith Lindsay, Hailong Liu, Matthew C. Long, Yoshiki Komuro, Simon J. Marsland, Simona Masina, Aleksi Nummelin, Jan Klaus Rieck, Yohan Ruprich-Robert, Markus Scheinert, Valentina Sicardi, Dmitry Sidorenko, Tatsuo Suzuki, Hiroaki Tatebe, Qiang Wang, Stephen G. Yeager, and Zipeng Yu
Geosci. Model Dev., 13, 3643–3708, https://doi.org/10.5194/gmd-13-3643-2020, https://doi.org/10.5194/gmd-13-3643-2020, 2020
Short summary
Short summary
The OMIP-2 framework for global ocean–sea-ice model simulations is assessed by comparing multi-model means from 11 CMIP6-class global ocean–sea-ice models calculated separately for the OMIP-1 and OMIP-2 simulations. Many features are very similar between OMIP-1 and OMIP-2 simulations, and yet key improvements in transitioning from OMIP-1 to OMIP-2 are also identified. Thus, the present assessment justifies that future ocean–sea-ice model development and analysis studies use the OMIP-2 framework.
Pedro Colombo, Bernard Barnier, Thierry Penduff, Jérôme Chanut, Julie Deshayes, Jean-Marc Molines, Julien Le Sommer, Polina Verezemskaya, Sergey Gulev, and Anne-Marie Treguier
Geosci. Model Dev., 13, 3347–3371, https://doi.org/10.5194/gmd-13-3347-2020, https://doi.org/10.5194/gmd-13-3347-2020, 2020
Short summary
Short summary
In the ocean circulation model NEMO, the representation of the overflow of dense Arctic waters through the Denmark Strait is investigated. In this
z-coordinate context, sensitivity tests show that the mixing parameterizations preferably act along the model grid slope. Thus, the representation of the overflow is more sensitive to resolution than to parameterization and is best when the numerical grid matches the local topographic slope.
Hideharu Sasaki, Shinichiro Kida, Ryo Furue, Hidenori Aiki, Nobumasa Komori, Yukio Masumoto, Toru Miyama, Masami Nonaka, Yoshikazu Sasai, and Bunmei Taguchi
Geosci. Model Dev., 13, 3319–3336, https://doi.org/10.5194/gmd-13-3319-2020, https://doi.org/10.5194/gmd-13-3319-2020, 2020
Short summary
Short summary
A quasi-global eddying ocean hindcast simulation using a new version of our model, called OFES2, was conducted to overcome several issues in its previous version. OFES2 simulated oceanic fields from 1958 to 2016 with improved global sea surface temperature and salinity, water mass properties in the Indonesian and Arabian seas, and Niño3.4 and Indian Ocean Dipole indexes. The output from OFES2 will be useful in studying various oceanic phenomena with broad spatiotemporal scales.
Katherine M. Smith, Skyler Kern, Peter E. Hamlington, Marco Zavatarelli, Nadia Pinardi, Emily F. Klee, and Kyle E. Niemeyer
Geosci. Model Dev. Discuss., https://doi.org/10.5194/gmd-2020-134, https://doi.org/10.5194/gmd-2020-134, 2020
Revised manuscript accepted for GMD
Short summary
Short summary
We present a newly developed reduced-order biogeochemical flux model that is complex and flexible enough to capture open-ocean ecosystem dynamics, but reduced enough to incorporate into highly resolved numerical simulations with limited additional computational cost. The model provides improved correlations between model output and field data, indicating that significant improvements in the reproduction of real-world data can be achieved with a small number of variables.
Dafydd Stephenson, Simon A. Müller, and Florian Sévellec
Geosci. Model Dev., 13, 2031–2050, https://doi.org/10.5194/gmd-13-2031-2020, https://doi.org/10.5194/gmd-13-2031-2020, 2020
Short summary
Short summary
Different water types are created at the sea surface with a signature based on the local conditions of the atmosphere. They then take these conditions with them into the deeper ocean, and so their journey is an important climate process to understand. In this study, we modify and repurpose a specialised model which simulates the ocean forward and backward in time to determine where new ocean water goes, where at the surface existing water comes from, and how old it is, by tracking it as a dye.
Violaine Piton, Marine Herrmann, Florent Lyard, Patrick Marsaleix, Thomas Duhaut, Damien Allain, and Sylvain Ouillon
Geosci. Model Dev., 13, 1583–1607, https://doi.org/10.5194/gmd-13-1583-2020, https://doi.org/10.5194/gmd-13-1583-2020, 2020
Short summary
Short summary
Consequences of tidal dynamics on hydro-sedimentary processes are a recurrent issue in estuarine and coastal processes studies, and accurate tidal solutions are a prerequisite for modeling sediment transport. This study presents the implementation and optimization of a model configuration in terms of bathymetry and bottom friction and assess the influence of these parameters on tidal solutions, in a macro-tidal environment: the Gulf of Tonkin (Vietnam).
Alexander Barth, Aida Alvera-Azcárate, Matjaz Licer, and Jean-Marie Beckers
Geosci. Model Dev., 13, 1609–1622, https://doi.org/10.5194/gmd-13-1609-2020, https://doi.org/10.5194/gmd-13-1609-2020, 2020
Short summary
Short summary
DINCAE is a method for reconstructing missing data in satellite datasets using a neural network. Satellite observations working in the optical and infrared bands are affected by clouds, which obscure part of the ocean underneath. In this paper, a neural network with the structure of a convolutional auto-encoder is developed to reconstruct the missing data based on the available cloud-free pixels in satellite images.
Jiangyu Li and Shaoqing Zhang
Geosci. Model Dev., 13, 1035–1054, https://doi.org/10.5194/gmd-13-1035-2020, https://doi.org/10.5194/gmd-13-1035-2020, 2020
Short summary
Short summary
Two assimilation systems developed using two nearly independent wave models are used to study the influences of various error sources including mode bias on wave data assimilation; a statistical method is explored to make full use of the merits of individual assimilation systems for bias correction, thus improving wave analysis greatly. This study opens a door to further our understanding of physical processes in waves and associated air–sea interactions for improving wave modeling.
Heiner Dietze, Ulrike Löptien, and Julia Getzlaff
Geosci. Model Dev., 13, 71–97, https://doi.org/10.5194/gmd-13-71-2020, https://doi.org/10.5194/gmd-13-71-2020, 2020
Short summary
Short summary
We present a new near-global coupled biogeochemical ocean-circulation model configuration of the Southern Ocean. The configuration features both a relatively equilibrated oceanic carbon inventory and an explicit representation of mesoscale eddies. In this paper, we document the model configuration and showcase its potential to tackle research questions such as the Southern Ocean carbon uptake dynamics on decadal timescales.
Jaap H. Nienhuis and Jorge Lorenzo-Trueba
Geosci. Model Dev., 12, 4013–4030, https://doi.org/10.5194/gmd-12-4013-2019, https://doi.org/10.5194/gmd-12-4013-2019, 2019
Short summary
Short summary
The response of barrier islands to sea level rise depends on their ability to move landward through the transport of sediment from the beach to the back barrier. The BRIE model simulates these processes and the resulting landward movement of barrier islands. The novelty of the BRIE model is the incorporation of tidal inlets (gaps between barrier islands) that can transport sediment landward and therefore help keep barrier islands above sea level.
Bjarne Büchmann
Geosci. Model Dev., 12, 3915–3922, https://doi.org/10.5194/gmd-12-3915-2019, https://doi.org/10.5194/gmd-12-3915-2019, 2019
Short summary
Short summary
Operational forecasting of the ocean state – e.g. used for ship route planning, sea rescue, and oil spill drift models – relies on data (forcing) obtained from weather forecasting. Unfortunately, the so-called meteorological analysis step introduces a discontinuity, which affects the ocean models adversely. In the present paper, a straightforward method to deal with the issue is introduced. Practical examples are given to illuminate the scale of the problem.
Philippe Delandmeter and Erik van Sebille
Geosci. Model Dev., 12, 3571–3584, https://doi.org/10.5194/gmd-12-3571-2019, https://doi.org/10.5194/gmd-12-3571-2019, 2019
Short summary
Short summary
Parcels is a framework to compute how ocean currents transport
stuffsuch as plankton and plastic around. In the latest version 2.0 of Parcels, we focus on more accurate interpolation schemes and implement methods to seamlessly combine data from different sources (such as winds and currents, possibly in different regions). We show that this framework is very efficient for tracking how microplastic is transported through the North Sea into the Arctic.
Franziska U. Schwarzkopf, Arne Biastoch, Claus W. Böning, Jérôme Chanut, Jonathan V. Durgadoo, Klaus Getzlaff, Jan Harlaß, Jan K. Rieck, Christina Roth, Markus M. Scheinert, and René Schubert
Geosci. Model Dev., 12, 3329–3355, https://doi.org/10.5194/gmd-12-3329-2019, https://doi.org/10.5194/gmd-12-3329-2019, 2019
Short summary
Short summary
A family of nested global ocean general circulation model configurations, the INALT family, has been established with resolutions of 1/10°, 1/20° and 1/60° in the South Atlantic and western Indian oceans, covering the greater Agulhas Current (AC) system. The INALT family provides a consistent set of configurations that allows to address eddy dynamics in the AC system and their impact on the large-scale ocean circulation.
Thomas Rackow, Dmitry V. Sein, Tido Semmler, Sergey Danilov, Nikolay V. Koldunov, Dmitry Sidorenko, Qiang Wang, and Thomas Jung
Geosci. Model Dev., 12, 2635–2656, https://doi.org/10.5194/gmd-12-2635-2019, https://doi.org/10.5194/gmd-12-2635-2019, 2019
Short summary
Short summary
Current climate models show errors in the deep ocean that are larger than the level of natural variability and the response to enhanced greenhouse gas concentrations. These errors are larger than the signals we aim to predict. With the AWI Climate Model, we show that increasing resolution to resolve eddies can lead to major reductions in deep ocean errors. AWI's next-generation (CMIP6) model configuration will thus use locally eddy-resolving computational grids for projecting climate change.
Keith J. Roberts, William J. Pringle, and Joannes J. Westerink
Geosci. Model Dev., 12, 1847–1868, https://doi.org/10.5194/gmd-12-1847-2019, https://doi.org/10.5194/gmd-12-1847-2019, 2019
Short summary
Short summary
Computer simulations can be used to reproduce the dynamics of the ocean near the coast. These simulations often use a mesh of triangles to represent the domain since they can be orientated and disparately sized in such a way to accurately fit the coastline shape. This paper describes a software package (OceanMesh2D v1.0) that has been developed in order to automatically and objectively design triangular meshes based on geospatial data inputs that represent the coastline and ocean depths.
Pearse J. Buchanan, Richard J. Matear, Zanna Chase, Steven J. Phipps, and Nathan L. Bindoff
Geosci. Model Dev., 12, 1491–1523, https://doi.org/10.5194/gmd-12-1491-2019, https://doi.org/10.5194/gmd-12-1491-2019, 2019
Short summary
Short summary
Oceanic sediment cores are commonly used to understand past climates. The composition of the sediments changes with the ocean above it. An understanding of oceanographic conditions that existed many thousands of years ago, in some cases many millions of years ago, can therefore be extracted from sediment cores. We simulate two chemical signatures (13C and 15N) of sediment cores in a model. This study assesses the model before it is applied to reinterpret the sedimentary record.
Carlos Rocha, Christopher A. Edwards, Moninya Roughan, Paulina Cetina-Heredia, and Colette Kerry
Geosci. Model Dev., 12, 441–456, https://doi.org/10.5194/gmd-12-441-2019, https://doi.org/10.5194/gmd-12-441-2019, 2019
Short summary
Short summary
Off southeast Australia, the East Australian Current (EAC) moves warm nutrient-poor waters towards the pole. In this region, the EAC and a large number of vortices pinching off it strongly affect phytoplankton’s access to nutrients and light. To study these dynamics, we created a numerical model that is able to solve the ocean conditions and how they modulate the foundation of the region’s ecosystem. We validated model results against available data and this showed that the model performs well.
Robinson Hordoir, Lars Axell, Anders Höglund, Christian Dieterich, Filippa Fransner, Matthias Gröger, Ye Liu, Per Pemberton, Semjon Schimanke, Helen Andersson, Patrik Ljungemyr, Petter Nygren, Saeed Falahat, Adam Nord, Anette Jönsson, Iréne Lake, Kristofer Döös, Magnus Hieronymus, Heiner Dietze, Ulrike Löptien, Ivan Kuznetsov, Antti Westerlund, Laura Tuomi, and Jari Haapala
Geosci. Model Dev., 12, 363–386, https://doi.org/10.5194/gmd-12-363-2019, https://doi.org/10.5194/gmd-12-363-2019, 2019
Short summary
Short summary
Nemo-Nordic is a regional ocean model based on a community code (NEMO). It covers the Baltic and the North Sea area and is used as a forecast model by the Swedish Meteorological and Hydrological Institute. It is also used as a research tool by scientists of several countries to study, for example, the effects of climate change on the Baltic and North seas. Using such a model permits us to understand key processes in this coastal ecosystem and how such processes will change in a future climate.
Hagen Radtke, Marko Lipka, Dennis Bunke, Claudia Morys, Jana Woelfel, Bronwyn Cahill, Michael E. Böttcher, Stefan Forster, Thomas Leipe, Gregor Rehder, and Thomas Neumann
Geosci. Model Dev., 12, 275–320, https://doi.org/10.5194/gmd-12-275-2019, https://doi.org/10.5194/gmd-12-275-2019, 2019
Short summary
Short summary
This paper describes a coupled benthic–pelagic biogeochemical model, ERGOM-SED. We demonstrate its use in a one-dimensional physical model, which is horizontally integrated and vertically resolved. We describe the application of the model to seven stations in the south-western Baltic Sea. The model was calibrated using pore water profiles from these stations. We compare the model results to these and to measured sediment compositions, benthopelagic fluxes and bioturbation intensities.
Dale Partridge, Tobias Friedrich, and Brian S. Powell
Geosci. Model Dev., 12, 195–213, https://doi.org/10.5194/gmd-12-195-2019, https://doi.org/10.5194/gmd-12-195-2019, 2019
Short summary
Short summary
This paper demonstrates the improvements made to an operational ocean forecast model around the Hawaiian Islands by performing a reanalysis of the model over a 10-year period. Using a number of different measurements we show the role a variety of observations play in producing the forecast, in particular the contribution of high-frequency radar.
Tuomas Kärnä, Stephan C. Kramer, Lawrence Mitchell, David A. Ham, Matthew D. Piggott, and António M. Baptista
Geosci. Model Dev., 11, 4359–4382, https://doi.org/10.5194/gmd-11-4359-2018, https://doi.org/10.5194/gmd-11-4359-2018, 2018
Short summary
Short summary
Unstructured meshes are attractive for coastal ocean modeling, as they allow more accurate representation of complex coastal topography. Unstructured mesh models are, however, often perceived as slow and inaccurate. We present a novel discontinuous Galerkin ocean model: Thetis. We demonstrate that the model is able to simulate baroclinic ocean flows with high accuracy on a triangular prismatic mesh. This work paves the way for highly accurate and efficient three-dimensional coastal ocean models.
Paul Sandery
Geosci. Model Dev., 11, 4011–4019, https://doi.org/10.5194/gmd-11-4011-2018, https://doi.org/10.5194/gmd-11-4011-2018, 2018
Short summary
Short summary
This article compares global mesoscale ocean forecasts with different data assimilation cycle lengths. Mean absolute increment is used to quantify differences in the overall impact of observations. Greater observation impact does not necessarily improve a forecast system. Experiments show a 1-day cycle generates improved 7-day forecasts when compared to a 3-day cycle. Cycle length is an important choice that influences system bias and predictability.
Daniel Boettger, Robin Robertson, and Gary B. Brassington
Geosci. Model Dev., 11, 3795–3805, https://doi.org/10.5194/gmd-11-3795-2018, https://doi.org/10.5194/gmd-11-3795-2018, 2018
Short summary
Short summary
This study focuses on the impact of the model vertical mixing parameterisation on the representation of the mixed layer depth (MLD) in ocean forecast models. We compare data from two recent versions of the OceanMAPS forecast system, and find that while there were large improvements in the later version of the model, the skill of each parameterisation varies with spatial location.
Marco van Hulten, Jean-Claude Dutay, and Matthieu Roy-Barman
Geosci. Model Dev., 11, 3537–3556, https://doi.org/10.5194/gmd-11-3537-2018, https://doi.org/10.5194/gmd-11-3537-2018, 2018
Short summary
Short summary
We present an ocean model of the natural radioactive isotopes thorium-230 and protactinium-231. These isotopes are often used to investigate past ocean circulation and particle transport. They are removed by particles produced by plankton and from uplifted desert dust that is deposited into the ocean. We approach observed dissolved and adsorbed Th-230 and Pa-231 activities. The Pa-231 / Th-230 sedimentation ratio is reproduced as well; this quantity can be used as a proxy for ocean circulation.
Dion Häfner, René Løwe Jacobsen, Carsten Eden, Mads R. B. Kristensen, Markus Jochum, Roman Nuterman, and Brian Vinter
Geosci. Model Dev., 11, 3299–3312, https://doi.org/10.5194/gmd-11-3299-2018, https://doi.org/10.5194/gmd-11-3299-2018, 2018
Short summary
Short summary
Well-performing, easy-to-use ocean models are a central ingredient to further the understanding of our Earth and climate. Veros, the versatile ocean simulator, is the first full-blown ocean model entirely written in the high-level programming language Python. It is considerably more approachable than traditional Fortran models and leverages modern best practices; at the same time, thanks to the Bohrium framework, Veros is about half as fast as a reference implementation in Fortran 90.
Christopher R. Sherwood, Alfredo L. Aretxabaleta, Courtney K. Harris, J. Paul Rinehimer, Romaric Verney, and Bénédicte Ferré
Geosci. Model Dev., 11, 1849–1871, https://doi.org/10.5194/gmd-11-1849-2018, https://doi.org/10.5194/gmd-11-1849-2018, 2018
Short summary
Short summary
Cohesive sediment (mud) is ubiquitous in the world's coastal regions, but its behavior is complicated and often oversimplified by computer models. This paper describes extensions to a widely used open-source coastal ocean model that allow users to simulate important components of cohesive sediment transport.
Knut-Frode Dagestad, Johannes Röhrs, Øyvind Breivik, and Bjørn Ådlandsvik
Geosci. Model Dev., 11, 1405–1420, https://doi.org/10.5194/gmd-11-1405-2018, https://doi.org/10.5194/gmd-11-1405-2018, 2018
Short summary
Short summary
We have developed a computer code with ability to predict how various substances and objects drift in the ocean. This may be used to, e.g. predict the drift of oil to aid cleanup operations, the drift of man-over-board or lifeboats to aid search and rescue operations, or the drift of fish eggs and larvae to understand and manage fish stocks. This new code merges all such applications into one software tool, allowing to optimise and channel any available resources and developments.
Qiang Wang, Claudia Wekerle, Sergey Danilov, Xuezhu Wang, and Thomas Jung
Geosci. Model Dev., 11, 1229–1255, https://doi.org/10.5194/gmd-11-1229-2018, https://doi.org/10.5194/gmd-11-1229-2018, 2018
Short summary
Short summary
For developing a system for Arctic research, we evaluate the Arctic Ocean simulated by FESOM. We use two global meshes differing in the horizontal resolution only in the Arctic Ocean (24 vs. 4.5 km). The high resolution significantly improves the model's representation of the Arctic Ocean. The most pronounced improvement is in the Arctic intermediate layer. The high resolution also improves the ocean surface circulation, mainly through a better representation of the Canadian Arctic Archipelago.
Jennifer A. Graham, Enda O'Dea, Jason Holt, Jeff Polton, Helene T. Hewitt, Rachel Furner, Karen Guihou, Ashley Brereton, Alex Arnold, Sarah Wakelin, Juan Manuel Castillo Sanchez, and C. Gabriela Mayorga Adame
Geosci. Model Dev., 11, 681–696, https://doi.org/10.5194/gmd-11-681-2018, https://doi.org/10.5194/gmd-11-681-2018, 2018
Short summary
Short summary
This paper describes the next-generation ocean forecast model for the European NW shelf, AMM15 (Atlantic Margin Model, 1.5 km resolution). The current forecast system has a resolution of 7 km. While this is sufficient to represent large-scale circulation, many dynamical features (such as eddies, frontal jets, and internal tides) can only begin to be resolved at 0–1 km resolution. Here we introduce AMM15 and demonstrate its ability to represent the mean state and variability of the region.
Julien Chauchat, Zhen Cheng, Tim Nagel, Cyrille Bonamy, and Tian-Jian Hsu
Geosci. Model Dev., 10, 4367–4392, https://doi.org/10.5194/gmd-10-4367-2017, https://doi.org/10.5194/gmd-10-4367-2017, 2017
Short summary
Short summary
This manuscript presents the development and validation of a two-phase flow Eulerian-Eulerian model based on OpenFOAM for sediment transport applications. The mathematical and numerical models are described in detail. The numerical implementation is demonstrated on four test cases: sedimentation of suspended particles, laminar bed load, sheet flow, and scour at an apron. These test cases illustrate the capabilities of SedFoam to deal with complex turbulent sediment transport problems.
Michael Lange and Erik van Sebille
Geosci. Model Dev., 10, 4175–4186, https://doi.org/10.5194/gmd-10-4175-2017, https://doi.org/10.5194/gmd-10-4175-2017, 2017
Short summary
Short summary
Here, we present version 0.9 of Parcels (Probably A Really Computationally Efficient Lagrangian Simulator). Parcels is an experimental prototype code aimed at exploring novel approaches for Lagrangian tracking of virtual ocean particles in the petascale age. The modularity, flexibility and scalability will allow the code to be used to track water, nutrients, microbes, plankton, plastic and even fish.
Cited articles
Ali, S. and Shah, M.: A Lagrangian Particle Dynamics Approach for Crowd Flow Segmentation and Stability Analysis, in: 2007 IEEE Conference on Computer Vision and Pattern Recognition, Minneapolis, MN, USA, 17–22 June 2007, IEEE, 1–6, https://doi.org/10.1109/CVPR.2007.382977, 2007. a, b
Allshouse, M. R., Ivey, G. N., Lowe, R. J., Jones, N. L., Beegle-Krause, C., Xu, J., and Peacock, T.: Impact of windage on ocean surface Lagrangian coherent structures, Environ. Fluid Mech., 17, 473–483, 2017. a
Barkan, R., McWilliams, J. C., Molemaker, M. J., Choi, J., Srinivasan, K., Shchepetkin, A. F., and Bracco, A.: Submesoscale dynamics in the northern Gulf of Mexico. Part II: Temperature–salinity relations and cross-shelf transport processes, J. Phys. Oceanogr., 47, 2347–2360, 2017. a
Beron-Vera, F. J., Olascoaga, M. J., and Goni, G. J.: Oceanic mesoscale eddies as revealed by Lagrangian coherent structures, Geophys. Res. Lett., 35, L12603, https://doi.org/10.1029/2008GL033957, 2008. a
Beron-Vera, F. J., Olascoaga, M. J., Brown, M. G., Koçak, H., and Rypina, I. I.: Invariant-tori-like Lagrangian coherent structures in geophysical flows, Chaos, 20, 017514, https://doi.org/10.1063/1.3271342, 2010. a
Bogacki, P. and Shampine, L. F.: A 3(2) pair of Runge-Kutta formulas, Appl. Math. Lett., 2, 321–325, 1989. a
Breivik, Ø. and Allen, A. A.: An operational search and rescue model for the Norwegian Sea and the North Sea, J. Marine Syst., 69, 99–113, 2008. a
de Boor, C.: A practical guide to splines, Springer-Verlag, New York Berlin Heidelberg, 2001. a
De Dominicis, M., Pinardi, N., Zodiatis, G., and Lardner, R.: MEDSLIK-II, a Lagrangian marine surface oil spill model for short-term forecasting – Part 1: Theory, Geosci. Model Dev., 6, 1851–1869, https://doi.org/10.5194/gmd-6-1851-2013, 2013. a
Dieci, L. and Lopez, L.: A survey of numerical methods for IVPs of ODEs with discontinuous right-hand side, J. Comput. Appl. Math., 236, 3967–3991, https://doi.org/10.1016/j.cam.2012.02.011, 2012. a
Dormand, J. and Prince, P.: A family of embedded Runge-Kutta formulae, J. Comput. Appl. Math., 6, 19–26, 1980. a
Dormand, J. and Prince, P.: A reconsideration of some embedded Runge-Kutta formulae, J. Comput. Appl. Math., 15, 203–211, https://doi.org/10.1016/0377-0427(86)90027-0, 1986. a
Dugstad, J., Fer, I., LaCasce, J., Sanchez de La Lama, M., and Trodahl, M.: Lateral Heat Transport in the Lofoten Basin: Near-Surface Pathways and Subsurface Exchange, J. Geophys. Res.-Oceans,124, 2992–3006, https://doi.org/10.1029/2018JC014774, 2019. a
Duran, R., Beron-Vera, F. J., and Olascoaga, M. J.: Extracting quasi-steady Lagrangian transport patterns from the ocean circulation: An application to the Gulf of Mexico, Scientific Reports, 8, 5218, https://doi.org/10.1038/s41598-018-23121-y, 2018. a, b
Enright, W., Jackson, K., Nørsett, S., and Thomsen, P.: Effective solution of discontinuous IVPs using a Runge-Kutta formula pair with interpolants, Appl. Math. Comput., 27, 313–335, https://doi.org/10.1016/0096-3003(88)90030-6, 1988. a
Farazmand, M. and Haller, G.: Computing Lagrangian coherent structures from their variational theory, Chaos, 22, 013128, https://doi.org/10.1063/1.3690153, 2012. a
García-Martínez, R. and Flores-Tovar, H.: Computer modeling of oil spill trajectories with a high accuracy method, Spill Sci. Technol. B., 5, 323–330, 1999. a
Gladwell, I., Shampine, L., and Thompson, S.: Solving ODEs with MATLAB, Cambridge University Press, New York, NY, USA, 2003. a
Gräwe, U.: Implementation of high-order particle-tracking schemes in a water column model, Ocean Model., 36, 80–89, 2011. a
Griffiths, D. F. and Higham, D. J.: Numerical methods for ordinary differential equations, Springer-Verlag, London, 2010. a
Hairer, E. and Wanner, G.: Solving Ordinary Differential Equations II: Stiff and Differential-Algebraic Problems, Springer-Verlag, Berlin, Heidelberg, 1996. a
Hairer, E., Nørsett, S. P., and Wanner, G.: Solving Ordinary Differential Equations I: Nonstiff Problems, 1st edn., Springer-Verlag Berlin Heidelberg, ISBN 978-3-662-12609-7, https://doi.org/10.1007/978-3-662-12607-3, 1987. a, b, c
Hairer, E., Wanner, G., and Lubich, C.: Geometric Numerical Integration, in:
Structure-Preserving Algorithms for Ordinary Differential Equations, Springer-Verlag, Berlin, Heidelberg, https://doi.org/10.1007/3-540-30666-8, 2006. a
Isaacson, E. and Keller, H. B.: Analysis of Numerical Methods, Dover Publications, New York, USA, 1994. a
Kloeden, P. E. and Platen, E.: Numerical Solution of Stochastic Differential
Equations, Springer-Verlag, Berlin, Heidelberg, 1992. a
Kress, R.: Numerical Analysis, in: Graduate Texts in Mathematics, Springer-Verlag, New York, https://doi.org/10.1007/978-1-4612-0599-9, 1998. a
Lekien, F. and Marsden, J.: Tricubic interpolation in three dimensions, Int. J. Numer. Meth. Eng., 63, 455–471, https://doi.org/10.1002/nme.1296, 2005. a, b
Lévy, M., Resplandy, L., Klein, P., Capet, X., Iovino, D., and Éthé, C.: Grid degradation of submesoscale resolving ocean models: Benefits for offline passive tracer transport, Ocean Model., 48, 1–9, 2012. a
Maslo, A., de Souza, J. M. A. C., Andrade-Canto, F., and Outerelo, J. R.: Connectivity of deep waters in the Gulf of Mexico, J. Marine
Syst., 203, 103267, https://doi.org/10.1016/j.jmarsys.2019.103267, 2020. a
Narváez, D. A., Klinck, J. M., Powell, E. N., Hofmann, E. E., Wilkin, J., and Haidvogel, D. B.: Modeling the dispersal of eastern oyster (Crassostrea virginica) larvae in Delaware Bay, J. Mar. Res., 70, 381–409, 2012. a
Nieto, R. and Gimeno, L.: A database of optimal integration times for Lagrangian studies of atmospheric moisture sources and sinks, Scientific Data, 6, 59, https://doi.org/10.1038/s41597-019-0068-8, 2019. a
Nordam, T.: nordam/ODE-integrators-for-Lagrangian-particles 0.9, Version 0.9, Zenodo, https://doi.org/10.5281/zenodo.4041979, 2020. a, b
Nordam, T., Brønner, U., Skancke, J., Nepstad, R., Rønningen, P., and Alver, M. O.: Numerical integration and interpolation in marine pollutant transport modelling, in: Proceedings of the 40th AMOP Technical Seminar, Calgary, AB, Canada, 3–5 October 2017, Environment and Climate Change Canada, Ottawa, 586–609, https://hdl.handle.net/11250/2652834, 2017. a
North, E. W., Adams, E. E., Schlag, Z., Sherwood, C. R., He, R., Hyun, K. H., and Socolofsky, S. A.: Simulating Oil Droplet Dispersal From the Deepwater Horizon Spill With a Lagrangian Approach, in: Monitoring and Modeling the Deepwater Horizon Oil Spill: A Record Breaking Enterprise, Wiley, 195, 217–226, 2011. a
Onink, V., Wichmann, D., Delandmeter, P., and van Sebille, E.: The Role of Ekman Currents, Geostrophy, and Stokes Drift in the Accumulation of Floating Microplastic, J. Geophys. Res.-Oceans, 124, 1474–1490, https://doi.org/10.1029/2018JC014547, 2019. a
Onu, K., Huhn, F., and Haller, G.: LCS Tool: A computational platform for Lagrangian coherent structures, J. Computat. Sci., 7, 26–36, 2015. a
Peng, J. and Dabiri, J. O.: Transport of inertial particles by Lagrangian coherent structures: application to predator–prey interaction in jellyfish feeding, J. Fluid Mech., 623, 75–84, https://doi.org/10.1017/S0022112008005089, 2009. a
Povinec, P., Gera, M., Holý, K., Hirose, K., Lujaniené, G., Nakano, M., Plastino, W., Sýkora, I., Bartok, J., and Gažák, M.: Dispersion of Fukushima radionuclides in the global atmosphere and the ocean, Appl. Radiat. Isotopes, 81, 383–392, https://doi.org/10.1016/j.apradiso.2013.03.058, 2013. a
Riuttanen, L., Hulkkonen, M., Dal Maso, M., Junninen, H., and Kulmala, M.: Trajectory analysis of atmospheric transport of fine particles, SO2, NOx and O3 to the SMEAR II station in Finland in 1996–2008, Atmos. Chem. Phys., 13, 2153–2164, https://doi.org/10.5194/acp-13-2153-2013, 2013. a
Rivas, D. and Samelson, R. M.: A Numerical Modeling Study of the Upwelling Source Waters along the Oregon Coast during 2005, J. Phys. Oceanogr., 41, 88–112, https://doi.org/10.1175/2010JPO4327.1, 2011. a, b
Rye, H., Reed, M., and Ekrol, N.: The ParTrack model for calculation of the spreading and deposition of drilling mud, chemicals and drill cuttings, Environ. Modell. Softw., 13, 431–441, 1998. a
Serra, M., Sathe, P., Rypina, I., Kirincich, A., Ross, S. D., Lermusiaux, P., Allen, A., Peacock, T., and Haller, G.: Search and rescue at sea aided by hidden flow structures, Nat. Commun., 11, 2525, https://doi.org/10.1038/s41467-020-16281-x, 2020. a
Shadden, S. C. and Taylor, C. A.: Characterization of Coherent Structures in the Cardiovascular System, Ann. Biomed. Eng., 36, 1152–1162, https://doi.org/10.1007/s10439-008-9502-3, 2008. a
Shadden, S. C., Astorino, M., and Gerbeau, J.-F.: Computational analysis of an aortic valve jet with Lagrangian coherent structures, Chaos: An Interdisciplinary J. Nonlinear Sci., 20, 017512, https://doi.org/10.1063/1.3272780, 2010. a
Siegel, D., Kinlan, B., Gaylord, B., and Gaines, S.: Lagrangian descriptions of marine larval dispersion, Mar. Ecol. Prog. Ser., 260, 83–96, 2003. a
Sirois, A. and Bottenheim, J. W.: Use of backward trajectories to interpret the 5-year record of PAN and O3 ambient air concentrations at Kejimkujik National Park, Nova Scotia, J. Geophys. Res., 100, 2867–2881, https://doi.org/10.1029/94JD02951, 1995. a
Spivakovskaya, D., Heemink, A. W., and Deleersnijder, E.: Lagrangian modelling of multi-dimensional advection-diffusion with space-varying diffusivities: theory and idealized test cases, Ocean Dynam., 57, 189–203, 2007. a
van Sebille, E., Scussolini, P., Durgadoo, J. V., Peeters, F. J., Biastoch, A., Weijer, W., Turney, C., Paris, C. B., and Zahn, R.: Ocean currents generate large footprints in marine palaeoclimate proxies, Nat. Commun., 6, 6521, https://doi.org/10.1038/ncomms7521, 2015. a
van Sebille, E., Griffies, S. M., Abernathey, R., Adams, T. P., Berloff, P., Biastoch, A., Blanke, B., Chassignet, E. P., Cheng, Y., Cotter, C. J., Deleersnijder, E., Döös, K., Drake, H. F., Drijfhout, S., Gary, S. F., Heemink, A. W., Kjellsson, J., Koszalka, I. M., Lange, M., Lique, C., MacGilchrist, G. A., Marsh, R., Mayorga Adame, C. G., McAdam, R., Nencioli, F., Paris, C. B., Piggott, M. D., Polton, J. A., Rühs, S., Shah, S. H., Thomas, M. D., Wang, J., Wolfram, P. J., Zanna, L., and Zika, J. D.: Lagrangian ocean analysis: Fundamentals and practices, Ocean Model., 121, 49–75, 2018. a
Visser, A. W.: Lagrangian modelling of plankton motion: From deceptively simple random walks to Fokker-Planck and back again, J. Marine Syst., 70, 287–299, 2008. a
Williams, J.: Bspline-Fortran: Multidimensional B-Spline Interpolation of Data on a Regular Grid (Version 5.4.0), Zenodo, https://doi.org/10.5281/zenodo.1215290, 2018. a
Woods, J.: The Lagrangian Ensemble metamodel for simulating plankton ecosystems, Prog. Oceanogr., 67, 84–159, 2005. a
Yang, Y., He, G.-W., and Wang, L.-P.: Effects of subgrid-scale modeling on Lagrangian statistics in large-eddy simulation, J. Turbul., 9, 1–24, https://doi.org/10.1080/14685240801905360, 2008. a
Zelenke, B., O'Connor, C., Barker, C., Beegle-Krause, C. J., and Eclipse, L. (Eds.): General NOAA Operational Modeling Environment (GNOME) Technical Documentation, U.S. Dept. of Commerce, NOAA Technical Memorandum NOS OR&R 40, Emergency Response Division, Seattle, WA, NOAA, 105 pp., available at: https://response.restoration.noaa.gov/sites/default/files/GNOME_Tech_Doc.pdf (last access: 17 November 2020), 2012. a
Short summary
In applied oceanography, a common task is to calculate the trajectory of objects floating at the sea surface or submerged in the water. We have investigated different numerical methods for doing such calculations and discuss the benefits and challenges of some common methods. We then propose a small change to some common methods that make them more efficient for this particular application. This will allow researchers to obtain more accurate answers with fewer computer resources.
In applied oceanography, a common task is to calculate the trajectory of objects floating at the...
