1 Introduction

According to the World Health OrganizationFootnote 1 (WHO), the severe acute respiratory syndrome corona virus 2 (SARS-CoV-2) previously known as new corona virus 2019 (2019-nCoV) is responsible of an infectious disease called Covid-19 [67]. In 1937, corona viruses were first identified as infectious bronchitis viruses with which birds suffered, that could devastate poultry stocks. Today, the viruses are the cause of the common cold in \( 15\% \) to \( 30\% \) of all cases. In the past 70 years, researchers have found camels, cattle, cats, dogs, horses, mice, pigs, rats and turkeys that were infected with corona viruses.Footnote 2 Older people, and those with underlying medical problems like cardiovascular diseases, diabetes, chronic respiratory diseases and cancer are more likely to develop serious illness. The Centers for Disease Control and PreventionFootnote 3 (CDC) affirm that the most common ways the virus spreads from an infected person to healthy people around them is when they cough or sneeze and release viral particles into the air and through touching, hands shaking and other forms of close personal contact. When healthy people touch objects or surfaces on which there are viral particles, then touch their eyes, nose or mouth before washing their hands, the virus can spread. In some rare cases, fecal contamination can cause the virus to spread as well. The best way to prevent and slow down transmission is being well informed about the SARS-CoV-2, the disease, its causes and how it spreads. Each person protects himself and others from infection by hands washing or using an alcohol-based rub frequently and not touching his face. The SARS-CoV-2 spreads primarily through droplets of saliva or discharge from the nose when an infected person coughs or sneezes, so it is important to also practice respiratory etiquette (for example, by coughing into a flexed elbow). Depending on its location, the SARS-CoV-2 can leave up to 7 days out of human body [17]. Several promising vaccines have recently emerged even though their efficiency is not fully established [43]. However, their adoption by populations is very partial, especially in the case of Cameroon where awareness is a major challenge. This slow uptake of vaccine solutions can be explained in part by the fact that vaccination is not a substitute for the application of barrier measures, even for those who are vaccinated. In this context, modeling, analyzing models, estimating risks and forecasting the potential spread of the disease in population remain very useful for decision makers [42].

Recently, numerous papers (published or not) appeared in order to contribute to the fight against the pandemic Covid-19. Those papers can be organized in three groups. The first group addressed the problem of forecasting the disease in order to help decision makers to better evaluate the logistic challenges they will face [3, 8, 9, 23, 24, 32, 34, 37, 38, 50, 52, 53, 55, 62, 65, 66, 69, 83, 84]. The second category of papers focused on evaluating the effectiveness of mitigations measures prescribed by the WHO and different governments in order to define better fighting strategies [10, 16, 19, 21, 22, 33, 38, 45, 46, 50, 55, 61,62,63, 68, 71,72,73,74,75, 77, 82, 84]. The last set of papers studied the social and the economical impacts of the pandemic [47, 81]. As we can observe in the literature, the wide part of models only considers person to person disease transmission. Of course it is important, but indirect transmission by environment may also be considered. On the other hand, the virus population dynamics appears also important and should be explicitly studied.

Indeed, when an individual gets in contact with the virus he is not directly exposed. According to the classic epidemiological literature, an exposed person refers to an infected person who cannot yet transmit the disease [7, 14, 51]. In the context of Covid-19, it refers to an individual whose virus is present in his body, but the level of reproduction has not yet reached the stage where the individual is shedding virus through various channels (coughing and sneezing). Hence, to become exposed the virus has to penetrate in the organism by oral or respiratory ways. Therefore, the exposition is mainly indirect in terms of human-to-human, or human-to-environment contacts. This consideration justifies the regular disinfection of individuals and the environment as recommended by public health agencies.

Following the authors in [27,28,29,30,31, 57], it is important to build and understand models that link the within-host dynamics and population-level dynamics of infectious. The prediction of the behavior of Covid-19 is very difficult because mitigation strategies are permanently changing. So it could be more efficient given an initial situation to propose optimal control strategies that will lead to disease extinction at a time to be determined. Hence, this work proposes efficient and low-cost control strategies against Covid-19.

The paper is organized as follows. A model is described in Sect. 2.2 according to the notations given in Sect. 2.1 and its asymptotic behavior is studied in Sect. 2.3 depending on some critical thresholds we define. In Sect. 3.1, we estimate the parameters of the model according to the real available data from Cameroon. A sensitivity analysis of critical thresholds depending on some control parameters is carried out in Sect. 3.2. Section 4 is devoted to the design and the computation of optimal control strategies according to different constraints. Illustrative simulations are carried out in Sect. 4.2. The algorithm which computes the optimal control is given in appendix B. The proofs of different theoretical result in “Appendix A.” The paper ends with a conclusion in Sect. 5.

2 The model and its general features

2.1 Definition of parameters

In this short section, different parameters are presented in order to make easier the understanding of the model to be proposed and studied later.

Table 1 Description of the parameters of the model (22)

2.2 Description of the model

We can summarize the dynamics of Covid-19 by Fig. 1.

Fig. 1
figure 1

Flowchart of the \(SEIQHTRDD_{I}B-V_{S}V_{E}V_{I}V_{Q}V_{H}V_{T}V_{R}V_{D}V_{D_{I}}V_{F}\) compartmental model

The disease dynamics can be viewed under two angles. Indeed, there are human dynamics on the first hand and the viruses dynamics on humans or in the environment on the other. According to the Cameroonian context, we consider ten compartments in the human population per unit area:

  • Susceptible \( \left( S \right) \): Public health institutions refer to a susceptible person as a person who is not a carrier of the virus and has no immunity to the disease. Such a person may become infected if he or she does not comply with the barrier measures prescribed by the government.

  • Exposed (Infected but not contagious) \( \left( E \right) \): public health institutions refer to an exposed person as a person who has the virus in his or her body but is not contagious in the sense that the virus has not yet reproduced sufficiently in the body to be present in the body’s secretions such as saliva and nasal secretions. Such a person has no symptoms and will test negative for Corona virus.

  • Infectious \( \left( I \right) \): public health institutions refer to an infectious person as anyone who is a carrier of the disease and contagious (the virus has reproduced sufficiently in the body to be present in the body’s secretions such as saliva and nasal discharge), but who has not yet been identified as such, probably because he or she does not manifest symptoms (asymptomatic case) or has not yet been tested. People who are infectious from this point of view are very often free to move around and represent a great risk of contamination.

  • Quarantined \( \left( Q \right) \): public health institutions refer to a person in quarantine as any person isolated (at home or in a specialized center) after having been identified or suspected as a carrier of the disease and contagious (infectious), but not presenting a serious form of the disease endangering his or her life (no respiratory assistance for example). The follow-up of these people being particular (at least because of their isolation, even it cannot be perfect due to contacts with health agents who are also in contact with the rest of the population), it seems important to distinguish them explicitly by a compartment. This distinction allows on the one hand to take into account the effect of isolation on the dynamics of the disease. On the other hand, it allows the economic impact of quarantine to be explicitly evaluated.

  • Hospitalized \( \left( H \right) \): public health institutions refer to a hospitalized person as any infectious person (infectious or in quarantine) who manifests severe forms of the disease to the point of requiring special medical attention (such as respiratory assistance and oxygen).

  • Treated \( \left( T \right) \): public health institutions refer to a treated person as any person who, after having been in the infectious phase (and possibly the quarantine or hospitalization phases), returns to a state where he or she is out of danger and is less contagious: even if the viruses remain in low concentration in the body, their reproduction rate decreases due to the body’s resistance or to any treatment. It is important to note that despite the controversy about the existence of a treatment, several treatment protocols are proposed either by conventional medicine, called modern medicine, or by local or traditional medicine using herbal teas (for those who do not go to the hospital and practice self-medication without having tested their serological status). It is common for a treated person to be sent back to his usual daily life (home) if he was followed in a health facility (out of quarantine or hospitalization). This is due to the limited resources available to monitor large numbers of patients.

  • Recovered \( \left( R \right) \): public health institutions refer to a cured person as anyone who, after having been in the treated phase, no longer has the virus in his or her body and has possibly developed a resistance (a priori temporary) to the virus. He/she is no longer contagious in the sense that he/she does not participate in the reproduction of the virus and his/her nasal or oral discharge is healthy.

  • Non-infectious dead \( \left( D \right) \): public health institutions refer to a non-infectious decedent as anyone who died for reasons unrelated to the Corona virus. The funeral procedures for such persons are classic and entrusted to their families (more or less long stay in the morgue, wake with the body, transfer of the body to its village and festive burialFootnote 4 [40]).

  • Infectious dead \( \left( D_{I} \right) \): Public health institutions refer to a deceased infectious person as a person who has died from Covid-19. Funeral procedures for such persons are very fast and ensured by public services: the body is sealed and is not given to the family to avoid any risk of secondary contamination in the handling of the corpse. Even if it is not proven that the corpses of these people are not contagious, no study to our knowledge shows the contrary. As recommended in [43], taking into account the effect of the deceased follow-up policy in the dynamics of the disease and the impacts of the disease and the economic aspects requires the explicit introduction of the compartments of the deceased in the modeling.

  • Buried \( \left( B \right) \): The buried compartment makes it possible to count burials and assess the impact of the rapid burial of infectious dead on the dynamics of the disease. It also makes it possible to assess the economic impact of deaths due to Covid-19.

WHO and the Government of Cameroon recommend several control measures that can be explicitly listed here as follows:

  • Limiting direct contacts with others through a control parameter \( \left( \tau _{3} \right) \): obviously, direct contact with a person who is contaminated with viruses leads to being contaminated as well.

  • Limiting movements through a control parameter \( \left( m\right) \): multiplying movements increases the risk of being in contact with the virus from environment or from a person after a direct contact.

  • Social distancing which is implemented through limitation of group size per unit area to a maximum value \( \left( \kappa \right) \): social distancing limits not mobility but the density of people per unit area. For example, guaranteeing a distance of one meter between individuals means that on a disk of radius 1 meter there is exactly one person. The implementation of social distancing does not affect the size of the total population but rather its spatial distribution. To our knowledge, no existing epidemiological model has implemented social distancing using a concept similar to “carrying capacity” which is rather common in ecological models.

  • Limiting the area of movements through a lock-down ratio \( \left( 1-\frac{{\mathcal {A}}}{{\mathcal {A}}_{G}}\right) \): the disease is better contained if its spatial spreading is limited by lock-down. Social distancing is not to be confused with lock-down. In the lock-down ratio \( 1-\frac{{\mathcal {A}}}{{\mathcal {A}}_{G}} \), \( {\mathcal {A}} \) denotes the “effective area” while \( {\mathcal {A}}_{G} \) is the global geographic area of study where people are living. The concept of effective area, which we also newly introduced, designates the surface of the spatial domain in which an individual is allowed to move in order to respect the rules of Lock-down. Indeed, lock-down results in small subpopulations that can be geographically delimited and do not interact. This is expected to facilitate the control of the disease.

  • Frequently hands washing and disinfection through a control parameter \( \left( \omega \right) \): by frequently hands washing and avoiding to touch the face (especially eyes, nose and mouth), the probability to get infected is clearly reduced.

  • Wearing mask and avoiding to touch the face through a control parameter \( \left( u \right) \): Wearing mask reduces both the risk of infection and the release of viruses by infectious, quarantined, hospitalized and treated individuals.

