Our material point method for ice shelf flow enables error-free advection of history variables, such as damage, and ice front tracking The method can be readily implemented into existing finite element software (here, Elmer/Ice), and is suitable for large-scale application The method is verified against analytical solutions for steady-state flow and front evolution, and tested on an idealized marine ice sheet