Staggered release policies for COVID-19 control: Costs and benefits of relaxing restrictions by age and risk.

Lockdown and social distancing restrictions have been widely used as part of policy efforts aimed at controlling the ongoing COVID-19 pandemic. Since these restrictions have a negative impact on the economy, there exists a strong incentive to relax these policies while protecting public health. Using a modified SEIR epidemiological model, this paper explores the costs and benefits associated with the sequential release of specific groups based on age and risk from lockdown and social distancing measures. The results in this paper suggest that properly designed staggered-release policies can do better than simultaneous-release policies in terms of protecting the most vulnerable members of a population, reducing health risks overall, and increasing economic activity.

Introduction

Between March 19 and April 7, various lockdown and social distancing measures were issued around the United States. Most of these restrictions were initially slated to last 30 days, after which the states reevaluated and decided on different paths forward. Since then, most states have begun some form of reopening process. While these policy choices vary from state to state, all of these restrictions apply to every member of the population uniformly. This simultaneous-release methodology ignores the importance of differential impacts of disease on various subsets of the population.

Severity and survival rates for COVID-19 infections vary significantly. They are a function of age and various comorbidities [1]. In this paper, we use documented infection–response variation to COVID-19 to design modified lockdown and social distancing restrictions that reduce overall death rates while increasing economic activity. Staggered-release policies are compared to the benchmark policy of simultaneously releasing all groups at a predetermined time. Policies are specifically evaluated on the reduction of overall death rates from the resulting projected outbreak.

Reported data show that infected individuals over the age of 65 face a much higher case mortality than individuals under the age of 44, while those between the ages 45 and 65 have an intermediate case mortality ([2], [3], see also Fig. 1). Thus far, lockdown and social distancing restrictions have kept the overall number of deaths much lower than they would have been otherwise, particularly among the elderly. However, there are deleterious effects on the economy as well. Since restrictions were put in place in America, visits to commercial venues are down two-thirds [4], a decrease in activity that has led to small business closure and layoff rates of around 50% in the Mid-Atlantic states [5]. Shutdown sectors represent over 20% of all US payroll employment [6] and the burden of these job losses has fallen primarily on the poor [7]. All told, the economic cost of closing non-essential businesses could total nearly $10,000 per household per quarter [8], which may cause persistent harm in the form of lower output and employment and associated higher overall morbidity and mortality even after the shutdown ends [9].

Fig. 1

COVID-19 data for New York City showing age-dependent numbers of cases (top), hospitalizations (middle), and deaths (bottom) per 100,000 people (see [2]. Accessed June 9, 2020.)

While saving lives must be the primary concern, judiciously increasing economic activity is an important secondary goal as well. Perhaps the most important result from the analysis in this paper is tied to the identification of policies that can responsibly mitigate the economic effects of a sustained shutdown without jeopardizing the most vulnerable members of the population. This paper analyzes the costs and benefits of policies that are predicated upon the timely release of younger individuals who face lower risk from COVID-19 infections from social distancing restrictions before the release of older individuals. Compared to a policy that simultaneously releases all sub-populations, carefully planned and executed staggered-release polices are likely to lead to:

Fewer infections over the entire duration of the outbreak

Lower disease mortality rate among the elderly

Fewer total disease deaths across the whole population

Increased economic activity at an earlier date

The intuition behind the merit of staggered-release policies is closely related to the ‘flattening the curve’ narrative that motivated initial public support for restrictions in the first place [1], [10]. The worst case scenario in terms of public health risk and lives lost transpires when infection rates are allowed to grow unchecked in a manner which leads to too many severe infections at once, overwhelming health care systems. Because these severe infections are much more likely among the elderly, reducing and managing elderly infection rates is a useful avenue toward improving health outcomes overall.

At the time of writing, it is believed that severe restrictions, even if maintained for an extended period of time, would be unlikely to prevent a second peak of infection since the population would not have built up significant levels of immunity before economic activity is resumed. It is strongly believed that releasing everybody at once will result in the population experiencing dramatic increases in infections, and we believe this risk may be well addressed through the usage of staggered-release policies.