One of the specificities of our modeling is based on the fact that each individual regardless of his compartment of belonging is potentially a reservoir since the outside part of his body (including his clothes) can carry viruses that he would have acquired in contact with the contaminated environment or a contaminated person. On the other hand, the amount of virus available in a human “reservoir” depends on its compartment since people who are contagious beyond transporting viruses taken elsewhere also produce them. This consideration makes it possible to explicitly assess the impact of barrier and disinfection measures. Again, we emphasize on the fact that getting in contact with viruses does not mean to be infected, but the infection occurs when the viruses penetrate the organism by oral or respiratory ways. Thus, we have ten compartments for virus population per individual status group located in a unitary surface or per unit area of environment \( \left( {\mathrm{virus}}\,\times \,{\mathrm{m}}^{-2}\right) \): global viral load in susceptible compartment \( \left( V_{S} \right) \), global viral load in infected compartment \( \left( V_{E} \right) \), global viral load in infectious compartment \( \left( V_{I} \right) \), global viral load in quarantined compartment \( \left( V_{Q} \right) \), global viral load in hospitalized compartment \( \left( V_{H} \right) \), global viral load in treated compartment \( \left( V_{T} \right) \), global viral load in recovered compartment \( \left( V_{R} \right) \), global viral load in non-infectious dead compartment \( \left( V_{D} \right) \), global viral load in infectious dead compartment \( \left( V_{D_{I}} \right) \), global free load of viruses in the environment \( \left( V_{F} \right) \). The viral load is transferred from an individual being in compartment i to the environment (transfer of type I: \( \tau _{1,i} \)) or conversely from the environment to an individual being in compartment i (transfer of type II: \( \tau _{2,i} \)). In addition, there are transfers of external viral load (available on the body and clothing) between individuals when there are direct contacts (transfer of type III: \( \tau _{3,i} \)). Under the simplifying assumption of mass action, it can be assumed that the amount of free virus lost by an individual from one compartment after contact is distributed to the other compartments in proportion to their different fractions in the population. It is likely that the balance of viruses received and lost is deficient in compartments I, Q and H since they participate in the reproduction of viruses. Infectious, quarantined, hospitalized and treated individuals produce new viruses at a rate \( \pi _{i} \) with i indexing a specific compartment. A susceptible individual becomes infected when a virus passes from the outside of his body to the inside (transfer of type IV: \( \tau _{4} \)). Infected persons die at a rate \( \mu \) or become infectious with a rate \( \alpha _{1} \). An infectious case is reported at a rate q and dies at a rate \( \mu +d \). The probability that a person identified as infectious does not immediately develop a severe form of the disease and is therefore quarantined is p, otherwise he is admitted to hospital. With a rate \( \alpha _{2} \) a quarantined person develop a severe form of the disease and is therefore admitted to hospital. Infectious, quarantined and hospitalized individuals die at a rate \( \mu +d \) or become treated (partially recovered) at the respective rates of \( \gamma _{1} \), \( \gamma _{2} \) and \( \gamma _{3} \). A treated individual is out of danger but is still infectious and can die at a rate \( \mu \). Treated persons recover completely and become immune at a rate \( \gamma _{4} \). A recovered person losses his immunity at a rate \( \ell \) or dies at a rate \( \mu \). Indeed, Recovered people may die from natural death or from causes unrelated to Covid-19. Potentially acquired immunity has certainly been considered in several modeling articles but no precise study on its effectiveness and duration exists to our knowledge. It is therefore more appropriate, as is the case in the literature, to assume that it is not definitive [42]. Moreover, despite the advent of vaccines, barrier measures remain strongly recommended for everyone, including those who have been vaccinated. With a rate \( \rho _{1} \), a non-infectious dead person is buried while an infectious dead individual is buried at a rate \( \rho _{2} \). Indeed, infectious dead persons are potentially very contagious and it is recommended to bury them rapidly in order to limit contacts.

Let \( X=\left[ \begin{matrix} S,&E,&I,&Q,&H,&T,&R,&D,&D_{I},&B \end{matrix} \right] ^{T} \) be the vector of human population densities per unit area in each compartment, \( Y=\left[ V_{S}, \quad V_{E}, \quad V_{I}, \quad V_{Q}, \quad V_{H}, V_{T}, V_{R}, \quad V_{D}, \quad V_{D_{I}}, \quad V_{F} \right] ^{T} \) be the vector of virus population densities per unit area in each compartment,

