A conservative semi-lagrangian discontinuous Galerkin scheme on the cubed sphere

DG methods coupled with strong-stability-preserving explicit Runge-Kutta discontinuous Galerkin (RKDG) time discretizations provide a robust numerical approach suitable for geoscience applications including atmospheric modeling. However, a major drawback of the RKDG method is its stringent Courant-Friedrichs-Lewy (CFL) stability restriction associated with explicit time stepping.

To address this issue, the authors adopt a dimension-splitting approach where a semi-Lagrangian (SL) time-stepping strategy is combined with the DG method. The resulting SLDG scheme employs a sequence of 1D operations for solving multidimensional transport equations. The SLDG scheme is inherently conservative and has the option to incorporate a local positivity-preserving filter for tracers.