This paper is organized as follows: Section 2 introduces the 3-group model and presents the setup of the benchmark and staggered-release policies. Analyses and results are presented in Section 3. Section 4 discusses the findings.

The three-group model and sequential release policies

Simple SIR and SEIR types of epidemic models have been used in most current studies of COVID-19 dynamics (see, for example, [11], [12], [13], [14], [15], [16], [17], [18], [19], [20]). The model introduced in this section is based on the standard Susceptible–Exposed–Infectious–Removed (SEIR) model with the incorporation of asymptomatic infections, disease-induced deaths, hospitalizations, and preferential mixing between different age groups. These modifications enable a flexible description of COVID-19 disease transmission dynamics within a model that accounts for special features associated with age-dependent groups. Since we are modeling a single outbreak, other considerations such as aging, migration, births, and unrelated deaths are ignored.

The population is divided into three groups: group 1 consisting of healthier younger people aged under 44 who have no medical conditions which may increase their risks to disease death, group 2 include people aged between 45 and 65 without health conditions, and group 3 consisting of the vulnerable members of the population including the elderly (65 years old and over) and people with higher risks than those in groups 1 and 2. These three sub-populations are labeled by and will be referred to as Group (or simply G1, G2, and G3). This separation of age groups is motivated by the age-specific COVID-19 data for New York presented in Fig. 1, which shows age-dependent numbers of cases, hospitalizations, and deaths per 100,000 people (see [2]), with the consideration that vulnerable people in the younger groups will be included in the high-risk group 3.

Based on available information about COVID-19, there is only a very small proportion of younger people who might have a higher risk to COVID-19 infections and deaths. Based on the US demographic data, there are about 13% of people above the age of 65, and about 25% of people between ages 45–64. Thus, the majority population is in the age group 0–44. Moving a very small proportion of the vulnerable young people into group 3 will not change much of the population size for Group 1. Each sub-population is further sub-divided into seven epidemiological classes: susceptible , exposed (), infectious but asymptomatic (), infectious and symptomatic (), hospitalized (), recovered (), and dead due to disease (). The total population is , where , with .

Let denote the per capita rate of progression to the infectious state ( represents the mean latent period), and denote the per-capita recovery rates ( and are the mean infectious periods), denote the rate of transition from to , denote the rate at which hospitalized individuals leave the class with proportion recovered and dead, and denote the disease related death rate for individuals in the class. Among the infectious people, assume that proportions and are symptomatic and asymptomatic, respectively. Asymptomatic and hospitalized individuals can also transmit the disease but possibly at lower rates than symptomatic individuals, which are denoted by the factors and , respectively. The parameters , and are assumed to be independent. Other assumptions on the dependence of these parameters might also be considered, but they would not affect the qualitative conclusions of this study (see Feng et al. [21], for more detailed discussions about various underlying biological assumptions on this topic in the context of Ebola models). A disease transmission diagram for each sub-group is depicted in Fig. 2.

Fig. 2

Depiction of disease transmission process. The infection rate includes transmissions from infectious people in all sub-groups. The sub-groups are connected through a mixing function given in Eq. (3).

The model is given by the following system of differential equations: where denotes the force of infection (FOI), the generator of new cases of infection among susceptible individuals in group . The functional form of the FOI for group is given by where describes the mixing among the groups and defined as the proportion of the th sub-group’s contacts that is with members of the th group. We will adopt the commonly used preferential mixing ([22] and later extended by [23]), in which case, the elements of have the following form: where describes the preference level of group and is the Kronecker delta (1 when and 0 otherwise). The function in (3) satisfies the required constraints for mixing functions (see [24]). In our previous study, we estimated the mixing parameters including the age-dependent contact rates and preference levels using the PolyMod data and other observations [25], which suggests that younger people tent to have higher contact rates. These results help inform our choice of model parameters in this study.

In the FOI function (2), is the per capita effective contact rate, i.e., contacts that can lead to infection [26]. In the absence of any intervention (and no hospitalization), the group-specific basic reproduction number (i.e., without interactions with other groups) for group is Thus, we can fix to get the baseline value for the effective contact rate for group : Note that the effective contact rate is often expressed as a product of the contact rate with the probability of infection on contact (e.g., ). If we assume to be the same for all then it can be canceled in the function in (3). Thus, the balance condition for mixing still holds. For ease of reference, we will refer to simply as a contact rate.