\( \Lambda =\left[ \begin{matrix} \Lambda _{1},&\Lambda _{2},&\Lambda _{3},&\Lambda _{4},&\Lambda _{5},&\Lambda _{6},&\Lambda _{7},&\Lambda _{8},&\Lambda _{9},&0 \end{matrix} \right] ^{T} \) be the vector of horizontal immigration rates in each compartment of the human population, \( \mathbf{b} =\left[ b_{1}, \quad b_{2}, b_{3}, b_{4}, \quad b_{5}, \quad b_{6}, \quad b_{7}, \quad 0, \quad 0, \quad 0 \right] ^{T} \) be the vector of vertical immigration rates (through births) in each compartment of the human population, \( \Pi = \left[ 0, \, 0, \, \varsigma _{3}\pi _{3}, \varsigma _{4}\pi _{4}, \, \varsigma _{5}\pi _{5}, \, \varsigma _{6}\pi _{6}, \, 0, \, 0, \, 0, \, 0\right] ^{T} \) be the vector of the release rates of viruses by contagious humans (via their nasal or oral discharge) into the environment, and \( {\overline{N}}=\sum _{i=1}^{7}X_{i} \text { and } N=\sum _{i=1}^{9}X_{i} \) be, respectively, the density of the living human population per unit area and the density of the living or deceased (but not yet buried) human population per unit area. A mathematical model corresponding to Fig. 1 can be described as

  • A mesoscopic model (intermediate scale of human gatherings):

    (1)

    The first term to the right of the equality in equation (1) corresponds to horizontal recruitment (immigration at rate \( \Lambda _{1} \)) and vertical recruitment (birth at rates \( b_{i} \)) in the susceptible human population. Population density per unit of surface is restricted by logistic expression \( \left( 1-\frac{N}{\kappa } \right) \) which reduce the population growth when N is near to \( \kappa \). The limiting capacity \( \kappa \) models social distancing which reduces the number of individuals per unit of area. The second term represents the exposition (\( \dfrac{a_{1}\tau _{4}Y_{1}}{a_{2}+Y_{1}} \)) and the natural death at rate \( \mu \). The last term corresponds to the lost of immunity at rate \( \ell \). Considering all those features in the same model are not common and it makes the analysis of the model more difficult.

    (2)

    The first term to the right of the equality in equation (2) corresponds to horizontal recruitment in the exposed human population. The second term represents the arrival of new exposed. The last term corresponds to the natural death at rate \( \mu \) and the worsening of the disease at rate \( \alpha _{1} \).

    $$\begin{aligned}&\overset{\bullet }{X_{3}}=\Lambda _{3}\left( 1-\frac{N}{\kappa } \right) +\alpha _{1}X_{2}-\left( \mu +d+q+\gamma _{1} \right) X_{3}. \end{aligned}$$
    (3)

    The first term to the right of the equality in equation (3) corresponds to horizontal recruitment in the infectious human population. The second term represents the arrival of new infectious. The last term corresponds to the death at rate \( \mu +d \), the detection of the disease by screening test (at rate q) and the treatment at rate \( \gamma _{1} \).

    $$\begin{aligned}&\overset{\bullet }{X_{4}}=\Lambda _{4}\left( 1-\frac{N}{\kappa } \right) +pqX_{3}\nonumber \\&\qquad -\left( \mu +d+\alpha _{2}+\gamma _{2} \right) X_{4}. \end{aligned}$$
    (4)

    The first term to the right of the equality in equation (4) corresponds to horizontal recruitment in the quarantined human population. The second term represents the arrival of new quarantined after testing an infectious at rate pq. The last term corresponds to the death at rate \( \mu +d \), the worsening of the disease (at rate \( \alpha _{2} \) and the treatment at rate \( \gamma _{2} \)).

    $$\begin{aligned}&\overset{\bullet }{X_{5}}=\Lambda _{5}\left( 1-\frac{N}{\kappa } \right) +\left( 1-p\right) qX_{3}+\alpha _{2}X_{4}\nonumber \\&\qquad \quad -\,\left( \mu +d+\gamma _{3} \right) X_{5}. \end{aligned}$$
    (5)

    The first term to the right of the equality in equation (5) corresponds to horizontal recruitment in the hospitalized human population. The second and third terms represent the arrival of new hospitalized after testing an infectious (at rate \( \left( 1-p\right) q\)) and disease worsening for quarantined (at rate \( \alpha _{2} \)). The last term corresponds to the death at rate \( \mu +d \) and the treatment at rate \( \gamma _{3} \).

    $$\begin{aligned}&\overset{\bullet }{X_{6}}=\Lambda _{6}\left( 1-\frac{N}{\kappa }\right) +\gamma _{1}X_{3}+\gamma _{2}X_{4}+\gamma _{3}X_{5}\nonumber \\&\qquad \quad -\,\left( \mu +\gamma _{4} \right) X_{6}. \end{aligned}$$
    (6)

    The first term to the right of the equality in equation (6) corresponds to horizontal recruitment in the treated human population. The second, the third and the fourth terms represent the respective treating rates from the infectious (\( \gamma _{1} \)), quarantined (\( \gamma _{2} \)) and hospitalized (\( \gamma _{3} \)). The last term corresponds to the death at rate \( \mu \) and the recovering at rate \( \gamma _{4} \).

    $$\begin{aligned}&\overset{\bullet }{X_{7}}=\Lambda _{7}\left( 1-\frac{N}{\kappa }\right) +\gamma _{4}X_{6}&\qquad -\left( \mu +\ell \right) X_{7}. \end{aligned}$$
    (7)

    The first term to the right of the equality in equation (7) corresponds to horizontal recruitment in the recovered human population. The second term represents the recovering rates of treated (\( \gamma _{4} \)). The last term corresponds to the death at rate \( \mu \) and the lost of immunity at rate \( \ell \).

    $$\begin{aligned}&\overset{\bullet }{X_{8}}=\Lambda _{8}\left( 1-\frac{N}{\kappa }\right) +\mu \left( X_{1}+X_{2}+X_{7}\right) -\rho _{1} X_{8}. \end{aligned}$$
    (8)

    The first term to the right of the equality in equation (8) corresponds to horizontal recruitment in the non-infectious dead human population. That recruitment could correspond to transportation or transfer of dead individuals for burial ceremonies like it is common in Cameroon or other valuable reasons. The second term represents the death of individuals in susceptible, exposed and treated compartment at rate \( \mu \). The last term corresponds to the burying at rate \( \rho _{1}\).

    $$\begin{aligned}&\overset{\bullet }{X_{9}}=\Lambda _{9}\left( 1-\frac{N}{\kappa }\right) \nonumber \\&\qquad +\left( \mu +d\right) \left( X_{3}+X_{4}+X_{5}\right) +\mu X_{6}-\rho _{2} X_{9}. \end{aligned}$$
    (9)

    The first term to the right of the equality in equation (9) corresponds to horizontal recruitment in the infectious dead human population. As in equation (8) that recruitment could correspond to transportation or transfer of dead individuals for burial ceremonies like it is common in Cameroon or other valuable reasons. The second term represents the death of individuals in infectious, quarantined, hospitalized at rate \( \mu +d \). The third term represents the death of individuals in recovered compartment \( \mu \). The last term corresponds to the burying at rate \( \rho _{2}\).

    $$\begin{aligned}&\overset{\bullet }{X_{10}}=\rho _{1} X_{8}+\rho _{2} X_{9}. \end{aligned}$$
    (10)

    The first and the second terms in Eq. (10) correspond the burying of dead individuals at respective rates \( \rho _{1}\) and \( \rho _{2}\). The horizontal migrations (represented by the \( \Lambda _ {i} \)s) in the different compartments can be explained by the movements of travelers who sometimes prefer to return to their countries of origin (as has been mentioned by the media about citizens of several African countries including Cameroonians residing in China), medical evacuations for more competent health centers, or repatriation of human remains. Thus, if the geographical area of study is limited, horizontal recruitments are possible in practically all the compartments.

  • A microscopic model (virus scale):

    $$\begin{aligned}&\overset{\bullet }{Y_{1}}=\left( \sum _{j=1}^{9}\tau _{3,j}Y_{j}+\tau _{2,1}Y_{10}\right) \frac{X_{1}}{N}\nonumber \\&\qquad \quad -\left( \mu +\nu +\omega _{1}+\tau _{1,1}+\tau _{3,1}+\frac{a_{1}\tau _{4}Y_{1}}{a_{2}+Y_{1}}\right) Y_{1}. \end{aligned}$$
    (11)

    The first term to the right of the equality in equation (11) represents the recruitment of viruses in susceptible population after a contact with environment (at rate \( \tau _{2,1} \)) and individuals of other compartments (at respective rates \( \tau _{3,j} \)). The second term corresponds to the disappearing of viruses by natural death (at rate \( \nu \)), by the death of the human host (at rate \( \mu \)), by disinfection measures applied by the host (at rate \( \omega _{1} \)), by the exchange of viruses with environment (at rate \( \tau _{1,1} \)) and other compartments (at rate \( \tau _{3,1} \)), and by the transition of the host from the susceptible compartment to the exposed compartment.

    $$\begin{aligned}&\overset{\bullet }{Y_{2}}=\left( \sum _{j=1}^{9}\tau _{3,j}Y_{j}+\tau _{2,2}Y_{10}\right) \frac{X_{2}}{N}+\frac{a_{1}\tau _{4}Y_{1}^{2}}{a_{2}+Y_{1}}\nonumber \\&\qquad \quad -\left( \mu +\alpha _{1}+\nu +\omega _{2}+\tau _{1,2}+\tau _{3,2} \right) Y_{2}. \end{aligned}$$
    (12)

    The first term to the right of the equality in equation (12) represents the recruitment of viruses in exposed population after a contact with environment (at rate \( \tau _{2,2} \)) and individuals of other compartments (at respective rates \( \tau _{3,j} \)). The second term corresponds to the update in population of viruses due to the transition from susceptible status to exposed status at rate \( \frac{a_{1}\tau _{4}Y_{1}}{a_{2}+Y_{1}} \). The third term represents the disappearing of viruses by natural death (at rate \( \nu \)), by the death of the human host (at rate \( \mu \)), by disinfection measures applied by the host (at rate \( \omega _{2} \)), by the exchange of viruses with environment (at rate \( \tau _{1,2} \)) and other compartments (at rate \( \tau _{3,2} \)), and by the transition of the host from the susceptible compartment to the exposed compartment at rate \( \alpha _{1} \).

    $$\begin{aligned}&\overset{\bullet }{Y_{3}}=\varsigma _{3}\pi _{3}X_{3}+\left( \sum _{j=1}^{9}\tau _{3,j}Y_{j}+\tau _{2,3}Y_{10}\right) \frac{X_{3}}{N}\nonumber \\&\qquad \quad +\alpha _{1}Y_{2}-\left( \mu +d+q+\gamma _{1}+\nu \right. \nonumber \\&\qquad \left. +\omega _{3}+\tau _{1,3}+\tau _{3,3} \right) Y_{3}. \end{aligned}$$
    (13)

    The first term to the right of the equality in equation (13) represents the release of internal viruses by infectious individuals on themselves. The second term corresponds to the recruitment of viruses in infectious population after a contact with environment (at rate \( \tau _{2,3} \)) and individuals of other compartments (at respective rates \( \tau _{3,j} \)). The third term corresponds to the update in population of viruses due to the transition from exposed status to infected status at rate \( \alpha _{1} \). The fourth term represents the disappearing of viruses by natural death (at rate \( \nu \)), by the death of the human host (at rate \( \mu +d \)), by disinfection measures applied by the host (at rate \( \omega _{3} \)), by the exchange of viruses with environment (at rate \( \tau _{1,3} \)) and other compartments (at rate \( \tau _{3,3} \)), by the transition of the host from the infectious compartment to the treated compartment at rate \( \gamma _{1} \), and by the transition of the host from the infectious compartment to the quarantined or the hospitalized compartment at screening rate q.

    $$\begin{aligned}&\overset{\bullet }{Y_{4}}=\varsigma _{4}\pi _{4}X_{4}+\left( \sum _{j=1}^{9}\tau _{3,j}Y_{j}+\tau _{2,4}Y_{10}\right) \frac{X_{4}}{N}+pqY_{3}\nonumber \\&\qquad \quad -\left( \mu +d+\alpha _{2}+\gamma _{2}+\nu +\omega _{4}+\tau _{1,4}+\tau _{3,4} \right) Y_{4}.\nonumber \\ \end{aligned}$$
    (14)

    The first term to the right of the equality in equation (14) represents the release of internal viruses by quarantined individuals on themselves. The second term corresponds to the recruitment of viruses in quarantined population after a contact with environment (at rate \( \tau _{2,4} \)) and individuals of other compartments (at respective rates \( \tau _{4,j} \)). The third term corresponds to the arrival of new quarantined from the infectious compartment at rate pq. The last term represents the disappearing of viruses by natural death (at rate \( \nu \)), by the death of the human host (at rate \( \mu +d \)), by disinfection measures applied by the host (at rate \( \omega _{4} \)), by the exchange of viruses with environment (at rate \( \tau _{1,4} \)) and other compartments (at rate \( \tau _{3,4} \)), by the transition of the host from the quarantined compartment to the hospitalized compartment at rate \( \alpha _{2} \), and by the transition of the host from the quarantined compartment to the treated compartment at rate \( \gamma _{2} \).

    $$\begin{aligned} \overset{\bullet }{Y_{5}}=&\varsigma _{5}\pi _{5}X_{5}+\left( \sum _{j=1}^{9}\tau _{3,j}Y_{j}+\tau _{2,5}Y_{10}\right) \frac{X_{5}}{N}\nonumber \\&+\left( 1-p\right) qY_{3}+\alpha _{2}Y_{4}\nonumber \\&-\left( \mu +d+\gamma _{3}+\nu +\omega _{5}+\tau _{1,5}+\tau _{3,5} \right) Y_{5}. \end{aligned}$$
    (15)

    The first term to the right of the equality in equation (15) represents the release of internal viruses by hospitalized individuals on themselves. The second term corresponds to the recruitment of viruses in hospitalized population after a contact with environment (at rate \( \tau _{2,5} \)) and individuals of other compartments (at respective rates \( \tau _{5,j} \)). The third term corresponds to the arrival of new hospitalized from the infectious compartment at rate \( \left( 1-p\right) q \). The fourth term materializes the arrival of new hospitalized from the quarantined compartment at rate \( \alpha _{2}\). The last term represents the disappearing of viruses by natural death (at rate \( \nu \)), by the death of the human host (at rate \( \mu +d \)), by disinfection measures applied by the host (at rate \( \omega _{5} \)), by the exchange of viruses with environment (at rate \( \tau _{1,5} \)) and other compartments (at rate \( \tau _{3,5} \)), and by the transition of the host from the hospitalized compartment to the treated compartment at rate \( \gamma _{3} \).

    $$\begin{aligned} \overset{\bullet }{Y_{6}}=&\varsigma _{6}\pi _{6}X_{6}+\left( \sum _{j=1}^{9}\tau _{3,j}Y_{j}+\tau _{2,6}Y_{10}\right) \frac{X_{6}}{N}\nonumber \\&\qquad \quad +\gamma _{1}Y_{3}+\gamma _{2}Y_{4}+\gamma _{3}Y_{5}\nonumber \\&\qquad \quad -\left( \mu +\gamma _{4}+\nu +\omega _{6}+\tau _{1,6}+\tau _{3,6} \right) Y_{6}. \end{aligned}$$
    (16)

    The first term to the right of the equality in equation (16) represents the release of internal viruses by treated individuals on themselves. The second term corresponds to the recruitment of viruses in treated population after a contact with environment (at rate \( \tau _{2,6} \)) and individuals of other compartments (at respective rates \( \tau _{6,j} \)). The third, the fourth and the fifth terms correspond to the update in viruses population due to the transition from infectious, quarantined and hospitalized compartments, respectively. The last term represents the disappearing of viruses by natural death (at rate \( \nu \)), by the death of the human host (at rate \( \mu +d \)), by disinfection measures applied by the host (at rate \( \omega _{6} \)), by the exchange of viruses with environment (at rate \( \tau _{1,6} \)) and other compartments (at rate \( \tau _{3,6} \)), and by the transition of the host from the treated compartment to the recovered compartment at rate \( \gamma _{4} \).

    $$\begin{aligned}&\overset{\bullet }{Y_{7}}=\left( \sum _{j=1}^{9}\tau _{3,j}Y_{j}+\tau _{2,7}Y_{10}\right) \frac{X_{6}}{N}+\gamma _{4}Y_{7}\nonumber \\&\qquad \quad -\left( \mu +\ell +\nu +\omega _{7}+\tau _{1,7}+\tau _{3,7} \right) Y_{7}. \end{aligned}$$
    (17)

    The first term to the right of the equality in equation (17) corresponds to the recruitment of viruses in recovered population after a contact with environment (at rate \( \tau _{2,7} \)) and individuals of other compartments (at respective rates \( \tau _{7,j} \)). The second term corresponds to the update in viruses population due to the transition from treated compartment. The last term represents the disappearing of viruses by natural death (at rate \( \nu \)), by the death of the human host (at rate \( \mu \)), by disinfection measures applied by the host (at rate \( \omega _{7} \)), by the exchange of viruses with environment (at rate \( \tau _{1,7} \)) and other compartments (at rate \( \tau _{3,7} \)), and by the transition of the host from the recovered compartment to the susceptible compartment at rate \( \ell \).

    $$\begin{aligned}&\overset{\bullet }{Y_{8}}=\left( \sum _{j=1}^{9}\tau _{3,j}Y_{j}+\tau _{2,8}Y_{10}\right) \frac{X_{8}}{N}+\mu \left( Y_{1}+Y_{2}+Y_{7}\right) \nonumber \\&\qquad \quad -\left( \rho _{1}+\nu +\omega _{8}+\tau _{1,8}+\tau _{3,8}\right) Y_{8}. \end{aligned}$$
    (18)

    The first term to the right of the equality in equation (18) corresponds to the recruitment of viruses in non-infectious dead population after a contact with environment (at rate \( \tau _{2,8} \)) and individuals of other compartments (at respective rates \( \tau _{8,j} \)). The second term corresponds to the update in viruses population due to the death of individuals belonging to susceptible, exposed and recovered compartments. The last term represents the disappearing of viruses by natural death (at rate \( \nu \)), by the burial of the dead host (at rate \( \rho _{1} \)), by disinfection measures applied on the dead host (at rate \( \omega _{8} \)), and by the exchange of viruses with environment (at rate \( \tau _{1,8} \)) and other compartments (at rate \( \tau _{3,8} \)).

    $$\begin{aligned}&\overset{\bullet }{Y_{9}}=\left( \sum _{j=1}^{9}\tau _{3,j}Y_{j}+\tau _{2,9}Y_{10}\right) \frac{X_{9}}{N}\nonumber \\&\qquad \quad +\left( \mu +d\right) \left( Y_{3}+Y_{4}+Y_{5}\right) +\mu Y_{6}\nonumber \\&\qquad \quad -\left( \rho _{2}+\nu +\omega _{9}+\tau _{1,9}+\tau _{3,9}\right) Y_{9}. \end{aligned}$$
    (19)

    The first term to the right of the equality in equation (19) corresponds to the recruitment of viruses in infectious dead population after a contact with environment (at rate \( \tau _{2,9} \)) and individuals of other compartments (at respective rates \( \tau _{9,j} \)). The second and the third terms correspond to the update in viruses population due to the death of individuals belonging to infectious, quarantined, hospitalized and treated compartments. The last term represents the disappearing of viruses by natural death (at rate \( \nu \)), by the burial of the dead host (at rate \( \rho _{2} \)), by disinfection measures applied on the dead host (at rate \( \omega _{9} \)), and by the exchange of viruses with environment (at rate \( \tau _{1,9} \)) and other compartments (at rate \( \tau _{3,9} \)).

    $$\begin{aligned}&\overset{\bullet }{Y_{10}}=\sum _{j=3}^{6}\left( 1-\varsigma _{j}\right) \pi _{j}X_{j}+\sum _{j=1}^{9}\tau _{1,j}Y_{j}\nonumber \\&\qquad \quad -\left( \nu +\omega _{10}+ \sum _{j=1}^{9}\frac{\tau _{2,j}X_{j}}{N}\right) Y_{10}. \end{aligned}$$
    (20)

    The first term to the right of the equality in equation (20) corresponds to the release of virus in the environment by infectious, quarantined, hospitalized and treated populations. The second terms represents the recruitment of viruses in environment after a contact with individuals of each compartment (at respective rates \( \tau _{1,j} \)). The last term represents the disappearing of viruses by natural death (at rate \( \nu \)), by disinfection measures (at rate \( \omega _{10} \)) and by the exchange of viruses from the environment to human hosts (at rate \(\sum _{j=1}^{9}\frac{\tau _{2,j}X_{j}}{N}\)).

The MonodFootnote 5 type infection force \( \uplambda =\frac{a_{1}\tau _{4}Y_{1}}{a_{2}+Y_{1}} \) is the product of the mass action factor \( a_{1} \), the adequate contact probability \( \frac{Y_{1}}{a_{2}+Y_{1}} \) (the probability that viruses reach critical parts of the body like nose, mouth and eyes), the transfer rate \( \tau _{4} \) and the susceptible body virus population size \( Y_{1} \).

The parameters of the model (1 )–(20) are described in Table 1. Although we focus explicitly on microscopic and mesoscopic scales, there is also a macroscopic scale of the global human population

$$\begin{aligned} X_{G}&=\left[ \begin{matrix} S_{G}&E_{G}&I_{G}&Q_{G}&H_{G}&T_{G}&R_{G}&D_{G}&D_{I,G}&B_{G} \end{matrix}\right] ^{T} \nonumber \\&= {\mathcal {A}}_{G}\times \left[ \begin{matrix} X_{1}&X_{2}&X_{3}&X_{4}&X_{5}&X_{6}&X_{7}&X_{8}&X_{9}&X_{10} \end{matrix}\right] ^{T}. \end{aligned}$$
(21)

Hence, equation at the macroscopic scale can be derived easily from the equation at the mesoscopic scale. Since the variables of the model (1 )–(20 ) are densities per unit area, the effect of lock-down is evaluated by multiplying these densities by the effective area \( {\mathcal {A}} \) instead of the total area \( {\mathcal {A}}_{G} \) like in (21). As mentioned before, lock-down allows the disease to be managed independently in small subpopulations. As we will see in Theorem 3, lock-down accelerates the extinction of the disease.

The model (1)–(20) can be written in the form

(22)

with M,Q Metzler and P having all its coefficients nonnegative.

Proposition 1

