FLEXPART version 11: Improved accuracy, efficiency, and flexibility

Abstract. Numerical methods and simulation codes are essential for the advancement of our understanding of complex atmospheric processes. As technology and computer hardware continue to evolve, the development of sophisticated code is vital for accurate and efficient simulations. In this paper, we present the recent advancements made in FLEXPART, a Lagrangian particle dispersion model, which has been used in a wide range of atmospheric transport studies over the past three decades, extending from tracing radionuclides from the Fukushima nuclear disaster, to inverse modelling of greenhouse gases, and to the study of atmospheric moisture cycles.

This version of FLEXPART includes notable improvements in accuracy and computational efficiency. 1) By leveraging the native vertical coordinates of European Centre for Medium Range Weather Forecasts (ECMWF) Integrated Forecasting System (IFS) instead of interpolating to terrain-following coordinates, we achieved an improvement in trajectory accuracy, leading to a ∼ 8−10 % reduction in conservation errors for semi-conserved quantities like potential vorticity. 2) The shape of aerosol particles is now accounted for in the gravitational settling and dry deposition calculation, increasing the simulation accuracy for non-spherical aerosol particles such as microplastic fibers. 3) Wet deposition has been improved by the introduction of a new below-cloud scheme, a new cloud identification scheme, and by improving the interpolation of precipitation. 4) Functionality from a separate version of FLEXPART, the FLEXPART-CTM model, is implemented, which includes linear chemical reactions. Additionally, the incorporation of Open Multi-Processing parallelisation makes the model better suited for handling large input data. Furthermore, we introduced novel methods for the input and output of particle properties and distributions. Users now have the option to run FLEXPART with more flexible particle input data, providing greater adaptability for specific research scenarios (e. g., effective backward simulations corresponding to satellite retrievals). Finally, a new user manual and restructuring of the source code into modules will serve as a basis for further development.