Abstract
Given that mass gathering events involve heterogeneous and time-varying contact between residents and visitors, we sought to identify possible measures to prevent the potential acceleration of the outbreak of an emerging infectious disease induced by such events. An individual-based simulator was built based on a description of the reproduction rate among people infected with the infectious disease in a hypothetical city. Three different scenarios were assessed using our simulator, in which controls aimed at reduced contact were assumed to be carried out only in the main event venue or at subsequent additional events, or in which behavior restrictions were carried out among the visitors to the main event. The simulation results indicated that the increase in the number of patients with COVID-19 could possibly be suppressed to a level equivalent to that if the event were not being held so long as the prevalence among visitors was only slightly higher than that among domestic residents and strict requirements were applied to the activities of visitors.
1. Introduction
Mass gathering events, such as the Olympic games, world expositions, and music festivals, involve large-scale mixing of people from different regions for extended durations (e.g., durations of a few days to several months). Motivated by the Tokyo 2020 Olympic Games, which were eventually held at the end of July 2021 amidst the coronavirus disease 2019 (COVID-19) pandemic, we aimed to understand the degree to which the reproduction rate can be accelerated by such heterogeneous mixing or be suppressed by behavioral changes from a dynamical perspective. Although an environmental exposure model can quantitatively provide infection risks for various situations (Jones, 2020; Nicas & Sun, 2006) and be applied to mass gathering events (Murakami, 2020), tracking potential infection events using a mechanistic model could be more effective for devising policies and allowing the general population to understand how such interventions work. For such a mechanistic description, an individual-based (or agent-based) model is known to be useful (Riley, 2007).
The application of agent-based models can be traced back to problems involving the containment of a potential pandemic caused by the novel influenza A virus subtype H5N1 in southeast Asia (Ferguson et al., 2005; Longini et al., 2005). Ferguson et al. (2020) also adjusted their framework to the COVID-19 epidemic in the UK and US in the early stage when non-pharmaceutical interventions were the only available measures and vaccines had yet to be developed. Their simulation results suggested that flattening the epidemic curve was not sufficient to avoid a surge of patients exceeding the capacity of hospitals, and thus, extensive non-pharmaceutical interventions combining home quarantine and social distancing should be administrated while closely monitoring their impacts on the economy and social activities, to reduce the effective reproduction number (Rt) below unity. Referring to mobility and census data, Aleta et al. (2020) constructed a detailed agent-based model of transmission in the Boston metropolitan area, and found that social distancing followed by testing and contact tracing could keep the number of cases below the limit of medical capacity. An individual-based model is considered more advantageous than a compartmental one (e.g., the susceptible, infected, and removed or recovered [SIR]-like model) when an infectious agent spreads among different communities, as it is highly stochastic in the early stage of the epidemic and keen to provide a correct description of the transmission dynamics. Son et al. (2020) successfully predicted the number of cases that would arise in Daegu, South Korea by using an agent-simulation model, where the majority of reported cases at that time (end of March 2020) were tied to a single church and the transition from the church to other communities was a key factor in the prediction.
Individual-based simulation has traditionally been used to target problems in epidemics involving the mixing of heterogenous populations and time-varying contact rates. Assessing the influence of a mass gathering event on an epidemic falls into this type of problem, where infectious cases from outside may have a non-small effect on the Rt even among local residents. In particular, the total number of visitors and the prevalence of infection among them, which may be comparable to the residents in the region, are key factors in such predictions. For this reason, we built an individual- (or agent)-based model with inflow and assessed how the growth of new cases could be controlled by different interventions under different prevalence rates among visitors.
2. Methods
2.1 Individual-based model
A mass gathering event yields additional gathering occasions, called secondary gatherings, that involve local and outside populations. For an appropriate description of transmission dynamics, the import of infected people between two types of temporarily formed populations should be considered. For this reason, we employed an individual-based description of the epidemic rather than a compartmental one such as the SIR model.
Specifically, we constructed an individual-based model by extending our study of the 2009 influenza pandemic (Saito et al., 2013) to deal with main and secondary gathering events. A simulation implementing the baseline model describes the activity of residents in a hypothetical city consisting of five towns that are tightly connected via commuter train, each of which being a collection of households, companies, schools, and shopping centers. In each local place, susceptible agents suffer infection risk P during the small time period Δt = 1min.
where N and I are the (temporal) numbers of total and infectious visitors, respectively, Rt,is the effective reproduction number in this place, and Tinf = 4day, is the mean infectious period. This city is a simplified model of central Tokyo; a schematic illustration of this city is shown in Figure 1, and the numbers of entities are listed in Table 1. In total, about 1 million residents live in the city and have the chance to visit the main event, which is set to have 50,000 attendees. The timing of attendance is random.
Schematic illustration of the simulated city.
Populations and number of places in the simulation.
Below, we briefly summarize the assumptions regarding what would happen during the event and its application to our simulator. The main mass gathering event consists of sub-venues, each of which randomly selected residents are scheduled to attend, and attendees in different sub-events are mutually segregated. It is likely that the epidemic in concern is not in an uptrend if a mass gathering event is allowable. For this reason, the transmissibility parameters in populations in the simulated city have been configured so that the Rt of the entire city is nearly equal to unity and the number of daily cases is kept almost constant until the event. This calibration was carried out in a heuristic manner because unlike compartmental models, this critical condition is not described in a simple formula. A non-small number of attendees would be expected to visit some other location after the event is over in order to meet with or socialize with one another, thereby creating many secondary gathering events. We further assume that inbound attendees come to the main event multiple times and always join an additional event after the main one, considering that such travelers would necessarily eat out, in contrast to local attendees. The abovementioned assumptions are reflected in the parameter setting shown in Table 2, where the infection risk is quantified in terms of the Rt. In households, schools, and corporations, the risk is small, unlike the 2009 swine flu pandemic. Therefore, the Rt values are set to 1. The attributed Rt value, 1.05, to places other than those for the main and sub-events was chosen heuristically to establish a sustained (non-growing) epidemic. The proportion of the main event attendees who are dropped (e.g., to meet their friends) is provided by a probabilistic distribution.
Parameters defining the activity of people in the simulation.
2.2 Scenarios
To examine the possibility of suppressing the additional impetus to the epidemic brought by visitors, we carried out simulations with the following three scenarios:
StadiumOnly: The control is administered only in the main event venue. In the simulation, the control is implemented by a reduction in the Rt from 25 to 1 (see Table 2).
Extra: In addition to StadiumOnly, secondary events are assumed to be controlled. This is implemented by reducing the proportion of highly risky places (Rt = 240) for secondary events from 10% to 1%.
WithoutDrop: In addition to StadiumOnly, most of the attendees (specifically 90% in the simulation) are assumed to go directly home in accordance with a corresponding governmental campaign. The risk distribution in secondary events is the same as that in the StadiumOnly scenario.
The abovementioned values are arbitrary, but reflect the following assessment. First, in the simulations, we set the latent period to 3.5 days and the infectious period to 3 days, as average. We here assign rather shorter values, considering these periods are exponentially distributed in our model while the reality is close to a gamma or Weibull distribution (Ferguson, 2020). Second, Murakami (2020) assessed the infection risk in the Opening Ceremonies of the Tokyo 2020 Olympic Games as follows. An infectious attendee would yield about 1.5 or 0.01 new cases in uncontrolled and controlled settings, respectively, during a 5-hour-long ceremony, which is translated into the Rt by:
As for the controlled case, we used the threshold value, Rt = 1, instead for this rather small risk. A similar assessment is available for another event, a 4-hour long soccer game: 0.74 or 0.022 cases are reproduced by one infectious case, respectively. A higher risk at some of the secondary gatherings is based on a survey targeting a large indoor convention (Sami, 2022). That study showed that respondents testing positive were more likely to have attended pubs, karaoke bars, or nightclubs. Using their reported values, attendees to these additional activities are at roughly three times higher risk of infection. Further assuming that this tripled risk is experienced within 1.5 hours, we take Rt = 24 × 5/1.5 × 3 = 240 to be the worst case scenario at additional meetings.
3. Results
Whether main mass gathering events accelerate epidemics depends on the baseline setting, the prevalence rate among visitors, and possible counter-measures to help prevent close contact. Figure 2 summarizes the outcome under each setting as an epidemic curve. As a baseline, the transmissibility parameters in respective small populations in the simulated city were configured so that the Rt of the entire city remained nearly equal to unity until the event opened. As was configured, the number of new cases remained almost constant until the main event opened on day 14, and then tended to increase with a different slope until the event closed on day 24.
Simulated epidemic curves for different assumptions on the prevalence rate among visitors and possible measures. Thick lines indicate the mean (blue) and median (red) trajectories, while thin lines correspond to respective quantiles (25th, 75th, 80th, 85th, 90th, and 95th percentiles). The main mass gathering event runs from day 14 through 23 and is bounded by vertical lines.
With the measure against only the event venue (scenario StadiumOnly) and zero infectious visitors, the new cases increased with a slope of Rt = 1.08, as shown at the top left of Figure 2. This increasing trend was sustained when the prevalence of infection among the visitors was mildly high, pinf,vis = 0.1%, but rapidly increased to Rt = 1.33 when the prevalence of infection was extremely high, at p_(inf, vis) = 1%. These results suggest that cases imported from the outside would have only a minor impact if pinf,vis < 0.1%.
Regarding secondary events under scenario Extra, it may be possible to erase the uptrend attributable to the main event when the visitors are fully noninfectious. Infectious visitors may reinstate an uptrend, but our findings suggest that it would be rather weak as long as the prevalence was mild: Rt = 1.04 for p_(inf,,vis) = 0.1%. A recommendation to visitors to return directly home (scenario WithoutDrop) would have a similar effect. A constant trend would be sustained after the start of the event start if no infectious visitors came in, or a shift to a weak uptrend if their prevalence rate was <0.1%.
So far, we have interpreted the impacts of infectious visitors and controls, viewing only the mean epidemic curve trajectories. However, these curves are highly stochastic, as shown in the gray shading in each panel of Figure 2, and the slope of the trend involves a greater uncertainty, as shown in Figure 3. Figure 4 shows the Rt with an uncertainty for all scenarios as a function of pinf,vis. The strongest measure of WithoutDrop suppresses the increase in the Rt based on the prevalence rate of visitors. On average, p_(inf,,vis),≤,0.1%,keeps the Rt almost at unity under this measure.
The effective reproduction number (Rt) in the entire city during the main event in each scenario. The uncertainty in each setting is evaluated from 1,000 runs with a different random number sequence.
Distributions of cumulative cases during the event corresponding to Figure 2. The uncertainty is based on 1,000 simulation runs.
4. Discussion and Conclusion
In this study, we demonstrated how the potential growth of an epidemic can be successfully or unsuccessfully suppressed in a mass gathering event involving a comparable number of local and outside visitors depending on the strength of the intervention and the prevalence rate among visitors. As a general remark, the growth rate potentially enhanced by hosting a mass gathering can be suppressed to a level identical or slightly higher than that of the baseline if the places for gathering after the main event are well controlled and the prevalence rate among visitors is mild relative to that among local residents. Reflecting the lack of qualitatively interpretable data for shops, restaurants, and other kinds of possible meeting places, we carried out simulations with scenario-based settings to assess risk in these places. Below, we discuss the availability of data from the parameters of our simulation.
The Japanese government enacted a strict quarantine at airports during the pandemic. The data collected from the quarantine (MHLW, 2021; Tokyo Metropolitan Government, 2020) could be expected to aid the estimation of prevalence rates among visitors, thereby making quarantine more effective in combination with testing before boarding and after arrival (Kamo et al., 2022).
In addition to the prevalence obtained along with the quarantine, the prevalence rate among the participants of the Tokyo 2020 Olympic Games is available (Tokyo Metropolitan Government, 2021a). Maintaining a sufficient level of control at secondary gathering places is more difficult than the main event venue. The Tokyo Metropolitan Government issued an “Infection Prevention Thorough Declaration Sticker” to shops, restaurants, and other meeting places to establish a certain level of risk reduction and provide information on well-sanitized places to local residents (Tokyo Metropolitan Government, 2021b). Hence, the ratio of safe places to ‘risky places’ (10% or 1% being assumed in the present study) may be determined via data. However, the requirements for this are mainly qualitative, and it is nontrivial to translate these data into the risk of infection computationally interpretable.
Observational facts that contribute to connecting the emission of droplets to the risk of infection have accumulated since the beginning of the COVID-19 pandemic and the abovementioned difficulties may be overcome by plugging these data into simulations. While some of this knowledge has yet to be fully authorized, a number of studies have already been published. Ando et al. (to appear) assessed infection risk at a restaurant using a dedicated hydrodynamical model on the Fugaku supercomputer. The effects of air ventilation in the venue have been monitored based on CO2 concentrations during a football game (Murakami et al., 2021), and droplet emissions have been measured during normal speech (Asadi et al., 2019), coughing (Milton et al., 2013), and singing (Alsveda et al., 2020). In addition, a COVID-19 outbreak event suspected to be induced by air conditioning in a restaurant has been reported (Lu et al., 2020). We expect that a reliable set of infection risk parameters will become available based on the outcomes of empirical studies interpreted inside hydrodynamical droplet simulations, and that individual-based models will contribute to improved risk assessments against new variants of severe acute respiratory syndrome coronavirus 2 and future pandemics caused by other infectious diseases.
Data Availability
All data produced in the present study are available upon reasonable request to the authors
References