Let \( {\underline{\omega }}=\underset{i=1,\ldots ,10}{\min }\left\{ \omega _{i}\right\} \), \( {\overline{\pi }}=\underset{i=1,\ldots ,10}{\max }\left\{ \pi _{i}\right\} \), \( C_{D}=\frac{\Lambda _{8}+\mu \kappa }{\rho _{1}} \), \( C_{D}=\frac{\Lambda _{9}+\left( \mu +d\right) \kappa }{\rho _{2}} \), \( C_{V}=\frac{\kappa {\overline{\pi }}}{\nu +{\underline{\omega }}} \) and \( N_{V}=\sum _{i=1}^{10} Y_{i}\). \( \Omega =\left\{ 0<N\le \kappa ,X_{8}\le C_{D}, X_{9}\le C_{D_{I}},N_{V}\le C_{V}\right\} \cap {\mathbb {R}}_{+}^{20} \) is an attractor of the invariant positive orthant \({\mathbb {R}}_{+}^{20}\) with respect to the system (22).

Proposition 1 shows the validity of the model ( 1)–(20). It shows the boundedness of the population density (\(\kappa \)) due to the social distancing. It also allows us to consider separately the microscopic dynamics (of viruses) and the macroscopic dynamics (of humans). Indeed, control strategies at the microscopic level concern the individual precautions to be applied (social distancing, mobility, disinfection of people and surfaces, mask wearing and treatment). At the macroscopic level, the global health policies are considered (social distancing, wearing mask, screening, rapid burying of infectious dead, isolation and treatment of infectious). Of course there are measures that act at both levels.

2.3 Equilibriums and asymptotic behaviors

Omitting the compartment of buried (B or \( X_{10} \)), we can address the issue of existence of equilibriums. We are particularly interested in the ’Disease Free Equilibrium’ (\(DFE= \left( X^{DFE},Y^{DFE} \right) \)) and other feasible ’Endemic Equilibrium’ (\( EE=\left( X^{EE},Y^{EE} \right) \)). It is clear that immigration in the infectious compartments with a constant positive rate makes impossible the existence of the DFE. It should be noted, however, that if horizontal immigration is proportional to the population, then a disease-free balance is still possible.

If the geographical area of study extends to the planet ( the model and the strategy of disease control are on a planet level) or a systematic control is applied at all the territorial frontiers (spatial isolation or lock-down). Of course, the latter possibility is less realistic than the former. Indeed, the concept of horizontal immigration is defined when there is a transition of material from one delimited geographical area to another delimited area. Therefore, if we consider the global population we will have only one area (global) and no horizontal immigration (ie \( \Lambda _{i}=0, i=2,\ldots ,6,8,9 \)).

If \( \Lambda _{i}=0, i=1,\ldots ,9 \) then there is no horizontal recruitment and the only possibility of recruitment of new individuals is only from new births. That last consideration will be adopted here for simplicity.

Proposition 2

Assume that \( \Lambda _{i}=0,\ i=1,\ldots ,9 \). Then the following statements hold.

  1. (i)

    The trivial null equilibrium always exists : \( X=Y=0 \),

  2. (ii)

    If \( b_{1}>\mu \), then the non-trivial disease-free equilibrium (DFE) exists :

    $$\begin{aligned}&X_{1}^{DFE}=N^{DFE}=\kappa \left( 1- \frac{\mu }{b_{1}}\right) ,\ X_{2}^{DFE}=\ldots ,\\&X_{7}^{DFE}{=}X_{9}^{DFE}{=}0, X_{8}^{DFE}{=}\frac{\mu X_{1}^{DFE}}{\rho _{1}},\ Y^{DFE}{=}0, \end{aligned}$$
  3. (iii)

    If \( b_{1}>\mu \) and there is an endemic equilibrium (EE) then \( c_{7}\ell +\sum _{i=1}^{7}c_{i}b_{i}>\mu +\uplambda >\mu \) where

    $$\begin{aligned}&\uplambda =\frac{a_{1}\tau _{4}Y_{1}^{EE} }{a_{2}+Y_{1}^{EE}},\\&X_{1}^{EE}=-\frac{\kappa \left( \mu +\uplambda -c_{7}\ell -\sum _{i=1}^{7}c_{i}b_{i} \right) }{\left( \sum _{i=1}^{7}c_{i} \right) \sum _{i=1}^{7}c_{i}b_{i} },\\&X_{i}^{EE}=c_{i}X_{1}^{EE} \text { and } \end{aligned}$$

    the \( c_{i} \)s are given as follows:

    $$\begin{aligned}&c_{1}=1,\ c_{2}=-\frac{\uplambda }{m_{2,2}},\ c_{3}=-\frac{m_{3,2}c_{2}}{m_{3,3}},\ \\&c_{4}=-\frac{m_{4,3}c_{3}}{m_{4,4}}, \quad c_{5}=-\frac{ m_{5,3}c_{3}+m_{5,4}c_{4}}{m_{5,5}}, \\&c_{6}=-\frac{ m_{6,3}c_{3}+m_{6,4}c_{4}+m_{6,5}c_{5}}{m_{6,6}} \text { and }\\&c_{7}=-\frac{ m_{7,6}c_{6}}{m_{7,7}}. \end{aligned}$$

Let us consider the Jacobian matrix of (22) given by

$$\begin{aligned} J= \left[ \begin{matrix} J_{1,1} &{} J_{1,2}\\ J_{2,1} &{} J_{2,2} \end{matrix} \right] \end{aligned}$$

with

$$\begin{aligned}&J_{1,1}\left( 1,j \right) =M_{1,j}-\sum _{i=1}^{7}\frac{b_{i}X_{i}}{\kappa },\ j=1,\ldots ,7 \text { and }\\&J_{1,1}\left( i,j \right) =M_{i,j} \text { in other cases,} \\&J_{1,2}\left( 1,1 \right) {=}-\frac{a_{1}a_{2}\tau _{4}X_{1}}{\left( a_{2}{+}Y_{1}\right) ^{2}},\ J_{1,2}\left( 2,1 \right) {=}\frac{a_{1}a_{2}\tau _{4}X_{1}}{\left( a_{2}+Y_{1}\right) ^{2}},\\&J_{1,1}\left( i,j \right) =0 \text { in other cases,} \\&J_{2,1}\left( i,i \right) {=}\varsigma _{i}\pi _{i}+\left( \sum _{j=1}^{9}\tau _{3,j}Y_{j}{+}\tau _{2,i}Y_{10}\right) \frac{N-X_{i}}{N^{2}} ,\\&i=1,\ldots ,9, \\&J_{2,1}\left( 10,j \right) =\left( 1-\varsigma _{j}\right) \pi _{j}-\frac{\tau _{2,j}Y_{10}}{N}\\&+Y_{10}\sum _{k=1}^{9} \frac{\tau _{2,k}X_{k} }{N^{2}} ,\ j=1,\ldots ,9, \\&J_{2,1}\left( i,j \right) =0 \text { in other cases,} \\&J_{2,2}\left( 1,1 \right) =Q_{1,1}+\frac{\tau _{3,1}X_{1}}{N}-\frac{a_{1}a_{2}\tau _{4}Y_{1}}{\left( a_{2}+Y_{1} \right) ^{2} },\\&J_{2,2}\left( i,i \right) =Q_{i,i}+\frac{\tau _{3,i}X_{i}}{N},\ i=2,\ldots ,9 \\&J_{2,2}\left( 2,1 \right) =\frac{\tau _{3,1}X_{2}}{N}+\frac{a_{1}\tau _{4}Y_{1}\left( 2a_{2}+Y_{1} \right) }{\left( a_{2}+Y_{1} \right) ^{2} },\\&J_{2,2}\left( i,j \right) =Q_{i,j} +\frac{\tau _{3,j}X_{i}}{N},\ i,j=1,\ldots ,9, \\&J_{2,2}\left( i,10 \right) =\frac{\tau _{2,i}X_{i}}{N},\ i=1,\ldots ,9, \text { and }\\&J_{2,2}\left( 10,j \right) =Q_{10,j},\ j=1,\ldots ,10. \end{aligned}$$

Proposition 3

Assume that \( \Lambda _{i}=0, i=1,\ldots ,10 \).

  1. (i)

    The null invariant set \(\left\{ 0\right\} ^{9}\times {\mathbb {R}}_{+}\times \left\{ 0\right\} ^{10}\) is locally asymptotically stable (LAS) if \(b_{1}\le \mu \).

  2. (ii)

    The null equilibrium \( X=Y=0 \) is globally asymptotically stable (GAS) if \(b_{i}\le \mu ,\ i{=}1,\ldots ,7\).

  3. (iii)

    The disease-free equilibrium is locally asymptotically stable (LAS) if and only if \(b_{1}>\mu \) and \( J_{1,1}-J_{1,2}J_{2,2}^{-1}J_{2,1} \) is Metzler stable.

  4. (iv)

    If \(b_{1}>\mu \) and \( \uplambda =0 \) then the disease-free equilibrium for the macroscopic model (\( X=X^{DFE}, \ Y\in \left( {\mathbb {R}}_{+}\right) ^{10}\)) is GAS.

Let \( {\mathcal {F}} \) denote the vector of flows due to infectious contacts between humans an viruses in the disease associated compartments (E, I, Q, H, T, \( D_{I} \), \( V_{S} \), \( V_{E} \), \( V_{I} \), \( V_{Q} \), \( V_{H} \), \( V_{T} \), \( V_{R} \), \( V_{D} \), \( V_{D,I} \) and \( V_{F} \)). Similarly, let \( {\mathcal {V}} \) denote the vector of flows in the disease associated compartments but which are not due to infectious contacts:

$$\begin{aligned}&{\mathcal {F}}_{1}=\frac{a_{1}\tau _{4}X_{1}Y_{1}}{a_{2}+Y_{1}},\ {\mathcal {F}}_{i}=0,\ i=2,\ldots , 7,\nonumber \\&{\mathcal {F}}_{8}=\frac{a_{1}\tau _{4}Y_{1}^{2}}{a_{2}+Y_{1}},\ {\mathcal {F}}_{i}=0,\ i=10,\ldots ,15,\nonumber \\&{\mathcal {V}}={\mathcal {G}}\left( X,Y\right) -{\mathcal {F}} \end{aligned}$$
(23)

with \( {\mathcal {G}} \) given in 22).

If F and V are the respective Jacobian matrix of \( {\mathcal {F}} \) and \( {\mathcal {V}} \) evaluated at the DFE, then the basic reproduction number \({\mathcal {R}}_{0}\) is given by the spectral radius of the next generation matrix (Van Den Driessche and Watmough method [78,79,80]). In this context, F is a \( 16\times 16 \) matrix given such as \( F\left( 1,1\right) =\frac{a_{1}\tau _{4}}{a_{2}} \), \( F\left( 1,7\right) =\frac{a_{1}\tau _{4}X_{1}^{DFE}}{a_{2}} \) and \( F\left( i,j\right) =0 \) otherwise. V is a \( 16\times 16 \) matrix given such as \( V_{1,7}=0 \), \( V_{i,j}=J\left( i+1,j+1\right) \) if \( i,j=1,\ldots ,5 \), \( V_{6,j}=J\left( 9,j+1\right) \) if \( j=1,\ldots ,5 \), \( V_{6,6}=J\left( 9,9\right) \), \( V_{6,j}=J\left( 9,j+4\right) \) if \( j=7,\ldots ,16 \), \( V_{i,6}=J\left( i+1,9\right) \) if \( i=1,\ldots ,5 \), \( V_{i,6}=J\left( i+4,9\right) \) if \( i=7,\ldots ,16 \), \( V_{i,j}=J\left( i+4,j+4\right) \) if \( i,j=7,\ldots ,16 \), \( V_{i,j}=J\left( i+1,j+4\right) \) if \( i=1,\ldots ,5,\ j=7,\ldots ,16,\ \left( i,j \right) \ne \left( 1,7\right) \) and \( V_{i,j}=J\left( i+4,j+1\right) \) if \( i=7,\ldots ,16,\ j=1,\ldots ,5\). Although the calculation of the inverse of V is nontrivial, we have the particular formula

$$\begin{aligned} {\mathcal {R}}_{0}=\rho \left( -FV^{-1}\right) =-F\left( 1,7\right) \times V^{-1}\left( 7,1 \right) . \end{aligned}$$
(24)

The effective reproduction number can also be computed by multiplying \({\mathcal {R}}_{0}\) by the proportion of susceptible individuals [70]. \({\mathcal {R}}_{0}\) can be biologically interpreted in this context as the average number of new viruses generated by an introduced virus during its lifespan.

The infection force \( \uplambda \) plays an important role on the disease spreading. Subsequently, regarding the expression of \( \uplambda \) in Table 1, we can see that the use of mask and avoiding to touch the face is crucial to limit disease spreading. The control u can be viewed as the fraction of time the mask is well used. \( \uplambda \) can vary with time depending on individual human habits. It is interesting to depict how the disease behaves if \( \uplambda \) is a decreasing function of time (\( \underset{t\rightarrow +\infty }{\lim }\uplambda \left( t\right) =0\)). We start by stating a useful result which is a particular case of very useful theorems in [36, 58, 76], about \( \omega \)-limit sets of asymptotic autonomous systems.

