Mathematical modeling of nonpharmaceutical interventions (NPIs) of coronavirus disease (COVID-19) in Kenya is presented. A susceptible-exposed-infected-recovered (SEIR) compartment model is considered with additional compartments of hospitalized population whose condition is severe or critical and the fatality compartment. The basic reproduction number () is computed by the next-generation matrix approach and later expressed as a time-dependent function so as to incorporate the NPIs into the model. The resulting system of ordinary differential equations (ODEs) is solved using fourth-order and fifth-order Runge–Kutta methods. Different intervention scenarios are considered, and the results show that implementation of closure of education institutions, curfew, and partial lockdown yields predicted delayed peaks of the overall infections, severe cases, and fatalities and subsequently containment of the pandemic in the country.

1. Introduction

On 7 January 2020, the World Health Organization (WHO) reported novel severe acute respiratory syndrome coronavirus (SARS-CoV-2) causative of the COVID-19 [1]. The WHO report details the chronology of how Wuhan in Hubei Province in China became a global epicenter of COVID-19 with the epidemiological link to the Huanan Seafood wet markets where there was sale of live animals. Over the period, COVID-19 has infected close to 26 million people worldwide with close to 1 million fatalities reported, necessitating the WHO to declare it a pandemic. With no viable vaccine at the moment, COVID-19 has initiated efforts focussed on containment and strategies aimed at reducing infections and fatalities. On the global platform, there has been a flurry of activity in determining the best strategy in the containment of the COVID-19 pandemic. A host of countries have instituted interventions not limited to social distancing, curfews, and total lockdowns. In [2], a model taking into consideration bats-hosts-reservoir-people transmission network was implemented, where reproduction numbers were computed to determine the transmissibility of the virus between people and reservoirs. The value of of reservoir to people was found to be 2.30 while that of people to people was 3.5. In [3], they delved into mathematical models of the COVID-19, and with a reproduction number of 2, they considered an SEIR model to model the effects of local social gatherings. With a 14-day infectious period, the model showed that, with 18-hour exposure, the attendees of the event had a protection efficacy of 70%. In [4], a conceptual SEIR model to gauge the effects of individual reactions and government reactions was investigated. This was in line with the Chinese government implementation of social distancing measures in areas considered as epicenters. The model incorporated zoonotic introductions and emigrations and showed that there were close relations with reported cases, even though asymptomatic transmission could not be established. In Ontario, Canada, Tuite et al. [5] implemented a mathematical model for transmission and mitigation strategies. They considered measures taken in the population which included quarantine, isolation of infectious cases, and hospitalization with intensive care unit (ICU) cases. It was reported that, in the absence of substantial physical social distancing and enhanced case detection and isolation, the ICU resources would be quickly overwhelmed. In [6], the authors considered a model for the prediction of new cases and made a comparison study of the United States of America, India, and Italy. Different scenarios under physical distancing measures such as lockdown were considered, and predictions were made using the available data. It was reported that, without mitigation of physical interaction of the populations, the peaks of infections would be realised much earlier in all the three countries.

Of late, the modeling efforts have advanced to cater for age-structured modeling to determine the transmission across different age strata. In [7], they use an age-structured SIR model with social contact matrices from surveys and Bayesian imputation to gauge the spread of the COVID-19 epidemic in India, where they use a generalization of the time-dependent case study data, age distribution, and social contact structure. Their predictions implore on the duration of social distancing mitigation on such lockdowns and give a scope of when periodic ease of the lockdowns can be implemented. Further age-structured morbidity and mortality are reduced on implementation of these measures. Ivorra et al. [8] investigated a more complex variant of the SEIR referred to as -SEIHRD model to take into account the undetected infections and based their study in China. Their study showed that mitigation measures of limiting physical contact decrease the basic reproduction number, and subsequently, there is a reduction in infections. Further brisk detection is needed in order to reduce the number of infections. In [9], results of social distancing strategies for curbing the COVID-19 epidemic are presented. They simulate the reduction of reproduction numbers as a result of social distancing measures. In [10], an important parameter referred to as the identification parameter is included in the simulations to predict the peaks of the epidemic in Japan. This parameter is the ratio between the positive cases to the number of tests performed. Further, they used available data and the least-square-based procedure with Poisson noise to estimate the infection rate as a function of the identification rate, and they considered an SEIR compartmental model. They report that interventions have a positive effect in delaying the peak of the epidemic, and they propose longer interventions to contain the epidemic.

