Editorial Type:
Article
Article Type:
Research Article

Modeling the Atmospheric General Circulation Using a Spherical Geodesic Grid: A New Class of Dynamical Cores

Todd D. Ringler Department of Atmospheric Science, Colorado State University, Fort Collins, Colorado

Search for other papers by Todd D. Ringler in
Current site
Google Scholar
PubMed Close
,
Ross P. Heikes Department of Atmospheric Science, Colorado State University, Fort Collins, Colorado

Search for other papers by Ross P. Heikes in
Current site
Google Scholar
PubMed Close
, and
David A. Randall Department of Atmospheric Science, Colorado State University, Fort Collins, Colorado

Search for other papers by David A. Randall in
Current site
Google Scholar
PubMed Close
Print Publication:
01 Jul 2000
DOI:
https://doi.org/10.1175/1520-0493(2000)128<2471:MTAGCU>2.0.CO;2
Page(s):
2471–2490
Restricted access

Abstract

This paper documents the development and testing of a new type of atmospheric dynamical core. The model solves the vorticity and divergence equations in place of the momentum equation. The model is discretized in the horizontal using a geodesic grid that is nearly uniform over the entire globe. The geodesic grid is formed by recursively bisecting the triangular faces of a regular icosahedron and projecting those new vertices onto the surface of the sphere. All of the analytic horizontal operators are reduced to line integrals, which are numerically evaluated with second-order accuracy. In the vertical direction the model can use a variety of coordinate systems, including a generalized sigma coordinate that is attached to the top of the boundary layer. Terms related to gravity wave propagation are isolated and an efficient semi-implicit time-stepping scheme is implemented. Since this model combines many of the positive attributes of both spectral models and conventional finite-difference models into a single dynamical core, it represents a distinctively new approach to modeling the atmosphere’s general circulation.

The model is tested using the idealized forcing proposed by Held and Suarez. Results are presented for simulations using 2562 polygons (approximately 4.5° × 4.5°) and using 10 242 polygons (approximately 2.25° × 2.25°). The results are compared to those obtained with spectral model simulations truncated at T30 and T63. In terms of first and second moments of state variables such as the zonal wind, meridional wind, and temperature, the geodesic grid model results using 2562 polygons are comparable to those of a spectral model truncated at slightly less than T30, while a simulation with 10 242 polygons is comparable to a spectral model simulation truncated at slightly less than T63.

In order to further demonstrate the viability of this modeling approach, preliminary results obtained from a full-physics general circulation model that uses this dynamical core are presented. The dominant features of the DJF climate are captured in the full-physics simulation.

In terms of computational efficiency, the geodesic grid model is somewhat slower than the spectral model used for comparison. Model timings completed on an SGI Origin 2000 indicate that the geodesic grid model with 10 242 polygons is 20% slower than the spectral model truncated at T63. The geodesic grid model is more competitive at higher resolution than at lower resolution, so further optimization and future trends toward higher resolution should benefit the geodesic grid model.

Corresponding author address: Dr. Todd D. Ringler, Department of Atmospheric Science, Colorado State University, Fort Collins, CO 80523-1371.

Abstract

This paper documents the development and testing of a new type of atmospheric dynamical core. The model solves the vorticity and divergence equations in place of the momentum equation. The model is discretized in the horizontal using a geodesic grid that is nearly uniform over the entire globe. The geodesic grid is formed by recursively bisecting the triangular faces of a regular icosahedron and projecting those new vertices onto the surface of the sphere. All of the analytic horizontal operators are reduced to line integrals, which are numerically evaluated with second-order accuracy. In the vertical direction the model can use a variety of coordinate systems, including a generalized sigma coordinate that is attached to the top of the boundary layer. Terms related to gravity wave propagation are isolated and an efficient semi-implicit time-stepping scheme is implemented. Since this model combines many of the positive attributes of both spectral models and conventional finite-difference models into a single dynamical core, it represents a distinctively new approach to modeling the atmosphere’s general circulation.

The model is tested using the idealized forcing proposed by Held and Suarez. Results are presented for simulations using 2562 polygons (approximately 4.5° × 4.5°) and using 10 242 polygons (approximately 2.25° × 2.25°). The results are compared to those obtained with spectral model simulations truncated at T30 and T63. In terms of first and second moments of state variables such as the zonal wind, meridional wind, and temperature, the geodesic grid model results using 2562 polygons are comparable to those of a spectral model truncated at slightly less than T30, while a simulation with 10 242 polygons is comparable to a spectral model simulation truncated at slightly less than T63.