Lemma 1

Let us consider a linear differential system

$$\begin{aligned}&\overset{\bullet }{U}=AU+C \end{aligned}$$
(25)
$$\begin{aligned}&\overset{\bullet }{V}=\left( A+v\left( t\right) B\right) V+C \end{aligned}$$
(26)

where A is a constant \( n \times n \) Metzler stable matrix, B is a constant \( n \times n \) Metzler matrix, C is a constant \( n \times 1 \) matrix and v is a positive numeric function of the time. Assume that the solution \( \left( U,V \right) \) of (25)–( 26) is bounded and . Then \( \underset{t\rightarrow +\infty }{\lim }\left( U\left( t\right) -V\left( t\right) \right) =0 \). This will happen in particular if \( \underset{t\rightarrow +\infty }{\lim }v\left( t\right) e^{\rho \left( A\right) t}=0 \), where \( \rho \left( A\right) \) denotes the spectral radius of A.

Instead of using the technical approaches in [36, 58, 76], Lemma 1 is proved here in a simple contextual manner in the appendix. The remaining results of the paper are original and so are their proofs. Although the latter combine fairly general results on dynamical systems, their implementation is non-trivial according to the structure of the model ( 1)–(20).

Theorem 1

If \( b_{1}>\mu \) and \( \underset{t\rightarrow +\infty }{\lim }\uplambda \left( t\right) \left\| e^{\rho \left( M_{1,1}|_{\uplambda =0} \right) t }\right\| =0\) then the DFE is GAS.

Theorem 1 gives us a first way to limit disease spreading.

Theorem 2

Assume that \( \Lambda _{i}=0, i=1,\ldots ,10 \) and \(b_{1}>\mu \). Let D denotes a diagonal matrix such that \( D_{i,i} \) is the sum of absolute values of only negative terms in the expression of \( Q_{i,i} \), \( R=D+Q \), \( r_{1}=\underset{j=1,\ldots ,10}{\max } \frac{R_{1,j}}{D_{j,j}} \) and \( r_{2}=\underset{j=1,\ldots ,10}{\max } \sum _{i=1}^{10}\left( RD^{-1}\right) _{i,j}\). Then the DFE is GAS if

$$\begin{aligned} {\mathcal {T}}_{0}=\frac{a_{1}\alpha _{1}r_{1}\tau _{4}\kappa \left( b_{1}-\mu \right) \left( \frac{\pi _{3}}{\mu +d+q+\gamma _{1}}+\frac{\pi _{4}}{\mu +d+\alpha _{2}+\gamma _{2}}+\frac{\pi _{5}}{\mu +d+\gamma _{3}}+\frac{\pi _{6}}{\mu +\gamma _{4}}\right) }{a_{2}b_{1}\left( 1-r_{2}\right) \left( \mu +\alpha _{1} \right) \left( \mu +\nu +\omega _{1}+\tau _{1,1}+\tau _{3,1}\right) }<1. \end{aligned}$$
(27)

Theorem 2 provides a sufficient condition that ensures the stability of the disease free equilibrium. Moreover, it shows that avoiding to touch the face and the frequency of wearing mask by susceptible population (\( u_{1} \) via \( \tau _{4} \)), the mobility (m via different contact rate among people \( \tau _{3} \), and with the environment \( \tau _{1} \)), the disinfection (washing of hands included in) by susceptible population (\( \omega _{1} \)) and social distancing (reducing \( \kappa \)) are crucial parameters in order to stop the spreading of the disease. The biological interpretation of \( {\mathcal {T}}_{0} \) is the following. During its lifespan \( \frac{1}{\left( \mu +\nu +\omega _{1}+\tau _{1,1}+\tau _{3,1}\right) } \) a virus infects a susceptible individual at the adequate rate \( \frac{a_{1}r_{1}\kappa \tau _{4}\left( b_{1}-\mu \right) }{a_{2}b_{1}\left( 1-r_{2}\right) } \). Once the susceptible individual is infected he becomes infectious with a probability \( \frac{\alpha _{1}}{\mu +\alpha _{1}} \) and will produce \( \frac{\pi _{3}}{\mu +d+q+\gamma _{1}}\), \(\frac{\pi _{4}}{\mu +d+\alpha _{2}+\gamma _{2}}\), \( \frac{\pi _{5}}{\mu +d+\gamma _{3}}\) and \(\frac{\pi _{6}}{\mu +\gamma _{4}} \) viruses during its respective stays in compartments of infectious, quarantined, hospitalized and treated.

By setting

$$\begin{aligned}&\delta {=}\max \left\{ 1{+}\frac{\left( \mu {+}\tau _{1,1}+\tau _{3,1}\right) \left( a_{2}\left( \nu +{\underline{\omega }} \right) +\kappa {\overline{\pi }}\right) +a_{1}\kappa \tau _{4}{\overline{\pi }}}{\left( \nu +\omega _{1} \right) \left( a_{2}\left( \nu +{\underline{\omega }} \right) +\kappa {\overline{\pi }}\right) },\right. \\&\left. 1+\underset{j=1,\ldots ,10}{\max }\frac{\tau _{2,j}}{\nu +\omega _{10}},\ \underset{j=2,\ldots ,9}{\max }\frac{D_{j,j}}{\nu +\omega _{j}} \right\} \text { and } \\&{\overline{\delta }}=\max \left\{ 1+\frac{\left( \mu +\tau _{1,1}+\tau _{3,1}\right) \left( a_{2}\nu +\kappa {\overline{\pi }}\right) +a_{1}\kappa \tau _{4}{\overline{\pi }}}{\nu \left( a_{2}\nu +\kappa {\overline{\pi }}\right) },\right. \\&\left. 1+\underset{j=1,\ldots ,10}{\max }\frac{\tau _{2,j}}{\nu },\ \underset{j=2,\ldots ,9}{\max }\frac{D_{j,j}}{\nu } \right\} , \end{aligned}$$

we can see using the boundedness of \( N_{V} \), that \( \frac{1}{1-r_{2}}\le \delta \le {\overline{\delta }} \). Hence, Theorem 2 also permits us to evaluate given \( \kappa \), u and m the minimum frequency at which the disinfection should occur. For example, if \( {\underline{\omega }}=0 \) probably due to the management of death persons, then the condition \( {\mathcal {T}}_{0}\le 1 \) is satisfied when

$$\begin{aligned}&\omega _{1}>\frac{a_{1}\alpha _{1}\tau _{4}\kappa {\overline{\delta }}\left( b_{1}-\mu \right) \left( \frac{\pi _{3}}{\mu +d+q+\gamma _{1}}+\frac{\pi _{4}}{\mu +d+\alpha _{2}+\gamma _{2}}+\frac{\pi _{5}}{\mu +d+\gamma _{3}}+\frac{\pi _{6}}{\mu +\gamma _{4}}\right) }{a_{2}b_{1}\left( \mu +\alpha _{1} \right) }\nonumber \\&\quad -\left( \mu +\nu +\tau _{1,1}+\tau _{3,1}\right) . \end{aligned}$$
(28)

Notice that by disinfection we mean a global action on the body, so it would be necessary to evaluate the specific contribution of hands washing.

Fig. 2
figure 2

Asymptotic behavior of infectious dynamics depending on \( {\mathcal {T}}_{0} \)

Figure 2 is an illustration of the asymptotic behavior of cumulative infectious individuals depending on the relative position of \( {\mathcal {T}}_{0} \) with respect to 1. The smaller \( {\mathcal {T}}_{0} \) is, the smaller is the amplitude of the peak and the earlier it occurs. If the mask is never adopted then the disease prevalence can pass above \( 30\% \) and stabilize around \( 3.5\% \).

We end the section with a theorem giving an upper bound to the time of extinction of the disease. By the time of extinction we mean the first time where there is neither infectious and nor virus. That is there is less than one infectious and less than one virus.

Theorem 3

Let \( {\mathcal {S}}\,{=}\,\frac{S}{S+E+I+Q+H+T+R} \), \({\mathcal {I}}\,{=}\frac{E+I+Q+H+T}{S+E+I+Q+H+T+R} \) and \( {\mathcal {A}}\) denotes the effective area of the spatial study zone in \( m^{2} \) (ie the area of the zone each individual is allowed to move in during the lock-down). Assume that \( \forall t\ge 0 \), \( \uplambda \left( t\right) =\uplambda _{0}e^{-\alpha t},\ \alpha >0 \). Then

$$\begin{aligned} \overset{\bullet }{{\mathcal {I}}}= & {} \uplambda {\mathcal {S}}+\frac{bN}{\kappa }{\mathcal {I}}-\left( b+\beta \left( 1-p_{d}\right) \right) {\mathcal {I}}\nonumber \\&-\left( \beta p_{d}-\mu \right) \left( 1-{\mathcal {I}}\right) {\mathcal {I}}, \end{aligned}$$
(29)

where \( \beta =\left( \sum _{i=0}^{4}T_{i} \right) ^{-1} \), \( p_{d} \) is the probability to die during the stay in infected or infectious compartments (E, I, Q, H, T) and the \( T_{i} \)s denote the average times spent in those compartments. Moreover, if \( \sigma _{1}=\min \left\{ \alpha ,\beta \left( 1-p_{d}\right) \right\} \) then

$$\begin{aligned} {\mathcal {I}}\left( t \right) \le {\mathcal {I}}\left( 0 \right) e^{-\left( \beta \left( 1-p_{d}\right) \right) t}+\uplambda _{0}te^{-\sigma _{1} t}. \end{aligned}$$
(30)

Furthermore, if \( \sigma _{2}=\min \left\{ \nu +{\underline{\omega }}\ ,\beta \left( 1-p_{d}\right) \right\} \), \( \sigma _{3}=\min \left\{ \nu +{\underline{\omega }}\ ,\sigma _{1}\right\} \) then

$$\begin{aligned}&N_{V}\left( t\right) \le N_{V}\left( 0\right) e^{-\left( \nu +{\underline{\omega }}\right) t} +\kappa {\overline{\pi }}{\mathcal {I}}\left( 0 \right) te^{-\sigma _{2}t}\nonumber \\&\qquad \quad +\kappa {\overline{\pi }}\uplambda _{0}t^{2}e^{-\sigma _{3}t} \end{aligned}$$
(31)

and the respective first times \( T^{\star }_{H}\) and \( T^{\star }_{V}\) such that, respectively, \( {\mathcal {A}}{\mathcal {I}}_{G}\left( t \right) {\overline{N}}\left( t \right) \le 1 \) and \( {\mathcal {A}}N_{V}\left( t\right) \le 1 \) satisfy

$$\begin{aligned} T^{\star }_{H} \le&\frac{\kappa {\mathcal {A}}{\mathcal {I}}\left( 0 \right) }{\beta \left( 1-p_{d}\right) e} +\frac{4\kappa \uplambda _{0}{\mathcal {A}}}{\sigma _{1}^{2}e^{2}} \text { and } \\ T^{\star }_{V} \le&\frac{{\mathcal {A}}N_{V}\left( 0\right) }{\left( \nu +{\underline{\omega }}\right) e}+\frac{\kappa {\overline{\pi }}{\mathcal {A}}\left( 8e\sigma _{3}^{3}{\mathcal {I}}\left( 0 \right) +27\sigma _{2}^{2}\uplambda _{0} \right) }{2\sigma _{2}^{2}\sigma _{3}^{3}e^{3}}. \end{aligned}$$

As we can see in the expressions of \( T^{\star }_{H} \) and \( T^{\star }_{V} \) that the higher \( \kappa \), \( \uplambda _{0} \), \( {\mathcal {I}}\left( 0\right) \), \( N_{V}\left( 0\right) \) and \( {\mathcal {A}}\) are, the longer the time of disease extinction \( \max \left\{ T^{\star }_{H}, T^{\star }_{V} \right\} \) is. If the lock-down is applied, then the effective global population to consider is \( {\overline{N}}_{G}={\mathcal {A}}{\overline{N}}\) which represents the number of living persons admitted in a given closed area. Hence, the lock-down level can be expressed by the ratio \( 1-\frac{{\mathcal {A}}}{{\mathcal {A}}_{G}} \), where \( {\mathcal {A}}_{G} \) denotes the real area of the physical study zone (town, district, department, region, country, continent, or the whole world).

3 A focus on Cameroonian context

3.1 Early estimation of the parameters of the disease dynamics

