COVID-19 in regions with low prevalence and low density of population. An uncertainty dynamic modeling approach =============================================================================================================== * José M. Garrido * David Martínez-Rodríguez * Fernando Rodríguez-Serrano * Javier Díez-Domingo * COVID-19_Granada Study Group * Rafael-J. Villanueva ## Abstract The coronavirus disease 2019 (COVID-19) that emerged in China at the end of 2019 has spread worldwide. In this article, we present a mathematical SEIR model focused on analysing the transmission dynamics of COVID-19, the patients circulating in the hospitals and evaluating the effects of health policies and vaccination on the control of the pandemic. We tested the model using registered cases and population data from the province of Granada (Spain), that represents a population size near 1 million citizens with low density of population and low prevalence. After calibrating the model with the data obtained from 15 March to 22 September 2020, we simulate different vaccination scenarios - including effectiveness and availability date - in order to study the possible evolution of the disease. The results show that: 1) infected will increase until 5.6% - 7.4% of the total population over next 3-4 months (2nd wave); 2) vaccination seems not to be enough to face the pandemic and other strategies should be used; 3) we also support the claim of the WHO about the effectiveness of the vaccine, that should be, at least, of 50% to represent a substantial progress against the COVID-19; 4) after the 2nd wave, the return to normal life should be controlled and gradual to avoid a 3rd wave. The proposed study may be a useful tool for giving insight into the transmission dynamics of SARS-CoV-2 and to design vaccination and health policies. Keywords * COVID-19 * mathematical model * transmission dynamics * vaccination strategies * low density population * low prevalence ## 1 Introduction Last December, the Chinese public health authorities informed about several cases of acute respiratory syndrome in Wuhan. By 7 Jan 2020, scientists had isolated a novel coronavirus from these patients not previously identified in humans, and referred as severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the agent responsible for the later designated coronavirus disease 2019 (COVID-19) in February 2020, and the World Health Organization (WHO) declared a pandemic [1, 2]. Much is known about the virus and the disease, but there is a crucial need to further knowledge about the virus transmission to control the spread and to implement measures to reduce the impact, like social distancing or the use of facemasks, while effective antivirals and vaccines are obtained. Different mathematical models have been proposed to analyse the evolution of this pandemic although with limitations due to the many uncertainties regarding to the disease. One of the main parameters in the epidemiological models is the effective reproductive number (*R*), which quantifies the average of secondary cases generated from each infected person in a susceptible environment during an outbreak. Different studies have estimated that the *R* value in China at the beginning of the pandemic was between 2 and 7.1, and also it has been observed that *R* value drops sharply after the implementation of control measures against COVID-19 [3]. Several approaches have tried to model mathematically the transmission dynamics of COVID-19 [4, 5, 6, 7]. Other references as [8, 9, 10, 11], moreover, consider asymptomatic infected individuals, most of the infected population, who play an important role in spreading the virus. In order to make comparable our study with other similar regions, we are going to give our results in absolute values and percentages. Also, under a demographic point of view, Granada is a province in southern Spain with a population of 921511 inhabitants, around 30% of the population live in the capital and the province has a population density of 73.5*/km*2 [12] Compared with the province of Madrid with 829.6*/km*2, the province of Granada has a low density of population. And the percentage of accumulated infected between 1 May 2020 and 1 June 2020 was [2.3% - 4.8%] in Granada and [10.3% - 13.3%] in Madrid [13]. Therefore, the province of Granada has a low prevalence. In this paper we propose a mathematical model to study the current transmission dynamics of SARS-CoV-2. It is a SEIR model in which the circuit of patients moving throughout the hospital dependencies is also considered, returning a more precise portrait of the use of hospital resources with the aim of predicting hospital’s overloads. This is a novelty that our model provides and any other of the models in the literature, to our knowledge, do. Also, we consider the asymptomatic infectious individuals who transmit the disease but not reported by the system. We proved the model using registered hospital’s cases and population data from the province of Granada (Spain). Once the model was implemented and calibrated, we performed calculations simulating possible future scenarios to evaluate the effect of the COVID-19 vaccination, that allow us to provide some public health recommendations. ## 2 Model building With respect to COVID-19, an individual can be considered as: * **(S)** susceptible, when the individual is healthy; * **(Q)** in lockdown, when the individual is at home to avoid the spread of the virus. In our model, lockdown only considers the lockdown time; * **(L)** latent or exposed, when the individual has been infected but it is not infectious yet; infectious, when the individual is capable of spreading SARS-CoV-2; * **(R)** recovered, when the individual recovers from the disease being asymptomatic or having mild symptoms; * **(H)** hospitalized at ward, when the individual has severe symptoms and needs to be hospitalized; * **(u)** in intensive care unit (ICU), when the individual has severe symptoms and needs to be treated in the intensive care unit; * **(F)** deceased, when the individual dies because of the disease; * **(HU)** after ICU, when an individual is transferred from ICU to other non-ICU department due to improvement in the evolution but still requires hospitalization; * **(A)** discharged, when the individual gets better and is discharged from hospital. We are going to use data of the Spanish province of Granada provided by the hospitals H.U. Virgen de las Nieves, H.U. Clínico San Cecilio, H.U. Santa Ana in Motril, H.U. of Baza, H. of San Rafael and H.L.A. Inmaculada, generating an inclusive model which represents the whole province. The province of Granada has *P**T* = 921 511 inhabitants, and we consider the *P**T* value as constant during the simulation time. Our simulations begin on 1 March 2020 and the time step *t* is one day. 1. When the alarm state is decreed and most of the people have to be in lockdown, that is, move from *S* to *Q*, it is modeled by the term *s**q*(*t*), where the model parameter *s**q*(*t*) determines the transit of people from *S* to *Q. s**q*(*t*) takes the value 0 except for 16 March 2020 and 31 March 2020, when the lockdown and the strict lockdown began in Spain, and 700 000 and 150 000 individuals in Granada, respectively, move from S to Q [14]. 2. When the lockdown finishes, the transit of individuals from *Q* to *S* is modeled by the term *q**s*(*t*) where the model parameter *q**s*(*t*) is 0 except for 13 April 2020 when the strict lockdown finishes and 150 000 individuals move from Q to S, and from 5 May to 21 June 2020, leaving the lockdown 8 750 people every day due to the gradual end of the confinement. 3. An individual moves to latent state (L) if he/she gets infected by contact with an infectious contagions. The transit is modeled by the non-linear term ![Graphic][1], where the transmission rate parameter *β* has to be calibrated. Furthermore, this parameter will change over time due to the global public health interventions. 4. A latent individual transits to infectious state after a while and this is modeled by the linear term *l**i**L*(*t*), where the latency period *t* = 1*/l**i* is the time from the moment in which an individual is infected until the moment in which is able to transmit the virus, *t*. This period is different to the typical incubation period time from infection to onset, and according to [15], it takes from 2 to 7 days. Even though there is limited evidence about the possibility of infection one or two days before onset [16], let us consider that *t* may take values from 1 to 6 days, and *l**i* = 1*/t*. 5. An infectious individual may become hospitalized (H), admitted to Intensive Care Unit (U) or get recovered (R), and these transits are modeled by the linear terms *i**h**I*(*t*), *i**u**I*(*t*) and *i**r**I*(*t*), respectively. Here we have 3 possible ways for infectious individuals. Every one takes its time and has its probability, that is, *i**h* = *p*1*/t*1, *i**u* = *p*2*/t*2 and *i**r* = (1 *− p*1 *− p*2)*/t*3. * *p*1 is the percentage of infected who become hospitalized and *p*2 is the percentage of infected people admitted directly at ICU. Then, the parameters of the transition from *I* to *H* and from *I* to *U* are *i**h* = *p*1*/t*1 and *i**u* = *p*2*/t*2, respectively. * Those infected not hospitalized require around *t*3 = 14 days to recover [17]. Thus, the transition from *I* to *R* is governed by the parameter *i**r* = (1 *− p*1 *− p*2)*/*14. 6. People hospitalized may move to ICU (U) if get worse, may decease (F) or may be discharged (A). These transits are modeled by the linear terms *h**u**H*(*t*), *h**f* *H*(*t*) and *h**a**H*(*t*), respectively. As before, here we have 3 possible ways for hospitalized individuals. * *p*4 is the percentage of hospitalized people who need to be admitted to ICU. Hence, the people in *H* move to *U* governed by the parameter *h**u* = *p*4*/t*4. * The transition parameter from *H* to *F* is *h**f* = *p*5*/t*5. * (1 *− p*4 *− p*5) of hospitalized are discharged after an average of *t*6 days in the hospital. Hence, the transition parameter from *H* to *A* is *h**a* = (1 *− p*4 *− p*5)*/t*6. 7. People in ICU (U) may decease (F) or may get better and be transferred to other non-ICU department (HU). These transits are modeled by the linear terms *u**f* *U* (*t*) and *u**hu**U* (*t*), respectively. Here we have 2 possible ways for individuals in ICU, die or get better. * The transition parameter from *U* to *F* is *u**f* = *p*7*/t*7. * The parameter that governs the transition from *U* to *HU* is *u**hu* = (1 *− p*7)*/t*8. 8. Finally, an individual in *HU* may get better and discharged. This transit is modeled by the linear term *hu**a**HU* (*t*), where *h**ua* = 1*/t*9. Since other known coronavirus have a seasonal behaviour [18], as a hypothesis, two main periods have been considered: virulence season (September to April) and non-virulence season (May to August). Due to this, some parameter values will vary depending on the season. Taking into account the above description of the populations and the model parameters, the following system of difference equations (1) describes the transmission dynamics of COVID-19 in the province of Granada (Spain) over the time. ![Formula][2] Fig. 1 shows a flow diagram about how individuals may move throughout the different groups with respect to the disease. Note that the right part of Fig. 1 (H, U, F, A, HU) represents the flow of patients within hospital departments. There is not movement between infectious (I) to deceased (F) as there is not data available. Also, we do not consider re-infection. ![Figure 1:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/10/25/2020.10.21.20215418/F1.medium.gif) [Figure 1:](http://medrxiv.org/content/early/2020/10/25/2020.10.21.20215418/F1) Figure 1: Flow diagram of the COVID-19 transmission dynamics. The boxes represent the groups of individuals respect to the disease, and the letters next to the arrows represent the transition rates between states. ## 3 Model calibration As we have mentioned before, we gathered data from different hospitals of the Spanish province of Granada with the aim at calibrating the model parameters, in such a way that the model may describe the evolution of the pandemics in this region. Hence, the model parameters *β, t*, *p*1, *t*1, *p*2, *t*2, *p*4, *t*4, *p*5, *t*5, *t*6, *p*7, *t*7, *t*8 and *t*9 have to be calibrated in the intervals described in Table 2. Also, the initial number of latent and infectious to establish an appropriate initial condition the 1 March 2020. According to the data, at 1 March there where 5 people in the hospital circuit and none at ICU or recovered. View this table: [Table 1:](http://medrxiv.org/content/early/2020/10/25/2020.10.21.20215418/T1) Table 1: Summary of the known mathematical model parameters. View this table: [Table 2:](http://medrxiv.org/content/early/2020/10/25/2020.10.21.20215418/T2) Table 2: Summary of the mathematical model parameters to be calibrate and their calibration ranges. The calibration ranges for the hospital’s data have been obtained calculating the percentiles 10 and 90 of the data. To calibrate the model, we consider four different *β*s: before and after the lockdown, and once the confinement has finished, non-virulence and virulence season behaviour. Also, we take into account that the basic reproduction number of our model *R* = *β/*(*i**r* + *i**u* + *i**h*) lies in reliable values as *t* goes on, time *t* in days. Our main source data (hospital’s data) collects the number of daily hospitalizations, daily people in ICU, accumulated deaths and accumulated discharges in the hospitals of the province of Granada from 1 March to 22 September 2020. These data allows us to determine intervals where to search the model parameter values during the calibration. Although the mean of the intervals determined by the hospital’s data seems to be the natural option for some of the parameter values, in this case, these mean values may not describe the behavior of the system, since the data do not follow centered distributions. In order to obtain suitable numerical values of the model parameters, the 10 and 90 percentile of the stay time in each state of the hospital circuit (H, U, F, A, HU) are obtained and settled as the calibration range where looking for the proper values of the parameters. The intervals to determine the parameter values, those given by the hospital’s data and the other, can be seen in Table 2 under the column *Calibration range*. Also, we are going to use the seroprevalence study [13], which determines the percentage of accumulated infected people from the whole population is also available. According to the prevalence study of COVID-19 in Spain, carried out between May and June 2020, the accumulated infected population in the province of Granada from the 27 April to 11 May is in the range [1.7%, 4.4%]. The range from 18 May to 1 June is [2.3%, 4.8%] and the range from 8 June to 22 June is [2.3%, 5.5%]. The model has been calibrated using the optimization algorithm NS for PSO [19] adapted for this problem. Difference between real data and the model outputs - hospitalized at ward, ICU, deceased, discharged and recovered - has been minimized in order to obtain the model parameter values. We have calibrated the model 600 times with 100 000 model evaluations each one, and for all the calibrated sets of model parameter values (all of them valid because lie on the calibration range intervals) we calculate the mean and the 95% confidence interval (CI95%) by the percentiles 2.5 and 97.5 in each time instant *t*. This way, we are able to capture and quantify the possible data uncertainty. The results of the calibration are good enough to represent the effects of the pandemic in Granada during 1 March - 22 September 2020 period (Fig. 2). ![Figure 2:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/10/25/2020.10.21.20215418/F2.medium.gif) [Figure 2:](http://medrxiv.org/content/early/2020/10/25/2020.10.21.20215418/F2) Figure 2: Model calibration. (A) Number of daily hospitalized at ward people. (B) Number of daily people in ICU. (C) Number of accumulated deceased people. (D) Number of accumulated discharged people. (E) Number of accumulated recovered people. (F) Transmission rate - *β*. The black bands are the mean and the 95%CI of the model outputs. The red points are the data gathered from the hospitals. Vertical red lines of (E) show the confidence intervals of the prevalence study [13] carried out in Spain last spring. ### 3.1 Some notes on the calibrated model parameter values The mean and confidence interval of the calibrated parameters is shown in Table 3. The availability of reliable hospital’s data allowed us to provide realistic intervals to calibrate the model parameters and, once the model has been calibrated, meaningful model parameter values. View this table: [Table 3:](http://medrxiv.org/content/early/2020/10/25/2020.10.21.20215418/T3) Table 3: Summary of the calibrated model parameter values. In 1 March 2020 the calibrated number of infectious 1190, *CI*95% 1165 *−* 1222, (0.13%, *CI*95% 0.125% *−*0.133%) and latent 573, *CI*95% 558 *−* 590, (0.062%, *CI*95% 0.060% *−* 0.064%) was low. Data provided by the hospitals in the province of Granada, from 1 March to 15 March, compared with the registry of the previous year did not show an increase of respiratory disease cases that could be attributed to COVID-19, what it is in accordance with our estimation of the low number of COVID-19 patients at the beginning of the pandemic (1 March 2020). Later, the number of infected patients experienced a great increase, what supports that the effective reproduction number *R* was initially very high, (5.42, *CI*95% 5.35 *−* 5.54), value which is in accordance with Sanche et al [20]. After the beginning of the lockdown declared by the Spanish government in 14 March, the number of infected patients kept growing as a consequence of the contagions produced in the previous weeks with high *β*. Nevertheless, the basic reproductive number *R* decreased to (0.55 CI95% 0.504 *−* 0.597), not only because of lockdown but also the social distancing, the use of facemasks and other population measures, reducing the *β* (lower right graphic of Fig. 2). During the non-virulence season after the lockdown is finished, the incidence of the disease is kept low, with a stabilization of the number of hospitalizations and deceases. The *R* value is low, (1.40 CI95% 1.367*−* 1.427), but it increases until slightly higher values than at the beginning of the virulence season (*R* = 1.73, CI95% 1.597*−* 1.879). This can be observed with the increment of hospitalized and intensive care unit people, together with the end of the stagnation in the number of deceased and discharged on Figure 2. Different values of the parameters can be appreciated, more precisely the percentage of infected people who needs hospitalization (hospitalization at ward *p*1, or ICU treatment *p*2). During the virulence season, about 2% of the infected people need hospitalization at ward, and 0.2% of the infected people enter the ICU directly without passing through the hospitalization ward. During the non-virulence season, about 1.4% of the infected people need hospitalization at ward, and 0.06% of the infected people enter the ICU directly without passing through the hospitalization ward. ## 4 Vaccination Once the model has been calibrated, we can simulate different possible scenarios for the future. In particular, taking into account that several vaccines candidates are being created, we are going to deal with vaccination strategies. In order to simulate the effect of the vaccine in the population of the Spanish region of Granada, a new state has to be considered in the mathematical model: * **(V)**vaccinated, when the individual is vaccinated and the vaccine is effective, protecting the individual. We suppose that only susceptible people are vaccinated. The vaccine protects against the most dangerous symptoms of the COVID-19, and stops the spread of the disease. However, the effect of the vaccine is not permanent and the vaccinated people may become susceptible again. The transit from *S* to *V* is modeled by the term *s**v*(*t*) and the transit from *V* to *S*, it is modeled by the term *v**s*(*t*). Both parameters describe the number of people who are moving from one state to other, at time *t*. The mathematical system of difference equations, which includes the vaccination of the population is shown in system (2). In black, the new terms. ![Formula][3] The administration capacity of vaccination is limited and it is not possible to vaccinate the whole population on the same day. We assume that, being this pandemic an extraordinary situation, it would be possible to vaccinate 1% of the population every day. Not every vaccinated person acquires protection against the disease. We consider the vaccine effectiveness as the percentage of vaccinated people who acquire immunity. Also, the duration effect of the vaccine is unknown, and we will assume an immunity of a year on all cases. To simulate the possible scenarios for the population vaccination, as it is still not known the possible effectiveness of the future COVID-19 vaccines, 100%, 75%, 50% and 25% of effectiveness have been compared with the non vaccination scenario. All the scenarios begin to administer the vaccine on 1 December 2020. The results of this simulation have been presented in Figure 3. ![Figure 3:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/10/25/2020.10.21.20215418/F3.medium.gif) [Figure 3:](http://medrxiv.org/content/early/2020/10/25/2020.10.21.20215418/F3) Figure 3: COVID-19 evolution depending on the vaccination effectiveness (starting on 1/12/2020): No vaccination, 25%, 50%, 75% 100%. (A) Number of infected people. (B) Number of daily hospitalized at ward people. (C) Number of daily people in ICU. The percentage of infected population on the peak is 7.4%, 6.5%, 6.1%, 5.7%, 5.6% for no vaccination and a vaccine effectiveness of 25%, 50%, 75% 100%, respectively. On Figure 3, the difference of the vaccine effectiveness effect is shown. As the effectiveness of the vaccine increases, the height of the peak decreases. It is important to remark the importance of achieving 50% of vaccine effectiveness, since the major reduction is seen from that point. The difference between the hospitalization peak from the scenario of non vaccination and the scenario of 50% effectiveness is 20.6% of reduction, meanwhile the difference on the hospitalization peak from the scenario of 50% effectiveness to the scenario of 100% effectiveness is only 9.6% of reduction. This results are in accordance with the WHO recommendation [21]. The availability of the vaccine is still unknown, but, according to some declarations of the Spanish government [22], it could be available at the beginning of December 2020. Anyway, we have simulated two availability scenarios: December 2020 and February 2021, assuming a vaccine effectiveness of 50%. The election of the vaccine effectiveness is based on the WHO recommendation where it is considered a substantial progress a vaccine effectiveness of 50% [21]. The results have been compared with no vaccination scenario in Figure 4. ![Figure 4:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/10/25/2020.10.21.20215418/F4.medium.gif) [Figure 4:](http://medrxiv.org/content/early/2020/10/25/2020.10.21.20215418/F4) Figure 4: COVID-19 evolution depending on the vaccination start date (with a 50% of effectiveness): No vaccination, vaccination from 1/12/2020 and vaccination from 1/2/2021. (A) Number of infected people. (B) Number of daily hospitalized at ward people. (C) Number of daily people in ICU. The difference between starting the vaccination campaign on 1 December 2020 or on 1 February 2021 is remarkable. As it can be seen in Figure 4, the difference between starting the vaccination campaign on 1 December 2020 or on 1 February 2021 is quite remarkable. There is 18.3% of reduction of infected people, 20.4% of reduction in hospitalized and 21.4% of reduction in ICU on the pandemics’ peak. There is also a time advance of 19 days on the infected and hospitalized peak and a time advance of 20 days on the ICU peak if the vaccination begins in December 1st instead of February 1st. It should be also noticed that, if the vaccination starts on 1 February 2021, the difference between the vaccination and non vaccination scenarios is only seen at the end of the pandemics, with no affectation on the peak. It is important to remark that, according the the date of release of the vaccine, the effect on this winter pandemic will be reduced. Herd immunity effect will be limited since the infection evolves faster than the vaccination campaign. Because of this fact, it would be necessary to increase and combine health measures in order to decrease the spread of the illness until a vaccine could be available to the population. ## 5 Discussion Coronaviruses cause respiratory and intestinal diseases in many animal species. In humans, cold-like symptoms and upper respiratory tract infections can be caused by 4 existing coronaviruses, including 2 alpha coronaviruses (229E and NL63) and 2 beta coronaviruses (OC43 and HKU1). In addition, there are 3 beta coronaviruses that cause severe respiratory syndromes in humans (SARS-CoV-1, MERS-CoV, and SARS-CoV-2) [23]. Main treatment options suggested against SARS-CoV-2 include chloroquine/hydroxychloroquine, favipiravir, remdesivir, lopinavir/ritonavir, umifenovir, steroids, cepharanthine, interferon, ivermectin, ribavirin, convalescent plasma, anticoagulants, mesenchymal stem cells and neutralizing antibodies, among others. However, there is an urgent need for effective rapid diagnostics, vaccines and therapeutics to detect, prevent, treat and contain SARS-CoV-2 [24, 25, 26]. Explanatory mathematical models of diseases must allow to predict the evolution of the disease considering a wide range of possible scenarios. It must be able to estimate the healthcare demands that will have to be covered by hospitals early enough as to implement management and operative adaptations regarding multiple aspects, including the extension of ICU to other places that could be equipped and be used for the management of the COVID-19 patients, as well as the reservation of out-of-hospital stays, such as hotels, for health personnel or even for patients with low severity who must be lockdown in a controlled context [27]. In this paper we propose a mathematical model to study the current transmission dynamics of SARS-CoV-2. It is a SEIR model where the circuit of patients moving inside the hospital dependencies, the particular features of the disease and the asymptomatic infectious, a principal source of contagiousness [28], are considered, returning a more precise portrait of the use of hospital resources. We tested the model using registered cases and population data from the hospitals of the Spanish province of Granada that represent a size of near 1 million citizens. In this province, there are low density of population and low prevalence of COVID-19 on the spring 2020 outbreak. A direct and permanent line of communication is important for the transmission of data in real time between the heads of hospital centres and the teams in charge of analysing the epidemiological evolution of the disease, especially through tools such as the one presented in this paper. This will ensure the optimization of future predictions on the basis of the augmentation of the knowledge database, contributing to face the fluctuations that may appear in the transmission of the virus depending on variables that are not sufficiently clarified at the moment. One of the strong aspects of our model is that it has been possible to calibrate and validate from data that has been collected and reported through an informative protocol agreed at the province of Granada, which ensures homogeneity as well as prevents biases that may appear on the count of each type of subject, consideration of infected by diagnostic, etc; also, the model can be fed with new data to update the predictions regularly. Another important aspect is that our model is able to distinguish patients recovered from COVID-19 and those who have died from the disease, since it will depend on the behaviour of the virus in the future, so that, the group of recovered can be considered as a group immune to infection, or they may again be considered as a susceptible population in the event of significant phenotypic changes in the virus or possible reinfections. Our simulations show that we have to keep population measures in order to control the pandemic and avoid a new concerning outbreak. After the beginning of the lockdown declared by the Spanish government in 14 March, the number of infected patients kept growing as a consequence of the contagions produced in the previous weeks of high *β*. However, the outbreak could be delayed if the disease *β* is reduced as much as possible, eluding the coincidence with the influenza’s seasonal peak. This will permit a better planning of the hospital resources, preparing a COVID-19 area without reducing the care attention to patients with other diseases, preventing the collapse of the health services, as happened last March-April 2020. The transmission capacity of the virus at the beginning of the pandemic was very high, with a *β* close to 0.4 and a *R* greater than 5. Numerous studies have stated the effectiveness of non-pharmaceutical measures in reducing the transmission of respiratory virus [29, 30, 31, 32, 33, 34, 35, 36, 37, 38]. Then, *β* and *R* decreased significantly, reaching an average of 1.40 and making it possible to control the pandemic by approaching to 1. The eradication of the virus in our environment is not to be expected in the short term, since it has not been possible to implement socio-sanitary measures that reduce the *R* and its confidence interval of potential values below 1. By 2 October 2020, there are 42 candidate vaccines for SARS-CoV-2 in clinical evaluation and 151 in preclinical models [21] that have been developed using a wide range of technologies and strategies, such as RNA, DNA, nanoparticles, recombinant proteins, live attenuated vaccines, inactivated vaccines, synthetic and modified virus-like particles [25, 39]. However, the impact of the COVID-19 vaccines on infection, re-infection, transmission, and disease severity has not been presented yet as conclusions of clinical trials [40]. Computational model of SARS-CoV-2 and vaccination in USA revealed that to extinguish the ongoing epidemic without other non-pharmaceutical interventions, the vaccine efficacy has to achieve at least 60% when coverage is 100%, and at least 80% when coverage is up to 75% to drop the peak by 85%-86%, 61%-62%, and 32% depending on the population previously exposed to the virus (5%, 15%, and 30%, respectively) [41]. WHO recommends a vaccine efficacy greater than 30% in 95% confidence interval estimation, and considers that an efficacy around 50% would represent substantial progress [40, 42]. The evolution of the pandemic with the current number of infected people and transmission rate, even though if an effective vaccine would be available, seems to be rising. More health policies must be implemented in order to decrease the transmission rate until the vaccine will be available for the whole population. Moreover, more vulnerable population groups (such as old people) should be vaccinated as soon as possible, since herd effect might be difficult to obtain before an outbreak, unless more drastic health policies are implemented. Finally, we have to take into account that, as shown in Figure 4, in summer 2021 when we reach low levels of infections, the social distancing, facemasks and other non-pharmaceutical measures cannot be removed suddenly. Otherwise, the transmission rate increases, also *R*, and a new outbreak may arise (as it can be seen in Figure 5) at the beginning of 2022. ![Figure 5:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/10/25/2020.10.21.20215418/F5.medium.gif) [Figure 5:](http://medrxiv.org/content/early/2020/10/25/2020.10.21.20215418/F5) Figure 5: COVID-19 evolution with the initial transmission rate from September 2021 due to the sudden end of non-pharmaceutical measures. Vaccine campaign 2020 starting on 1 December 2020 with 50% effectiveness and a year effect and no vaccination campaign in 2021. ## 6 Conclusion The modified SEIR model presented in this paper can be a useful tool for providing insight into the transmission dynamics of SARS-CoV-2. It considers the circuit of patients moving inside the hospital dependencies, returning a more precise portrait of the use of hospital resources, in a region with low density of population and low COVID-19 prevalence. Once the model is calibrated, vaccination simulations have been performed to evaluate the evolution of COVID-19. The simulations show that the number of infected is going to increase reaching a peak in mid December - mid January. Also, the vaccination should be combined with other health measures to face COVID-19 and avoid overloading in the health system. Furthermore, our model supports the WHO recommendation about the effectiveness of the vaccine, that should be, at least, of 50%. Moreover, when the next outbreak declines, we will not yet be out of danger, and we will have to resume life as we knew it before the pandemic relaxing measures of social situation and masks in a gradual and controlled way to avoid a new outbreak. The predictions may provide enough time to the health systems to establish the appropriate measures to face outbreaks that could collapse human and material resources. Moreover, future pharmaceutical and non-pharmaceutical interventions, as well as aspects related to SARS-CoV-2 biology and behaviour, can be represented through the model parameters to act as modulator of the *β* what provide the required flexibility to adapt the model throughout the time. ### Limitations of the model The proposed model is a classical system of difference equations. For its building, we assume usual hypothesis as the homogeneous mixing of the population (any individual may infect any individual). No age groups are considered, and the hospitalization and decease rate may vary depending on the age of the individuals. Also, we assume that all lockdown people are susceptible, when latent, asymptomatic infectious and recovered are also in lockdown. Furthermore, mobility and other spatial aspects are not considered. Moreover, it is important to consider that the model must be particularized for each population, either provincial or larger, since the transmission rate must respond to the effect of health policy measures and the response of the population to the recommendations, such as opening educational centers, social distancing, use of face masks, adaptation of public and private services to establish physical barriers that hinder the spread of the virus, etc. ## Supporting information Ethics committee approval [[supplements/215418_file02.pdf]](pending:yes) ## Data Availability All data used in the study is provided in the manuscript. ## Funding This publication was made possible by grants from * the Spanish Ministerio de Economía, Industria y Competitividad (MINECO), the Agencia Estatal de Investigacióon (AEI) and Fondo Europeo de Desarrollo Regional (FEDER UE) grant MTM2017-89664-P; * the European Union through the Operational Program of the [European Regional Development Fund (ERDF) / European Social Fund (ESF)] of the Valencian Community 2014-2020. Files: GJIDI/2018/A/010 and GJIDI/2018/A/009; * the Ramóon Areces Foundation, Madrid, Spain (CIVP18A3920). ## Conflict of interest None. ## Footnotes * * COVID-19_Granada Study Group: José Manuel Garrido Jiménez, Antonio Cansino Osuna, Juan Carlos Carrillo Santos, Cristina Carvajal Pedrosa, Ma Angeles García Rescalvo, Francisco José Guerrero García, Sebastián Manzanares Galán, Francisco Marti Jiménez, Manuel Enrique Reyes Nadal, Pedro Manuel Ruiz Lorenzo, José Luis Salcedo Lagullón, Indalecio Sánchez-Montesinos García, Fernando Rodríguez-Serrano. * ***emails:*** josemgarrido{at}ugr.es, damarro3{at}upv.es, fernrs{at}ugr.es, javier.diez{at}fisabio.es, rjvillan{at}imm.upv.es * Received October 21, 2020. * Revision received October 21, 2020. * Accepted October 25, 2020. * © 2020, Posted by Cold Spring Harbor Laboratory This pre-print is available under a Creative Commons License (Attribution-NonCommercial-NoDerivs 4.0 International), CC BY-NC-ND 4.0, as described at [http://creativecommons.org/licenses/by-nc-nd/4.0/](http://creativecommons.org/licenses/by-nc-nd/4.0/) ## References 1. [1]. T. Zhang, Q. Wu, Z. Zhang, Probable Pangolin Origin of SARS-CoV-2 Associated with theCOVID-19 Outbreak, Current Biology 30 (7) (2020) p1346 – 1351.e2. doi:[https://doi.org/10.1016/j.cub.2020.03.022](https://doi.org/10.1016/j.cub.2020.03.022). URL [http://www.sciencedirect.com/science/article/pii/S0960982220303602](http://www.sciencedirect.com/science/article/pii/S0960982220303602) 2. [2]. F. Zhou, T. Yu, R. Du, G. Fan, Y. Liu, Z. Liu, J. Xiang, Y. Wang, B. Song, X. Gu, L. Guan, Y. Wei, H. Li, X. Wu, J. Xu, S. Tu, Y. Zhang, H. Chen, B. Cao, Clinical course and risk factorsfor mortality of adult inpatients with COVID-19 in Wuhan, China: a retrospective cohort study, The Lancet 395 (10229) (2020) 1054–1062. doi:10.1016/s0140-6736(20)30566-3. URL [https://doi.org/10.1016/s0140-6736(20)30566-3](https://doi.org/10.1016/s0140-6736(20)30566-3) [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/s0140-6736(20)30566-3&link_type=DOI) 3. [3]. E. Shim, A. Tariq, W. Choi, Y. Lee, G.Chowell, Transmission potential and severity of COVID-19 in South Korea, International Journal of Infectious Diseases 93 (2020) 339–344. doi:[https://doi.org/10.1016/j.ijid.2020.03.031](https://doi.org/10.1016/j.ijid.2020.03.031). URL [http://www.sciencedirect.com/science/article/pii/S1201971220301508](http://www.sciencedirect.com/science/article/pii/S1201971220301508). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.ijid.2020.03.031&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F10%2F25%2F2020.10.21.20215418.atom) 4. [4]. M. Gatto, E. Bertuzzo, L. Mari, S. Miccoli, L. Carraro, R. Casagrandi, A. Rinaldo, Spread and dynamics of the COVID-19 epidemic in Italy: Effects of emergency containment measures, Proceedings of the National Academy of Sciences 117 (19) (2020) 10484–10491. arXiv:[https://www.pnas.org/content/117/19/10484.full.pdf](https://www.pnas.org/content/117/19/10484.full.pdf), doi:10.1073/pnas.2004978117. URL [https://www.pnas.org/content/117/19/10484](https://www.pnas.org/content/117/19/10484) [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoicG5hcyI7czo1OiJyZXNpZCI7czoxMjoiMTE3LzE5LzEwNDg0IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjAvMTAvMjUvMjAyMC4xMC4yMS4yMDIxNTQxOC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 5. [5]. B. Tang, X. Wang, Q. Li, N. L. Bragazzi, S. Tang, Y. Xiao, J. Wu, Estimation of the transmission risk of the 2019-nCoV and its implication for public health interventions, Journal of Clinical Medicine 9 (2) (2020) 462. doi: 10.3390/jcm9020462.URL [https://doi.org/10.3390/jcm9020462](https://doi.org/10.3390/jcm9020462) [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3390/jcm9020462&link_type=DOI) 6. [6]. Z. Liu, P. Magal, O. Seydi, G. Webb, A COVID-19 epidemic model with latency period, Infectious Disease Modelling 5 (2020) 323–337. doi: 10.1016/j.idm.2020.03.003. URL [https://doi.org/10.1016/j.idm.2020.03.003](https://doi.org/10.1016/j.idm.2020.03.003) [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.idm.2020.03.003&link_type=DOI) 7. [7]. K. Prem, Y. Liu, T. W. Russell, A. J. Kucharski, R. M. Eggo, N. Davies, M. Jit, P. Klepac, S. Flasche, S. Clifford, C. A. B. Pearson, J. D. Munday, S. Abbott, H. Gibbs, A. Rosello, B. J. Quilty, T. Jombart, F. Sun, C. Diamond, A. Gimma, K. van Zandvoort, S. Funk, C. I. Jarvis, W. J. Edmunds, N. I. Bosse, J. Hellewell, The effect of control strategies to reduce social mixing on outcomes of the COVID-19 epidemic in wuhan, china: a modelling study, The Lancet Public Health 5 (5) (2020) e261–e270. doi:10.1016/s2468-2667(20)30073-6. URL [https://doi.org/10.1016/s2468-2667(20)30073-6](https://doi.org/10.1016/s2468-2667(20)30073-6) [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/s2468-2667(20)30073-6&link_type=DOI) 8. [8]. M. Serhani, H. Labbardi, Mathematical modeling of COVID-19 spreading with asymptomatic infected and interacting peoples, Journal of Applied Mathematics and Computing (Aug. 2020). doi:10.1007/s12190-020-01421-9. URL [https://doi.org/10.1007/s12190-020-01421-9](https://doi.org/10.1007/s12190-020-01421-9) [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s12190-020-01421-9&link_type=DOI) 9. [9]. U. Avila-Ponce de León, A. G. Pérez, E. Avila-Vales, An SEIARD epidemic model for COVID19 in Mexico: Mathematical analysis and state-level forecast, Chaos, Solitons & Fractals 140 (2020) 110165. doi:[https://doi.org/10.1016/j.chaos.2020.110165](https://doi.org/10.1016/j.chaos.2020.110165). URL [http://www.sciencedirect.com/science/article/pii/S0960077920305610](http://www.sciencedirect.com/science/article/pii/S0960077920305610) 10. [10]. K. Sarkar, S. Khajanchi, J. J. Nieto, Modeling and forecasting the COVID-19 pandemic in India, Chaos, Solitons & Fractals 139 (2020) 110049.doi:[https://doi.org/10.1016/j.chaos.2020.110049](https://doi.org/10.1016/j.chaos.2020.110049). URL [http://www.sciencedirect.com/science/article/pii/S096007792030446X](http://www.sciencedirect.com/science/article/pii/S096007792030446X) 11. [11]. G. Giordano, F. Blanchini, R. Bruno, P. Colaneri, A. D. Filippo, A. D. Matteo, M. Colaneri, Modelling the COVID-19 epidemic and implementation of population-wide interventions in Italy, Nature Medicine 26 (6) (2020) 855–860. oi:doi:10.1038/s41591-020-0883-7.URL [https://doi.org/10.1038/s41591-020-0883-7](https://doi.org/10.1038/s41591-020-0883-7) [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41591-020-0883-7&link_type=DOI) 12. [12].Province of Granada. URL [https://en.wikipedia.org/wiki/Province\_of\_Granada](https://en.wikipedia.org/wiki/Province_of_Granada) 13. [13]. M. Pollán, B. Pérez-Gómez, R. Pastor-Barriuso, J. Oteo, M. A. Hernán, M. Pérez-Olmeda, J. L. Sanmartín, A. Fernández-García, I. Cruz, N. F. de Larrea, M. Molina, F. Rodríguez- Cabrera, M. Martín, P. Merino-Amador, J. L. Paniagua, J. F. Munõz-Montalvo, F. Blanco, R. Yotti, F. Blanco, R. G. Fernández, M. Martín, S. M. Navarro, M. Molina, J. F. Munõz- Montalvo, M. S. Hernández, J. L. Sanmartín, M. Cuenca-Estrella, R. Yotti, J. L. Paniagua, N. F. de Larrea, P. Fernández-Navarro, R. Pastor-Barriuso, B. Pérez-Gómez, M. Pollán, A. Avellón, G. Fedele, A. Fernández-García, J. O. Iglesias, M. T. P. Olmeda, I. Cruz, M. E. F. Martinez, F. D. Rodríguez-Cabrera, M. A. Hernán, S. P. Fernández, J. M. R. Aguirre, J. M. N. Marí, B. P. Borrás, A. B. P. Jiménez, M. Rodríguez-Iglesias, A. M. C. Gascón, M. L. L. Alcaine, I. D. Suárez, O. S. Álvarez, M. R. Pérez, M. C. Sanchís, C. J. V. Gomila, L. C. Saladrigas, A. H. Fernández, A. Oliver, E. C. Feliciano, M. N. G. Quintana, J. M. B. Fernández, M. A. H. Betancor, M. H. Febles, L. M. Martín, L.-M. L. López,T. U. Miota, I. D. B. Población, M. S. C. Pérez, M. N. V. Fernández, T. M. Enríquez, M. V. Arranz, M. D.-G. González, I. Fernández-Natal, G. M. Lobón, J. L. M. Bellido, P. Ciruela, A. M. i Casals, M. D. Botías, M. A. M. Maeso, D. P. del Campo, A. F. de Castro, R. L. Ramírez, M. F. E. Retamosa, M. R. González, M. S. B. Lobeiras, A. F. Losada, A. Aguilera, G. Bou, Y. Caro, N. Marauri, L. M. S. Blanco, I. del Cura González, M. H. Pascual, R. A. Fernández, P. Merino-Amador, N. C. Castro, A. T. Lizcano, C. R. Almagro, M. S. Hernández, N. A. Elizaga, M. E. Sanz, C. E. Baquedano, A. B. Bascaran, S. I. Tamayo, L. E. Otazua, R. B. Benarroch, J. L. Flores, A. V. de la Villa, Prevalence of SARS-CoV-2 in Spain (ENE-COVID): a nationwide, population-based seroepidemiological study, The Lancet 396 (10250) (2020) 535–544. doi:10.1016/s0140-6736(20)31483-5.URL [https://doi.org/10.1016/s0140-6736(20)31483-5](https://doi.org/10.1016/s0140-6736(20)31483-5) [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/s0140-6736(20)31483-5&link_type=DOI) 14. [14].Google, Google movility report. (accessed 5 april 2020). URL [https://www.google.com/covid19/mobility/](https://www.google.com/covid19/mobility/) 15. [15]. Wei-jie Guan and Zheng-yi Ni and Yu Hu and Wen-hua Liang and Chun-quan Ou and Jian-xing He and Lei Liu and Hong Shan and Chun-liang Lei and David S.C. Hui and Bin Du and Lanjuan Li and Guang Zeng and Kwok-Yung Yuen and Ru-chong Chen and Chun-li Tang and Tao Wang and Ping-yan Chen and Jie Xiang and Shi-yue Li and Jin-lin Wang and Zi-jing Liang and Yi-xiang Peng and Li Wei and Yong Liu and Ya-hua Hu and Peng Peng and Jian-ming Wang and Ji-yang Liu and Zhong Chen and Gang Li and Zhi-jian Zheng and Shao-qin Qiu and Jie Luo and Chang-jiang Ye and Shao-yong Zhu and Nan-shan Zhong, Clinical Characteristics of Coronavirus Disease 2019 in China, New England Journal of Medicine 382 (18) (2020) 1708–1720. doi:10.1056/nejmoa2002032. URL [https://doi.org/10.1056/nejmoa2002032](https://doi.org/10.1056/nejmoa2002032) [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1056/nejmoa2002032&link_type=DOI) 16. [16]. R. M. Anderson, H. Heesterbeek, D. Klinkenberg, T. D. aHollingsworth, How will country-based mitigation measures influence the course of the COVID-19 epidemic?, The Lancet 395 (10228) (2020) 931–934. doi:10.1016/s0140-6736(20)30567-5. URL [https://doi.org/10.1016/s0140-6736(20)30567-5](https://doi.org/10.1016/s0140-6736(20)30567-5) [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/s0140-6736(20)30567-5&link_type=DOI) 17. [17].World Health Organization (WHO), Q&A on coronaviruses (COVID-19). URL [https://www.who.int/emergencies/diseases/novel-coronavirus-2019/question-and-answers-hub/q-a-detail/q-a-coronaviruses](https://www.who.int/emergencies/diseases/novel-coronavirus-2019/question-and-answers-hub/q-a-detail/q-a-coronaviruses) 18. [18]. A. W. D. Edridge, J. Kaczorowska, A. C. R. Hoste, M. Bakker, M. Klein, K. Loens, M. F. Jebbink, A. Matser, C. M. Kinsella, P. Rueda, M. Ieven, H. Goossens, M. Prins, P. Sastre, M. Deijs, L. van der Hoek, Seasonal coronavirus protective immunity is short-lasting, Nature Medicine (aSep. 2020). doi:10.1038/s41591-020-1083-1. URL [https://doi.org/10.1038/s41591-020-1083-1a](https://doi.org/10.1038/s41591-020-1083-1a) [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41591-020-1083-1&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32929268&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F10%2F25%2F2020.10.21.20215418.atom) 19. [19]. D. Martínez-Rodríguez, J. M. Colmenar, J. I. Hidalgo, R.-J. Villanueva Micó, S. Salcedo-Sanz, Particle swarm grammatical evolution for energy demand estimation, Energy Science & Engineering 8 (4) (2020) 1068–1079. arXiv:[https://onlinelibrary.wiley.com/doi/pdf/10.1002/ese3.568](https://onlinelibrary.wiley.com/doi/pdf/10.1002/ese3.568), doi:10.1002/ese3.568. URL [https://onlinelibrary.wiley.com/doi/abs/10.1002/ese3.568](https://onlinelibrary.wiley.com/doi/abs/10.1002/ese3.568) [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/ese3.568&link_type=DOI) 20. [20]. S. Sanche, Y. T. Lin, C. Xu, E. Romero-Severson, N. Hengartner, R. Ke, High Contagiousness and Rapid Spread of Severe Acute Respiratory Syndrome Coronavirus 2, Emerging Infectious Diseases 26 (7) (aJul. 2020). doi:10.3201/eid2607.200282. URL [https://doi.org/10.3201/eid2607.200282](https://doi.org/10.3201/eid2607.200282) [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3201/eid2607.200282&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32255761&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F10%2F25%2F2020.10.21.20215418.atom) 21. [21].World Health Organization (WHO), Draft landscape of COVID-19 candidate vaccines. 2[https://www.who.int/publications/m/item/draft-landscape-of-covid-19-candidate-vaccines](https://www.who.int/publications/m/item/draft-landscape-of-covid-19-candidate-vaccines) October 2020. (accessed 8 October 2020).URL [https://www.who.int/publications/m/item/draft-landscape-of-covid-19-candidate-vaccines](https://www.who.int/publications/m/item/draft-landscape-of-covid-19-candidate-vaccines) 22. [22]. El País, Sánchez espera que espan£ empiece a suministrar la vacuna contra la covid-19 en diciembre (accessed 8 october 2020). URL [https://elpais.com/sociedad/2020-09-07/sanchez-espera-que-espana-empiece-a-suministrar-la-vacuna-contra-la-covid-19-en-diciembre.html](https://elpais.com/sociedad/2020-09-07/sanchez-espera-que-espana-empiece-a-suministrar-la-vacuna-contra-la-covid-19-en-diciembre.html) 23. [23]. B. M. Bennet, J. Wolf, R. Laureano, R. S. Sellers, Review of current vaccine development strategies to prevent coronavirus disease 2019 (COVID-19), Toxicologic Pathology (2020) 019262332095909 doi: 10.1177/0192623320959090. URL [https://doi.org/10.1177/0192623320959090](https://doi.org/10.1177/0192623320959090) [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1177/0192623320959090&link_type=DOI) 24. [24]. G. Aktas, A comprehensive review on rational and effective treatment strategies against an invisible enemy; SARS-Cov-2 infection, Experimental Biomedical Research 3 (4) (2020) 293–311. doi:10.30714/j-ebr.2020463629. URL [https://doi.org/10.30714/j-ebr.2020463629](https://doi.org/10.30714/j-ebr.2020463629) [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.30714/j-ebr.2020463629&link_type=DOI) 25. [25]. J. Pang, M. X. Wang, I. Y. H. Ang, S. H. X. Tan, R. F. Lewis, J. I.-P. Chen, R. A. Gutierrez, S. X. W. Gwee, P. E. Y. Chua, Q. Yang, X. Y. Ng, R. K. Yap, H. Y. Tan, Y. Y. Teo, C. C. Tan, A. R. Cook, J. C.-H. Yap, L. Y. Hsu, Potential rapid diagnostics, vaccine and therapeutics for 2019 novel coronavirus (2019-nCoV): A systematic review, Journal of Clinical Medicine 9 (3) (2020) 623. doi:10.3390/jcm9030623. URL[https://doi.org/10.3390/jcm9030623a](https://doi.org/10.3390/jcm9030623a) [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3390/jcm9030623&link_type=DOI) 26. [26]. G. Mihaescu, M. C. Chifiriuc, C. Iliescu, C. O. Vrancianu, L.-M. Ditu, L. G. Marutescu, R. Grigore, S. Bertesteanu, M. Constantin, G. G. Pircalabioru, SARS-CoV-2: From structure to pathology, host immune response and therapeutic management, Microorganisms 8 (10) (2020) 1468. doi:10.3390/microorganisms8101468. URL [https://doi.org/10.3390/microorganisms8101468](https://doi.org/10.3390/microorganisms8101468) [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3390/microorganisms8101468&link_type=DOI) 27. [27]. G. E. Weissman, A. Crane-Droesch, C. Chivers, T. Luong, A. Hanish, M. Z. Levy, J. Lubken, M. Becker, M. E. Draugelis, G. L. Anesi, P. J. Brennan, J. D. Christie, C. W. Hanson, M. E. Mikkelsen, S. D. Halpern, Locally informed simulation to predict hospital capacity needs during the COVID-19 pandemic, Annals of Internal Medicine 173 (1) (2020) 21–28. doi:10.7326/m20-1260. URL [https://doi.org/10.7326/m20-1260](https://doi.org/10.7326/m20-1260) [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.7326/M20-1260&link_type=DOI) 28. [28]. Y. Bai, L. Yao, T. Wei, F. Tian, D.-Y. Jin, L. Chen, M. Wang, Presumed Asymptomatic Carrier Transmission of COVID-19, JAMA 323 (14) (2020) 1406–1407. arXiv: [https://jamanetwork.com/journals/jama/articlepdf/2762028/jama\_bai\_2020\_ld\_200013.pdf](https://jamanetwork.com/journals/jama/articlepdf/2762028/jama\_bai_2020_ld_200013.pdf), doi:10.1001/jama.2020.2565.URL [https://doi.org/10.1001/jama.2020.2565](https://doi.org/10.1001/jama.2020.2565) [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1001/jama.2020.2565&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32083643&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F10%2F25%2F2020.10.21.20215418.atom) 29. [29]. S.-C. Chen, C.-M. Liao, Modelling control measures to reduce the impact of pandemic influenza among schoolchildren, Epidemiology and Infection 136 (8) (2007) 1035–1045. doi:10.1017/s0950268807009284. URL [https://doi.org/10.1017/s0950268807009284](https://doi.org/10.1017/s0950268807009284) [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1017/s0950268807009284&link_type=DOI) 30. [30]. J. Cui, Y. Zhang, Z. Feng, S. Guo, Y. Zhang, Influence of asymptomatic infections for the effectiveness of facemasks during pandemic influenza, Mathematical Biosciences and Engineering 16 (mbe-16-05-194) (2019) 3936. doi:[http://dx.doi.org/10.3934/mbe.2019194](http://dx.doi.org/10.3934/mbe.2019194). 31. [31]. S.-C. Chen, C.-M. Liao, Cost-effectiveness of influenza control measures: a dynamic transmission model-based analysis, Epidemiology and Infection 141 (12) (2013) 2581–2594. URL [http://www.jstor.org/stable/24475987](http://www.jstor.org/stable/24475987) 32. [32]. N. Arinaminpathy, N. Raphaely, L. Saldana, C. Hodgekiss, J. Dandridge, K. Knox, N. D. McCaarthy, Transmission and control in an institutional pandemic influenza A(H1N1) 2009 outbreak, Epidemiology and Infection 140 (6) (2011) 1102–1110. doi:10.1017/s0950268811001518. URL [https://doi.org/10.1017/s0950268811001518](https://doi.org/10.1017/s0950268811001518) [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1017/s0950268811001518&link_type=DOI) 33. [33]. B. J. Cowling, S. T. Ali, T. W. Y. Ng, T. K. Tsang, J. C. M. Li, M. W. Fong, Q. Liao, M. Y. Kwan, S. L. Lee, S. S. Chiu, J. T. Wu, P. Wu, G. M. Leung, Impact assessment of nonpharmaceutical interventions against coronavirus disease 2019 and influenza in Hong Kong: an observational study, The Lancet Public Health 5 (5) (2020) e279–e288. doi: [https://doi.org/10.1016/S2468-2667(20)30090-6](https://doi.org/10.1016/S2468-2667(20)30090-6). URL [http://www.sciencedirect.com/science/article/pii/S2468266720300906](http://www.sciencedirect.com/science/article/pii/S2468266720300906) 34. [34]. S. M. Tracht, S. Y. Del Valle, J. M. Hyman, Mathematical Modeling of the Effectiveness Facemasks in Reducing the Spread of Novel Influenza A (H1N1), PLOS ONE 5 (2) (2010) 1–12. doi:10.1371/journal.pone.0009018. URL [https://doi.org/10.1371/journal.pone.0009018](https://doi.org/10.1371/journal.pone.0009018) [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pone.0010840&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20824070&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F10%2F25%2F2020.10.21.20215418.atom) 35. [35]. A. E. Aiello, G. F. Murray, V. Perez, R. M. Coulborn, B. M. Davis, M. Uddin, D. K. Shay, S. H. Waterman, A. S. Monto, Mask use, hand hygiene, and seasonal influenza-like illness among young adults: A randomized intervention trial, The Journal of Infectious Diseases 201 (4) (2010) 491–498. doi:10.1086/650396. URL[https://doi.org/10.1086/650396](https://doi.org/10.1086/650396) [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1086/650396&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20088690&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F10%2F25%2F2020.10.21.20215418.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000273843900004&link_type=ISI) 36. [36]. S.-A. Lee, S. A. Grinshpun, T. Reponen, Respiratory performance offered by n95 respirators and surgical masks: Human subject evaluation with NaCl aerosol representing bacterial and viral particle size range, The Annals of Occupational Hygiene 52 (3) (2008) 177–185. doi: 10.1093/annhyg/men005. URL [https://doi.org/10.1093/annhyg/men005](https://doi.org/10.1093/annhyg/men005) [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/annhyg/men005&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=18326870&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F10%2F25%2F2020.10.21.20215418.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000254877300004&link_type=ISI) 37. [37]. N. C. J. Brienen, A. Timen, J. Wallinga, J. E. Van Steenbergen, P. F. M. Teunis, The effect of mask use on the spread of influenza during a pandemic, Risk Analysis 30 (8) (2010) 1210–1218. arXiv:[https://onlinelibrary.wiley.com/doi/pdf/10.1111/j.1539-6924.2010.01428.x](https://onlinelibrary.wiley.com/doi/pdf/10.1111/j.1539-6924.2010.01428.x), doi: 10.1111/j.1539-6924.2010.01428.x. URL [https://onlinelibrary.wiley.com/doi/abs/10.1111/j.1539-6924.2010.01428.x](https://onlinelibrary.wiley.com/doi/abs/10.1111/j.1539-6924.2010.01428.x) [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/j.1539-6924.2010.01428.x&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20497389&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F10%2F25%2F2020.10.21.20215418.atom) 38. [38]. S. Sim, K. Moey, N. Tan, The use of facemasks to prevent respiratory infection: a literature review in the context of the health belief model, Singapore Medical Journal 55 (3) (Mar. 2014). doi:10.11622/smedj.2014037. URL [https://doi.org/10.11622/smedj.2014037](https://doi.org/10.11622/smedj.2014037) [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.11622/smedj.2014037&link_type=DOI) 39. [39]. F. Amanat, F. Krammer, SARS-CoV-2 vaccines: Status report, Immunity 52 (4) (2020) 583–589. doi:10.1016/j.immuni.2020.03.007. URL [https://doi.org/10.1016/j.immuni.2020.03.007](https://doi.org/10.1016/j.immuni.2020.03.007) [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.immuni.2020.03.007&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32259480&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F10%2F25%2F2020.10.21.20215418.atom) 40. [40]. M. Peiris, G. M. Leung, What can we expect from first-generation COVID-19 vaccines?, The Lancet (aSep. 2020).doi:10.1016/s0140-6736(20)31976-0. URL [https://doi.org/10.1016/s0140-6736(20)31976-0](https://doi.org/10.1016/s0140-6736(20)31976-0) [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/s0140-6736(20)31976-0&link_type=DOI) 41. [41]. S. M. Bartsch, K. J. O. ÂShea, M. C. Ferguson, M. E. Bottazzi, P. T. Wedlock, U. Strych, J. A. McKinnell, S. S. Siegmund, S. N. Cox, P. J. Hotez, B. Y. Lee, Vaccine efficacy needed for a COVID-19 coronavirus vaccine to prevent or stop an epidemic as the sole intervention, American Journal of Preventive Medicine 59 (4) (2020) 493–503. doi:10.1016/j.amepre.2020.06.011. URL [https://doi.org/10.1016/j.amepre.2020.06.011](https://doi.org/10.1016/j.amepre.2020.06.011) [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.amepre.2020.06.011&link_type=DOI) 42. [42]. P. Krause, T. R. Fleming, I. Longini, A. M. Henao-Restrepo, R. Peto, N. Dean, M. Halloran, Y. Huang, T. Fleming, P. Gilbert, V. DeGruttola, H. Janes, P. Krause, I. Longini, M. Nason, R. Peto, P. Smith, A. Riveros, P. Gsell, A. Henao-Restrepo, COVID-19 vaccine trials should seek worthwhile efficacy, The Lancet 396 (10253) (2020) 741–743. doi:10.1016/s0140-6736(20)31821-3. URL [https://doi.org/10.1016/s0140-6736(20)31821-3](https://doi.org/10.1016/s0140-6736(20)31821-3) [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/s0140-6736(20)31821-3&link_type=DOI) [1]: /embed/inline-graphic-1.gif [2]: /embed/graphic-1.gif [3]: /embed/graphic-7.gif