The wind industry has experienced a rapid growth in Europe over the last decades. The early generation turbines were designed for a life of 20–25 years. Decommissioning of offshore wind turbines is becoming more important since many of these installed assets are approaching their end of lifetime. In this study, using vibratory extraction of monopile foundations instead of current practice of cutting them is investigated numerically. Correct estimation of extraction force helps operators to choose suitable vibro-hammer and vessel (or crane-barge), which leads to reduction of decommissioning costs. A Coupled Eulerian-Lagrangian (CEL) approach of ABAQUS/Explicit combined with a modified Mohr-Coulomb (MMC) model are used to find the pile shaft resistance during total removal operation under saturated dense sand condition. The MMC model captures the nonlinear pre-peak hardening and post-peak softening of the dense sand which is not modelled by conventional Mohr-Coulomb model. The VUSDFLD subroutine, which is a user-defined framework, has been used to implement the MMC model into CEL analysis. A parametric study is conducted to analyze how the characteristics of the vibro-hammer, such as its frequency, eccentric moment, and the extraction rate influence the results. The present numerical results show that using proper frequency results in reduction of soil resistance to less than 25% of the initial resistance. However, appropriate hammer with enough eccentric moment and suitable extraction rate, are vital to ensure soil degradation. The results show that the proposed methodology is both robust and straightforward and it has the potential to reduce computational time which is efficient for engineering applications.