The results of the research studies outlined provide techniques knowledge on implementing a mathematical model for the epidemic in Kenya. Recently, in [11], the authors outlined the forecast of the scale of the COVID-19 epidemic in Kenya. In the study, they investigated the dissimilarities between China and Kenya in terms of demographics and geography. They forecasted the potential incidence rate and magnitude of the epidemic using the observed growth rate and age distribution of confirmed cases in China. They reported that the number of infected cases is likely to be high and that isolation of the asymptomatic infectious population will not be an adequate measure. Further, they proposed exceptional social distancing to “flatten the curve” and cautioned on the expected rebound of infections when the restrictions are lifted. In [12], an SEIR model was implemented taking into consideration the population’s interaction with the environment. This was achieved by considering the environment as the host of SARS-CoV-2 virus. In the absence of social distancing, transmission and spread of COVID-19 was rife. As at 7 May 2020, Kenya had performed 28,002 tests, with 607 confirmed positive cases of COVID-19, 197 recoveries, and 29 fatalities [13]. The country’s peak is not yet predicted and, in this paper, we investigate a variant of the SEIR model to predict the peaks of the infections, severity of the illness, and fatalities with and without the nonpharmaceutical measures.

2. Model Formulation

To model the mitigation efforts to curb the transmission and spread of COVID-19, we use the compartmental susceptible-exposed-infected-recovered model. We consider human-human transmission and divide the human population at a time into eight compartments. The susceptible population is denoted by , the exposed population by , the asymptomatic infectious population by , the mild symptomatic population by , the population who are severe and hospitalized by , the critically ill population in the intensive care unit (ICU) by , the recovered population by , and the fatalities by . Hence, the total human population is given as follows:

is the human population in Kenya that may be infected with (SARS-CoV-2) virus. The susceptible individuals move to after the incubation period to the onset of the disease. The exposed individuals move to either asymptomatic or the mild symptomatic . Assuming the population in the asymptomatic have “immunocompetency,” they move to . The population from can move either to or if their conditions deteriorate warranting hospitalization, they move to and are considered severe. The severe population on assessment can move to critical compartment , and if their condition improves and becomes less critical, they are referred back to the general ward as severe cases. Further, on improved conditions, they recover and move to . The population in who succumb move to .

On 13 March 2020, after the first confirmed case, the Government of Kenya (GOK) instituted measures towards addressing the spread and transmission of COVID-19. The measures and interventions initially include a dusk-dawn curfew and massive campaign and sensitization measures of social and physical distancing. On 6 April 2020, the GOK extended and upscaled these measures with partial lockdowns in towns considered as hotspots in Kenya. These measures constituted enabled the formulation of assumptions for the development of the mathematical model. We assume the following: (i)The disease is transmitted through human-human transmissions. There are no cross-infections occurring from pathogens in the environment nor human-animal transmissions.(ii)Susceptible (S) individuals are exposed/infected through contact with infectious individuals. Each infectious individual causes an average secondary infection.(iii)After an average incubation period of 5.1 days, exposed (E) individuals either become asymptomatic (A) or exhibit mild infections (M), but not all infected persons exhibit symptoms.(iv)The virus-infected person is not infectious during the incubation period.(v)Individuals with mild infection either recover (R) or worsen to a severe case (H).(vi)Individuals with severe infection either recover (R) or worsen to a critical case (C).(vii)Critically ill individuals either return to regular hospital or die.(viii)Demographics such as birth and death rates and immigration were not considered.(ix)Only a fraction of infective individuals can be identified by diagnosis.(x)An infected individual acquires immunity upon recovery.

The flowchart below illustrates the compartmental SEIR model under consideration.

