Background Density dependence in helminth establishment and heterogeneity in exposure to infection are known to drive resilience to interventions based on mass drug administration (MDA). However, the interaction between these processes is poorly understood. We developed a novel individual-based model for onchocerciasis transmission, EPIONCHO-IBM, which accounts for both processes. We fit the model to pre-intervention epidemiological data and explore parasite dynamics during MDA with ivermectin. Methodology/Principal findings Density dependence and heterogeneity in exposure to blackfly (vector) bites were estimated by fitting the model to matched pre-intervention microfilarial prevalence, microfilarial intensity and vector biting rate data from savannah areas of Cameroon and Cô te d'Ivoire/Burkina Faso using Latin hypercube sampling. Transmission dynamics during 25 years of annual and biannual ivermectin MDA were investigated. Density dependence in parasite establishment within humans was estimated for different levels of (fixed) exposure heterogeneity to understand how parametric uncertainty may influence treatment dynamics. Stronger overdispersion in exposure to blackfly bites results in the estimation of stronger density-dependent parasite establishment within humans, consequently increasing resilience to MDA. For