Submitted as: development and technical paper
21 Nov 2023
Submitted as: development and technical paper |  | 21 Nov 2023
Status: this preprint is currently under review for the journal GMD.

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

Jan Clemens, Lars Hoffmann, Bärbel Vogel, Sabine Grießbach, and Nicole Thomas

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 Atmosphere (CLaMS), show low computational performance 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 effectively exploits modern HPC architectures, it was previously limited to kinematic trajectories on 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.

The evaluation of the new transport scheme in MPTRAC shows that after 90-day 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 one day, internal uncertainties (e.g., due to interpolation or the numerical integration method) in the Lagrangian transport calculations are at least one 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.

Jan Clemens et al.

Status: open (until 16 Jan 2024)

Comment types: AC – author | RC – referee | CC – community | EC – editor | CEC – chief editor | : Report abuse

Jan Clemens et al.

Data sets

Implementation and evaluation of diabatic advection in the Lagrangian transport model MPTRAC 2.6 - code and data for review Creators Jan Clemens

Model code and software

Massive-Parallel Trajectory Calculations (MPTRAC) L. Hoffmann, J. Clemens, S. Griessbach, K. Haghighi Mood, F. Khosrawi, M. Liu, Y.-S. Lu, J. Sonnabend, and L. Zou

Jan Clemens et al.


Total article views: 50 (including HTML, PDF, and XML)
HTML PDF XML Total BibTeX EndNote
36 11 3 50 2 0
  • HTML: 36
  • PDF: 11
  • XML: 3
  • Total: 50
  • BibTeX: 2
  • EndNote: 0
Views and downloads (calculated since 21 Nov 2023)
Cumulative views and downloads (calculated since 21 Nov 2023)

Viewed (geographical distribution)

Total article views: 50 (including HTML, PDF, and XML) Thereof 50 with geography defined and 0 with unknown origin.
Country # Views %
  • 1
Latest update: 01 Dec 2023
Short summary
Lagrangian transport models simulate the transport of air masses in the atmosphere. For example, one model (CLaMS), is well suited for calculating transport, as it uses a special coordinate system and special vertical wind. However, it only runs inefficiently on modern supercomputers. Hence, we have implemented the benefits of CLaMS into a new model (MPTRAC), which is already highly efficient on modern supercomputers. Finally, in extensive tests, we showed that CLaMS and MPTRAC agree very well.