In the recent COVID-19 pandemic, mathematical modeling constitutes an important tool to evaluate the prospective effectiveness of non-pharmaceutical interventions (NPIs) and to guide policy-making. Most research is, however, centered around characterizing the epidemic based on point estimates like the average infectiousness or the average number of contacts. In this work, we use stochastic simulations to investigate the consequences of a population's heterogeneity regarding connectivity and individual viral load levels. Therefore, we translate a COVID-19 ODE model to a stochastic multi-agent system. We use contact networks to model complex interaction structures and a probabilistic infection rate to model individual viral load variation. We observe a large dependency of the dispersion and dynamical evolution on the population's heterogeneity that is not adequately captured by point estimates, for instance, used in ODE models. In particular, models that assume the same clinical and transmission parameters may lead to different conclusions, depending on different types of heterogeneity in the population. For instance, the existence of hubs in the contact network leads to an initial increase of dispersion and the effective reproduction number, but to a lower herd immunity threshold (HIT) compared to homogeneous populations or a population where the heterogeneity stems solely from individual infectivity variations.

Human mobility is the fuel of global pandemics. In this simulation study, we analyze how mobility restrictions mitigate epidemic processes and how this mitigation is influenced by the epidemic’s degree of dispersion.We find that (even imperfect) mobility restrictions are generally efficient in mitigating epidemic spreading. Notably, the effectiveness strongly depends on the dispersion of the offspring distribution associated with the epidemic. We also find that mobility restrictions are useful even when the pathogen is already prevalent in the whole population. However, also a delayed implementation is more efficient in the presence of overdispersion. Conclusively, this means that implementing green zones is easier for epidemics with overdispersed transmission dynamics (e.g., COVID-19). To study these relationships at an appropriate level of abstraction, we propose a spatial branching process model combining the flexibility of stochastic branching processes with an agent-based approach allowing a conceptualization of locality, saturation, and interaction structure.