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.
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
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
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
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 I): Evolution of ecosystem composition under limited light and nutrient conditions
Evaluation of the CMCC global eddying ocean model for the Ocean Model Intercomparison Project (OMIP2)
Reproducible and relocatable regional ocean modelling: fundamentals and practices
Barents-2.5km v2.0: An operational data-assimilative coupled ocean and sea ice ensemble prediction model for the Barents Sea and Svalbard
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
A nonhydrostatic oceanic regional model, ORCTM v1, for internal solitary wave simulation
How does 4DVar data assimilation affect the vertical representation of mesoscale eddies? A case study with observing system simulation experiments (OSSEs) using ROMS v3.9
An ensemble Kalman filter-based ocean data assimilation system improved by adaptive observation error inflation (AOEI)
GULF18, a high-resolution NEMO-based tidal ocean model of the Arabian/Persian Gulf
The Baltic Sea Model Intercomparison Project (BMIP) – a platform for model development, evaluation, and uncertainty assessment
An ensemble Kalman filter system with the Stony Brook Parallel Ocean Model v1.0
Wind work at the air-sea interface: a modeling study in anticipation of future space missions
Improved upper-ocean thermodynamical structure modeling with combined effects of surface waves and M2 internal tides on vertical mixing: a case study for the Indian Ocean
The bulk parameterizations of turbulent air–sea fluxes in NEMO4: the origin of sea surface temperature differences in a global model study
NeverWorld2: an idealized model hierarchy to investigate ocean mesoscale eddies across resolutions
Observing system simulation experiments reveal that subsurface temperature observations improve estimates of circulation and heat content in a dynamic western boundary current
Parallel implementation of the SHYFEM (System of HydrodYnamic Finite Element Modules) model
Block-structured, equal-workload, multi-grid-nesting interface for the Boussinesq wave model FUNWAVE-TVD (Total Variation Diminishing)
Evaluation of an emergent feature of sub-shelf melt oscillations from an idealized coupled ice sheet–ocean model using FISOC (v1.1) – ROMSIceShelf (v1.0) – Elmer/Ice (v9.0)
GNOM v1.0: an optimized steady-state model of the modern marine neodymium cycle
Implementation and evaluation of open boundary conditions for sea ice in a regional coupled ocean (ROMS) and sea ice (CICE) modeling system
ROMSPath v1.0: offline particle tracking for the Regional Ocean Modeling System (ROMS)
DINCAE 2.0: multivariate convolutional neural network with error estimates to reconstruct sea surface temperature satellite and altimetry observations
RADIv1: a non-steady-state early diagenetic model for ocean sediments in Julia and MATLAB/GNU Octave
IBI-CCS: a regional high-resolution model to simulate sea level in western Europe
Empirical Lagrangian parametrization for wind-driven mixing of buoyant particles at the ocean surface
Improving ocean modeling software NEMO 4.0 benchmarking and communication efficiency
Improvements in the regional South China Sea Operational Oceanography Forecasting System (SCSOFSv2)
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.
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.
Lucille Barré, Frédéric Diaz, Thibaut Wagener, France Van Wambeke, Camille Mazoyer, Christophe Yohia, and Christel Pinazo
Geosci. Model Dev. Discuss., https://doi.org/10.5194/gmd-2023-33, https://doi.org/10.5194/gmd-2023-33, 2023
Revised manuscript accepted for GMD
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.
Doroteaciro Iovino, Pier Giuseppe Fogli, and Simona Masina
EGUsphere, https://doi.org/10.5194/egusphere-2023-469, https://doi.org/10.5194/egusphere-2023-469, 2023
Short summary
Short summary
The 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 OMIP-2 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.
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.
Johannes Röhrs, Yvonne Gusdal, Edel Rikardsen, Marina Duran Moro, Jostein Brændshøi, Nils Melsom Kristensen, Sindre Fritzner, Keguang Wang, Ann Kristin Sperrevik, Martina Idžanović, Thomas Lavergne, Jens Debernard, and Kai H. Christensen
Geosci. Model Dev. Discuss., https://doi.org/10.5194/gmd-2023-20, https://doi.org/10.5194/gmd-2023-20, 2023
Revised manuscript accepted for GMD
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 – 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.
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.
Hao Huang, Pengyang Song, Shi Qiu, Jiaqi Guo, and Xueen Chen
Geosci. Model Dev., 16, 109–133, https://doi.org/10.5194/gmd-16-109-2023, https://doi.org/10.5194/gmd-16-109-2023, 2023
Short summary
Short summary
The Oceanic Regional Circulation and Tide Model (ORCTM) is developed to reproduce internal solitary wave dynamics. The three-dimensional nonlinear momentum equations are involved with the nonhydrostatic pressure obtained via solving the Poisson equation. The validation experimental results agree with the internal wave theories and observations, demonstrating that the ORCTM can successfully describe the life cycle of nonlinear internal solitary waves under different oceanic environments.
David E. Gwyther, Shane R. Keating, Colette Kerry, and Moninya Roughan
Geosci. Model Dev., 16, 157–178, https://doi.org/10.5194/gmd-16-157-2023, https://doi.org/10.5194/gmd-16-157-2023, 2023
Short summary
Short summary
Ocean eddies are important for weather, climate, biology, navigation, and search and rescue. Since eddies change rapidly, models that incorporate or assimilate observations are required to produce accurate eddy timings and locations, yet the model accuracy is rarely assessed below the surface. We use a unique type of ocean model experiment to assess three-dimensional eddy structure in the East Australian Current and explore two pathways in which this subsurface structure is being degraded.
Shun Ohishi, Takemasa Miyoshi, and Misako Kachi
Geosci. Model Dev., 15, 9057–9073, https://doi.org/10.5194/gmd-15-9057-2022, https://doi.org/10.5194/gmd-15-9057-2022, 2022
Short summary
Short summary
An adaptive observation error inflation (AOEI) method was proposed for atmospheric data assimilation to mitigate erroneous analysis updates caused by large observation-minus-forecast differences for satellite brightness temperature around clear- and cloudy-sky boundaries. This study implemented the AOEI with an ocean data assimilation system, leading to an improvement of analysis accuracy and dynamical balance around the frontal regions with large meridional temperature differences.
Diego Bruciaferri, Marina Tonani, Isabella Ascione, Fahad Al Senafi, Enda O'Dea, Helene T. Hewitt, and Andrew Saulter
Geosci. Model Dev., 15, 8705–8730, https://doi.org/10.5194/gmd-15-8705-2022, https://doi.org/10.5194/gmd-15-8705-2022, 2022
Short summary
Short summary
More accurate predictions of the Gulf's ocean dynamics are needed. We investigate the impact on the predictive skills of a numerical shelf sea model of the Gulf after changing a few key aspects. Increasing the lateral and vertical resolution and optimising the vertical coordinate system to best represent the leading physical processes at stake significantly improve the accuracy of the simulated dynamics. Additional work may be needed to get real benefit from using a more realistic bathymetry.
Matthias Gröger, Manja Placke, H. E. Markus Meier, Florian Börgel, Sandra-Esther Brunnabend, Cyril Dutheil, Ulf Gräwe, Magnus Hieronymus, Thomas Neumann, Hagen Radtke, Semjon Schimanke, Jian Su, and Germo Väli
Geosci. Model Dev., 15, 8613–8638, https://doi.org/10.5194/gmd-15-8613-2022, https://doi.org/10.5194/gmd-15-8613-2022, 2022
Short summary
Short summary
Comparisons of oceanographic climate data from different models often suffer from different model setups, forcing fields, and output of variables. This paper provides a protocol to harmonize these elements to set up multidecadal simulations for the Baltic Sea, a marginal sea in Europe. First results are shown from six different model simulations from four different model platforms. Topical studies for upwelling, marine heat waves, and stratification are also assessed.
Shun Ohishi, Tsutomu Hihara, Hidenori Aiki, Joji Ishizaka, Yasumasa Miyazawa, Misako Kachi, and Takemasa Miyoshi
Geosci. Model Dev., 15, 8395–8410, https://doi.org/10.5194/gmd-15-8395-2022, https://doi.org/10.5194/gmd-15-8395-2022, 2022
Short summary
Short summary
We develop an ensemble-Kalman-filter-based regional ocean data assimilation system in which satellite and in situ observations are assimilated at a daily frequency. We find the best setting for dynamical balance and accuracy based on sensitivity experiments focused on how to inflate the ensemble spread and how to apply the analysis update to the model evolution. This study has a broader impact on more general data assimilation systems in which the initial shocks are a significant issue.
Hector S. Torres, Patrice Klein, Jinbo Wang, Alexander Wineteer, Bo Qiu, Andrew F. Thompson, Lionel Renault, Ernesto Rodriguez, Dimitris Menemenlis, Andrea Molod, Christopher N. Hill, Ehud Strobach, Hong Zhang, Mar Flexas, and Dragana Perkovic-Martin
Geosci. Model Dev., 15, 8041–8058, https://doi.org/10.5194/gmd-15-8041-2022, https://doi.org/10.5194/gmd-15-8041-2022, 2022
Short summary
Short summary
Wind work at the air-sea interface is the scalar product of winds and currents and is the transfer of kinetic energy between the ocean and the atmosphere. Using a new global coupled ocean-atmosphere simulation performed at kilometer resolution, we show that all scales of winds and currents impact the ocean dynamics at spatial and temporal scales. The consequential interplay of surface winds and currents in the numerical simulation motivates the need for a winds and currents satellite mission.
Zhanpeng Zhuang, Quanan Zheng, Yongzeng Yang, Zhenya Song, Yeli Yuan, Chaojie Zhou, Xinhua Zhao, Ting Zhang, and Jing Xie
Geosci. Model Dev., 15, 7221–7241, https://doi.org/10.5194/gmd-15-7221-2022, https://doi.org/10.5194/gmd-15-7221-2022, 2022
Short summary
Short summary
We evaluate the impacts of surface waves and internal tides on the upper-ocean mixing in the Indian Ocean. The surface-wave-generated turbulent mixing is dominant if depth is < 30 m, while the internal-tide-induced mixing is larger than surface waves in the ocean interior from 40
to 130 m. The simulated thermal structure, mixed layer depth and surface current are all improved when the mixing schemes are jointly incorporated into the ocean model because of the strengthened vertical mixing.
Giulia Bonino, Doroteaciro Iovino, Laurent Brodeau, and Simona Masina
Geosci. Model Dev., 15, 6873–6889, https://doi.org/10.5194/gmd-15-6873-2022, https://doi.org/10.5194/gmd-15-6873-2022, 2022
Short summary
Short summary
The sea surface temperature (SST) is highly influenced by the transfer of energy driven by turbulent air–sea fluxes (TASFs). In the NEMO ocean general circulation model, TASFs are computed by means of bulk formulas. Bulk formulas require the choice of a given bulk parameterization, which influences the magnitudes of the TASFs. Our results show that parameterization-related SST differences are primarily sensitive to the wind stress differences across parameterizations.
Gustavo M. Marques, Nora Loose, Elizabeth Yankovsky, Jacob M. Steinberg, Chiung-Yin Chang, Neeraja Bhamidipati, Alistair Adcroft, Baylor Fox-Kemper, Stephen M. Griffies, Robert W. Hallberg, Malte F. Jansen, Hemant Khatri, and Laure Zanna
Geosci. Model Dev., 15, 6567–6579, https://doi.org/10.5194/gmd-15-6567-2022, https://doi.org/10.5194/gmd-15-6567-2022, 2022
Short summary
Short summary
We present an idealized ocean model configuration and a set of simulations performed using varying horizontal grid spacing. While the model domain is idealized, it resembles important geometric features of the Atlantic and Southern oceans. The simulations described here serve as a framework to effectively study mesoscale eddy dynamics, to investigate the effect of mesoscale eddies on the large-scale dynamics, and to test and evaluate eddy parameterizations.
David E. Gwyther, Colette Kerry, Moninya Roughan, and Shane R. Keating
Geosci. Model Dev., 15, 6541–6565, https://doi.org/10.5194/gmd-15-6541-2022, https://doi.org/10.5194/gmd-15-6541-2022, 2022
Short summary
Short summary
The ocean current flowing along the southeastern coast of Australia is called the East Australian Current (EAC). Using computer simulations, we tested how surface and subsurface observations might improve models of the EAC. Subsurface observations are particularly important for improving simulations, and if made in the correct location and time, can have impact 600 km upstream. The stability of the current affects model estimates could be capitalized upon in future observing strategies.
Giorgio Micaletto, Ivano Barletta, Silvia Mocavero, Ivan Federico, Italo Epicoco, Giorgia Verri, Giovanni Coppini, Pasquale Schiano, Giovanni Aloisio, and Nadia Pinardi
Geosci. Model Dev., 15, 6025–6046, https://doi.org/10.5194/gmd-15-6025-2022, https://doi.org/10.5194/gmd-15-6025-2022, 2022
Short summary
Short summary
The full exploitation of supercomputing architectures requires a deep revision of the current climate models. This paper presents the parallelization of the three-dimensional hydrodynamic model SHYFEM (System of HydrodYnamic Finite Element Modules). Optimized numerical libraries were used to partition the model domain and solve the sparse linear system of equations in parallel. The performance assessment demonstrates a good level of scalability with a realistic configuration used as a benchmark.
Young-Kwang Choi, Fengyan Shi, Matt Malej, Jane M. Smith, James T. Kirby, and Stephan T. Grilli
Geosci. Model Dev., 15, 5441–5459, https://doi.org/10.5194/gmd-15-5441-2022, https://doi.org/10.5194/gmd-15-5441-2022, 2022
Short summary
Short summary
The multi-grid-nesting technique is an important methodology used for modeling transoceanic tsunamis and coastal effects. In this study, we developed a two-way nesting interface in a multi-grid-nesting system for the Boussinesq wave model FUNWAVE-TVD. The interface acts as a
backboneof the nesting framework, handling data input, output, time sequencing, and internal interactions between grids at different scales.
Chen Zhao, Rupert Gladstone, Benjamin Keith Galton-Fenzi, David Gwyther, and Tore Hattermann
Geosci. Model Dev., 15, 5421–5439, https://doi.org/10.5194/gmd-15-5421-2022, https://doi.org/10.5194/gmd-15-5421-2022, 2022
Short summary
Short summary
We use a coupled ice–ocean model to explore an oscillation feature found in several contributing models to MISOMIP1. The oscillation is closely related to the discretized grounding line retreat and likely strengthened by the buoyancy–melt feedback and/or melt–geometry feedback near the grounding line, and frequent ice–ocean coupling. Our model choices have a non-trivial impact on mean melt and ocean circulation strength, which might be interesting for the coupled ice–ocean community.
Benoît Pasquier, Sophia K. V. Hines, Hengdi Liang, Yingzhe Wu, Steven L. Goldstein, and Seth G. John
Geosci. Model Dev., 15, 4625–4656, https://doi.org/10.5194/gmd-15-4625-2022, https://doi.org/10.5194/gmd-15-4625-2022, 2022
Short summary
Short summary
Neodymium isotopes in seawater have the potential to provide key information about ocean circulation, both today and in the past. This can shed light on the underlying drivers of global climate, which will improve our ability to predict future climate change, but uncertainties in our understanding of neodymium cycling have limited use of this tracer. We present a new model of neodymium in the modern ocean that runs extremely fast, matches observations, and is freely available for development.
Pedro Duarte, Jostein Brændshøi, Dmitry Shcherbin, Pauline Barras, Jon Albretsen, Yvonne Gusdal, Nicholas Szapiro, Andreas Martinsen, Annette Samuelsen, Keguang Wang, and Jens Boldingh Debernard
Geosci. Model Dev., 15, 4373–4392, https://doi.org/10.5194/gmd-15-4373-2022, https://doi.org/10.5194/gmd-15-4373-2022, 2022
Short summary
Short summary
Sea ice models are often implemented for very large domains beyond the regions of sea ice formation, such as the whole Arctic or all of Antarctica. In this study, we implement changes in the Los Alamos Sea Ice Model, allowing it to be implemented for relatively small regions within the Arctic or Antarctica and yet considering the presence and influence of sea ice outside the represented areas. Such regional implementations are important when spatially detailed results are required.
Elias J. Hunter, Heidi L. Fuchs, John L. Wilkin, Gregory P. Gerbi, Robert J. Chant, and Jessica C. Garwood
Geosci. Model Dev., 15, 4297–4311, https://doi.org/10.5194/gmd-15-4297-2022, https://doi.org/10.5194/gmd-15-4297-2022, 2022
Short summary
Short summary
ROMSPath is an offline particle tracking model tailored for use with output from Regional Ocean Modeling System (ROMS) simulations. It is an update to an established system, the Lagrangian TRANSport (LTRANS) model, including a number of improvements. These include a modification of the model coordinate system which improved accuracy and numerical efficiency, and added functionality for nested grids and Stokes drift.
Alexander Barth, Aida Alvera-Azcárate, Charles Troupin, and Jean-Marie Beckers
Geosci. Model Dev., 15, 2183–2196, https://doi.org/10.5194/gmd-15-2183-2022, https://doi.org/10.5194/gmd-15-2183-2022, 2022
Short summary
Short summary
Earth-observing satellites provide routine measurement of several ocean parameters. However, these datasets have a significant amount of missing data due to the presence of clouds or other limitations of the employed sensors. This paper describes a method to infer the value of the missing satellite data based on a convolutional autoencoder (a specific type of neural network architecture). The technique also provides a reliable error estimate of the interpolated value.
Olivier Sulpis, Matthew P. Humphreys, Monica M. Wilhelmus, Dustin Carroll, William M. Berelson, Dimitris Menemenlis, Jack J. Middelburg, and Jess F. Adkins
Geosci. Model Dev., 15, 2105–2131, https://doi.org/10.5194/gmd-15-2105-2022, https://doi.org/10.5194/gmd-15-2105-2022, 2022
Short summary
Short summary
A quarter of the surface of the Earth is covered by marine sediments rich in calcium carbonates, and their dissolution acts as a giant antacid tablet protecting the ocean against human-made acidification caused by massive CO2 emissions. Here, we present a new model of sediment chemistry that incorporates the latest experimental findings on calcium carbonate dissolution kinetics. This model can be used to predict how marine sediments evolve through time in response to environmental perturbations.
Alisée A. Chaigneau, Guillaume Reffray, Aurore Voldoire, and Angélique Melet
Geosci. Model Dev., 15, 2035–2062, https://doi.org/10.5194/gmd-15-2035-2022, https://doi.org/10.5194/gmd-15-2035-2022, 2022
Short summary
Short summary
Climate-change-induced sea level rise is a major threat for coastal and low-lying regions. Projections of coastal sea level changes are thus of great interest for coastal risk assessment and have significantly developed in recent years. In this paper, the objective is to provide high-resolution (6 km) projections of sea level changes in the northeastern Atlantic region bordering western Europe. For that purpose, a regional model is used to refine existing coarse global projections.
Victor Onink, Erik van Sebille, and Charlotte Laufkötter
Geosci. Model Dev., 15, 1995–2012, https://doi.org/10.5194/gmd-15-1995-2022, https://doi.org/10.5194/gmd-15-1995-2022, 2022
Short summary
Short summary
Turbulent mixing is a vital process in 3D modeling of particle transport in the ocean. However, since turbulence occurs on very short spatial scales and timescales, large-scale ocean models generally have highly simplified turbulence representations. We have developed parametrizations for the vertical turbulent transport of buoyant particles that can be easily applied in a large-scale particle tracking model. The predicted vertical concentration profiles match microplastic observations well.
Gaston Irrmann, Sébastien Masson, Éric Maisonnave, David Guibert, and Erwan Raffin
Geosci. Model Dev., 15, 1567–1582, https://doi.org/10.5194/gmd-15-1567-2022, https://doi.org/10.5194/gmd-15-1567-2022, 2022
Short summary
Short summary
To be efficient on supercomputers, software must be high-performance at computing many concurrent tasks. Communications between tasks is often necessary but time consuming, and ocean modelling software NEMO 4.0 is no exception.
In this work we describe approaches enabling fewer communications, an optimization to share the workload more equally between tasks and a new flexible configuration to assess NEMO's performance easily.
Xueming Zhu, Ziqing Zu, Shihe Ren, Miaoyin Zhang, Yunfei Zhang, Hui Wang, and Ang Li
Geosci. Model Dev., 15, 995–1015, https://doi.org/10.5194/gmd-15-995-2022, https://doi.org/10.5194/gmd-15-995-2022, 2022
Short summary
Short summary
SCSOFS has provided daily updated marine forecasting in the South China Sea for the next 5 d since 2013. Comprehensive updates have been conducted to the configurations of SCSOFS's physical model and data assimilation scheme in order to improve its forecasting skill. The three most sensitive updates are highlighted. Scientific comparison and accuracy assessment results indicate that remarkable improvements have been achieved in SCSOFSv2 with respect to the original version SCSOFSv1.
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...
