Gravitational mass movements may erode and/or entrain a significant amount of bed material that can strongly affect the flow dynamics until the moving mass eventually deposits and comes to rest. Snow avalanches generally release on slopes covered by a metastable and thus potentially erodible snow cover that can have a wide range of strength – or cohesion – depending on the type of snow and its physical properties. As the avalanche flows, the snow cover is fully or partially entrained at the front and at the base of the flow, increasing the mass of the avalanche. Conversely, at the tail, snow may be deposited along the track, reducing the overall flowing mass. The balance between entrainment and deposition therefore determines the growing or decaying of the avalanche in terms of mass. To date, it remains unclear how cohesion influences these processes and what consequences it has for avalanche dynamics and run-out. Here, we perform simulations based on the Discrete Element Method (DEM) to analyze the influence of cohesion and slope angle on the erosion, entrainment, mixing and deposition processes. This method makes it possible to follow the dynamics of the particles within the flow very precisely, something that cannot be done in real experiments. In the model, the cohesion is represented as the combined effects of a fragmentation potential associated with the strength of the bonds, and an aggregation potential associated with the stickiness of the particles. For various combinations of input parameters and material properties, we release a heap of particles over an erodible bed and simulate the entrainment and deposition mechanisms. Our results show on the one hand that a low strength (< 3 kPa) promotes a ploughing entrainment mechanism and a large entrainment velocity, up to 3 m/s. On the other hand a high strength (> 3 kPa) favors basal abrasion as the flow front is not able to destabilize the erodible bed at once. In this case, the entrainment velocity decreases typically below 1 m/s. This has important consequences on mixing: the front in granular flows with low strength and adhesion is typically made of freshly entrained material coming from the whole depth of the bed, while remains of the released material can be found at the front of highly cohesive avalanches. Finally, the deposition process is analyzed by evaluating the relationship between the deposit thickness h_stop and slope angle θ which extends the framework of the model of h_stopθ from cohesionless to cohesive granular flows. We find that a higher bond strength of the flowing material increases the deposition height. Our work improves our understanding of the mechanics of cohesive granular flows and may contribute to improving parameterizations in depth-averaged models used to simulate geophysical mass flows such as rock, ice, snow avalanches, debris flows and landslides.