The variant of the SEIR compartmental model [14, 15] illustrated in Figure 1 culminates to an eight-dimensional dynamical system of ordinary differential equations (ODEs) given bywith the initial conditions:

The parameters used in the model are given in Table 1.

The terms and describe the rate at which the susceptible population are infected by the asymptomatic population and the mild symptomatic, respectively.

In order to gauge the extent of the spread of COVID-19, experts have recommended extensive testing in the populations, but due to economic constraints in the country and globally, this testing has been below par. Therefore, many countries have resorted to introducing a wide range of NPIs to at least slow down the epidemic spread as they gather resources for mass testing and isolation of confirmed cases. In order to simulate and predict infection cases from the model, we use the framework developed by [10]. In this case, it is assumed only a fraction of the infected individuals are identified by the diagnosis. This fraction is given by and it is referred to as the identification parameter which is . For this study, is calculated as the ratio of the confirmed cases to the tests performed. Its usage in this study is as follows: suppose , and one person is infected in Kenya of population at time . Then, the infectious individuals are given by , and hence, . At any time , the number of infectious individuals is given by . Further, at , if we assume , then .

2.1. Existence of a Disease-Free Equilibrium (DFE)

To determine the DFE of the COVID-19 in Kenya, we solve the system of equations (2) after equating the right-hand side of the system to zeros. In this case, , and subsequently, , and , which yields a DFE point as

We compute the basic reproduction number at DFE using the next-generation approach. Let and , where

Obtaining the derivatives of and at DFE point yields F and V matrices as follows:

The is the spectral radius of the product of :where is given by

This consists of two terms showing transmission by the asymptomatic population and the mild symptomatic population.

2.2. Existence of Endemic Equilibrium (EE) Point

The endemic equilibrium point of model (2) is determined as positive steady state in which the COVID-19 is said to persist in a given population.

Theorem 1. The COVID-19 model has a unique endemic equilibrium point when ; otherwise, the endemic equilibrium does not exist.

Proof. Suppose that is a nontrivial equilibrium point of (2).From (11) and (12), we getFrom (14), we getFrom (15), we getFrom (13), we getUsing (18) in (19) yieldsFrom (15), we getBoth and should satisfyUsing (23) in (22) yieldsUsing (23) in (16) yieldsUsing (24) in (21) yieldsAn endemic scenario , where is the attack rate, since almost everyone gets infected.
Therefore, solving the following equations simultaneously yields and :Hence,

2.3. Prediction of Infection Cases

In the quest to determine the level of preparedness to curb the spread of COVID-19, it is important to predict future numbers of daily infection cases and cumulative infections. This provides a basis for policymakers to know when the available healthcare infrastructure is likely to be overwhelmed. In this section, we used a Metropolis–Hastings Markov chain Monte Carlo (MCMC) algorithm for nonlinear Gaussian functions to estimate the daily and cumulative cases from Kenya’s data of daily confirmed cases. Using 20,000 iterations of the MCMC procedure, we are able to predict the daily infection cases and the cumulative infections for a period of 90 days.

Using the available daily data of the COVID-19 cases [13], we are able to use these procedures to predict the daily infection cases as shown in Figure 2 and when the country is likely to have reported the first 1,000 infection cases from Figure 3.

3. Simulation and Results

In this section, we report on the results of simulations of the COVID-19 model. The model’s equation (2) is solved using a fourth- and fifth-order Runge–Kutta method, which is implemented in MATLAB. We consider people which is Kenya’s total population as per the 2019 census. The simulations are done using the parameters listed in Table 1. The human-human contact reduction measures we consider are school closedown, dusk-to-dawn curfew, and partial lockdowns in towns and cities perceived as COVID-19 hotspots. The essence of contact reduction is to mitigate the interaction between the susceptible and infectious populations. Based on studies of COVID-19 spread in China, Italy, and Spain, we suppose that if no mitigation measures are in place, then an infectious individual would infect three secondary cases, in his/her interaction sphere. Thus, we consider the as 3, for the unmitigated scenario. The first case of the novel coronavirus was confirmed in Kenya on 13 March 2020. On 16 March 2020, schools were closed, which was the first NPI that was implemented by the Kenyan government. Then 14 days into the school closedown, curfew measures were implemented. Thereafter 24 days into school closedown, partial lockdown was implemented. In our simulations, we implement the school closure for 210 days, curfew for 196 days, and partial lockdown for 186 days. This is done via a time-dependent R0, which was taken as a cosine function (see [9]). We presume that school closure yields a reduction in contact while adding a curfew on to the school closure results in a reduction, and introduction of the partial lockdown on to the other two NPIs yields a reduction in contacts as shown in Figure 4.