Lockdown and social distancing restrictions are introduced to reduce the normal contact rates. Here, it is assumed that is a function of time and that their values are influenced by the selected policy. Specifically, we model them via step functions defined as follows: where , , and represent the times when policies may change. In the case of (5), is the time when the initial restriction on everyone begins, and are the numbers of days from when Group 1 and Group 2 may be released from restrictions, respectively, and is the day from when the restriction for all groups is reduced to a residual level until the outbreak is over. The factor 0.95 on the severity of restrictions is meant to reflect the fact that complete isolation of every single member of a population is infeasible and likely inhumane.

The parameters in (5) represent the reductions of contacts for group (or severity of restrictions imposed on group from business as usual (0) to no contact (1)) during the period of time , with and (). We remark that the parameters also allow us to take into consideration that a proportion of Group may not change their contacts for any reason. For example, is a weighted average of restriction levels of people in Group 1 who will or will not change their contacts during the time period . To simplify notation, let denote the reduction of contacts during the initial restrictions, i.e., , and denote the residual restriction after the strict policies are lifted for all groups, i.e., for . Throughout this paper, the policy corresponding to for and for , which is the case where all groups are held under the same severity of restrictions for the period and then all released at once at the time , will be referred to as the simultaneous-release policy. We use this policy as a baseline or benchmark for evaluating the efficacy of other release policies. A staggered-release policy will be represented by relaxed restrictions for G1 at time and for G2 at with levels . Examples of the schedules of the benchmark simultaneous-release policy and a staggered-release policy are depicted in Fig. 3(a) and (b), respectively. Staggered-release scenarios involving releasing Group 2 at will also be considered.

Fig. 3

Depiction of examples showing the policy switching for the three groups. The top plot represents the benchmark simultaneous-release scenario in which all groups have restriction level for and all groups have relaxed restriction level for . The bottom diagram corresponds to one staggered-release scenario in which the times for policy switching are at and : all groups have the restriction level for and the level for , but Group 1 has a relaxed level of restriction for .

We will explore the effects of different policies by simulating Model (1) with various schedules of relaxing restrictions prescribed by chosen paths determined by , , and .

Analysis and results

In this section, whenever we state results about staggered-release policies, they will be in terms of their performance against the benchmark delayed simultaneous-release policy. One such scenario is depicted in Fig. 4. It illustrates that when all groups are released after a long period of strong restrictions (larger ), there could be a second wave with a high peak if the residual restriction level is not sufficiently high and there are no other intervention measures in place to prevent it. This figure is produced by a single-group model using parameter values within the ranges listed in Table A.1.

Fig. 4

Illustration of a potential scenario where a second wave can occur with high peak after lockdown and social distancing restrictions are relaxed. is the time when the restriction started and is the time when the restrictions are lifted.

Table A.1

Definition of the symbols used in Model (1) and their values used in the examples.

Symbol	Description	Values & range
k_i	Rate of progression from E_i to I_i, i.e., 1∕k_i is the mean latent period	(1/10, 1/7)
γ_ai	Rate of recovery, i.e., 1∕γ_ai is the mean infectious period for A class	(1/10, 1/5)
γ_i	Rate of recovery, i.e., 1∕γ_i is the mean infectious period for I class	(1/10, 1/5)
R_0i	Basic reproduction number for group i	(1.7, 3.6)
a_0i	Effective contact rate in the absence of intervention	See Eq. (4)
a_i(t)	Effective contact rate under intervention such as social distancing	See Eq. (5)
θ_i	Infectivity ratio of A_i to I_i individuals, 0≤θ_i≤1	(0.3, 0.6)
χ_i	Infectivity ratio of H_i to I_i individuals, 0≤χ_i≤1	(0, 0.2)
ϵ_i	Level of preference for contacting one’s own group, 0≤ϵ_i≤1	(0.5, 0.9)
c_ij	Proportion of contacts a member of group i has with group j	See Eq. (3)
λ_i(t)	Force of infection for susceptibles in group i at time t	See Eq. (2)
p_i	Proportions of infectious that are symptomatic. p_1≤p_2≤p_3	(0.4,0.8)
q_i	Proportions of disease deaths from the H class	(0.00064, 0.032)
ϕ_i	1∕ϕ_i is the mean duration in H_i before disease death	(0.05, 0.1)
δ_i	Rate of disease death from the I_i class	(0,0.0008)
η_i	Rate of moving from the I_i class into the H_i class	(0.2, 0.5)
T_j	Time points when a new policy starts. T_0≤T_1≤T_2≤T_3=T_end	Vary
sj(i)	Reduction of contacts for group i in time interval (T_j−1,T_j]	Vary
s_b=s1(i)	Reduction of contacts during the initial restrictions, i.e., in (T_0,T_1]	0.8
s_r	Reduction of contacts in the residual period t>T_end	0.2

