What pushed Israel out of herd immunity? Modeling COVID-19 spread of Delta and Waning immunity

Following a successful vaccination campaign at the beginning of 2021 in Israel, where approximately 60% of the population were vaccinated with an mRNA BNT162b2 vaccine, it seemed that Israel had crossed the herd immunity threshold (HIT). Nonetheless, Israel has seen a steady rise in COVID-19 morbidity since June 2021, reaching over 1,000 cases per million by August. This outbreak is attributed to several events that came together: the temporal decline of the vaccine's efficacy (VE); lower efficacy of the vaccine against the current Delta (B.1.617.2) variant; highly infectiousness of Delta; and temporary halt of mandated NPIs (non-pharmaceutical interventions) or any combination of the above. Using a novel spatial-dynamic model and recent aggregate data from Israel, we examine the extent of the impact of the Delta variant on morbidity and whether it can solely explain the outbreak. We conclude that both Delta infectiousness and waning immunity could have been able to push Israel above the HIT independently, and thus, to mitigate the outbreak effective NPIs are required. Our analysis cautions countries that once vaccines' will wane a highly infectious spread is expected, and therefore, the expected decline in the vaccine's effectiveness in those countries should be accompanied by another vaccination campaign and effective NPIs

. CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. replacing all other variants of concern (VOC) [1]. The rapid dominance of Delta is attributed to its high infectiousness, and is predicted to have a basic reproduction number, R 0 , of 3.2 to 8, with a mean of 5.08, as opposed to Alpha, which was estimated to be approximately 4 [3]. A recent study confirmed that the viral load of individuals infected with Delta is significantly higher than Alpha [4].
Vaccination roll-out in Israel of the COVID-19 vaccines started on December 20, 2020. By the end of February 2021, over a third of the whole population and over 80% of adults over the age of sixty, were already fully vaccinated by the BNT162b2 vaccine developed by BioNTech and Pfizer [5]. This successfully conducted vaccination campaign resulted in a drop in cases in Israel from over 1,000 per million per day to less than 3 cases per million per day in May 2021. By May, Israel was thought to have reached the herd immunity threshold (HIT), which prompted the lifting of non-pharmaceutical interventions (NPIs), including halt of mandated masking by the end of that month.
Nevertheless, a new outbreak began in mid-June, first, with dozens of infections coming from abroad, triggering some outbreaks within local communities that went on to spread across the country. By the end of August, Israel went back to over 1,000 confirmed cases per million per day. This outbreak, from its very beginning, was completely dominated by infections with the Delta variant. Additionally, since the beginning of the outbreak, a high increase in breakthrough infections in vaccinated individuals was reported. Later, the waning of immunity was confirmed in multiple studies worldwide, showing that the BNT162b2 vaccine effectiveness (VE) is reduced to below 50% after five months of the second dose [6][7][8][9]. Finally, the Israeli government reinstated some NPIs in late June, including indoor masking, green passports and mandatory isolations.
What caused this surge to start in Israel in June 2021 or what propelled Israel above the HIT is a hot topic of debate. In the HIT formula, both the R 0 and the fraction of protected individuals are taken into account [10]. Consequently, both the higher infectiousness of Delta and the reduced effectiveness of the vaccine in those vaccinated in January and February may have pushed Israel above the HIT. Using a novel spatialdynamic model of disease spread, which has shown high success in modeling the spread in Israel throughout the COVID-19 pandemic [11,12], we compared models of infections that combine an R 0 of Alpha or Delta and decline and no decline in the effectiveness. Our analysis suggests that the current outbreak in Israel is caused by a combination of both the Delta variant being more contagious than the Alpha variant and the . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. (which was not certified by peer review) preprint The copyright holder for this this version posted September 15, 2021. ; https://doi.org/10.1101/2021.09.12.21263451 doi: medRxiv preprint decline of VE as time passes. However, we find that even in the presence of the Alpha variant, the decline in VE would have pushed Israel above the HIT. Moreover, we found that the Delta variant itself would push Israel above the HIT even without waning immunity, but with a different age mixture of confirmed cases compared to the scenario in which there is waning in the efficacy but the dominant variant is Alpha.

