key: cord-1005134-esqq615d authors: Ortega-Quijano, Daniel; Ortega-Quijano, Noe title: Impact of age-selective vs non-selective physical-distancing measures against coronavirus disease 2019: a mathematical modelling study date: 2021-03-12 journal: Int J Epidemiol DOI: 10.1093/ije/dyab043 sha: e82c088f7610010df60817927a90698bea30b514 doc_id: 1005134 cord_uid: esqq615d BACKGROUND: There is a real possibility of successive COVID-19-epidemic waves with devastating consequences. In this context, it has become mandatory to design age-selective measures aimed at achieving an optimal balance between protecting public health and maintaining a viable economic activity. METHODS: We programmed a Susceptible, Exposed, Infected, Removed (SEIR) model in order to introduce epidemiologically relevant age classes into the outbreak-dynamics analysis. The model was fitted to the official death toll and calculated age distribution of deaths in Wuhan using a constrained linear least-squares algorithm. Subsequently, we used synthetic location-specific and age-structured contact matrices to quantify the effect of age-selective interventions both on mortality and on economic activity in Wuhan. For this purpose, we simulated four different scenarios ranging from an absence of measures to age-selective interventions with stronger physical-distancing measures for older individuals. RESULTS: An age-selective strategy could reduce the death toll by >30% compared with the non-selective measures applied during Wuhan’s lockdown for the same workforce. Moreover, an alternative age-selective strategy could allow a 5-fold increase in the population working on site without a detrimental impact on the death toll compared with the Wuhan scenario. CONCLUSION: Our results suggest that age-selective-distancing measures focused on the older population could have achieved a better balance between COVID-19 mortality and economic activity during the first COVID-19 outbreak in Wuhan. However, the implications of this need to be interpreted along with considerations of the practical feasibility and potential wider benefits and drawbacks of such a strategy. On 23 January 2020, with <400 accumulated confirmed diagnoses of coronavirus disease (COVID- 19) , an unprecedented lockdown that lasted for >10 weeks was established in the city of Wuhan, China. 1 Control measures enacted by local and national authorities included a travel ban; non-selective social-distancing policies; the obligation to wear masks; the extension of the Chinese Lunar New Year holiday; and the closure of schools, universities, factories, government offices and, ultimately, of all non-essential services. 2, 3 As a result of these measures, on 8 April, with just one new case of local transmission reported over the preceding 2 weeks, local authorities began to ease restrictions in Wuhan. 4 Meanwhile, the disease spread globally to become a pandemic and >50 000 people were infected and another 2500 died in Wuhan according to local authorities. 1, 5 Non-selective physical-distancing measures have substantially halted transmission of COVID-19 worldwide, giving extra time to the unprecedented global effort that is underway to find treatments and vaccines in a context of successive epidemic waves. [6] [7] [8] However, as various mathematical models predicted, uncontrolled community transmission might make it necessary to alternate periods of relaxation of measures with repeated periods of lockdown in the near and mid-term, with potentially devastating consequences for the public health and economy of countries. 9 In fact, Israel was the world's first country to reimpose lockdown in September 2020. 10 To put the economic consequences into perspective, the monthly losses during Wuhan's closure have been estimated at >52% of the 2019 first-quarter municipal gross domestic product. 11 The purpose of this study is to present an alternative approach to the design of social-distancing measures. In particular, we investigate the impact of age-selective physical-distancing measures focused on older age groups on the cumulative number of COVID-19-related deaths and the number of people in the active workforce in the setting of Wuhan. It is an exercise of a mathematical nature and, although the practical limitations are outlined in the discussion, the specific design and application of the measures are beyond the scope of this study. To simulate Wuhan's outbreak and the effect of alternative lockdown scenarios, we programmed a deterministic agedependent and location-specific Susceptible, Exposed, Infected, Removed (SEIR) model in order to introduce epidemiologically relevant age classes into the outbreak-dynamics analysis. Initial parameters, including the Wuhan-population structure and Prem's synthetic age-dependent and locationspecific contact matrices for modelling social interactions among the population in Wuhan, 8 were obtained from the literature (Table 1) . All the data used to calculate these parameters were ultimately extracted from official publications from the National Health Commission, China CDC, World Bank, United Nations Statistics Division, International Labor Organization, POLYMOD study and other publicly available sources. 1, 18, 19 Procedures Compartments In deterministic compartmental models, individuals are divided into classes solely based on their infection status. • Age is an independent risk factor for mortality in patients with coronavirus disease (COVID-19). • Non-selective physical-distancing measures enacted by the Chinese authorities effectively contributed to controlling Wuhan's outbreak. • No published article has modelled and compared the effect of age-selective vs non-selective physical-distancing interventions against COVID-19 on both mortality and economic workforce, nor proposed age-structured control measures based on both aspects. • We assessed the effect of age-selective physical-distancing interventions compared with non-selective interventions on mortality and workforce reduction in Wuhan. • We found that age-selective social restrictions specifically balanced for overprotecting vulnerable age groups would likely reduce the total number of deaths while allowing a higher percentage of the workforce to remain at their workplace. Susceptible (S) individuals have not been exposed to the pathogen yet and are not immunized against it, so they are prone to be infected. Susceptible individuals become Exposed (E) when they are infected but not yet infectious, and they eventually become Infected (I) when they can further transmit the infection. Finally, infected individuals are Removed (R) when the infectious period is finished. Under the close-system assumption ignoring population dynamics, the total population is N ¼ S þ E þ I þ R. In view of the strongly asymmetric impact of COVID-19 on different age groups, with a much higher mortality for the elderly compared with the younger population, this work incorporates age groups into the SEIR model described above in order to include the age-structured results into the simulations outcome. 17 We have grouped the population into 10-year bands except for the last group, which includes all individuals >70 years old, which results in eight age categories. This decision arises from the convenience to harmonize our agedependent SEIR model with the age categories used by Chinese authorities to disseminate the official mortality data of the COVID-19 outbreak. 1, 19 As a result, the number of individuals and the total population in each compartment are no longer quantified by a single number, but by a column vector of eight values corresponding to the age groups, where index i denotes the age group. Our SEIR model and the transitions of individuals across the model compartments are presented in Figure 1 . The number of days for which exposed individuals remain in the incubation period before becoming infected follows an Erlang-probability density function of shape factor 2 and the rate of becoming infected r ¼ 1=d L , where d L is the mean incubation period. Likewise, the infectious period is modelled by an Erlang distribution of shape factor 2 and mean infectious period d I , which determines the removal rate c ¼ 1=d I that describes the transition of individuals from the Infected to the Removed compartment. The interaction between susceptible and infectious individuals generates exposed individuals at a rate denoted as k and better known as force of infection, defined as the rate at which susceptible individuals become infected but not yet infectious. In conventional age-independent SEIR models, the per-capita infection rate is k ¼ bI=N, which means that it is governed by the transmission rate b and by the infection prevalence expressed as the fraction of infectious individuals I=N, where N is the total population. Importantly, the transmission rate encompasses both the underlaying social-contact dynamics and the probability of infection transmission. 20 The rate of becoming infected (r) and the removal rate (c) are considered to be age-independent. Regarding the force of infection (k), the transmission rate can now be written as b ¼ sC, where s is the probability of infection transmission and C is the contact matrix. Each element C i;j of the contact matrix in a given scenario quantifies the average number of daily contacts made by an individual in age group i with an individual in age group j. 21 The basic reproduction number R 0 is defined as the average number of secondary cases produced by a primary infectious individual in a susceptible population. 20 It is calculated as the ratio between the transmission rate and the recovery rate, which can in turn be expressed as: where C is the average number of contacts of the overall population, which can be readily calculated by the weightaveraging the contact matrix C by the age-structured population vectorÑ. The median basic reproduction number in Wuhan before the lockdown has been previously estimated. 13 As for the infection period, it can easily be derived from clinical data and approximate values are known from the beginning of the outbreak (see Table 1 ). Knowing all these values, the calculation of the mean probability of transmission prior to lockdown yields a value of 7.08%. Assuming that the probability of transmission per contact is a biologically stable parameter, we introduced it into the model as a constant. With the aim of quantifying the number of deceases produced by the COVID-19 outbreak, as well as the effect of the physical-distancing measures proposed in this work, we used age-specific case-fatality ratios adjusted by the method proposed by Verity et al. 17 As a result, the SEIR model describing the number of individuals in each compartment for age group i is described by the following system of coupled ordinary differential equations: All parameters intervening in the equations and their description are listed in Table 1 . The subscript of the contact-matrix location has been omitted for the sake of simplicity. Due to the fact that COVID-19 is a new disease, it is assumed that the entire Wuhan population was initially susceptible. Following Kucharski and colleagues, it is assumed that the outbreak was generated by a single infectious individual on 22 November 2019. 16 The final dates for each scenario are specified in Table 1 and the end-date criteria are detailed in the Supplementary Material, pp. 8-9, 14, 16 and 17-18, available as Supplementary data at IJE online. The initial infectious individual is assigned to age group 5 (40-49 years) in our model, since it is the group that includes the median age of the first 41 reported cases that were reported in Wuhan. 22 The system of ordinary differential equations governing the SEIR model is numerically solved by the Runge-Kutta method using the Dormand-Price pair. Several scenarios have been considered. In order to gain an insight into the outbreak dynamics and to provide a meaningful approach to the impact of the dominant factors on the efficiency of the physical-distancing measures, we grouped the location-specific contact matrices into two ( Figure 2 ): one for work and another one for the rest, the latter encompassing Prem's school, home and others. 8 The next paragraphs provide an overview of each simulated scenario. A detailed description can be found in the Supplementary Material, pp. 8-18, available as Supplementary data at IJE online. Scenario 1 (no physical distancing) models the theoretical evolution of a COVID-19 outbreak in Wuhan without any control measure, taking into account the school Winter Break from 15 January to 10 February 2020 and the Lunar New Year holidays from 24 to 30 January 2020. 23, 24 Scenario 2 (Wuhan's lockdown) simulates the effect on population mixing and the mortality of non-selective socialdistancing interventions taken in Wuhan. We assumed that 5% of the workforce remained at their workplace throughout the lockdown to provide essential services, so the workcontact matrix was weighted accordingly. Regarding the contact matrix for the rest, we used a constrained linear leastsquares algorithm using the minimum norm criterion to fit the estimated total number of deceased individuals to the official death toll ( Figure 3 and Supplementary Material, p. 6, available as Supplementary data at IJE online). This method provides a plausible contact matrix that best fits the Wuhan situation. This scenario will be used as the baseline for comparatively analysing the outcome of the Wuhan lockdown with the proposed alternative physical-distancing measures. Scenario 3 (public-health lockdown) introduced additional selective, age-specific physical-distancing measures consisting of an additional 30% reduction in the number of contacts for people >60 years old with the aim of minimizing the number of deaths compared with the Wuhan lockdown (Scenario 2) while keeping a similar workforce. Finally, Scenario 4 (health-economy trade-off lockdown) exploited the efficiency of the 30% reduction in the number of contacts for people >60 years old of Scenario 3 for balancing the trade-off between public health and economic activity, and quantified the additional workforce that could be maintained at their workplaces while keeping the same death toll as in Wuhan's scenario. All analyses were done with R software (version 3.6.2) and Matlab (version R2020a Update 2). A sensitivity analysis on the key assumptions can be found in the Supplementary Material, pp. 18-20, available as Supplementary data at IJE online. We found in Scenarios 3 and 4 that age-selective physicaldistancing measures preferentially targeting the older population can have an impact on the total number of deaths and make it possible to design alternative control measures with a tailored balance between public health and economic activity in order to minimize the overall stress produced by the outbreak on society as a whole. In Figures 4 and 5 , we represent the daily and cumulative infected individuals per age group, the aggregated daily and cumulative infections for all age groups and the cumulative deaths per age group between January and April/May 2020 for the Wuhan lockdown (Scenario 2) and alternative social-distancing (Scenarios 1, 3 and 4). Our deterministic age-structured SEIR model enables epidemicdynamics analyses in two different axes. The first axis represents the epidemic curve as a function of time for either a specific age group or the population as a whole, whereas the second axis provides information about the epidemic distribution as a function of age at a given time, which provides an insight into the outbreak dynamics and its impact on different population age groups. The school Winter Break and the Lunar New Year holidays alone (Scenario 1) would result in the smallest reduction in the incidence and number of deaths of all the scenarios considered (first scenario, Figure 4A -E). By contrast, Wuhan's lockdown enacted by the Chinese authorities (Scenario 2) effectively reduced the number of contacts ( Figure 2H ), the incidence, the number of deaths and the size of the outbreak compared with the 'no physical distancing' scenario ( Figures 4F-J and 6-figures 6A-B) . We found that non-selective measures in Scenario 2 apparently decreased contacts more effectively in the intermediate age groups than in the extreme age groups ( Figure 2H ) and that contacts were reduced by !80% to produce a number of deaths and an epidemic curve similar to that of Wuhan ( Figure 3) . Additional age-selective physical-distancing interventions targeting vulnerable age groups, modelled in Scenario 3 as an additional restriction of 30% further reduction in the number of contacts for people >60 years old in the rest of locations, achieved a reduction of 33% in the death toll in the simulated period ( Figure 5A-E) . The importance lies in the fact that this reduction in contacts represents only 5% of the total, but its impact on the number of deaths is disproportionate due to the high mortality for the elderly. Scenario 4 revealed that additional specific interventions would have allowed an increase in the workforce returning to work by 800% in the age group of 20 years, 600% in the age group 30 years, 400% in the age group 40 years, 100% in the age group 50 years and 100% in the age group 60 years without an increase in the number of accumulated deaths compared with Scenario 2. Although this increase in workforce would have prolonged the duration of the epidemic by roughly an additional month ( Figure 5 , date axes), Figure 6 shows that the effective workforce is 400% relative to Scenario 2 even taking this factor into account. These findings have been summarized and plotted in Figure 6D , which specifically shows that, on the one hand, measures in Scenario 3 would have resulted in the same number of people at work but a lower death toll compared with the Wuhan lockdown (Scenario 2), whereas, on the other hand, Scenario 4 would have had a comparable death toll with Wuhan lockdown but a large increase in the workforce. Changing initial parameters in the sensitivity analysis had an impact on several aspects such as the size of the epidemic, but the overall trends were maintained and the age-selective measures of Scenarios 3 and 4 were shown to have beneficial effects independently of the parameter variations (Supplementary Material, pp. 18-20, available as Supplementary data at IJE online). We have presented a deterministic age-structured SEIR model that enables the investigation of alternative control measures for effectively protecting the most vulnerable age groups while keeping a balanced active workforce. The relevance of this approach lies in the fact that one of the main prognostic factors regarding hospitalization and mortality of COVID-19 is age. 22 Compartmentalized evaluation of the outbreak impact in each age group allowed us to design and quantify the effects of hypothetical age-selective physical-distancing interventions (Scenarios 3 and 4) on the number of deaths and workforce compared with the Wuhan lockdown (Scenario 2) between January and April/ May 2020. In order to exemplify the possibilities of the model, we decided to use the Wuhan setting (Scenario 2) and then simulate the effect of different social-mixing-pattern interventions on the outbreak dynamics and workforce in a defined period of time. 1, 17 From a public-health point of view, the ideal objective of interventions would be to achieve the lowest number of infected patients and no deaths. We decided to fit our model to the official death toll, as it was considered more stable and reliable than the number of confirmed cases. 25, 26 Similarly, for a given number of cumulative deaths, the percentage of the workforce that could return to work as a result of the establishment of age-selective social-mixing measures directed at vulnerable age groups was considered a relevant measure of the economic efficacy of interventions. Wuhan's lockdown has been the model followed by many countries to control COVID-19 outbreaks. 27, 28 Other countries have chosen to carry out massive tests and established measures with less economic impact. 29, 30 However, for many reasons, not all countries have these capabilities. 31 Under the context of successive epidemic waves and without widespread access to vaccination to date, governments need to find a suitable balance and adopt measures adjusted to the health and economic reality of their country. We have found in Scenario 3 that an additional 30% reduction in the number of contacts in people >60 years old would reduce by >30% the accumulated number of deaths at the end of the simulation period. Similarly, these measures would allow governments to increase the population at work by >400% while containing the impact on public health. According to our calculations, the non-selective measures of physical distance in Wuhan caused an 80% decrease in contacts and were precisely less effective in extreme age groups, which is consistent with the available evidence in China and the UK. 32, 33 In fact, attempting to shield the elderly may have important mental-and physical-health consequences. 34 National health systems are not designed to combat pandemics. In an extreme situation, the relevance of the opportunity cost becomes palpable and the lowest opportunity cost in terms of saved lives is achieved when the available resources are allocated primarily to the most vulnerable groups. It is essential to design protocols with the best available evidence to guarantee individual protection and asepsis in interactions within age groups and vulnerable groups. 35 In the Wuhan lockdown (Scenario 2), our model accurately predicts the timescale of the outbreak, the number of deaths and their age distribution. 1, 19 However, we have found a higher number of infected individuals than reported by the Chinese authorities. 1 This issue has been previously highlighted by other authors and may be affected by changes in the definition of cases, health-system saturation, undetected subclinical or asymptomatic cases and lack of diagnostic tests. 16, 26 The study has several important limitations. On the one hand, with regard to practical issues, it is worth highlighting that this model focuses on a specific location (Wuhan) and does not explore the restrictions necessary to curtail the geographic spread of the virus. There are also social difficulties arising from keeping the elderly separated from the rest of the population for long time periods (i.e. difficulties in childcare if parents return to the workforce but schools remain closed), in addition to the direct negative medical consequences discussed above. On the other hand, there are additional technical limitations. First, despite the biological plausibility of the parameters we used, they may undergo changes as new evidence accumulates. Second, as it happened in the models on which we based our analyses, we assume that the latency period is equivalent to the incubation period. 16 Third, we assumed no heterogeneity in susceptibility and transmissibility between children and adults, although there are conflicting results regarding this issue. 32, 36 We also did not explicitly include individuallevel heterogeneity in contact, clustering of household transmission and nosocomial transmission. Fourth, we did not model the limited capacity of the health system/intensive care units (ICUs), as the official data from the Chinese government did not record daily or global ICU admissions/ COVID-case severity in Wuhan, and it should be kept in mind that long-term effects beyond May 2020 were not investigated. Finally, contact matrices were synthetic and, even though the official death toll is probably more stable and reliable than the official number of confirmed cases, it may also underestimate the actual number of deaths. As a result, there is a high level of uncertainty in the parameter values (including contact matrices), deaths data and the chosen model structure. Therefore, we carried out a sensitivity analysis for different values of R 0 , the start date of the outbreak, the number of initial patients and the ages of the first infected patients. Although contact matrices are synthetic, the constrained linear least-squares method used to fit the simulate scenarios to the available data yields more rigorous and realistic results than simply multiplying contact matrices by a fixed coefficient and provides added value over previous COVID-19 studies modelling socialdistancing or age-selective interventions. [7] [8] [9] 37 To conclude, our results suggest that age-selective interventions could achieve either a reduction on the death toll or an increase in the population working on site without a detrimental impact on the death toll for a comparable time period. However, the practical feasibility of such measures as well as their long-term effects remain to be investigated. Supplementary data are available at IJE online. None. National Health Commission of the People's Republic of China. Updates on the Epidemic The State Council of the People's Republic of China. The Announcement from Wuhan's Headquarter on the Novel Coronavirus Prevention and Control The State Council of the People's Republic of China. The State Council's announcement on extending the Lunar New Year Holiday in 2020 Wuhan reopens after 76-day lockdown speeches/detail/who-director-general-s-opening-remarks-at-the-me dia-briefing-on-covid Impact assessment of nonpharmaceutical interventions against coronavirus disease 2019 and influenza in Hong Kong: an observational study First-wave COVID-19 transmissibility and severity in China outside Hubei after control measures, and second-wave scenario planning: a modelling impact assessment The effect of control strategies to reduce social mixing on outcomes of the COVID-19 epidemic in Wuhan, China: a modelling study Impact of non-pharmaceutical interventions (NPIs) to reduce COVID-19 mortality and healthcare demand The government of Israel. Ministry of Health guidelines for lockdown Assessment of monthly economic losses in Wuhan under the lockdown against COVID-19 Early transmissibility assessment of a novel coronavirus in Wuhan Early transmission dynamics in Wuhan, China, of novel coronavirus-infected pneumonia Transmission dynamics of 2019 novel coronavirus (2019-nCoV) Early dynamics of transmission and control of COVID-19: a mathematical modelling study. The Lancet Infectious Diseases Estimates of the severity of coronavirus disease 2019: a model-based analysis Projecting social contact matrices in 152 countries using contact surveys and demographic data The Novel Coronavirus Pneumonia Emergency Response Epidemiology Team. The epidemiological characteristics of an outbreak of 2019 novel coronavirus diseases (COVID-19)-China Modeling Infectious Diseases in Humans and Animals, 1st edn Social contact patterns and control strategies for influenza in the elderly Clinical features of patients infected with 2019 novel coronavirus in Wuhan, China Curriculum Plan for the Primary and Middle School Kindergartens in the 2019-2020 School Year The State Council of the People's Republic of China. Notice of partial holiday arrangements Wrong but useful-what Covid-19 epidemiologic models can and cannot tell us Effect of changing case definitions for COVID-19 on the epidemic curve and transmission parameters in mainland China: a modelling study COVID-19 in Italy: momentous decisions and many uncertainties Covid-19-the law and limits of quarantine Spread of SARS-CoV-2 in the Icelandic population Coronavirus cases have dropped sharply in South Korea The South African Response to the Pandemic Changes in contact patterns shape the dynamics of the COVID-19 outbreak in China CMMID COVID-19 working group et al. Quantifying the impact of physical distance measures on the transmission of COVID-19 in the UK COVID-19 and the consequences of isolating the elderly Physical distancing, face masks, and eye protection to prevent person-to-person transmission of SARS-CoV-2 and COVID-19: a systematic review and meta-analysis. The Lancet Epidemiology and transmission of COVID-19 in 391 cases and 1286 of their close contacts in Shenzhen, China: a retrospective cohort study Effects of nonpharmaceutical interventions on COVID-19 cases, deaths, and demand for hospital services in the UK: a modelling study The data underlying this article are available in GitHub repository IJE_2021_Covid19AgeStructuredSEIR at https://github.com/nortega quijano/IJE_2021_Covid19AgeStructuredSEIR and can be publicly accessed. Ethics approval was not required for this mathematicalmodelling study. None declared.