Editorial Type:
Article
Article Type:
Research Article

Characteristics of a Spectral Inverse of the Laplacian Using Spherical Harmonic Functions on a Cubed-Sphere Grid for Background Error Covariance Modeling

Hyo-Jong Song Korea Institute of Atmospheric Prediction Systems, Seoul, South Korea

Search for other papers by Hyo-Jong Song in
Current site
Google Scholar
PubMed Close
,
In-Hyuk Kwon Korea Institute of Atmospheric Prediction Systems, Seoul, South Korea

Search for other papers by In-Hyuk Kwon in
Current site
Google Scholar
PubMed Close
, and
Junghan Kim Korea Institute of Atmospheric Prediction Systems, Seoul, South Korea

Search for other papers by Junghan Kim in
Current site
Google Scholar
PubMed Close
Print Publication:
01 Jan 2017
DOI:
https://doi.org/10.1175/MWR-D-16-0134.1
Page(s):
307–322
Restricted access

Abstract

In this study, a spectral inverse method using spherical harmonic functions (SHFs) represented on a cubed-sphere grid (SHF inverse) is proposed. The purpose of the spectral inverse method studied is to help with data assimilation. The grid studied is the one that results from a spectral finite element decomposition of the six faces of the cubed sphere on Gauss–Legendre–Lobatto (GLL) points with equiangular gnomonic projection. For a given discretization of the cube in this form, as the total wavenumber of the test functions increases, there comes a point at which the cube’s eigenstructure fails to be able to replicate the spherical harmonic functions. The authors call this point a limit wavenumber in using the SHF inverse. In common with the authors’ previous research, the allowable total wavenumber of the SHF inverse increases more effectively with an enhanced polynomial order. The use of the eigenvectors and eigenvalues of the Laplacian, discretized on the grid spacing used in this study, to the Poisson equation is compared with the benchmark set by using the spherical harmonics solution to the problem. In terms of accuracy, the SHF inverse is superior to a direct inverse of the Laplacian using eigendecomposition. The feasibility of SHF inverse in operational implementation is examined under a massive computational environment.

Corresponding author e-mail: Hyo-Jong Song, hj.song@kiaps.org

Abstract

In this study, a spectral inverse method using spherical harmonic functions (SHFs) represented on a cubed-sphere grid (SHF inverse) is proposed. The purpose of the spectral inverse method studied is to help with data assimilation. The grid studied is the one that results from a spectral finite element decomposition of the six faces of the cubed sphere on Gauss–Legendre–Lobatto (GLL) points with equiangular gnomonic projection. For a given discretization of the cube in this form, as the total wavenumber of the test functions increases, there comes a point at which the cube’s eigenstructure fails to be able to replicate the spherical harmonic functions. The authors call this point a limit wavenumber in using the SHF inverse. In common with the authors’ previous research, the allowable total wavenumber of the SHF inverse increases more effectively with an enhanced polynomial order. The use of the eigenvectors and eigenvalues of the Laplacian, discretized on the grid spacing used in this study, to the Poisson equation is compared with the benchmark set by using the spherical harmonics solution to the problem. In terms of accuracy, the SHF inverse is superior to a direct inverse of the Laplacian using eigendecomposition. The feasibility of SHF inverse in operational implementation is examined under a massive computational environment.

Corresponding author e-mail: Hyo-Jong Song, hj.song@kiaps.org
Save
Cover Monthly Weather Review
Monthly Weather Review

  • Amdahl, G. M., 1967: Validity of the single processor approach to achieving large-scale computing capabilities. AFIPS Conf. Proc., 30, 483–485, doi:10.1145/1465482.1465560.

  • Cheong, H.-B., and H.-G. Kang, 2015: Eigensolutions of the spherical Laplacian for the cubed-sphere and icosahedral-hexagonal grids. Quart. J. Roy. Meteor. Soc., 141, 3383–3398, doi:10.1002/qj.2620.

    • Search Google Scholar
    • Export Citation

  • Courant, R., and D. Hilbert, 1962: Methods of Mathematical Physics. Vol. 1. Wiley-Interscience, 560 pp.

  • Daley, R., 1991: Atmospheric Data Analysis. Cambridge University Press, 457 pp.

  • Fisher, M., 2003: Background error covariance modeling. Proc. Seminar on Recent Developments in Data Assimilation for Atmosphere and Ocean, Reading, United Kingdom, ECMWF, 45–63.

  • Fournier, A., M. A. Taylor, and J. J. Tribbia, 2004: The Spectral Element Atmosphere Model (SEAM): High-resolution parallel computation and localized resolution of regional dynamics. Mon. Wea. Rev., 132, 726–748, doi:10.1175/1520-0493(2004)132<0726:TSEAMS>2.0.CO;2.

    • Search Google Scholar
    • Export Citation

  • Haurwitz, B., 1940: The motion of atmospheric disturbances on the spherical earth. J. Mar. Res., 3, 254–267.

  • Holton, J. R., 2004: An Introduction to Dynamic Meteorology. Elsevier Academic Press, 535 pp.

  • Hortal, M., and A. J. Simmons, 1991: Use of reduced Gaussian grids in spectral models. Mon. Wea. Rev., 119, 1057–1074, doi:10.1175/1520-0493(1991)119<1057:UORGGI>2.0.CO;2.

    • Search Google Scholar
    • Export Citation

  • Levy, M. N., R. D. Nair, and H. M. Tufo, 2009: A high-order element-based Galerkin method for the barotropic vorticity equation. Int. J. Numer. Methods Fluids, 59, 1369–1387, doi:10.1002/fld.1874.

    • Search Google Scholar
    • Export Citation

  • Lorenc, A. C., and Coauthors, 2000: The Met. Office global three-dimensional variational data assimilation scheme. Quart. J. Roy. Meteor. Soc., 126, 2991–3012, doi:10.1002/qj.49712657002.

    • Search Google Scholar
    • Export Citation

  • Nair, R. D., S. J. Thomas, and R. D. Loft, 2005: A discontinuous Galerkin transport scheme on the cubed sphere. Mon. Wea. Rev., 133, 814–828, doi:10.1175/MWR2890.1.

    • Search Google Scholar
    • Export Citation

  • Song, H.-J., and I.-H. Kwon, 2015: Spectral transformation using a cubed-sphere grid for a three-dimensional variational data assimilation system. Mon. Wea. Rev., 143, 2581–2599, doi:10.1175/MWR-D-14-00089.1.

    • Search Google Scholar
    • Export Citation

  • Taylor, M., J. Tribbia, and M. Iskandarani, 1997: The spectral element method for the shallow water equations on the sphere. J. Comput. Phys., 130, 92–108, doi:10.1006/jcph.1996.5554.

    • Search Google Scholar
    • Export Citation

  • Williamson, D. L., J. B. Drake, J. J. Hack, R. Jakob, and P. N. Swarztrauber, 1992: A standard test set for numerical approximations to the shallow water equations in spherical geometry. J. Comput. Phys., 102, 211–224, doi:10.1016/S0021-9991(05)80016-6.

    • Search Google Scholar
    • Export Citation

  • Wu, W.-S., R. J. Purser, and D. F. Parrish, 2002: Three-dimensional variational analysis with spatially inhomogeneous covariances. Mon. Wea. Rev., 130, 2905–2916, doi:10.1175/1520-0493(2002)130<2905:TDVAWS>2.0.CO;2.

    • Search Google Scholar
    • Export Citation
All Time Past Year Past 30 Days
Abstract Views 0 0 0
Full Text Views 1078 586 46
PDF Downloads 363 50 4