In this study, a mathematical integrated model is developed to investigate the wave-induced sloping seabed response in the vicinity of breakwater. In the present model, the wave model is based on the Volume-Averaged/Reynolds Averaged Navier–Stokes (VARANS) equations, while Biot's consolidation equation is used to govern the soil model. The influence of turbulence fluctuations on the mean flow with respect to the complicated interaction between wave, sloping seabed and breakwater are obtained by solving the Volume-Averaged k − ϵ model. Unlike previous investigations, the phase-resolved absolute shear stress is used as the source of accumulation of residual pore pressure, which can link the oscillatory and residual mechanisms simultaneously. Based on the proposed model, parametric studies regarding the effects of wave and soil characteristics as well as bed slopes on the wave-induced soil response in the vicinity of breakwater are investigated. Numerical results indicate that wave-induced seabed instability is more likely to occur in a steep slope in the case of soil with low relative density and low permeability under large wave loadings. It is also found that, the permeability of breakwater significantly affect the potential for liquefaction, especially in the region below the breakwater.