Inert Trace Constituent Transport in Sigma and Hybrid Isentropic–Sigma Models. Part II: Twelve Semi-Lagrangian Algorithms

Fred M. Reames Space Science and Engineering Center, University of Wisconsin—Madison, Madison, Wisconsin

Search for other papers by Fred M. Reames in
Current site
Google Scholar
PubMed
Close
and
Tom H. Zapotocny Space Science and Engineering Center, University of Wisconsin—Madison, Madison, Wisconsin

Search for other papers by Tom H. Zapotocny in
Current site
Google Scholar
PubMed
Close
Restricted access

We are aware of a technical issue preventing figures and tables from showing in some newly published articles in the full-text HTML view.
While we are resolving the problem, please use the online PDF version of these articles to view figures and tables.

Abstract

Part I of this paper examined nine trace constituent advection algorithms as applied in channel versions of the University of Wisconsin hybrid isentropic–sigma (θ–σ) and sigma (σ) models. This paper examines the performance of 12 semi-Lagrangian transport (SLT) algorithms in the same models. The interpolants of trace constituent include second- through twelfth-order Lagrangian, an overlapping polynomial, a quasi-monotonic scheme, and two sixth-order schemes that employ derivative estimates. Additional experiments are performed that emulate the SLT algorithm in the NCAR Community Climate Model 2. As in Part I, these three-dimensional simulations are under adiabatic conditions so that conservation of the initial trace constituent maximum on an isentropic surface, and conservation of the areas between any two constituent contours, can be used as objective measures of SLT algorithm accuracy. Further, the experiments provide comparisons not only between the different SLT formulations but also between their performance in models using θ and σ coordinates.

Similar to Part I, an important result of the experiments is that comparison of algorithms is most revealing under three-dimensional transport within a baroclinic wave in which vertical transport is important. The experiments also show 1) that the “cascade” interpolation scheme is a reasonable method of greatly reducing computation time in SLT without affecting accuracy; 2) that “shape-preserving” interpolation schemes reduce accuracy in both the θ–σ and σ models; and 3) that Lagrangian interpolants of tenth and twelfth order do not significantly improve results. Comparisons to results in Part I suggest that the conservation of second-order moments advection scheme is the most consistent of all options tested.

* Current affiliation: Department of Mechanical Engineering, University of Wisconsin—Madison, Madison, Wisconsin

Corresponding author address: Dr. Fred M. Reames, SSEC, University of Wisconsin—Madison, 1225 West Dayton Street, Madison, WI 53706-1695.

Abstract

Part I of this paper examined nine trace constituent advection algorithms as applied in channel versions of the University of Wisconsin hybrid isentropic–sigma (θ–σ) and sigma (σ) models. This paper examines the performance of 12 semi-Lagrangian transport (SLT) algorithms in the same models. The interpolants of trace constituent include second- through twelfth-order Lagrangian, an overlapping polynomial, a quasi-monotonic scheme, and two sixth-order schemes that employ derivative estimates. Additional experiments are performed that emulate the SLT algorithm in the NCAR Community Climate Model 2. As in Part I, these three-dimensional simulations are under adiabatic conditions so that conservation of the initial trace constituent maximum on an isentropic surface, and conservation of the areas between any two constituent contours, can be used as objective measures of SLT algorithm accuracy. Further, the experiments provide comparisons not only between the different SLT formulations but also between their performance in models using θ and σ coordinates.

Similar to Part I, an important result of the experiments is that comparison of algorithms is most revealing under three-dimensional transport within a baroclinic wave in which vertical transport is important. The experiments also show 1) that the “cascade” interpolation scheme is a reasonable method of greatly reducing computation time in SLT without affecting accuracy; 2) that “shape-preserving” interpolation schemes reduce accuracy in both the θ–σ and σ models; and 3) that Lagrangian interpolants of tenth and twelfth order do not significantly improve results. Comparisons to results in Part I suggest that the conservation of second-order moments advection scheme is the most consistent of all options tested.

* Current affiliation: Department of Mechanical Engineering, University of Wisconsin—Madison, Madison, Wisconsin

Corresponding author address: Dr. Fred M. Reames, SSEC, University of Wisconsin—Madison, 1225 West Dayton Street, Madison, WI 53706-1695.

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

  • Carpenter, R. L., Jr., K. K. Droegemeier, P. R. Woodward, and C. E. Hane, 1990: Application of the Piecewise Parabolic Method (PPM) to meteorological modeling. Mon. Wea. Rev.,118, 586–612.

  • Hack, J. J., B. A. Boville, B. P. Briegleb, J. T. Kiehl, P. J. Rasch, and D. L. Williamson, 1993: Description of the NCAR Community Climate Model (CCM2). NCAR Tech. Note NCAR/TN-382+STR, 108 pp. [Available from NCAR, P.O. Box 3000, Boulder, CO 80307.].

  • Holnicki, P., 1995: A shape-preserving interpolation: Applications to semi-Lagrangian advection. Mon. Wea. Rev.,123, 862–870.

  • Johnson, D. R., T. H. Zapotocny, F. M. Reames, B. J. Wolf, and R. B. Pierce, 1993: A comparison of simulated precipitation by hybrid isentropic–sigma and sigma models. Mon. Wea. Rev.,121, 2088–2114.

  • Junker, N. W., and J. E. Hoke, 1990: An examination of nested grid model precipitation forecasts in the presence of moderate-to-strong low-level southerly inflow. Wea. Forecasting,5, 333–344.

  • Kalnay, E., M. Kanamitsu, and W. E. Baker, 1990: Global numerical weather prediction at the National Meteorological Center. Bull. Amer. Meteor. Soc.,71, 1410–1428.

  • Kurihara, V., and J. L. Holloway, 1967: Numerical integration of a nine level global primitive equation model formulated by the box method. Mon. Wea. Rev.,95, 509–529.

  • Matsuno, T., 1966: Numerical integrations of the primitive equations by a simulated backward difference method. J. Meteor. Soc. Japan,44, 76–84.

  • Pierce, R. B., D. R. Johnson, F. M. Reames, T. H. Zapotocny, and B. J. Wolf, 1991: Numerical investigations with a hybrid isentropic–sigma model. Part I: Normal-mode characteristics. J. Atmos. Sci.,48, 2005–2024.

  • Prather, M. J., 1986: Numerical advection by conservation of second-order moments. J. Geophys. Res.,91, 6671–6681.

  • Purser, R. J., and L. M. Leslie, 1991: An efficient interpolation procedure for high-order three-dimensional semi-Lagrangian models. Mon. Wea. Rev.,119, 2492–2498.

  • Reames, F. M., and T. H. Zapotocny, 1999: Inert trace constituent transport in sigma and hybrid isentropic–sigma models: Part I: Nine advection algorithms. Mon. Wea. Rev.,127, 173–187.

  • Whittaker, T. M., and R. A. Petersen, 1977: Objective cross-sectional analyses incorporating thermal enhancement of observed winds. Mon. Wea. Rev.,105, 147–153.

  • Williamson, D. L., 1990: Semi-Lagrangian transport in the NMC spectral model. Tellus,42A, 412–428.

  • Zapotocny, T. H., A. J. Lenzen, D. R. Johnson, T. K. Schaack, and F. M. Reames, 1997: A comparison of inert trace constituent transport between the University of Wisconsin isentropic–sigma model and the NCAR Community Climate Model. Mon. Wea. Rev.,125, 120–142.

All Time Past Year Past 30 Days
Abstract Views 0 0 0
Full Text Views 657 576 279
PDF Downloads 30 20 1