We describe the first application of two mixed finite element methods to equations glacier evolution under different simplifying assumptions, along with a framework for implicit solution coupled velocity-thickness equations. The method uses Raviart-Thomas elements velocity and piecewise constants thickness is reframing classic staggered-grid difference case unstructured triangular meshes. secon...