Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Modelling the impact of interventions on the progress of the COVID-19 outbreak including age segregation

  • Jorge Rodríguez ,

    Roles Data curation, Formal analysis, Funding acquisition, Investigation, Methodology, Project administration, Resources, Software, Supervision, Validation, Visualization, Writing – original draft, Writing – review & editing

    jorge.rodriguez@ku.ac.ae

    Affiliation Department of Chemical Engineering, College of Engineering, Khalifa University, Abu Dhabi, United Arab Emirates

  • Mauricio Patón,

    Roles Data curation, Formal analysis, Investigation, Software, Validation, Visualization, Writing – original draft, Writing – review & editing

    Affiliation Department of Chemical Engineering, College of Engineering, Khalifa University, Abu Dhabi, United Arab Emirates

  • Joao M. Uratani,

    Roles Formal analysis, Investigation, Writing – original draft, Writing – review & editing

    Affiliation Department of Chemical Engineering, College of Engineering, Khalifa University, Abu Dhabi, United Arab Emirates

  • Juan M. Acuña

    Roles Conceptualization, Formal analysis, Funding acquisition, Resources, Writing – review & editing

    Affiliation Department of Epidemiology and Public Health, College of Medicine, Khalifa University, Abu Dhabi, United Arab Emirates

Abstract

In this work, a SEIR-type mathematical model of the COVID-19 outbreak was developed that describes individuals in compartments by infection stage and age group. The model assumes a close well-mixed community with no migrations. Infection rates and clinical and epidemiological information govern the transitions between stages of the disease. The impact of specific interventions (including the availability of critical care) on the outbreak time course, the number of cases and the outcome of fatalities were evaluated. Data available from the COVID-19 outbreak from Spain as of mid-May 2020 was used. Key findings in our model simulation results indicate that (i) universal social isolation measures appear effective in reducing total fatalities only if they are strict and the number of daily interpersonal contacts is reduced to very low numbers; (ii) selective isolation of only the elderly (at higher fatality risk) appears almost as effective as universal isolation in reducing total fatalities but at a possible lower economic and social impact; (iii) an increase in the number of critical care capacity directly avoids fatalities; (iv) the use of personal protective equipment (PPE) appears to be effective to dramatically reduce total fatalities when adopted extensively and to a high degree; (v) extensive random testing of the population for more complete infection recognition (accompanied by subsequent self-isolation of infected aware individuals) can dramatically reduce the total fatalities only above a high percentage threshold that may not be practically feasible.

Introduction

