Articles | Volume 13, issue 12
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.
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
| 01 Dec 2020
Numerical integrators for Lagrangian oceanography
SINTEF Ocean, Trondheim, Norway
Department of Physics, Norwegian University of Science and Technology, Trondheim, Norway
Rodrigo Duran
National Energy Technology Laboratory, Albany, OR 97321, USA
Theiss Research, San Diego, CA 92037, USA
Related authors
Tor Nordam, Ruben Kristiansen, Raymond Nepstad, Erik van Sebille, and Andy M. Booth
Geosci. Model Dev., 16, 5339–5363, https://doi.org/10.5194/gmd-16-5339-2023, https://doi.org/10.5194/gmd-16-5339-2023, 2023
Short summary
Short summary
We describe and compare two common methods, Eulerian and Lagrangian models, used to simulate the vertical transport of material in the ocean. They both solve the same transport problems but use different approaches for representing the underlying equations on the computer. The main focus of our study is on the numerical accuracy of the two approaches. Our results should be useful for other researchers creating or using these types of transport models.
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.
Luca Kunz, Alexa Griesel, Carsten Eden, Rodrigo Duran, and Bruno Sainte-Rose
EGUsphere, https://doi.org/10.5194/egusphere-2024-1215, https://doi.org/10.5194/egusphere-2024-1215, 2024
Short summary
Short summary
TRansient Attracting Profiles (TRAPs) indicate the most attractive regions of the flow and have the potential to facilitate offshore cleanups in the Great Pacific Garbage Patch. We study the characteristics of TRAPs and the prospects for predicting debris transport from a mesoscale permitting dataset. Our findings show the relevance of TRAP lifetime estimations to an operational application and our algorithm to track TRAPs may benefit even more challenges that are related to the search at sea.
Tor Nordam, Ruben Kristiansen, Raymond Nepstad, Erik van Sebille, and Andy M. Booth
Geosci. Model Dev., 16, 5339–5363, https://doi.org/10.5194/gmd-16-5339-2023, https://doi.org/10.5194/gmd-16-5339-2023, 2023
Short summary
Short summary
We describe and compare two common methods, Eulerian and Lagrangian models, used to simulate the vertical transport of material in the ocean. They both solve the same transport problems but use different approaches for representing the underlying equations on the computer. The main focus of our study is on the numerical accuracy of the two approaches. Our results should be useful for other researchers creating or using these types of transport models.
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
DELWAVE 1.0: deep learning surrogate model of surface wave climate in the Adriatic Basin
StraitFlux – precise computations of water strait fluxes on various modeling grids
Comparison of the Coastal and Regional Ocean COmmunity model (CROCO) and NCAR-LES in non-hydrostatic simulations
Intercomparisons of Tracker v1.1 and four other ocean particle-tracking software packages in the Regional Ocean Modeling System
CAR36, a regional high-resolution ocean forecasting system for improving drift and beaching of Sargassum in the Caribbean archipelago
Implementation of additional spectral wave field exchanges in a three-dimensional wave–current coupled WAVEWATCH-III (version 6.07) and CROCO (version 1.2) configuration: assessment of their implications for macro-tidal coastal hydrodynamics
Comparison of 4-dimensional variational and ensemble optimal interpolation data assimilation systems using a Regional Ocean Modeling System (v3.4) configuration of the eddy-dominated East Australian Current system
LOCATE v1.0: numerical modelling of floating marine debris dispersion in coastal regions using Parcels v2.4.2
An optimal transformation method applied to diagnosing the ocean carbon sink
New insights into the South China Sea throughflow and water budget seasonal cycle: evaluation and analysis of a high-resolution configuration of the ocean model SYMPHONIE version 2.4
MQGeometry-1.0: a multi-layer quasi-geostrophic solver on non-rectangular geometries
Parameter estimation for ocean background vertical diffusivity coefficients in the Community Earth System Model (v1.2.1) and its impact on El Niño–Southern Oscillation forecasts
Towards a real-time modeling of global ocean waves by the fully GPU-accelerated spectral wave model WAM6-GPU
Great Lakes wave forecast system on high-resolution unstructured meshes
Development of a Novel Storm Surge Inundation Model Framework for Efficient Prediction
Skin Sea Surface Temperature schemes in coupled ocean-atmosphere modeling: the impact of chlorophyll-interactive e-folding depth
A simple approach to represent precipitation-derived freshwater fluxes into nearshore ocean models: an FVCOM4.1 case study
Impact of increased resolution on Arctic Ocean simulations in Ocean Model Intercomparison Project phase 2 (OMIP-2)
A high-resolution physical–biogeochemical model for marine resource applications in the northwest Atlantic (MOM6-COBALT-NWA12 v1.0)
A flexible z-layers approach for the accurate representation of free surface flows in a coastal ocean model (SHYFEM v. 7_5_71)
Implementation and assessment of a model including mixotrophs and the carbonate cycle (Eco3M_MIX-CarbOx v1.0) in a highly dynamic Mediterranean coastal environment (Bay of Marseille, France) – Part 1: Evolution of ecosystem composition under limited light and nutrient conditions
Ocean wave tracing v.1: a numerical solver of the wave ray equations for ocean waves on variable currents at arbitrary depths
Design and evaluation of an efficient high-precision ocean surface wave model with a multiscale grid system (MSG_Wav1.0)
Evaluation of the CMCC global eddying ocean model for the Ocean Model Intercomparison Project (OMIP2)
Barents-2.5km v2.0: an operational data-assimilative coupled ocean and sea ice ensemble prediction model for the Barents Sea and Svalbard
Open-ocean tides simulated by ICON-O, version icon-2.6.6
Using Probability Density Functions to Evaluate Models (PDFEM, v1.0) to compare a biogeochemical model with satellite-derived chlorophyll
Data assimilation sensitivity experiments in the East Auckland Current system using 4D-Var
Using the COAsT Python package to develop a standardised validation workflow for ocean physics models
Improving Antarctic Bottom Water precursors in NEMO for climate applications
Formulation, optimization, and sensitivity of NitrOMZv1.0, a biogeochemical model of the nitrogen cycle in oceanic oxygen minimum zones
Waves in SKRIPS: WAVEWATCH III coupling implementation and a case study of Tropical Cyclone Mekunu
Adding sea ice effects to a global operational model (NEMO v3.6) for forecasting total water level: approach and impact
Enhanced ocean wave modeling by including effect of breaking under both deep- and shallow-water conditions
An internal solitary wave forecasting model in the northern South China Sea (ISWFM-NSCS)
The 3D biogeochemical marine mercury cycling model MERCY v2.0 – linking atmospheric Hg to methylmercury in fish
Global seamless tidal simulation using a 3D unstructured-grid model (SCHISM v5.10.0)
Arctic Ocean simulations in the CMIP6 Ocean Model Intercomparison Project (OMIP)
ChemicalDrift 1.0: an open-source Lagrangian chemical-fate and transport model for organic aquatic pollutants
The Met Office operational wave forecasting system: the evolution of the regional and global models
Implementation and assessment of a model including mixotrophs and the carbonate cycle (Eco3M_MIX-CarbOx v1.0) in a highly dynamic Mediterranean coastal environment (Bay of Marseille, France) (Part. II): Towards a better representation of total alkalinity when modelling the carbonate system and air-sea CO2 fluxes
4DVarNet-SSH: end-to-end learning of variational interpolation schemes for nadir and wide-swath satellite altimetry
Development and validation of a global 1∕32° surface-wave–tide–circulation coupled ocean model: FIO-COM32
Reproducible and relocatable regional ocean modelling: fundamentals and practices
Barotropic tides in MPAS-Ocean (E3SM V2): impact of ice shelf cavities
Using the two-way nesting technique AGRIF with MARS3D V11.2 to improve hydrodynamics and estimate environmental indicators
Multidecadal and climatological surface current simulations for the southwestern Indian Ocean at 1∕50° resolution
The tidal effects in the Finite-volumE Sea ice–Ocean Model (FESOM2.1): a comparison between parameterised tidal mixing and explicit tidal forcing
HIDRA2: deep-learning ensemble sea level and storm tide forecasting in the presence of seiches – the case of the northern Adriatic
Moana Ocean Hindcast – a > 25-year simulation for New Zealand waters using the Regional Ocean Modeling System (ROMS) v3.9 model
Peter Mlakar, Antonio Ricchi, Sandro Carniel, Davide Bonaldo, and Matjaž Ličer
Geosci. Model Dev., 17, 4705–4725, https://doi.org/10.5194/gmd-17-4705-2024, https://doi.org/10.5194/gmd-17-4705-2024, 2024
Short summary
Short summary
We propose a new point-prediction model, the DEep Learning WAVe Emulating model (DELWAVE), which successfully emulates the Simulating WAves Nearshore model (SWAN) over synoptic to climate timescales. Compared to control climatology over all wind directions, the mismatch between DELWAVE and SWAN is generally small compared to the difference between scenario and control conditions, suggesting that the noise introduced by surrogate modelling is substantially weaker than the climate change signal.
Susanna Winkelbauer, Michael Mayer, and Leopold Haimberger
Geosci. Model Dev., 17, 4603–4620, https://doi.org/10.5194/gmd-17-4603-2024, https://doi.org/10.5194/gmd-17-4603-2024, 2024
Short summary
Short summary
Oceanic transports shape the global climate, but the evaluation and validation of this key quantity based on reanalysis and model data are complicated by the distortion of the used modelling grids and the large number of different grid types. We present two new methods that allow the calculation of oceanic fluxes of volume, heat, salinity, and ice through almost arbitrary sections for various models and reanalyses that are independent of the used modelling grids.
Xiaoyu Fan, Baylor Fox-Kemper, Nobuhiro Suzuki, Qing Li, Patrick Marchesiello, Peter P. Sullivan, and Paul S. Hall
Geosci. Model Dev., 17, 4095–4113, https://doi.org/10.5194/gmd-17-4095-2024, https://doi.org/10.5194/gmd-17-4095-2024, 2024
Short summary
Short summary
Simulations of the oceanic turbulent boundary layer using the nonhydrostatic CROCO ROMS and NCAR-LES models are compared. CROCO and the NCAR-LES are accurate in a similar manner, but CROCO’s additional features (e.g., nesting and realism) and its compressible turbulence formulation carry additional costs.
Jilian Xiong and Parker MacCready
Geosci. Model Dev., 17, 3341–3356, https://doi.org/10.5194/gmd-17-3341-2024, https://doi.org/10.5194/gmd-17-3341-2024, 2024
Short summary
Short summary
The new offline particle tracking package, Tracker v1.1, is introduced to the Regional Ocean Modeling System, featuring an efficient nearest-neighbor algorithm to enhance particle-tracking speed. Its performance was evaluated against four other tracking packages and passive dye. Despite unique features, all packages yield comparable results. Running multiple packages within the same circulation model allows comparison of their performance and ease of use.
Sylvain Cailleau, Laurent Bessières, Léonel Chiendje, Flavie Dubost, Guillaume Reffray, Jean-Michel Lellouche, Simon van Gennip, Charly Régnier, Marie Drevillon, Marc Tressol, Matthieu Clavier, Julien Temple-Boyer, and Léo Berline
Geosci. Model Dev., 17, 3157–3173, https://doi.org/10.5194/gmd-17-3157-2024, https://doi.org/10.5194/gmd-17-3157-2024, 2024
Short summary
Short summary
In order to improve Sargassum drift forecasting in the Caribbean area, drift models can be forced by higher-resolution ocean currents. To this goal a 3 km resolution regional ocean model has been developed. Its assessment is presented with a particular focus on the reproduction of fine structures representing key features of the Caribbean region dynamics and Sargassum transport. The simulated propagation of a North Brazil Current eddy and its dissipation was found to be quite realistic.
Gaetano Porcile, Anne-Claire Bennis, Martial Boutet, Sophie Le Bot, Franck Dumas, and Swen Jullien
Geosci. Model Dev., 17, 2829–2853, https://doi.org/10.5194/gmd-17-2829-2024, https://doi.org/10.5194/gmd-17-2829-2024, 2024
Short summary
Short summary
Here a new method of modelling the interaction between ocean currents and waves is presented. We developed an advanced coupling of two models, one for ocean currents and one for waves. In previous couplings, some wave-related calculations were based on simplified assumptions. Our method uses more complex calculations to better represent wave–current interactions. We tested it in a macro-tidal coastal area and found that it significantly improves the model accuracy, especially during storms.
Colette Gabrielle Kerry, Moninya Roughan, Shane Keating, David Gwyther, Gary Brassington, Adil Siripatana, and Joao Marcos A. C. Souza
Geosci. Model Dev., 17, 2359–2386, https://doi.org/10.5194/gmd-17-2359-2024, https://doi.org/10.5194/gmd-17-2359-2024, 2024
Short summary
Short summary
Ocean forecasting relies on the combination of numerical models and ocean observations through data assimilation (DA). Here we assess the performance of two DA systems in a dynamic western boundary current, the East Australian Current, across a common modelling and observational framework. We show that the more advanced, time-dependent method outperforms the time-independent method for forecast horizons of 5 d. This advocates the use of advanced methods for highly variable oceanic regions.
Ivan Hernandez, Leidy M. Castro-Rosero, Manuel Espino, and Jose M. Alsina Torrent
Geosci. Model Dev., 17, 2221–2245, https://doi.org/10.5194/gmd-17-2221-2024, https://doi.org/10.5194/gmd-17-2221-2024, 2024
Short summary
Short summary
The LOCATE numerical model was developed to conduct Lagrangian simulations of the transport and dispersion of marine debris at coastal scales. High-resolution hydrodynamic data and a beaching module that used particle distance to the shore for land–water boundary detection were used on a realistic debris discharge scenario comparing hydrodynamic data at various resolutions. Coastal processes and complex geometric structures were resolved when using nested grids and distance-to-shore beaching.
Neill Mackay, Jan Zika, Taimoor Sohail, Richard Williams, Oliver Andrews, and Andrew Watson
EGUsphere, https://doi.org/10.5194/egusphere-2023-2448, https://doi.org/10.5194/egusphere-2023-2448, 2024
Short summary
Short summary
The ocean absorbs carbon dioxide from the atmosphere, mitigating climate change, but estimates of the uptake do not always agree. There is a need to reconcile these differing estimates, and to improve our understanding of ocean carbon uptake. We present a new method for estimating ocean carbon uptake and test it with model data. The method effectively diagnoses the ocean carbon uptake from limited data, and therefore shows promise for reconciling different observational estimates.
Ngoc B. Trinh, Marine Herrmann, Caroline Ulses, Patrick Marsaleix, Thomas Duhaut, Thai To Duy, Claude Estournel, and R. Kipp Shearman
Geosci. Model Dev., 17, 1831–1867, https://doi.org/10.5194/gmd-17-1831-2024, https://doi.org/10.5194/gmd-17-1831-2024, 2024
Short summary
Short summary
A high-resolution model was built to study the South China Sea (SCS) water, heat, and salt budgets. Model performance is demonstrated by comparison with observations and simulations. Important discards are observed if calculating offline, instead of online, lateral inflows and outflows of water, heat, and salt. The SCS mainly receives water from the Luzon Strait and releases it through the Mindoro, Taiwan, and Karimata straits. SCS surface interocean water exchanges are driven by monsoon winds.
Louis Thiry, Long Li, Guillaume Roullet, and Etienne Mémin
Geosci. Model Dev., 17, 1749–1764, https://doi.org/10.5194/gmd-17-1749-2024, https://doi.org/10.5194/gmd-17-1749-2024, 2024
Short summary
Short summary
We present a new way of solving the quasi-geostrophic (QG) equations, a simple set of equations describing ocean dynamics. Our method is solely based on the numerical methods used to solve the equations and requires no parameter tuning. Moreover, it can handle non-rectangular geometries, opening the way to study QG equations on realistic domains. We release a PyTorch implementation to ease future machine-learning developments on top of the presented method.
Zheqi Shen, Yihao Chen, Xiaojing Li, and Xunshu Song
Geosci. Model Dev., 17, 1651–1665, https://doi.org/10.5194/gmd-17-1651-2024, https://doi.org/10.5194/gmd-17-1651-2024, 2024
Short summary
Short summary
Parameter estimation is the process that optimizes model parameters using observations, which could reduce model errors and improve forecasting. In this study, we conducted parameter estimation experiments using the CESM and the ensemble adjustment Kalman filter. The obtained initial conditions and parameters are used to perform ensemble forecast experiments for ENSO forecasting. The results revealed that parameter estimation could reduce analysis errors and improve ENSO forecast skills.
Ye Yuan, Fujiang Yu, Zhi Chen, Xueding Li, Fang Hou, Yuanyong Gao, Zhiyi Gao, and Renbo Pang
EGUsphere, https://doi.org/10.5194/egusphere-2024-169, https://doi.org/10.5194/egusphere-2024-169, 2024
Short summary
Short summary
Accurate and timely forecasting of ocean waves is of great importance to the safety of marine transportation and offshore engineering. In this study, Graphics Processing Unit (GPU)-accelerated computing is introduced in the WAM (Cycle 6). With this effort, global high-resolution wave simulations now can run on GPUs up to tens of times faster than currently available models on a CPU node with just as accurate results.
Ali Abdolali, Saeideh Banihashemi, Jose Henrique Alves, Aron Roland, Tyler J. Hesser, Mary Anderson Bryant, and Jane McKee Smith
Geosci. Model Dev., 17, 1023–1039, https://doi.org/10.5194/gmd-17-1023-2024, https://doi.org/10.5194/gmd-17-1023-2024, 2024
Short summary
Short summary
This article presents an overview of the development and implementation of Great Lake Wave Unstructured (GLWUv2.0), including the core model and workflow design and development. The validation was conducted against in situ data for the re-forecasted duration for summer and wintertime (ice season). The article describes the limitations and challenges encountered in the operational environment and the path forward for the next generation of wave forecast systems in enclosed basins like the GL.
Xuanxuan Gao, Shuiqing Li, Dongxue Mo, Yahao Liu, and Po Hu
Geosci. Model Dev. Discuss., https://doi.org/10.5194/gmd-2024-12, https://doi.org/10.5194/gmd-2024-12, 2024
Revised manuscript accepted for GMD
Short summary
Short summary
Storm surges generate coastal inundation and expose populations and properties in danger. We developed a novel storm surge inundation model for efficient prediction. Estimates compare well with in-situ measurements and results from a numerical model. The new model significantly improves over the existing numerical models with much higher computational efficiency and stability, which allows timely disaster prevention and mitigation.
Vincenzo de Toma, Daniele Ciani, Yassmin Hesham Essa, Chunxue Yang, Vincenzo Artale, Andrea Pisano, Davide Cavaliere, Rosalia Santoleri, and Andrea Storto
Geosci. Model Dev. Discuss., https://doi.org/10.5194/gmd-2024-13, https://doi.org/10.5194/gmd-2024-13, 2024
Revised manuscript accepted for GMD
Short summary
Short summary
This study explores methods to reconstruct diurnal variations in skin sea surface temperature in a model of the Mediterranean Sea. Our new approach, considering chlorophyll concentration, enhances spatial and temporal variations in the warm layer. Comparative analysis shows context-dependent improvements. The proposed "chlorophyll-interactive" method brings the surface net total heat flux closer to zero annually, despite a net heat loss from the ocean to the atmosphere.
Krysten Rutherford, Laura Bianucci, and William Floyd
EGUsphere, https://doi.org/10.5194/egusphere-2023-3064, https://doi.org/10.5194/egusphere-2023-3064, 2024
Short summary
Short summary
Nearshore ocean models often lack complete information about freshwater fluxes due to numerous ungauged rivers and streams. We tested a simple rain-based hydrological model as inputs into an ocean model of Quatsino Sound, B.C., Canada with the aim of improving the representation of the land-ocean connection in the nearshore model. Through multiple tests, we found that the performance of the ocean model improved when providing 60 % or more of the freshwater inputs from the simple runoff model.
Qiang Wang, Qi Shu, Alexandra Bozec, Eric P. Chassignet, Pier Giuseppe Fogli, Baylor Fox-Kemper, Andy McC. Hogg, Doroteaciro Iovino, Andrew E. Kiss, Nikolay Koldunov, Julien Le Sommer, Yiwen Li, Pengfei Lin, Hailong Liu, Igor Polyakov, Patrick Scholz, Dmitry Sidorenko, Shizhu Wang, and Xiaobiao Xu
Geosci. Model Dev., 17, 347–379, https://doi.org/10.5194/gmd-17-347-2024, https://doi.org/10.5194/gmd-17-347-2024, 2024
Short summary
Short summary
Increasing resolution improves model skills in simulating the Arctic Ocean, but other factors such as parameterizations and numerics are at least of the same importance for obtaining reliable simulations.
Andrew C. Ross, Charles A. Stock, Alistair Adcroft, Enrique Curchitser, Robert Hallberg, Matthew J. Harrison, Katherine Hedstrom, Niki Zadeh, Michael Alexander, Wenhao Chen, Elizabeth J. Drenkard, Hubert du Pontavice, Raphael Dussin, Fabian Gomez, Jasmin G. John, Dujuan Kang, Diane Lavoie, Laure Resplandy, Alizée Roobaert, Vincent Saba, Sang-Ik Shin, Samantha Siedlecki, and James Simkins
Geosci. Model Dev., 16, 6943–6985, https://doi.org/10.5194/gmd-16-6943-2023, https://doi.org/10.5194/gmd-16-6943-2023, 2023
Short summary
Short summary
We evaluate a model for northwest Atlantic Ocean dynamics and biogeochemistry that balances high resolution with computational economy by building on the new regional features in the MOM6 ocean model and COBALT biogeochemical model. We test the model's ability to simulate impactful historical variability and find that the model simulates the mean state and variability of most features well, which suggests the model can provide information to inform living-marine-resource applications.
Luca Arpaia, Christian Ferrarin, Marco Bajo, and Georg Umgiesser
Geosci. Model Dev., 16, 6899–6919, https://doi.org/10.5194/gmd-16-6899-2023, https://doi.org/10.5194/gmd-16-6899-2023, 2023
Short summary
Short summary
We propose a discrete multilayer shallow water model based on z-layers which, thanks to the insertion and removal of surface layers, can deal with an arbitrarily large tidal oscillation independently of the vertical resolution. The algorithm is based on a two-step procedure used in numerical simulations with moving boundaries (grid movement followed by a grid topology change, that is, the insertion/removal of surface layers), which avoids the appearance of very thin surface layers.
Lucille Barré, Frédéric Diaz, Thibaut Wagener, France Van Wambeke, Camille Mazoyer, Christophe Yohia, and Christel Pinazo
Geosci. Model Dev., 16, 6701–6739, https://doi.org/10.5194/gmd-16-6701-2023, https://doi.org/10.5194/gmd-16-6701-2023, 2023
Short summary
Short summary
While several studies have shown that mixotrophs play a crucial role in the carbon cycle, the impact of environmental forcings on their dynamics remains poorly investigated. Using a biogeochemical model that considers mixotrophs, we study the impact of light and nutrient concentration on the ecosystem composition in a highly dynamic Mediterranean coastal area: the Bay of Marseille. We show that mixotrophs cope better with oligotrophic conditions compared to strict auto- and heterotrophs.
Trygve Halsne, Kai Håkon Christensen, Gaute Hope, and Øyvind Breivik
Geosci. Model Dev., 16, 6515–6530, https://doi.org/10.5194/gmd-16-6515-2023, https://doi.org/10.5194/gmd-16-6515-2023, 2023
Short summary
Short summary
Surface waves that propagate in oceanic or coastal environments get influenced by their surroundings. Changes in the ambient current or the depth profile affect the wave propagation path, and the change in wave direction is called refraction. Some analytical solutions to the governing equations exist under ideal conditions, but for realistic situations, the equations must be solved numerically. Here we present such a numerical solver under an open-source license.
Jiangyu Li, Shaoqing Zhang, Qingxiang Liu, Xiaolin Yu, and Zhiwei Zhang
Geosci. Model Dev., 16, 6393–6412, https://doi.org/10.5194/gmd-16-6393-2023, https://doi.org/10.5194/gmd-16-6393-2023, 2023
Short summary
Short summary
Ocean surface waves play an important role in the air–sea interface but are rarely activated in high-resolution Earth system simulations due to their expensive computational costs. To alleviate this situation, this paper designs a new wave modeling framework with a multiscale grid system. Evaluations of a series of numerical experiments show that it has good feasibility and applicability in the WAVEWATCH III model, WW3, and can achieve the goals of efficient and high-precision wave simulation.
Doroteaciro Iovino, Pier Giuseppe Fogli, and Simona Masina
Geosci. Model Dev., 16, 6127–6159, https://doi.org/10.5194/gmd-16-6127-2023, https://doi.org/10.5194/gmd-16-6127-2023, 2023
Short summary
Short summary
This paper describes the model performance of three global ocean–sea ice configurations, from non-eddying (1°) to eddy-rich (1/16°) resolutions. Model simulations are obtained following the Ocean Model Intercomparison Project phase 2 (OMIP2) protocol. We compare key global climate variables across the three models and against observations, emphasizing the relative advantages and disadvantages of running forced ocean–sea ice models at higher resolution.
Johannes Röhrs, Yvonne Gusdal, Edel S. U. Rikardsen, Marina Durán Moro, Jostein Brændshøi, Nils Melsom Kristensen, Sindre Fritzner, Keguang Wang, Ann Kristin Sperrevik, Martina Idžanović, Thomas Lavergne, Jens Boldingh Debernard, and Kai H. Christensen
Geosci. Model Dev., 16, 5401–5426, https://doi.org/10.5194/gmd-16-5401-2023, https://doi.org/10.5194/gmd-16-5401-2023, 2023
Short summary
Short summary
A model to predict ocean currents, temperature, and sea ice is presented, covering the Barents Sea and northern Norway. To quantify forecast uncertainties, the model calculates ensemble forecasts with 24 realizations of ocean and ice conditions. Observations from satellites, buoys, and ships are ingested by the model. The model forecasts are compared with observations, and we show that the ocean model has skill in predicting sea surface temperatures.
Jin-Song von Storch, Eileen Hertwig, Veit Lüschow, Nils Brüggemann, Helmuth Haak, Peter Korn, and Vikram Singh
Geosci. Model Dev., 16, 5179–5196, https://doi.org/10.5194/gmd-16-5179-2023, https://doi.org/10.5194/gmd-16-5179-2023, 2023
Short summary
Short summary
The new ocean general circulation model ICON-O is developed for running experiments at kilometer scales and beyond. One targeted application is to simulate internal tides crucial for ocean mixing. To ensure their realism, which is difficult to assess, we evaluate the barotropic tides that generate internal tides. We show that ICON-O is able to realistically simulate the major aspects of the observed barotropic tides and discuss the aspects that impact the quality of the simulated tides.
Bror F. Jönsson, Christopher L. Follett, Jacob Bien, Stephanie Dutkiewicz, Sangwon Hyun, Gemma Kulk, Gael L. Forget, Christian Müller, Marie-Fanny Racault, Christopher N. Hill, Thomas Jackson, and Shubha Sathyendranath
Geosci. Model Dev., 16, 4639–4657, https://doi.org/10.5194/gmd-16-4639-2023, https://doi.org/10.5194/gmd-16-4639-2023, 2023
Short summary
Short summary
While biogeochemical models and satellite-derived ocean color data provide unprecedented information, it is problematic to compare them. Here, we present a new approach based on comparing probability density distributions of model and satellite properties to assess model skills. We also introduce Earth mover's distances as a novel and powerful metric to quantify the misfit between models and observations. We find that how 3D chlorophyll fields are aggregated can be a significant source of error.
Rafael Santana, Helen Macdonald, Joanne O'Callaghan, Brian Powell, Sarah Wakes, and Sutara H. Suanda
Geosci. Model Dev., 16, 3675–3698, https://doi.org/10.5194/gmd-16-3675-2023, https://doi.org/10.5194/gmd-16-3675-2023, 2023
Short summary
Short summary
We show the importance of assimilating subsurface temperature and velocity data in a model of the East Auckland Current. Assimilation of velocity increased the representation of large oceanic vortexes. Assimilation of temperature is needed to correctly simulate temperatures around 100 m depth, which is the most difficult region to simulate in ocean models. Our simulations showed improved results in comparison to the US Navy global model and highlight the importance of regional models.
David Byrne, Jeff Polton, Enda O'Dea, and Joanne Williams
Geosci. Model Dev., 16, 3749–3764, https://doi.org/10.5194/gmd-16-3749-2023, https://doi.org/10.5194/gmd-16-3749-2023, 2023
Short summary
Short summary
Validation is a crucial step during the development of models for ocean simulation. The purpose of validation is to assess how accurate a model is. It is most commonly done by comparing output from a model to actual observations. In this paper, we introduce and demonstrate usage of the COAsT Python package to standardise the validation process for physical ocean models. We also discuss our five guiding principles for standardised validation.
Katherine Hutchinson, Julie Deshayes, Christian Éthé, Clément Rousset, Casimir de Lavergne, Martin Vancoppenolle, Nicolas C. Jourdain, and Pierre Mathiot
Geosci. Model Dev., 16, 3629–3650, https://doi.org/10.5194/gmd-16-3629-2023, https://doi.org/10.5194/gmd-16-3629-2023, 2023
Short summary
Short summary
Bottom Water constitutes the lower half of the ocean’s overturning system and is primarily formed in the Weddell and Ross Sea in the Antarctic due to interactions between the atmosphere, ocean, sea ice and ice shelves. Here we use a global ocean 1° resolution model with explicit representation of the three large ice shelves important for the formation of the parent waters of Bottom Water. We find doing so reduces salt biases, improves water mass realism and gives realistic ice shelf melt rates.
Daniele Bianchi, Daniel McCoy, and Simon Yang
Geosci. Model Dev., 16, 3581–3609, https://doi.org/10.5194/gmd-16-3581-2023, https://doi.org/10.5194/gmd-16-3581-2023, 2023
Short summary
Short summary
We present NitrOMZ, a new model of the oceanic nitrogen cycle that simulates chemical transformations within oxygen minimum zones (OMZs). We describe the model formulation and its implementation in a one-dimensional representation of the water column before evaluating its ability to reproduce observations in the eastern tropical South Pacific. We conclude by describing the model sensitivity to parameter choices and environmental factors and its application to nitrogen cycling in the ocean.
Rui Sun, Alison Cobb, Ana B. Villas Bôas, Sabique Langodan, Aneesh C. Subramanian, Matthew R. Mazloff, Bruce D. Cornuelle, Arthur J. Miller, Raju Pathak, and Ibrahim Hoteit
Geosci. Model Dev., 16, 3435–3458, https://doi.org/10.5194/gmd-16-3435-2023, https://doi.org/10.5194/gmd-16-3435-2023, 2023
Short summary
Short summary
In this work, we integrated the WAVEWATCH III model into the regional coupled model SKRIPS. We then performed a case study using the newly implemented model to study Tropical Cyclone Mekunu, which occurred in the Arabian Sea. We found that the coupled model better simulates the cyclone than the uncoupled model, but the impact of waves on the cyclone is not significant. However, the waves change the sea surface temperature and mixed layer, especially in the cold waves produced due to the cyclone.
Pengcheng Wang and Natacha B. Bernier
Geosci. Model Dev., 16, 3335–3354, https://doi.org/10.5194/gmd-16-3335-2023, https://doi.org/10.5194/gmd-16-3335-2023, 2023
Short summary
Short summary
Effects of sea ice are typically neglected in operational flood forecast systems. In this work, we capture these effects via the addition of a parameterized ice–ocean stress. The parameterization takes advantage of forecast fields from an advanced ice–ocean model and features a novel, consistent representation of the tidal relative ice–ocean velocity. The new parameterization leads to improved forecasts of tides and storm surges in polar regions. Associated physical processes are discussed.
Yue Xu and Xiping Yu
Geosci. Model Dev., 16, 2811–2831, https://doi.org/10.5194/gmd-16-2811-2023, https://doi.org/10.5194/gmd-16-2811-2023, 2023
Short summary
Short summary
An accurate description of the wind energy input into ocean waves is crucial to ocean wave modeling, and a physics-based consideration of the effect of wave breaking is absolutely necessary to obtain such an accurate description, particularly under extreme conditions. This study evaluates the performance of a recently improved formula, taking into account not only the effect of breaking but also the effect of airflow separation on the leeside of steep wave crests in a reasonably consistent way.
Yankun Gong, Xueen Chen, Jiexin Xu, Jieshuo Xie, Zhiwu Chen, Yinghui He, and Shuqun Cai
Geosci. Model Dev., 16, 2851–2871, https://doi.org/10.5194/gmd-16-2851-2023, https://doi.org/10.5194/gmd-16-2851-2023, 2023
Short summary
Short summary
Internal solitary waves (ISWs) play crucial roles in mass transport and ocean mixing in the northern South China Sea. Massive numerical investigations have been conducted in this region, but there was no systematic evaluation of a three-dimensional model about precisely simulating ISWs. Here, an ISW forecasting model is employed to evaluate the roles of resolution, tidal forcing and stratification in accurately reproducing wave properties via comparison to field and remote-sensing observations.
Johannes Bieser, David J. Amptmeijer, Ute Daewel, Joachim Kuss, Anne L. Soerensen, and Corinna Schrum
Geosci. Model Dev., 16, 2649–2688, https://doi.org/10.5194/gmd-16-2649-2023, https://doi.org/10.5194/gmd-16-2649-2023, 2023
Short summary
Short summary
MERCY is a 3D model to study mercury (Hg) cycling in the ocean. Hg is a highly harmful pollutant regulated by the UN Minamata Convention on Mercury due to widespread human emissions. These emissions eventually reach the oceans, where Hg transforms into the even more toxic and bioaccumulative pollutant methylmercury. MERCY predicts the fate of Hg in the ocean and its buildup in the food chain. It is the first model to consider Hg accumulation in fish, a major source of Hg exposure for humans.
Y. Joseph Zhang, Tomas Fernandez-Montblanc, William Pringle, Hao-Cheng Yu, Linlin Cui, and Saeed Moghimi
Geosci. Model Dev., 16, 2565–2581, https://doi.org/10.5194/gmd-16-2565-2023, https://doi.org/10.5194/gmd-16-2565-2023, 2023
Short summary
Short summary
Simulating global ocean from deep basins to coastal areas is a daunting task but is important for disaster mitigation efforts. We present a new 3D global ocean model on flexible mesh to study both tidal and nontidal processes and total water prediction. We demonstrate the potential for
seamlesssimulation, on a single mesh, from the global ocean to a few estuaries along the US West Coast. The model can serve as the backbone of a global tide surge and compound flooding forecasting framework.
Qi Shu, Qiang Wang, Chuncheng Guo, Zhenya Song, Shizhu Wang, Yan He, and Fangli Qiao
Geosci. Model Dev., 16, 2539–2563, https://doi.org/10.5194/gmd-16-2539-2023, https://doi.org/10.5194/gmd-16-2539-2023, 2023
Short summary
Short summary
Ocean models are often used for scientific studies on the Arctic Ocean. Here the Arctic Ocean simulations by state-of-the-art global ocean–sea-ice models participating in the Ocean Model Intercomparison Project (OMIP) were evaluated. The simulations on Arctic Ocean hydrography, freshwater content, stratification, sea surface height, and gateway transports were assessed and the common biases were detected. The simulations forced by different atmospheric forcing were also evaluated.
Manuel Aghito, Loris Calgaro, Knut-Frode Dagestad, Christian Ferrarin, Antonio Marcomini, Øyvind Breivik, and Lars Robert Hole
Geosci. Model Dev., 16, 2477–2494, https://doi.org/10.5194/gmd-16-2477-2023, https://doi.org/10.5194/gmd-16-2477-2023, 2023
Short summary
Short summary
The newly developed ChemicalDrift model can simulate the transport and fate of chemicals in the ocean and in coastal regions. The model combines ocean physics, including transport due to currents, turbulence due to surface winds and the sinking of particles to the sea floor, with ocean chemistry, such as the partitioning, the degradation and the evaporation of chemicals. The model will be utilized for risk assessment of ocean and sea-floor contamination from pollutants emitted from shipping.
Nieves G. Valiente, Andrew Saulter, Breogan Gomez, Christopher Bunney, Jian-Guo Li, Tamzin Palmer, and Christine Pequignet
Geosci. Model Dev., 16, 2515–2538, https://doi.org/10.5194/gmd-16-2515-2023, https://doi.org/10.5194/gmd-16-2515-2023, 2023
Short summary
Short summary
We document the Met Office operational global and regional wave models which provide wave forecasts up to 7 d ahead. Our models present coarser resolution offshore to higher resolution near the coastline. The increased resolution led to replication of the extremes but to some overestimation during modal conditions. If currents are included, wave directions and long period swells near the coast are significantly improved. New developments focus on the optimisation of the models with resolution.
Lucille Barré, Frédéric Diaz, Thibaut Wagener, Camille Mazoyer, Christophe Yohia, and Christel Pinazo
Geosci. Model Dev. Discuss., https://doi.org/10.5194/gmd-2023-34, https://doi.org/10.5194/gmd-2023-34, 2023
Revised manuscript accepted for GMD
Short summary
Short summary
Carbonate system is typically studied using measurements, yet modelling can contribute valuable insights. Using a biogeochemical model, we propose a new representation of total alkalinity, dissolved inorganic carbon, pCO2 and pH in a highly dynamic Mediterranean coastal area: the Bay of Marseille, a useful addition to measurements. Through a detailed analysis of pCO2 and air-sea CO2 fluxes we show that their variations are strongly impacted by the hydrodynamic processes that affect the bay.
Maxime Beauchamp, Quentin Febvre, Hugo Georgenthum, and Ronan Fablet
Geosci. Model Dev., 16, 2119–2147, https://doi.org/10.5194/gmd-16-2119-2023, https://doi.org/10.5194/gmd-16-2119-2023, 2023
Short summary
Short summary
4DVarNet is a learning-based method based on traditional data assimilation (DA). This new class of algorithms can be used to provide efficient reconstructions of a dynamical system based on single observations. We provide a 4DVarNet application to sea surface height reconstructions based on nadir and future Surface Water and Ocean and Topography data. It outperforms other methods, from optimal interpolation to sophisticated DA algorithms. This work is part of on-going AI Chair Oceanix projects.
Bin Xiao, Fangli Qiao, Qi Shu, Xunqiang Yin, Guansuo Wang, and Shihong Wang
Geosci. Model Dev., 16, 1755–1777, https://doi.org/10.5194/gmd-16-1755-2023, https://doi.org/10.5194/gmd-16-1755-2023, 2023
Short summary
Short summary
A new global surface-wave–tide–circulation coupled ocean model (FIO-COM32) with a resolution of 1/32° × 1/32° is developed and validated. Both the promotion of the horizontal resolution and included physical processes are shown to be important contributors to the significant improvements in FIO-COM32 simulations. It is time to merge these separated model components (surface waves, tidal currents and ocean circulation) and start a new generation of ocean model development.
Jeff Polton, James Harle, Jason Holt, Anna Katavouta, Dale Partridge, Jenny Jardine, Sarah Wakelin, Julia Rulent, Anthony Wise, Katherine Hutchinson, David Byrne, Diego Bruciaferri, Enda O'Dea, Michela De Dominicis, Pierre Mathiot, Andrew Coward, Andrew Yool, Julien Palmiéri, Gennadi Lessin, Claudia Gabriela Mayorga-Adame, Valérie Le Guennec, Alex Arnold, and Clément Rousset
Geosci. Model Dev., 16, 1481–1510, https://doi.org/10.5194/gmd-16-1481-2023, https://doi.org/10.5194/gmd-16-1481-2023, 2023
Short summary
Short summary
The aim is to increase the capacity of the modelling community to respond to societally important questions that require ocean modelling. The concept of reproducibility for regional ocean modelling is developed: advocating methods for reproducible workflows and standardised methods of assessment. Then, targeting the NEMO framework, we give practical advice and worked examples, highlighting key considerations that will the expedite development cycle and upskill the user community.
Nairita Pal, Kristin N. Barton, Mark R. Petersen, Steven R. Brus, Darren Engwirda, Brian K. Arbic, Andrew F. Roberts, Joannes J. Westerink, and Damrongsak Wirasaet
Geosci. Model Dev., 16, 1297–1314, https://doi.org/10.5194/gmd-16-1297-2023, https://doi.org/10.5194/gmd-16-1297-2023, 2023
Short summary
Short summary
Understanding tides is essential to accurately predict ocean currents. Over the next several decades coastal processes such as flooding and erosion will be severely impacted due to climate change. Tides affect currents along the coastal regions the most. In this paper we show the results of implementing tides in a global ocean model known as MPAS–Ocean. We also show how Antarctic ice shelf cavities affect global tides. Our work points towards future research with tide–ice interactions.
Sébastien Petton, Valérie Garnier, Matthieu Caillaud, Laurent Debreu, and Franck Dumas
Geosci. Model Dev., 16, 1191–1211, https://doi.org/10.5194/gmd-16-1191-2023, https://doi.org/10.5194/gmd-16-1191-2023, 2023
Short summary
Short summary
The nesting AGRIF library is implemented in the MARS3D hydrodynamic model, a semi-implicit, free-surface numerical model which uses a time scheme as an alternating-direction implicit (ADI) algorithm. Two applications at the regional and coastal scale are introduced. We compare the two-nesting approach to the classic offline one-way approach, based on an in situ dataset. This method is an efficient means to significantly improve the physical hydrodynamics and unravel ecological challenges.
Noam S. Vogt-Vincent and Helen L. Johnson
Geosci. Model Dev., 16, 1163–1178, https://doi.org/10.5194/gmd-16-1163-2023, https://doi.org/10.5194/gmd-16-1163-2023, 2023
Short summary
Short summary
Ocean currents transport things over large distances across the ocean surface. Predicting this transport is key for tackling many environmental problems, such as marine plastic pollution and coral reef resilience. However, doing this requires a good understanding ocean currents, which is currently lacking. Here, we present and validate state-of-the-art simulations for surface currents in the southwestern Indian Ocean, which will support future marine dispersal studies across this region.
Pengyang Song, Dmitry Sidorenko, Patrick Scholz, Maik Thomas, and Gerrit Lohmann
Geosci. Model Dev., 16, 383–405, https://doi.org/10.5194/gmd-16-383-2023, https://doi.org/10.5194/gmd-16-383-2023, 2023
Short summary
Short summary
Tides have essential effects on the ocean and climate. Most previous research applies parameterised tidal mixing to discuss their effects in models. By comparing the effect of a tidal mixing parameterisation and tidal forcing on the ocean state, we assess the advantages and disadvantages of the two methods. Our results show that tidal mixing in the North Pacific Ocean strongly affects the global thermohaline circulation. We also list some effects that are not considered in the parameterisation.
Marko Rus, Anja Fettich, Matej Kristan, and Matjaž Ličer
Geosci. Model Dev., 16, 271–288, https://doi.org/10.5194/gmd-16-271-2023, https://doi.org/10.5194/gmd-16-271-2023, 2023
Short summary
Short summary
We propose a new fast and reliable deep-learning architecture HIDRA2 for sea level and storm surge modeling. HIDRA2 features new feature encoders and a fusion-regression block. We test HIDRA2 on Adriatic storm surges, which depend on an interaction between tides and seiches. We demonstrate that HIDRA2 learns to effectively mimic the timing and amplitude of Adriatic seiches. This is essential for reliable HIDRA2 predictions of total storm surge sea levels.
Joao Marcos Azevedo Correia de Souza, Sutara H. Suanda, Phellipe P. Couto, Robert O. Smith, Colette Kerry, and Moninya Roughan
Geosci. Model Dev., 16, 211–231, https://doi.org/10.5194/gmd-16-211-2023, https://doi.org/10.5194/gmd-16-211-2023, 2023
Short summary
Short summary
The current paper describes the configuration and evaluation of the Moana Ocean Hindcast, a > 25-year simulation of the ocean state around New Zealand using the Regional Ocean Modeling System v3.9. This is the first open-access, long-term, continuous, realistic ocean simulation for this region and provides information for improving the understanding of the ocean processes that affect the New Zealand exclusive economic zone.
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...
