The discontinuous Galerkin (DG) method is effective and promising tool for the hyperbolic equations with discontinuities. In recent years, the DG method has been applied and shown to be effective in the analysis of shallow water equations (SWEs), which are widely used for river flooding, storm surge, tsunami, and other natural disaster simulations. In this study, the DG method is applied to the governing equations for riverbed variation, which are coupling of the SWEs and equilibrium sediment transport equation, and wetting-drying treatment and numerical flux are discussed for our purpose. The proposed method has been verified with numerical analysis such as erodible dam-break flow and sandbar migration simulation. The numerical results have shown that the proposed method is effective for riverbed variation analysis. © 2024 by the Japan Society for Computational Engineering and Science.