Deriving the Error of Time Filtered Leapfrog Scheme via Modified Equations

Year 2026, Volume: 15 Issue: 1 , 480 - 491 , 24.03.2026

Ahmet Güzel

https://doi.org/10.17798/bitlisfen.1840686

https://izlik.org/JA24WC87WE

Abstract

The leapfrog (LF) scheme is a cornerstone of numerical weather prediction and large-scale atmospheric modeling due to its computational efficiency and ability to preserve the amplitude of pure oscillations during long integrations. However, the three-time-level nature of the LF method introduces a parasitic computational mode that can grow over time and contaminate physical solutions. Traditionally, the Robert-Asselin (RA) filter has been employed to suppress this mode, but it inadvertently damps the physical mode, reducing the LF scheme's formal accuracy from second to first order. This research provides a rigorous mathematical analysis of modern time filters—specifically the Robert Asselin (RA), Robert Asselin Williams (RAW), and higher-order Robert Asselin (hoRA) filters— using the method of modified equations to evaluate phase and amplitude errors. By solving the linear system for each filtered scheme, we derive equivalent linear multistep methods and their corresponding two-term modified equations. Our findings confirm that the RAW filter significantly mitigates the physical mode damping of the RA filter, recovering second-order accuracy when parameters are optimally tuned (e.g., █(α=0.53)). Furthermore, the hoRA filter demonstrates even higher performance, attaining second-order accuracy generally and third-order accuracy for the specific choice of █(β=0.4). Numerical tests on the oscillation equation validate these theoretical derivations, showing that the hoRA filter yields the lowest amplitude and phase error magnitudes compared to the RA and RAW alternatives.

Keywords

Time Filter , Modified Equations , Leapfrog Method , Error Analysis.

Ethical Statement

The study is complied with research and publication ethics.

References

• R. Bleck, “Short-range prediction in isentropic coordinates with filtered and unfiltered numerical models,” Mon. Wea. Rev., vol. 102, no. 12, pp. 813–829, 1974.
• G. J. Haltiner and J. M. McCollough, “Experiments in the initialization of a global primitive equation model,” J. Appl. Meteor., vol. 14, no. 3, pp. 281–288, 1975.
• M. Déqué and D. Cariolle, “Some destabilizing properties of the Asselin time filter,” Mon. Wea. Rev., vol. 114, no. 5, pp. 880–884, 2012/11/15 1986. [Online]. Available: https://doi.org/10.1175/1520-0493(1986)114<0880:SDPOTA>2.0.CO;2
• A. Robert and M. L´epine, “An anomaly in the behaviour of the time filter used with the leapfrog scheme in atmospheric models,” Atmosphere-Ocean, vol. 35, no. sup2, pp. S3–S15, 1997.
• R. E. Schlesinger, L. W. Uccellini, and D. R. Johnson, “The effects of the Asselin time filter on numerical solutions to the linearized shallow-water wave equations,” Mon. Wea. Rev., vol. 111, no. 3, pp. 455–467, 1983.
• E. Cordero and A. Staniforth, “A problem with the Robert–Asselin time filter for three- time-level semi-implicit semi-lagrangian discretizations,” Mon. Wea. Rev., vol. 132, no. 2, pp. 600–610, 2004.
• P. D. Williams, “A proposed modification to the Robert–Asselin time filter,” Mon. Wea. Rev., vol. 137, no. 8, pp. 2538–2546, 2009. [Online]. Available: http://dx.doi.org/10.1175/2009MWR2724.1
• N. Hurl, W. Layton, Y. Li, and C. Trenchea, “Stability analysis of the Crank–Nicolson- Leapfrog method with the Robert–Asselin–Williams time filter,” BIT, vol. 54, no. 4, pp. 1009–1021, 2014. [Online]. Available: http://dx.doi.org/10.1007/s10543-014-0493-1
• P. D. Williams, “Achieving seventh-order amplitude accuracy in leapfrog integrations,” Mon. Wea. Rev., vol. 141, pp. 3037–3051, 2013. [Online]. Available: http://dx.doi.org/10.1175/MWR-D-12-00303.1
• J. Amezcua, E. Kalnay, and P. D. Williams, “The effects of the RAW filter on the climatology and forecast skill of the SPEEDY model,” Mon. Wea. Rev., vol. 139, no. 2, pp. 608–619, 2012/11/09 2011. [Online]. Available: http://dx.doi.org/10.1175/2010MWR3530.1
• Y. Li and C. Trenchea, “Analysis of time filters used with the leapfrog scheme,” in Proceedings of the VI Conference on Computational Methods for Coupled Problems in Science and Engineering, Venice, Italy, May 2015, pp. 1261–1272.
• C.-C. Young, Y.-C. Liang, Y.-H. Tseng, and C.-H. Chow, “Characteristics of the RAW- filtered leapfrog time-stepping scheme in the ocean general circulation model,” Mon. Wea. Rev., vol. 142, no. 1, pp. 434–447, 2013.
• Y. Li and C. Trenchea, “A higher-order Robert–Asselin type time filter,” J. Comput. Phys.,vol.259,pp.23–32,2014.[Online].Available: https://doi.org/10.1016/j.jcp.2013.11.022
• D. R. Durran, “The third-order Adams-Bashforth method: An attractive alternative to leapfrog time differencing,” Mon. Wea. Rev., vol. 119, no. 3, pp. 702–720, 1991. [Online]. Available: https://doi.org/10.1175/1520-0493(1991)119<0702:TTOABM>2.0.CO;2
• A. Quarteroni, R. Sacco, and F. Saleri, Numerical mathematics, 2nd ed., ser. Texts in Applied Mathematics. Berlin: Springer-Verlag, 2007, vol. 37.
• D. R. Durran, Numerical methods for fluid dynamics, 2nd ed., ser. Texts in Applied Mathematics. New York: Springer, 2010, vol. 32, with applications to geophysics. [Online]. Available: http://dx.doi.org/10.1007/978-1-4419-6412-0
• E. Hairer, S. P. Nørsett, and G. Wanner, Solving ordinary differential equations. I, 2nd ed., ser. Springer Series in Computational Mathematics. Berlin: Springer-Verlag, 1993, vol. 8, nonstiff problems.
• Y. Li and C. Trenchea, “Analysis of time filters used with the leapfrog scheme,” in Coupled Problems in Science and Engineering VI COUPLED PROBLEMS 2015, E. O. n. B. Schrefler and M. Papadrakakis, Eds. Barcelona, Spain: International Center for Numerical Methods in Engineering (CIMNE), april 2015, pp. 1261–1272.
• D. F. Griffiths and D. J. Higham, Numerical methods for ordinary differential equations, ser. Springer Undergraduate Mathematics Series. Springer-Verlag London, Ltd., London, 2010, initial value problems. [Online]. Available: http://dx.doi.org/10.1007/978-0-85729-148-6
• R. Asselin, “Frequency filter for time integrations,” Mon. Wea. Rev., 100(6), pp. 487-490, 1972. https://doi.org/10.1175/1520-0493(1972)100<0487:FFFTI>2.3.CO;2
• A. Guzel, “Higher Order Time Filters for Evolution Equations,” D-Scholarship@Pitt, Jun. 28, 2018. https://d-scholarship.pitt.edu/concern/etds/77efa379-a048-4032-9a81-2daf72c37334