We propose and analyze two novel decoupled numerical schemes for solving the Cahn-Hilliard-Stokes-Darcy (CHSD) model for two-phase flows in karstic geometry. In the first numerical scheme, we explore a fractional step method (operator splitting) to decouple the phase-field (Cahn-Hilliard equation) from the velocity field (Stokes-Darcy fluid equations). To further decouple the Stokes-Darcy system, we introduce a first order pressure stabilization term in the Darcy solver in the second numerical scheme so that the Stokes system is decoupled from the Darcy system and hence the CHSD system can be solved in a fully decoupled manner. We show that both decoupled numerical schemes are uniquely solvable, energy stable, and mass conservative. Ample numerical results are presented to demonstrate the accuracy and efficiency of our schemes.Keywords Cahn-Hilliard-Stokes-Darcy system · two phase flow · karstic geometry · interface boundary conditions · diffuse interface modelMany natural and engineering applications involve multiphase flows in karstic geometry, i.e., geometry with both conduit (or vug) and porous media [25]. Such kind of problems are intrinsically difficulty due to the multi-scale multi-physics nature. In [25], the authors utilized Onsager's extremum principle to derive a diffuse interface model, the so-called Cahn-Hilliard-Stokes-Darcy system (CHSD), for two-phase incompressible flows with matched densities in the karstic geometry. Existence and weak-strong uniqueness of weak solutions for the CHSD system have been proved recently in [28]. For complex systems like the CHSD model, efficient and accurate numerical schemes are highly desirable. There are several challenges associated with the system.First, due to the relative slow motion of fluid in porous media, long time simulations are needed in order to capture physically important phenomena. In particular, we would like to have numerical schemes that inherit, with modification if needed, the energy law of the continuous model. Second, the CHSD model, similar to all phase field model with a sharp interface limit, is stiff due to the existence of relatively steep transition regions.This stiffness leads to a severe time-step restriction if one adopts classical explicit time stepping. Third, the CHSD system involves at least three coupled physical processes: the dynamics of the phase field variable (governed by the Cahn-Hilliard equation), the fluid flow in the conduit (governed by the Stokes equation), and the fluid flow in the porous media (governed by the Darcy system). Efficient and accurate numerical schemes for each of the sub-models do exist. Therefore, it would be highly desirable to have numerical schemes that decouple these subsystems while maintaining the long time stability. Such decoupled schemes would reduce the complexity of the computation and allow the possibility of the utilization of legacy codes.In this paper, we introduce and analyze two novel decoupled schemes for the CHSD system. In particular, we show that both schemes are uniquely solvab...