A reduction in implies that fewer individuals are being infected, and hence, the disease spread is slowed down. In Figure 5, predictions of the daily confirmed cases and the prediction of the cumulative infections are given which shows that, with a basic reproduction number of greater than 1, the infections shall continue increasing.

In Figures 6(a) and 6(b), we report on the simulation of infections and the cumulative infections under different scenarios. These different scenarios are when there is no mitigation, when the schools are closed, when schools are closed and dusk-dawn curfew is instituted, and when schools are closed, curfew is in place, and partial lockdown is in place.

From the figures, we deduce that the NPIs instituted aid in containment of the epidemic spread, resulting in fewer people being infected as compared with the unmitigated scenario. In particular, we see a much higher peak of infections in the unmitigated scenario, which arrives much faster due to the high rate of infections when  = 3.0. For the first two NPIs, the infection peaks are lower and come at a later date. When all the three NPIs are combined, resulting in of 1.2, the infections peak at a much lower value and the epidemic begins to subside. However, this leaves a large population in the susceptible compartment such that when the value of begins to rise, towards the end of the NPI duration, a rebound of infections occurs. Due to the effectiveness of the combined NPIs, the new peak is much lower and the epidemic is eradicated. The measures result in a reduction in the cumulative number of infections (see Figure 6(b)). Combining all the three NPIs leads to a significant delay of the epidemic spread, a time the health system can utilize to put mechanisms and facilities to respond to the epidemic appropriately.

In Figures 7(a) and 7(b), we report on the number of severe cases of COVID-19 who are hospitalized. Reduction of the severe cases is important as this represents the proportion of the population who are directly utilizing the healthcare facilities at a given time. Due to limited hospital infrastructure and resources, delay of this peak is vital for the health systems to prepare for a guaranteed influx of patients in need of medical attention. Furthermore, any reduction of the hospitalized results in significant reduction of severe cases transfer to critical cases that require ICU services.

From Figures 7(a) and 7(b), we deduce that stringent implementation of the NPIs is key in ensuring the available hospital infrastructure is not overwhelmed during the epidemic. The national government and county governments under the Health Ministry directive have boosted the bed capacity at all hospital levels, including establishing isolation centers. If schools are closed, curfew is implemented and lockdown is instituted, the simulation shows that the hospitalized population may not overwhelm the established healthcare infrastructure capacity.

In Figures 8(a) and 8(b), we present the results of simulations of fatalities under different intervention scenarios.

Obviously, fatalities from COVID-19 will have the highest peak if NPIs are not instituted. From Figure 8(a), delayed peaks are realised if the NPIs are in place and the best scenario is if the considered measures are combined, for the considered duration of implementation. A rebound of infections also results in a resurgence of fatalities due to the severity of the disease, especially among the older population.

4. Conclusion

Since COVID-19 vaccine development will take a significantly longer period of time before it becomes globally available, it is imperative on all medical sectors in the country to adhere to mitigation measures and intervention policies set by the government. The implementation of school closure, dusk-to-dawn curfew, and partial lockdown is the measure in curbing the spread of COVID-19 that are effective and lead to the flattening of the curve within the duration of implementation. Relaxation of the mitigation measures on or before September 2020 will likely lead to a resurgence and the country may experience a new wave of infections. Devoid of the aftershocks of these interventions, the country will have no infections as at 31 December 2020. This model developed at the advent of the COVID-19 pandemic demonstrates the projections of infections and fatalities as a result of the pandemic in Kenya. It forms a basis of determination of hospital needs estimates of the country which raises the level of preparedness in case of a second wave of infections.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that they have no conflicts of interest.