Note: . Time unit is days. Most parameter ranges are based on [1], [16], [19], [28].

Fig. 4 is a special case of Model (1) when all groups are identical so there is only one group. We selected a value for to be between the lower estimates of 2–3 (assumed in many modeling studies for COVID-19) and the high estimates such as those reported in [19], which is 5.7 (95% CI 3.8, 8.9). The value provided a better match of the data of new cases shown in Fig. A.1 (see Fig. A.1 in the Appendix). Other parameter values used in this figure are , , , , , , , , .

Fig. A.1

(a) Data of new cases for NYC. (b) New cases from our model.

When three groups are considered, because of the higher activity level in the younger group, would be much higher than other groups. For given (), the overall basic reproduction number for the entire population can be calculated using the formula derived in the Appendix (see Appendix A.1). In most of the numerical simulations in this section, we used the values , , and , for which . To demonstrate that the results are not sensitive to the choice of , examples with other lower or higher values are also included in the Appendix. Some other parameters are also assumed to have different values for the three groups. For example, (i.e., Group 1 has a much lower probability of dying from COVID-19 than the older groups, as shown in Fig. 1).

We will compare the total number of disease deaths in Group from to the end of the outbreak under a staggered-release policy and the benchmark policy, denoted by and , respectively, which can be calculated by (or equivalently ) for a given policy. A staggered-release policy is evaluated based on the reductions in the cumulative deaths of group in comparison with the benchmark policy, which is termed “efficacy” and denoted by , i.e., Elderly and overall deaths are the primary outcomes. We will generally also present the effect of a staggered-release policy on the projected overall number of deaths in all groups, since this provides context for the possible trade-offs that occur in order to save elderly lives. Let denote the overall efficacy. Another consideration in the comparison of policies is the peak size of the total symptomatic infections.

Releasing the low-risk young group first saves lives

In this section, we analyze the costs and benefits of releasing Group 1 from restrictions before other groups. In this case, . Examples of releasing all three groups at different time points are considered in Section 3.2. To identify the effects, we hold constant the total duration of restrictions and the severity of restrictions on Group 2 and elderly Group 3. At , there is a choice of whether (and to what degree) to release Group 1 from restrictions; a choice of .

For ease of presentation, in the remainder of this paper we will shift the time axis to have (the time when the initial restriction starts) and let and represent the times from . The values at of the state variables of the original system (1) without intervention will be taken as the initial values at for the shifted system. The comparisons of outcomes between staggered- and simultaneous-release scenarios will be considered based on changes after .

Fig. 5 presents how the efficacy for the elderly group () and the overall efficacy () vary with the policy parameters and . Several staggered-release policies are represented by and , while the delayed simultaneous-release policy corresponds to or . Other parameter values are: , , , , , , , , , , , . The efficacies for Group 3 and overall are illustrated in Fig. 5(a) and (b), respectively.

Fig. 5

Dependence of the efficacy for elderly () and the overall efficacy () on the timing and restriction level .