This section aims to determine values of the parameters of the model ( 1)–(20) in such a way that its predicted dynamics is close to the Corona’s observed data in Cameroon, from March 2 to June 7, 2020. That is estimation of parameters also known in the literature as fitting the model. We used the daily data available on the website https://www.data.gouv.fr/fr/datasets/coronavirus-covid19-evolution-par-pays-et-dans-le-monde-maj-quotidienne/www.data.gouv.fr for our estimates. We also used the situation reports on Covid-19 in Cameroon available on the website https://www.humanitarianresponse.info/en/operations/cameroon/health/documentshttps://www.humanitarianresponse.info. It corresponds to the beginning of the epidemic, when the natural dynamics of the disease were less complex to predict since control measures were simpler and secondary approaches (e.g. periodic containment and vaccine) had not yet taken place. The data consist of daily counts of infections, recoveries, deaths, quarantines and hospitalizations. Following the reference [52, 53], we set \( CRI\left( t\right) \), \( CRQ\left( t\right) \), \( {\mathrm{CRH}}\left( t\right) \), \( CRR\left( t\right) \) and \( CRD\left( t\right) \) to be, respectively, the cumulative numbers of reported infectious cases, of reported quarantined, of reported hospitalized (active cases), of reported recovery and of reported dead at time \( t\ge t_{0} \), where \( t_{0} \) denotes the date of the first introduction of an infectious to be evaluated. Note that these values will be divided by the cumulative area of populated areas (\( {\mathcal {A}}_{G} \)) to match the model ( 1–(20) for which variables are expressed in density per unit area. There are several approaches for the estimation of the parameters of a model. A first method is to discretize the differential equation model with respect to available time series and to fit unknown parameters by reducing the distance between the numerical solution and observed time series [13]. Another method consists of empirical estimating function linking the state of the system with the time using times series and after in estimation through the considered model [12, 24, 54, 65].

Here we use a semi-empirical method using statistical regressions based on empirical models according to data (the upcoming model ( 40) for example) and key relations derived from the main model (1)–(20). In the global population (as defined in (21)), those relations are the following :

$$\begin{aligned}&CRI\left( t\right) =q\int _{t_{0}}^{t}I_{G}\left( s\right) ds , \end{aligned}$$
(32)
$$\begin{aligned}&\int _{t_{0}}^{t}Q_{G}\left( s\right) ds=CRQ\left( t\right) =pq \int _{t_{0}}^{t} I_{G}\left( s\right) ds , \end{aligned}$$
(33)
$$\begin{aligned}&\int _{t_{0}}^{t}H_{G}\left( s\right) ds=CRH\left( t\right) \nonumber \\&=\int _{t_{0}}^{t}\left( \left( 1-p \right) q I_{G}\left( s\right) +\alpha _{2}Q_{G}\left( s\right) \right) ds , \end{aligned}$$
(34)
$$\begin{aligned}&\int _{t_{0}}^{t}R_{G}\left( s\right) ds=CRR\left( t\right) \nonumber \\&=\int _{t_{0}}^{t}\left( \gamma _{2}Q_{G}\left( s\right) +\gamma _{3}H_{G}\left( s\right) \right) ds , \text { and } \end{aligned}$$
(35)
$$\begin{aligned}&\int _{t_{0}}^{t}D_{G}\left( s\right) ds=CRD\left( t\right) =\left( \mu +d \right) \nonumber \\&\int _{t_{0}}^{t}\left( Q_{G}\left( s\right) +H_{G}\left( s\right) \right) ds . \end{aligned}$$
(36)

According to certain cultural practices in CameroonFootnote 6 (See [40]), we will assume that the average time that elapses between the death of an individual not identified as infectious (i.e., he or she dies without being quarantined or hospitalized) and his or her burial is 3 days (\( \rho _{1}=\frac{1}{3}\)), while burial takes place within 24 hours of the death of an individual considered infectious (\( \rho _{2}=1\)). According to the literature, we will assume that the transition rate from the infected compartment (E) to the infectious compartment (I) is \( \alpha _{1}=\frac{1}{14} \), the transition rates from the infectious compartments (I,Q,H) to the treated compartment (T) are \( \gamma _{i}=\frac{1}{7}, \ i=1,\ldots ,3 \) and the transition rate from the treated compartment (T) to the recovered/immune compartment (R) is \( \gamma _{4}=\frac{1}{7}\) [52, 53]. It is not clearly established in the literature that there is a real immunity, but according to [42], we will suppose that \( \ell =2.74\times 10^{-3} \). The recommended social distancing parameter \( \kappa \) in Cameroon is \( 1 {\mathrm{person}}\,\times \,{\mathrm{m}}^{-2}\) but we will estimate the real value of \( \kappa \) according to available data and the model (1)–( 20).

The data on CRQ, CRQ being available in situation reports only for the 30 firsts days contrarily to the data on CRI, CRD and CRR which cover the whole period of study. Hence we divide the estimation in two steps :

  • Step 1: Use of data from March 2 to April 1, 2020. Summing equations (32), (33) and (34), we have

    $$\begin{aligned} \alpha _{2}\approx 1+{\mathbb {E}}\left[ \frac{CRH\left( t\right) -{\mathrm{CRI}}\left( t\right) }{{\mathrm{CRQ}}\left( t\right) }\right] , \end{aligned}$$
    (37)

    \( t_{0} \) being the greater time where \(I=Q=H=T=R=0\). The mathematical expectation \( {\mathbb {E}}\left[ \bullet \right] \) approximated by the empirical average is used here to contain uncertainties causing variation on \( \frac{{\mathrm{CRH}}\left( t\right) -{\mathrm{CRI}}\left( t\right) }{{\mathrm{CRQ}}\left( t\right) } \). Such as statistical approach could also permit to evaluate a confidence interval assuming adequate hypotheses on the distribution of the noises. Commonly, large numbers law and central limit theorem are used for large dataset [48]. Due to their relatively rapid growth, we suppose as in the references [52, 53] that \(CRD\left( t\right) \) and \({\mathrm{CRH}}\left( t\right) \) can be approximated through log-linear models as follows

    $$\begin{aligned}&{\mathrm{CRQ}}\left( t\right) \approx e^{\left( k_{Q}\left( t-t_{1}\right) \right) }-1, \end{aligned}$$
    (38)
    $$\begin{aligned}&{\mathrm{CRH}}\left( t\right) \approx e^{\left( k_{H}\left( t-t_{2}\right) \right) }-1, \end{aligned}$$
    (39)
  • Step 2: Use of data from March 2 to June 7, 2020. Regarding the rapid growth of reported infectious cases, we assume that CRI can be approximated through the following empirical model

    $$\begin{aligned} {\mathrm{CRI}}\left( t\right) \approx e^{\left( k_{I}\left( t-t_{0}\right) \right) }-1, \end{aligned}$$
    (40)

    According to the equations (32)–(36) we have

    $$\begin{aligned}&T_{G}\left( t\right) =e^{-\left( \mu +\gamma _{4} \right) t }\nonumber \\&\int _{t_{0}}^{t}e^{\left( \mu {+}\gamma _{4} \right) s}\left( \gamma _{1}I_{G}\left( s\right) {+}\gamma _{2}Q_{G}\left( s\right) {+}\gamma _{3}H_{G}\left( s\right) \right) \nonumber \\&{\mathrm{d}}s \text { and} \end{aligned}$$
    (41)
    $$\begin{aligned}&R_{G}\left( t\right) =\gamma _{4}e^{-\left( \mu +\ell \right) t }\int _{t_{0}}^{t}e^{\left( \mu +\ell \right) s}T_{G}\left( s\right) ds. \end{aligned}$$
    (42)

    The following approximations hold

    $$\begin{aligned}&p\approx {\mathbb {E}}\left[ \frac{{\mathrm{CRQ}}\left( t\right) }{{\mathrm{CRI}}\left( t\right) }\right] , \end{aligned}$$
    (43)
    $$\begin{aligned}&q=\frac{1}{p\left( t_{1}-t_{0}\right) }, \end{aligned}$$
    (44)
    $$\begin{aligned}&d\approx {\mathbb {E}}\left[ \frac{{\mathrm{CRD}}\left( t\right) }{\int _{t_{0}}^{t}\left( Q_{G}\left( s\right) +H_{G}\left( s\right) \right) ds}\right] -\mu \text { and} \end{aligned}$$
    (45)
    $$\begin{aligned}&E_{G}\left( t\right) =\frac{{\mathrm{d}}}{{\mathrm{d}}t}\left[ \frac{qI_{G}\left( t\right) +\left( \mu +d+q+\gamma _{1}\right) {\mathrm{CRI}}\left( t\right) }{q\alpha _{1}}\right] .\nonumber \\ \end{aligned}$$
    (46)

The coefficients \( k_{I}=0.1717935\), \( k_{Q}=0.202587\), \( k_{H}=0.128817\), \( t_{0}=-2.646517\), \( t_{1}=1.484067\), \( t_{2}=-5.746704\) were estimated through log-linear regression with the software R. Among the above coefficient only \( t_{2}\) was found not statistically different from zero during the regression analysis of variance. The adopted empirical models were highly significant with all the adjusted \( R^{2}\)s between \(82.98\% \) and \( 96.03\% \), and the probability-values (commonly called p-values) between \( 8.096\times 10^{-8} \) and \( 2.716\times 10^{-3} \). A negative \( t_{0} \) means that the first infectious case occurs probably before the origin date of the study. The same holds for \( t_{1} \) and \( t_{2} \). Since CRI was given by (32) it was easy to determine the expression of \( E_{G}\left( t\right) \), in particular \( E_{G}\left( 0\right) \). The statistical approach we adopted is justified by the uncertainties. \( \alpha _{2} \), p, q, d and other estimated values are available in Table 2.

It appeared more difficult to explicitly estimate \( a_{1} \), \( a_{2} \), \( \pi _{i} \), \( \tau _{i,j}^{\min } \), \( \tau _{i,j}^{\max } \), \( \tau _{4}^{\min } \), \( \tau _{4}^{\max } \), \( \varsigma _{i}^{\min } \), \( \varsigma _{i}^{\max } \), \( m_{i} \) and \( u_{i} \) using available data. Regarding public transport conditions in Cameroon the maximum value of \( \kappa \) can be fixed to \( {\overline{\kappa }}=8\,{\mathrm{person}}\,\times \,{\mathrm{m}}^{-2} \). According to the density of population per unit of surfaceFootnote 7 the minimum value of \(\kappa \) is fixed to \( {\underline{\kappa }}=5.003\times 10^{-5}\ {-2} \). Without social distancing at the disease outbreak, we have \( \kappa ={\overline{\kappa }} \) and the infection term \( \uplambda S \) is proportional to the global population (say \( S_{G} \)) due to the mass action. S being bounded from above by \( {\overline{\kappa }} \) according to the model, we have \( a_{1}\left( t\right) =\frac{S_{G}\left( t\right) }{{\overline{\kappa }}\left( t\right) }\). Notice that \( {\underline{\kappa }}\) and \( {\overline{\kappa }} \) are non decreasing functions of \( S_{G} \). Thus, we can simplify by assuming that \( a_{1} \) is definitely a constant given by \( \frac{S_{G}\left( 0 \right) }{{\overline{\kappa }}\left( 0 \right) } \). It is possible to estimate \( S_{G} \) and \( N_{G} \) (global population of alive or dead except buried individuals) similarly to \( T_{G} \) and \( R_{G} \). Under the hypothesis of homogeneity of the population the values of S, E and I can be estimated by multiplying the respective proportions in the global population by the threshold \( \frac{\kappa \left( b_{1}-\mu \right) }{b_{1}} \) (for example \( S\left( 0 \right) = \frac{\kappa \left( b_{1}-\mu \right) S_{G}\left( 0 \right) }{b_{1}N_{G}\left( 0 \right) } \)). The term \( \uplambda S=\frac{a_{1}\tau _{4}Y_{1}S}{a_{2}+Y_{1}} \) can be estimated using the relations

$$\begin{aligned} \uplambda S&=\overset{\bullet }{E}+\left( \mu +\alpha _{1}\right) E \nonumber \\&=\frac{\kappa \left( b_{1}-\mu \right) }{b_{1}}\left( \frac{\overset{\bullet }{E}_{G}}{N_{G}}-\frac{E_{G}}{N_{G}^{2}}\overset{\bullet }{N}_{G}+\left( \mu +\alpha _{1}\right) \frac{E_{G}}{N_{G}}\right) \nonumber \\&\text { and } \end{aligned}$$
(47)
$$\begin{aligned}&Y_{1}=\frac{a_{2}\uplambda S}{a_{1}\tau _{4}-\uplambda S}. \end{aligned}$$
(48)
Fig. 3
figure 3

Cumulative reported cases from the 2nd of March to the 07th of June 2020

Figure 3 shows the cumulative reported infectious cases in Cameroon from the 2nd of March to the 7th of June 2020. When studying the dynamics of cumulative reported infectious cases we noticed a change approximatively around the thirtieth days. Indeed, the growth of the cumulative number of reported infectious appeared to be higher during the first thirty days. This is likely due to the delay in the implementation of the barrier measures (which were announced by the Cameroonian government only as of March 20, 2020) and the delay in their effects on the disease dynamics.Footnote 8 Thus, to capture this apparent change, we split the empirical modeling into two periods: those of the first 30 days and those of the days after. This splitting also allowed for a better fit of the model. Parameters concerning the disease outbreak (starting time, initial values of state variables) have been estimated using the first dataset. The other parameters have been fitted using the whole sample. The estimates are given in Table 2.

