Article Type:
Research Article

The Use of Spline Interpolation in Semi-Lagrangian Transport Models

L. P. Riishøjgaard Joint Center for Earth Systems Technology, University of Maryland Baltimore County, Baltimore, Maryland

Search for other papers by L. P. Riishøjgaard in
Current site
Google Scholar
PubMed Close
,
S. E. Cohn NASA/Goddard Space Flight Center, Greenbelt, Maryland

Search for other papers by S. E. Cohn in
Current site
Google Scholar
PubMed Close
,
Y. Li Joint Center for Earth Systems Technology, University of Maryland Baltimore County, Baltimore, Maryland

Search for other papers by Y. Li in
Current site
Google Scholar
PubMed Close
, and
R. Ménard Joint Center for Earth Systems Technology, University of Maryland Baltimore County, Baltimore, Maryland

Search for other papers by R. Ménard in
Current site
Google Scholar
PubMed Close
Print Publication:
01 Jul 1998
DOI:
https://doi.org/10.1175/1520-0493(1998)126<2008:TUOSII>2.0.CO;2
Page(s):
2008–2016
Restricted access

Abstract

The accuracy of interpolating semi-Lagrangian (SL) discretization methods depends on the choice of the interpolating function. Results from barotropic transport simulations on the sphere are presented, using either bicubic Lagrangian or bicubic spline SL discretization. The spline-based scheme is shown to generate excessively noisy fields in these simulations. The two methods are then tested in a one-dimensional advection problem. The damping and dispersion relations for the schemes are examined. The analysis and numerical experiments suggest that the excessive noise found in the spline-based simulations is a consequence of insufficient damping of the small scales for small and near-integer values of the Courant number. Inspection of the local Courant number for the two-dimensional spline-based simulation confirms this hypothesis. This noise can be controlled by adding a scale-selective diffusion term to the spline-based scheme, while retaining its excellent dispersion characteristics.

Corresponding author address: Dr. L. P. Riish=ojgaard, Joint Center for Earth Systems Technology, University of Maryland Baltimore County, 1000 Hilltop Circle, Baltimore, MD 21250.

Abstract

The accuracy of interpolating semi-Lagrangian (SL) discretization methods depends on the choice of the interpolating function. Results from barotropic transport simulations on the sphere are presented, using either bicubic Lagrangian or bicubic spline SL discretization. The spline-based scheme is shown to generate excessively noisy fields in these simulations. The two methods are then tested in a one-dimensional advection problem. The damping and dispersion relations for the schemes are examined. The analysis and numerical experiments suggest that the excessive noise found in the spline-based simulations is a consequence of insufficient damping of the small scales for small and near-integer values of the Courant number. Inspection of the local Courant number for the two-dimensional spline-based simulation confirms this hypothesis. This noise can be controlled by adding a scale-selective diffusion term to the spline-based scheme, while retaining its excellent dispersion characteristics.

Corresponding author address: Dr. L. P. Riish=ojgaard, Joint Center for Earth Systems Technology, University of Maryland Baltimore County, 1000 Hilltop Circle, Baltimore, MD 21250.

Save
Cover Monthly Weather Review
Monthly Weather Review

  • Bates, J. R., Y. Li, A. Brandt, S. F. McCormick, and J. Ruge, 1995:A global shallow-water numerical model based on the semi-Lagrangian advection of potential vorticity. Quart. J. Roy. Meteor. Soc.,121, 1981–2005.

  • Bermejo, R., and A. Staniforth, 1992: The conversion of semi-Lagrangian advection schemes to quasi-monotone schemes. Mon. Wea. Rev.,120, 2622–2632.

  • Davis, P. J., 1979: Circulant Matrices. Wiley, 250 pp.

  • Huang, C. Y., 1994: Semi-Lagrangian advection schemes and Eulerian WKL algorithms. Mon. Wea. Rev.,122, 1647–1658.

  • Li, Y., S. Moorthi, and J. R. Bates, 1994: Direct solution of the implicit formulation of fourth order horizontal diffusion for gridpoint models on the sphere. NASA Tech. Memo. 104606, 41 pp. [Available from NASA CASI, Linthicum Heights, MD 21090-2934.].

  • Makar, P. A., and S. R. Karpik, 1996: Basis-spline interpolation on the sphere: Applications to semi-Lagrangian advection. Mon. Wea. Rev.,124, 182–199.

  • McDonald, A., 1984: Accuracy of multiply-upstream, semi-Lagrangian advective schemes. Mon. Wea. Rev.,112, 1267–1275.

  • Moorthi, S., and R. W. Higgins, 1992: A fast direct solver for a class of two-dimensional separable elliptic equations on the sphere. NASA Tech. Memo. 104567. [Available from NASA CASI, Linthicum Heights, MD 21090-2934.].

  • Polavarapu, S., M. Tanquay, R. Ménard, and A. Staniforth, 1996: The tangent linear model for semi-Lagrangian schemes: Linearizing the process of interpolation. Tellus,48A, 74–95.

  • Press, W. H., S. A. Teukolsky, W. T. Vetterling, and B. P. Flannery, 1992: Numerical Recipes. 2d ed. Cambridge University Press, 963 pp.

  • Pudykiewicz, J., and A. Staniforth, 1984: Some properties and comparative performance of the semi-Lagrangian method of Robert in the solution of the advection–diffusion equation. Atmos.–Ocean,22, 283–308.

  • Purnell, D. K., 1976: Solution of the advective equation by upstream interpolation with a cubic spline. Mon. Wea. Rev.,104, 42–48.

  • Riishøjgaard, L. P., 1996: On four-dimensional variational assimilation of ozone data in weather prediction models. Quart. J. Roy. Meteor. Soc.,122, 1545–1571.

  • Sod, G. A., 1985: Numerical Methods in Fluid Dynamics: Initial and Boundary Value Problems. Cambridge University Press, 446 pp.

All Time Past Year Past 30 Days
Abstract Views 0 0 0
Full Text Views 176 48 5
PDF Downloads 127 36 1