Understanding the role of mantle convection in driving continental deformations remains an unresolved problem of geodynamics. The degree to which viscous mantle flow interact with the lithosphere can be different form one place to another depending on vigor of flow beneath a given are which can be tested using fluid-dynamic models. In this work we use Advanced Solver for Problems in Earth's ConvecTion (ASPECT), a numerical code based on finite element method, to solve for surface velocities induced by edge driven convection in the continental island, Madagascar. Our numerical model incorporates three-dimensional ellipsoid chunk geometry, an initial temperature constrained by lithospheric structure using the global model Litho.1 and elastic-viscous material behavior govern by diffusion and dislocation creep flow law. Preliminary solution shows that major down welling flow occur mostly beneath thick cratonic root at central-west of the island. Upward flow of mantle material is found at the corner of thick lithosphere at central Madagascar and thinner lithosphere at the south.