Articles | Volume 11, issue 2
https://doi.org/10.5194/gmd-11-575-2018
© Author(s) 2018. This work is distributed under
the Creative Commons Attribution 3.0 License.
the Creative Commons Attribution 3.0 License.
https://doi.org/10.5194/gmd-11-575-2018
© Author(s) 2018. This work is distributed under
the Creative Commons Attribution 3.0 License.
the Creative Commons Attribution 3.0 License.
Development and technical paper
| 
08 Feb 2018
Development and technical paper |
| 08 Feb 2018
Trajectory errors of different numerical integration schemes diagnosed with the MPTRAC advection module driven by ECMWF operational analyses
Thomas Rößler
CORRESPONDING AUTHOR
Forschungszentrum Jülich, Jülich Supercomputing Centre,
Jülich, Germany
Department of Mathematics,
University of Wisconsin–Milwaukee, Milwaukee, WI, USA
Olaf Stein
Forschungszentrum Jülich, Jülich Supercomputing Centre,
Jülich, Germany
Yi Heng
School of Chemical Engineering and Technology, Sun
Yat-sen University, Guangzhou, China
Paul Baumeister
Forschungszentrum Jülich, Jülich Supercomputing Centre,
Jülich, Germany
Lars Hoffmann
Forschungszentrum Jülich, Jülich Supercomputing Centre,
Jülich, Germany
Related authors
No articles found.
Phoebe Noble, Haruka Okui, Joan Alexander, Manfred Ern, Neil P. Hindley, Lars Hoffmann, Laura Holt, Annelize van Niekerk, Riwal Plougonven, Inna Polichtchouk, Claudia C. Stephan, Martina Bramberger, Milena Corcos, William Putnam, Christopher Kruse, and Corwin J. Wright
Atmos. Chem. Phys., 26, 7607–7630, https://doi.org/10.5194/acp-26-7607-2026, https://doi.org/10.5194/acp-26-7607-2026, 2026
Short summary
Short summary
Gravity waves are small-scale processes that drive the circulation in the middle and upper atmosphere. In this work, we assess 3 new high-resolution (3-5km horizontal resolution) models against satellite data. Generally, models capture the spatial patterns and represent stratospheric northern hemisphere mountain generated waves well. However, they still underestimate amplitudes globally and struggle with the representation of southern hemispheric convective waves.
Yen-Sen Lu, Lars Hoffmann, Xue Wu, Corwin J. Wright, and Neil P. Hindley
EGUsphere, https://doi.org/10.5194/egusphere-2026-1469, https://doi.org/10.5194/egusphere-2026-1469, 2026
Short summary
Short summary
We studied how computer models represent atmospheric waves triggered by intense tropical cyclones. We used a real typhoon as a case study to test how different mathematical descriptions of the atmospheric model affect these wave patterns. We found that small changes in model settings significantly impact the predicted waves in the upper atmosphere. This research helps scientists choose better model configurations to improve weather forecasts and our understanding of atmospheric science.
Zhen Yang, Bärbel Vogel, Felix Plöger, Zhixuan Bai, Dan Li, Sabine Griessbach, Lars Hoffmann, Frank G. Wienhold, Elizabeth Asher, Alexandre A. Baron, Katie R. Smith, Troy Thornberry, Jianchun Bian, and Michaela I. Hegglin
Atmos. Chem. Phys., 26, 4749–4769, https://doi.org/10.5194/acp-26-4749-2026, https://doi.org/10.5194/acp-26-4749-2026, 2026
Short summary
Short summary
Balloon measurements over Lhasa, China, combined with satellite guided modelling, tracked aerosol from the 2019 Raikoke eruption through the Asian summer monsoon anticyclone. We find two altitude dependent routes into the anticyclone. It partly blocks transport but also allows entry and mixing, explaining the gradual weakening of the layers. Balloon measurements over Boulder, United States, outside the anticyclone, provide an independent check outside the region.
Peter G. Berthelemy, Corwin J. Wright, Neil P. Hindley, Phoebe E. Noble, and Lars Hoffmann
Atmos. Chem. Phys., 25, 17595–17611, https://doi.org/10.5194/acp-25-17595-2025, https://doi.org/10.5194/acp-25-17595-2025, 2025
Short summary
Short summary
Atmospheric gravity waves are one of the key mechanisms for moving energy upwards through the atmosphere. We use temperature data to see them from a satellite, and here have made a new method to automatically detect them. This works by seeing if points next to each other are from the same wave. This is useful for creating larger gravity wave datasets without noise, which can then be used by climate forecasters to improve their understanding of the atmosphere.
Arno Keppens, Daan Hubert, José Granville, Oindrila Nath, Jean-Christopher Lambert, Catherine Wespes, Pierre-François Coheur, Cathy Clerbaux, Anne Boynard, Richard Siddans, Barry Latter, Brian Kerridge, Serena Di Pede, Pepijn Veefkind, Juan Cuesta, Gaelle Dufour, Klaus-Peter Heue, Melanie Coldewey-Egbers, Diego Loyola, Andrea Orfanoz-Cheuquelaf, Swathi Maratt Satheesan, Kai-Uwe Eichmann, Alexei Rozanov, Viktoria F. Sofieva, Jerald R. Ziemke, Antje Inness, Roeland Van Malderen, and Lars Hoffmann
Atmos. Meas. Tech., 18, 6893–6916, https://doi.org/10.5194/amt-18-6893-2025, https://doi.org/10.5194/amt-18-6893-2025, 2025
Short summary
Short summary
The first Tropospheric Ozone Assessment Report (TOAR) encountered discrepancies between several satellite sensors’ estimates of the distribution and change of ozone in the free troposphere. Therefore, contributing to the second TOAR, we harmonise as much as possible the observational perspective of sixteen tropospheric ozone products from satellites. This only partially accounts for the observed discrepancies, with a reduction of 10–40 % of the inter-product dispersion upon harmonisation.
Farahnaz Khosrawi and Lars Hoffmann
EGUsphere, https://doi.org/10.5194/egusphere-2025-3147, https://doi.org/10.5194/egusphere-2025-3147, 2025
Preprint archived
Short summary
Short summary
Computer performance has increased immensely in recent years, but the ability to store data has only increased slightly. This presents scientists with major challenges. Many compression methods have been developed in recent years with which data can be stored either lossless or lossy. Here we test three of these methods: two lossy compression methods and one lossless compressor. Our study shows that compression is a valuable tool to cope with the high demand of disk space from these data sets.
Mingzhao Liu, Lars Hoffmann, Jens-Uwe Grooß, Zhongyin Cai, Sabine Grießbach, and Yi Heng
Atmos. Chem. Phys., 25, 4403–4418, https://doi.org/10.5194/acp-25-4403-2025, https://doi.org/10.5194/acp-25-4403-2025, 2025
Short summary
Short summary
We studied the transport and chemical decomposition of volcanic SO2, focusing on the 2019 Raikoke event. By comparing two different chemistry modeling schemes, we found that including complex chemical reactions leads to a more accurate prediction of how long SO2 stays in the atmosphere. This research helps improve our understanding of volcanic pollution and its impact on air quality and climate, providing better tools for scientists to track and predict the movement of these pollutants.
Ling Zou, Reinhold Spang, Sabine Griessbach, Lars Hoffmann, Farahnaz Khosrawi, Rolf Müller, and Ines Tritscher
Atmos. Chem. Phys., 24, 11759–11774, https://doi.org/10.5194/acp-24-11759-2024, https://doi.org/10.5194/acp-24-11759-2024, 2024
Short summary
Short summary
This study provided estimates of the occurrence of ice polar stratospheric clouds (PSCs) observed by the Michelson Interferometer for Passive Atmospheric Sounding (MIPAS) and their connection with temperatures above the frost point (Tice) using a Lagrangian model derived from ERA5. We found that ice PSCs above Tice with temperature fluctuations along the backward trajectory are 33 % in the Arctic and 9 % in the Antarctic. This quantitative assessment enhances our understanding of ice PSCs.
Jan Clemens, Lars Hoffmann, Bärbel Vogel, Sabine Grießbach, and Nicole Thomas
Geosci. Model Dev., 17, 4467–4493, https://doi.org/10.5194/gmd-17-4467-2024, https://doi.org/10.5194/gmd-17-4467-2024, 2024
Short summary
Short summary
Lagrangian transport models simulate the transport of air masses in the atmosphere. For example, one model (CLaMS) is well suited to 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.
Lars Hoffmann, Kaveh Haghighi Mood, Andreas Herten, Markus Hrywniak, Jiri Kraus, Jan Clemens, and Mingzhao Liu
Geosci. Model Dev., 17, 4077–4094, https://doi.org/10.5194/gmd-17-4077-2024, https://doi.org/10.5194/gmd-17-4077-2024, 2024
Short summary
Short summary
Lagrangian particle dispersion models are key for studying atmospheric transport but can be computationally intensive. To speed up simulations, the MPTRAC model was ported to graphics processing units (GPUs). Performance optimization of data structures and memory alignment resulted in runtime improvements of up to 75 % on NVIDIA A100 GPUs for ERA5-based simulations with 100 million particles. These optimizations make the MPTRAC model well suited for future high-performance computing systems.
Jan Clemens, Bärbel Vogel, Lars Hoffmann, Sabine Griessbach, Nicole Thomas, Suvarna Fadnavis, Rolf Müller, Thomas Peter, and Felix Ploeger
Atmos. Chem. Phys., 24, 763–787, https://doi.org/10.5194/acp-24-763-2024, https://doi.org/10.5194/acp-24-763-2024, 2024
Short summary
Short summary
The source regions of the Asian tropopause aerosol layer (ATAL) are debated. We use balloon-borne measurements of the layer above Nainital (India) in August 2016 and atmospheric transport models to find ATAL source regions. Most air originated from the Tibetan plateau. However, the measured ATAL was stronger when more air originated from the Indo-Gangetic Plain and weaker when more air originated from the Pacific. Hence, the results indicate important anthropogenic contributions to the ATAL.
Abhiraj Bishnoi, Olaf Stein, Catrin I. Meyer, René Redler, Norbert Eicker, Helmuth Haak, Lars Hoffmann, Daniel Klocke, Luis Kornblueh, and Estela Suarez
Geosci. Model Dev., 17, 261–273, https://doi.org/10.5194/gmd-17-261-2024, https://doi.org/10.5194/gmd-17-261-2024, 2024
Short summary
Short summary
We enabled the weather and climate model ICON to run in a high-resolution coupled atmosphere–ocean setup on the JUWELS supercomputer, where the ocean and the model I/O runs on the CPU Cluster, while the atmosphere is running simultaneously on GPUs. Compared to a simulation performed on CPUs only, our approach reduces energy consumption by 45 % with comparable runtimes. The experiments serve as preparation for efficient computing of kilometer-scale climate models on future supercomputing systems.
Bärbel Vogel, C. Michael Volk, Johannes Wintel, Valentin Lauther, Jan Clemens, Jens-Uwe Grooß, Gebhard Günther, Lars Hoffmann, Johannes C. Laube, Rolf Müller, Felix Ploeger, and Fred Stroh
Atmos. Chem. Phys., 24, 317–343, https://doi.org/10.5194/acp-24-317-2024, https://doi.org/10.5194/acp-24-317-2024, 2024
Short summary
Short summary
Over the Indian subcontinent, polluted air is rapidly uplifted to higher altitudes during the Asian monsoon season. We present an assessment of vertical transport in this region using different wind data provided by the European Centre for Medium-Range Weather Forecasts (ECMWF), as well as high-resolution aircraft measurements. In general, our findings confirm that the newest ECMWF reanalysis product, ERA5, yields a better representation of transport compared to the predecessor, ERA-Interim.
Xue Wu, Lars Hoffmann, Corwin J. Wright, Neil P. Hindley, M. Joan Alexander, Silvio Kalisch, Xin Wang, Bing Chen, Yinan Wang, and Daren Lyu
EGUsphere, https://doi.org/10.5194/egusphere-2023-3008, https://doi.org/10.5194/egusphere-2023-3008, 2024
Preprint archived
Short summary
Short summary
This study identified a noteworthy time-lagged correlation between hurricane intensity and stratospheric gravity wave intensities during hurricane intensification. Meanwhile, the study reveals distinct frequencies, horizontal wavelengths, and vertical wavelengths in the inner core region during hurricane intensification, offering essential insights for monitoring hurricane intensity via satellite observations of stratospheric gravity waves.
Mingzhao Liu, Lars Hoffmann, Sabine Griessbach, Zhongyin Cai, Yi Heng, and Xue Wu
Geosci. Model Dev., 16, 5197–5217, https://doi.org/10.5194/gmd-16-5197-2023, https://doi.org/10.5194/gmd-16-5197-2023, 2023
Short summary
Short summary
We introduce new and revised chemistry and physics modules in the Massive-Parallel Trajectory Calculations (MPTRAC) Lagrangian transport model aiming to improve the representation of volcanic SO2 transport and depletion. We test these modules in a case study of the Ambae eruption in July 2018 in which the SO2 plume underwent wet removal and convection. The lifetime of SO2 shows highly variable and complex dependencies on the atmospheric conditions at different release heights.
Zbigniew P. Piotrowski, Jaro Hokkanen, Daniel Caviedes-Voullieme, Olaf Stein, and Stefan Kollet
EGUsphere, https://doi.org/10.5194/egusphere-2023-1079, https://doi.org/10.5194/egusphere-2023-1079, 2023
Preprint withdrawn
Short summary
Short summary
The computer programs capable of simulation of Earth system components evolve, adapting new fundamental science concepts and more observational data on more and more powerful computer hardware. Adaptation of a large scientific program to a new type of hardware is costly. In this work we propose cheap and simple but effective strategy that enable computation using graphic processing units, based on automated program code modification. This results in better resolution and/or longer predictions.
Lars Hoffmann, Paul Konopka, Jan Clemens, and Bärbel Vogel
Atmos. Chem. Phys., 23, 7589–7609, https://doi.org/10.5194/acp-23-7589-2023, https://doi.org/10.5194/acp-23-7589-2023, 2023
Short summary
Short summary
Atmospheric convection plays a key role in tracer transport in the troposphere. Global meteorological forecasts and reanalyses typically have a coarse spatiotemporal resolution that does not adequately resolve the dynamics, transport, and mixing of air associated with storm systems or deep convection. We discuss the application of the extreme convection parameterization in a Lagrangian transport model to improve simulations of tracer transport from the boundary layer into the free troposphere.
Reimar Bauer, Jens-Uwe Grooß, Jörn Ungermann, May Bär, Markus Geldenhuys, and Lars Hoffmann
Geosci. Model Dev., 15, 8983–8997, https://doi.org/10.5194/gmd-15-8983-2022, https://doi.org/10.5194/gmd-15-8983-2022, 2022
Short summary
Short summary
The Mission Support System (MSS) is an open source software package that has been used for planning flight tracks of scientific aircraft in multiple measurement campaigns during the last decade. Here, we describe the MSS software and its use during the SouthTRAC measurement campaign in 2019. As an example for how the MSS software is used in conjunction with many datasets, we describe the planning of a single flight probing orographic gravity waves propagating up into the lower mesosphere.
Paul Konopka, Mengchu Tao, Marc von Hobe, Lars Hoffmann, Corinna Kloss, Fabrizio Ravegnani, C. Michael Volk, Valentin Lauther, Andreas Zahn, Peter Hoor, and Felix Ploeger
Geosci. Model Dev., 15, 7471–7487, https://doi.org/10.5194/gmd-15-7471-2022, https://doi.org/10.5194/gmd-15-7471-2022, 2022
Short summary
Short summary
Pure trajectory-based transport models driven by meteorology derived from reanalysis products (ERA5) take into account only the resolved, advective part of transport. That means neither mixing processes nor unresolved subgrid-scale advective processes like convection are included. The Chemical Lagrangian Model of the Stratosphere (CLaMS) includes these processes. We show that isentropic mixing dominates unresolved transport. The second most important transport process is unresolved convection.
Zhongyin Cai, Sabine Griessbach, and Lars Hoffmann
Atmos. Chem. Phys., 22, 6787–6809, https://doi.org/10.5194/acp-22-6787-2022, https://doi.org/10.5194/acp-22-6787-2022, 2022
Short summary
Short summary
Using AIRS and TROPOMI sulfur dioxide retrievals and the Lagrangian transport model MPTRAC, we present an improved reconstruction of injection parameters of the 2019 Raikoke eruption. Reconstructions agree well between using AIRS nighttime and TROPOMI daytime retrievals, showing the potential of our approach to create a long-term volcanic sulfur dioxide inventory from nearly 20 years of AIRS retrievals.
Ling Zou, Sabine Griessbach, Lars Hoffmann, and Reinhold Spang
Atmos. Chem. Phys., 22, 6677–6702, https://doi.org/10.5194/acp-22-6677-2022, https://doi.org/10.5194/acp-22-6677-2022, 2022
Short summary
Short summary
Ice clouds in the stratosphere (SICs) greatly affect the water vapor balance and radiation budget in the upper troposphere and lower stratosphere (UTLS). We quantified the global SICs and analyzed their relationships with tropopause temperature, double tropopauses, UTLS clouds, gravity waves, and stratospheric aerosols. The correlations between SICs and all abovementioned processes indicate that the occurrence of and variability in SICs are spatiotemporally dependent on different processes.
Lars Hoffmann, Paul F. Baumeister, Zhongyin Cai, Jan Clemens, Sabine Griessbach, Gebhard Günther, Yi Heng, Mingzhao Liu, Kaveh Haghighi Mood, Olaf Stein, Nicole Thomas, Bärbel Vogel, Xue Wu, and Ling Zou
Geosci. Model Dev., 15, 2731–2762, https://doi.org/10.5194/gmd-15-2731-2022, https://doi.org/10.5194/gmd-15-2731-2022, 2022
Short summary
Short summary
We describe the new version (2.2) of the Lagrangian transport model MPTRAC, which has been ported for application on GPUs. The model was verified by comparing kinematic trajectories and synthetic tracer simulations for the free troposphere and stratosphere from GPUs and CPUs. Benchmarking showed a speed-up of a factor of 16 of GPU-enabled simulations compared to CPU-only runs, indicating the great potential of applying GPUs for Lagrangian transport simulations on upcoming HPC systems.
Lars Hoffmann and Reinhold Spang
Atmos. Chem. Phys., 22, 4019–4046, https://doi.org/10.5194/acp-22-4019-2022, https://doi.org/10.5194/acp-22-4019-2022, 2022
Short summary
Short summary
We present an intercomparison of 2009–2018 lapse rate tropopause characteristics as derived from ECMWF's ERA5 and ERA-Interim reanalyses. Large-scale features are similar, but ERA5 shows notably larger variability, which we mainly attribute to UTLS temperature fluctuations due to gravity waves being better resolved by ECMWF's IFS forecast model. Following evaluation with radiosondes and GPS data, we conclude ERA5 will be a more suitable asset for tropopause-related studies in future work.
Paul F. Baumeister and Lars Hoffmann
Geosci. Model Dev., 15, 1855–1874, https://doi.org/10.5194/gmd-15-1855-2022, https://doi.org/10.5194/gmd-15-1855-2022, 2022
Short summary
Short summary
The efficiency of the numerical simulation of radiative transport is shown on modern server-class graphics cards (GPUs). The low-cost prefactor on GPUs compared to general-purpose processors (CPUs) enables future large retrieval campaigns for multi-channel data from infrared sounders aboard low-orbit satellites. The validated research software JURASSIC is available in the public domain.
Cited articles
Bowman, K. P., Lin, J. C., Stohl, A., Draxler, R., Konopka, P., Andrews, A., and Brunner, D.: Input Data Requirements for Lagrangian Trajectory Models, B. Am. Meteorol. Soc., 94, 1051–1058, 2013.
Brioude, J., Angevine, W. M., McKeen, S. A., and Hsie, E.-Y.: Numerical uncertainty at mesoscale in a Lagrangian model in complex terrain, Geosci. Model Dev., 5, 1127–1136, https://doi.org/10.5194/gmd-5-1127-2012, 2012.
Buizza, R., Houtekamer, P. L., Toth, Z., Pellerin, G., Wei, M., and Zhu, Y.: A comparison of the ECMWF, MSC, and NCEP Global ensemble prediction systems, Mon. Weather Rev., 133, 1076–1097, 2005.
Butcher, J. C.: Numerical methods for ordinary differential equations, John Wiley & Sons, 2008.
CDO: Climate Data Operators, available at: http://www.mpimet.mpg.de/cdo (last access: 3 May 2017), 2015.
Davis, L. S. and Dacre, H. F.: Can dispersion model predictions be improved by increasing the temporal and spatial resolution of the meteorological input data?, Weather, 64, 232–237, 2009.
Dee, D. P., Uppala, S. M., Simmons, A. J., Berrisford, P., Poli, P., Kobayashi, S., Andrae, U., Balmaseda, M. A., Balsamo, G., Bauer, P., Bechtold, P., Beljaars, A. C. M., van de Berg, L., Bidlot, J., Bormann, N., Delsol, C., Dragani, R., Fuentes, M., Geer, A. J., Haimberger, L., Healy, S. B., Hersbach, H., Hólm, E. V., Isaksen, L., Kãllberg, P., Köhler, M., Matricardi, M., McNally, A. P., Monge-Sanz, B. M., Morcrette, J.-J., Park, B.-K., Peubey, C., de Rosnay, P., Tavolato, C., Thépaut, J.-N., and Vitart, F.: The ERA-Interim reanalysis: configuration and performance of the data assimilation system, Q. J. Roy. Meteor. Soc., 137, 553–597, 2011.
Draxler, R. R. and Hess, G. D.: An overview of the HYSPLIT_4 modeling system of trajectories, dispersion, and deposition, Aust. Meteorol. Mag., 47, 295–308, 1998.
ECMWF: CY40R1 operational analyses and forecasts, available at: http://www.ecmwf.int/en/forecasts (last access: 3 May 2017), 2013.
ECMWF: CY41R1 operational analyses and forecasts, available at: http://www.ecmwf.int/en/forecasts (last access: 3 May 2017), 2015.
Engl, H. W., Hanke, M., and Neubauer, A.: Regularization of Inverse Problems, Kluwer Academic Publishers, Dordrecht, 1996.
Harris, J. M., Draxler, R. R., and Oltmans, S. J.: Trajectory model sensitivity to differences in input data and vertical transport method, J. Geophys. Res., 110, D14109, https://doi.org/10.1029/2004JD005750, 2005.
Heng, Y., Hoffmann, L., Griessbach, S., Rößler, T., and Stein, O.: Inverse transport modeling of volcanic sulfur dioxide emissions using large-scale simulations, Geosci. Model Dev., 9, 1627–1645, https://doi.org/10.5194/gmd-9-1627-2016, 2016.
Hoffmann, L., Xue, X., and Alexander, M. J.: A global view of stratospheric gravity wave hotspots located with Atmospheric Infrared Sounder observations, J. Geophys. Res., 118, 416–434, 2013.
Hoffmann, L., Rößler, T., Griessbach, S., Heng, Y., and Stein, O.: Lagrangian transport simulations of volcanic sulfur dioxide emissions: impact of meteorological data products, J. Geophys. Res., 121, 4651–4673, https://doi.org/10.1002/2015JD023749, 2016.
Hoffmann, L., Hertzog, A., Rößler, T., Stein, O., and Wu, X.: Intercomparison of meteorological analyses and trajectories in the Antarctic lower stratosphere with Concordiasi superpressure balloon observations, Atmos. Chem. Phys., 17, 8045–8061, https://doi.org/10.5194/acp-17-8045-2017, 2017.
Hoppe, C. M., Hoffmann, L., Konopka, P., Grooß, J.-U., Ploeger, F., Günther, G., Jöckel, P., and Müller, R.: The implementation of the CLaMS Lagrangian transport core into the chemistry climate model EMAC 2.40.1: application on age of air and transport of long-lived trace species, Geosci. Model Dev., 7, 2639–2651, https://doi.org/10.5194/gmd-7-2639-2014, 2014.
Jones, A., Thomson, D., Hort, M., and Devenish, B.: The UK Met Office's next-generation atmospheric dispersion model, NAME III, in: Air Pollution Modeling and its Application XVII, Springer, 580–589, 2007.
Jülich Supercomputing Centre: JURECA: General-purpose supercomputer at Jülich Supercomputing Centre, Journal of Large-scale Research Facilities, 2, A62, https://doi.org/10.17815/jlsrf-2-121, 2016.
Kalnay, E., Kanamitsu, M., Kistler, R., Collins, W., Deaven, D., Gandin, L., Iredell, M., Saha, S., White, G., Woollen, J., Zhu, Y., Chelliah, M., Ebisuzaki, W., Higgins, W., Janowiak, J., Mo, K. C., Ropelewski, C., Wang, J., Leetmaa, A., Reynolds, R., Jenne, R., and Joseph, D.: The NCEP/NCAR 40-year reanalysis project, B. Am. Meteorol. Soc., 77, 437–471, 1996.
Kuo, Y.-H., Skumanich, M., Haagenson, P. L., and Chang, J. S.: The accuracy of trajectory models as revealed by the observing system simulation experiments, Mon. Weather Rev., 113, 1852–1867, 1985.
Legras, B., Joseph, B., and Lefèvre, F.: Vertical diffusivity in the lower stratosphere from Lagrangian back-trajectory reconstructions of ozone profiles, J. Geophys. Res., 108, 4562, https://doi.org/10.1029/2002JD003045, 2003.
Lin, J., Gerbig, C., Wofsy, S., Andrews, A., Daube, B., Davis, K., and Grainger, C.: A near-field tool for simulating the upstream influence of atmospheric observations: The Stochastic Time-Inverted Lagrangian Transport (STILT) model, J. Geophys. Res., 108, 4493, https://doi.org/10.1029/2002JD003161, 2003.
Lin, J., Brunner, D., Gerbig, C., Stohl, A., Luhar, A., and Webley, P. (Eds.): Lagrangian modeling of the atmosphere, vol. 200 of Geophysical Monograph Series, American Geophysical Union, Washington DC, 2012.
Manney, G. L., Lawrence, Z. D., Santee, M. L., Read, W. G., Livesey, N. J., Lambert, A., Froidevaux, L., Pumphrey, H. C., and Schwartz, M. J.: A minor sudden stratospheric warming with a major impact: Transport and polar processing in the 2014/2015 Arctic winter, Geophys. Res. Lett., 42, 7808–7816, https://doi.org/10.1002/2015GL065864, 2015.
Maryon, R. and Heasman, C.: The accuracy of plume trajectories forecast using the UK Meteorological Office operational forecasting models and their sensitivity to calculation schemes, Atmos. Environ., 22, 259–272, 1988.
Petterssen, S.: Weather analysis and forecasting, McGraw-Hill, New York, 1940.
Pisso, I., Real, E., Law, K. S., Legras, B., Bousserez, N., Attié, J. L., and Schlager, H.: Estimation of mixing in the troposphere from Lagrangian trace gas reconstructions during long-range pollution plume transport, J. Geophys. Res., 114, D19301, https://doi.org/10.1029/2008JD011289, 2009.
Pisso, I., Marécal, V., Legras, B., and Berthet, G.: Sensitivity of ensemble Lagrangian reconstructions to assimilated wind time step resolution, Atmos. Chem. Phys., 10, 3155–3162, https://doi.org/10.5194/acp-10-3155-2010, 2010.
Press, W. H., Teukolsky, S. A., Vetterling, W. T., and Flannery, B. P.: Numerical Recipes in C, The Art of Scientific Computing, vol. 1, Cambridge University Press, 2nd Edn., 2002.
Preusse, P., Eckermann, S. D., Ern, M., Oberheide, J., Picard, R. H., Roble, R. G., Riese, M., Russell III, J. M., and Mlynczak, M. G.: Global ray tracing simulations of the SABER gravity wave climatology, J. Geophys. Res., 114, D08126, https://doi.org/10.1029/2008JD011214, 2009.
Rabier, F., Jærvinen, H., Klinker, E., Mahfouf, J.-F., and Simmons, A.: The ECMWF operational implementation of four-dimensional variational assimilation, Q. J. Roy. Meteor. Soc., 126, 1143–1170, 2000.
Rienecker, M. M., Suarez, M. J., Gelaro, R., Todling, R., Bacmeister, J., Liu, E., Bosilovich, M. G., Schubert, S. D., Takacs, L., Kim, G.-K., Bloom, S., Chen, J., Collins, D., Conaty, A., da Silva, A., Gu, W., Joiner, J., Koster, R. D., Lucchesi, R., Molod, A., Owens, T., Pawson, S., Pegion, P., Redder, C. R., Reichle, R., Robertson, F. R., Ruddick, A. G., Sienkiewicz, M., and Woollen, J.: MERRA: NASA's Modern-Era Retrospective Analysis for Research and Applications, J. Climate, 24, 3624–3648, 2011.
Rolph, G. D. and Draxler, R. R.: Sensitivity of three-dimensional trajectories to the spatial and temporal densities of the wind field, J. Appl. Meteorol., 29, 1043–1054, 1990.
Seibert, P.: Convergence and accuracy of numerical methods for trajectory calculations, J. Appl. Meteorol., 32, 558–566, 1993.
Stohl, A.: Computation, accuracy and applications of trajectories – a review and bibliography, Atmos. Environ., 32, 947–966, 1998.
Stohl, A. and Seibert, P.: Accuracy of trajectories as determined from the conservation of meteorological tracers, Q. J. Roy. Meteor. Soc., 124, 1465–1484, 1998.
Stohl, A., Wotawa, G., Seibert, P., and Kromp-Kolb, H.: Interpolation errors in wind fields as a function of spatial and temporal resolution and their impact on different types of kinematic trajectories, J. Appl. Meteorol., 34, 2149–2165, 1995.
Stohl, A., Haimberger, L., Scheele, M., and Wernli, H.: An intercomparison of results from three trajectory models, Meteorol. Appl., 8, 127–135, 2001.
Stohl, A., Forster, C., Frank, A., Seibert, P., and Wotawa, G.: Technical note: The Lagrangian particle dispersion model FLEXPART version 6.2, Atmos. Chem. Phys., 5, 2461–2474, https://doi.org/10.5194/acp-5-2461-2005, 2005.
Stohl, A., Prata, A. J., Eckhardt, S., Clarisse, L., Durant, A., Henne, S., Kristiansen, N. I., Minikin, A., Schumann, U., Seibert, P., Stebel, K., Thomas, H. E., Thorsteinsson, T., Tørseth, K., and Weinzierl, B.: Determination of time- and height-resolved volcanic ash emissions and their use for quantitative ash dispersion modeling: the 2010 Eyjafjallajökull eruption, Atmos. Chem. Phys., 11, 4333-4351, https://doi.org/10.5194/acp-11-4333-2011, 2011.
Walmsley, J. L. and Mailhot, J.: On the numerical accuracy of trajectory models for long-range transport of atmospheric pollutants, Atmos. Ocean, 21, 14–39, 1983.
Woollings, T., Czuchnicki, C., and Franzke, C.: Twentieth century North Atlantic jet variability, Q. J. Roy. Meteor. Soc., 140, 783–791, https://doi.org/10.1002/qj.2197, 2014.
Wu, X., Griessbach, S., and Hoffmann, L.: Equatorward dispersion of a high-latitude volcanic plume and its relation to the Asian summer monsoon: a case study of the Sarychev eruption in 2009, Atmos. Chem. Phys., 17, 13439–13455, https://doi.org/10.5194/acp-17-13439-2017, 2017.
Short summary
In this study, we performed an assessment of truncation errors and computational efficiency of trajectory calculations using six popular numerical integration schemes of the Runge–Kutta family. More than 5000 transport simulations for different seasons and regions of the free troposphere and stratosphere were conducted, driven by the latest version of ECMWF operational analyses and forecasts. The study provides guidelines to achieve the most accurate and efficient trajectory calculations.
In this study, we performed an assessment of truncation errors and computational efficiency of...
