“…Eqs 4, 6 were numerically solved using a finite element (FE) code (COMSOL Multiphysics ® , version 5.5) within the framework of CFD. Many earlier studies have used this CFD code to deal with different geodynamic problems, such as mantle convection (He, 2014), magma upwelling (Shahraki and Schmeling, 2012), Rayleigh-Taylor instabilities (Ruffino et al, 2016), mid-ocean ridge development (Montési and Behn, 2007), wedge melting (Lee and Kim, 2021), and plate subduction (Carminati and Petricca, 2010;Rodriguez-Gonzalez et al, 2012). To track the evolving surface of model surface topography, we implemented an arbitrary Lagrangian Eulerian (ALE) scheme, the detailed mathematical formulation of which can be seen in earlier publications Mandal, 2022 andDasgupta et al,2021a).…”