COVID-19 is the worst pandemic the world has seen in a century. The devastating consequences of the pandemic, both in number of fatalities as in economic harm, have brought unprecedented attention from nearly all fields. The large number of cases in a short period of time have forced governments to implement interventions to slow the spread of the virus. Such interventions consisted of the isolation of individuals (e.g. via curfews or full lockdowns, the banning of people in public events (e.g. concerts, stadiums), the enforcement to use protective equipment masks in public places, the building of field hospitals and the implementation of mass testing and contact tracing among others.

Epidemic modelling has a long history [19], as the prediction of the course and possible outcomes of a pandemic is crucial for society. The field has flourished during the pandemic period, with several mathematical models for the COVID-19 disease published from March 2020 to date. Some of these modelling efforts have evaluated the evolution of the disease spread in several countries such as in Italy [10, 11], France [12], Australia [13], the United Kingdom, the United States [14, 15], Ukraine [16], India [17] or Africa [18]. Other modelling efforts have focused on the earlier dynamics of transmission from initial cases, to the potential of the implementation of interventions to limit the disease spread, such as international travel restrictions [19, 20], contact tracing and isolation of infected individuals at onset [21], the use of masks [22], different scales of social distancing and isolation or the impact on the healthcare system [23].

Other works have focused on the estimation and/or the understanding of fundamental characteristics (i.e. potential model parameters) of the disease, such as the incubation period [24], how the transmission might occur, as well as to assess short-term [25] and long term forecasts [26].

In this work a modified SEIR-type model that evaluates intervention scenarios for the COVID-19 outbreak is presented. The model considers an approach based on the segregation by age groups and known disease stages. The model aims at retaining mechanistic meaning of all variables and parameters while capturing the relevant phenomena at play. The model aims also at contributing to the understanding by broader audiences (researchers, public health authorities, and the general public) of what to expect on the propagation of infectious diseases and how specific interventions may help. The complexity level was intentionally limited to maintain accessibility to non-experts and policy makers to comprehend the model results such that expert advice and decision making can be brought closer together to help guide interventions for immediate and longer-term needs.

Materials and methods

Model description

The model presented is based on balances of individuals, segregated by age group, transitioning between infection stages. All individuals are placed in a common single domain or compartment (e.g. a well-mixed city or town), no geographical clustering nor separation of any type is considered, nor is any form of migration in or out of the community. Large cities with ample use of public transportation are thought to be settings best described by the model.

The model also provides a direct estimation of the effective reproduction number(Rt) [27, 28] under different circumstances of individual characteristics (such as use of personal protection or awareness) as well as under population-based interventions (such as imposed social isolation). The basic reproduction number (R0) provides an estimation of the number of cases expected to be generated by one case in the population when all individuals are considered to be susceptible [28, 29]. Rt however, provides the average number of secondary infections per primary case at a given time [30] and it is used to evaluate the effectiveness of interventions [31]. The ability to estimate the Rt for different times of the outbreak (given the interventions), outbreak settings and interventions is considered to be a valuable characteristic of the model. Rt is predicted to change over time with interventions not related to increased immunity (isolation or use of personal protection equipment (PPE), as opposed to vaccination).

Model constituents.

The model solves dynamic variables or states. Every individual belongs, in addition to their age group (which she/he never leaves), to only one of the possible stages of the infection in terms of infectiousness and severity of symptoms, namely: healthy susceptible (H); asymptomatic non-infectious (NI); asymptomatic infectious (AS); symptomatic (S); in need of hospitalisation (SH); in need of critical care (SC); deceased (D) and recovered immune (R).

Definitions of the model state variables are shown in Table 1. Each variable is a state vector with the number of individuals in that stage per age group, which is defined per decade from 0–9 to 80+ year olds (9 age groups). Therefore, each state is a vector of dimensions 1x9, and the total number of states is a matrix of dimensions 8x9. Note that vector variables and parameters are represented in bold font and scalar ones in regular font.

thumbnail
Table 1. Model states in vectors (1x9) of number of individuals in each infection stage.

https://doi.org/10.1371/journal.pone.0248243.t001

Fig 1 shows the population progress through stages as described in the model. An additional schematic representation of the model approach with the population groups considered for the infection stages, rates of infection and transition between groups and showing possible interactions between population groups is shown in S5 Fig in S1 File.

thumbnail
Fig 1. Schematic representation of the population compartments by infection stages and age.

Only interactions between infectious individuals (both asymptomatic and symptomatic) and healthy susceptible ones can increase the rate of infections.

https://doi.org/10.1371/journal.pone.0248243.g001

Rates of transition between infection stages.

The transitions between stages are governed by the rates of infection and transition shown in Table 2. Notation is defined as ri,j, where the individual moves from the compartment j to the compartment i.

thumbnail
Table 2. Rates of infection and transition between states, vectors (1x9).

https://doi.org/10.1371/journal.pone.0248243.t002

The average rates of transition between states are defined such that the latest epidemiological and clinical data can be used to determine and continuously update the parameters as more knowledge of the disease emerges. These parameters include the proportion of individuals that transition to a more severe stage or recover (see Table 3) and the average times reported at each stage before transition or recovery (see Table 4).

thumbnail
Table 3. Epidemiological parameters (all in vectors per age group).

https://doi.org/10.1371/journal.pone.0248243.t003

thumbnail
Table 4. Clinical average times in each infection stage (all in vectors per age group).

https://doi.org/10.1371/journal.pone.0248243.t004

The rates of transition between stages healthy to critical (in number of individuals per day) are described in Eqs 1.A1.E. All rates are vectors per age group of dimensions (1x9). Note that point operators between vectors indicate an operation element-by-element. (1.A) (1.B) (1.C) (1.D) (1.E)

The rates of individuals recovering from the different infected stages (in number of individuals per day) are described in Eqs 2.A2.E. (all rates in vectors per age group).

(2.A)(2.B)(2.C)(2.D)(2.E)

The rate of transition from critical to deceased is the sum of that of those critical receiving intensive care (rd_scic) plus that of those critical without available care (rd_scnc) as per Eqs 3.A3.C. All critical individuals not receiving intensive care (Nsc_ncc) are assumed to become fatalities after a time (td_nc). The allocation of critical care is described below.

(3.A)(3.B)(3.C)

Fig 2 provides a representation of the transition between stages, the proportions of individuals recovering, or worsening are as per Fig 1.

thumbnail
Fig 2. Schematic representation of the transitions rates between stages.

Individuals spend in each stage an amount of time depending on their transition path towards recovery or increased severity.

https://doi.org/10.1371/journal.pone.0248243.g002

Rates of infection.

The infection of healthy susceptible individuals (H) is modelled as occurring only via their interaction with infectious either asymptomatic (AS) or symptomatic (S) individuals. Hospitalised (SH) and critical (SC) individuals are assumed not available for interactions neither are those deceased (D).

Two rates of infection of healthy susceptible individuals (in number of infections per day) are defined, one from each one of the two possible infecting groups (AS and S). The rates of infection (in vectors per age group) result from the product of (i) the fraction of interactions occurring with AS (or S) individuals among the total interactions (fias or fis) (Eqs 8.E and 8.F) times (ii) the likelihood of contagion in an interaction with AS (or S) (pi_as or pi_s, Eqs 8.A and 8.B) (per age group), (iii) the average number of daily interpersonal contacts that H individuals have (nih) and (iv) the number of H individuals themselves (per age group) (see Fig 3 and Eqs 4.A and 4.B).

(4.A)(4.B)
thumbnail
Fig 3. Schematic representation of the terms of the infection rate.

https://doi.org/10.1371/journal.pone.0248243.g003

Stage transition equations.

The dynamic variation on the number of individuals in each stage over time and per age group is governed by the population balance equations described in Eqs 5.A5.H. (all in vectors by age group).

(5.A)(5.B)(5.C)(5.D)(5.E)(5.F)(5.G)(5.H)

The state transitions as governed by these rates are represented in a matrix form in S5a Fig in S1 File.

Calculation of the effective reproduction number (Rt).

The reproduction number refers to the potential infection of susceptible individuals from infected individuals [29]. Since the model produces, for each parameter set used, a set of values for its outputs over at any given time, an instantaneous estimation of the reproduction number (Rt) is also obtained. Several parameters and variables influence the Rt including the duration of infectious stages; the likelihoods of infection per social contact as well as the percentages of individuals transitioning to more severe stages.

The effective reproduction number (Rt) during the outbreak is computed over time according to Eq 6 from the current values of the model state variables. Under this approach, infectious individuals can only infect others while they are in asymptomatic (AS) and symptomatic (S) stages. Although it has been speculated that post-symptomatic recovered individuals may be infectious for some period of time, this has not been considered in the model at this time due to lack of data. Hospitalised and critical individuals are assumed to be well isolated and also not able to infect others. The provided dynamic output of the reproduction number Rt can be used to guide and interpret the impact of interventions in terms of Rt.

Modelled infected individuals can take only three possible infectious paths, namely: (i) AS → R; (ii) AS → S → R and (iii) AS → S → SH. These paths are made of combinations of four possible infectious stage intervals in which infected individuals spend time and infect at their corresponding rate (see Table 5).

thumbnail
Table 5. Possible infectious stages intervals for Rt computation.

https://doi.org/10.1371/journal.pone.0248243.t005

The dynamic computation of Rt consists of adding the total infection contributions of every stage interval as shown in Eq 6. (6) in which, the age group weighted average rates of infection by AS and S are given as per Eqs 7.A and 7.B.

(7.A)(7.B)

Model limitations.

The model shares many of the fundamental characteristics of compartment SEIR-type models and is based on dynamic balances of individuals in compartments classified by their stage of infection and age groups only. No other differentiation within those groups is captured by the model. This characteristic allows for the model application only to single, densely populated clusters. The model has a low complexity and requires a small number of parameters that are also mechanistic and meaningful. Most of these parameters can be directly estimated from epidemiological and clinical data and are not recommended for calibration against data. Analogous to all SEIR-type models this model carries limitations since all variables and parameters refer to representative averages for each compartment of stage and age group. This may limit the model representation of the non-linear interactions that occur in reality. Phenomena like so-called super spread events or any location specific phenomena are not reproduced by these types of models. Any quantitative interpretation of results for prediction purposes should therefore be always accompanied by a critical discussion against these limitations.

Modelling interventions

The model presented above was used to evaluate different possible interventions to slow the spread of the COVID-19 outbreak. Four main interventions that have been or may be applied were evaluated namely (i) the degree of social isolation of the individuals in the population; (ii) the level of personal protection and awareness that individuals apply to protect themselves and others against contagion during interactions; (iii) the implementation of mass testing and (iv) the intensive care capacity available. These interventions can be stratified by age groups. The key parameters that define the interventions as modelled are described in Table 6.

Social isolation.

The degree of social isolation is described through its impact on the parameter (nih) (vector per age group) corresponding to the representative average number of interpersonal contacts that healthy susceptible individuals have with others per day. Different nih values can be applied to different age groups to describe age selective isolation strategies such as e.g. isolation of the elderly and/or the young alone.

Level of protection and awareness.

The level of protection and awareness describes the use of personal and protective equipment (e.g. masks). The level of protection impacts the likelihood of infection (see Eq 4) and is described through the parameters (lpah) for healthy and (lpaas and lpas) for infectious AS and S individuals (both in vectors per age group). Values of the lpa parameters can vary between 0 and 1, with 1 corresponding to the use of complete protective measures and zero to no use of any protection equipment. Analogously to the nih parameter, different values can be assigned per age group e.g. elders having a higher lpah value due to being more prone to wearing masks compared to teenagers.

The likelihoods of infection per interaction are calculated as per Eqs 8.A and 8.B. (8.A) (8.B) where lpah reflects a level of protection and awareness related to interventions and defined below; lpaasav and lpasav are scalars corresponding to the weighted averages over all age groups of the pool of AS and S with which H individuals can interact (Eqs 8.C and 8.D). (8.C) (8.D) where NasT and NsT are the total numbers of AS and S individuals of all ages respectively. The Σ symbol indicates summation across all age groups.

Awareness of infection by mass testing.

The awareness of infection after a positive test is assumed to lead to a full quarantine and removes those individuals from regular interaction with others. This awareness of infection is described through a reduction factor of interactions (rfi). Symptomatic individuals are already assumed to have a self-imposed quarantine (rfis) even if unaware of a positive test.

The fractions of individuals (without and with symptoms) aware of infection is therefore equal to the product of the fraction of total individuals randomly tested from the entire population (ptas for AS) and from the entire pool of individuals showing symptoms (pts for S) times the corresponding test sensitivity (tsns_as and tsns_s respectively). Test sensitivity might differ between each of the two groups (symptomatic and asymptomatic). In asymptomatic and mild disease individuals tested, a sensitivity of the tests of 80% was assumed.

The fractions of infectious AS and S individuals that remain in interaction with others (fias and fis) are therefore calculated as per Eqs 8.E and 8.F. Hospitalised, critical and deceased are considered excluded from the pool of interacting individuals. (8.E) (8.F) where pts is the proportion of symptomatic individuals tested (randomly) and ptas is the proportion of randomly tested non-symptomatic individuals of all types. The parameters tsns_s and tsns_as refer to the sensitivity of the tests for both groups, respectively. Table 3 shows the definitions and units of all the parameters used in the modelling of interventions.

Critical care capacity.

The impact of available critical care capacity is modelled by a specific function to allocate critically ill individuals as per the available ICU. The function allocates critically ill individuals in two possible groups, namely those admitted to ICU (Nsc_ic) and those not admitted to ICU due to lack of capacity or for other medical reasons (Nsc_ncc). At each simulation time point the allocation function is computed for the total Nsc per age group.

The function allocates ICU resources with priority to age groups with higher ICU survival rate (fr_sc) until the maximum number of intensive care units is reached leaving any remaining individuals without care, in this way Nsc_ic and Nsc_ncc are computed.

As the COVID-19 outbreak has progressed, data indicate that not all patients in critical condition have been admitted into intensive care units (ICU). Data show that many individuals with very poor prognosis, particularly those of oldest age may have never been referred to ICU due to capacity limitations or other medical humanitarian reasons. This is based on data from Spain as of May 2020 [32] showing that for individuals over 70, only a fraction of the reported fatalities previously hospitalised was ever admitted. In order to maintain consistency with the reported data as of May 2020 [32] the parameters of fd_sc and fsc_sh have been estimated such that the product of fd_sc * fsc_sh (fatality ratios over hospitalised individuals) is consistent with reported numbers for all ages irrespective of reported ICU admissions.

Results and discussion

Impact of interventions on a COVID-19 outbreak case study

A case study based on a scenario of propagation of the COVID-19 pandemic using data available from Spain as of May 2020 [32, 33] is presented below. The results obtained are intended to be interpreted qualitatively and to be contextualized to the specific setting characteristics. They serve also as a demonstration of the model potential if applied with higher confidence parameter values. Several selected scenarios were simulated aimed at illustrating the impact of different interventions.

Default reference epidemiological and clinical parameter values were extracted from different information sources on the COVID-19 outbreak as available in May 2020 [24, 32, 34, 35]. Details of values and sources are provided in S1 and S2 Tables in S1 File respectively, with indications of the level of confidence. A population with an age distribution as that of the region of Madrid (Spain) in 2019 was used [36] for the simulations.

Default reference values for interventions-related parameters (number of interactions per age group, level of protection awareness, and the isolation due to testing) were selected in order to describe a situation prior to the outbreak and without any specific intervention (see values and rationale in S3 Table in S1 File). The value for available intensive care beds per million people (capICpM) of 261 has been used by default in all case studies. The dynamic simulation results of the default outbreak scenario under no intervention is shown in S4 Fig in S1 File.

All scenarios are simulated for 365 days and evaluated in terms of (i) final total number of fatalities at outbreak termination and (ii) final number of fatalities per age group. In addition, the scenarios are presented also in terms of dynamic profiles over time for (iii) number of active cases; (iv) reproduction number; (iv) number of critical cases; (v) number of fatalities.

Intervention #1. Social isolation.

Under this scenario, the impact of the degree of social isolation (nih) was evaluated by age brackets. Four cases are shown, namely (i) universal social isolation, (ii) isolation for elders only (for individuals over 60 years old), (iii) isolation for youngsters (below 20 years old), and (iv) isolation for both elders and youngsters (only people 20 years or older and younger than 60 years old were not isolated).

The comparison between the impacts of the four types of isolations is shown in Fig 4. For the complete results of the four isolation intervention cases, please refer to the S1 File, Section VI.

thumbnail
Fig 4.

Impact of the degree of social isolation under the four different strategies (universal and per age groups) on the total number of fatalities at the end of the outbreak (left) and on the maximum values of critically ill individuals ever reached during the outbreak (right). Numbers are as percentage of the total population.

https://doi.org/10.1371/journal.pone.0248243.g004

It appears that social isolation starts to be effective in terms of significantly reducing total deaths after the number of daily contacts is placed below a threshold number (approximately two to three daily contacts in this example). Above this number of contacts per day, it does not appear to significantly modify the total final fatalities. This represents a decrease in number of contacts from 70–80%, which falls in line with the observations reported in literature that social distancing required over 80% to be effective [13].

Full details of the results shown in Fig 4 can be seen in (see S6 Fig in S1 File). The model is capable of capturing the effect of the number of interactions partly due to its description of the saturation of the healthcare capacity. S6a Fig in S1 File (middle left) illustrates the “flatten the curve” concept, now globally popular. In line with other publications [10, 11], if interactions are not decreased, the number of cases grows rapidly and explosively. S6a Fig in S1 File (middle right) shows the estimate of Rt over time. This illustrates how the level of social isolation can define the infectability and the number of cases each infected individual will infect (Rt) showing how factors such as interventions can impact Rt. The number of critical cases increases throughout time as the social interaction increases S6a Fig in S1 File (bottom left).

As shown in Fig 4, the selective social isolation of the elderly has a potentially very significant impact on final total fatalities at an almost comparable level than for universal isolation. This is consequence of the much higher mortality by the disease in the elder than in the younger (S1 Table in S1 File, Section I) [32]. This is a result with potentially significant consequences as it indicates that a sustained isolation selective only to the elderly and not to the other age groups could be applied for economic and social reasons at a small cost in number of increased total fatalities. The isolation of the young individuals (< 20 years old) produces no effect in the overall number of total fatalities. The detailed results of imposed social isolation selective to age groups is shown in S6b-S6d Fig in S1 File.

Many of the early interventions during the COVID-19 outbreak started by isolating the elderly and the young (no schools, no colleges, or universities for students), decreasing the number of interactions of the two subpopulations substantially. The isolation of these population groups together results in similar effects to that of the isolation of the elderly alone with no significant added value in isolating the young respect only the elderly as shown in Fig 4.

Intervention #2. Level of protection and awareness.

The impact of the level of self-protection and awareness was evaluated in this scenario. The parameter that describes this intervention is the level of personal protection and awareness (lpa) of the healthy and infectious population groups (lpah, lpaas and lpas) (see S3 Table in S1 File). Increases in these parameters decrease the likelihood of infection per interaction (see Eqs 8.A and 8.B) and subsequently the rates of infection (Eqs 4.A and 4.B). Fig 5 presents the main impacts of this intervention on the total final fatalities and on the maximum number of critical cases. For the full descriptive results, including the predictions of the effective reproduction number Rt, please refer to S7 Fig in S1 File, Section VII.

thumbnail
Fig 5.

Impact of the levels of protection and awareness on the total number of fatalities at the end of the outbreak (left) and on the maximum values of critically ill individuals ever reached during the outbreak (right). Numbers are as percentage of the total population.

https://doi.org/10.1371/journal.pone.0248243.g005

Fig 5 shows that high levels of protection and awareness appear as having a potentially major impact, leading to a very significant reduction on total outbreak fatalities at the highest levels of protection. This is in line of the substantial impact of the masks previously reported in other models [22, 37], however in their case they point that even low protection of masks might deem important decreases in number of cases. The height of the peak of critical cases decreases the higher the level of protection is, Fig 5 (right). In addition, if the peak does not exceed the critical care capacity (0.026% of the population), the total number of fatalities reaches much smaller values. Although moderate values for the levels of protection and awareness of 60% or lower decrease the speed of the outbreak they eventually yield the similar number of fatalities due to the exceeded ICU capacity.

Intervention #3. Awareness of infection by testing.

The impact of widespread extensive random testing was evaluated in this scenario. Fig 6 shows the comparison of the impact of testing only symptomatic individuals, testing randomly only the non-symptomatic population (therefore the same fraction of asymptomatic) and testing everyone both with and without symptoms.

thumbnail
Fig 6. Impact of the infection detection level (product of the proportion of individuals tested times the test the sensitivity, 80%) for symptomatic alone, asymptomatic alone and for both.

The impact on the total number of fatalities at the end of the outbreak (left) and on the maximum values of critically ill individuals ever reached during the outbreak (right) are shown. Numbers are as percentage of the total population.

https://doi.org/10.1371/journal.pone.0248243.g006

Only when both symptomatic and asymptomatic groups are extensively tested to levels allowing for a detection of near 90% of total infections, a meaningful impact is predicted. The current test sensitivities (around 80%) imply that even at 100% of population tested those high levels of detections required will be unachievable and a large number of infections will remain undetected. Therefore, the intervention appears to be ineffective if only applied by itself (for the complete simulation results of the impact of this intervention, please refer to S8a-S8c Fig in S1 File, Section VIII). The number of fatalities (and also number of cases) could be decreased substantially only when mass testing is combined with other interventions (social isolation, increased level of protection and/or increase in ICU beds) [11]. In addition, Fig 6 shows clearly that testing only symptomatic individuals leads to barely no impact at all. The results show that extensive testing as a sole strategy (not combined with other interventions such as social isolation, higher lpa and increase in ICU beds) appears as unfeasible to decrease the number of fatalities due to the high percentage of detections, sustained over time, required.

Intervention #4. Critical care capacity.

The impact of the availability of intensive care beds was evaluated in this scenario. The parameter that describes this intervention is the number of available intensive care beds per million population. Fig 7 illustrates the model predictions for this scenario, in terms of total final fatalities and numbers of critical cases over time. Note that once the ICU beds capacity is exceeded the critically ill patients become fatalities in one day.

thumbnail
Fig 7.

Impact of the availability of intensive care beds on the final number of fatalities total and also per age group (left) as well as for the different time course profiles of the number of critical cases (right). Numbers are in percentage of the total population of all ages.

https://doi.org/10.1371/journal.pone.0248243.g007

From Fig 7 (for more complete details, refer to S9 Fig in S1 File, Section IX), the large impact that the increase in critical care resources can have in decreasing total fatalities becomes evident. The higher the availability of critical beds, the lower number of fatalities. These appears to be in accordance with the interventions in several countries for building field hospitals [38]. The trend applies until there is no shortage of IC beds and all remaining fatalities are only the unavoidable ones. This intervention avoids those deaths that are preventable by the availability of critical care support for those that need it.

Conclusions

The impact of specific interventions on the outbreak time course, number of cases and outcome of fatalities were evaluated. Data available from Spain for the COVID-19 outbreak as of May 2020 was used. Our results on the impact and mechanisms of the various interventions indicate that:

  1. Universal social isolation measures may be effective in reducing total fatalities only if they are strict and the average number of daily social interactions is reduced below a low threshold number.
  2. Selective isolation of only the age groups most vulnerable to the disease (i.e. older than 60) appears almost as effective as the universal isolation in reducing total fatalities but at a possible much lower economic damage. The compared impacts of different social isolation interventions on the final total number of fatalities (Fig 4, left) shows that the isolation of only the elderly can achieve equivalent impact to that of all.
  3. An increase in the number of critical care beds could save a significant number of lives. Using our current parameters values, an estimate of 8 fatalities could be avoided per extra available critical care bed unit.
  4. The use of protective equipment (PPE) appears capable of reducing very significantly the total fatalities if implemented extensively and to a high degree.
  5. Extensive random testing of the population leading to isolation of the infected individuals, appears to be an ineffective intervention (if applied only by itself) due to the required (unreachable with existing test sensitivities) high percentage of infection detections and the logistic difficulty for its maintenance over time.

Any quantitative interpretation of the above results must be accompanied with a critical discussion in terms of the model limitations and its frame of application. The sensitivity analyses provided can be used to help such analysis.

Supporting information

S1 File. The Matlab® source code and Excel file containing all parameter values used as well as a non-age segregated version of the model are available at https://github.com/EnvBioProM/COVID_Model.

https://doi.org/10.1371/journal.pone.0248243.s001

(PDF)

Acknowledgments

All authors wish to thank Khalifa University and the Government of Abu Dhabi for the support.

References

  1. 1. Chowell G, Sattenspiel L, Bansal S, Viboud C. Mathematical models to characterize early epidemic growth: A review. Physics of Life Reviews. 2016;18: 66–97. pmid:27451336
  2. 2. Grassly NC, Fraser C. Mathematical models of infectious disease transmission. Nat Rev Microbiol. 2008;6: 477–487. pmid:18533288
  3. 3. Hethcote HW. The Mathematics of Infectious Diseases. SIAM Rev. 2000;42: 599–653.
  4. 4. Keeling MJ, Eames KTD. Networks and epidemic models. Journal of The Royal Society Interface. 2005;2: 295–307. pmid:16849187
  5. 5. Kermack WO, McKendrick AG. A contribution to the mathematical theory of epidemics. Proceedings of the Royal Society of London Series A, Containing Papers of a Mathematical and Physical Character. 1927;115: 700–721.
  6. 6. May RM, Anderson RM. Population biology of infectious diseases: Part II. Nature. 1979;280: 455–461. pmid:460424
  7. 7. Pastor-Satorras R, Castellano C, Van Mieghem P, Vespignani A. Epidemic processes in complex networks. Rev Mod Phys. 2015;87: 925–979.
  8. 8. Riley S. Large-Scale Spatial-Transmission Models of Infectious Disease. Science. 2007;316: 1298–1301. pmid:17540894
  9. 9. Berezovskaya F, Karev G, Song B, Castillo-Chavez C. A Simple Epidemic Model with Surprising Dynamics. Mathematical Biosciences & Engineering. 2005;2: 133. pmid:20369916
  10. 10. Gatto M, Bertuzzo E, Mari L, Miccoli S, Carraro L, Casagrandi R, et al. Spread and dynamics of the COVID-19 epidemic in Italy: Effects of emergency containment measures. Proc Natl Acad Sci USA. 2020; 202004978. pmid:32327608
  11. 11. Giordano G, Blanchini F, Bruno R, Colaneri P, Filippo AD, Matteo AD, et al. Modelling the COVID-19 epidemic and implementation of population-wide interventions in Italy. Nat Med. 2020; 1–6. pmid:31932805
  12. 12. Salje H, Kiem CT, Lefrancq N, Courtejoie N, Bosetti P, Paireau J, et al. Estimating the burden of SARS-CoV-2 in France. Science. 2020;369: 208–211. pmid:32404476
  13. 13. Chang SL, Harding N, Zachreson C, Cliff OM, Prokopenko M. Modelling transmission and control of the COVID-19 pandemic in Australia. Nature Communications. 2020;11: 5710. pmid:33177507
  14. 14. Wang L, Wang G, Gao L, Li X, Yu S, Kim M, et al. Spatiotemporal Dynamics, Nowcasting and Forecasting of COVID-19 in the United States. arXiv:200414103 [stat]. 2020 [cited 25 Nov 2020]. Available: http://arxiv.org/abs/2004.14103
  15. 15. Zou D, Wang L, Xu P, Chen J, Zhang W, Gu Q. Epidemic Model Guided Machine Learning for COVID-19 Forecasts in the United States. medRxiv. 2020; 2020.05.24.20111989.
  16. 16. Kyrychko YN, Blyuss KB, Brovchenko I. Mathematical modelling of the dynamics and containment of COVID-19 in Ukraine. Scientific Reports. 2020;10: 19662. pmid:33184338
  17. 17. Sarkar K, Khajanchi S, Nieto JJ. Modeling and forecasting the COVID-19 pandemic in India. Chaos, Solitons & Fractals. 2020;139: 110049. pmid:32834603
  18. 18. van Zandvoort K, Jarvis CI, Pearson CAB, Davies NG, Nightingale ES, Munday JD, et al. Response strategies for COVID-19 epidemics in African settings: a mathematical modelling study. BMC Medicine. 2020;18: 324. pmid:33050951
  19. 19. Chinazzi M, Davis JT, Ajelli M, Gioannini C, Litvinova M, Merler S, et al. The effect of travel restrictions on the spread of the 2019 novel coronavirus (COVID-19) outbreak. Science. 2020;368: 395–400. pmid:32144116
  20. 20. Tang B, Wang X, Li Q, Bragazzi NL, Tang S, Xiao Y, et al. Estimation of the Transmission Risk of the 2019-nCoV and Its Implication for Public Health Interventions. Journal of Clinical Medicine. 2020;9: 462. pmid:32046137
  21. 21. Hellewell J, Abbott S, Gimma A, Bosse NI, Jarvis CI, Russell TW, et al. Feasibility of controlling COVID-19 outbreaks by isolation of cases and contacts. The Lancet Global Health. 2020;8: e488–e496. pmid:32119825
  22. 22. Eikenberry SE, Mancuso M, Iboi E, Phan T, Eikenberry K, Kuang Y, et al. To mask or not to mask: Modeling the potential for face mask use by the general public to curtail the COVID-19 pandemic. Infectious Disease Modelling. 2020;5: 293–308. pmid:32355904
  23. 23. Davies NG, Kucharski AJ, Eggo RM, Gimma A, Edmunds WJ, Jombart T, et al. Effects of non-pharmaceutical interventions on COVID-19 cases, deaths, and demand for hospital services in the UK: a modelling study. The Lancet Public Health. 2020;5: e375–e385. pmid:32502389
  24. 24. Lauer SA, Grantz KH, Bi Q, Jones FK, Zheng Q, Meredith HR, et al. The Incubation Period of Coronavirus Disease 2019 (COVID-19) From Publicly Reported Confirmed Cases: Estimation and Application. Annals of Internal Medicine. 2020 [cited 29 Apr 2020]. pmid:32150748
  25. 25. Roosa K, Lee Y, Luo R, Kirpich A, Rothenberg R, Hyman JM, et al. Real-time forecasts of the COVID-19 epidemic in China from February 5th to February 24th, 2020. Infectious Disease Modelling. 2020;5: 256–263. pmid:32110742
  26. 26. Kissler SM, Tedijanto C, Goldstein E, Grad YH, Lipsitch M. Projecting the transmission dynamics of SARS-CoV-2 through the postpandemic period. Science. 2020 [cited 30 Apr 2020]. pmid:32291278
  27. 27. De Serres G, Gay NJ, Farrington CP. Epidemiology of Transmissible Diseases after Elimination. Am J Epidemiol. 2000;151: 1039–1048. pmid:10873127
  28. 28. Adam D. A guide to R—the pandemic’s misunderstood metric. Nature. 2020;583: 346–348. pmid:32620883
  29. 29. Delamater PL, Street EJ, Leslie TF, Yang YT, Jacobsen KH. Complexity of the Basic Reproduction Number (R0). Emerging Infectious Diseases journal. 2019;25: 1–4. pmid:30560777
  30. 30. Nishiura H, Chowell G. The Effective Reproduction Number as a Prelude to Statistical Estimation of Time-Dependent Epidemic Trends. In: Chowell G, Hyman JM, Bettencourt LMA, Castillo-Chavez C, editors. Mathematical and Statistical Estimation Approaches in Epidemiology. Dordrecht: Springer Netherlands; 2009. pp. 103–121. https://doi.org/10.1007/978-90-481-2313-1_5
  31. 31. Gostic KM, McGough L, Baskerville EB, Abbott S, Joshi K, Tedijanto C, et al. Practical considerations for measuring the effective reproductive number, Rt. PLoS Comput Biol. 2020;16: e1008409. pmid:33301457
  32. 32. ISCIII. Situación de COVID-19 o Coronavirus en España. 30 Apr 2020 [cited 10 May 2020]. Available: https://covid19.isciii.es/
  33. 33. Pollán M, Pérez-Gómez B, Pastor-Barriuso R, Oteo J, Hernán MA, Pérez-Olmeda M, et al. Prevalence of SARS-CoV-2 in Spain (ENE-COVID): a nationwide, population-based seroepidemiological study. The Lancet. 2020;396: 535–544. pmid:32645347
  34. 34. Nishiura H, Kobayashi T, Miyama T, Suzuki A, Jung S, Hayashi K, et al. Estimation of the asymptomatic ratio of novel coronavirus infections (COVID-19). International Journal of Infectious Diseases. 2020;94: 154–155. pmid:32179137
  35. 35. Zhou F, Yu T, Du R, Fan G, Liu Y, Liu Z, et al. Clinical course and risk factors for mortality of adult inpatients with COVID-19 in Wuhan, China: a retrospective cohort study. The Lancet. 2020;395: 1054–1062. pmid:32171076
  36. 36. INE. Instituto Nacional de Estadística. In: INE [Internet]. 2020 [cited 4 Nov 2020]. Available: https://www.ine.es/
  37. 37. Ngonghala CN, Iboi E, Eikenberry S, Scotch M, MacIntyre CR, Bonds MH, et al. Mathematical assessment of the impact of non-pharmaceutical interventions on curtailing the 2019 novel Coronavirus. Mathematical Biosciences. 2020;325: 108364. pmid:32360770
  38. 38. Wallis N, Gust C, Porter E, Gilchrist N, Amaral A. Implementation of field hospital pharmacy services during the COVID-19 pandemic. Am J Health Syst Pharm. 2020;77: 1547–1551. pmid:32539128