According to the values in Table 2, relation (28 ) and Theorem 2, if people never wear a mask, no social distancing is applied and only susceptible individuals disinfect themselves, then the rate of complete disinfection (\( \omega _{1} \)) should not be less than \( 20\ {\mathrm{day}}^{-1} \) (ie about once every 72 minutes) in order to have \( {\mathcal {T}}_{0}=7.3163071\) and \( {\mathcal {R}}_{0}=0.887382\). If at least the social distancing of \( 1\,{\mathrm{person}}\,\times \,{\mathrm{m}}^{-2} \) is respected then the lower bound value of \( \omega _{1} \) drops to \( 3 \ {\mathrm{day}}^{-1}\) (ie about once every 8 hours) in order to have \( {\mathcal {T}}_{0}=5.7754491\) and \( {\mathcal {R}}_{0}=0.7064687\). Thus, a complementary application of all mitigation measures is mandatory.

Table 2 Estimation and baseline values of parameters of the model (22)

3.2 Effects of controls on disease dynamics

This section aims at evaluating first the sensitivity of different control parameters with respect to the thresholds \( {\mathcal {R}}_{0} \) and \( {\mathcal {T}}_{0} \). The second fold of this section is to survey the effect of some control parameters on the asymptotic behavior of the disease dynamics. Following Theorem 1, we define the convergence index

$$\begin{aligned} {\mathcal {T}}_{1}=\underset{t\rightarrow +\infty }{\lim }\frac{-1}{t}\int _{0}^{t}\ln \left( \frac{\uplambda \left( s \right) }{\uplambda \left( 0 \right) }\right) {\mathrm{d}}s. \end{aligned}$$
(49)

\( {\mathcal {T}}_{1} \) is an indicator of the rate of the potential convergence of the disease dynamics toward an equilibrium. Indeed, if \( {\mathcal {T}}_{1}>0\) then there is an exponential convergence to the disease free equilibrium. If \( {\mathcal {T}}_{1}\le 0\) then there is an endemicity of the disease, in particular if \( {\mathcal {T}}_{1}\) is near to zero then we are in the neighborhood of an endemic equilibrium. \( {\mathcal {T}}_{1}\) is an estimation of the key parameter \( \alpha \) which permits according to Theorem 3, to bound the average time the disease is expected to disappear.

There are several methods for sensitivity analysis in the literature (see [18, 56] for example). The normalized forward sensitivity index (of variable u with respect to parameter p) \( \Upsilon _{p}^{u}=\frac{\partial u}{\partial p}\frac{p}{u} \) defined in [18] provides an idea on effect of a local perturbations around a given situation in terms of the values of the parameter p. Paper [56] presents several alternative methods of assessing the global effect of variation on coefficients: Pearson correlation coefficient (CC), Spearman rank correlation coefficient (RCC), partial correlation coefficient (PCC), standardized correlation coefficient (SCC), standardized rank correlation coefficient (SRCC), partial rank correlation coefficient (PRCC) and Extended Fourier amplitude sensitivity test (eFAST). The authors in [56] recommended PRCC and eFAST as robust sensitivity indexes to be used, respectively, in cases of monotonic and non-monotonic sensitivities. These indexes are based on samples obtained through latin hyper cube sampling, Monte-Carlo sampling or any other appropriate experimental design [11, 20, 35]. Hence, there are statistics and confidence intervals or region can be provided in order to decide on their significance.

We adopt the PRCC with a three levelsFootnote 9 (\( -1 \), 0, \( +1 \)) fractional experimental design with the factors \( \kappa \), q, \(m_{i}\), \(u_{i}\) and \(\omega _{i}\), \(i=1,\ldots ,10\). Such an experimental design allows us to get a statistical approximation of optimal parameters based on a quadratic model obtained by regression. Ideally, strategies should be specific to different compartments. However, since it is difficult to know the status of people, it is reasonable to apply the same strategies for similar” compartments. Thus, we reduce the factors by considering the classes \( \left[ 1\right] \equiv \left\{ 1,2,3,6,7 \right\} \), \( \left[ 4\right] \equiv \left\{ 4,5 \right\} \), and \( \left[ 8\right] \equiv \left\{ 8,9,10 \right\} \). We also assume that \( u_{8}=u_{9}=u_{10}=0 \). Finally, we have 10 factors with two responses \( {\mathcal {T}}_{0} \) and \( {\mathcal {T}}_{1} \). We consider a \(3^{10-6}\) fractional experimental design with a sample size 81. The experimental matrix is computed using Box calculus for the aliases with second order interactions. Table 3 provides us with the codification of factors values that have been arbitrarily chosen but appear realistic.

Table 3 Codification of factors values
Fig. 4
figure 4

Effects of control parameters on \( R_{0}\), \( T_{0}\) and \( T_{1}\)

Figure 4 shows graphically the effects of control parameters on \( {\mathcal {R}}_{0} \), \( {\mathcal {T}}_{0} \) and \( {\mathcal {T}}_{1} \). \( u_{1} \), \( m_{1} \), \( \omega _{1} \) and \( \omega _{10} \) are in decreasing order the more sensitive parameters with respect to \( {\mathcal {R}}_{0} \). Thus, an efficient control strategy should prioritize the usage of mask, the reduction of mobility and the disinfection in class \( \left[ 1\right] \) compartment, and the disinfection of the environment. \( u_{1} \), \( \omega _{1} \), \( \omega _{4} \) and \( \omega _{10} \) are in decreasing order the more sensitive parameters with respect to \( {\mathcal {T}}_{0} \). Thus, an efficient control strategy based on \( {\mathcal {T}}_{0} \) should prioritize the usage of mask in class \( \left[ 1\right] \), and disinfection in classes \( \left[ 1\right] \), \( \left[ 2\right] \) and the environment. \( m_{1} \), \( u_{1} \), \(\omega _{1} \), and \( \omega _{10} \) are in decreasing order the more sensitive parameters with respect to \( {\mathcal {T}}_{1} \). Thus, to accelerate the convergence toward the disease free equilibrium the fighting policies should prioritize the reduction of mobility in class \( \left[ 1\right] \), the frequent usage of mask, the disinfection of people in class \( \left[ 1\right] \) and the disinfection of the environment.

Globally, the usage of mask in class \( \left[ 1\right] \) compartment is the most sensitive control followed by the reduction of mobility and the disinfecting in class \( \left[ 1\right] \), and the disinfection of the environment. As we can observe, the effects of the social distancing and the screening, the effect of using masks in class \( \left[ 4\right] \), and the mobility in classes \( \left[ 4\right] \) and \( \left[ 8\right] \) are not very significant. However, as we will observe in Sect. 4, those controls have influence on the economical point of view.

4 An optimal control analysis

4.1 Design of the control strategy

There are several works tackling the issue of optimality of fighting strategies against Covid-19 [19, 46, 47, 61, 68, 72]. The objective function to optimize is commonly expressed in the literature by the sum of a positive linear combination of disease states and a positive quadratic combination of control parameters. Such a formulation is a penalization approach to solve the multiple objective problem which is to reduce both the population size of disease compartments and the control energy. The quadratic dependence of the objective functional with respect to control parameters makes easier the optimization procedure since the problem has better chances to be convex [4, 49]. Unfortunately, the choice of penalization coefficients is very often not justified and there is no specific interpretation. Only few papers tried a practical optimization in terms of social and the economical impacts of the pandemic [47, 81]. We follow the references [47, 81] to define a realistic new economical objective function since the main issue of the fight against Covid-19 pandemic seems to eradicate it with limited existing resources. We highlight the contradictory costs, respectively, due to disease spreading and disease containment strategies.

The control parameters we consider here are mainly the limit density of population \( \kappa \), the probability of mobility m, the proportion of time spent on wearing a mask u, the disinfection rate \( \omega \) and the reporting rate of infectious (through screening test especially) q. Indeed, those controls appear explicitly in the expression of the threshold \( {\mathcal {T}}_{0} \) which will play a similar role to the basic reproduction number. Since the threshold \( {\mathcal {T}}_{0} \) appears important for the asymptotic behavior of the disease, it could be interesting to minimize an analog to the effective reproduction number : \( {\mathcal {T}}_{0}\frac{X_{1}}{N} \). Unfortunately, the expression of \( {\mathcal {T}}_{0}\) is very complex and it is not differentiable with respect to all its parameters. This irregularity involves an additional complexity in the optimization procedure. Thus, we will consider the infection rate \( \uplambda X_{1} \) instead of \( {\mathcal {T}}_{0}\frac{X_{1}}{N} \). \( \uplambda X_{1} \) takes into consideration the dynamics of viruses.

The medical monitoring of a patient during his infectious period has an economical cost c that can be evaluated depending on the time spent before recovery and on the expenses during each stay in a disease compartment. The expression of c should depend on \(c_{g}\), the daily revenue of a person in health which will be lost due to the disease. The application of each control involves an economical cost proportional to the size of the population and to the time. We will note the cost per unit of time and per individual \( c_{u} \) for u and \( c_{\omega } \) for \( \omega \). Reducing m and \( \kappa \) involves a relative reduction in economical activities and subsequently a reduction of gross domestic product assumed proportional to \( m\kappa \) (\(GDP\overset{\text {def}}{=} g =1527\,{\mathrm{USD}}\,\times \,{\mathrm{preson}}^{-1}\times \,{\mathrm{year}}^{-1}\) in 2018.Footnote 10) A screening test has a cost \( c_{s} \) which is not clearly defined according to the media. We will assume that \( c_{s}=10\ USD \) in average. Since it is difficult to really estimate unreported infectious, it seems better to define general measures to be applied by everyone independently of the status. Thus, we will assume here that \( \tau _{1,i}=\tau _{1} \), \( \tau _{2,i}=\tau _{2} \), \( \tau _{3,i}=\tau _{3} \), \( u_{i}=u \), \( \omega _{i}=\omega \), \( \varsigma _{i}=\varsigma \), \( m_{i}=m \). Let \( C=\left( \kappa ,q,u,m,\omega \right) \). We consider the following functional cost per area \( {\mathcal {A}} \):

$$\begin{aligned}&Z\left( T_{f},X,Y,C\right) = \int _{0}^{T_{f}}c\uplambda \left( s\right) X_{1}\left( s\right) ds \nonumber \\&\qquad +\int _{0}^{T_{f}}\left( c_{u}u\left( s\right) +\frac{c_{\omega }\omega \left( s\right) }{{\overline{\omega }}+\omega \left( s\right) }\right) {\overline{N}}\left( s\right) ds\nonumber \\&\qquad +\int _{0}^{T_{f}}c_{s}q\left( X_{1}\left( s\right) +X_{2}\left( s\right) +X_{3}\left( s\right) +X_{6}\left( s\right) \right. \nonumber \\&\qquad \left. +X_{7}\left( s\right) \right) {\mathrm{d}}s \nonumber \\&\qquad -\int _{0}^{T_{f}}g\frac{\kappa \left( s\right) m\left( s\right) \left( X_{1} \left( s\right) +X_{2} \left( s\right) +X_{3} \left( s\right) \right) }{{\overline{\kappa }}}{\mathrm{d}}s \end{aligned}$$
(50)

where \({\overline{\kappa }}=10\,{\mathrm{person}}\,\times \,{\mathrm{m}}^{-2}\) denotes the usual average density of population per unit of surface, and \({\overline{\omega }}=1440s^{-1}\) is the rate of hands washing corresponding to spending the half of time for disinfection. The negative term in the expression of Z represents the reduction of economic losses when people are free to move around in their various economic activities. A negative value of the cost means that globally the losses due to Covid-19 are not greater than the production of wealth.

We assume that an efficient hands washing takes about 30s. Thus, if \( \omega ^{-1}=60s \) the time spent on hands washing is half of the time between two consecutive washes. On the other hand, when \( \omega \) tends to \( +\infty \) then the proportion of time spent to wash hands is 1. A function of \( \omega \) displaying the behavior explained above is \( \frac{\omega \left( s\right) }{{\overline{\omega }}+\omega \left( s\right) } \) with \({\overline{\omega }}=1440s^{-1}\). We suppose that the disinfection costs is \( c_{\omega }=1.736\times 10^{-5}\,{\mathrm{USD}}\,\times \,{\mathrm{person}}^{-1}\,\times \,{\mathrm{s}}^{-1}\) in average. We also assume that after using a mask for a given time (4 h) it is systematically changed for a new mask having a certain unitary cost (\( 1\,{\mathrm{USD}}\)). Hence, \( c_{u}=0.25\,{\mathrm{USD}}\,\times \,{\mathrm{person}}^{-1}\,\times \,{\mathrm{h}}^{-1} \).

