Dynamical Analysis of a Stochastic Ebola Model with Nonlinear Incidence Functions

The stochastic Ebola epidemic model with a combination of nonlinear incidence functions is analyzed to gain deeper insights into the complex dynamics of disease transmission. We incorporate the effects of media coverage through a saturated incidence function and consider the efficacy of protective equipment in the force of infection. Additionally, we include the impact of insufficient medication in the density-independent recovery rate to reflect its effect on treatment. Initially, we investigate the global existence and uniqueness of solutions, as well as the asymptotic behavior of solutions for both disease-free and endemic equilibrium, by constructing an appropriate Lyapunov function. We analyze the model's long-term behavior, focusing on persistence and eradication by defining two new threshold values, R-0(s) and R-0(e). When R-0(s )> 1, the presence of a stationary distribution is confirmed. Conversely, when R-0(e )< 1, the exponential extinction of disease is obtained. We derive the exact expression of the probability density function around the quasi-endemic equilibrium by solving the related Fokker-Planck equation. Finally, we perform numerical simulations to demonstrate the theoretical findings. Key results indicate that strong noise intensities can drive disease extinction, even in scenarios where deterministic models predict persistence. Additionally, the analysis highlights that media coverage, protective efficacy, safe burial practices, and recovery rate are crucial for significantly reducing the spread of infection.