Mathematical model and simulations of MERS outbreak: Predictions and implications for control measures

Abstract

The Middle East Respiratory Syndrome (MERS) has been identified in 2012 and since then outbreaks have been reported in various localities in the Middle East and in other parts of the world. To help predict the possible dynamics of MERS, as well as ways to contain it, this paper develops a mathematical model for the disease. It has a compartmental structure similar to SARS models and is in the form of a coupled system of nonlinear ordinary differential equations (ODEs). The model predictions are fitted to data from the outbreaks in Riyadh (Saudi Arabia) during 2013-2016. The results reveal that MERS will eventually be contained in the city. However, the containment time and the severity of the outbreaks depend crucially on the contact coefficients and the isolation rate constant. When randomness is added to the model coefficients, the simulations show that the model is sensitive to the scaled contact rate among people and to the isolation rate. The model is analyzed using stability theory for ODEs and indicates that when using only isolation, the endemic steady state is locally stable and attracting. Numerical simulations with parameters estimated from the city of Riyadh illustrate the analytical results and the model behavior, which may have important implications for the disease containment in the city. Indeed, the model highlights the importance of isolation of infected individuals and may be used to assess other control measures. The model is general and may be used to analyze outbreaks in other parts of the Middle East and other areas.