Issue 
Math. Model. Nat. Phenom.
Volume 15, 2020
Coronavirus: Scientific insights and societal aspects



Article Number  24  
Number of page(s)  6  
DOI  https://doi.org/10.1051/mmnp/2020006  
Published online  25 March 2020 
On a quarantine model of coronavirus infection and data analysis
^{1}
Institut Camille Jordan, UMR 5208 CNRS, University Lyon 1,
69622
Villeurbanne, France.
^{2}
INRIA Team Dracula, INRIA Lyon La Doua,
69603
Villeurbanne, France.
^{3}
Peoples Friendship University of Russia (RUDN University),
6 MiklukhoMaklaya St,
Moscow
117198, Russia.
^{4}
Department of Mathematics & Statistics, IIT Kanpur,
Kanpur
208016, India.
^{5}
School of Mathematics & Actuarial Science, University of Leicester,
Leicester
LE1 7RH, UK.
^{*} Corresponding author: volpert@math.univlyon1.fr
Received:
16
March
2020
Accepted:
16
March
2020
Attempts to curb the spread of coronavirus by introducing strict quarantine measures apparently have different effect in different countries: while the number of new cases has reportedly decreased in China and South Korea, it still exhibit significant growth in Italy and other countries across Europe. In this brief note, we endeavour to assess the efficiency of quarantine measures by means of mathematical modelling. Instead of the classical SIR model, we introduce a new model of infection progression under the assumption that all infected individual are isolated after the incubation period in such a way that they cannot infect other people. Disease progression in this model is determined by the basic reproduction number R_{0} (the number of newly infected individuals during the incubation period), which is different compared to that for the standard SIR model. If R_{0} > 1, then the number of latently infected individuals exponentially grows. However, if R_{0} < 1 (e.g. due to quarantine measures and contact restrictions imposed by public authorities), then the number of infected decays exponentially. We then consider the available data on the disease development in different countries to show that there are three possible patterns: growth dynamics, growthdecays dynamics, and patchy dynamics (growthdecaygrowth). Analysis of the data in China and Korea shows that the peak of infection (maximum of daily cases) is reached about 10 days after the restricting measures are introduced. During this period of time, the growth rate of the total number of infected was gradually decreasing. However, the growth rate remains exponential in Italy. Arguably, it suggests that the introduced quarantine is not sufficient and stricter measures are needed.
Mathematics Subject Classification: 92D30
Key words: Coronavirus infection / quarantine model / modes of infection development
© The authors. Published by EDP Sciences, 2020
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
1 SIR model
The classical susceptibleinfectedrecovered (SIR) model in epidemiology (see, e.g., [2]) allows the determination of critical condition of disease development in the population irrespective of total population size. Among many versions of the model, let us consider the following ODE system (1.1) (1.2) (1.3)
for the sizes (numbers) of the susceptible subpopulation S, infected I, and recovered R. The term kIS describes the disease transmission rate due to the contact with infected, βI characterizes the rate of recovery of infected, and σI mortality due to infection only. It is assumed here that recovered do not return to susceptible class, that is, recovered individuals have immunity against new infection, and they cannot infect susceptible.
Analysis of this model is well known and we will not discuss it here. Let us only determine the condition of the disease progression in the case where the number of infected/recovered/dead is much less than the number of susceptible, and S in the model above can be approximate by a constant, S ≈ S_{0}. The equation (1.2) becomes as follows: (1.4)
We obtain an ordinary differential equation with constant coefficients whose solution can be readily found, I(t) = I_{0}e^{μt}, where I_{0} is the number of infected at the initial moment of detection of the infection/disease. Substituting this expression into equation (1.4), we get
where the new parameter is called the basic reproduction number. If , then the number of infected will grow, if it will decay which plays the key role behind disease propagation.
In particular, one can change the course of the disease dynamics by changing the basic reproduction number. For instance, if measures are introduced that push it below one (e.g. by making k to decrease), the exponential growth changes to exponential decay. Interestingly, this is exactly what has happened in South Korea (see Sect. 3) after restrictions on daily life were introduced.
Thus, the classical SIR model determines the condition of the disease development from the comparison of the disease transmission rate with the sum of the recover and death rates. In the other words, we compare the rates of adding and removal of infected.
The model does not account for the incubation period of the disease that was shown to be important in the case of coronavirus spread: individuals can become infective before showing any symptoms. A model that allows for the incubation period is considered in the next section.
2 Quarantine model
2.1 Model
Exceptional measures adopted for the coronavirus infection suggest to introduce another model of infection development. We consider the subpopulation of latently infected individuals who are already infected but do not show any symptom during the incubation period. When the incubation period is over, the disease manifests itself with its symptoms, and the individual is isolated in the quarantine where he/she cannot infect the others. Under these assumptions, instead of system (1.1)–(1.3) we have (2.1) (2.2)
where I is the subpopulation of latently infected individuals, τ is the incubation period. The second term in the righthand side of equation (2.2) corresponds to the individuals infected at time t − τ. Their incubation period is finished at time t, and they are put to quarantine. As a result they cannot spread infection anymore.
2.2 Solution
As before, we approximate S(t) ≈ S_{0}, S(t −τ) ≈ S_{0} as τ is not large, and replace equation (2.2) by the equation (2.3)
Substituting I(t) = I_{0}e^{μt} in above equation, we get (2.4)
This equation has solution μ = 0 for all values of parameters. Besides this solution, it can have a positive or a negative solution μ. The function f(μ) = kS_{0}(1 − e^{−μτ}) is an increasing function with the asymptotic limit kS_{0} at infinity. The sign of the solution is determined by the derivative f′ (0). If f′ (0) > 1, then there is a positive solution (Fig. 1), if f′(0) < 1, the solution is negative. In terms of the basic reproduction number , the condition is similar as for the SIR model. However, the meaning of this parameter is different. It characterizes the total infection rate during the incubation period and not its ratio with the rates of recover and death.
With the notation , we can write equation (2.4) as follows: (2.5)
Its solution depends on the single parameter .
Figure 1 Graphical solution of equation (2.4), μ = f(μ) with the values of parameters: τ = 1, kS_{0} = 2 (dotted curve) and kS_{0} = 0.4 (dashed curve). 
2.3 Public health
In the absence of vaccination and lack of effective treatment, the only way to influence the disease development is to act on the basic reproduction number , that is to decrease the value of the parameter k characterizing the disease transmission rate by the infected individuals (among which the symptoms are not prominent yet). Note that S_{0} is a characteristic of the given population and τ is the property of the disease, therefore neither of them is controllable. However, parameter k depends on the human mobility and the social distance and hence it can be changed. If in the beginning of the disease development k > 1∕(S_{0}τ), that is each infected individual infects on average more than one susceptible individual during the incubation period, then the number of newly infected individuals will grow exponentially. If, at some moment of time, there are measures restricting potential contacts of infected and susceptible individuals, then k can become less than the critical value, and the number of newly infected individuals will exponentially decay. There is a transition between the growing and decaying exponential.
2.4 Limitations of the model
The model has a number of obvious limitations. It takes into account only the beginning of the disease development where the number of suspected can be considered as constant. It is a good approximation if the disease propagation is stopped/regulated when the number of infected is relatively small (compared to the total population).
Next, the model does not take into account the spatial distribution of infected and their displacement. In this case, instead of equation (2.3) one has to consider a model with an explicit space, for instance the following equation: (2.6)
where I and S_{0} now have the meaning of the corresponding densities instead of sizes and the diffusion term describes random motion of the individuals with intensity δ. In the case where neither of S_{0}, k and τdepend on space, equation (2.6) can be reduced to the nonspatial model. If we consider, for simplicity, the whole axis (or a bounded interval with noflux boundary conditions), then we introduce the variable . Integrating equation (2.6) with respect to x, we obtain an equation similar to (2.3).
However, the spatially averaged model does not describe the dynamics well in case some of the parameters depend on space. If there are two different patches of the disease development with different basic reproduction numbers, then the disease can be eradicated in the first patch due to the imposed restrictions but it can give a new outbreak in another patch if the restrictions are not adopted there or they are not sufficient. This situation is observed in Europe where the disease progresses exponentially while it is already slowed down in China. Hence, in some cases, a multipatch model should be considered with some connectivity between the patches at least for some time period. Also, the human movement can follow a pattern more complicated than the one given by the Fickian diffusion, e.g. to follow the network made by longdistance connections between large airports.
Another important limitation of the model is related to the assumption of a single incubation period. According to the data, it is possible that there are different incubation periods or maybe the whole spectrum of incubation periods from several days up to four weeks. The distributed delay models are appropriate in this case.
Finally, let us mention virus mutations that can have a strong influence on the disease progression and treatment [1]. At the moment, there are no available data on mutations of coronavirus, it may take some time to have confirmation of this aspect.
3 Data
We used the data from Worldometer^{1} showing the total number of infected individuals and daily cases in different countries and worldwide. Across the world daily cases clearly show twopatch dynamics basically due to China (growthdecay) and Europe (growth). There are periods of exponential growth, decay and regrowth. The origin of the outbreak on February 12 (in China) is not clear yet. It may be related to the method of data collection. Daily reported cases curves in China and South Korea correspond to the growthdecay dynamics. The daily new cases curves for France, Germany, Italy, Spain, USA correspond to the growth mode.
Let us note that the data on daily reported cases do not exactly correspond to the variable I(t). The latter represents a sum of daily cases during the incubation period. Taking into consideration that the incubation period to be of 5 days, we obtain the graphs for the latently infected individuals. We observe growthdecay dynamics for South Korea (Fig. 2) and growth dynamics in Italy (Fig. 3).
Fitting the data allows us to determine the basic reproduction number from equation (2.5). For South Korea, on the growing branch and on the decaying branch. For Italy, . Let us note that the basic reproduction number depends on the incubation period.
Figure 2 Latently infected individuals in South Korea obtained as a sum of daily cases during the incubation period taken equal 5 days. The growing branch is approximated by the exponential 200e^{0.21x}, decaying exponential by 35 ⋅ 10^{3}e^{−0.15x}. 
Figure 3 Latently infected individuals in Italy obtained as a sum of daily cases during the incubationperiod taken equal 5 days. The data curve is approximated by the exponential 230e^{0.18x}. 
4 Discussion
The model confirms the efficiency of the approach to stop the disease by the limitation of contacts between the individuals through quarantine of infected individuals. This is quite obvious in theory, with the help of simplest model formulation,but difficult in practice. Experience of China and South Korea shows that the peak of infection (maximum of newly reported cases on daily basis) is reached about 10 days after adopting serious restrictive measures. The number of infection increased during this time in 10–20 times. In Italy 10 days after the universities and schools were closed (March 4) the peak of infection does not seem to be reached, and exponential growth continues.
Moreover, the exponential growth rate of the total number of infected in China and in South Korea observed before the adoptedmeasures (January 25 and February 22, respectively) rapidly changed to a slower growth rate afterward (Fig. 4). Similar situation is observed in Iran though the information about adopted measures is not clear yet. However, in Italy the exponential growth rate does not change up to March 4. This can bean indication that these measures are not sufficient (or that they are not respected by local people).
Figure 4 Total number of infected in the logarithmic scale. In South Korea (left), the exponential growth rate (linear growth in this scale) before February 22 is replaced by a slower growth after this date. In Italy (right), the exponential growth rate (linear growth in this scale) after March 4 does not change. 
Acknowledgements
The first author acknowledges the IHES visiting program during which this work was done. The work was supported by the Ministry of Science and Education of Russian Federation, project number FSSF20200018, and by the FrenchRussian program PRC2307.
References
 N. Bessonov, G. Bocharov, C.C. Leon, V. Popov and V. Volpert, Genotype dependent virus distribution and competition of virus strains. To appear in Memocs (2020). [Google Scholar]
 J. Murray, Vol.1 of Mathematical biology, third edition. SpringerVerlag, Heidelberg (2002). [Google Scholar]
All Figures
Figure 1 Graphical solution of equation (2.4), μ = f(μ) with the values of parameters: τ = 1, kS_{0} = 2 (dotted curve) and kS_{0} = 0.4 (dashed curve). 

In the text 
Figure 2 Latently infected individuals in South Korea obtained as a sum of daily cases during the incubation period taken equal 5 days. The growing branch is approximated by the exponential 200e^{0.21x}, decaying exponential by 35 ⋅ 10^{3}e^{−0.15x}. 

In the text 
Figure 3 Latently infected individuals in Italy obtained as a sum of daily cases during the incubationperiod taken equal 5 days. The data curve is approximated by the exponential 230e^{0.18x}. 

In the text 
Figure 4 Total number of infected in the logarithmic scale. In South Korea (left), the exponential growth rate (linear growth in this scale) before February 22 is replaced by a slower growth after this date. In Italy (right), the exponential growth rate (linear growth in this scale) after March 4 does not change. 

In the text 
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.