Results shown in Fig. 5 suggest that there is a significant reduction in elderly deaths and overall deaths as the low-risk young group faces less severe restrictions () and as this relaxation of restrictions takes place at an earlier date (), and these responses are strictly monotonic in both and . Particularly, the timely release of Group 1 with and can decrease elderly and overall disease deaths by 7% to 28% of their projected values under the simultaneous-release policy, for the parameter values used. Results for different sets of parameter values are discussed in later sections.

Fig. 6 illustrates the influence of the policy parameters and on the total symptomatic infections (per 100,000). All parameter values are the same as in Fig. 5. It shows that in all 4 plots the peak size decreases when decreased from 0.8 to 0.5, but in (a) the peak increased when is decreased further to 0.4. For , the early releases of Group 1 reduce the peak size by 19%–44% in all four plots.

Fig. 6

Dependence of the total symptomatic infections on the policy parameters and . In all four plots, the benchmark simultaneous-release scenario corresponds to the curve with the highest peak, which corresponds to the strict restriction level .

For the scenarios presented in Fig. 5, Fig. 6, we observe that the restriction level would provide the most beneficial policy for all values of , and that and 30 are better choices than or 50. For ease of reference, we refer to these two policies as:

Policy I: and ;

Policy II: and .

If we compare Policy I and Policy II, both have advantages and disadvantages. Policy I has a relatively lower peak size (see the thicker-red curve in Fig. 6(a) and (b)) while Policy II leads to higher efficacies (see the red column in Fig. 5(a) and (b)). Examples of policies with will be discussed later.

Additional simulations of Model (1) with a wider range of parameter values show similar qualitative behavior, some of which are included in the following sections and in Appendix. Particularly, the main results are not sensitive to the choices of group-specific and overall basic reproduction numbers ( and ), proportions of symptomatic infections (), mean infectious period (), and . The following results summarize how the reductions in disease deaths () for the elderly group and overall respond to the early release of Group 1 from restrictions.

Result 1. The following results are based on a large number of simulations of Model (1) with a broader range of reasonable parameter values for COVID-19.

Both the overall efficacy  and the efficacy for the most vulnerable and elderly group  are strictly increasing as the severity of secondary restriction level decreases.

The peak size of the total symptomatic infections may be the lowest at an intermediate value of severity .

For smaller , the efficacy for Group 1 may become negative if the disease mortality for Group 1 is high.

In all simulated results, staggered-release policies are more beneficial than the simultaneous-release benchmark policy in the sense that they help reduce disease deaths in the most vulnerable population and lower the peak size, while allowing increased economic activities.

Results presented in Fig. 5, Fig. 6 suggest that, in the evaluation of staggered-release policies, balanced considerations are needed between reductions of disease deaths and the peak size of symptomatic infections. We may also consider the efficacies for Groups 1 and 2. This is presented in Fig. 7. We observe that the behavior of as shown in Fig. 7(b) is very similar to that of Group 3 as shown in Fig. 5(a), while Fig. 7(a) shows that the changes of disease deaths for Group 1 are much smaller. This is because the disease mortality for Group 1 is much lower. Parameter values are the same as in Fig. 5.

Fig. 7

Similar to Fig. 5 but for Groups 1 and 2. All parameter values are the same as in Fig. 5.

Fully staggered releases require proper timing

Results in the previous sections present the positive impacts of releasing Group 1 before the other two groups are released at the same time. In this section, we explore the impact of relax restrictions for Group 2 at some point before the final releases of all groups at . The intuition would be that since Group 2 is less vulnerable than Group 3, staggering their releases may have similar benefits to the early release of Group 1 only. We will observe a similar result here, but with some caveats.

Let represent the timing for releasing Group 2 with a restriction level . As before, the restriction level for Group 3 will be kept at 0.8 for all , but now is another policy variable. Assume that the severity level for Group 1 does not change in , i.e., . As more detailed analysis of the degree of Group 1 release, , and the timing have already been done in the previous sections, in this section we will fix and . We will vary the cutoff time and the severity . Note that the staggered policy considered in the previous section corresponds to . We refer to this as the “G1-release” policy. The baseline simultaneous-release policy will be the same idea as before, with all three groups being held under the same initial restrictions until . This is represented by with .