In order to further demonstrate the viability of this modeling approach, preliminary results obtained from a full-physics general circulation model that uses this dynamical core are presented. The dominant features of the DJF climate are captured in the full-physics simulation.

In terms of computational efficiency, the geodesic grid model is somewhat slower than the spectral model used for comparison. Model timings completed on an SGI Origin 2000 indicate that the geodesic grid model with 10 242 polygons is 20% slower than the spectral model truncated at T63. The geodesic grid model is more competitive at higher resolution than at lower resolution, so further optimization and future trends toward higher resolution should benefit the geodesic grid model.

Corresponding author address: Dr. Todd D. Ringler, Department of Atmospheric Science, Colorado State University, Fort Collins, CO 80523-1371.

Save
Cover Monthly Weather Review
Monthly Weather Review

  • Arakawa, A., and V. R. Lamb, 1977: Computational design of the basic dynamical process of the UCLA general circulation model. Methods in Computational Physics, Vol. 17, Academic Press, 173–265.

  • Baumgardner, J. R., and P. O. Frederickson, 1985: Icosahedral discretization of the two-sphere. SIAM J. Numer. Anal.,22, 1107–1115.

  • Bourke, W., 1974: A multilevel spectral model. I. Formulation and hemispheric integrations. Mon. Wea. Rev.,102, 687–701.

  • Charney, J. G., R. Fjortoft, and J. von Neumann, 1950: Numerical integration of the barotropic vorticity equation. Tellus,2, 237–254.

  • Cooley, J. W., and J. W. Tukey, 1965: An algorithm for the machine computation of complex Fourier series. Math. Comput.,19, 297–301.

  • Ding, P., and D. A. Randall, 1998: A cumulus parameterization with multiple could base levels. J. Geophys. Res.,103, 11 341–11 354.

  • Durran, D. R., 1991: The third-order Adams–Bashforth method: An attractive alternative to leapfrog time differencing. Mon. Wea. Rev.,119, 702–720.

  • Eliasen, E., B. Machenhauer, and E. Rasmussen, 1970: On a numerical method for integration of the hydrodynamical equations with a spectral representation of the horizontal fields. Institut for Teoretisk Meteorologi Rep. 2, University of Copenhagen, 37 pp.

  • Fowler, L. D., D. A. Randall, and S. A. Rutledge, 1996: Liquid and ice cloud microphysics in the CSU general circulation model. Part I: Model description and simulated microphysical processes. J. Climate,9, 489–529.

  • Greenberg, M. D., 1978: Foundations of Applied Mathematics. Prentice-Hall, 636 pp.

  • Haltiner, G. J., and R. T. Williams, 1980: Numerical Prediction and Dynamic Meteorology. John Wiley and Sons, 477 pp.

  • Harshvardhan, D. A. Randall, T. G. Corsetti, and D. A. Dazlich, 1989:Earth radiation budget and cloudiness simulations with a general circulation model. J. Atmos. Sci.,46, 1922–1942.

  • Heikes, R., 1993: The shallow water equations on a spherical geodesic grid. Department of Atmospheric Science Tech. Rep. 524, Colorado State University, Fort Collins, CO, 175 pp.

  • ——, and D. A. Randall, 1995a: Numerical integration of the shallow-water equations on a twisted icosahedral grid. Part I: Basic design and results of tests. Mon. Wea. Rev.,123, 1862–1880.

  • ——, and ——, 1995b: Numerical integration of the shallow-water equations on a twisted icosahedral grid. Part II: A detailed description of the grid and an analysis of numerical accuracy. Mon. Wea. Rev.,123, 1881–1887.

  • Held, I. M., and M. J. Suarez, 1994: A proposal for the intercomparison of the dynamical cores of atmospheric general circulation models. Bull. Amer. Meteor. Soc.,73, 1825–1830.

  • Hoskins, B. J., and A. J. Simmons, 1975: A multi-layer spectral model and the semi-implicit method. Quart. J. Roy. Meteor. Soc.,101, 636–655.

  • Hsu, Y.-J. G., and A. Arakawa, 1990: Numerical modeling of the atmosphere with an isentropic vertical coordinate. Mon. Wea. Rev.,118, 1933–1959.

  • Jaluria, Y., and K. E. Torrance, 1986: Computational Heat Transfer. Hemisphere Publishing, 366 pp.

  • Jarrud, M., and A. J. Simmon, 1983: The spectral technique. Proc. Seminar on Numerical Methods for Weather Prediction, Reading, United Kingdom, European Centre for Medium-Range Weather Forecasts, 1–60.

  • Jones, P. W., 1999: First- and second-order conservative remapping schemes for grids in spherical coordinates. Mon. Wea. Rev.,127, 2204–2210.

  • Kalnay, E., and Coauthors, 1996: The NCEP/NCAR 40-Year Reanalysis Project. Bull. Amer. Meteor. Soc.,77, 437–471.

  • Konor, C. S., and A. Arakawa, 1997: Design of an atmospheric model based on a generalized vertical coordinate. Mon. Wea. Rev.,125, 1649–1673.

  • Masuda, Y., and H. Ohnishi, 1986: An integration scheme of the primitive equations model with an icosahedral-hexagonal grid system and its application to the shallow water equations. Short- and Medium-Range Numerical Weather Prediction, T. Matsuno, Ed., Japan Meteorological Society, 317–326.

  • Orlanski, I., 1985: Advances in the theory of atmospheric fronts. Advances in Geophysics, Vol. 28, Academic Press, 223–252.

  • Orszag, S. A., 1970: Transform method for the calculation of vector-coupled sums: Application to the spectral form of the vorticity equation. J. Atmos. Sci.,27, 890–895.

  • Palmer, T. N., G. J. Shutts, and R. Swinbank, 1986: Alleviation of a systematic westerly bias in general circulation and numerical weather prediction models through an orographic gravity wave drag parameterization. Quart. J. Roy. Meteor. Soc.,112, 1001–1039.

  • Panton, R. L., 1984: Incompressible Flow. John Wiley and Sons, 780 pp.

  • Phillips, N. A., 1957: A coordinate system having some special advantages for numerical forecasting. J. Meteor.,14, 184–185.

  • Phillips, T. J., 1996: Documentation of the AMIP models on the World Wide Web. Bull. Amer. Meteor. Soc.,77, 1191–1196.

  • Randall, D. A., 1994: Geostrophic adjustment and the finite-difference shallow water equations. Mon. Wea. Rev.,122, 1371–1377.

  • Richardson, L. F., 1922: Weather Prediction by Numerical Process. Cambridge University Press, 236 pp.

  • Ringler, T. D., R. P. Heikes, and D. A. Randall, 1998: Solving the primitive equations on a spherical geodesic grid: A technical report to a new class of dynamical cores. Department of Atmospheric Science Tech. Rep. 665, Colorado State University, Fort Collins, CO, 62 pp.

  • Royer, J. F., 1986: Correction of negative mixing ratios in spectral models by global horizontal borrowing. Mon. Wea. Rev.,114, 1406–1410.

  • Sadourny, R., A. Arakawa, and Y. Mintz, 1968: Integration of the nondivergent barotropic vorticity equation with an icosahedal-hexagonal grid for the sphere. Mon. Wea. Rev.,96, 351–356.

  • Sellers, P. J., and Coauthors, 1996: A revised land surface parameterization (SiB2) for atmospheric GCMs. Part I: Model formulation. J. Climate,9, 676–705.

  • Stuhne, G. R., and W. R. Peltier, 1996: Vortex erosion and amalgamation in a new model of large scale flow on the sphere. J. Comput. Phys.,128, 58–81.

  • ——, and ——, 1999: New icosahedral grid-point discretizations of the shallow water equations on the sphere. J. Comput. Phys.,144, 23–58.

  • Suarez, M. J., A. Arakawa, and R. A. Randall, 1983: Parameterization of the planetary boundary layer in the UCLA general circulation model: Formulation and results. Mon. Wea. Rev.,111, 2224–2243.

  • Swarztrauber, P. N., D. L. Williamson, and J. B. Drake, 1997: The Cartesian method for solving partial differential equations in spherical geometry. Dyn. Atmos. Oceans,27, 679–706.

  • Tennekes, H., and J. L. Lumley, 1972: A First Course in Turbulence. The MIT Press, 300 pp.

  • Thuburn, J., 1997: A PV-based shallow-water model on a hexagonal–icosahedral grid. Mon. Wea. Rev.,125, 2328–2347.

  • Williamson, D. L., 1968: Integration of the barotropic vorticity equation on a spherical geodesic grid. Tellus,20, 642–653.

  • ——, 1976: Linear stability of finite-difference approximations on a uniform latitude–longitude grid with Fourier filtering. Mon. Wea. Rev.,104, 31–41.

  • Xie, P. P., and P. A. Arkin, 1996: Analyses of global monthly precipitation using gauge observations, satellite estimates, and numerical model predictions. J. Climate,9, 840–858.

  • Zalesak, S., 1979: Fully multidimensional flux-corrected transport algorithms for fluids. J. Comput. Phys.,31, 335–362.

All Time Past Year Past 30 Days
Abstract Views 0 0 0
Full Text Views 1181 661 81
PDF Downloads 493 110 3