“…As it is well known, the standard mixed formulations for the Stokes equations and Darcy equations earn different compatibility conditions, thus a straightforward application of the existing solvers for the Stokes equations and Darcy equations may not be feasible. To this end, a great amount of effort has been devoted to developing accurate and efficient numerical schemes for the coupled Stokes-Darcy problem, and a non-exhaustive list of these approaches include Lagrange multiplier methods [21,17,32,18], weak Galerkin method [9,22], strongly conservative methods [20,16], stabilized mixed finite element method [28,24], discontinuous Galerkin (DG) methods [26,34], virtual element method [23,33], a lowest-order staggered DG method [37] and penalty methods [38]. The coupled Stokes-Darcy problem describes multiphysics phenomena, and involves a Stokes subproblem and a Darcy subproblem, it is thus natural to resort to domain decomposition methods, which allows one to solve the coupled system sequentially with a low computational cost.…”