METHODS
In this work, we use a diffusive Monte-Carlo (MC) model introduced by De-Leon and Pederiva [11,12] to model the spread of the disease. We used this model to examine the expected effects of vaccinations on both confirmed cases (CC) and severe hospitalizations, as well as to estimate the vaccine's current efficacy.
As opposed to other infection models, such as the Susceptible Infected Removed (SIR) [13][14][15], this particle model enables us to distinguish between different age groups and treat each one separately, assuming that the infection occurs throughout the population simultaneously. Furthermore, a particle model can be adjusted to adapt the rate of population immunization to accurately examine the different effects of the vaccine on subgroups of the vaccinated population as well as the entire population.
It is necessary to differentiate between R t (the theoretical replication number of the virus) and R e (the effective replication number of the virus) in order to estimate the spread of a pandemic in the presence of vaccines. R t estimates the number of encounters between carriers and healthy individuals that would have resulted in an infection if the vaccines had not been available and is defined by how contagious the current variant is, whereas, R e is affected by vaccination rates and protection offered by the vaccines. In this work, based on the easing of social restrictions in Israel in February 2021 and Delta's increasing dominance since June, we estimate the theoretical R t in Israel in July 2021 at 3 for the Delta variant and 2.2 for the Alpha variant (we perform the modelling on additional R t levels). It is important to note, that in contrast to R 0 , the basic reproduction rate, R t takes into account the NPIs that are still customary in Israel, especially the duty of isolation for verified patients and those who have been exposed to a verified patient and the duty of wearing masks indoors.
Recent studies from Israel, Qatar, UK and [6,9,[16][17][18] have shown reduction in effectiveness in preventing COVID-19 infection and transmission blocking in individuals vaccinated more than 5 months before exposure. While the exact effectiveness after this waning humoral immunity, reports from Israel have estimated the effectiveness of the vaccine against infection to be between 20%-50% for those vaccinated six months ago, compared to 80%-90% for those vaccinated recently.
Data on deaths and vaccination rates was download on September 1st from OurWorldInData.org [? ].
Only countries with a population over 1 million and with available information are presented.
. CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. (which was not certified by peer review) preprint

RESULTS
After a wide vaccination campaign in Israel, the number of new confirmed and hospitalized cases declined from mid-January 2021, first among adults over 60 years of age, and later in the entire population [19,20]. After almost complete elimination of cases, a new outbreak began in mid-June 2021. We applied our spatial-dynamic MC algorithm with information on current vaccination rates, COVID-19 recoveries, and the number of confirmed cases in Israel from July 2020 to June 2021 to model difference scenarios of the expected daily confirmed cases in the outbreak that started in June 2021 (Figure 1). We first ran our Interestingly, due to the fact that in the Waning model the reduced VE is expected first in the elderly population who were vaccinated first, the number of CC in 60 year old and over individuals is predicted to be higher in the Waning model compared to the Delta model (250 cases vs. 150 cases in mid-August) ( Figure 1b). In the Delta+Waning model, this number grows to 700 cases (in the real-world there were 716 cases). By examining the age distribution of cases at the beginning of the outbreak, we found that our prediction matched what actually occurred. While during the winter outbreak the fraction of cases among adults 60 and over was significantly lower than their fraction in the total population ( 10% of cases vs. 15% of the total population, One-sample proportional test p-value = 0.009) due to the higher vaccination rates in older populations, in the June outbreak we observed an increased percentage in the adult population that were affected ( 17% vs 15%, one-sample proportional test p-value = 0.7), more or less as expected from random spread in the population (Supplementary Figure 2).
. CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.  (Figure 2).
. CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. .

DISCUSSION
In this study, we used a spatio-dynamic Monte Carlo algorithm to estimate the effect of the Delta variant on the COVID-19 outbreak that started in Israel in June 2021. We compared four different models with a combination of varying VE and theoretical reproduction rate. Based on the analysis, both the introduction of Delta and its higher reproduction rate as well as a time-decline in VE can explain this outbreak. Interestingly, we observe that the time-decline VE is responsible for the fact that the first to be infected are those vaccinated first, which are most prominently older populations, who are also at higher risk for severe disease. It is therefore recommended that counties that have vaccinated their elderly population in recent months monitor not just the number of confirmed cases, but also their age distribution and immunization status, and be aware of when this population will lose some of its immune defenses. Such monitoring will allow to detect in real-time the protective attenuation of the vaccine and provide this population another dose of the vaccine on time, which according to recent studies, has been found to be very effective against the delta variant [21,22].
. CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. (which was not certified by peer review) preprint It is important to keep in mind that our models attempt to predict events that did not appear in reality.
In the real world all the parameters of the model are inter-dependent. Increased reproduction rates might cause the government to impose restrictions sooner and cause a reduction in social interactions. A high reproduction rate may also result in more people getting vaccinated. As such, our analysis focuses on factors contributing to the June outbreak in Israel and its global implications, not on the prediction of cases in alternative scenarios.
Interestingly, when we examine the effect of R t on the spread of the pandemic (Supplementary Figure 1) we find that in case of decline in the effectiveness of the vaccine against infection (which is the case today) the theoretical reproduction number, R t should be relatively low (R t < 2), in order not to push Israel above the HIT, which could be the case with a less contagious variant and/or presence of effective NPIs.
In summary, we argue that both the Delta variant higher infectiousness and the reduction of effectiveness of the vaccines contributed to pushing Israel above the HIT. If this is the case, the booster shots provided to all Israelis since August, are expected to reduce severe cases and the burden on hospitals, but not to bring Israel back to herd immunity without reinstating effective NPIs.  . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. (which was not certified by peer review) preprint The copyright holder for this this version posted September 15, 2021. ; https://doi.org/10.1101/2021.09.12.21263451 doi: medRxiv preprint