A novel ALE scheme with the internal boundary for true free surface simulation in geodynamic models

Abstract. The accurate simulation of Earth's surface is essential for understanding lithospheric and mantle dynamics, especially in processes such as subduction and surface deformation. Traditional boundary conditions, such as free-slip or no-slip, do not fully capture the complex interactions occurring at the surface. The commonly used 'Sticky Air' method, while practical, suffers from several limitations, including increased computational cost and marker fluctuation issues. In this study, we propose a novel scheme within the finite element framework that integrates the 'Sticky Air' concept into an Arbitrary Lagrangian-Eulerian (ALE) formulation by employing an internal boundary to simulate a true free surface, referred to as the ALE-IB. This approach effectively addresses the limitations of existing methods, notably by reducing marker fluctuation issues and enhancing numerical stability. Moreover, it maintains a true surface in the computational domain that can be further reshaped by surface processes such as erosion and deposition, provides a foundational scheme for further coupling framework of tectonic modelling and landscape evolution modelling. We detail the theoretical formulation, implementation strategies, and validation through a series of numerical experiments. The results demonstrate that our method achieves higher accuracy and broader applicability compared to conventional techniques. Ultimately, this framework provides a more realistic and robust tool for geodynamic modelling of the Earth's free surface.