A stable second-order finite volume method was proposed to predict sediment transport under rapidly varied flow conditions such as transcritical flow. For the use under unsteady flow conditions, a sediment transport model was coupled with shallow water equations. HLLC approximate Riemann solver based on a monotone upstream-centered schemes for conservation laws (MUSCL) reconstruction was used for the computation of the flux terms. From the comparisons of dam break flow experiments on erodible beds in one- and two-dimensional channels, good agreements were obtained when proper parameters were provided. Lastly, dam surface erosion problem by overtopped water was simulated. Overall, the numerical solutions showed reasonable results, which demonstrated that the proposed numerical scheme could provide stable and physical results in the cases of subcritical and supercritical flow conditions.