Two natural comparisons are the “G1-release” policy and the “G2-release” policy with . One example is illustrated in Fig. 8, in which , , , with various values of . All other parameter values are the same as in Fig. 5, Fig. 6.

Fig. 8

Comparison of the G1-release and G2-release policies represented by the solid and dashed curves, respectively, in (a) and (b). The efficacies for Group 3 and overall are shown in (c) for four values of , where and 0.8 correspond to the G1-release and G2-release, respectively.

We observe from Fig. 8(a) that the G2-release () has a higher peak size than the G1-release (), although both peaks are much lower than that of the simultaneous-release. Fig. 8(b) shows that the increased infections of the G2-release resulted in an increased number of disease deaths in Group 2, particularly during an earlier period of the outbreak. Plot (c) shows that the early release of Group 2 generated slightly higher efficacy in the elderly group 3 but decreased the overall efficacy, which is due to the increased deaths in Group 2. The changes in peak size, number of deaths, and efficacy vary with the timing . This suggests that a careful selection of and a more detailed cost-benefit analysis should be carried out to determine the most beneficial policies for releasing Group 2.

Robustness of the simulation results

In most of our simulations of Model (1) using a broad range of parameter values, the qualitative behavior concerning the efficacies and the peak sizes of symptomatic infections are similar. Several examples are provided in this section and in the Appendix.

Fig. 9 illustrates results for reduced group reproduction numbers: , and . All other parameters are the same as in Fig. 5, Fig. 6 with and . As decreases from 0.8 to 0.4, plot (a) shows that the peak size is the lowest at ; plot (b) shows that the cumulative numbers of deaths decrease monotonically in Groups 2 and 3 but remains almost the same in Group 1; and plot (c) presents the corresponding efficacies for all groups. These plots show a similar qualitative behavior as in Fig. 5, Fig. 6 (see the case of and ).

Fig. 9

Influence of the policy parameter on (a) symptomatic infections, (b) cumulative disease deaths, and (c) Efficacies .

To explore the influence of , Fig. 10 illustrates the results for . All other parameters are the same as in Fig. 5, Fig. 6 with and . It shows that, as decreases, the peak size increases while the efficacy decreases within each group. One of the reasons for the larger peak size for smaller values (e.g., 80–100) is because the number of infections are still in the increasing phase at the time of releasing all groups, whereas for the infection curve has already started decreasing. This suggests that the timing of group releases can be critically important.

Fig. 10

Effects of on the efficacy with fixed days and . All other parameters are the same as in Fig. 5, Fig. 6.

The significantly lower case mortality for the younger groups shown in Fig. 1 was considered in the choice for values of for the figures in previous sections. If, however, in a different population is not so much lower than and , a more detailed cost and benefit analysis might be needed in the examination of policies for releasing Group 1 earlier. In Fig. 11, the values , , are used. All other parameters are the same as in Fig. 5, Fig. 6 with , and an additional of is included in the comparison. One key difference between the results shown in Fig. 11 and others presented earlier is that, the least strict restriction () may shift the peak earlier in time with a much higher peak size than the case of (see (a)). On one hand, this can greatly increase the population immunity, which may reduce the loss of life in the elderly group 3 (see (b)) and increase the efficacies for both Group 3 and overall (see (c)). On the other hand, however, the efficacy for Group 1 is negative for (see (c)), which represents increased deaths in Group 1, although by a very small number. This can be avoided by restricting to be higher than 0.2 to ensure non-negative , as done in the earlier examples.

Fig. 11

Similar to Fig. 9 but for , . The values of are the same as in Fig. 5 and an additional value of is included.

Discussion

The findings generated by a careful study of various scenarios support the view that releasing some appropriately defined groups earlier from lockdown and social distancing restrictions may be highly beneficial under the premise that reducing overall disease deaths is the priority. The use of a staggered-release policy built on the immunological strength of sub-groups of a population may indeed prove to be highly effective in terms of protecting vulnerable individuals while mitigating economic setbacks.

