In parallel, an innovative internal erosion numerical modeling approach, based on Discrete Element-DEM and Lattice Boltzmann methods-LBM, has been developed through a partnership between the French Geological Survey (BRGM) and the Environment and Agriculture National Research Institute (IRSTEA). The upward propagation of cavities within the cover were successfully simulated. The role of different parameters (soft cover cohesion, hydraulic head, system geometry, etc.) in the sinkhole occurrence were tested by a parametric analysis.