Abstract
After one year of stop-and-go COVID-19 mitigation, some European countries still experience sustained viral circulation due to the B.1.1.7 variant. As the prospect of phasing out this stage through vaccination draws closer, it is critical to balance the efficacy of long-lasting interventions and their impact on the quality of life. Focusing on the current situation in France, we show that moderate interventions require a much longer time to achieve the same result as high intensity lockdowns, with the additional risk of deteriorating control as adherence wanes. Integrating intensity and duration of social distancing in a data-driven “distress” index, we show that shorter strict lockdowns are largely more performant than longer moderate lockdowns, for similar intermediate distress and infringement on individual freedom. Our study shows that favoring milder interventions over more stringent short approaches on the basis of perceived acceptability could be detrimental in the long term, especially with waning adherence.
The emergence of SARS-CoV-2 B.1.1.7 variant in December 20201,2 disrupted the management of COVID-19 pandemic in Europe. The alert arrived as some governments were lifting interventions that had been applied to curb the second wave. Some countries, such as the UK and Ireland, were forced to rapidly implement strict lockdowns to control the explosion of cases due to the variant. Others maintained or strengthened their restrictions because of concerns over the new variant3.
Few months after, with vaccination lagging behind (25% of the population of the European Union with a first dose on May 1 vs. 44% in the US, 51% in the UK, and 62% in Israel4) and a third wave due to the B.1.1.7 variant, continental Europe faced the challenge of relying once again on heavy restrictions to reduce sustained viral circulation and improve the epidemic situation approaching the summer. But what is the optimal strategy, given vaccination rollouts, current epidemic conditions, and the sustainability of perduring restrictive policies? On one side, there remain limited available options beyond high intensity interventions, once milder layers of social distancing have been accumulated, strengthened, and extended over time (e.g. curfew, closure of restaurants and bars, while maintaining schools open). On the other side, the efficacy and long-term sustainability of the adopted policies are potentially threatened by loss of adherence and policy-induced fatigue5,6, affecting the quality of life of the population.
Building on observed adherence waning and introducing a data-driven measure capturing the limitations on individual freedom resulting from restrictions, we compared intervention scenarios of varying intensity and duration, and examined the role of adherence and sustainability on optimal epidemic control. The study is applied to the third wave in Île-de-France – the Paris region, the most populated of France and heavily hit by the pandemic – accounting for vaccination rollout plans, seasonality, and plans for the phasing out of restrictions.
RESULTS
Adherence and impact of interventions of varying stringency and duration
We used an age-stratified and regionally structured transmission model, fitted to hospital admission data and validated against serological studies’ estimates7. The model integrates mobility data to parameterize social contacts over time accounting for social distancing measures. It describes the co-circulation of the historical strains and B.1.1.7 variant, matching virological and genomic surveillance data8 and considering the variant’s estimated advantage in transmissibility1,2,9 and hospitalization rate10,11. B.1.1.7 variant was estimated to become dominant in the region by mid-February 20218.
We fitted the model to daily hospital admission data till week 11, when the region was still under curfew before strengthened measures were applied on March 20 (end of week 11) to control the rise of the third wave. We compared projected trends based on unchanged curfew conditions with the trajectories resulting from the application of more stringent interventions in week 12. Stringency and efficacy of these interventions were based on past mobility data12 and modeling estimates7,13 (Methods). Figure 1 shows overall mobility variations compared to pre-pandemic levels extracted from mobile phone data12 during the three pandemic waves (Methods). During the first wave (March-May 2020), a full lockdown was enforced with strict stay-at-home orders, schools were closed, and only essential services were allowed to operate. During the 8 weeks of lockdown, mobility showed a reduction of about 70% in the region (65% at the national level12) that remained fairly constant over time. During the second wave (October-December 2020), a milder lockdown was implemented for 6 weeks. Schools remained open and a larger number of job sectors were allowed to operate. Recorded mobility and estimated presence at work decreased but remained almost two times higher compared to the first lockdown (average mobility reduction of about 40% in the first three weeks compared to pre-pandemic levels) and showed a rapid and marked increase over time. This loss of adherence occurred remarkably faster (after the third week) and more substantially during the second lockdown compared to the first, and it was estimated to increase the effective reproductive number by about 11% (Methods). Higher mobility was registered in the last two weeks due to the reopening of all retail for Christmas shopping. Meanwhile, fear of contracting COVID-19 decreased over time since the second wave (average −0.31% per week), whereas anxiety increased (average +0.13% per week; Figure 1) suggesting a progressive buildup of distress with time under any social distancing measure.
Top row: weekly hospital admissions in Île-de-France during the first (left), second (center), and third (right) pandemic wave. Dots correspond to data, curves and shaded areas correspond to the fitted trajectories and 95% probability range. Horizontal dashed lines refer to the peak of the first and second wave in the region. Central row: mobility reduction in Île-de-France during the first (left), second (center), and third (right) pandemic wave. Yellow histograms represent the variation of mobility with respect to pre-pandemic levels, based on the number of trips extracted from mobile phone data12. Blue curves show the estimated change in presence at workplace locations over time with respect to pre-pandemic levels based on Google location history data14. Shaded rectangles in the plots of the first two rows correspond to social distancing measures (strict lockdown in the first wave, moderate lockdown in the second wave, strengthened measures in the third wave). The second week of the second lockdown and the third week of the strengthened measures against the third wave have lower mobility and presence at workplaces due to bank holidays in the week. Vertical dotted grey lines correspond to school holiday periods. Bottom row: percentage of individuals avoiding crowded public places15 (left), percentage of individuals scared to contract COVID-1915 and prevalence of anxiety in the context of COVID-19 epidemic (center)16 as functions of time; scattered plot between the percentage of individuals scared to contract COVID-19 and the percentage of individuals avoiding crowded places (right) in the time period October 2020 - April 2021, with the results of a Pearson correlation test. Results for these indicators refer to the national scale.
We modeled the effect of a strict lockdown and a moderate lockdown starting week 12 for a duration ranging between 2 and 8 weeks. These scenarios were based on measured mobility reductions and estimated transmissibility conditions during the first and second lockdowns, respectively, before relaxation emerged. For lockdowns longer than 2 weeks, we compared scenarios assuming full adherence with situations characterized by a loss of adherence over time, equal in magnitude as the one observed during the second lockdown, but occurring earlier (starting two weeks after the implementation of restrictions). We considered two scenarios: a limited loss of adherence with a one-time drop, and a continuous loss with repeated drops by the same amount every two weeks. More details are reported in the Methods section.
Vaccination prioritized to older age classes and with efficacies based on available estimates17–19 was simulated according to observed and planned daily rhythms (from 100,000 doses/day at the national level in February to 200,000 first doses/day in March, further increased to 300,000 first doses/day from April, see Methods), following government announcements.
Regardless of adherence, the strict lockdown was predicted to be the only measure able to achieve a rapid decrease of the epidemic trajectories (Figure 2, Figure 3), in line with observations in the UK and Ireland following similar interventions. It would outperform moderate lockdowns of any duration, on both short- and longer-term epidemic impacts. Starting from about 3,000 weekly hospitalizations at the time measures were applied, admissions would be reduced to less than 400 when exiting a strict lockdown of at least 1 month vs. more than 2,000 after a moderate lockdown (Figure 3a). Even with adherence waning, a strict lockdown would reduce the epidemic to the levels recorded at the exit of the first lockdown in May 2020 (670 weekly admissions in w20, 433 in w21, 330 in w22 in 2020), and it would be maintained low by increasing immunization rates. These levels would also enable a better control of viral circulation through test-trace-isolate when partially alleviating restrictions7,20. Importantly, a short circuit-breaker21 of 2 weeks, after which curfew was restored, was predicted to be already enough to rapidly reduce hospitalizations to levels below the ones of February 2020. Rebounds at the end of the short lockdown would be prevented by maintaining a certain degree of social distancing (curfew) and increasing immunization, with stronger reductions over time for increasing vaccination rhythms (from 300,000 to 500,000 first doses/day since April; Figure 2).
From left to right: increasing lockdown duration, expressed in weeks. Top row: vaccination pace accelerated to 300,000 first doses/day since the start of April; bottom row: 500,000 first doses/day. Interventions are applied in w12 and assume a delay of one week to the peak in hospital admissions. Dots refer to data; filled dots correspond to the data used to fit the model and to provide the trajectory for the curfew scenario; void dots correspond to more recent data. Curves refer to the median trajectory; shaded areas around the curves correspond to the 95% probability range obtained from 250 stochastic simulations. The type of intervention is coded by different line colors. Horizontal dashed lines refer to the peak of the first and second wave in the region. Results for strict lockdown scenarios with full adherence or loss of adherence overlap; for this reason, we do not show the scenario with limited loss of adherence.
Weekly hospital admissions at the end of the lockdown (a), weekly hospital admissions at the end of June (w26) (b), cumulative hospital admissions computed over the time period w12-w26 (c), hospital pressure, defined as the number of weeks in which hospital admissions remain above the peak level achieved during the second wave, in the period w12-w26 (d) as functions of the adherence level - full adherence over time, limited loss of adherence, continuous loss of adherence over time. The point with the curfew (gray circle) represents the estimate under the curfew scenario with no additional intervention, and is shown for comparison. Results refer to a vaccination rhythm accelerated to 300,000 first doses/day since April. Symbol types refer to the stringency of intervention (squares representing a strict lockdown scenario, diamonds representing a moderate lockdown scenario). Color shades of the symbol contour refer to the duration (weeks) of the lockdown intervention (from the lightest shade corresponding to 2 weeks, to the darkest one corresponding to 8 weeks). Adherence is coded with the fill color (filled symbols with the color of the scenario correspond to scenarios with full adherence, void symbols represent scenarios with limited loss of adherence, blue filled-in symbols correspond to scenarios with continuous loss of adherence). Horizontal dashed lines in panel (a) refer to the peak of the first and second wave in the region; horizontal dashed line in panel (b) refers to the level of mid-May 2020 at the exit of the first lockdown.
Obtaining results equivalent to a short strict lockdown would require moderate interventions to last longer than 2 months, and could potentially be compromised by loss of adherence to restrictions (Figure 3a). This could slow down and stop the decrease in hospital admissions, leading to a plateau or a rise in hospitalizations after several weeks of moderate lockdown, potentially higher than the peak of the third wave (Figure 2, Figure 3a). This occurs in our scenarios as repeated drops in adherence over time may reduce the efficacy of a lockdown to values lower than a simple curfew after a few weeks (Methods). Since moderate lockdowns would not be able to considerably reduce viral circulation, they would entail a larger impact on the hospital system (median hospitalizations in the period w12-w26 around 38,000-50,000 compared to 10,000-23,000 for strict lockdowns, Figure 3c) for a longer time (median 6-10 weeks with hospitalization incidence above the peak of the second wave compared to at most 2 weeks for a strict lockdown of any duration, Figure 3d). This impact would be more substantial if adherence waned, leaving the hospital system under high pressure for twice the amount of time (median 12 weeks above the peak of the second wave assuming continuous adherence loss, compared to 6 weeks for full adherence, corresponding to 80% of the time period under study).
Despite different trajectories, our model anticipates that moderate lockdowns could reach at the end of June the hospitalization levels measured in May 2020 (670 weekly admissions) if adherence was maintained over time, similarly to a short strict lockdown, and with no advantage compared to curfew measures at this stage (Figure 3b). Adherence loss would lead to higher hospital admission levels.
Assuming the slowdown of the epidemic estimated during winter holidays in February8, we found that the two weeks of school closure for spring holidays (w15-16) would have a significant impact on the efficacy of moderate lockdowns, otherwise with schools open. They would allow flattening the curve and avoiding even longer plateaus in critical conditions before the accrued effect of immunization would decrease the epidemic.
Optimizing interventions’ sustainability by minimizing policy-induced distress
Another critical dimension associated with the nature of interventions – besides their stringency, duration, and adherence – is their sustainability over time, which is a combination of intensity of restrictions and how long they last. To account for this aspect, we introduced a distress index, integrating the intensity, duration, and adherence level in each scenario, and providing a quantitative measure of the policy-induced distress perceived on average by an individual. The higher the distress index and the lower is the sustainability of the measure. The index is computed by multiplying the recorded population-level mobility reduction of each intervention and its duration (Methods). For a given duration, this index provides a lower level of distress associated to smaller reductions of mobility, corresponding to a higher individual freedom given by milder restrictions and/or lower adherence. The index is defined on a range from 0 to 10, with 10 representing the highest distress level, associated to 2 months of strict lockdown with full adherence, i.e. the level experienced in the first wave in March-May 2020, and 0 representing absence of restrictions.
Moderate lockdowns of less than 6 weeks are all characterized by low levels of distress (<4), similar to those of a curfew and of a 2-week strict lockdown (Figure 4). In this range of distress values, a net advantage was observed for the short strict lockdown that substantially reduced the total number of hospitalizations (23,000 vs. an average of 47,000) and hospital pressure (2 weeks vs. more than 8 weeks). High values of the distress index (>7) were associated exclusively to strict and long lockdowns (of 6 weeks, with full adherence, or longer, also with adherence waning over time), which correspond to the most performant measures in suppressing viral circulation and reducing the healthcare impact, but also the least sustainable.
Weekly hospital admissions at the end of June (w26) (a), cumulative hospital admissions (computed in the time period w12-w16) (b), hospital pressure, defined as the number of weeks in which hospital admissions remain above the peak level achieved during the second wave, in the period w12-w26 (c) as functions of the distress index. Results refer to the accelerated vaccination pace of 300,000 first doses/day since April. Symbol types refer to the stringency of intervention (squares representing a strict lockdown scenario, diamonds representing a moderate lockdown scenario, void circle representing the projection under the curfew scenario). Color shades of the symbol contour refer to the duration (weeks) of the lockdown intervention (from the lightest shade corresponding to 2 weeks, to the darkest one corresponding to 8 weeks). Adherence to moderate and strict lockdowns is coded with the fill color (filled symbols with the color of the scenario correspond to scenarios with full adherence, void symbols represent scenarios with limited loss of adherence, blue filled-in symbols correspond to scenarios with continuous loss of adherence). Horizontal dashed line in panel (a) refers to the level of mid-May 2020 at the exit of the first lockdown.
There exists, however, quite a diversified range of intervention options that, for moderate distress (index between 4 and 7), would achieve better control of the epidemic than curfews and moderate lockdowns. One-month strict lockdowns would largely outperform moderate interventions in terms of health metrics while inducing similar distress, as the latter must be maintained for a longer duration. Six- or eight-week moderate lockdowns would lead to about three times as many patients hospitalized and about three to six times the hospitalization incidence at the end of June compared to interventions that exhibit a similar distress level as a strict lockdown of 1 month. If moderate lockdown restrictions were less respected over time, epidemiological and healthcare indicators would considerably worsen, for relatively small gains in lowering the policy-induced distress.
Vaccination and seasonality while managing reopening plans
Further acceleration of vaccination pace coupled with the potential effect of seasonality22 may act in synergy to (i) counteract the deterioration of the epidemic due to the waning of adherence over time, or (ii) - if stronger-bring down the epidemic faster than what is expected from moderate interventions. Figure 5 shows the interplay of these factors assuming that seasonality acts on reducing viral transmission starting May. Keeping the planned vaccination rhythm at 300,000 first injections per day since April, a 5 to 10% reduction of transmission induced by seasonality would be necessary to absorb the potential loss of adherence against moderate interventions by the end of June (label (1)). Without counting on seasonal effects, vaccination rollout should increase by 33% (i.e. from 300,000 to 400,000 first doses per day). Larger seasonality (>20%) or accelerations in vaccination rollouts (up to 800,000 first doses per day) would be able to compensate for the larger cumulated number of patients requiring hospitalization due to adherence loss (label (2)). Reaching by the start of the summer the weekly admissions achieved by an imperfectly adhered 1-month strict lockdown would require substantial seasonality coupled with large increases in vaccination rhythms (contour line at 300 in the top row of Figure 5).
Heatmaps show median values of weekly hospital admissions at the end of June (w26) (top row) and cumulative hospital admissions in the time period w12-w26 (bottom row), as functions of vaccination rhythm (y axis) and seasonal reduction in transmission (x axis) for moderate lockdowns of 8 weeks with full adherence (left), limited loss of adherence (center), continuous loss of adherence (right). Contour lines indicate reference values of specific scenarios defined in the legends. Circled numbers refer to a subset of conditions of interventions (heatmaps from left to right), vaccination and seasonality (variables in the y and x axis in each heatmap) achieving the same outcome, identified by the contour lines (see legend) and discussed in the main text. Plots in the third row show projections of the weekly hospital admissions under different hypotheses for the reopening conditions, assumed right after lifting the moderate lockdown (lighter line color), or through a progressive transition (darker line color): conditions experienced in mid-July 2020, but with schools in session (left); curfew scenarios with 40% fewer individuals respecting physical distancing (center); curfew scenarios with 15% fewer individuals respecting physical distancing (right). Scenarios assume a 10% reduction in transmissibility due to seasonality (except for the mid-July 2020 conditions that already embed seasonal aspects) and a vaccination rhythm of 300,000 first doses per day starting April.
In all situations, a certain degree of social distancing is required to accompany the gradual lifting of lockdown to avoid slowdowns or rebounds (bottom row of Figure 5). Even the summer conditions estimated in mid-July 2020, but considering schools in session, may lead to an epidemic resurgence if incidence is high, despite the growth in population immunity and summer seasonality. Results show that a progressive transition in phasing out restrictions is essential, and they further support the importance to lower the incidence level to better manage potential rebounds while reopening.
DISCUSSION
Managing sustained viral circulation after long periods of social distancing measures of varying intensity faces the challenge to reduce the strain on the healthcare system and to limit perduring or stringent interventions affecting the quality of life of the population. Moreover, with accelerating vaccination campaigns and the prospects of reopening the society, adherence waning may represent a threat to phasing out restrictions. Using Île-de-France as a case study, we compared the efficacy of different measures against their sustainability and potential for case resurgence due to imperfect adherence of the population. Given the high incidence levels reached by the epidemic in the region by mid-March, exceeding the peak of the second wave8, only high intensity interventions would have been able to rapidly curb viral circulation, allowing the region to reduce significantly the burden of hospitalization after only 2 weeks and despite loss of adherence. Once incidence substantially declined, the management of the epidemic could largely benefit from test-trace-isolate strategies7,20 and the large-scale availability of self-test kits for iterative screening23, while immunization due to vaccination builds up in the population. Hospitals could more rapidly restore routine care beyond COVID-19. Moreover, rapidly reaching low incidence levels would also lower the potential for SARS-CoV-2 evolution conferring fitness advantages, and allow a better control of the possible emergence or importation of variants of concern24.
Moderate interventions as the strategy adopted in November-December 2020 to curb the second wave constitute suboptimal options for the management of the epidemic through summer. Their efficacy remains limited because of B.1.1.7 higher transmissibility1,2,9 and still low immunization levels (15.7% population vaccinated with a first dose in the region by April 20). Our results show that these measures should be maintained for much longer (>2 months) to reach incidence values similar to the result of a short and strict circuit-breaker21, at the expense of a large number of severe cases requiring hospital care, a continuously high pressure on the hospital system, and high levels of distress cumulated over time. The strengthened measures currently in place, based on closure of non-essential businesses, ban on gatherings and recommendations to telework, are similar in intensity to the moderate scenarios considered in this study, as also confirmed by the similarity of the stringency index25. These were however accompanied by the advanced closure of schools just before the 2-weeks regular school holidays in April that provided an extra break on the epidemic evolution showing in the current trend.
Despite different trajectories, epidemic conditions by the time summer starts are predicted to be similar across intervention scenarios (with the exception of high intensity interventions lasting one month or longer, largely suppressing the epidemic) and with the curfew scenario in absence of additional restrictions. Differences in how the epidemic is managed throughout spring are absorbed over time thanks to vaccination. However, large disparities remain for cumulative epidemiological and public health indicators, depending on whether early suppression or mitigation were achieved by the interventions. This would have an impact not only throughout the third wave (by increasing the overall number of hospitalizations, patients requiring critical care, and deaths), but also on the medium-to-long term due to the rising number of individuals who are likely to suffer from long-term health consequences following a COVID-19 infection (long COVID)26,27. Early estimates indicate that about 10% of individuals testing positive for COVID-19 exhibit symptoms after 4 months26, and about 3/4 of hospitalized patients report at least one symptom after 6 months27 – mainly fatigue, muscle weakness, sleep difficulties, anxiety, depression. Choosing a 2-week strict lockdown against an 8-week moderate lockdown would correspond to estimated 30,000 avoided long COVID cases among detected infections from mid-March to the end of June.
The choice of interventions inevitably also impacts the quality of life of the population due to imposed restrictions, leading to possible spontaneous relaxation. The shorter the measure’s application, the less likely it is to observe adherence waning over time. Interventions of high intensity but short duration may therefore constitute an optimal approach to reduce both epidemic and healthcare burdens, while minimizing possible loss of adherence as well as policy-induced distress. Evidence from OECD countries after one year of COVID-19 pandemic show that swift lockdown measures were overall less restrictive of civil liberties, thanks to achieved control, compared to recurrent mitigation policies severely impacting individual freedom28. Indeed, moderate or mild (curfew) interventions cumulate distress over time, as they need to be implemented for much longer to achieve the reduction of health indicators, with the potential risk of losing population adherence. This would considerably worsen both incidence and cumulative indicators, slowing down or stopping the decrease in incidence obtained with restrictions, thus remaining on a long plateau at sustained viral circulation, as occurred after the second wave. If relaxation against measures is left uncontrolled, epidemic rebounds can also be expected. At the same time, loss of adherence would correspond to a limited gain in personal freedom, when averaged over all individuals (−26% in distress index by continuously losing adherence in moderate lockdowns lasting 8 weeks), compared to interventions of lower stringency (−60% in distress index from an 8-week moderate lockdown to a 2-week strict lockdown).
In our study, loss of adherence occurs over time and is informed from observed increases in mobility during the second lockdown and corresponding estimated impact on the epidemic, during unchanged restrictions and recommendations. We did not consider initial adherence to restrictions different from what was measured in the first and second lockdown. While lower initial adherence may be expected as stringent social distancing measures are being applied for the third time, this may also depend on the acceptability of new measures, clarity of restriction and recommendations. For example, a recent survey showed that about 70% of individuals approved the strengthened measures recently applied in France, however almost half of them planned to disobey the rules29. Also, adherence loss should not be confused with population response to restrictions induced by socio-economic conditions and life circumstances6,12,30,31. Prior work showed that this response – despite numerical differences depending on the stringency of measures (first lockdown, second lockdown, localized curfew at 8pm, nationwide curfew at 6pm) – is associated to the composition of the population, with blue-collar jobs and household crowding emerging as determinants of higher mobility during restrictions in France31.
We introduced an index integrating mobility reduction and duration of restrictions to provide a quantitative measure of policy-induced distress along the spectrum of varying stringency. This is meant to integrate the impact of restrictions infringing on individual freedom, as well as psychosocial effects of prolonged measures, linked for example to isolation, uncertainty, loss of purpose, and lack of social contacts32,33. While both distress index and adherence were informed from data, we did not consider an explicit relation between distress and adherence loss, potentially leading to feedback mechanisms reinforcing relaxation for increasingly long durations. Related to “pandemic fatigue”, a concept often introduced as the presumed cause to limited adherence, this relation remains highly debated6. Some behavioral scientists warned against an ill-defined concept used to justify avoiding strict and/or early interventions6,34. Different features and origins of fatigue are likely at play – including for example life constraints independent of motivation, as discussed above – that would require a range of definitions, data, and frameworks for analysis. A study on data from 14 countries showed that adherence to physical distancing evolved following a U-shape between March and December 20205. However, in France this drop would correspond to the summer period, between the first two waves, during which restrictions were lifted and only recommendations on the use of personal preventive measures were in place. As such, it does not relate to the adherence loss throughout interventions considered here. Different indicators obtained from surveys show that fear of contracting COVID-19 decreased over time after the second wave in France, while anxiety continued to increase in the population. We found a positive association between fear and social distancing (expressed by the percentage of individuals avoiding crowded places), confirming the role that risk perception has in shaping health-related behaviors35. However, we did not find a significant association between increasing anxiety, concurrent with lasting restrictions, and decreasing social distancing. So far there exists little evidence on the mechanisms of action of behavioral interventions that could improve our understanding and be leveraged to boost policy observance.
Available evidence indicates that interventions implemented in 2020 largely reduced the incidence of COVID-1913,36–40, in the absence of effective treatments and prior to vaccination. Substantial differences were observed between analyses aiming to assess the efficacy of single social distancing measures (e.g. closure of schools, businesses, all but essential services, ban on mass gatherings and public events, restrictions on movements and stay-at-home orders). Our study did not focus on isolated measures, but considered the estimated efficiency of policy packages that were deployed during the first and second wave in France, along with observed policy observance and wane in time. A lockdown as strict as the first one is unlikely to reach nowadays the efficiency observed last year, and for this reason we considered reduced adherence, which we show would marginally affect the results. The two lockdowns implemented in 2020 did not differ exclusively for the closure or opening of schools, but also for the mobility levels and presence at workplace estimated from data. Behaviors related to mobility, presence at work and school are not independent and we currently lack enough data to parameterize their relationship. In addition, alternative versions of interventions allowing time outdoor where risk of transmission is reduced41 – such as recommendations currently in place – may reshape mobility, contacts and associated risk in ways different than previously observed, preventing their assessment on the basis of historical data. Open questions remain on the combination and sequence of restrictions to be progressively lifted after the lockdown, as specific measures are too detailed for mathematical models to quantify (e.g. reopening of restaurants). Strategic prioritization will likely depend on countries’ interests.
Vaccination is key to exit the health crisis, however our numerical evidence shows that epidemic management still needs to rely on social distancing to curb viral transmission, confirming prior work42–44. Increasing vaccination rollout coupled with 5-30% reduction in transmission due to seasonal effects would be able to compensate for the slowdown or rebound effects of adherence waning or fast reopening. Multiple studies have investigated the relationship between SARS-CoV-2 transmission and weather. Results suggest that warm and humid conditions, and high UV radiation levels, are less favorable to disease spread22. Based on previous estimates45,46, we derived that the average increase in UV radiation and temperature reported in Île-de-France from March to June corresponds to approximately 10% reduction in transmission (Methods). Additional mitigating effects are expected due to seasonal behavior, with individuals spending more time outdoor than indoor, and aerating indoor settings more compared to winter time. But misconceptions on seasonality may generate excessive trust in the public altering their risk perception, and in authorities affecting their decision-making47. Despite a building literature on the topic, there remain aspects that are difficult to measure and include a strong behavioral component. A large second wave started last year in the United States during summer because of early reopening, and cases started to rise in France from mid-July, paving the way to the second wave in the fall. Lifting restrictions with the conditions experienced in mid-July 2020 is expected to lead to an epidemic rebound if incidence is high. We did not consider here the situation at the end of the first lockdown in spring 2020 because it was characterized by the maintenance of cautious behaviors, and additional levers existed that continued curbing transmission after lockdown was lifted (e.g., the increase in mask use, from 45% in mid-May 2020 to >70% at the end of the summer16, also due to mask mandates). Managing the epidemic while gradually releasing non-pharmaceutical interventions through the summer should mainly rely on the speed of vaccination rollout.
Our study has a set of limitations. It is applied to a region only, as indicators for France hide a variable situation at the local level, limiting the accuracy of modeling approaches extended to the whole country. Geographical heterogeneity depends on the current epidemic situation, population immunity due to natural infection, and variant frequency, so that results are not directly generalizable to other regions. We did not consider waning of immunity48 or reinfections over the time frames modeled. We assumed the transmissibility advantage of the B.1.1.7 variant from early estimates in France9, in agreement with other studies1,2, however this may be altered over time by social distancing and competition with other strains. Assuming a smaller transmissibility advantage for the variant would lead to lower incidence projections, however it would not be able to capture the evolution in time of B.1.1.7 frequency in the region. We did not consider the interaction with other variants, such as 20H/501Y.V2 or 20J/501Y.V3, that are already present in the country and show so far limited diffusion. If these variants can at least partially escape natural or vaccine-induced immunity49, they may pose a challenge for the management of the epidemic as population immunity increases. Our approach is not suited to account for contacts in low-risk and high-risk conditions, e.g. in closed ill aerated settings vs. open settings, but seasonal reductions effectively account for these aspects. Modeled vaccination rhythms according to authorities’ plans were slightly faster than observed. By May 4, 23.6% of the population was vaccinated with a first dose in the model, compared to 20.3% according to data, however this is not expected to affect our findings. We did not consider slowdowns that were recently observed after the temporary stop of AstraZeneca vaccine administration, undermining demand relatively to other vaccines. We considered 50% coverage in the adult population, following the declared intentions to get vaccinated of this age class in France50, but we did not consider that after a certain portion of adults is vaccinated, vaccine hesitancy may rise beyond expectations, also in relation to the successful reduction of epidemic incidence, as currently seen in some countries51. Our findings and prior work show that relaxing social distancing with limited immunization may result in epidemic rebounds42–44. We did not consider the economic impact of social distancing measures, as our study focused on the epidemiological, healthcare, and behavioral components. There is increasing evidence, however, that economic growth, public health, and civil liberties do not need to be in opposition in the management of the COVID-19 crisis, with countries aiming for elimination faring largely better than countries adopting mitigation strategies28. Also, we did not consider health impacts beyond COVID-19 that can result from a high pressure on the hospital system. Psychosocial impact was instead introduced through a simplified empirically-driven indicator based on restricted mobility, the core of many social distancing measures. However, this indicator is an average, therefore it hides the effects on vulnerable populations who may experience disproportionately higher distress6,30,31. Also, being informed by mobility only, it aims at providing a measure of infringement of personal freedoms, but without explicitly capturing other elements associated with the quality of life33.
Control of the epidemic in the next months depends, in non-linear ways, on the interplay between the stringency of restrictions, vaccination roll-out plans, population adherence to measures and vaccination, and plans for gradually lifting social distancing measures towards summer. Mathematical models help to unravel the complexity of these interactions, accounting for the uncertainties characterizing some of these aspects, and to quantitatively inform on the optimal solutions for epidemic control. Our study shows that favoring milder interventions over more stringent approaches limited in time on the basis of perceived acceptability could be detrimental in the long term, especially with waning adherence.
METHODS
Data
Hospital surveillance data
The model is calibrated to regional daily hospital admission data, collected in the SIVIC database52. The database includes the number of admissions of COVID-19 confirmed patients to regular hospital or intensive care units. Hospital data are corrected for notification delays and do not suffer changes in detection or sampling, unlike the number of detected cases. As such, they provide a robust data source and have been used throughout 2020 in France for pandemic assessment and response7,8,13,40.
Mobility data
Mobility reductions shown in Figure 1 were extracted from two different data sources. Overall mobility was reconstructed from mobile phone data provided by Orange Business Service Flux Vision12,31. Data included origin-destination travel flows of mobile phone users among 1,436 geographical areas in France. Each area corresponds to a group of municipalities, defined according to the 2018 EPCI level (Établissements Publics de Coopération Intercommunale53). Mobility reduction in a given week was computed as the relative variation of the number of trips with respect to the pre-pandemic baseline. Estimated presence at workplaces was obtained from Google Mobility Reports14. This dataset provides the relative change in the daily number of visitors to places of work compared to a pre-pandemic baseline, based on Google location-history data.
Indicators of social distancing, risk perception, mental health
Several initiatives collect data over time through surveys to explore individual behaviours in response to COVID-19 pandemic. Here we use data from YouGov15 and Santé publique France16. Surveys gather self-reported data, tracking compliance with preventive measures (e.g. avoiding social gatherings or contacts with other people, frequency of the use of masks), as well as risk perception and mental health indicators (e.g. fear to contract the virus, anxiety, depression). Indicators for specific social distancing behaviors (avoiding gatherings, use of masks) are used in addition to mobility data described above. YouGov surveys cover multiple countries and provide data at least every two weeks. Santé publique France polls collect data at the national level at least every month.
Indicators obtained from surveys report that fear of contracting COVID-19 showed an overall decrease over time after the second wave in France, whereas prevalence of anxiety in the population showed an increasing tendency. Performing a linear regression in this time window (i.e. October 2020 - April 2021), we found a weekly average reduction of −0.31% for individuals scared to contract the virus, and −0.39% for individuals avoiding crowded places; and an average weekly increase of +0.13% in the prevalence of anxiety in the population. Fear of contracting COVID-19 showed a positive correlation with the behavior of avoiding crowded places (Pearson r=0.71, p<10−4, in the time period from w40, 2020 to w15, 2021, Figure 1; results are robust when extending the timeframe of analysis). We observed a non-significant association between the prevalence of anxiety in the population and adoption of social distancing (Pearson r=0.2, p=0.46, in the time period from w11, 2020 to w11, 2021).
SARS-CoV-2 two-strain transmission model
We used a stochastic discrete age-stratified two-strain transmission model, integrating data on demography, age profile, social contacts, mobility and adoption of preventive measures. The model accounts for the co-circulation of two strains and vaccination. Four age classes are considered: [0–11), [11–19), [19–65) and 65+ years old (children, adolescents, adults and seniors respectively). Transmission dynamics follows a compartmental scheme specific for COVID-19 where individuals are divided into susceptible, exposed, infectious, hospitalized and recovered. The infectious phase is divided into two steps: a prodromic phase (Ip) and a phase where individuals may remain either asymptomatic (Ias, with probability pa=40%54) or develop symptoms. We distinguished between different degrees of severity of symptoms, ranging from paucisymptomatic (Ips), to individuals with mild symptoms (Ims), or severe symptoms (Iss) requiring hospitalization40,55. The duration of the infectious period was computed from the estimated mean generation time of 6.6 days56. Prodromic, asymptomatic and pauci-symptomatic individuals have a reduced transmissibility57. A reduced susceptibility is considered for children and adolescents, along with a reduced relative transmissibility for children, based on available evidence58–61. We assume that infectious individuals with severe symptoms reduce of 75% their number of contacts because of the illness they experience. Sensitivity analysis on the probability of being asymptomatic, the susceptibility of younger age classes and transmissibility of children was performed in previous work7,13,62.
Contact matrices are parametrized over time to account for behavioral response to social distancing interventions and adoption of preventive measures. Contacts at school, work and on transports are considered according to the French school calendar, school closures, and presence at workplaces estimated by Google. Physical contacts are reduced based on data from regular large-scale surveys conducted by Santé Publique France7. Contacts engaged by seniors are subject to an additional reduction of 30%, to account for evidence of a higher risk aversion behavior of the older age class compared to other age classes7.
B.1.1.7 variant
Genomic and virological surveillance to identify specific mutations are in place in France since the start of 2021 to monitor variants over time. The first large-scale genome sequencing initiative (called Flash1 survey) was conducted on January 7-8 and analyzed all positive samples provided by participating laboratories9. The proportion of B.1.1.7 in Île-de-France was estimated to be 6.9%, compared to the national estimate of 3.3%, making Île-de-France the region with the highest penetration registered in the country. Flash surveys are performed on average every two weeks, currently on a sample of sequences. Starting week 6 a new protocol for virological surveillance was implemented to provide more timely estimates on the weekly frequency of detected viruses with specific mutations. It was based on second-line RT-PCR tests with specific primers that allow the detection of the main mutations that characterize the variants of concern. They must include at least the N501Y mutation and allow to distinguish the 20I/501Y.V1 variant from the 20H/501Y.V2 or 20J/501Y.V3 variants.
We considered the co-circulation of B.1.1.7 variant together with the historical strains, assuming complete cross-immunity. An increase in transmissibility of 59% (95% credible interval: 54-65%)9 was considered for B.1.1.7 variant compared to the historical strains. This early estimate was obtained from the Flash1 and Flash2 survey in France, and it is in line with other estimates1,2. To account for uncertainty in the transmission advantage and possible changes due to restrictions, we also show for sensitivity (Supplement) the results assuming 40% of increase in transmissibility, i.e. the lower estimate provided by Ref.2. We considered a 64% increase in hospitalization rates, following evidence of an increased risk of hospitalization after infection due to B.1.1.7 compared with other lineages (adjusted OR of 1.64 (95%CI, 1.32-2.04))10, in line with recent estimates11. B.1.1.7 variant was initialized in the model on January 7, 2021 using the estimates of the first large-scale nationwide genomic surveillance survey (Flash1).
Vaccination rollout campaign
Administration of vaccines was included in the model according to the vaccination rhythm adopted in France starting January 2020. We considered the administration of 100,000 doses per day (including first and second doses) at the national level from the end of January (w04), accelerated to 200,000 first doses per day starting the beginning of March (w10), and 300,000 first doses per day starting April (w13). Rollout plans were expressed in terms of first administrations from March on to follow the objectives of authorities, delaying the administration of the second dose to reach a higher coverage in a smaller timeframe, as done in the UK. Higher vaccination paces (400,000 to 800,000 doses/day) were also tested. Paces are defined at the national level, and the number of doses is proportionally distributed to the region according to the population eligible for the vaccine. Vaccination is prioritized to the older age class, assuming 80% coverage, and then shifted to adults considering 50% coverage, according to surveys on vaccine hesitancy50. Vaccination to healthcare personnel and patients in long-term care facilities, performed at the start of the vaccination program, could not be explicitly included.
We considered 75% vaccine efficacy against susceptibility17 and 65% vaccine efficacy against transmission63, estimated after the first injection. We further considered 80% vaccine efficacy against symptoms given infection, computed from the estimated vaccine reduction of symptomatic disease63,64 estimated at 95% after the second dose, and found to be similar after the first dose19. As the landscape for vaccine efficacy rapidly evolves, we also tested vaccine efficacy against transmission equal to 40%65. We assumed efficacy to start 3 weeks after the first injection, and tested a delay of 2 weeks for sensitivity.
Inference framework
The model is fitted to daily hospital admission data through a maximum likelihood procedure, by estimating the transmission rate in each pandemic phase. More precisely, prior to the first lockdown and in absence of intervention (period January-March 2020), we estimated {β,t0}where β is the transmission rate per contact and t0 is the date of the start of the simulation. Then, in each phase we estimated αphase, i.e. the scaling factor of the transmission rate per contact specific to the pandemic phase under study (e.g. lockdown, exit from lockdown, summer, start of second wave, second lockdown, etc.). The transmission rate per contact in each phase is then defined as the transmission rate per contact in the pre-lockdown phase β multiplied by the scaling factor αphase. A pandemic phase is defined by the interventions implemented (e.g. lockdown, curfew, and other restrictions) and activity of the population (school holidays, summer holidays, etc.). The likelihood function is of the form
where Θ indicates the set of parameters to be estimated, Hobs(t) is the observed number of hospital admissions on day t, Hperd(t) is the number of hospital admissions predicted by the model on day t, Poiss (· |Hpred(t))is the probability mass function of a Poisson distribution with mean Hpred(t), and [t1, tn] is the time window considered for the fit.
For Île-de-France, we seeded the model with 140 infected individuals to reduce the strong fluctuations associated with fitting the rapid increase and the high peak of hospitalizations observed in the first wave (the region was one of the areas mostly affected by the epidemic in early 2020). Simulations progress throughout 2020 to build immunity in the population. The model was validated against the estimates of three independent serological surveys conducted in France7. We used 250 stochastic simulations to compute median values and associated 95% probability range for all quantities of interest.
First lockdown, second lockdown, curfew
French authorities implemented two national lockdowns in 2020 to face the rapid surge of COVID-19 cases observed in the first and second wave. The first lockdown started on March 17 and lasted 8 weeks. It involved strict mobility restrictions outside home, together with closure of schools and non-essential activities. The associated effective reproductive number was estimated to be 0.73 [95% confidence interval: 0.72, 0.74]. The delay from the date of implementation of lockdown and the peak of hospitalizations was estimated to be 9 days in the region, and varied between 7 and 12 days across regions7. In our scenarios we assumed a 7-day delay, and tested 10 days for sensitivity. The second lockdown started on October 30 and lasted till mid-December, with a progressive lift of measures while approaching Christmas holidays. It was a milder lockdown with schools in session, corresponding to an estimated effective reproductive number of 0.88 [95% CI: 0.86, 0.90] in the first few weeks of implementation, before relaxation occurred (see next subsection). The second lockdown was lifted in mid-December with the application of a curfew starting at 8pm, then anticipated in January to 6pm to face increasing SARS-CoV-2 spread. The resulting effective reproductive number was estimated to be 0.90 [95% CI: 0.86–0.93] for the historical strains and 1.43 [95% CI: 1.37–1.48] for the B.1.1.7 variant at the end of January8. Starting March 20, strengthened measures were additionally put in place in the region of Île-de-France to curb the third wave. These measures included mobility restrictions for trips exceeding 10 kilometers, closure of business and of schools (1 week for primary schools, 2 weeks for middle and high schools in addition to 2-week school holidays in April).
Loss of adherence
Analysis of mobility data during the second lockdown and estimates of its efficacy over time show that adherence to adopted policy waned over time. We compared the estimated mobility reduction and reproductive number in the first three weeks of implementation (w45-47, November 2 – November 22) with respect to the following week. We considered the average over the three-week period to smooth out the effect of the national holiday on November 11, altering mobility and presence at work with respect to a regular week. We observed the mobility reduction with respect to the pre-pandemic phase to loosen from 42.6% in the first three weeks to 34.3% in the fourth week of the lockdown, corresponding to a relative change of 19%, and associated to a relative increase of 10.9% in the reproductive number. The epidemic conditions estimated in w48, while under the second lockdown and after adherence was partially lost, are equivalent to the ones estimated at the end of January 2021, and corresponding to curfew conditions (previous subsection).
In this study, we modeled the loss of adherence throughout interventions with a 10.9% relative increase in the reproductive number, according to estimates from the second lockdown. We applied it after 2 weeks from implementation of interventions (to model a faster dynamics of adherence waning compared to the one observed in the second lockdown), and considered it limited in time (one drop) or continuous (repeated drops every two weeks). Continuous loss of adherence under a moderate lockdown scenario would become less efficient than curfew for long enough interventions, because of the small difference between the estimated efficacies of the second lockdown (before relaxation) and curfew conditions.
Distress index
In order to quantify the infringement on individual freedom associated with lockdowns and provide a measure of the policy impact on the quality of life, we introduced a quantity called distress index. This measure takes into account both the duration and the intensity of restrictions. It is defined as the sum of the absolute values of weekly mobility reductions, over the number of weeks in which each restriction is maintained, and normalized to a scale from 0 to 10 (10 representing a strict 8-weeks lockdown and 0 the absence of restrictions). In case of a strict or moderate lockdown without loss of adherence, we considered the mobility reductions of 68.9% or 42.6% recorded during the two interventions in 2020, respectively, and varied durations from 2 to 8 weeks. Loss of adherence is computed with a 19% variation of the mobility after two weeks (limited loss) and repeated every two weeks (continuous loss), according to estimates from the second lockdown. We took the end of January 2021 (w04) as reference for the mobility reduction associated with curfew. The resulting distress index is reported in Table 1.
Distress index corresponding to curfew conditions lasting 2 to 8 weeks is shown for comparison.
Seasonality
Multiple studies have investigated the relationship between SARS-CoV-2 transmission and weather factors, including temperature, humidity, ultraviolet radiation22, suggesting that summer conditions may help in reducing transmission of the virus. Seasonal factors and simultaneous social distancing interventions are difficult to disentangle, however containment measures are estimated to have a larger impact on the epidemic compared to seasonal effects only. Considering the estimated dependence of the reproductive number on UV radiation45 and temperature46, we extracted data on downward UV radiation at the surface and daily temperature recorded in Paris, in Île-de-France, in the last three years (2018-2020)66. An average increase of 23 kJ/m2 in UV radiation and of 11 °C in temperature were registered from March to June. These would correspond to estimated 6.7% and 7.3% reduction in viral transmission, respectively.
Data Availability
The mobility data supporting the findings of this study were available to authors from the Orange Business Service Flux Vision within the framework of the research project ANR EVALCOVID-19 (ANR-20-COVI-0007). Restrictions apply to the availability of these data, which were used under license, and so are not publicly available. All other data are available online at the references cited.
DATA AVAILABILITY
The mobility data supporting the findings of this study were available to authors from the Orange Business Service Flux Vision within the framework of the research project ANR EVALCOVID-19 (ANR-20-COVI-0007). Restrictions apply to the availability of these data, which were used under license, and so are not publicly available. All other data are available online at the references cited.
FUNDING STATEMENT
This study is partially funded by: ANR projects DATAREDUX (ANR-19-CE46-0008-03), EVALCOVID-19 (ANR-20-COVI-0007); EU H2020 grants MOOD (H2020-874850) and RECOVER (H2020-101003589).
Footnotes
↵⋄ co-first authors
References
