Abstract:Abstract. Diabatic transport schemes with hybrid zeta coordinates, which follow isentropes in the stratosphere, are known to greatly improve Lagrangian transport calculations compared to the kinematic approach. However, some Lagrangian transport calculations with a diabatic approach, such as the Chemical Lagrangian Transport Model of the Stratosphere (CLaMS), are not well prepared to run on modern high-performance computing (HPC) architectures. Here, we implemented and evaluated a new diabatic transport scheme in the Massive-Parallel Trajectory Calculations (MPTRAC) model. While MPTRAC can be used either with shared-memory multiprocessing on CPUs or with GPUs to offload computationally intensive calculations, making it flexible for many HPC applications, it has been limited to kinematic trajectories in pressure coordinates. The extended modelling approach now enables the use of either kinematic or diabatic vertical velocities and the coupling of different MPTRAC modules based on pressure or hybrid zeta coordinates. This study focus on the accuracy of the implementation in comparison to the CLaMS model. The evaluation of the new transport scheme in MPTRAC shows that, after 90 d of forward calculations, distributions of air parcels in the upper troposphere and lower stratosphere (UTLS) are almost identical for MPTRAC and CLaMS. No significant bias between the two Lagrangian models was found. Furthermore, after 1 d, internal uncertainties (e.g. due to interpolation or the numerical integration method) in the Lagrangian transport calculations are at least 1 order of magnitude smaller than external uncertainties (e.g. from reanalysis selection or downsampling of ERA5). Differences between trajectories using either CLaMS or MPTRAC are on the order of the combined internal uncertainties within MPTRAC. Since the largest systematic differences are caused by the reanalysis and the vertical velocity (diabatic vs. kinematic), the results support the development efforts for trajectory codes that can access the full resolution of ERA5 in combination with diabatic vertical velocities. This work is part of a larger effort to adapt Lagrangian transport in state-of-the-art models such as CLaMS and MPTRAC to current and future HPC architectures and exascale applications.

A Long-Time-Step-Permitting Tracer Transport Model on the Regular Latitude–Longitude Grid

Improved Local Time Step for 2D Shallow-Water Modeling Based on Unstructured Grids

Computationally Efficient Modeling of Hydro-Sediment-morphodynamic Processes Using a Hybrid Local Time Step/global Maximum Time Step

Three-Dimensional Simulation of Scalar Transport in Large Shallow Water Systems Using Flux-Form Eulerian–Lagrangian Method

Stabilized bases for high-order, interpolation semi-Lagrangian, element-based tracer transport

A Well-Balanced Numerical Model for the Simulation of Long Waves over Complex Domains

Geometry of Non-Diffusive Tracer Transport in Gridded Atmospheric Models

A Lagrangian Advection Scheme with Shape Matrix (LASM) for Solving Advection Problems

2d Well-Balanced Positivity Preserving Shallow Flow Model With Wetting And Drying On Triangular Grids

Assessing the contribution of groundwater to catchment travel time distributions through integrating conceptual flux tracking with explicit Lagrangian particle tracking

Flux-Form Eulerian-Lagrangian Method for Solving Advective Transport of Scalars in Free-Surface Flows

Two-Dimensional Finite-Volume Eulerian-Lagrangian Method on Unstructured Grid for Solving Advective Transport of Passive Scalars in Free-Surface Flows

Accelerating the Lagrangian particle tracking of residence time distributions and source water mixing towards large scales

An unstructured grid, nonhydrostatic, generalized vertical coordinate ocean model

Algorithm Development for Vertical Lagrangian Coordinate Based Non-Hydrostatic and Hydrostatic Models

Implementation and evaluation of diabatic advection in the Lagrangian transport model MPTRAC 2.6

A high order semi-Lagrangian discontinuous Galerkin method for the two-dimensional incompressible Euler equations and the guiding center Vlasov model without operator splitting

Variational Inequality Transport Model on the Sphere by the Active-Set Reduced-Space Algorithm

Cross-Scale Modeling of Shallow Water Flows in Coastal Areas with an Improved Local Time-Stepping Method

High Resolution Simulations to Transport and Diffusion over Complex Terrain

An efficient forward semi-Lagrangian model