Let \( T_{i-2} \) and \( c_{i-2} \) denote, respectively, the time spent by an individual in the compartment \( i=2,\ldots ,6 \) and the cost generated by unit of time. According to probabilistic properties of Poisson processes [15, 25, 44], we have

$$\begin{aligned} T_{0}= & {} \dfrac{1}{\mu +\alpha _{1}}, \end{aligned}$$
(51)
$$\begin{aligned} T_{1}= & {} T_{1}\left( q \right) =\dfrac{\alpha _{1}T_{0}}{ q+\mu +d+\gamma _{1}}, \end{aligned}$$
(52)
$$\begin{aligned} T_{2}= & {} T_{2}\left( q \right) =\dfrac{pqT_{1}}{\alpha _{2}+\mu +d+\gamma _{2}}, \end{aligned}$$
(53)
$$\begin{aligned} T_{3}= & {} T_{3}\left( q \right) =\dfrac{\left( 1-p\right) qT_{1}+\alpha _{2}T_{2}}{\mu +d+\gamma _{3}}\text { and } \end{aligned}$$
(54)
$$\begin{aligned} T_{4}= & {} T_{4}\left( q \right) =\dfrac{\gamma _{1}T_{1}+\gamma _{2}T_{2}+\gamma _{3}T_{3}}{\mu +\gamma _{4}}. \end{aligned}$$
(55)

One can associate a cost \( c_{d} \) per person which is related to his burial ceremony. The probability to die during the stay in disease compartments I, Q and H is given by

$$\begin{aligned} p_{d}=\left( \mu +d \right) \left( T_{1}+T_{2}+T_{3}\right) . \end{aligned}$$
(56)

Hence,

$$\begin{aligned}&c=c\left( q \right) =\left( c_{1}+c_{g}\right) T_{1}+\left( c_{2}+c_{g}\right) T_{2}\nonumber \\&\quad +\,\left( c_{3}+c_{g}\right) T_{3}+c_{d}p_{d}. \end{aligned}$$
(57)

Table 4 summarizes the different unit costs we mentioned above.

Table 4 Baseline values of elementary costs

Let us consider the Hamiltonian

(58)

where the adjoint state \( \left[ \begin{matrix} X^{*} \\ Y^{*} \end{matrix} \right] \) satisfies

(59)

and

$$\begin{aligned} \left[ \begin{matrix} X^{*}\left( T_{f}\right) \\ Y^{*}\left( T_{f}\right) \end{matrix} \right] = \left[ \begin{matrix} 0 \\ 0 \end{matrix} \right] . \end{aligned}$$
(60)

The optimal control \( C^{*} \) will be given by the optimality condition

$$\begin{aligned} \nabla _{C} H =0. \end{aligned}$$
(61)

Notice that \( C^{*} \) is projected into the admissible set \({\mathcal {C}}_{Ad}=] 0,{{\overline{\kappa }}}]\times [0,{\overline{q}}] \times [0,{u_{\max }}]\times [m_{\min },1]\times [0,{\omega _{\max }}]\) where \( {\overline{q}}=48 \) (persons are tested in average every 30 minutes), \( u_{\max }=0.9 \), \( m_{\min }=0 \), \( \omega _{\max }=24 \). The existence of a unique optimal strategy is guaranteed by the fact that the strict convexity of the Hamiltonian with respect to C and the regularity of the optimal feedback strategy \( C^{*}\left( X,Y,X^{*},Y^{*}\right) \). For more details on the determination of optimal control we refer to the books [4, 49] and references therein.

Theorem 4

There is a unique piecewise continuous optimal control \(C^{*}\) for the model (1)–(20 ) characterized by Algorithm 1–.

4.2 Numerical simulations and discussion

In order to numerically determine the optimal control strategy, we implement the forward-backward sweep method described in Chapter 3, page 101 of [4]. It is an iterative approach that consists of replacing the feedback by a guess value first, solving the forward state problem, solving the backward adjoint problem and then updating the feedback using the gradient direction with a step length obtained by the golden section method. The algorithm stops when the updates of the control does not significantly change its values according to a chosen threshold. That method is similar to the one given into the Chapter 3 of [49], but it permits to better control the convergence. Notice that one could also envisage a method which minimizes the distance between the final condition of the adjoint problem obtained for a guess initial condition and the reference final condition. That method could be a multidimensional variant of the secant algorithm [26, 39].

We found after several simulations that, if we globally consider the economic losses and gains then, no specific control effort is needed additionally to take care of hospitalized patients and putting in quarantine of infectious reported individuals as it is already done. Probably that is a consequence of the high recovery rate of patients and the relatively small death probability. However, since the logistic resources for quarantine and hospitalization are limited, the economical losses in terms of gross domestic product (due to social distancing, reduction of mobility and lock-down) can be penalized in the cost functional Z. That is done in the rest of the section.

Figure 5 is an illustration of the optimal control strategy against Covid-19 and the evolution of the infection force \( \uplambda \). We can observe that social distancing and disinfection frequency are the main elements of the optimal control strategy. The frequent oscillations of \(\omega \) can be explained by the sensitivity of the controlled system behavior with respect to \(\omega \). Social distancing permits to have no restriction on mobility. The use of mask is useful but not necessary when social distancing and regular disinfection are respected. Notice that if the mask is misused it is useless, and it can promote other respiratory complications. As the sensitivity analysis showed the screening has no significant effect on disease behavior and it is economically better to keep q at its minimal value corresponding to self decision to be tested or the hospitalization of severe disease cases.

Fig. 5
figure 5

Optimal control dynamics of Covid-19 prevalence

Figure 6 shows the progression of the SARS-CoV-2 population (in terms of number of viruses per \( m^{-2} \)) and the disease prevalence either the optimal control is applied or not. Relative prevalence of disease is preferred over absolute prevalence because it is more appropriate for scaling up. The peak of the disease occurs earlier and with smaller values under optimal control policy. We also observe that without control there is a possibility to observe several peaks of the disease. Without control, the peak of the disease (\( 32\% \) prevalence) occurs after 57 days (about 8.5 weeks) while under control it occurs after 26.37 days (about 4 weeks) with a weaker magnitude (\( 0.4478\% \) prevalence). With the initial compartmental population sizes given in Table 2 and under the optimal control strategy, the disease extinction in the human population is reached after \( 9.571{\mathcal {A}}\times 10^{-3}\,{\mathrm{day}}\,\times \,{\mathrm{m}}^{-2}\) while it occurs in the virus population after \( 2.023{\mathcal {A}}\times \,10^{-1}\,{\mathrm{day}}\,\times \,{\mathrm{m}}^{-2}\). If we consider the maximal disease prevalence \( {\mathcal {I}} 1 \) and the virus population size (\( N_{V} \)) is equal to its upper bound \( C_{V} \) given in Proposition 1, then the disease extinction in the human population is reached after \( 5.499{\mathcal {A}}\times \,10^{3}\,{\mathrm{day}}\,\times \,{\mathrm{m}}^{-2} \). Similarly, the extinction in virus population is reached after \( 9.8321 {\mathcal {A}}\times 10^{4}\,{\mathrm{day}}\,\times \,{\mathrm{m}}^{-2} \). Thus, the sooner the control measures are taken, the sooner the disease can disappear. The dependence of the disease extinction time with respect to the surface \({\mathcal {A}} \) clearly highlights the importance of lock-down. As we can observe, the extinction can be very slower in the virus population than in the human population. The cost of the optimal strategy has been computed for 180 days corresponding to half of the year. It has been evaluated at \( {\mathcal {A}}Z^{\star }=3.978\,{\mathrm{USD}}\,\times \,{\mathcal {A}}\,\times \,{\mathrm{m}}^{-2}\) which corresponds on average to \(0.0221\,{\mathrm{USD}}\,\times \,{\mathcal {A}}\,\times \,{\mathrm{m}}^{-2}\times {\mathrm{day}}^{-1}\).

Fig. 6
figure 6

Controlled dynamics of Covid-19 prevalence

5 Conclusion

This work is concerned with the mathematical modeling of the severe acute respiratory corona virus 2. According to governmental data and World Health Organization, the corona virus made several millions of dead in the world. It spreads very rapidly and affects not only the public health sector, but also the global economy. The transmission of the disease can occur through direct contact with an infectious individual or through indirect contacts via the environment. The high complexity of the disease spreading suggested numerous modeling activities. Indeed, through a mathematical model it is possible to rigorously study fighting policies against corona. Hence, the goals we aimed to achieve were to propose a realistic model for the Covid-19 dynamics, to estimate the parameters of the proposed model according to the available real data, to study asymptotic behaviors depending on some featured control parameters related to WHO and governmental recommendations, and to design control strategies that are optimal with regards to socio-economical constraints.

We have constructed a two scale compartmental model with 20 classes corresponding to 10 human population states (macroscopic scale) and 10 viruses locations (microscopic scale). Indeed, we considered Susceptible \( \left( S \right) \), Infected \( \left( E \right) \), Infectious \( \left( I \right) \), Quarantined \( \left( Q \right) \), Hospitalized \( \left( H \right) \), Treated \( \left( T \right) \), Recovered \( \left( R \right) \), Non-Infectious dead \( \left( D \right) \), Infectious dead \( \left( D_{I} \right) \), Buried \( \left( B \right) \). Each individual was assumed to have an external viral load on his body while there were free viruses in the environment. The model has been labeled \( SEIQHTRDD_{I}B-V_{S}V_{E}V_{I}V_{Q}V_{H}V_{T}V_{R}V_{D}V_{D_{I}}V_{F} \) as described graphically in Fig. 1. To better access the effects of social distancing, the units of state variables have been taken in terms of person per unit area and virus per unit of surface. The parameters of the model have been successfully fitted using a Cameroonian dataset and a semi-empiric procedure. Additionally to the quantitative study of the system, we carried out a qualitative analysis of the model.

The existence of equilibriums has been studied. Contrarily to the endemic equilibrium, the disease free equilibrium (DFE) exists unconditionally. The DFE is globally asymptotically stable if a threshold \( {\mathcal {T}}_{0} \) is less or equal to one. We also provided a formula for the basic reproduction number \( {\mathcal {R}}_{0} \) which has a biological interpretation regarding the model: the average number of new viruses generated by an introduced virus during its lifespan. \( {\mathcal {T}}_{0} \) has a similar biological interpretation. A convergence index \( {\mathcal {T}}_{1} \) has been defined. Under an exponential convergence hypothesis, the times of disease extinction respectively in human population (\( T_{H}^{\star } \)) and virus population (\( T_{V}^{\star } \)) are determined. A sensitivity analysis of \( {\mathcal {R}}_{0} \), \( {\mathcal {T}}_{0} \) and \( {\mathcal {T}}_{1} \) have been done using the partial rank correlation coefficient (PRCC). We found that the sensitive parameters are in descending order the proportion of time of wearing mask, the proportion of time spent on disinfecting people and the environment, and the mobility probability. Through a concept of specific area \( {\mathcal {A}} \), the importance of lock-down is highlighted. Indeed, the times of extinction are proportional to \( {\mathcal {A}} \) knowing that \( 1-\frac{{\mathcal {A}}}{{\mathcal {A}}_{G}} \) is the lock-down level. Since the disease and the application of the different control strategies generate economical costs, an optimal control analysis has been already carried out.

In order to guarantee effectiveness of control strategies and to reduce the expenses generated by Covid-19, we defined a functional cost to be minimized. The existence and the uniqueness of the optimal control is established and its characterization is given. By the sweep decent method the optimal strategies are computed according to some practical constraints materialized by an admissible control set. We found that if one globally considers the economical losses and gains, then no specific control effort is needed additionally to take care of hospitalized patients and putting in quarantine of infectious reported individuals as it is already done. However, since the logistic resources for quarantine and hospitalization are limited, the economical losses in terms of gross domestic product (due to social distancing, reduction of mobility and lock-down) can be penalized in the functional cost. That has been done and we computed new optimal strategies.

According to the simulations, we observed that social distancing and disinfection frequency are the main elements of the optimal control strategy. Social distancing permits to have no restriction on mobility. The use of mask is useful but not necessary when social distancing and regular disinfection are respected. Screening has no significant effect on disease behavior and it is economically better to keep q at its minimal value corresponding to self decision to be tested or the hospitalization of severe disease cases. Under the optimal control policy, the peak of the disease occurs earlier and with smaller values.

The current work has several perspectives. The first one is to take into consideration the effects of the available vaccine against Covid-19 both in the asymptotic behavior analysis and in the optimization of the control strategy. Another aspect is to consider the spatial heterogeneity in terms of available logistic resources, population densities, attractiveness of specific places at given periods of a day, cultural and religious habits, and others significant factors. An ongoing work is trying to tackle the problem regarding those aspects. The model can also be extended using the promising theory of fractional calculus [1, 2, 5, 6, 59, 60, 64].