A Stochastic Model with Jumps for the COVID-19 Epidemic in the Greater Abidjan Region during Public Health Measures

Our goal in this paper is to adapt to the context of the Greater Abidjan Region (Côte d'Ivoire) an existing model reflecting the evolution of the COVID-19 epidemic in Wuhan (China). This model is a deterministic compartmental model which is translated by a system of ordinary differential equations. We are reinvesting this model to obtain the parameters of the epidemic in the great Abidjan. We study some mathematical aspects of this deterministic model. Using data from the evolution of the COVID-19 epidemic in Abidjan, we obtain constants corresponding to the Greater Abidjan Region epidemic. We are performing numerical simulations using these constants to predict the behavior of the epidemic in the Greater Abidjan Region. The initial model does not take into account certain disturbances and sudden shocks which could disturb its behavior. In the Greater Abidjan Region, some continuous or sudden events disrupt the behavior of the epidemic. To stick to the realities of the Greater Abidjan Region, we introduce a white noise and jumps that correspond to the different disturbances that can occur. We obtain a stochastic model of COIVD-19 with jumps. We prove that our new model with jumps has a positive global solution. In certain conditions, this solution oscillates in the set of the disease-free equilibrium of the deterministic model. We perform numerical simulations which corroborate our theoretical results.

od of public measures application. For this reasons, we have worked on data related only to that period. For the next step, we plan to make predictions to account the evolution of this epidemic.
In [1], Pierre Magal, et al. develop a mathematic model for COVID-19 epidemic in Wuhan (China) in which the population was composed of four compartments: • S(t) is the number of susceptible individuals at time t i.e., people who are not infected yet but might become infectious individuals in the future, • I(t) is the number of asymptomatic infectious individuals at time t i.e., people who have contracted the disease but have not yet developed it at time t, • R(t) is the number of reported symptomatic infectious individuals i.e., symptomatic infectious with severe symptoms at time t, • U(t) is the number of unreported symptomatic infectious individuals i.e., symptomatic infectious with mild symptoms at time.

Introduction
This study aims at analyzing the evolution of the COVID-19 epidemic in the great Abidjan over the peri-   (1.1) τ , ν , 1 ν , 2 ν and η are real constants on [0;1]. τ is the transmission rate. ν is the rate at which the asymptomatic infectious individuals become symptomatic. They assume 1 ν is the average time during which asymptomatic infectious are asymptomatic.
f is the fraction of symptomatic infectious that become reported symptomatic infectious. So 1 f ν ν = is the rate at which asymptomatic infectious become reported symptomatic infectious.
1 f − is obviously the fraction of symptomatic infectious that become unreported symptomatic infectious. So is the rate at which asymptomatic infectious become unreported symptomatic infectious.
η is the rate at which the symptomatic infectious individuals lose their symptoms by death or healing. They get 1 η the average time symptomatic infectious have symptoms.
Pierre Magal, et al. calculated the parameters of their model from the Wuhan epidemic data before public confinement and isolation measures.
A nonlinear adjustment method allowed them to obtain an expression of the CR function giving the cumulative number of reported symptomatic infectious individual at time t. They obtain that With numerical simulations they show the evolution of the epidemic in Wuhan without public confinement and isolation measures. By varying the parameters of the model, they use simulations to predict the peak of the epidemic.
Since January 23, public confinement and isolation measures have been taken in Wuhan. Pierre Magal, et al. considered the effectiveness of its measurements after January 25. They therefore modified their model with the following non constant transmission rate: With these public measures, numerical simulations show that the epidemic in Wuhan will be eradicated. We quote them: "As a consequence of our study, we note that public health measures, such as isolation, quarantine, and public closings, greatly reduce the final size of this epidemic, and make the turning point much earlier than without these measures." In our work, we use the first deterministic model of Pierre Magal, et al. (1.1) because we consider that it could reflect the dynamics of the epidemic in the Greater Abidjan Region as we will see later with the data of this region. We rephrase some starting notations. We study the system from an instant which we will define later. At this instant 0 t = we note 0 0 0 0 ( , , , ) S I R U the state of the system. N is the total population size before the epidemic. In this model, there is no birth or recruitment. So, we have 0 In section 1, we prove that the solution of the system (1.1) is bounded in is the set of the different disease-free equilibrium for the system  [2,3], we show that all the solutions of (1.3), bounded for 0 t ≥ in ∆ , converge in M when t → +∞ . These mathematical results comfort us in the possibility of reaching a state without the disease.
In Côte d'Ivoire, the Greater Abidjan Region includes the city of Abidjan and six other districts which are peripheral.
In Côte d'Ivoire, the first patient of COVID-19 was detected in Abidjan on March 11, 2020. At the start of the epidemic, the rest of the country was unaffected by the disease. So the public authorities took the following measures [4].
• The suspension of air flights to countries with more than 100 reported cases (decided on March 20, 2020).
• The closure of all land, sea and air borders from midnight on March 22, 2020 until further notice.
• On March 23, 2020 the state of emergency is established. A curfew is in place from 9 p.m. to 5 a.m. for two weeks.
• From March 29, 2020, no transport is authorized between the Greater Abidjan region and the rest of the country.
• On April 9, 2020, it was decided: the wearing of the mandatory mask in the Greater Abidjan; compulsory confinement of fragile people; reduction of non-essential travels and the number of people in vehicles.
The Greater Abidjan Region is the epicenter of the epidemic with 98 percent of the confirmed cases. Since March 29 the big Abidjan has been isolated from the rest of the country. The rest of the country has no more than 5 reported cases. Since at least April 21, 2020, the rest of the country has no COVID-19 patients. National data in [4] can therefore be assimilated to that of the Big Abidjan (Table 1, Table 2, Table 3, Table 4, Table 5 and Table 6).
The majority of the economic activities of the Ivorian population are activities of the informal sector (small trade, rustic restaurant, street vendors, illegal taxi drivers,...). These activities mainly focus on the city of Abidjan. This did not allow the establishment of total containment in Côte d'Ivoire. In addition to conventional barrier measures    (washing hands regularly, coughing in the elbow...), the measures above have been gradually implemented to curb the epidemic. As our study is limited to the Greater Abidjan Region, we will analyze the data from the isolation of the Greater Abidjan Region (March 30, 2020). We will therefore study our model with the eradication measures implemented.
On the other hand, from this isolation date, the epidemic data are more reliable because of the multiplication of screening centers.
These measures, which do not correspond to total containment, favor the transmission of the virus. The transmission rate of the Greater Abidjan Region model is therefore not zero.
In section 2, as was done in [1,5,6], we will determine the constants of the system (1.1) corresponding to the epidemic in the Great Abidjan Region with the effect of the public measures. We consider the period of the effect of the public measures for two reasons: daily statistical data is more significant during this period because of the multiplication of screening centers and we want to know the behavior of the epidemic with these measures. We use the software CurveExpert Pro to fit our data. That permits us to compute the different parameters of the Greater Abidjan Region model and the basic reproductive number ℜ 0 . Simulation with MatLab software permits us to find the peak of the epidemic in the Greater Abidjan Region in the case of this deterministic model. Simulations with other parameter values make it possible to observe a more rapid eradication of the epidemic and to propose measures to achieve this.
In Côte d'Ivoire, continuously, certain protective measures such as minimum physical distance of one meter, regular hand washing etc. are not respected in markets and public transport vehicles. Abidjan has many large markets in each of its districts. Our region has many uncontrolled urban transport that does not allow compliance with health measures. Also, the practice of economic activities in the informal sector promotes many contacts between individuals. All these factors continually disturb the transmission rate of the model and could be assimilated to white noise to be introduced into the system (1.1). That's why we add a white noise to the model (1.1) to disrupt the transmission rate.
The isolation of the Greater Abidjan Region is not fully respected despite the efforts of the police. Sudden inflows or outflows of people take place formally with administrative authorizations or illegally by village roads.
Sudden mid-term measures such as the obligation to wear a mask and a curfew extension caused sudden shocks to the model at the first time when they were applied.
The model described by the system (1.1) does not take these factors into account. We need to introduce jumps in the system (1.1) to represent these sudden changes which disturb the transmission rate.
( X t -) means the left limit of ( ) X t , ( , ) N dt dz  is a Poisson counting measure with the stationary compensator In Section 3, we prove that the system (1.2) has a unique global and positive solution. We show how the jumps influence the behavior of the solution around the set of the equilibrium states. Inspired by [10,11] and using MaLab R2019b software, we perform numerical simulations which corroborate our theoretical results.

Some Mathematical Properties of the Deterministic Model
Here we justify simple mathematical properties that confer biological significance to our model. N is the size of the population before the epidemic. We have 0 0 0 0 S I R U N + + + ≤ . We prove that for any initial state in , the solution of the system (1.1) are bounded in

Lemma
For any given initial value in , the set ∆ is positively invariant for the system (1.1).

Proof
With the system (1.1), we obtain by summing the equation of the system: If we prove that for will be decreasing and therefore less than or equal to that N .
On the other hand, initiating the solution of system (1.1) in positive orthant 4 +  , we have :

Proof
We define : Since all the model parameters are non-negative, 0 S S ≤ and 1 1 ν ≤ , we get :

Remark
In practical terms, when 1 ℜ ≤ 0 , any solution of the system (1.1) converges towards a disease-free state.

Determination of the constants in the model of the COVID-19 epidemic in the Greater Abidjan Region with the effect of public measures
In this section, we will try to determine the parameters of the deterministic model corresponding to the dynamics of the epidemic in the Greater Abidjan Region after the implementation of all the public measures that led to its isolation.  All this brings us to consider our data 14 days after the isolation of the Greater Abidjan Region to tend towards the dynamics of the epidemic after all the public measures have been taken.
The initial instant 0 t = of our study of the epidemic in the big Abidjan with the effect of public measures therefore corresponds to April 12, 2020. In our case we will name 0 0 0 0 ( , , , ) S I R U the state of the system at this time. This is the initial state of the model that we are going to study.
We will not use the value of f indicated in [1] by Pierre Magal, et al. We will consider that 0.05 f = . It is assumed in [12] and [13] that the young age of an individual favors the absence of symptoms or the development of very weak symptoms after COVID-19 contraction. The last population census of Côte d'Ivoire [14] shows that children (0-14 years-old) represent 41.8 percent of the total population and young people (15-34 years-old completed) constitute 35.5 percent of the total population. 77.3 percent of the total population is under 35 years of age. This extreme youth of the population could explain the low number of symptomatic infectious because the symptoms of the disease rarely appear. On the other hand, the low attendance of health facilities (around 47 percent, [15]) may not favor the detection of COVID-19. Screening is only systematic in the event of serious symptoms and symptomatic infectious are hardly detected. There are not enough screening centers, although their number has been increased. For all these reasons, we will assume that around 5 percent of infectious individuals become confirmed cases.

N =
is the size of the population of the big Abidjan before the epidemic.
As in [1], we define the CR function of the cumulative number of reported symptomatic infectious cases at time t for 0 t ≥ with the function The statistical data of the COVID-19 epidemic show that the variations in the number of susceptible people due to the epidemic in the beginning of the epidemic are so small so we can consider 0 ( ) S t N = and 0 (0) 5816060 S S = = because as we will see later the distance between 0 and 0 t is small. That will permit us to evaluate the other parameters of the model. We can later adjust the value of 0 S .
As we look at the epidemic from the effect of the measures at date t = 0 then, for 0 0 t t ≤ < , ( ) CR t exists but has no interesting biological significance and is used just for calculations ( Figure 1).
To exploit this adjustment function, we consider that For When we use the methods to estimate the parameters presented and justified in [1] in section 4 we obtain that :  Under the conditions we have adopted, the epidemic in the Greater Abidjan lives around 1400 days after April 12, 2020, with a peak reaching 230 days after this date.
We can see that at the end of the epidemic, the number of susceptible people in the population of the big Abidjan is around 5,300,000. This drop should not be seen as deaths but rather as cures because the death rate from the disease is very weak (it is less than 1 percent as the statistics indicate. Refer also to the linear graphs). At the end of the epidemic, we have approximately 300,000 cured which does not influence the dynamics of the model.

Tracks for rapid eradication of the epidemic in the great Abidjan
In support of the public measures taken, we agree with the two new directions taken by the government of Côte d'Ivoire: An increase in screening tests results in an increase in the parameter f and therefore in the rate 1 ν .
We suggest a large increase from f to 0.8. We hope that 80 percent of the infected will be reported cases in order to get them out of the dynamics of the epidemic.
The degree of repression and sanction can result in a coefficient k between 0 and 1 by which the current transmission rate is multiplied. If the repression and the sanctions are high then k is small and the rate of transmission decreases. For our simulation, we opt for k = 0.5, i.e. a reduction by half of the current transmission rate. With the previous values of the other parameters, using the same methods of parameter approximations the following simulation gives 0.5131 ℜ = 0 ( Figure 4).
The following simulations shows in the same conditions the evolution of the number of individuals susceptible which knows a small decrease to stabilize at 5,815,600 at the end of the epidemic ( Figure 5).

The Stochastic Model with Jumps
In this section, we study the model (1.2) for 0 t ≥ where 0 t = represents the date of April 12, 2020, the date from which we assume that the public measures are effective.
Lyapunov analysis [7,8,16], permits us to prove that the solution of model (1.2) is positive and global.
We assume that the jump coefficient satisfies the following conditions: We will first prove that the system (1.2) has a unique solution which is global and positive. In a second point we will study the asymptotic behavior of this solution in the case

Global positive solution
As in [16], we show that jump processes can suppress the explosion. To prove this local solution is global, we need to show that e τ = ∞ a.s. We will also prove that for any given At first we prove that the solution is positive. In second part, we prove that the solution checks is in this set. So we define the following stopping time Clearly, m τ is increasing as m ↑ ∞ a.s. Set

Z C z S t I t U t C z S t I t U t S t I t I t S t z
We define for each u ≥ 0, Due to the property of the function $ F $ defined by We know that if X is a Levy process then, for each 0 t ≥ we have S , R , I and $ R $ are Levy processes. So we obtain almost surely from (3.1) : Since T is arbitrary, we must have ( ) 1 τ ∞ = +∞ = Let ω be fixed in Ω .
is a deterministic function such that ω ω ω η ω ω is a decreasing function.
Therefore for By varying ω in the universe of probability , we thus show that: This completes our proof and for any given initial value Proof: We know that

t I t R t U t
LGdt We obtain Integrating both sides of (3.2) from 0 to t, then taking expectation, yields the solution of (1.2) oscillates more closely in the set M of the disease-free equilibrium for the system (1.1) as the intensity of the noise and the jumps decrease. The case fatality rate is less than 1 percent (see linear graphs). Convergence towards a disease-free state is an encouraging sign. It is a disease-free situation with few deaths.

Numerical Simulation of the Stocastic Model
As in [7,10] using MatLab R2019b software, we add white noise and Poisson jumps to the first example seen in the deterministic case which reflects the reality in the big Abidjan. To be close to the difficulties of lowering the transmission rate in greater Abidjan, we lowered this rate by 10 percent in the first two figures and by 25 percent in the last two. We get basic reproductive numbers less than 1.
We obtain the following figures which corroborate the results of our last theorem ( Figure 6). Now we present the evolution of susceptible individuals (Figure 7 and Figure 8). Now we present the evolution of susceptible individuals (Figure 9). Peak of epidemic appears to be around 600 days after April 12, 2020.

Conclusion and Perspectives
In this study, we tried to describe the behavior of the COVID-19 epidemic in the great Abidjan during the period of public measures. We use a stochastic model with white noise and jumps which reflects the disturbances that hamper the implementation of these measures. Total containment which would reduce the transmission rate to zero is impossible in the current context. The application of controls and actions restriction, would greatly reduce the transmission rate τ . Massive screening corresponds to a high value of the parameter f . This means a large value of 1 ν and a small value of 2 ν . This permits to reduce the propagation of the epidemic because this allows a large number of infectious individuals to escape from the dynamics of the epidemic. Therefore, the basic reproductive number