A Potential Enstrophy and Energy Conserving Numerical Scheme for Solution of the Shallow-Water Equations on a Geodesic Grid

Abstract Using the shallow water equations, a numerical framework on a spherical geodesic grid that conserves domain-integrated mass, potential vorticity, potential enstrophy, and total energy is developed. The numerical scheme is equally applicable to hexagonal grids on a plane and to spherical geodesic grids. This new numerical scheme is compared to its predecessor and it is shown that the new scheme does considerably better in conserving potential enstrophy and energy. Furthermore, in a simulation of geostrophic turbulence, the new numerical scheme produces energy and enstrophy spectra with slopes of approximately K−3 and K−1, respectively, where K is the total wavenumber. These slopes are in agreement with theoretical predictions. This work also exhibits a discrete momentum equation that is compatible with the Z-grid vorticity-divergence equation.

[1]  F. Winninghoff.  O the Adjustment Toward a Geostrophic Balance in a Simple Primitive Equation Model with Application to the Problems of Initialization and Objective Analysis. , 1968 .

[2]  Akio Arakawa,et al.  Integration of the Nondivergent Barotropic Vorticity Equation with AN Icosahedral-Hexagonal Grid for the SPHERE1 , 1968 .

[3]  David L. Williamson,et al.  Integration of the barotropic vorticity equation on a spherical geodesic grid , 1968 .

[4]  P. Morel,et al.  A FINITE-DIFFERENCE APPROXIMATION OF THE PRIMITIVE EQUATIONS FOR A HEXAGONAL GRID ON A PLANE , 1969 .

[5]  David L. Williamson,et al.  NUMERICAL INTEGRATION OF FLUID FLOW OVER TRIANGULAR GRIDS , 1969 .

[6]  Akio Arakawa,et al.  Computational Design of the Basic Dynamical Processes of the UCLA General Circulation Model , 1977 .

[7]  A. Arakawa,et al.  A Potential Enstrophy and Energy Conserving Scheme for the Shallow Water Equations , 1981 .

[8]  J. Holton.  Geophysical fluid dynamics. , 1983, Science.

[9]  P. Frederickson,et al.  Icosahedral Discretization of the Two-Sphere , 1985 .

[10]  Yoshinobu Masuda,et al.  An Integration Scheme of the Primitive Equation Model with an Icosahedral-Hexagonal Grid System and its Application to the Shallow Water Equations , 1986 .

[11]  David A. Randall,et al.  Geostrophic Adjustment and the Finite-Difference Shallow-Water Equations , 1994 .

[12]  R. Heikes,et al.  Numerical Integration of the Shallow-Water Equations on a Twisted Icosahedral Grid , 1995 .

[13]  David A. Randall,et al.  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 , 1995 .

[14]  G. R. Stuhne,et al.  Vortex Erosion and Amalgamation in a New Model of Large Scale Flow on the Sphere , 1996 .

[15]  J. Thuburn.  A PV-Based Shallow-Water Model on a Hexagonal-Icosahedral Grid , 1997 .

[16]  R. Salmon,et al.  Geophysical Fluid Dynamics , 2019, Classical Mechanics in Geophysical Fluid Dynamics.

[17]  G. R. Stuhne,et al.  New icosahedral grid-point discretization of the shallow water equation on the sphere , 1999 .

[18]  Todd D. Ringler,et al.  Modeling the Atmospheric General Circulation Using a Spherical Geodesic Grid: A New Class of Dynamical Cores , 2000 .

[19]  David A. Randall,et al.  Global Atmospheric Modeling Using a Geodesic Grid with an Isentropic Vertical Coordinate , 2000 .

[20]  Francis X. Giraldo,et al.  Lagrange—Galerkin methods on spherical geodesic grids: the shallow water equations , 2000 .

[21]  Todd D. Ringler,et al.  The ZM Grid: An Alternative to the Z Grid , 2002 .