Topography-based local spherical Voronoi grid refinement on classical and moist shallow-water finite-volume models

Abstract

Locally refined grids for global atmospheric models are attractive since they are expected to provide an alternative to solve local phenomena without the requirement of a global high-resolution uniform grid, whose computational cost may be prohibitive. Spherical centroidal Voronoi tessellation (SCVT), as used in the atmospheric Model for Prediction Across Scales (MPAS), allows a flexible way to build and work with local refinement. In addition, the Andes Range plays a key role in the South American weather, but it is hard to capture its fine-structure dynamics in global models. This paper describes how to generate SCVT grids that are locally refined in South America and that also capture the sharp topography of the Andes Range by defining a density function based on topography and smoothing techniques. We investigate the use of the mimetic finite-volume scheme employed in the MPAS dynamical core on this grid considering the nonlinear classic and moist shallow-water equations on the sphere. We show that the local refinement, even with very smooth transitions from different resolutions, generates spurious numerical inertia–gravity waves that may even numerically destabilize the model. In the moist shallow-water model, wherein physical processes such as precipitation and cloud formation are included, our results show that the local refinement may generate spurious rain that is not observed in uniform-resolution SCVT grids. Fortunately, the spurious waves originate from small-scale grid-related numerical errors and can therefore be mitigated using fourth-order hyperdiffusion. We exploit a grid geometry-based hyperdiffusion that is able to stabilize spurious waves and has very little impact on the total energy conservation. We show that, in some cases, the clouds are better represented in a variable-resolution grid when compared to a respective uniform-resolution grid with the same number of cells, while in other cases, grid effects can affect the cloud and rain representation.