Tracing contacts to evaluate the transmission of COVID-19 from highly exposed individuals in public transportation ================================================================================================================== * Caio Ponte * Humberto A. Carmona * Erneson A. Oliveira * Carlos Caminha * Antonio S. Lima Neto * José S. Andrade, Jr. * Vasco Furtado ## Abstract We investigate, through a data-driven contact tracing model, the transmission of COVID-19 inside buses during distinct phases of the pandemic in a large Brazilian city. From this microscopic approach, we recover the networks of close contacts within consecutive time windows. A longitudinal comparison is then performed by upscaling the traced contacts with the transmission computed from a mean-field compartmental model for the entire city. Our results show that the effective reproduction numbers inside the buses, *Re**bus*, and in the city, *Re**city*, followed a compatible behavior during the first wave of the local outbreak. Moreover, by distinguishing the close contacts of healthcare workers in the buses, we discovered that their transmission, *Re**health*, during the same period, was systematically higher than *Re**bus*. This result reinforces the need for special public transportation policies for highly exposed groups of people. Keywords * COVID-19 * Public Transportation * Contact Tracing * Complex Networks * Compartmental Models ## I. INTRODUCTION Human mobility is crucial to understanding the COVID-19 pandemic since the Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) is disseminated individual to individual via droplet and airborne transmissions [1]. Considering that non-pharmacological interventions, such as social distancing and isolation, still represent fundamental measures to control the COVID-19 outbreak, the nature of SARS-CoV-2 dissemination also unveils the need to understand the role of the space where interaction between people occurs. There is consensus that superspreading events, which are usually investigated through contact tracing models [2–7], are more likely to happen in indoor environment, as substantiated by previous studies of indoor contagion in hospitals [8, 9], restaurants [10], offices [11], and even on cruise ships [12, 13]. However, the relation between the microscopic level of contagion in indoor environments and the macroscopic observables, such as the numbers of cases and deaths at the city scale remains unclear. Public transportation is one of the main forms of commuting, playing an important role in the pace of life in cities [14], specially in epidemics [15–20]. In spite of the fact that some cities have adopted social distancing and sanitary protocols on public transportation to control the COVID-19 outbreak, it is common that buses or subways get crowded at rush hour, mainly in the developing countries. In recent months, some studies have been proposed to establish the safety of public transportation regarding the indoor COVID-19 contagion [21–27]. To the best of our knowledge, however, none of these studies have considered the possibility of a comparative analysis based on a two-fold perspective, namely, the dynamic of people’s movement in a city and the dynamic of the virus dissemination within vehicles of public transport. Here, using data about people’s movement on buses and COVID-19 infection in a large metropolis, we define two data-driven mathematical models based on concepts of complex networks and non-linear dynamics in order to foster the understanding of the role public transportation plays in the COVID-19 pandemic. At the microscopic scale, we define a contact tracing model to estimate the transmission within city buses and, at the macroscopic scale, a compartmental model is employed to estimate the transmission in the entire city. The main contribution of our study is a comparative analysis between these two distinct modeling approaches through the combination of daily epidemiological and mobility data during the first 9 months of the local COVID-19 outbreak, and through different social distancing restriction regimes. One specially relevant aspect of this work is the fact that we are able to trace within the public transportation vehicles (*i*.*e*., indoor environments) two groups of people, one of them with a higher exposure to the virus in comparison to the other. This allows us to shed light on potentially dangerous superspreading events in public transportation. ## II. MODELING APPROACH ### A. Contact tracing model We propose a contact tracing model using two datasets that relate bus validations to COVID-19 confirmed cases during the periods of social isolation, lockdown, and economic reopening in the city of Fortaleza, Ceará, Brazil (see Methods). Our model is a network based on Potentially Infectious Contacts (PICs), in which bus passengers during their infectious period -according to subsequent diagnosis of COVID-19 -have shared the transport for a certain amount of time with other passengers, the latter in their exposed period -also according to subsequent COVID-19 diagnosis. Precisely, the proposed network is composed of vertices *p**i* that represent the passengers diagnosed with COVID-19, and weighted directed edges *c**k* = (*p**i*, *p**j*, *τ**ij*) that represent PICs. For each edge, the direction is assigned from an infectious passenger *p**i* to an exposed passenger *p**j*, and the weight *τ**ij* is defined as the estimated value of the ride time shared by *p**i* and *p**j* on the same bus, as shown in Fig. 1a. We calculate *τ**ij* by superimposing the estimated ride times from *p**i* and *p**j*, considering the different moments of their boarding. Here, the epidemiological profile for COVID-19 transmission is characterized by the dates of the passengers’ Onset of Symptoms (OS). The infectious period corresponds to the days in which a passenger diagnosed with COVID-19 can transmit the virus, initiating 2 days before OS and ending 12 days after OS. The exposed period refers to the time window during which the passenger can get the virus and maintain it latent until the infectious period. In this context, the exposed period begins 14 days before OS and ends 2 days before OS, *i*.*e*., the infectious and the exposed periods have a width of 14 and 12 days, respectively, and they do not overlap [28–30]. Furthermore, if there is more than one PIC related to an exposed passenger *p**j*, we consider solely the edge with the largest value of *τ**ij*. It is important to notice that, by crossing the datasets of bus validations and confirmed cases of COVID-19 in Fortaleza during the period from March to December 2020, we are able to identify 5, 159 pairs of infectious and exposed passengers that rode the same bus on the same day. However, their associated values of *τ**ij* could only be computed for 3, 023 (58.6%), due to missing information in the dataset of bus validations. From these pairs, we obtain that the network of PICs corresponds to a forest composed of 213 trees with a total of 530 vertices (infectious passengers) and 317 edges (PICs). From all vertices found, 97 were identified as healthcare workers (see Methods). The Centers for Disease Control and Prevention (CDC) recommends that any contact tracing strategy for COVID-19 should consider the concept of Close Contacts (CCs) [31], *i*.*e*., anybody who has been for at least 15 minutes within 6 feet (≈ 2 meters) of an infectious person. Since buses are small, enclosed, and they have a great tendency to get crowded at rush hours, we define the CCs in the network of PICs only considering the time condition *τ**ij* *> τ**c*, where the threshold *τ**c* = 15 minutes. Applying this criterion to the network of PICs, we find that the network of CCs is composed of 154 trees with a total of 360 vertices (infectious passengers) and 206 edges (CCs). In this case, 75 vertices were identified as healthcare workers. In order to understand the COVID-19 spreading in public transportation, we define the effective reproduction number for the contact tracing model, *Re**bus*, as the expected number of secondary cases produced by a single (typical) infection. Precisely, it accounts for two contributions in relation to who is spreading the disease: one due to reported infectious individuals, ![Graphic][1], and another due to unreported infectious individuals ![Graphic][2]. Here, we assume that the fraction of newly reported to newly unreported cases generated by a typical reported infectious individual remains invariant during time. This is equivalent to consider the value of ![Graphic][3] proportional to the average number of outdegrees from the vertices in the network of CCs during a given time window, ![Graphic][4], ![Fig. 1.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/06/03/2021.05.31.21257245/F1.medium.gif) [Fig. 1.](http://medrxiv.org/content/early/2021/06/03/2021.05.31.21257245/F1) Fig. 1. Proposed Models for COVID-19 and Spreading Scenarios. **(a)** Potentially Infectious Contacts (PICs). We define a PIC when an infectious passenger *p**i* (in red) and an exposed passenger *p**j* (in yellow) share the same bus. The weight *τ**ij* is the estimated value of the ride time shared by *p**i* and *p**j*. The time lines show the infectious (in red) and the exposed periods (in yellow) of each passenger, where each square represents one day. The time lines are built based on the Onset of Symptoms (OS). Precisely, the infectious period begins 2 days before OS and ends 12 days after OS, while the exposed period begins 14 days before OS and ends 2 days before OS. Other passengers (in gray), even though they have shared the same bus with *p**i*, either were not notified as COVID-19 cases or, however notified, they were not considered as PICs because they were not in their exposed period. **(b)** The SEIIR model. The total population of size *N* provides the susceptible population *S* (in blue). The susceptible individuals become exposed *E* (in yellow) at a time-dependent rate *λ*(*t*). The exposed individuals become infectious at a time rate *σ*. A fraction *α* of the infectious population is reported *I**r* (in red), while a fraction (1 − *α*) is unreported *I**u* (in purple). The infectious individuals that recover, reported or not, become recovered *R* (in green) at a time rate *γ*. Finally, it is assumed that a fraction *ϕ* of the removed population *γI**r* deceases *D**r* (in dark gray). ![Formula][5] The constant of proportionality *χ* involved in this relation will be explicitly computed through the calibration between the contact tracing and the compartmental models. Each consecutive time window has a width of 22 days and a step size of 5 days. We emphasize that our model has an intrinsic time delay regarding the consolidation of ![Graphic][6] that can reach ≈ 53 days. This value is associated to the time delay in the consolidation of COVID-19 dataset (≈ 15 days) and to the superposition of the maxima of two infectious periods and one exposed period. ### B. Compartmental Model We also adopt a compartmental model to describe the transmission of COVID-19 in order to estimate the levels of infection of the pathogen in Fortaleza. Here, we propose a SEIIR model that distinguishes the populations of Susceptible, Exposed, Infectious (reported or unreported), and Removed (recovered or deceased) individuals, as shown in Fig. 1b. Our model is inspired by the SEIIR model proposed by Li *et al*. [32]. The reported infectious population *I**r* corresponds to the number of individuals that had the SARS-CoV-2 infection confirmed by the health system. The unreported infectious population *I**u* comprises the complement of *I**r*, *i*.*e*., individuals that were infected with COVID-19 but remained unknown to health authorities. We assume that the large majority of the reported infectious individuals are symptomatic cases, in contrast to the population of unreported infectious individuals -of which the large majority is assumed to be of asymptomatic cases. Given this fundamental assumption and considering the recent finding that asymptomatic people are 42% less likely to transmit the SARS-CoV-2 than symptomatic ones [33], we define that the transmission rate for the unreported infectious population *I**u* is reduced by a factor of *µ* in relation to the parameter *β* that represents the transmission rate for the reported infectious population *I**r*. In this context, the time-dependent rate at which the susceptible population *S* becomes the exposed population *E* is given by ![Formula][7] where *N* is the total population of Fortaleza, taken as constant, being approximately equal to 2.67 million people. A fraction *α* of the exposed individuals is presumed to become reported infectious at a rate *σ*, and the complementary fraction (1 − *α*) to evolve to unreported infected at the same rate. Also, both reported and unreported infectious population are assumed to become part of the removed population at the same rate *γ*. We also keep track of the fraction *ϕ* of the removed reported infectious population evolving to death, so that the reported deceased population *D**r* increases at a rate of *ϕγI**r*. The following system of coupled differential equations rules our model: ![Formula][8] ![Formula][9] ![Formula][10] ![Formula][11] ![Formula][12] ![Formula][13] The total population *N* = *S* + *E* + *I**r* + *I**u* + *R* + *D**r* is conserved. Furthermore, it can be readily shown [32] that the effective reproduction number *Re**city* is given by ![Formula][14] From Eq. (9), we can identify ![Graphic][15] as the average number of secondary infections due to contagion with reported infectious individuals, while ![Graphic][16] is the effective reproduction number due to contagion with unreported infectious individuals. Finally, the SEIIR model is used here as a core model within the Iterative Ensemble Kalman Filter (IEnKF) framework (see Methods). This approach allows us to investigate the time evolution of the effective reproduction number *Re**city* by inferring the SEIIR model parameters and their populations (see Figs. S1 and S2 of the Supplementary Information). The IEnKF framework is systematically applied to running windows of 22 days, with step size of 5 days, starting from March 24 to November 9, 2020. We use as observable the cumulative number of deaths by SARS-CoV-2 reported daily by the health authorities. For the first and subsequent windows, the guesses for the initial populations of exposed, *E*, and deceased individuals, *D*, are obtained from the daily number of COVID-19 confirmed cases and the daily cumulative number of COVID-19 confirmed deaths. As the cumulative number of deaths by SARS-CoV-2, both quantities are calculated from the dataset of COVID-19 confirmed cases and deaths (see Methods). In the particular case of the initial guess for the exposed population, *E* = *C**i**/*(*ασ*), where *C**i* is the reported number of daily cases. This corresponds, for example, to 4, 982 individuals in the first window. After using IEnKF to estimate the values of all model parameters for the first window, the factor *Re**city* is calculated at its center. These parameters and all populations obtained by numerical integration of Eqs. (3)-(8), except for *E* and *D*, as previously explained, are used as initial guesses for the second window. The same procedure is then repeated for the third and subsequent windows. ## III. RESULTS AND DISCUSSION Figure 2a shows the normalized moving averages of bus validations of individuals that got COVID-19, including healthcare workers. We note that healthcare workers that came into contact with SARS-CoV-2 during the studied period did not reduce their bus rides as much as other passengers. In addition, their normalized moving averages of bus validations are getting closer to each other again as the economic reopening progresses. The inset of Fig. 2a shows the daily bus validations, which gradually started to increase in the economic reopening. Figures 2b and 2c show the daily numbers of cases and deaths, respectively, following the same previous normalization and stratification. ![Fig. 2.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/06/03/2021.05.31.21257245/F2.medium.gif) [Fig. 2.](http://medrxiv.org/content/early/2021/06/03/2021.05.31.21257245/F2) Fig. 2. Time series and moving averages of bus validations, COVID-19 cases, and COVID-19 deaths. Time evolutions of the normalized moving averages of **(a)** bus validations, COVID-19 cases, and **(c)** deaths, for healthcare workers (in red) and all individuals (in black). The insets show their corresponding daily numbers. In **(a)**, we note that healthcare workers that came into contact with SARS-CoV-2 during the studied period did not reduce their bus rides as much as other passengers. In addition, the normalized moving averages of bus validations of healthcare workers and of all individuals are getting closer to each other again as the economic reopening progresses. In **(b)** and **(c)**, we find that both the normalized moving averages of cases and of deaths, respectively, for healthcare workers increased before those of all individuals until the lockdown regime. The windows of moving averages have 5 days of width for all curves. We normalized each moving average by its maximum. The vertical dotted lines represent the beginning of social isolation (State Decree 33,519), lockdown (State Decree 33,574), and economic reopening (State Decree 33,608) regimes imposed on March 20, May 8, and June 1, 2020, respectively. We also highlight, in light red, the lockdown period in the city of Fortaleza. The representativeness of the dataset of COVID-19 confirmed cases on buses is assessed comparing the daily numbers of infectious individuals within those vehicles and in the entire city, as shown in Fig. 3. The daily number of infectious individuals is computed taking into account the 14 days that the individuals remain infectious, *i*.*e*., each individual who tested positive for COVID-19 counts up to 14 times, once per day, for the infectious curve. Figures 3a and 3b show the daily numbers regarding all infectious passengers and those infectious passengers who are healthcare workers, respectively. Similarly, the daily numbers of infectious individuals and infectious healthcare workers of the entire city are shown in Figs. 3c and 3d, respectively. While the first and second waves of the epidemic can be clearly identified in both curves shown in Fig. 3a (infectious passengers) and Fig.3c (infectious individuals), only highly attenuated peaks during the second wave period can be visualized in the corresponding curves for healthcare workers, as shown in Figs. 3b and 3d. We conjecture that the explanation for this behavior may be twofold. First, due to the high contagion of healthcare workers during the first wave, this group of people may have achieved a large percentage of immunity, as compared to the rest of the population. Second, efficient Personal Protective Equipment (PPE) became more available in hospitals after the first wave. The results in Fig. 3e show that the percentage of infectious passengers with respect to all infectious individuals in Fortaleza was higher than 1% during most of the epidemic period. Finally, the evolution in time of the fraction between infectious passengers and infectious individuals in the city who are both healthcare workers is shown in Fig. 3f. ![Fig. 3.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/06/03/2021.05.31.21257245/F3.medium.gif) [Fig. 3.](http://medrxiv.org/content/early/2021/06/03/2021.05.31.21257245/F3) Fig. 3. Representativeness of the dataset of COVID-19 confirmed cases on buses. The daily numbers of **(a)** all infectious passengers, **(b)** infectious passengers who are healthcare workers, all infectious individuals in the entire city, and **(d)** infectious individuals who are healthcare workers in the entire city. For healthcare workers, we highlight an unexpected emergence of a single peak in the daily numbers of infectious individuals within buses and in the entire city, which contrasts to the first and the second waves of COVID-19. We conjecture that the explanation for this behavior may be the lack of Personal Protective Equipment (PPE) in hospitals during the first wave or the herd immunity of healthcare workers during the second wave. **(e)** The fraction of all infectious passengers. **(f)** The fraction of infectious passengers who are healthcare workers. These results show that the percentage of infectious passengers with respect to all infectious individuals in Fortaleza was higher than 1% during most of the epidemic period (dashed gray line). All solid lines represent moving averages with windows of 7 days. The histogram of the values of *τ**ij* for the network of PICs is shown in Fig. 4a. The obtained distribution is characterized by the average ⟨*τ**ij* ⟩ *PICs* ≈ 28 minutes. We find that CCs, defined by *τ**ij* *> τ**c* = 15 minutes, represent about 62% of the PICs, as shown by the Complementary Cumulative Distribution Function (CCDF) in the inset of Fig. 4a. For the network of CCs, the average of the shared ride times is ⟨*τ**ij* ⟩*CCs* ≈ 39 minutes. Fig. 4b shows the network of CCs taking into consideration the periods of social isolation, lockdown, and economic reopening. As depicted, it is composed of several trees, where the vertices represent bus passengers that were diagnosed with COVID-19 and the edges correspond to CCs. Bus passengers identified as healthcare workers in the network are highlighted in red. The size of the vertices is proportional to their outdegrees. ![Fig. 4.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/06/03/2021.05.31.21257245/F4.medium.gif) [Fig. 4.](http://medrxiv.org/content/early/2021/06/03/2021.05.31.21257245/F4) Fig. 4. Shared Ride Time Histogram *τ**ij* and Network of Close Contacts (CCs). (**a**) We show the weight distribution *τ**ij* of the network of Potentially Infectious Contacts (PICs) in 15-minute-width bins. The average of the shared ride times of PICs is ⟨ *τ**ij* ⟩ *PICs* ≈ 28 minutes. Applying the threshold *τ**c* = 15 minutes in the network of PICs, we define a network of Close Contacts (CCs). For the network of CCs, the average of the shared ride times is ⟨ *τ**ij* ⟩ *CCs* ≈ 39 minutes. The inset shows the Complementary Cumulative Distribution Function (CCDF) of *τ**ij*. We find that the percentage of the edges with *τ**ij* greater than *τ**c* = 15 minutes is ≈ 62% (dashed line), *i*.*e*., most part of PICs are CCs. (**b**) The vertices represent the bus passengers that were diagnosed with COVID-19 and the edges corresponds to the CCs each passenger had using the public transportation system in Fortaleza. The healthcare and non-healthcare workers are represented by red and gray vertices, respectively. The size of the vertices is proportional to their outdegrees. At this point, we show that it is possible to perform a direct comparison between the computed values of ![Graphic][17] obtained from the contact tracing model for different time windows and the corresponding effective reproduction numbers *Re**city* estimated from the compartmental model. First, it is reasonable to assume that ![Graphic][18], as long as the population traveling by public buses can be considered as statistically equivalent, from an epidemiologic point of view, to the rest of the city. As a consequence of this assumption and using Eq. (9), we can write that ![Formula][19] where the parameter *ψ* = [*α* + (1 − *α*)*µ*]−1 depends on the time window used for model inference with the IEnKF technique. We now proceed with the comparison between contact tracing and compartmental models. In practical terms, this is achieved by upscaling ![Graphic][20] to the numerical values obtained for *Re**city* during the early period of the SARS-CoV-2 epidemic, before the restrictions of isolation and social distancing imposed by the State Government took effect. Considering Eqs. (1) and (10), we use the relation ![Graphic][21] and the numerical values of ![Graphic][22], and *Re**city* on the day that corresponds to the maximum of *Re**city* during the first wave (April 8, 2020) to calculate *χ* ≈ 37. This constant combined with the values of *ψ* from the inference with the compartmental model, and the values of ![Graphic][23] from the contact tracing, both calculated for all time windows, are then used to obtain the entire curve of ![Graphic][24] (see the time evolution of *χ/ψ* in Fig. S3 of the Supplementary Information). The value of *χ* can be understood as the product of two factors, *χ* = *χ**rr**χ**ru*. Assuming the equality between the proportions of CCs in the pairs of infectious and exposed passengers with existing and missing values of *τ**ij*, we estimate *χ**rr* ≈ 1*/*0.586 ≈ 1.70 as a balance factor for possible missing CCs. The value of the remaining factor *χ**ru* ≈ 21.76 expresses the sub-notification of the confirmed cases as well as our lack of knowledge on the transmission from reported to unreported infectious passengers, for which the factor (1 − *α*)*/α* could be a lower bound (see Fig. S2 of the Supplementary Information for a sensitivity analysis of the model with parameter *α*), approximately between 5 and 9 [32]. In an entirely similar fashion, by considering only reported infectious passengers that can be identified as healthcare workers, we can estimate their particular effective reproduction number as ![Graphic][25] with the same upscaling factor *χ/ψ* used for all infectious passengers and ![Graphic][26] is the average of the vertices outdegrees for healthcare workers. In Figure 5, we show the comparison between the estimates of *Re**bus* and *Re**city* from March to November 2020. Although the contact tracing and compartmental models are defined on different scales, the former on a microscopic scale and the latter on a macroscopic scale, the two curves capture the same decreasing trend associated to both social isolation and lockdown periods. We note that *Re**bus* consistently follows *Re**city* during the local COVID-19 outbreak, except for a three-month period between the first and the second waves of daily cases. In this period, the *Re**bus* decayed to undetectable standards despite the fact that the number of daily bus validations has increased (see Fig. 2a). As also shown in Fig. 5, *Re**health* was systematically higher than *Re**bus*, which unveils that the healthcare workers played an important role in the transmission within buses during the first wave of COVID-19 in Fortaleza. Furthermore, *Re**health* remained undetectable even in the beginning of the second wave, in contrast to *Re**bus* and notwithstanding the increase of the number of daily bus validations of healthcare workers, as shown in Fig. 2a. As shown in the inset of Fig. 5, the maximum ratio *Re**health**/Re**bus* occurred soon after the lockdown period, since the hospitals were still overloaded due to the peak of cases at the beginning of May and the new daily infections were low in the beginning of the reopening period. We emphasize that the complement of *Re**health*, due to non-healthcare workers, behaves similar to *Re**bus*. ![Fig. 5.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/06/03/2021.05.31.21257245/F5.medium.gif) [Fig. 5.](http://medrxiv.org/content/early/2021/06/03/2021.05.31.21257245/F5) Fig. 5. Time evolution of the effective reproduction numbers. Moving averages of the effective reproduction number for the entire city, *Re**city* (blue ∘), for buses, *Re*r*bus* (gray *×*), and for healthcare workers in the buses, *Re**health* (red +). We find that *Re**bus* consistently follows *Re**city* during the local COVID-19 outbreak, except for a three-month period between the first and the second waves of daily cases. We also show that *Re**health* was systematically higher than *Re**bus*, which unveils that the healthcare workers played an important role in the transmission within buses during the first wave of COVID-19 in Fortaleza. The inset shows that the maximum ratio *Re**health**/Re**bus* occurred soon after the lockdown period. The windows of moving averages have 22 days of width with step size of 5 days for all curves. In the period indicated by the shaded regions in the main plot and its inset, both *Re**bus* and *Re**health* decayed to undetectable standards, *i*.*e*., no CCs could be identified under the framework of our contact tracing approach. The vertical dotted lines represent the beginning of social isolation (State Decree 33,519), lockdown (State Decree 33,574), and economic reopening (State Decree 33,608) regimes imposed on March 20, May 8, and June 1, 2020, respectively. We also highlight, in light red, the lockdown period in the city of Fortaleza. ## IV. CONCLUSIONS In summary, two epidemiological models have been used in this work to understand the transmission on public transportation during the COVID-19 outbreak in Fortaleza, Ceará, Brazil. Whilst the compartmental model accounts for the transmission in the entire city (macroscopic scale), the contact tracing model has been used to estimate the transmission inside city buses (microscopic scale) through the concept of CCs. Both models were fed with real data of bus validations and of COVID-19 confirmed cases and deaths. Our results show that *Re**bus* consistently follows *Re**city* during the local COVID-19 outbreak, except for a three-month period between the first and the second waves of daily cases. Furthermore, the transmission from healthcare workers within buses until the end of July is characterized by a value of *Re**health* persistently greater than *Re**bus*. Healthcare workers, even the non-frontline professionals, are more likely to get and, consequently, spread the pathogen because their social network distances to individuals that tested positive for COVID-19 are very short compared to non-highly exposed workers. Despite being more tested, healthcare workers may not even know that they are infectious when they board a bus due to eventual time delays of the result of a COVID-19 test. Other groups of highly exposed people may affect the dynamics of dissemination of the virus in a similar way, *e*.*g*., education workers and police officers. Therefore, our results reinforce the worldwide claim that it is imperative to propose special policies to support displacement (or to avoid it) of highly exposed groups of people. Finally, we suggest that the intensity and the necessity of using public transportation by highly exposed groups must be seriously considered as a criterion to prioritize their vaccination. ## V. METHODS ### A. Datasets #### 1. Bus validations Most part of bus passengers in Fortaleza (≈ 94%) pay their bus fares with a smart card. Every time a passenger passes their card on a ticket gate of a bus, a validation record is created. The Fortaleza City Hall compiled and made available an anonymized dataset of bus validations with the following information: a citizen’s ID (a hash code), a vehicle ID (another hash code), the date and time of the validation record and the estimated ride time. The dataset ranges from March to December 2020, totaling 107, 488, 528 validation registers that refers to 1, 426, 569 different passengers. #### 2. COVID-19 confirmed cases and deaths The dataset of COVID-19 confirmed cases and deaths is an anonymized list of all individuals diagnosed with the disease in Fortaleza from March to December 2020. These data were also processed and made available by the Fortaleza City Hall. Such dataset is organized in columns as follows: a citizen’s ID (the same hash code used in the previous dataset), the date of OS, a confirmed death flag, the date of death and a healthcare worker flag. In the period of time ranged by the data, there are 85, 553 confirmed cases (5, 960 of healthcare workers) and 3, 075 confirmed deaths (227 of healthcare workers). We emphasize that these healthcare workers are not only the frontline professionals but also people whose jobs are related to the health field. Finally, we found that 9, 032 people (721 healthcare workers) were diagnosed with COVID-19 and used their smart card on buses at least once from March to December 2020. ### B. Iterated Ensemble Kalman Filter We use the Iterated Ensemble Kalman Filter (IEnKF) framework [32, 35–37] to infer the compartmental model parameters and initial subpopulations. The algorithm is based on comparing predictions of the model *f* (.) obtained by the numerical integration of Eqs. (3)-(7) of the main text with a set of *T* observations *𝒪*1, …, *𝒪**T* taken at discrete times *t*1, …, *t**T* within an observation window (see Fig. S4 of the Supplementary Information). The inference framework starts from an initial state vector *X*(0) = *{S, E, I**r*, *I**u*, *R, D**r**}*(0), and an initial parameter vector *θ*(0) = *{β, µ, σ, γ, α, ϕ}*(0). To these vectors, uncertainties are attributed in terms of the variance matrices *σ**X* and *σ**θ*, respectively. For each iteration *m*, an ensemble of *P* “particles” is generated such that each particle has the initial state at time *t* drawn from a multivariate normal distribution with mean *X*(*m*−1) and variance *a*(*m*−1)*σ**X*, where < *a* < 1 is a “cooling factor”. The initial state vector for particle *i* is denoted by *X*(*t*, *i*) = *𝒩* (*X*(*m*−1), *a*(*m*−1)*σ**X*). These states are also used to set *X**F* (*t*, *i*) which define the posterior distribution at time *t*. Analogously, each particle *i* has an initial parameter vector *θ*(*t*, *i*) = *𝒩* (*θ*(*m*−1), *b*(*m*−1)*σ**θ*), where 0 < *b* < 1 is another cooling factor. The inference proceeds by numerically integrating the model from these initial conditions, such that the predicted vector state for each particle *i* at time *t**n* is obtained from the following distribution, *X**P* (*t**n*, *i*) = *f* (*X**F* (*t**n*−1, *i*), *θ*(*t**n*−1, *i*)). Based on these predictions, a weight *W* (*t**n*, *i*) is assigned to each particle *i*, such that ![Formula][27] where *𝒪* (*t**n*, *i*) is the predicted value for the observed quantities at time *t**n* for particle *i*, and Θ is a “temperature”. In our case, *𝒪* (*t**n*, *i*) is the prediction for the cumulative number of daily reported deaths *D**r*(*t**n*, *i*). The filtering process is accomplished by keeping the particles with the largest weights with probability *𝒫* = *W* (*t**n*, *i*)*/* ∑ *j* *W* (*t**n*, *j*). The states of the filtered particles will set the posterior distribution at time *t**n*, *X**F* (*t**n*, *i*) = *X**P* (*t**n*, *i**best*), where *i**best* is the index of the filtered particles [36]. The parameter vector is updated at time *t**n* using *θ*(*t**n*, *i*) = *𝒩* (*θ*(*t**n*−1, *i**best*), *b*(*m*−1)*σ**θ*). This filtering process continues until all the observations *𝒪*1, …, *𝒪**N* are compared. The iterative process continues by setting the initial state vector *X*(*m*) and parameter vector *θ*(*m*) for the next iteration. The next parameter vector is given by [35]: ![Formula][28] where ![Graphic][29] is the sample mean of *θ*(*t**n*, *i**best*) and *V* (*t**n*) is the variance [35, 36]. The next state vector is given by the sample mean, ![Formula][30] After each iteration *m*, the initial state vector *X*(*m*) and parameter vector *θ*(*m*) are used to compute the evolution of the model for the whole observation window 1, …, *t**T*. The performance of the inferred model is computed by evaluating the error ![Formula][31] The iteration continues until | *ε*(*m*) − *ε*(*m*−1) | < *ϵ**max*, where the threshold used here is *ϵ**max* = 0.01. In Table S1 of the Supplementary Information, we report the values of the estimated parameters for consecutive time windows. ## Supporting information Supplementary Information [[supplements/257245_file03.pdf]](pending:yes) ## Data Availability Two datasets were used with the approval and consent obtained by the Fortaleza City Hall, Ceará, Brazil. The first is a list of COVID-19 confirmed cases and deaths of patients in Fortaleza and the second consists of bus validations records from smart cards of passengers, both collected during the period from March to December, 2020. In the context of the ongoing health crisis, we make available these anonymized datasets under request to E. Oliveira at erneson{at}unifor.br. ## VII. CONTRIBUTIONS C.P., H.A.C., E.A.O., C.C., A.S.L., J.S.A. and V.F. designed research; C.P., H.A.C., E.A.O., C.C., A.S.L., J.S.A. and V.F. performed research; C.P., H.A.C., E.A.O., C.C., A.S.L., J.S.A. and V.F. analyzed data; and C.P., H.A.C., E.A.O., C.C., A.S.L., J.S.A. and V.F. wrote the paper. All authors reviewed the manuscript. ## VIII. COMPETING INTERESTS The authors declare no competing financial interests. ## IX. DATA AVAILABILITY This study was approved by the Institutional Review Board (IRB) at Universidade de Fortaleza (UNIFOR). Two datasets were used with the approval and consent obtained by the Fortaleza City Hall, Ceará, Brazil. The first is a list of COVID-19 confirmed cases and deaths of patients in Fortaleza and the second consists of bus validations records from smart cards of passengers, both collected during the period from March to December, 2020. In the context of the ongoing health crisis, we make available these anonymized datasets under request to E. Oliveira at erneson{at}unifor.br. ## VI. ACKNOWLEDGMENTS We gratefully acknowledge CNPq, CAPES, FUNCAP, the National Institute of Science and Technology for Complex Systems in Brazil and the Edson Queiroz Foundation for financial support. * Received May 31, 2021. * Revision received May 31, 2021. * Accepted June 2, 2021. * © 2021, Posted by Cold Spring Harbor Laboratory This pre-print is available under a Creative Commons License (Attribution 4.0 International), CC BY 4.0, as described at [http://creativecommons.org/licenses/by/4.0/](http://creativecommons.org/licenses/by/4.0/) ## References 1. [1].COVID-19 transmission—up in the air. The Lancet. Respiratory Medicine (2020). [https://doi.org/10.1016/S2213-2600(20)30514-2](https://doi.org/10.1016/S2213-2600(20)30514-2) 2. [2].Lloyd-Smith, J. O., Schreiber, S. J., Kopp, P. E., & Getz, W. M. Superspreading and the effect of individual variation on disease emergence. Nature 438, 355–359 (2005). [https://doi.org/10.1038/nature04153](https://doi.org/10.1038/nature04153) [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nature04153&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=16292310&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F06%2F03%2F2021.05.31.21257245.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000233300200048&link_type=ISI) 3. [3].Lossio-Ventura, J. A. et al. DYVIC: DYnamic VIrus Control in Peru. In 2020 IEEE International Conference on Bioinformatics and Biomedicine (BIBM) 2264-2267 (2020). [https://doi.org/10.1103/PhysRevResearch.3.013163](https://doi.org/10.1103/PhysRevResearch.3.013163) 4. [4].Serafino, M. et al. Superspreading k-cores at the center of COVID-19 pandemic persistence. arXiv preprint arxiv:2103.08685 (2021). [https://doi.org/10.1101/2020.08.12.20173476](https://doi.org/10.1101/2020.08.12.20173476) 5. [5].Reyna-Lara, A. et al. Virus spread versus contact tracing: Two competing contagion processes. Physical Review Research, 3, 013163 (2021). [https://doi.org/10.1103/PhysRevResearch.3.013163](https://doi.org/10.1103/PhysRevResearch.3.013163) 6. [6].Hamner, L. et al. High SARS-CoV-2 Attack Rate Following Exposure at a Choir Practice - Skagit County, Washington, March 2020. MMWR Morb Mortal Wkly Rep 69, 606–610 (2020). [http://doi.org/10.15585/mmwr.mm6919e6](http://doi.org/10.15585/mmwr.mm6919e6) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F06%2F03%2F2021.05.31.21257245.atom) 7. [7].Majra, D., Benson, J., Pitts, J. & Stebbing, J. SARS-CoV-2 (COVID-19) superspreader events. Journal of Infection 82, 36–40 (2021). [https://doi.org/10.1016/j.jinf.2020.11.021](https://doi.org/10.1016/j.jinf.2020.11.021) 8. [8].Liu, Y. et al. Aerodynamic analysis of SARS-CoV-2 in two Wuhan hospitals. Nature 582, 557–560 (2020). [https://doi.org/10.1038/s41586-020-2271-3](https://doi.org/10.1038/s41586-020-2271-3) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F06%2F03%2F2021.05.31.21257245.atom) 9. [9].Lednicky, J. A. et al. Viable SARS-CoV-2 in the air of a hospital room with COVID-19 patients International Journal of Infectious Diseases 100, 476–482 (2020). [https://doi.org/10.1016/j.ijid.2020.09.025](https://doi.org/10.1016/j.ijid.2020.09.025) [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.ijid.2020.09.025&link_type=DOI) 10. [10].Kwon, K. et al. Evidence of Long-Distance Droplet Transmission of SARS-CoV-2 by Direct Air Flow in a Restaurant in Korea. Journal of Korean Medical Science 35, e415 (2020). [https://doi.org/10.3346/jkms.2020.35.e415](https://doi.org/10.3346/jkms.2020.35.e415) 11. [11].Böhmer, M. M. et al. Investigation of a COVID-19 outbreak in Germany resulting from a single travel-associated primary case: a case series. The Lancet Infectious Diseases bf 20, 920–928 (2020). [https://doi.org/10.1016/S1473-3099(20)30314-5](https://doi.org/10.1016/S1473-3099(20)30314-5) 12. [12].Kakimoto, K. et al. Initial Investigation of Transmission of COVID-19 Among Crew Members During Quarantine of a Cruise Ship - Yokohama, Japan, February 2020. MMWR Morb Mortal Wkly Rep 69, 312–313 (2020). [http://doi.org/10.15585/mmwr.mm6911e2](http://doi.org/10.15585/mmwr.mm6911e2) [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.15585/mmwr.mm6911e2&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F06%2F03%2F2021.05.31.21257245.atom) 13. [13].Mizumoto, K. & Chowell, G. Transmission potential of the novel coronavirus (COVID-19) onboard the diamond Princess Cruises Ship, 2020. Infectious Disease Modelling 5, 264–270 (2020). [https://doi.org/10.1016/j.idm.2020.02.003](https://doi.org/10.1016/j.idm.2020.02.003) 14. [14].Vuchic, V. R. Urban transit: operations, planning, and economics (John Wiley & Sons, Hoboken, 2017). 15. [15].Stoddard, S. T. et al. The role of human movement in the transmission of vector-borne pathogens. Plos Neglected Tropical Diseases 3, e481 (2009). [https://doi.org/10.1371/](https://doi.org/10.1371/) journal.pntd.0000481 16. [16].Edelson, P. J. & Phypers, M. TB transmission on public transportation: A review of published studies and recommendations for contact tracing. Travel Medicine and Infectious Disease 9, 27–37 (2011). [https://doi.org/10.1016/j.tmaid.2010.11.001](https://doi.org/10.1016/j.tmaid.2010.11.001) 17. [17].Bomfim, R. et al. Predicting dengue outbreaks at neighbourhood level using human mobility in urban areas. Journal of the Royal Society Interface 17, 20200691 (2020). [https://doi.org/10.1098/rsif.2020.0691](https://doi.org/10.1098/rsif.2020.0691) 18. [18].Kraemer, M. U. G. et al. The effect of human mobility and control measures on the COVID-19 epidemic in China. Science 368 493–497 (2020). [https://doi.org/10.1126/science](https://doi.org/10.1126/science). abb4218 [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Mzoic2NpIjtzOjU6InJlc2lkIjtzOjEyOiIzNjgvNjQ5MC80OTMiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMS8wNi8wMy8yMDIxLjA1LjMxLjIxMjU3MjQ1LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 19. [19].Schlosser, F. et al. COVID-19 lockdown induces disease-mitigating structural changes in mobility networks. Proceedings of the National Academy of Sciences 117, 32883–32890 (2020). [https://doi.org/10.1073/pnas.2012326117](https://doi.org/10.1073/pnas.2012326117) [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoicG5hcyI7czo1OiJyZXNpZCI7czoxMjoiMTE3LzUyLzMyODgzIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjEvMDYvMDMvMjAyMS4wNS4zMS4yMTI1NzI0NS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 20. [20].Melo, H. P. et al. Heterogeneous impact of a lockdown on inter-municipality mobility. Physical Review Research 3, 013032 (2021). [https://doi.org/10.1103/PhysRevResearch.3.013032](https://doi.org/10.1103/PhysRevResearch.3.013032) 21. [21].Di Carlo, P. et al. Air and surface measurements of SARS-CoV-2 inside a bus during normal operation. Plos One 15, e0235943 (2020). [https://doi.org/10.1371/journal.pone.0235943](https://doi.org/10.1371/journal.pone.0235943) 22. [22].Goscé, L. & Johansson, A. Analysing the link between public transport use and airborne transmission: mobility and contagion in the London underground. Environmental Health 17, 1–11 (2018). [https://doi.org/10.1186/s12940-018-0427-5](https://doi.org/10.1186/s12940-018-0427-5) 23. [23].Jenelius, E. & Cebecauer, M. Impacts of COVID-19 on public transport ridership in Sweden: Analysis of ticket validations, sales and passenger counts. Transportation Research Interdisciplinary Perspectives 8, 100242 (2020). [https://doi.org/10.1016/j.trip.2020.100242](https://doi.org/10.1016/j.trip.2020.100242) 24. [24].Shen, Y. et al. Community outbreak investigation of SARS-CoV-2 transmission among bus riders in eastern China. JAMA internal medicine 180, 1665–1671 (2020). [https://doi.org/10.1001/jamainternmed.2020.5225](https://doi.org/10.1001/jamainternmed.2020.5225) 25. [25].Shen, J. et al. Prevention and control of COVID-19 in public transportation: experience from China. Environmental pollution 266, 115291 (2020). [https://doi.org/10.1016/j.envpol.2020.115291](https://doi.org/10.1016/j.envpol.2020.115291) 26. [26].Zhang, Z. et al. Disease transmission through expiratory aerosols on an urban bus. Physics of Fluids 33, 015116 (2021). [https://doi.org/10.1063/5.0037452](https://doi.org/10.1063/5.0037452) 27. [27].Hu, M. et al. Risk of coronavirus disease 2019 transmission in train passengers: An epidemiological and modeling study. Clinical Infectious Diseases 72, 604–610 (2021). [https://doi.org/10.1093/cid/ciaa1057](https://doi.org/10.1093/cid/ciaa1057) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F06%2F03%2F2021.05.31.21257245.atom) 28. [28].He, X. et al. Temporal dynamics in viral shedding and transmissibility of COVID-19. Nature Medicine 26, 672–675 (2020). [https://doi.org/10.1038/s41591-020-0869-5](https://doi.org/10.1038/s41591-020-0869-5) [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.7326/M20-3012&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F06%2F03%2F2021.05.31.21257245.atom) 29. [29].Li, Q. et al. Early transmission dynamics in Wuhan, China, of novel coronavirus-infected pneumonia. N. Engl. J. Med. 382, 1199–1207 (2020). [https://doi.org/10.1056/NEJMoa2001316](https://doi.org/10.1056/NEJMoa2001316) [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1056/NEJMoa2001316&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F06%2F03%2F2021.05.31.21257245.atom) 30. [30].Sanche, S. et al. High Contagiousness and Rapid Spread of Severe Acute Respiratory Syndrome Coronavirus 2. Emerging Infectious Diseases 26, 1470–1477 (2020). [https://doi.org/10.3201/eid2607.200282](https://doi.org/10.3201/eid2607.200282) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F06%2F03%2F2021.05.31.21257245.atom) 31. [31].Centers for Disease Control and Prevention. Available at [https://www.cdc.gov/coronavirus/2019-ncov/php/contact-tracing/contact-tracing-plan/contact-tracing.html](https://www.cdc.gov/coronavirus/2019-ncov/php/contact-tracing/contact-tracing-plan/contact-tracing.html) Accessed May 14, 2021. 32. [32].Li, R. et al. Substantial undocumented infection facilitates the rapid dissemination of novel coronavirus (SARS-CoV-2). Science 368, 489–493 (2020). [https://doi.org/10.1126/science.abb3221](https://doi.org/10.1126/science.abb3221) [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Mzoic2NpIjtzOjU6InJlc2lkIjtzOjEyOiIzNjgvNjQ5MC80ODkiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMS8wNi8wMy8yMDIxLjA1LjMxLjIxMjU3MjQ1LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 33. [33].Byambasuren, O. et al. Estimating the extent of asymptomatic COVID-19 and its potential for community transmission: Systematic review and meta-analysis. Journal of the Association of Medical Microbiology and Infectious Disease Canada 5, 223–234 (2020). [https://doi.org/10.3138/jammi-2020-0030](https://doi.org/10.3138/jammi-2020-0030) 34. [34].Le Marshall, J., Rea, A., Leslie, L., Seecamp, R. & Dunn, M. Error characterisation of atmospheric motion vectors. Aust. Meteorol. Mag. 53, 123–131 (2004). 35. [35].Ionides, E. L., Bretó, C. & King, A. A. Inference for nonlinear dynamical systems. Proc. Natl. Acad. Sci. U. S. A. 103, 18438–18443 (2006). [https://doi.org/10.1073/pnas.0603181103](https://doi.org/10.1073/pnas.0603181103) [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoicG5hcyI7czo1OiJyZXNpZCI7czoxMjoiMTAzLzQ5LzE4NDM4IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjEvMDYvMDMvMjAyMS4wNS4zMS4yMTI1NzI0NS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 36. [36].King, A. A., Ionides, E. L., Pascual, M. & Bouma, M. J. Inapparent infections and cholera dynamics. Nature 454, 877–880 (2008). [https://doi.org/10.1038/nature07084](https://doi.org/10.1038/nature07084) [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nature07084&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=18704085&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F06%2F03%2F2021.05.31.21257245.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000258398600034&link_type=ISI) 37. [37].Sakov, P., Oliver, D. S. & Bertino, L. An Iterative EnKF for Strongly Nonlinear Systems. Monthly Weather Review 140, 1988–2004 (2012). [https://doi.org/10.1175/MWR-D-11-00176.1](https://doi.org/10.1175/MWR-D-11-00176.1) [1]: /embed/inline-graphic-1.gif [2]: /embed/inline-graphic-2.gif [3]: /embed/inline-graphic-3.gif [4]: /embed/inline-graphic-4.gif [5]: /embed/graphic-2.gif [6]: /embed/inline-graphic-5.gif [7]: /embed/graphic-3.gif [8]: /embed/graphic-4.gif [9]: /embed/graphic-5.gif [10]: /embed/graphic-6.gif [11]: /embed/graphic-7.gif [12]: /embed/graphic-8.gif [13]: /embed/graphic-9.gif [14]: /embed/graphic-10.gif [15]: /embed/inline-graphic-6.gif [16]: /embed/inline-graphic-7.gif [17]: /embed/inline-graphic-8.gif [18]: /embed/inline-graphic-9.gif [19]: /embed/graphic-14.gif [20]: /embed/inline-graphic-10.gif [21]: /embed/inline-graphic-11.gif [22]: /embed/inline-graphic-12.gif [23]: /embed/inline-graphic-13.gif [24]: /embed/inline-graphic-14.gif [25]: /embed/inline-graphic-15.gif [26]: /embed/inline-graphic-16.gif [27]: /embed/graphic-16.gif [28]: /embed/graphic-17.gif [29]: /embed/inline-graphic-17.gif [30]: /embed/graphic-18.gif [31]: /embed/graphic-19.gif