Staggered-release policies are designed to avoid dangerous levels of severe infections among the most vulnerable by redistributing the timing and likelihood of infection across groups. The release of low-risk and young from restrictions earlier might increase their disease exposure, but most infections would be asymptomatic or mild. The second group would be then released after some level of population immunity has been achieved among the first group. The related concept of deploying recovered individuals with antibodies into society to restart the economy while keeping infectious contacts low has been put forward in Weitz et al. [27], who term it ‘shield immunity.’

The effectiveness of staggered-release policies is crucially tied to timing choices. We show that the mistimed release of groups may lead to an increase in overall disease deaths. While the early release of the youngest group is generally going to be significantly helpful, the response of the epidemic curve to specific choices of when to release older groups is complicated, as pointed out more generally by Morris et al. [17].

Simulations show that if the second group’s release is timed to take place just before the peak of infections, then the release may lead to a higher peak of infections and more disease deaths. On the other hand, releasing this second group after the peak of infections generated by the first group has passed produces only a small or no bump in infections that does not increase significantly the risks of those involved, a result consistent with the findings in Morris et al. [17]. Morris et al. observe that optimal control in this context responds to changing levels of susceptibility and infection rates in the population.

Compared to the delayed simultaneous-release benchmark, staggered-release policies that release Group 1 first can save a significant portion of lives among the elderly. In most of the examples presented in this paper, we see commonly 10%–20% reduction in elderly disease deaths compared to the simultaneous-release policy. The estimated size of this effect is robust to alternate parameter specifications.

Results shown in Fig. 10, Fig. A.3 reveal a stark change in outcomes if the release of older groups is mistimed, even if it is carried out a few weeks too early. Releasing these individuals from restrictions while infections are still rising results in a significantly higher peak size of symptomatic infections compared to those from a policy of waiting to release them only after infections have begun to fall.

Fig. A.3

Symptomatic infections under four different policies corresponding to with (a) and (b) and . Both plots are for and . The benchmark simultaneous-release policy () is also plotted. Plot (a) shows a much higher peak size than plot (b) for .

Importantly, while economic considerations were mentioned as a motivating factor behind the earlier removal of restrictions for some groups, there were no actual economic components in the model we used in this paper. Our conclusion that staggered-release policies outperform simultaneous-release policies is shown through a direct comparison of health risks and aggregate loss of life estimated from projections of the outbreak under these policies. The fact that earlier release of certain sub-groups enhances economic activity on the whole is a further benefit to staggered-release policies, but this is a benefit that exists outside of our model.

It would be insightful to incorporate the adverse health consequences associated with lengthy lockdown restrictions, such as increased risk of depression and loss of health insurance from employment to cover unrelated medical issues, within the comparison structure used in this paper. However, we wanted the focus to be on a direct comparison of health risks associated with the disease itself, so that it is clear that these additional consequences do not drive the result. We leave the explicit consideration of these additional health risks to future work.

Ultimately, it was our goal to analyze the use of staggered-release policies and compare the outcomes to those resulting from a policy of the delayed simultaneous release of everyone. This is meant to be an illustrative paper, not a prescriptive one. Implementing staggered-release policies would require a more thorough analysis of the specific population under consideration in order to appropriately structure the schedule of release.

Conclusion

Lockdown and social distancing restrictions have been an important tool in controlling the spread of COVID-19 and keeping infection rates low enough to avoid exceeding health care capacity. Now that the initial phase of the outbreak has been addressed, the next step is to properly schedule the relaxation of constraints in order to optimally control the entirety of the outbreak. The decision to lift restrictions on everyone simultaneously is perhaps the natural policy, but this can have dire consequences if a second peak arises as a result of insufficient immunity in the population. Such an outcome would defeat the purpose of installing restrictions in the first place, so it is crucially important to consider more effective policies in terms of when to release people from restrictions.

To summarize, the early release of low-risk young individuals provides an interim period that allows them to build up some degree of immunity without facing significant health risks themselves, in order to protect older and more vulnerable individuals from severe health risks by reducing their likelihood of infection and death. As an added benefit, the staggered-release policy allows for more activity in the population at an earlier date than the delayed simultaneous-release policy would. Thus staggered-release policies can positively impact everyone involved by allowing people to return to their lives without risking them.

Declaration of Competing Interest

The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.