Local mesh refinement significantly in uences the performance of explicit time-stepping methods for numerical wave propagation. Local time-stepping (LTS) methods improve the efficiency by using smaller time-steps precisely where the smallest mesh elements are located, thus permitting a larger time-step in the coarser regions of the mesh without violating the stability condition. However, when the mesh contains nested patches of refinement, any local time-step will be unnecessarily small in some regions. To allow for an appropriate time-step at each level of mesh refinement, multi-level local time-stepping (MLTS) methods have been proposed. Starting from the Runge{Kutta-based LTS methods derived by Grote et al. [17], we propose explicit MLTS methods of arbitrarily high accuracy. Numerical experiments with finite difference and continuous finite element spatial discretizations illustrate the usefulness of the novel MLTS methods and show that they retain the high accuracy and stability of the underlying Runge{Kutta methods
