Probable cross-corridor transmission of SARS-CoV-2 due to cross airflows and its control Pan Cheng1, Wenzhao Chen1, Shenglan Xiao2, Fan Xue3, Qun Wang1, Pak Wai Chan4, Ruoyu You6, Zhang Lin5, Jianlei Niu6, and Yuguo Li1,7* This is the peer-reviewed post-print version of the paper: Cheng, P., Chen, W., Xiao, S., Xue, F., Wang, Q., Chan, P. W., You, R., Lin, Z., Niu, J., & Li, Y. (2022). Probable cross-corridor transmission of SARS-CoV-2 due to cross airflows and its control. Building and Environment, 218, 109137. Doi: 10.1016/j.buildenv.2022.109137 The final version of this paper is available at: https://doi.org/10.1016/j.buildenv.2022.109137. The use of this file must follow the Creative Commons Attribution Non-Commercial No Derivatives License, as required by Elsevier’s policy. Graphical abstract 1 Department of Mechanical Engineering, The University of Hong Kong, Hong Kong SAR, China 2 School of Public Health (Shenzhen), Sun Yat-sen University, Shenzhen, China 3 Department of Real Estate and Construction, The University of Hong Kong, Hong Kong SAR, China 4 Hong Kong Observatory, Kowloon, Hong Kong SAR, China 5 Division of Building Science and Technology, City University of Hong Kong, Hong Kong SAR, China 6 Department of Building Services Engineering, The Hong Kong Polytechnic University, Hong Kong SAR, China 7 School of Public Health, The University of Hong Kong, Hong Kong SAR, China *Corresponding author. Yuguo LI – Address: Department of Mechanical Engineering, The University of Hong Kong, Hong Kong SAR, China; School of Public Health, The University of Hong Kong, Hong Kong SAR, China; Email: liyg@hku.hk 1 Highlights • Cross-corridor airflow probably explains a horizontal cluster of COVID-19. • Downstream flats under prevailing wind directions were at highest risk. • Had northerly wind direction prevailed, the outbreak would not occur. • We provide a probable explanation for outbreaks in quarantine hotels. • Positive pressure ventilation in corridor can reduce risk. Abstract 5 10 15 20 A COVID-19 outbreak occurred in May 2020 in a public housing building in Hong Kong – Luk Chuen House, located in Lek Yuen Estate. The horizontal cluster linked to the index case’ flat (flat 812) remains to be explained. Computational fluid dynamics simulations were conducted to obtain the wind-pressure coefficients of each external opening on the eighth floor of the building. The data were then used in a multi-zone airflow model to estimate the airflow rate and aerosol concentration in the flats and corridors on that floor. Apart from flat 812 and corridors, the virus-laden aerosol concentrations in flats 811, 813, 815, 817 and 819 (opposite to flat 812, across the corridor) were the highest on the eighth floor. When the doors of flats 813 and 817 were opened by 20%, the hourly-averaged aerosol concentrations in these two flats were at least four times as high as those in flats 811, 815 and 819 during the index case’s home hours or the suspected exposure period of secondary cases. Thus, the flats across the corridor that were immediately downstream from flat 812 were at the highest exposure risk under a prevailing easterly wind, especially when their doors or windows that connected to the corridor were open. Given that the floorplan and dimension of Luk Chuen House are similar to those of many hotels, our findings provide a probable explanation for COVID-19 outbreaks in quarantine hotels. Positive pressure and sufficient ventilation in the corridor would help to minimise such cross-corridor infections. Keywords: COVID-19; aerosol transmission; ventilation; airflow; quarantine hotel; SARSCoV-2 2 25 30 35 40 45 50 55 60 1. Introduction Several Coronavirus Disease 2019 (COVID-19) outbreaks involving probable airborne transmission of severe respiratory syndrome coronavirus 2 (SARS-CoV-2) in poorly ventilated spaces have been documented and investigated (Li Y. et al., 2021; Miller et al., 2021; Ou et al., 2022). The airborne transmission of SARS-CoV-2 has also been recognised by the World Health Organization since July 2020 (WHO, 2021). Theoretically, if airborne transmission occurs within poorly ventilated indoor spaces, then transmission between two or more directly connected and poorly ventilated small indoor spaces is also possible. In this case, direct connection refers to the direct airflow between separated spaces without significant filtration and deposition during air exchange. Transmission between guest rooms and the corridor or between guest rooms sharing the same corridor has been reported in quarantine hotels in many countries such as Australia (Leong et al., 2021), Canada (CityNews, 2021), New Zealand (Fernando and McPhee, 2020), China and the UK (BBC, 2021), and Hong Kong (Li X. et al., 2021; Wong et al., 2021). By 15 June 2021 Australia experienced 22 COVID-19 outbreaks in quarantine hotels, while New Zealand reported 10 (Grout et al., 2021). An outbreak at Rydges Hotel and another at Stamford Plaza Hotel in Melbourne led to the second wave of infections in Victoria, Australia, with 99% of the over 19,800 cases linked to infections at these two hotels. This wave of infections resulted in more than 800 deaths and a lockdown of the state that lasted 112 days (PHE, 2021). A lockdown in Victoria may cost AU$1 billion a week (Fernando and McPhee, 2020). In Australia and New Zealand, an estimated 5.0 quarantine system failures per 100,000 travellers and 6.1 quarantine system failures per 1,000 positive cases of SARS-CoV-2 were estimated (Fernando and McPhee, 2020). In Hong Kong, one guest in a quarantine hotel that had ventilation that was ‘not up to standard’ caught the virus from two people who were staying in the room that was directly opposite across the corridor (RTHK, 2021). Environmental samples taken from the rooms tested positive for the virus (Jiang et al., 2020), although fomite transmission was deemed to be unlikely (CDC, 2021). The COVID-19 pandemic has seen the emergence of ‘quarantine hotels’, venues specifically dedicated to isolate newly arrived travellers and contain the spread of the virus. Guests of quarantine hotels are typically confined to their rooms for fixed periods (e.g., 14 days or 21 days). During this period, the guests are strictly required to stay in their rooms; they are only allowed to open their doors occasionally to receive deliveries of food and other essential items, as well as provide samples for swab tests. In some countries, no other people are allowed to enter guests’ rooms throughout the entire duration of the stay. As there is no close contact or sharing of items between guests or between guests and hotel workers, an occurrence of transmission would indicate that transmission is airborne. 3 65 70 75 80 85 90 95 100 With the exception of facilities that were specially built for quarantine recently (CNN, 2021), most quarantine hotels were neither designed nor constructed for quarantine purposes. Thus, there is an urgent need to understand the conditions under which the transmission of infections occur and to develop intervention strategies. As with investigations of other airborne outbreaks, quantifying the airflow between rooms and the ventilation rate within rooms is important for revealing the virus-laden aerosol transport and possibility of airborne transmission. However, there is currently limited access to quarantine hotels to conduct such investigations due to the continued use of the facilities for quarantine during the ongoing pandemic, even after outbreaks are identified. Given that COVID-19 outbreaks have been reported in ordinary hotels (e.g., Hoefer et al., 2020), it is unsurprising that SARS-CoV-2 can also be transmitted between guests and workers in quarantine hotels. In 2003, a major outbreak of SARS occurred in a hotel in Hong Kong (CDC, 2003). It resulted in 13 cases of infection, 9 of whom (including the index case) stayed on the ninth floor of the hotel. The infected international travellers travelled elsewhere and initiated outbreaks in numerous other countries, including Singapore and Canada (CDC, 2003). In June 2020, we carried out field measurements of the dispersion of tracer gas in Luk Chuen House in Hong Kong (Wang et al., 2022). A COVID-19 outbreak had occurred here at the end of May, 2020, and a cross-corridor virus transmission was suspected (Figure 1a). Built in 1975, Luk Chuen House is a multistorey public housing building with a floorplan similar to those of typical multistorey hotels, in which separate flats (or rooms in the case of a hotel) are located along both sides of a long corridor (Figure 1b, Bahadori-Jahromi et al., 2017). The outbreak included two spatial clusters (Wang et al., 2022): a horizontal cluster on the eighth floor, including flat 812 where the index case lived and flats 813 and 817 where secondary infections occurred (Figure 1a); and a vertical cluster connected by the building’s shared drainage system (including flats 710, 810, 1012 and 1112). Wang et al. (2022) conducted a detailed investigation of the role of the unusual two-stack drainage system that connected the flats in spreading the virus in the vertical cluster. The potential role of airflow between flats in facilitating cross-corridor virus transmission in the horizontal cluster was discussed in this study. Here we present a focused and detailed investigation of the horizontal cluster in the COVID19 outbreak in Luk Chuen House. Specifically, using multi-zone modelling of airflows driven by winds and the buoyancy in the building, we explain how airborne transmission contributed to the outbreak. We discuss the implications of the findings for virus transmission between guest rooms in quarantine hotels that have a similar floorplan to that of Luk Chuen House, although the corridors are naturally ventilated in Luk Chuen House and mechanically ventilated in most quarantine hotels as discussed later. To achieve the above objectives, we 4 estimate the hourly variation in aerosol concentrations in flats and the corridor on the eighth floor of Luk Chuen House during the suspected exposure period. We used the relevant field measurements and data from the literature to validate the modelling study. 105 2. Methods 2.1 The outbreak 110 115 120 125 130 135 The index patient was a 34-year-old female nightshift worker who worked in a warehouse where fruit and vegetables imported from a European country were labelled. She developed a cough and fatigue on 22 May, 2020 and was later hospitalised on 30 May, 2020. Through contact tracing, two co-workers in the warehouse with earlier symptom onset – on 25 and 27 April, 2020 – were identified. The index case lived in flat 812 of Luk Chuen House (Figure 1a). Between 31 May and 13 June, 2020, SARS-CoV-2 infections were confirmed in nine residents living in Luk Chuen House. Among the flats with infected residents, flat 812 (where the husband and sister of the index case were infected), flats 813 (one senior female infected), and flat 817 (one senior female infected) formed a horizontal cluster. We obtained the hourly weather data recorded at Shatin station from 21 May to 24 May, 2020 from the Hong Kong Observatory (Figure 2). The possible exposure period was suspected to begin from 00:00 on 22 May to 24:00 on 24 May, 2020, when the index case was believed to have a peak viral load (i.e., within 2 days before and after the onset of her symptoms). Due to unfavourable background wind conditions, the date 21 May was excluded from the study period. During the suspected 72-hour exposure period from 22 May to 24 May, 2020, the wind direction fluctuated between a northerly and south-westerly direction. Wind angle in the vertical axis of Figure 2 is the angle between wind incoming direction and the north direction, which ranges from 0 - 360 degree. When the angle is within the range of 22.5 - 67.5 degree, the wind direction is said to be northeasterly. Similarly, 67.5 -112.5 degree, easterly; 112.5 157.5 degree, southeasterly; 157.5 - 202.5 degree, southerly; 202.5 - 247.5, southwesterly; 247.5 - 292.5, westerly; 292.5 - 337.5, northwesterly; 337.5 - 360 and 0 - 22.5, northerly. The wind speed varied between 0.2 m/s and 5.5 m/s, with an average value of 2.07 m/s. Temperature ranged between 24.8 and 29.9 °C, with an average of 27.1 °C. The outbreak occurred during a time when there were no other locally transmitted cases in Hong Kong (Wang et al., 2022). The building has 15 storeys, with residential flats located on 2F–15F. Each floor contains 40 flats, with 20 flats lining each side; this corresponds to 40 columns of flats in the building. A flat is defined here as a residential unit and referred to as ‘flat ##&&’, in which ‘##’ refers to the floor level, and ‘&&’ refers to flat number. For instance, flat 0812 – or its shortened version 812 – refers to the twelfth flat (i.e., flat number 12) on the eighth floor. Each column of flats is referred to as ‘-&& flats’; for instance, ‘-12 flats’ refers to the column of flats sharing the flat number 12 that are aligned vertically across 5 140 145 150 155 160 all floors. Columns of flats that are adjacent to one another are connected by a shared twostack drainage pipe. For example, the -10 flats and the -12 flats, which are adjacent to each other, share a two-stack drainage pipe; this was found to be the likely cause of the vertical spread of the virus from flat 812 (and possibly flat 810) to flats 710, 1012 and 1112 (Wang et al., 2022). For two corridors on each floor, each corridor connects 20 flats, with 10 flats lining on both sides. Each corridor is divided into 10 corridor zones, i.e., corridor A divided into corridor zones 1-10 and corridor B divided into corridor zones 11-20. Thereafter, each corridor zone will be referred to by their zone number. For example, corridor zone 10 is the corridor zone next to the lift lobby and connects flats 819 and 820. (a) (b) Figure 1. Floorplan of the eighth floor of Luk Chuen House showing the flats with infected cases, and a typical floorplan of a hotel. (a) Layout of the eighth floor of Luk Chuen House. Spaces with line shading – that is, the lifts, refuse room and switch room – were not included in the airflow analysis. There are 40 flats (801–840) on the eighth floor, three stairwells (ST1, ST2 and ST3), two corridors (Corridor A and Corridor B, and each corridor is divided into 10 corridor zones ), and one lift lobby. Flat 812, where the index case lived, is shown in red, while flats 810, 813 and 817 where secondary infections occurred are shown in blue. (b) A typical floorplan of a hotel, re-drawn based on the hotel depicted in Bahadori-Jahromi et al. (2017). Both layouts are not drawn to scale. 6 165 170 Figure 2. Hourly temperature, wind speed and wind direction from 21 to 24 May 2020 at Shatin weather station (proximate to Luk Chuen House). Note that each light blue box on the x-axis indicates an hour for which the index case was at home, and one light red box on the wind direction plot indicates an hour in the suspected infection period for the secondary cases in flats 813 and 817 (that is, when the index case was at home and the hourly prevailing wind direction was easterly). 2.2 Computer modelling 175 180 185 We first conducted a CFD simulation to investigate the patterns of airflow around the building under wind directions that possibly occurred during the suspected exposure period. The CFD simulation encompassed a domain that was 2 km long (streamwise direction), 1 km wide and approximately 500 m high; this included Luk Chuen House and its 126 buildings in the surrounding area (Figure 3a). The commercial software ANSYS Fluent was used for CFD simulations. A pressure-based solver was used, and only a steady-state solution was obtained. The thermal buoyancy effect on wind flows was not considered, and the realisable k-epsilon model with standard wall functions was used to model turbulence. The second-order upwind scheme were used for gradient, pressure, momentum, turbulent kinetic energy, and turbulent dissipation rate discretisation. The inlet velocity is an exponential function of height (Equations 1, 2; Hsu et al., 1994). An outflow condition was applied on the downwind boundary, opposite to the inlet boundary. The top and sides of the computation domain were symmetrical. 𝑧𝑧 0.11 (1) 𝑈𝑈(𝑧𝑧) = 𝑈𝑈𝑓𝑓𝑓𝑓𝑓𝑓𝑓𝑓 ∙ �40� when z is less than or equal to 40 meters; 𝑈𝑈(𝑧𝑧) = 𝑈𝑈𝑓𝑓𝑓𝑓𝑓𝑓𝑓𝑓 190 (2) when z is larger than 40 meters, where 𝑈𝑈𝑓𝑓𝑓𝑓𝑓𝑓𝑓𝑓 is the average wind speed during the investigation period multiplied by 1.11. 7 195 200 All buildings were assumed to have sealed walls; this meant that the cross ventilation inside each building was not included in the CFD modelling. In accordance with the wind directions detected during 22-24 May, 2020 by a proximate weather station (Shatin Station) of the Hong Kong Observatory, a total of five major wind directions were simulated: northerly, northeasterly, easterly, south-easterly and south-westerly. Two pressure contours on the surface of Luk Chuen House under an easterly wind are shown in Figures 3b and 3c. We conducted model validation and grid independence tests (Figures S1, S2). 2 The wind coefficients were calculated as (𝑝𝑝 − 𝑝𝑝𝑟𝑟𝑟𝑟𝑟𝑟 )/(0.5𝜌𝜌𝑈𝑈𝑟𝑟𝑟𝑟𝑟𝑟 ) (Confluence, 2021), in which 𝑈𝑈𝑟𝑟𝑟𝑟𝑟𝑟 employs the incoming wind speed at a height of 10 m, (i.e., 1.70 m/s), and 10 m 205 210 215 220 is the height of the anemometer in Shatin station. The predicted distributions of wind pressure coefficients on the external walls of the eighth floor of Luk Chuen House for different wind directions are displayed in Figure 4. These wind pressure coefficients were subsequently used as input data for our multi-zone airflow simulation. The multi-zone airflow program (MIX) calculates airflow rates across the external wall of a building and between each zone inside the building based on mass conservation (Li et al., 2000). Specifically, MIX is based on the concept of mass balance and the modified Bernoulli’s equation and can predict the airflow rate across the envelope and between the flats in a building. When it was used to predict the airflow between different flats in Block E of Amoy Gardens, where a major SARS outbreak occurred in 2003, MIX allowed for interpreting the associations between spatial infection patterns and the dispersion of SARSCoV-containing aerosols (Li et al., 2005). More recently, MIX was used to calculate the concentrations of virus-laden aerosols and to estimate the infection risk via long-range airborne transmission in the first nosocomial outbreak of MERS-CoV in the Republic of Korea (Xiao et al., 2018). Note that Xiao et al. (2018) estimated the infection risk employing Wells-Riley equation, based on the indoor virus-laden aerosol concentration, viral load, and exposure time of the susceptible occupants. In this study, we concern only the relative magnitude of virus-laden aerosol concentration when we use the terms “infection risk” or “exposure risk”, implicitly assume that all the bio-aerosol have equal viral load and all the residents have the same exposure time. 8 (a) (b) (c) 9 225 230 Figure 3. Building layout in CFD simulation and the pressure contour on the surface of Luk Chuen House. (a) The layout of Luk Chuen House and 126 surrounding buildings in the CFD simulation (left: two-dimensional view from the top; right: a three-dimensional view), and the pressure contour on the surface of Luk Chuen House under an easterly wind from (b) the southern view and (c) the northern view. The detailed zone partitioning of the eighth floor of Luk Chuen House for MIX prediction is shown in Figure S3. Besides the wind pressure coefficients, the air temperature and geometry of each zone at the time of exposure are also needed as input of the MIX program. A constant set of air temperatures was used for different zones, that is, 24.15 °C in the flat (27.15 °C in the kitchen and 27.25 °C in the toilet of flat 812), and 25.4–26.15 °C in the corridor, with a 235 difference of 0.1–0.2 °C between neighbouring corridor zones (Table S1). Detailed definitions of the zones and their geometry are described in Supplementary Material S3. 240 245 250 255 The floors, ceilings, and walls inside the building were assumed to be perfectly sealed. We only considered leakage and openings to occur through windows and doors, which were modelled as porous openings uniformly distributed over a relatively large area (Li et al., 2000). The leakage coefficients were set at 0.4% for closed windows and 1.8% for closed doors (Orme et al., 1998; Li et al., 2005). Based on our on-site investigation, there was no exhaust fan in any of the accessible toilets as it was naturally ventilated; there was also no mechanical ventilation in each flat except for the air conditioning at the window between the kitchen and the main flat, which was assumed to be off during the simulation period. Aside from mechanical ventilation, particle deposition is also a mechanism by which airborne droplets are removed. The measured deposition loss rate coefficients (Thatcher et al., 2002) in a fully decorated room and an empty room were employed for residence flats and common areas (Table S4), respectively. By common areas, we mean three stairwells, two corridors (Corridor A and Corridor B), and one lift lobby. Note that the airborne bio-aerosols in this study have a diameter less than 10 microns and are emitted by coughing. The cough frequency, the number of droplets expelled per cough and size distribution of expelled droplets are the same as that in Xiao et al. (2018). We noticed that several residents would leave the doors of their flats (that were connected to the corridor) ajar for cross ventilation. Thus, we evaluated the possible effects of dooropening on airflow in three different scenarios (Table 1). In each scenario, the windows connected to the corridors were assumed to be closed. We denote Scenario 1 [All doors closed], Scenario 2 [813 and 817 ajar] and Scenario 3 [803 and 807 ajar]. 260 10 Table 1 Three simulation scenarios for the eighth floor of Luk Chuen House. Scenario 1 2 3 265 270 275 Description All doors of flats (that are connected to the corridor) are closed Doors of flats 813 and 817 are ajar Doors of flats 803 and 807 are ajar The MIX-estimated inter-zonal flow rates were used to estimate the zonal concentrations of aerosols originating from the index case, based on the mass balance of the aerosols, and assuming that the aerosols were instantaneously and uniformly distributed in each flat and in each corridor zone. The infectious bio-aerosol was released by the index case in her flat, which was composed of a bedroom, a kitchen and a toilet. We assumed that during the exposure period (from 00:00 on 22 May, 2020 to 24:00 on 24 May, 2020), the index case stayed in her flat during her home hours, i.e., between her shifts (05:00–19:30 on 22 and 23 May) and on her day off (05:00–24:00 on 24 May). We also assumed that the rate at which infectious aerosols were generated by the index case was constant when she was home, and zero when she was out of home. We suspected that during the period when the index case was at home, the secondary cases in flats 813 and 817 became infected when the hourly prevailing wind direction was easterly. 3 Results 3.1 Predicted hourly aerosol concentration profiles 280 285 Figure 5 displays the hourly normalised bio-aerosol concentration in flat 812 under the three scenarios, as well as the hourly normalised bio-aerosol concentrations in the corridor zone and the flats opposite and adjacent to flat 812 in Scenario 1. Note that the concentration used for normalisation is the average bio-aerosol concentration in flat 812 during the index case’s home hours in Scenario 1, and that all the normalised bio-aerosol concentration profiles are presented as percentages. Figure 6 illustrates the net airflow through the internal openings of the 68 zones (flats and corridor zones) that were investigated, according to five different wind directions. 11 Flat Window Fan N NE E SE SW N NE E SE SW Window 290 840 1.24 -0.07 -0.58 -0.31 -0.32 1.34 -0.10 -0.63 -0.42 -0.31 838 0.98 -0.05 -0.52 -0.18 -0.32 0.91 -0.05 -0.51 -0.17 -0.32 836 0.75 -0.07 -0.49 -0.14 -0.31 0.82 -0.06 -0.50 -0.15 -0.31 834 0.55 -0.10 -0.47 -0.13 -0.29 0.48 -0.10 -0.47 -0.12 -0.29 832 0.35 -0.11 -0.44 -0.11 -0.29 0.40 -0.11 -0.45 -0.12 -0.29 830 0.22 -0.11 -0.40 -0.10 -0.29 0.20 -0.10 -0.40 -0.10 -0.29 828 0.12 -0.07 -0.36 -0.09 -0.29 0.15 -0.09 -0.37 -0.09 -0.29 826 0.03 0.00 -0.30 -0.08 -0.29 0.01 0.05 -0.29 -0.07 -0.29 824 -0.01 0.13 -0.24 -0.06 -0.29 -0.00 0.09 -0.26 -0.07 -0.28 822 -0.03 0.24 -0.18 -0.04 -0.29 -0.02 0.27 -0.16 -0.04 -0.29 Lba 0.05 0.30 -0.10 -0.03 -0.29 - 820 -0.80 0.26 -0.14 -0.11 -0.29 -0.99 0.27 -0.20 -0.14 -0.30 818 -0.33 0.21 -0.05 -0.07 -0.29 -0.30 0.20 -0.05 -0.06 -0.29 816 -0.25 0.17 -0.03 -0.06 -0.29 -0.26 0.17 -0.03 -0.06 -0.29 814 -0.21 0.16 -0.00 -0.06 -0.28 -0.22 0.16 0.00 -0.07 -0.28 812 -0.19 0.15 0.01 -0.07 -0.28 -0.20 0.15 0.01 -0.07 -0.28 810 -0.19 0.14 0.02 -0.08 -0.27 -0.19 0.15 0.02 -0.09 -0.27 808 -0.19 0.13 0.02 -0.10 -0.27 -0.19 0.14 0.02 -0.09 -0.27 ST3a 1.38 -0.12 -0.65 -0.23 -0.31 ST3b -1.25 -0.22 -0.47 -0.21 -0.39 Flat Fan dir. N NE E SE SW N NE E SE SW dir. N NE E SE SW N NE E SE SW 806 -0.19 0.09 0.02 -0.13 -0.26 -0.20 0.09 0.02 -0.14 -0.26 804 -0.21 0.06 0.00 -0.16 -0.26 -0.20 0.06 0.01 -0.15 -0.26 ST1a -0.37 -0.26 0.25 -0.15 -0.24 ST1b -0.28 -0.22 -0.37 0.06 -0.29 839 -0.60 -0.20 -0.47 -0.20 0.16 -0.60 -0.20 -0.47 -0.22 0.06 837 -0.61 -0.20 -0.48 -0.14 0.27 -0.62 -0.20 -0.48 -0.13 0.30 835 -0.63 -0.20 -0.48 -0.13 0.31 -0.63 -0.20 -0.48 -0.13 0.29 833 -0.64 -0.19 -0.47 -0.12 0.31 -0.64 -0.19 -0.47 -0.12 0.31 831 -0.64 -0.19 -0.47 -0.11 0.30 -0.64 -0.19 -0.47 -0.11 0.32 829 -0.62 -0.19 -0.48 -0.10 0.30 -0.61 -0.19 -0.48 -0.10 0.29 827 -0.59 -0.19 -0.48 -0.10 0.28 -0.60 -0.19 -0.48 -0.10 0.31 825 -0.53 -0.19 -0.48 -0.09 0.28 -0.52 -0.19 -0.48 -0.09 0.27 823 -0.48 -0.19 -0.48 -0.09 0.26 -0.49 -0.19 -0.48 -0.09 0.27 821 -0.42 -0.18 -0.48 -0.08 0.25 -0.42 -0.18 -0.48 -0.08 0.25 Lbb -0.30 -0.18 -0.47 -0.06 0.17 - 819 -0.23 -0.18 -0.47 -0.04 0.12 -0.25 -0.18 -0.47 -0.04 0.15 817 -0.20 -0.18 -0.47 -0.03 0.09 -0.18 -0.18 -0.47 -0.04 0.08 815 -0.17 -0.18 -0.46 -0.03 0.06 -0.18 -0.18 -0.46 -0.02 0.07 813 -0.15 -0.18 -0.46 -0.02 0.03 -0.14 -0.18 -0.45 -0.02 0.01 811 -0.14 -0.18 -0.45 -0.02 -0.01 -0.14 -0.18 -0.45 -0.01 0.01 809 -0.13 -0.19 -0.45 -0.01 -0.04 -0.12 -0.19 -0.45 -0.02 -0.06 807 -0.13 -0.19 -0.44 -0.01 -0.07 -0.13 -0.19 -0.44 0.00 -0.07 805 -0.12 -0.19 -0.46 0.00 -0.10 -0.12 -0.19 -0.46 0.00 -0.12 803 -0.12 -0.20 -0.50 0.00 0.14 -0.12 -0.20 -0.49 0.00 -0.13 802 -0.25 -0.06 -0.04 -0.24 -0.25 -0.27 -0.16 -0.09 -0.22 -0.24 ST2 -0.34 -0.18 -0.48 -0.07 0.21 - Figure 4. Predicted distributions of wind pressure coefficients at the external openings of the eighth floor of Luk Chuen House for different wind directions from the CFD. The wind directions (dir.) N, NE, E, SE and SW represent the northerly, north-easterly, easterly, south-easterly and south-westerly wind directions, respectively. Positive values are shown in grey. 12 801 -0.15 -0.20 -0.63 0.00 -0.19 -0.17 -0.21 -0.68 0.02 -0.25 (a) (b) (c) 13 (d) (e) (f) 295 300 (g) Figure 5. Normalised bio-aerosol concentration profiles. Plots show the normalised bioaerosol concentration profiles (a) in flat 812 under the three different scenarios, (b) in adjacent corridor zones 46–55 (for flats 801-820) in Scenario 1 [All doors closed], (c) in remote corridor zones 58–67 (for flats 821-840) in scenario 1, (d) in the flats opposite to flat 812 in scenario 1, (e) in the flats adjacent to flat 812 in scenario 1, (f) in the flats distant from flat 812 and on the opposite side of the corridor in scenario 1, (g) in the flats distant from flat 812 and on the same side of the corridor in scenario 1. Note that the vertical axes of the seven plots differ in their ranges. Each light blue box and light red box indicates an hour for which the index case was at home and an hour during the period for potential secondary infections 14 305 310 315 320 325 330 335 340 to occur, respectively (as Figure 2). Note that for each scenario (Table 1), we “normalise” the estimated aerosol concentration using the average concentration for Scenario 1 in flat 812 during the period when the index case was at home, and some “normalised” values can be greater than 1. It would not be appropriate to present the absolute aerosol concentration profiles as the strength of the source of the aerosols and its variations are unknown. Hence, we specify a constant source of aerosols (for the index case of flat 812) in our simulation. The aerosol concentration in flat 812 is affected by ventilation, which varies according to background wind speed and direction. For each scenario in flat 812 (Table 1), we normalise the estimated aerosol concentration using the average aerosol concentration for Scenario 1 in flat 812 during the period when the index case was at home (Table S5). Note that the reference for normalization is an average value, instead of a maximum value, which justify that he normalised concentration fluctuates around 100% in flat 812 (Figure 5a). In order for crosscorridor transmission to occur, the infectious aerosols must first escape into the corridor. As shown in Figure 5b, when the wind direction is easterly, the corridor area adjacent to flat 812 (i.e., zones 51–55; note that these are the zones defined in the multi-zone MIX simulation as shown in Figure S3) experiences the highest aerosol concentrations excluding flat 812 (normalised concentrations, 6-12%) because this area is immediately downwind from flat 812. It is interesting that corridor zones adjacent to flat 812 remain uncontaminated or experience far lower levels of contamination when the wind blows in other directions. This finding confirms our suspected infectious periods. When the wind blows in other directions, such as in a south-easterly direction, the adjacent corridor zones 46–50 may also be contaminated, but at lower levels. This explains why flats 811, 813, 815, 817, and 819 experienced the second highest aerosol concentrations (approximately 0.4–1.2%), while flats 801, 803, 805, 807, and 809 had lower concentrations (i.e., less than 0.2% or substantially lower) (Figure 5d). Flat 811 is directly opposite to flat 812 and therefore reports the highest concentrations among all the flats that are opposite to flat 812. The corridor areas distant from flat 812 are also contaminated when the wind is easterly. However, these areas record lower concentrations of aerosols (mostly less than 0.6%, Figure 5c). Thus, the aerosol concentrations in the flats that are distantly located from flat 812 are very low (mostly less than 0.06%) (Figure 5f, 5g). The hourly aerosol concentration in each flat is a result of the pattern of airflow as well as the ventilation conditions. For example, between 14:00 and 16:00 on 22 May, the wind direction is easterly (Figure 2). Accordingly, the aerosol concentration profiles for flats 811, 813, 815, 817 and 819 increase as the contaminated airflow from flat 812 is directed into these flats 15 345 350 355 360 365 through the corridor (Figure 6c). In the meantime, there is no inflow of contaminated air from the corridor for flats 801, 803, 805, 807 and 809; hence, the bio-aerosol concentration in each of these flats is almost zero. However, between 11:00 and 12:00 on 23 May, the wind direction is north-easterly (Figure 2), and the inter-zonal airflow in corridor A is reversed and the aerosol concentration in zones 46-50 increases (Figure 6b). Due to a lower inflow rate, the peak profile for flats 801, 803, 805, 807 and 809 during these hours is not as high as that for flat 811, 813, 815, 817 and 819 between 14:00 and 16:00 on 22 May. In comparison with the flats on the opposite side of corridor A, the flats adjacent to flat 812 display considerably lower exposure levels (< 0.2%) and shorter risk periods (Figure 5e). Aerosols from the corridor enter flats 814, 816, 818 and 820 only when the wind shifts from an easterly direction to other directions. As for flat 810, an increase in aerosol concentration is observed for a brief period from 20:00 to 24:00 on 22 May and from 21:00 to 24:00 on 24 May; this indicates a substantially lower possibility of horizontal airborne transmission for residents of this flat than for the residents of flats 813 and 817. Among the 20 corridor zones (i.e., 46–55 and 58–67), zones 51–55 present the highest bio-aerosol concentrations, excluding flat 812. These are followed by zones 46–50 and 58–67, with zones 58–67 recording bio-aerosol concentrations trend that are relatively similar to those in zones 51–55 (Figure 5b, 5c). Comparison of the predicted inter-zonal airflow patterns from the MIX (Figure 6) and the spatial infection pattern (Figure 1) shows that the likely infection period for the secondary cases in flats 813 and 817 was the period when the wind direction was easterly and when the index case in flat 812 was at home (Figure 6c). Furthermore, the hours in which the prevailing wind direction was northerly would have been deemed the safest period, in which the least number of residents in their flats (flat 802) were exposed to infection risk, and the region of exposure was almost confined to the corridor (Figure 6a). 16 370 (a) (b) (c) (d) (e) Figure 6. Illustration of net airflow through the main internal openings and the normalised bio-aerosol concentrations (in %) at times characterised by different wind directions. Plots show the patterns of net airflow on the eighth floor of Luk Chuen House 17 380 (under Scenario 1 [All doors closed].) with (a) a northerly wind at 11:00 on 24 May, (b) a north-easterly wind at 19:30 on 22 May, (c) an easterly wind at 15:00 on 22 May, (d) a southeasterly wind at 17:30 on 23 May, and (e) a south-westerly wind at 24:00 on 24 May in the year of 2020. The degree of transparency of the red square in each zone is inversely (but not proportionally) related to its normalised bio-aerosol concentration [refer to the colour scale in (a)]. The direction of net airflow is indicated by the black arrows, with a solid arrow indicating unidirectional airflow and a dashed arrow indicating bidirectional airflow. 385 Figure 7. Normalised tracer gas concentrations (%) in flats 810 and 812 and the corridor zones near the meshed window of flats 808 and 817 monitored at 17:47 on 17 June 2020. 375 390 395 400 To verify the MIX predictions, we compared the model data with those in a related set of tracer gas experiment in Wang et al. (2022); see Figure 7. At a background concentration of 0.017 ppm, the normalised tracer gas concentrations in the corridor 9, in flat 810, and in the corridor 4 were 0.0023%, 0.13%, and 0.50%, respectively; the relative magnitude of the relationship agrees with that modelled by the MIX program (Figure 6e). The concentration used for normalisation was that in the bedroom of flat 812 at the same time (6.196 ppm). Moreover, between 17:47 and 18:08 when the corridor door was open, the tracer gas concentration increased in flat 810 and declined in the two monitored corridor zones (Wang et al., 2022); this was in agreement with the result of the MIX program (Figure 5b, 5e). It is observed that the corridor zones with high concentrations are on the windward side of the source room 812 under a southwesterly wind (Figure 6e). The bio-aerosol concentration distribution at 24:00 on 24 May, was preceded with a time period when the wind direction is easterly (22:00 - 23:00 on 24 May). During 22:00-23:00 on 24 May, the bio-aerosols accumulate in the corridor zones 51-55. Subsequently at 23:00-24:00 on 24 May, the bioaerosol concentrations in corridor zones 51-55 decrease. The airflow direction between flat 812 and corridor 51 during 23:00-24:00 on 24 May was bi-directional, which means virusladen bio-aerosols still penetrated into corridor zones. 18 3.2 Spatial characteristics of exposure risk and the effects of door-opening 405 To elucidate the spatial characteristics of exposure risk, we group all the flats on the eighth floor (excluding flat 812) into eight areas in a clockwise direction and beginning from flat 801. Each area includes five flats or zones, except for area 7, which consists of four flats (flats 820, 818, 816 and 814) (Table S3 and Figure S5). Accordingly, 20 zones in the two corridors are also grouped into four areas. 410 415 For the three scenarios, we calculate the average normalised bio-aerosol concentration in the 12 areas during the index case’s home hours (Figure 8a) and during the suspected infection hours for the secondary cases in flats 813 and 817 (Figure 8b). In all three scenarios, the average concentration for the flats in area 2 is distinctly higher than that for flats in all other areas; this may explain why the residents of the flats in area 2 became infected while the others did not. Interestingly, the corridor zones show higher concentrations than the flats in area 2. We expect that the residents spent most of their time in their own flats. All other flat areas did not register any infections. (a) (b) 19 (c) 420 425 (d) Figure 8. Comparison of the average normalised bio-aerosol concentrations between the 12 areas, and between the flats in area 1 and area 2, under the three scenarios. (Scenario 1 [All doors closed], Scenario 2 [813 and 817 ajar] and Scenario 3 [803 and 807 ajar]). Plots show the average normalised bio-aerosol concentrations for (a) the 12 areas during index case’s home hours, (b) the 12 areas during suspected infection period for the secondary cases in flats 813 and 817, (c) the flats in area 1 and area 2 during the index case’s home hours, and (d) the flats in area 1 and area 2 during the suspected infection period for the secondary cases in flats 813 and 817. Note that values on the y-axis are presented on a log scale and the blue, brown, and green dashed lines notifies the corresponding normalized concentration in flats of area 2 ([a] and [b]) or flat 813 ([c] and [d]) in Scenario 1, 2, 3, respectively. 430 435 We also obtain the average bio-aerosol concentrations for flats directly opposite to flat 812 under the three scenarios (Figures 8c, 8d). This allows for examining the impact of dooropening on a resident’s exposure to aerosols. In contrast to Scenarios 1 and 3, in which the bio-aerosol concentrations of flats in area 2 decrease from flat 811 to 819, under Scenario 2, the average bio-aerosol concentrations in flats 813 and 817 exceed those in flats 811, 815 and 819 (Figures 8c, 8d), due to an increased effective leakage area and higher ventilation rate. The air change rates in flats 813 and 817 increase by 0.9-13 (mean: 3.53) in scenario 2 [813 and 817 ajar] compared to scenario 1 [All doors closed], which brings in more virus-laden 20 440 445 aerosols (Figure 9). We speculate that the residents of flats 813 and 817 may have opened the doors or windows that connected their flats to the corridor and were thus exposed to higher exposure risks than the residents of flats 811, 815 and 819. Nonetheless, this remains a speculation because we have no information on door-opening behaviour and the exact virus release period in this study. In Scenario 3, when the residents of flats 803 and 807 open the doors connecting their flats to the corridor for 20% (effective leakage area, ELA), the air change rates in flats 803 and 807 increase by 0.9-7.6 (mean: 3.12) h-1 compared to scenario 1 [All doors closed] and the aerosol concentrations in flats 803, 807 and 809 increase but remain substantially lower than those in flats in area 2, which also increase due to dooropening of flats 803 and 807. 450 Figure 9. Estimated hourly ventilation rates (ACH) for flats 803, 807, 813, and 817 in 3 scenarios, derived from the MIX program. 4. Discussion 455 460 465 4.1 Horizontal transmission is likely caused by cross-corridor airflow and insufficient ventilation The secondary infections in flats 813 and 817 were probably caused by the cross-corridor airflow from flat 812, where the index case lived. Using the hourly wind data, we inferred the most likely periods of exposure for residents in Luk Chuen House. The local health authorities could not identify any possible encounters between the secondary infected individuals and the index case. Although residents may have shared the same lifts or walked past each other in the corridor, the index case’s work activity patterns (i.e., working the nightshift) would have reduced likelihood of such events happening. Except for the resident of flat 810, no other residents of flats that were on the same side of the corridor as flat 812 experienced secondary infections. This pattern of spatial infection agrees with our modelled pattern of airflow (Figure 6). 21 470 475 480 485 490 495 500 505 It may be useful to reiterate the roles of two connected drainage stacks in facilitating the vertical spread of the virus (Wang et al., 2022). There are 20 pairs of connected drainage stacks in the building: for example, the pair of stacks for -10 and -12 flats discharge wastewater from all -10 and -12 flats. The viral aerosols in the vertical stack may have originated from the washbasin or toilet of flat 812, before it was transported by the stack pair, and leaked into flats 810, 710 and 1112. No other flats that were connected by other stackpairs were contaminated, except for flats 813 and 817, which form a horizontal spatial cluster with flat 812. The results of our simulation suggest that flats in area 2 that were located downstream of flat 812 when the wind direction was easterly during 22–24 May had the highest exposure risk, while all other flats had a low exposure risk (Figure 8a, 8b). This aligns with the observed spatial pattern of infections. Two flats in area 2 had secondary infections, while all flats in the other areas did not (except for flat 810, which belonged to the vertical cluster was the result of contamination via the shared connected drainage stacks; Wang et al., 2022). We cannot rule out the possibility that the two infected neighbours once physically encountered the index case or touched common surfaces that were contaminated. The two residents who contracted COVID-19 in flats 813 and 817 were 68 and 78 years old, respectively. They might have spent more time at home; this would result in a longer exposure duration and a higher probability of overlapping with periods during which the index case was in her flat. Our results showed that door-opening may have also played a role in the cross-corridor transmission of SARS-CoV-2. If they opened their corridor doors by 20% (effective leakage area, ELA), the residents of flats 813 and 817 would have had an exposure risk that was greater than four times that of the residents of other flats in area 2 (Figure 8c, 8d). The index case was the only virus-laden aerosol source in this outbreak, and it was assumed that she remained in her bedroom throughout the entire period for which she was home. Our results suggest that when the wind direction was south-easterly (e.g., 18:00–19:00 on 23 May) or south-westerly (23:00–24:00 on 24 May), the airflow through the door and window of flat 812 would have been a unidirectional inflow. This would have been the desirable situation for the flats opposite to flat 812; the aerosol concentration in these flats would also have diminished rapidly during the hours for which the prevailing wind direction was southeasterly or south-westerly (Figure 5d). Nevertheless, this would have been an undesirable situation for the flats adjacent to flat 812 because the residual aerosols in the corridor would have been transmitted into these flats during these hours (Figure 5e). The coexistence of vertical and horizontal clusters as observed in Luk Chuen House is not new. Both a vertical and horizontal spread was documented during the 2003 Amoy Garden 22 SARS-CoV outbreak (Yu et al., 2004), which involved a transmission across the corridors on each floor in Block E of the estate (Li et al., 2005). 510 515 520 525 530 535 540 Besides the studies reviewed in the Introduction, several studies have reported outbreaks of respiratory infections that were transmitted between rooms. For instance, Hutton et al. (1990) described a tuberculosis (TB) outbreak in a hospital. They found that the droplet nuclei of a draining abscess from a TB patient in a positive pressure room dispersed into the corridor and subsequently to other rooms. The corridor ventilation rate was low, and there was a likelihood that the intake of air from outdoors decreased or ceased due to the cold outdoor weather at the time. Similarly, Gustafson et al. (1982) described a chickenpox outbreak among children in a hospital, with all exposure occurring within one afternoon. The index case remained in strict isolation but was unfortunately placed in a positive pressure room. These studies have demonstrated the combined roles that ventilation rates and room air pressure play in shaping the outbreaks of respiratory infections. Specifically, low ventilation in the index patient’s room led to an accumulation of highly infectious air, while positive pressure – due to either the setup of a room or wind conditions at the time – leaked the infectious air into a poorly ventilated corridor. At the time of revision (early March 2022), the Omicron variants have been rapidly spreading in Hong Kong, particularly in public housing with both horizontal and vertical transmission are suspected to play a role (Wong et al., SCMP, 2022). 4.2 The dual roles of airflow Strictly speaking, still air does not transport aerosols. As airflow transports aerosols that are suspended in the air, it is airflow that transmits or carries airborne particles from one location to another. This is referred to as the ‘transport role’ of the airflow. In addition, the inflow of clean air into a space will dilute any polluted air, and the outflow of polluted air will remove aerosols from the space. This is referred to as the ‘removal role’ of the airflow. For any given airflow, the removal role is associated with the transport role. If polluted air is removed not by transportation of the air to the external environment, but to a neighbouring indoor space, then the effect of removal on a space will bring an effect of pollution on another space. This phenomenon involving the expelling of pollutants from one room and the transporting of the same pollutants to other rooms, for which we refer to as the ‘dual roles of airflow’. Understanding the dual roles of airflow is crucial to infection control. In a corridor-based hotel setting, the settings shown in Figures 10a (left) and 10b (left) will be ideal for quarantine hotels in which such naturally driven flows exist. Figures 10a and 10b (middle and right) show possible pollution situations depending on natural wind directions that drive cross-corridor ventilation; this occurs in Luk Chuen House. Obviously, such flows of infiltration or natural ventilation exist when the building envelope is not airtight. For simplicity, the bi-directional flows across doorways or windows (which exist when there is a 23 545 sufficient buoyancy effect) are not shown in Figure 10. Note that the difference between the middle and right plots of Figures 10a and 10b is the ventilation rate that is represented by the arrow width. (a) (b) 550 555 (c) Figure 10. Schematic diagram of airflow in typical hotel floorplans under different ventilation scenarios. Each small square represents a guest room. The guest room with the index case is always the second room in the top row of each plot. A corridor or guest room that is subsequently polluted is shown in lighter red. The uncontaminated space is shown in white. Green arrows represent the rate and direction of airflow, with arrow widths corresponding to airflow rates. The individual plots show airflow for (a) a floorplan featuring a side corridor, (b) a floorplan featuring a central corridor, and (c) a floorplan where sufficient positive pressure in the corridor prevents the escape of polluted air from the room of the index case into the corridor. 24 560 565 570 575 The dual roles of airflow discussed above may provide at least two suggested approaches for controlling the spread of infections in hotel-like settings. These approaches correspond to two different scenarios – when the building envelope is airtight, and when it is not. When the building envelope is not airtight, sufficient dilution can be used for infection control. In the case of air pollution, the levels of airborne pollutants can be diluted to sufficiently low levels with sufficient ventilation (Figure 10b right). In theory, there is no need to worry about the transport of the aerosols by airflows if sufficient ventilation is available. The challenge lies in determining the exact ventilation rate that suffices for infection control; this question remains to be answered. When the building envelope is airtight, we can use the isolation principle to minimise the negative effect of the transport role, that is, by isolating the contaminated air. This can be achieved by creating a positive pressure in the corridor and making it well ventilated (Figure 10c). In addition to introducing outdoor air to establish positive pressure, the corridor air may be cleaned by filtration and Ultraviolet germicidal irradiation (UVGI) devices. The removal of infectious aerosols differs from the removal of gaseous air pollutants, as aerosols can also be removed by filtration and disinfection such as UVGI. Therefore, each individual guest room may also be equipped with filtration and UVGI. 4.3 Ventilation in quarantine hotels 580 585 590 595 The infection scenario discussed in the present study should be similar to instances of guestto-worker or between-guest infections documented in quarantine hotels that have floorplans similar to that of Luk Chuen House. The inter-zonal flows driven by winds, buoyancy and mechanical ventilation into each room (e.g., the bathroom fans in a hotel room) have the potential to spread viral particles. More infections can occur when the corridor is poorly ventilated (e.g., Figure 10a and 10b, middle) than when it is sufficiently ventilated (Figure 10a and 10b, right). If both the guest rooms and the corridor are well ventilated, the infection risk is low. If the corridor has a sufficiently strong positive pressure, then in theory no infection should occur (i.e., assuming airborne transmission is the predominant transmission route) (Figure 10c). The fact that some items such as quilts, pillows and towels can be contaminated suggests that hotel service personnel may be exposed to infection risk due to the resuspension of particles, followed by inhalation. McKendrick and Emond (1976) presented an intriguing study of the cross-infection risks of chickenpox and measles in seven hospitals with different floorplans. These authors recommended floorplans featuring lateral corridors and cautioned against those featuring central corridors. They also proposed that air from a patient’s room should not be discharged into any corridors and found that the movement of patients and door-opening increased cross25 600 605 610 615 620 625 630 infection risk. Anderson et al. (1985) showed that floorplans featuring patient rooms with negative pressure (i.e., corridors with positive pressure) helped to reduce infections in hospitals. Another study showed that higher ventilation rates led to fewer common cold infections in crowded student dormitories in China (Sun et al., 2011). Many student dormitories in China as well as those in other countries also feature long central corridors like those found in hotels. These studies in hotels and student dormitories have revealed the design characteristics and human behaviours that can help to minimise cross infections: lateral corridors, airflow from corridors to guest rooms, sufficient ventilation, the prohibition or reduction of human movements outside of the guest room, and the minimisation of periods in which the doors of opposite or adjacent rooms are open. Two users of a quarantine hotel in the UK (Warner and Warner, 2021) reported that while they could participate in outdoor exercises, there was a possibility of meeting other quarantined occupants in the corridors. The occupants also opened their doors within seconds of each other to receive delivered food. This was not the experience of one of the present study’s authors (Li), who stayed in quarantine hotels in mainland China and Hong Kong. There were, however, occasions on which test samples were taken, and the doors of the opposite rooms were open within the same period (as the medical team had multiple staff taking samples from different rooms simultaneously). On these occasions, the occupants had to remove their masks for oral or nasal sampling. There are several methods for minimising the leak of air from a room into a corridor. The necessary condition for preventing a leak is that the air pressure of the corridor is always higher than that of the room, that is, a positive pressure is maintained at the corridor. This can be achieved by providing a sufficient supply of air to the corridor. The volume of additional air required may be estimated based on the worst possible pressure profile, which is determined by the air temperatures in the guest rooms and the corridor (e.g., Figure 11) as well as the wind direction. In summer, a cooler temperature may be set in the guest room for human comfort (Figures 11c, 11d). In such a situation, the worst pressure profile is created when positive wind acts on the window (Figure 11c). The wind pushes the neutral plane in the opposite direction, allowing more air to escape into the corridor. A greater effort is needed to push the neutral plane down to the floor level. Under a windier condition, a cooler temperature in the corridor may be the more suitable choice. The worst condition is also shown in Figure 11a when the room temperature is higher than that in the corridor. The required over-supply of airflow to the corridor is also a function of the volume of air that is leaked from the door and the window. A multi-zone airflow program such as MIX (which was used in the present study) can be used for designing room or corridor temperature and opening geometry. 635 26 (a) 640 645 650 655 (b) (c) (d) Figure 11. Estimating the need for an imbalance in airflow or the over-supply of air to the corridor to avoid the escape of air from a guest room to the corridor. The ∆𝑃𝑃 as shown in all plots is the minimum required pressurization in the corridor. The goal of design is to push the neutral plane at the doorway up or down such that no air is leaked from the room to the corridor (desired pressure profile: indicated by the solid red line). The neutral plane is located at mid-height when only a buoyancy force from differences in air density exists (buoyancy-driven pressure profile: indicated by the dotted line). When wind pressure also acts on the window, the neutral plane moves (wind- and buoyancy-driven pressure profile: indicated by the solid black line). Individual plots illustrate the movement of the neutral plane when (a) the guest room is warmer than the corridor, and a positive wind acts on the window, (b) the guest room is warmer than the corridor, and a negative wind acts on the window, (c) the guest room is cooler than the corridor, and a positive wind acts on the window, and (d) the guest room is cooler than the corridor, and a negative wind acts on the window. 4.4 Challenges in improving ventilation in public housing in Hong Kong Although the focus of discussion in this paper centres on quarantine hotels, the studied COVID-19 outbreak occurred in a public housing building in Hong Kong. Thus, our findings also provide insights for improving ventilation in public housing in Hong Kong. Furthermore, many older buildings intended for public housing in Hong Kong were not originally designed to accommodate air conditioning (as this technology was not available at the time). In 2020, 45% of the population in Hong Kong live in public permanent housing (Slotta, 2021). According to Lu et al. (2021), 69.4% of COVID-19 patients in Hong Kong were residents of public housing. 27 660 665 670 675 680 685 690 695 Luk Chuen House is an elongated rectangular block with a long corridor that connects the flats on both sides. Each of these flats has a floor area and single-entry design that is similar to that of a single hotel room. Within each flat, there are no built-in partitions to separate the living room and the bedrooms; some families may make their own basic, flexible partitions, or install solid partitions. Each flat includes a narrow balcony where both the kitchen and a naturally ventilated toilet are located. The balcony is semi-open to the ambient air, with a door that connects to the main flat area, that is, living room or bedroom. An upper meshed window that functions as an air-vent, along with a door, connects the main flat to the corridor. These meshed windows are usually sealed up to facilitate storage or privacy. The layout of Luk Chuen House is the same as that of Wing Shui House, which was studied by Niu and Tung (2007). Their study measured the ventilation rate of the building and revealed the roles of both transport and removal of the corridor air at both ends of the corridor, which were well ventilated, as well as at the middle of the corridor, where the lift lobby was located. A similar effect of ventilation in a cross-corridor floorplan ventilation was documented by Mu et al. (2016, 2017). The positive effect of the balcony on the ventilation and contaminants dispersion in a naturally ventilated building was studied by Cui et al (2014). It was unfortunate that the flat of the index case at Luk Chuen House and the two flats of the secondary infected cases were located at the middle of the exact section of the corridor where ventilation was at its worst when the wind blew in an easterly direction. Ventilation rates in an almost identical floor plan on the thirteenth floor in a 16-storey public housing building were measured by Wu et al. (2016, 2019). Interestingly, the same authors (Wu et al., 2016) previously predicted that flats downstream from a source of contaminants would experience a significant infection risk. Our findings for the horizontal spatial cluster in the COVID-19 outbreak at Luk Chuen House confirm these predictions. There is an urgent need to improve ventilation in public housing. At Luk Chuen House, we observed that many residents had placed objects that obstructed the airflow from the meshed windows of their flats to the corridor. Such old public housing buildings were designed and constructed during the period when air conditioners were not yet widely used. Air conditioners had since been installed in almost all flats, and to save energy when their air conditioning was switched on, the residents would often shut the windows and any openings of their flats. However, we found that the openings in the walls at the two ends of the corridor and at the central lift lobby area provided insufficient ventilation for the corridor. When airflow was obstructed at the doors and meshed windows of the flats upwind, the ventilation in the corridor would become poor, and the flats downwind would not be well ventilated. It would be very difficult to ensure adequate ventilation in such situations. A well-ventilated corridor is essential for minimising the cross-corridor transmission of infectious viruses and other airborne pollutants. One possible solution is to seal the corridor (i.e., by sealing the openings at the two ends of the corridor and the centre lift lobby with an operable airtight window), and to install a positive pressure ventilation system for the corridor. 28 700 4.5 Limitations 705 We recognise several major limitations of the study. First, the field experiment was conducted in a very short period; the research team was only provided with access to the premises for two days (Wang et al., 2022). The team had no access to the flats opposite to flat 812 on the eighth floor. It was not possible to take direct measurements of the tracer gas concentrations in flats 813 and 817, although we placed a sensor near the meshed window of flat 813, which was accessible from the corridor. In addition, the prevailing wind directions during the field experiment differed from the easterly wind direction that prevailed during the suspected period of exposure. Thus, the estimates of exposure risk were based on a multi-zone airflow model. 710 715 720 725 730 735 Second, our study lacked information on the behaviour of the residents, including the exact times at which they stayed in their flats during the suspected exposure period, whether and how often they opened the doors and windows of their flats, and their use of air-conditioning. The latter two factors would have influenced the air temperatures in their flats, which may have affected the buoyancy effect and inter-zonal flows. Given that the affected flats on the eighth floor of the building are in relatively open surroundings, their airflow is likely predominantly driven by wind, and the buoyancy effect is not expected to be significant. Third, the CFD-predicted pressure coefficients are for a sealed building, which may lead to errors in the predicted airflow rates by MIX (Gautam et al., 2019). However, an early study showed that when the wall porosity is less than 10-20% of the wall area (Seifert et al., 2006), such as sealed building approach can provide reasonably accurate results.. 5. Conclusion In the COVID-19 outbreak at Luk Chuen House, cross-corridor airflows may explain the spread of SARS-CoV-2-laden aerosols from the flat where the index case lived to two flats on the other side of the central corridor. During the suspected period of exposure – that is, the day of and two days following the symptoms onset in the index case – a prevailing easterly wind resulted in a high concentration of viral aerosols in flats 811, 813, 815, 817 and 819. Flats 814, 816, 818 and 820, which were on the same side of the corridor, had aerosol concentrations that were approximately 2 orders of magnitude lower than that in the flat of the index case. Other flats on the same floor that were distantly located from the flat of the index case (e.g., 822, 824, 826, and 828) were almost risk-free. When the doors of flats 813 and 817 that connected to the corridor were open, the residents of these flats experienced the highest exposure risk, which was higher than that experienced by the residents of flats 811, 815 and 819. The results from our computer modelling align with the observed patterns of spatial infection in the horizontal cluster of the outbreak. It will be a challenge to improve ventilation conditions in similar public housing buildings. However, portable air cleaners 29 with high efficiency particulate air (HEPA) filtration or ultraviolet germicidal irradiation (UVGI) may be applied. 740 745 750 755 Following a detailed consideration of the dual roles of airflows in long-corridor buildings – that is, the removal and transport roles – we identify two distinct strategies for controlling the spread of airborne infections in such buildings. The first is to provide sufficient ventilation in all areas, such that the transport of infectious aerosols is no longer an issue. This may be feasible with well-designed and naturally ventilated buildings. The second strategy is to control the transport of infectious aerosols, such as by controlling the direction of airflow between rooms. The cross-corridor transmission investigated in the present study appears to reflect the outbreak patterns of COVID-19 infections in quarantine hotels worldwide. Most of these quarantine hotels were not originally designed for quarantine purposes. We therefore recommend that, if available, single-sided corridor hotels with the right prevailing wind directions are more suited to the purposes of quarantining cases of airborne infections. In hotels with central corridors, the maintenance of positive pressure and sufficient ventilation in the corridor must be ensured to help reduce infections. Acknowledgements This work was supported by Research Grants Council (RGC) of Hong Kong, China [grant number 1720 3321]. 760 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. Author contributions 765 Pan Cheng: Methodology, Formal analysis, Investigation, Data Curation, Writing - Original Draft, Visualisation; Wenzhao Chen, Shenglan Xiao, Fan Xue, Qun Wang, and Pak Wai Chan: Methodology, Writing - Review & Editing; Jianlei Niu and Zhang Lin: Data Curation, Writing - Review & Editing; Yuguo Li: Conceptualisation, Methodology, Writing Original Draft, Review & Editing, Funding acquisition 770 References • Anderson, J. D., Bonner, M., Scheifele, D. W., & Schneider, B. C. (1985). Lack of nosocomial spread of varicella in a pediatric hospital with negative pressure ventilated patient rooms. Infection Control, 6(03), 120–121. 30 775 • Bahadori-Jahromi, A., Rotimi, A., Mylona, A., Godfrey, P., & Cook, D. (2017). Impact of window films on the overall energy consumption of existing UK hotel buildings. Sustainability, 9(5), 731–753. • British Broadcasting Corporation (BBC). (2021). Covid: ‘Urgent’ call to shut Reading quarantine hotel after outbreak. https://www.bbc.com/news/uk-england-berkshire57417585 (Accessed 20 December 2021). • Centers for Disease Control and Prevention (CDC). (2003). Update: Outbreak of Severe Acute Respiratory Syndrome — Worldwide, 2003. Morbidity and Mortality Weekly Report, 52(12), 241-248. • Cui, D. J., Mak, C. M., & Niu, J. L. (2014). Effect of balconies and upper–lower vents on ventilation and indoor air quality in a wind-induced, naturally ventilated building. Building Services Engineering Research and Technology, 35(4), 393–407. • CityNews. (2021). 2 quarantine hotels in the GTA hit with COVID-19 outbreaks. https://toronto.citynews.ca/2021/05/11/quarantine-hotel-covid-19-outbreak/ (Accessed 31 October 2021). • Confluence. (2021). https://confluence.cornell.edu/display/SIMULATION/FLUENT++Flow+over+an+Airfoil-+Step+6 (Accessed on 8 January 2022). • Cable News Network (CNN). (2021). China has built a 5,000-room quarantine center for overseas arrivals. It could be the first of many. https://edition.cnn.com/2021/09/29/china/guangzhou-covid-quarantine-center-mic-intlhnk/index.html (Accessed 2 November 2021). • Centers for Disease Control and Prevention (CDC). (2021). Science brief: SARS-CoV-2 and surface (fomite) transmission for indoor community environments. https://www.cdc.gov/coronavirus/2019-ncov/more/science-and-research/surfacetransmission.html (Accessed 19 September 2021). • Fernando, G., & McPhee, S. (2020). Coronavirus: Victoria’s new lockdown could cost $1 billion a week, Treasurer says. https://www.news.com.au/finance/economy/australianeconomy/coronavirus-victorias-new-lockdown-could-cost-1-billion-a-week-treasurersays/news-story/9d26f9986565aee82cfeb9623c0b0ee6 (Accessed 31 October 2021). • Gustafson, T. L., Lavely, G. B., Brawner Jr, E. R., Hutcheson Jr, R. H., Wright, P. F. & Schaffner, W. (1982). An outbreak of airborne nosocomial varicella. Pediatrics, 70(4), 550–556. • Gautam, K. R., Rong, L., Zhang, G., & Abkar, M. (2019). Comparison of analysis methods for wind-driven cross ventilation through large openings. Building and Environment, 154, 375-388. • Grout, L., Katar, A., Ait Ouakrim, D., Summers, J. A., Kvalsvig, A., Baker, M. G., Blakely, T., & Wilson, N. (2021). Failures of quarantine systems for preventing COVID‐ 780 785 790 795 800 805 810 19 outbreaks in Australia and New Zealand. Medical Journal of Australia, 215(7), 320–324. 31 • Hoefer, A., Pampaka, D., Wagner, E. R., Herrera, A. A., Alonso, E. G. R., López-Perea, N., Portero, R. C., Herrera-León, L., Herrera-León, S., & Gallo, D. N. (2020). Management of a COVID-19 outbreak in a hotel in Tenerife, Spain. International Journal of Infectious Diseases, 96, 384–386. • Hutton, M. D., Stead, W. W., Cauthen, G. M., Bloch, A. B., & Ewing, W. M. (1990). Nosocomial transmission of tuberculosis associated with a draining abscess. Journal of Infectious Diseases, 161(2), 286–295. • Jiang, F. C., Jiang, L., Wang, Z. G., Meng, Z. H., Shao, S. F., Anderson, B. D., & Ma, M. J. (2020). Detection of severe acute respiratory syndrome coronavirus 2 RNA on surfaces in quarantine rooms. Emerging Infectious Diseases, 26, 2162–2164. • Li, Y., Delsante, A., & Symons, J. (2000). Prediction of natural ventilation in buildings with large openings. Building and Environment, 35(3), 191–206. • Li, Y., Duan, S., Yu, I. T. S., & Wong, T. W. (2005). Multi-zone modeling of probable SARS virus transmission by airflow between flats in Block E, Amoy Gardens. Indoor Air, 15(2), 96–111. • Li, Y., Qian, H., Hang, J., Chen, X., Cheng, P., Ling. H., … & Kang, M. (2021). Probable airborne transmission of SARS-CoV-2 in a poorly ventilated restaurant. Build and Environment, 196, 107788. • Li, X., Chen, H., Lu, L., Chen, L. L., Chan, B. P. C., Wong, S. C., Cheng, V. C.C., Yuen, K. Y., Chan, K. H., & To, K. K. W. (2021). High compliance to infection control measures prevented guest-to-staff transmission in COVID-19 quarantine hotels. Journal of Infection, S0163-4453(21), 00533-8. • Leong, L. E., Soubrier, J., Turra, M., Denehy, E., Walters, L., Kassahn, K., ... & Lim, C. K. (2021). Whole-genome sequencing of SARS-CoV-2 from quarantine hotel outbreak. Emerging Infectious Diseases, 27(8), 2219–2221. • Lu, H., Li, J., & Ni, X. (2021). Promoting resilience in the age of covid-19 pandemic: a new era of strategic foresight. Aging and Disease, 12(1), 1–2. • McKendrick, G. D. W., & Emond, R. T. D. (1976). Investigation of cross-infection in isolation wards of different design. Epidemiology & Infection, 76(1), 23–31. • Mu, D., Gao, N., & Zhu, T. (2016). Wind tunnel tests of inter-flat pollutant transmission characteristics in a rectangular multi-storey residential building, part A: Effect of wind direction. Building and Environment, 108, 159–170. • Mu, D., Shu, C., Gao, N., & Zhu, T. (2017). Wind tunnel tests of inter-flat pollutant transmission characteristics in a rectangular multi-storey residential building, part B: Effect of source location. Building and Environment, 114, 281–292. • Miller S. L., Nazaroff, W. W., Jimenez, J. L., Boerstra, A., Buonanno, G., Dancer, S. J., & Noakes, C. (2021). Transmission of SARS-CoV-2 by inhalation of respiratory aerosol in the Skagit Valley Chorale superspreading event. Indoor Air, 31(2), 314–323. 815 820 825 830 835 840 845 850 32 • Niu, J. & Tung, T. C. (2007). On-site quantification of re-entry ratio of ventilation exhausts in multi-family residential buildings and implications. Indoor Air, 18(1), 12–26. • Orme, M., Liddament, M. & Wilson, A. (1998). Numerical data for air infiltration and natural ventilation calculations. International Energy Agency. Paris. • Ou, C., Hu, S., Luo, K., Yang, H., Hang, J., Cheng, P., … & Li, Y. (2022). Insufficient ventilation led to a probable long-range airborne transmission of SARS-CoV-2 on two buses. Building and Environment, 207, 108414. • Public Health Expert (PHE). (2021). Australia’s Quarantine Systems Failures: Lessons for NZ. https://sciblogs.co.nz/public-health-expert/2021/05/12/australias-quarantine-systemsfailures-lessons-for-nz/ (Accessed 7 January 2022). • Radio Television Hong Kong (RTHK). (2021). Experts suspect cross-infection at quarantine hotel. https://news.rthk.hk/rthk/en/component/k2/1606261-20210817.htm (Accessed 31 October 2021). • Sun, Y., Wang, Z., Zhang, Y., & Sundell, J. (2011). In China, students in crowded dormitories with a low ventilation rate have more common colds: evidence for airborne transmission. PloS One, 6(11), e27140. • Slotta, D. (2021). Population distribution in Hong Kong 2020 by type of housing. https://www.statista.com/statistics/629802/hong-kong-population-by-type-of-housing/ (Assessed 13 January 2021). • Thatcher, T. L., Lai, A. C., Moreno-Jackson, R., Sextro, R. G., & Nazaroff, W. W. (2002). Effects of room furnishings and air speed on particle deposition rates indoors. Atmospheric Environment, 36(11), 1811-1819. • Wu, Y., Tung, T. C., & Niu, J. L. (2016). On-site measurement of tracer gas transmission between horizontal adjacent flats in residential building and cross-infection risk assessment. Building and Environment, 99, 13–21. • Wu, Y., Tung, T. C., & Niu, J. L. (2019). Experimental analysis of driving forces and impact factors of horizontal inter-unit airborne dispersion in a residential building. Building and Environment, 151, 88–96. • World Health Organization (WHO). (2021). Coronavirus disease (COVID-19): How is it transmitted? Updated 23 December 2021. https://www.who.int/news-room/questions-andanswers/item/coronavirus-disease-covid-19-how-is-it-transmitted (Accessed on 10 March 2022). • Warner, J. O., & Warner, J. A. (2021). The UK’s hotel quarantine system is not fit for purpose. The bmj opinion. https://blogs.bmj.com/bmj/2021/06/04/the-uks-hotelquarantine-system-is-not-fit-for-purpose/ (Accessed on 2 November 2021). • Wong, S. C., Chen, H., Lung, D. C., Ho, P. L., Yuen, K. Y. & Cheng, V. C. C. (2021). To prevent SARS‐CoV‐2 transmission in designated quarantine hotel for travellers: Is the 855 860 865 870 875 880 885 ventilation system a concern? Indoor Air, 31(5), 1295–1297. 890 • Wang, Q., Lin, Z., Niu, J., Choi, G. K. Y., Fung, J. C., Lau, A. K., ... & Li, Y. (2022). 33 Spread of SARS-CoV-2 aerosols via two connected drainage stacks in a high-rise housing outbreak of COVID-19. Journal of Hazardous Materials, 128475. • Wong N, Cheng : and Cheung, E. 2022. Coronavirus: what kind of lockdown can Hong Kong handle ... or is it even ready at all? South China Morning Post. 2 March, 2022. https://www.scmp.com/news/hong-kong/healthenvironment/article/3168870/coronavirus-what-kind-lockdown-can-hong-kong. Last accessed 12 March 2022. • Xiao, S., Li, Y., Sung, M., Wei, J., & Yang, Z. (2018). A study of the probable transmission routes of MERS‐CoV during the first hospital outbreak in the Republic of 895 Korea. Indoor Air, 28(1), 51–63. 900 • Yu, I. T., Li, Y., Wong, T. W., Tam, W., Chan, A. T., Lee, J. H., ... & Ho, T. (2004). Evidence of airborne transmission of the severe acute respiratory syndrome virus. New England Journal of Medicine, 350(17), 1731-1739. 34 905 Supplementary Materials Probable cross-corridor transmission of SARS-CoV-2 due to cross airflows and its control 910 915 920 Pan Cheng, Wenzhao Chen, Shenglan Xiao, Fan Xue, Qun Wang, Pak Wai Chan, Jianlei Niu, Zhang Lin, and Yuguo Li* S1 The index case’s daily routine Her daily shift begins at 20:30 and ends at 04:00 the next day. According to Google Maps, it takes around 50 min for her to reach her workplace from home or to reach home from her workplace (Kerry Logistics Network Limited) by bus or Mass Transit Railway (MTR) (Google, 2021). One hour is deemed to be a sufficient time for the index patient to get to her workplace from home or get home from work. The index case stopped going to work on May 24. Therefore, she was home from 05:00 to 19:30 on 22 and 23 May, 2020 and from 05:00 to 24:00 on 24 May, 2020. We assumed that the index case stayed in her flat (bedroom) during her time at home. 35 925 S2 Wind pressure coefficients obtained by CFD simulation To validate the model, a uniform velocity inlet (7.675 m/s) case with a 0.2-m cube mounted on the floor was used to compare our results to those from wind tunnel experiments by Castro and Robins (1977); see (Figure S1). The validation case was facilitated by a mesh size of 1,290,366 cells, with a density box around the cube for local refinement. In the studied area, the mean absolute error between the CFD simulation results in this study and the experiment results in Castro and Robins (1977) is around 30% with a minimum value near 4%. (a) 930 935 940 (b) Figure S1. Schematic display of a comparison between the results of the CFD simulation and the wind tunnel experiment by Castro and Robins (1977). In this study, unstructured mesh was used. For the grid independence test, three tested cases were performed; each used a different maximum cell size around the buildings (1.8 m, 2 m or and 3 m). Away from the buildings, the maximum cell size was the same for the three tested cases above. The total numbers of cell elements were 5,749,814 (fine), 4,834,083 (median) and 3244707 (coarse), respectively. Note that only the scenario for an easterly wind was tested for grid independence. The other four wind inlet scenarios adopted the median grid setting, as the median grid produced results that were almost identical to those obtained from the fine grid (Figure S2). 36 (a) (b) Figure S2. Results of the grid independence test, with the wind-pressure coefficient along vertical lines AB (windward) and EF (leeward) on the facade of Luk Chuen House when the incoming wind is easterly (a) and the layout of the employed grid (b). 37 945 950 955 S3 Geometry of the eighth floor of Luk Chuen House and the multi-zone model On the eighth floor of Luk Chuen House, there are 40 rooms (zones 01–44), three stairwells (zones 45, 57 and 68), two corridors (zones 46–55 and 58–67) and one lift lobby (zone 56). The entire floor is divided into 68 zones in MIX program. Flats 810 and 812 are separated into three zones (zones 35–37 and 38–40) each, and the two corridors are split into 10 zones each (Figure S3). In each flat, there is one exhaust fan and one window (Window 1 in Figure S3) on the external wall [Figure S4a], as well as one meshed window and one door (Window 3 and Door 3 in Figure S3) connected to the corridor [Figure S4d]. In flats 810 and 812, there is a door and a void space (‘Door 1/Void space’ shown in Figure S3) between the kitchen (zones 35 and 39) and the toilet (zones 36 and 38) [Figure S4b]; a door and a window (‘Door 2’ and ‘Window 2’ in Figure S3) connecting the kitchen and the bedroom (zone 37 and 40) [Figure S4c]. The dimensions of all of the 68 zones are listed in Table S1. There are two openings on the external walls at the lift lobby and the two stairwells (zones 56, 45 and 68). The stairwell (zone 57) connected to the lift lobby has one external opening. The dimensions and effective leakage area (ELA) coefficients for all of the openings are listed in Table S2. 960 965 There are three entrances to the building: one entrance is located on each end of the building (each with stairwell access), while another entrance at the middle of the building serves as the main entry portral (with both lift and stairwell access). A thorough investigation by the local health authority did not reveal any evidence to suggest that transmission occurred via the sharing of lifts (the building has separate lifts for odd and even floors) or other common areas. After the vertical cluster was identified, all other residents living in -10 and -12 flats were evacuated and placed in a 2-week quarantine beginning from 4 June, 2020. 970 Figure S3. Division of the eighth floor into 68 zones for multi-zone airflow modelling. The elevators, refuse room and switch room (filled with oblique lines) are excluded from the 38 975 computational domain. Note that the void area of the toilet is above the door, and that Window 1 on the facade is indicated by a yellow bar if it is composed of irons bars only, and a blue bar if glass windows exist but are kept open. (a) (b) (c) (d) 39 Figure S4. Schematic display of eight openings for a flat. Table S1 Dimensions and air temperature for all of the 68 zones. 980 Zone No. Length, m Flats 1, 20 5.613 2–9, 12–19 5.613 10, 11 5.613 21 7.186 22–29 7.186 30 7.186 31 8.763 32–34, 41–43 8.763 35, 39 3.0 36, 38 1.42 37, 40 7.463 44 8.763 Stair cases 45, 68 8.315 57 5.613 Corridor zones 46 4.509 47–54 4.42 Width, m Height, m Temperature, °C 55 58 59–66 4.509 4.42 4.623 4.509 4.42 4.623 4.623 4.42 1.3 1.3 4.42 4.623 2.713 2.713 2.713 2.713 2.713 2.713 2.713 2.713 2.713 2.713 2.713 2.713 24.15 24.15 24.15 24.15 24.15 24.15 24.15 24.15 27.15 27.25 24.15 24.15 3.353 3.353 2.713 2.713 26.15 26.10 1.81 1.81 2.713 2.713 4.623 4.623 4.42 1.81 1.829 1.829 2.713 2.713 2.713 67 4.509 1.829 2.713 26.05 25.85, 25.75, 25.65, 25.4, 25.5, 25.65, 25.75, 25.85 26.05 26.05 25.85, 25.75, 25.65, 25.5, 25.4, 25.65, 25.75, 25.85 26.05 Lift lobby 56 12 7.065 2.713 26.15 Table S2 Dimensions and effective leakage area (ELA) coefficients for all of the openings of flats. Structure (Position) Exhaust fan (kitchen) Width (m) 0.4 Height (lower–higher bound) 1.3 (0.913–2.213) 40 ELA coefficient (%) 20 985 Window 1 (kitchen) 2.1 1.3 (0.913–2.213) Door 1 (toilet) Void space (toilet) Door 2 (bedroom) Window 2 (bedroom) Door 3 (flat) Window 3 (flat) 0.65 0.70 0.70 1.30 1.9 (0–1.9) 0.5 (1.9–2.4) 2.0 (0–2.0) 1.0 (1.0–2.0) 95/75 (without/with glass) 1.8 (closed) 100 (fully open) 1.8 (closed) 20 (half open) 0.75 2.0 2.0 (0–2.0) 0.5 (2.0–2.5) 1.8/20 (closed/ajar) 0.4 (closed) The virus-laden bio-aerosols are simulated as passive tracers (Li et al., 2005) and have no influence on the airflow, as small particles evaporate rapidly after being released. Deposition loss rate coefficients for particles under different furnishing conditions (Thatcher et al., 2002) are adopted here to consider bio-aerosol removal caused by particle deposition inside the flats and in the corridors. Table S3 Areas; see also Figure S5. Area no. 1 2 3 4 5 6 7 8 9 10 11 12 Flats or corridor zones included 801, 803, 805, 807, 809 811, 813, 815, 817, 819 821, 823, 825, 827, 829 831, 833, 835, 837, 839 832, 834, 836, 838, 840 822, 824, 826, 828, 830 814, 816, 818, 820 802, 804, 806, 808, 810 1, 2, 3, 4, 5 6, 7, 8, 9, 10 11, 12, 13, 14, 15 16, 17, 18, 19, 20 990 Table S4 Deposition loss rate coefficient (h-1) Particle diameter range (𝝁𝝁𝒎𝒎) 0 - 0.6 0.6 - 0.73 0.73 - 0.905 0.905 - 1.12 1.12 - 1.39 Residence flats Common areas 0.23 0.24 0.27 0.33 0.47 0.09 0.10 0.11 0.15 0.25 41 1.39 - 1.725 1.725 - 2.14 2.14 - 2.655 2.655 - 3.295 3.295 - 4.04 4.04 - 5.075 5.075 - 6.30 6.30 - 7.84 7.84 - 10 995 0.67 0.93 1.32 1.93 3.39 4.71 5.73 7.78 10.5 0.39 0.61 0.92 1.45 2.54 3.79 4.88 6.48 8.84 Figure S5. Schematic display of 12 areas. The flat in red (812) is where the index case lived. The two flats with secondary infections are indicated by Blue (813 and 817). Light green and yellow are used to distinguish zones. (See also Table S3). 42 1000 S4 Weather condition on the day of tracer gas measurement and measurement results Figure S6. Hourly temperatures, wind speeds and wind directions recorded at Shatin weather station (located proximally to Luk Chuen House) on 16–17 June, 2020, when the tracer gas measurements were conducted. 1005 43 S5 Bio-aerosol concentration normalisation 1010 Table S5 Average bio-aerosol concentration in flat 812 during the index case’s home hours on each day and all days. Scenario No. 1 2 3 05:00–19:30 on 22 May 34.27 #/m3 33.37 #/m3 33.70 #/m3 05:00–19:30 on 23 May 34.16 #/m3 33.83 #/m3 33.82 #/m3 05:00–24:00 on 24 May 35.62 #/m3 35.05 #/m3 35.55 #/m3 All home hours on 22–24 May 34.77 #/m3 34.36 #/m3 34.47 #/m3 From Table S5, it is apparent that the average bio-aerosol concentration in flat 812 is the highest in Scenario 1. Therefore, we use the average bio-aerosol concentration in flat 812 during all of the index case’s home hours to normalise and report the bio-aerosol concentration data for other zones (34.77 #/m3). 1015 44 S6 Inflow and outflow rates for flat 812 (a) 1020 (b) Figure S7. Hourly inflow and outflow rates for flat 812, derived from the MIX program. Plots show the (a) inflow rate and (b) outflow rate. EF, W1, D3 and W3 indicate the exhaust fan, Window 1, Door 3 and Window 3, respectively, in Figure S3. The estimated average flow rate during the 72-h computation period by our multi-zone program (MIX) is 4.28 ACH, which agrees with the value measured by Tung et al. (2021) in a similar apartment under strong wind conditions (i.e., 4.18 ACH). 1025 References • Castro, I. P., & Robins, A. G. (1977). The flow around a surface-mounted cube in uniform and turbulent streams. Journal of fluid Mechanics, 79(2), 307–335. 45 • Google. (2021). https://www.google.com/search?q=From+Luk+Chuen+House+to+Kerry+Logistics+Netw ork+Limited&oq=from+Luk+&aqs=chrome.0.69i59j69i57j69i59.6902j0j15&sourceid=c hrome&ie=UTF-8 (Accessed on 8 January 2022). • Hsu, S. A., Meindl, E. A., & Gilhousen, D. B. (1994). Determining the power-law windprofile exponent under near-neutral stability conditions at sea. Journal of Applied Meteorology and Climatology, 33(6), 757–765. • Li, Y., Duan, S., Yu, I. T. S., & Wong, T. W. (2005). Multi-zone modeling of probable SARS virus transmission by airflow between flats in Block E, Amoy Gardens. Indoor Air, 15(2), 96–111. • Thatcher, T. L., Lai, A. C., Moreno-Jackson, R., Sextro, R. G., & Nazaroff, W. W. (2002). Effects of room furnishings and air speed on particle deposition rates indoors. Atmospheric Environment, 36(11), 1811-1819. • Tung, C. W., Mak, C. M., Niu, J. L., Hung, K., Wu, Y., Tung, N., & Wong, H. M. (2021). Enlightenment of re-entry airflow: The path of the airflow and the airborne pollutants transmission in buildings. Building and Environment, 195, 107760. • Wang, Q., Lin, Z., Niu, J., Choi, G. K. Y., Fung, J. C., Lau, A. K., ... & Li, Y. (2022). Spread of SARS-CoV-2 aerosols via two connected drainage stacks in a high-rise housing outbreak of COVID-19. Journal of hazardous materials, 128475. 1030 1035 1040 1045 46