Critical Role of the Subways in the Initial Spread of SARS-CoV-2 in New York City

Jeffrey E. Harris
doi: https://doi.org/10.1101/2021.07.03.21259973

Abstract

We studied the possible role of the subways in the spread of SARS-CoV-2 in New York City during late February and March 2020. Data on cases and hospitalizations, along with phylogenetic analyses of viral isolates, demonstrate rapid community transmission throughout all five boroughs within days. The near collapse of subway ridership during the second week of March was followed within 1-2 weeks by the flattening of COVID-19 incidence curve. We observed persistently high entry into stations located along the subway line serving a principal hotspot of infection in Queens. We used smartphone tracking data to estimate the volume of subway visits originating from each zip code tabulation area (ZCTA). Across ZCTAs, the estimated volume of subway visits on March 16 was strongly predictive of subsequent COVID-19 incidence during April 1-8. In a spatial analysis, we distinguished between the conventional notion of geographic contiguity and a novel notion of contiguity along subway lines. We found that the March 16 subway-visit volume in subway-contiguous ZCTAs had an increasing effect on COVID-19 incidence during April 1-8 as we enlarged the radius of influence up to 5 connected subway stops. By contrast, the March 31 cumulative incidence of COVID-19 in geographically-contiguous ZCTAs had an increasing effect on subsequent COVID-19 incidence as we expanded the radius up to 3 connected ZCTAs. The combined evidence points to the initial citywide dissemination of SARS-CoV-2 via a subway-based network, followed by percolation of new infections within local hotspots.

Main

An accurate, thorough understanding of the rapid, widespread propagation of SARS-CoV-2 infection during the early phase of the massive outbreak in New York City is crucial to the successful control of future pandemic threats.

To that end, we test three main hypotheses here. First, New York City’s extensive public transport system, particularly its subways, played a critical role in the widespread dissemination of SARS-CoV-2 infection throughout the city during the end of February and the beginning of March 2020. Second, the ensuing marked decline in subway use was an important vehicle by which the public’s growing perception of risk was translated into reduced community transmission of the virus. Third, those areas with an attenuated decline in subway use, we posit, subsequently became the loci for high-density clusters of viral infection in late March 2020.

The Metropolitan Transportation Authority (MTA), a network of subways, buses and commuter rail cars serving the NYC area, is larger than all of other metropolitan transport systems in the United States combined. While nearly 85 percent of U.S. workers drive to their jobs, according to the MTA, 80 percent of rush-hour commuters to the city’s central business districts use transit.1 The MTA’s subway system is particularly unique, with a total of 1697.8 million turnstile entries during the calendar year 2019,2 compared to 157.2 million entries into the Washington DC metro,3 the next largest subway system in the country.

Our hypotheses are hardly novel. The role of transportation networks in the spread of SARS-CoV-2 has been supported by recent studies of the initial outbreak in Wuhan, China.4-6 One study of the NYC epidemic found an association between continued subway use among essential workers and a delayed flattening of the epidemic curve.7 Another study based in part on NYC subway ridership data found a link between mobility and COVID-19 risk.8 Yet another study found strong correlations between NYC subway turnstile entries and COVID-19 cases and deaths.9

What sets our study apart is its comprehensive, multidisciplinary approach. We rely on such diverse lines of evidence as phylogenetic analysis of early viral samples, public health data on confirmed COVID-19 cases, public transport data on turnstile entries, location-tracking data on the movements of smartphones, and census data on the prevalence of at-risk multi-generational households. Our spatial analysis of emerging case clusters distinguishes critically between the conventional notion of geographic contiguity and what we call subway contiguity.

Early Rapid, Widespread Community Transmission

Assessment of the extent of infection during the earliest days of the NYC outbreak has been hampered by the initial lack of adequate testing materials. Still, Fig. 1a shows that, despite the narrow testing criteria initially imposed on February 28 by the Centers for Disease Control (CDC),10 positive tests had been detected in residents of every borough of the city by March 6. Fig. 1b further demonstrates that by March 1, hospitals had already admitted patients residing in every borough. The incubation period between infection and first symptoms of COVID-19 is 5 days on average,11 with a range of up to 2 weeks.12 Add to that elapsed time an extra 4–10 more days before a symptomatic individual becomes sick enough to be hospitalized.13 Accordingly, in all likelihood, SARS-CoV-2 infections were already occurring by mid-February in every one of the five boroughs of a city of over 8 million inhabitants. This pattern of early rapid, widespread dispersion is sharply distinguishable from the gradual radial geographic expansion of COVID-19 cases observed in the earliest days of epidemic in Los Angeles County,14 a comparable sized jurisdiction with 10 million inhabitants.

Fig. 1:
Fig. 1: Evidence of Early Rapid, Widespread Community Transmission

a. Counts of the earliest cases of test-confirmed COVID-19 reported by the NYC health department, starting on February 29, 2020.20 The counts represent individuals initially identified through targeted testing of symptomatic persons in accordance with restricted criteria issued on February 28 by the U.S. Centers for Disease Control (CDC).10 The horizontal scale indicates the dates that the cases were diagnosed over the ensuing 8 days. b. Timeline of the numbers of individuals ultimately diagnosed with COVID-19 in connection with their inpatient hospitalizations, derived from the same data source.20 The counts of these hospitalization are graphed according to each individual’s date of admission during the same 9-day interval. c. Timing and locations of 78 viral isolates from dominant clade A2a collected from patients of the Mount Sinai Health System (MSHS) in New York.15 In addition to four of the New York City boroughs (Brooklyn, Bronx, Manhattan, and Queens), two of the MSHS A2a patients were from Westchester County (colored cyan) and five patients had unknown residence (colored white). Pink bubbles a cluster of 17 samples sharing a common point mutation, A1844V in open reading frame (ORF) 1a. d. Map of all subway lines and stops in NYC, distinguishing 129 zip code tabulation areas (ZCTAs) containing a subway station, 30 ZCTAs geographically contiguous with a ZCTA containing a subway station, and 29 other ZCTAs. The Jamaica – 179th Street station at the end of the F Line connects to the 43 bus-route running along Hillside Avenue, which terminates in ZCTA 11004.

The data in Fig. 1c help to distinguish between two alternative explanations for this pattern of early rapid, widespread dispersion of SARS-CoV-2 infections: parallel, contemporaneous importation from multiple outside sources; and rapid mixing via community transmission. The figure describes the timing and locations of 78 viral isolates belonging to phylogenetic clade A2a that were collected from patients of the Mount Sinai Health System (MSHS) in New York15 soon after the CDC liberalized its testing criteria.10 Within this dominant clade, the investigators identified a local transmission cluster with a signature mutation in samples drawn from residents of Brooklyn, Manhattan, Queens, and Westchester County over a 5-day period. This observation goes against exclusively parallel seeding from distinct sources.

The evidence from Figs. 1a–c alone does not identify the distinct mechanisms underlying such widespread community transmission in so short an interval. Despite a large body of investigation attempting to retrospectively track down super-spreader events,16 the only such documented occurrence is an outbreak of COVID-19 among MTA front-line workers.17,18 If only by exclusion, we are left with NYC’s unique subway system,1 which, in combination with the MTA’s extensive bus routes,19 covers virtually every corner of the city (Fig. 1d).

Subway Volume and COVID-19 Cases

The Collapse of Subway Travel and The Flattening of the Epidemic Curve

For the city as a whole, Fig. 2a compares daily subway turnstile entries21 to daily numbers of confirmed COVID-19 diagnoses.20 Counts of confirmed cases based on voluntary testing of symptomatic individuals are known to have significantly understated actual numbers of SARS-CoV-2 infections.22,23 Still, once the CDC liberalized its testing criteria,10 one can see the rapid growth in daily confirmed cases, from 21 on March 8 to 1,038 on March 15.

Fig. 2:
Fig. 2: Subway Volume and COVID-19 Cases

a. COVID-19 case counts and subway volume during February 23 – April 19, 2020. The dark green-colored data points show the numbers of daily, city-wide confirmed COVID-19 cases reported by the NYC health department,20 measured on a logarithmic scale at the left. The lilac-colored bars show the daily volume of trips on the city’s subway system, computed from the Metropolitan Transportation Authority (MTA) turnstile data21 and measured on a linear scale at the right. b. COVID-19 case counts in Manhattan and Queens during March 1 – April 5, 2020, with common, initial exponential growth at a doubling time of 1.1 during the week of March 8– 15, as estimated by Poisson regression, followed by divergence of the epidemic paths in the two boroughs. c. Zip code tabulation areas (ZCTAs) in New York City, color coded according to cumulative case incidence as of March 31, 2020, showing a high-incidence hot spot in the Queens-Elmhurst area. d. Section of Fig. 2c, overlaid by the locations of the 22 stations of the 7 (Flushing) subway line, including those in Manhattan (sky blue), the hot spot (yellow), and the remainder of Queens (pink).27 The 82nd Street – Jackson Heights station within the yellow group is identified for reference. The pink-colored Mets-Willets station within ZCTA 11368 is on the other side of Grand Central Parkway. e. Relative numbers of daily turnstile entries for each of the three zones of the 7 (Flushing) line identified in Fig. 2d. The daily turnstile entries, likewise derived from the MTA turnstile data 21, were normalized so that the volume on Monday, March 2 was equal to 100 for each zone. As of March 16, the subway entries into the yellow hotspot stations were 63.2 percent of their March 2 level, while entries into the remaining Queens stations and Manhattan stations were, respectively, 47.7 and 32.2 percent of their March 2 baseline

During that same week from March 8–15, subway volume was already declining from its prior average of 5.6 million turnstile entries per weekday. By the end of that week, daily COVID-19 case counts had begun to deviate from their exponential trend. By the time subway rides had fallen to less than one-quarter of their regular volume in the third week of March, the epidemic curve had flattened out. The data are compatible with a causal relation between the drop in subway demand and the deceleration of the epidemic curve, with a delay of 2 ± 1 weeks between the two time-series.

The flattening of the epidemic curve cannot be wholly attributable to official government actions to restrict mobility and reduce interpersonal contact. The decline in subway turnstile entries in Fig. 2a occurred before the mayor closed entertainment venues and limited restaurants, bars and cafes to food take-out and delivery on March 17.24 While the mayor indeed shut down nightclubs, movie theaters, and concert halls, no one ordered the subways closed. To the contrary, state and local officials attempted to quell the public’s rising fears about the risks of coronavirus transmission on public transit.25,26 A more plausible explanation is that voluntary action motivated by fear of contagion – and not a response government coercion – precipitated the collapse of subway demand, which at least in part contributed to the subsequent flattening of curve.

The Attenuated Decline in Subway Use and the Emergence of Hotspots

If the decline in subway use in fact caused the observed deceleration of the epidemic in Fig. 2a, then those areas of the city with a more rapid decline would experience a greater deceleration, while those areas with an attenuated decline would experience continued epidemic growth. This prediction is tested in Figs. 2b through 2e, where we focus on an emerging hotspot in the Elmhurst area of Queens and the specific subway line running through it.

Fig. 2b plots daily confirmed COVID-19 cases over time in two boroughs: Manhattan and Queens.20 During the week starting March 8, the case counts from both boroughs followed an exponential path with a slope of 0.63/day, which, based on a generation time of 5.5 days,28 implies a basic reproductive number ℛ0= 0.63 × 5.5 = 3.47. This estimate of ℛ0 is comparable to that estimated for the outbreak in Wuhan,29 a city with its own massive subway system,30 but higher than that estimated for Italy.31 By the third week in March, however, the two incidence curves began to diverge significantly. By the last full week of the month, weekday reported cases in Manhattan were down to about 600, while weekday reported cases in Queens exceeded 1,500.

Fig. 2c maps the cumulative incidence of confirmed COVID-19 cases according to zip code tabulation area (ZCTA) as of March 31, 2020. While there are isolated high-incidence ZCTAs in Brooklyn and the Bronx, there is a notable cluster in the Elmhurst area in Queens, especially ZCTAs 11369 and 11370, where the cumulative incidence of confirmed cases had already exceeded 1 percent of the population. Manhattan, by contrast, shows no foci of cumulative COVID-19 incidence in excess of 0.75 percent of the population. Comparison of the borough-level data in Fig. 2b suggests that the Queens-Elmhurst hotspot seen in Fig. 2c may have begun to emerge in the third week of March.

Fig. 2d displays the 22 stations of the 7 (Flushing) subway line27 overlaid on a section of the map of Fig. 2c. Fig. 2e shows that turnstile entries into the six yellow-colored stations within the Queens-Elmhurst hotspot remained significantly higher than the remaining pink Queens and blue Manhattan stations, especially from the week of March 15 onward. The divergence in the decline in subway volume among these three groups is consistent with the prediction that the attenuated decline in turnstile volume promoted continued epidemic spread of SARS-CoV-2.

Smartphone Device Movements and COVID-19 Incidence

Smartphone Device Movements as a Proxy for Subway Turnstile Entries

The turnstile volume data in Figs. 2a and 2e show how many riders entered the subway system at various stations, but not where these subway riders originated. To fill this data gap, we relied on data on the movements of smartphones equipped with location-tracking software.32

Each smartphone movement (or “visit”) had a recorded origin and a destination census block group (CBG). Fig. 3a illustrates how we reproduced the daily pattern of turnstile entries into each subway station by adding up those smartphone visits whose destination was the CBG where that station was located. That finding allowed us to rely upon smartphone visits to station CBGs as a proxy for turnstile entries, and thus to study the origins of subway visitors.

Fig. 3:
Fig. 3: Smartphone Device Movements and COVID-19 Cases

a. Daily turnstile entries into the 82nd Street – Jackson Heights Station (blue vertical bars, right axis) and numbers of smartphone device visits to the census block group (CBG) containing that station (red and green lines, left axis). The green data series shows the number of visits from devices originating within the same CBG, while the red data series shows red line shows the number of visits from devices originating in other CBGs. b. Section of Queens showing CBG boundaries within ZCTA boundaries, overlaid with locations of stations along the 7 (Flushing) subway line. Light-dark blue shading identifies those origin CBGs with the highest number of combined visits to a pair of destination CBGs along the 7 (Flushing) line: one of the yellow hotspot stations and one station in the Queensboro Plaza-Court Square commercial complex. In addition, those ZCTAs within the Queens-Elmhurst hot spot have been shaded light-dark green according to the two-color scheme of Figs. 2c and 2d. c. Number of device visits to subway CBGs on March 16, expressed as a percent of visits during March 1–7, 2020. d. Incidence of newly diagnosed COVID-19 cases during April 1–8, 2020. e. Incremental COVID-19 incidence during April 1–8 (vertical axis) related to the number of visits to subway CBGs on March 16, 2020 (horizontal axis). Each point in the log-log plot is an individual ZCTA. As in Fig. 3c, visit counts are normalized so that average volume during the first week of March equaled 100. The superimposed line is the ordinary least squares fit. (See Supplement.) f. Prevalence of at-risk multi-generational households, measured as the proportion of households in each ZCTA with at least 4 persons, of whom at least one person was 18–34 years of age and at least one other person was at least 50 years of age.14 The map shows census tract boundaries within ZCTA boundaries. Color scheme reflects quartiles of prevalence.

Fig. 3b illustrates how smartphones entering the 7 (Flushing) Line at the yellow stations (already identified in Fig. 2d) originated not only from the ZCTA where the station was located, but also from the adjacent high-incidence ZCTAs. This finding suggested that we could reliably estimate the number of subway visits originating from each ZCTA by adding up subway-station smartphone visits that originated from that ZCTA. To that end, Fig. 3c displays the resulting map of the estimated distribution of subway visits by ZCTA as of March 16, 2020, expressed as a percentage of the corresponding baseline volume during the first week of March.

Relation Between Subway Visits in Mid-March and COVID-19 Incidence in Early April

Fig. 3d maps the incidence of newly diagnosed COVID-19 cases per 10,000 during April 1–8, 2020. If smartphone visits are in fact a reliable proxy for subway visits, and if an attenuated decline in subway visits in certain areas of the city in mid-March resulted in the subsequent emergence of high-incidence hotspots in those areas by early April, then we should expect to observe a strong correlation between the visit volume mapped in Fig. 3c and the incidence mapped in Fig. 3d. This prediction is in fact borne out in the bivariate plot of Fig. 3e. The slope of the ordinarily least squares-fitted line was significant at the level p < 0.001.

The significant bivariate association in Fig. 3e held up in multivariate models that took into account of two additional ZCTA-specific covariates: (i) cumulative COVID-19 incidence through March 31 (already mapped in Fig. 2c); and (ii) the prevalence of multi-generational households (mapped in Fig. 3f), a well-established ecological determinant of transmission rates.14,33-35. In all such multivariate models, the estimated parameters were significantly different from zero at the level p = 0.006 or lower. (See Supplement Table A.)

Spatial Analysis

Geographic and Subway Contiguity

The foregoing multivariate models of COVID-19 incidence during the first week of April do not account for the possibility of contagion across ZCTAs. While models of the spatial propagation of SARS-CoV-2 across geographic units have been proposed and tested,14,36 New York City presents a potentially unique example of contagion in subway space, as opposed to geographic space.

To that end, consider two distinct ZCTAs, abstractly labeled i and j. We say that ZCTAs i and j are geographically contiguous, or simply g-contiguous, when they share at least one common boundary point. By contrast, the same two ZCTAs are subway contiguous, or simply s-contiguous, when ZCTA j is the next stop after ZCTA i on some subway line in some direction. As detailed in the section “Contiguity in Geographic and Subway Space” in the Supplement, these elemental relations between ZCTAs can be compounded. For example, two ZCTAs i and j are g2-contiguous if there is a third distinct ZCTA labeled k, such that ZCTA i is g-contiguous with ZCTA k and ZCTA k is in turn is g-contiguous with ZCTA j.

As a further extension of the concept of compound contiguity, we say that two ZCTAs are (g+g2)-contiguous if they are either g-contiguous or g2-contiguous. This situation is illustrated in Fig. 4a, which shows ZCTA 11415 (colored orange) in Queens, surrounded by a total of 18 ZCTAs (colored peach) that are (g+g2)-contiguous with ZCTA 11415. Within this group, four ZCTAs are g-contiguous with ZCTA 11145, while the remaining 14 ZCTAs are situated effectively within a radius of 2 from the reference ZCTA 11415.

Fig. 4:
Fig. 4: Spatial Analysis

a. Map of ZCTA 11415 (colored orange), surrounded by 18 ZCTAs (colored peach) within a geographic radius of 2 ZTCAs. b. Map of ZCTA 11415 (colored orange), along with 12 ZCTAs either within a geographic radius of 1 ZCTA or a subway radius of 5 station stops. c. Estimated spatial effects of cumulative incidence through March 31 and subway volume on March 16 in relation to radius in subway space, as the contiguity criterion varied from g, to g+s, to g+s+s2, up to g+s+s2+s3+s4+s5. Subway volume exhibited a significantly increasing trend. (In a 2-sided z-test comparing a radius of 3 with radius of 0, p = 0.026.) Cumulative incidence did not. d. Estimated spatial effects of cumulative incidence through March 31 and subway volume on March 16 in relation to radius in geographic space, as the contiguity criterion was varied from g to g+g2 to, g+g2+g3. Cumulative incidence exhibited a significantly increasing trend. (In a 2-sided z-test comparing a radius of 3 with radius of 1, p < 0.001.) Subway volume did not. (In an analogous 2-sided z test, p = 0.268.)

Fig. 4b, by contrast, displays 12 ZCTAs (again colored peach) that are (g+s+s2+s3+s4+s5)-contiguous with the reference ZCTA 11415 (again colored orange). ZCTA 11367 is exclusively g-contiguous with the reference ZCTA 1145. The remaining ZCTAs are accessible within 5 subway stops along the same or a connecting line. Thus, ZCTA 11101 in Queens is accessible via four stops on the E Line, while ZCTA 10065 in Manhattan is further accessible after a transfer at the Queens Plaza station to the R, N or W Lines.

Spatial Regressions: Subways, Networks, and Percolation

Our non-spatial regression models permitted us to measure how prior conditions in a particular ZCTA (including March 16 subway volume and March 31 cumulative cases) influenced subsequent COVID-19 incidence during April 1–8 within the same ZCTA. By contrast, our spatial regression models (detailed in the Supplement) permitted us to measure how COVID-19 incidence during the first week of April was influenced by prior conditions in other ZCTAs. To implement these spatial models, we did not arbitrarily allow each ZCTA to be influenced by all other ZCTAs, but instead restricted the radius of potential contagion in both geographic and subway space. Thus, for a particular ZCTA in Queens, Fig. 4a illustrates a limited geographic radius of 2 ZCTAs, while Fig. 4b illustrates a limited subway radius of 5 stops.

We repeatedly estimated such between-ZCTA spatial effects as we varied the allowable radius of influence – from 1 to 3 in geographic space, and from 0 to 5 in subway space. As we enlarged the allowable radius in subway space, as shown in Fig. 4c, we found that the volume of subway visits originating in other ZCTAs on March 16 was an increasingly strong predictor of subsequent COVID-19 incidence during the first week in April. Cumulative incidence in other ZCTAs as of March 31, however, showed no such trend in relation to the allowable radius in subway space. By contrast, as we enlarged the allowable radius of influence in geographic space, as shown in Fig. 4d, we found the reverse. Cumulative incidence in other ZCTAs as of March 31 became an increasingly strong predictor of subsequent COVID-19 incidence, whereas subway visits in other ZCTAs as of March 16 did not.

Our finding that subway volume as of March 16 exhibited increasingly contagious effects in subway space supports the conclusion that SARS-CoV-2 was being propagated via a subway-based network at least through March 16. Our finding that March 31 cumulative incidence exhibited increasingly contagious effects in geographic space supports the conclusion that percolation of new cases through local geographic spread had subsequently become the dominant mode of propagation by the end of March. Once local clusters developed, further percolation of new cases via transmission within multi-generational households (Fig. 3f) became dominant.

Discussion

The evidence presented here supports three distinct but not mutually exclusive hypotheses. First, the subway system played a critical role in the rapid, widespread community transmission of SARS-CoV-2 infection throughout New York City during late February and early March 2020. Second, the ensuing marked decline in subway travel was an important mechanism by which the public’s growing perception of risk was translated into reduced community transmission of the virus. Third, those areas with an attenuated decline in subway use subsequently became hotspots of viral infection in late March and early April 2020.

One alternative interpretation is that subway travel was no more than a proxy for other determinants of vulnerability to COVID-19. In higher-risk communities, so the argument goes, many residents had service jobs that could not be performed remotely. Such an interpretation, however, does not square with the spatial-effect findings in Fig. 4c, which imply some mechanism of contagion running along subway lines. A more responsive counterargument would have to assign at least an indirect role to the subway system. Thus, the decline in turnstile entries seen in Figs. 2a and 2e (and Supplement Fig. B) could have reflected employees’ responses to their employers’ requests to work from home, which in turn reduced workplace exposure, where contagion would in fact have taken place. This version does not require that infected individuals transmitted their infections inside subway cars or on station platforms. It concedes only that public transport was an efficient vehicle for moving infected individuals from the periphery of the city to its commercial centers and back again many times a day.

This last counterargument, however, does not square with the evidence on the known mechanisms of SARS-CoV-2 transmission. An infected person exhales moist air containing very small droplets loaded with the virus.37 A passenger without a mask standing two feet away from an infected rider without a mask for just 15 minutes would almost certainly have inhaled virus particles, even if the infected rider never coughed or sneezed.38 An infected person constantly sheds virus particles in the form of fomites on almost every surface he touches, such as glasses, keys and phones.39 That would include the stainless-steel poles shared by standing passengers. Social distancing can be difficult if not impossible in crowded subway cars and platforms, as well as in public transportation conveyances and transportation hubs generally.40 A crowded subway train or platform would thus have been an ideal incubator for coronavirus transmission. In a study of outbreaks involving three or more cases in municipalities in China outside Hubei Province, transport-based transmission was second only to home-based transmission.41 The extensive outbreak among MTA front-line workers (and later on, their family members) has no alternative explanation.17,18

Yet another counterargument draws upon conflicting studies of the transmission of other respiratory viruses in public transport. A study of the London Underground offered supporting evidence of the transmission of influenza-like illness.42,43 But a simulation study calibrated to the 1957–1958 flu epidemic in New York City estimated only a small contribution from subway travel.44 A cross-sectional study of 121 cities found a negative association between public transit use and mortality from pneumonia and influenza during 2006-2015.45 In contrast to a basic reproductive number of ℛ0= 3.47 (95% confidence interval, 3.16–3.78) for SARS-CoV-2 in New York City estimated here, seasonal influenza has an ℛ0 in the range of 1.2–1.4, while pandemic influenza has an ℛ0 in the range of 1.4–1.8, with the high end representing the 1918 pandemic.46 While a wave of COVID-19 cases swept through the U.S. during October 2020– January 2021, reported diagnoses of influenza A and B were way down.47 The relevance of studies of influenza in public transport is, at the very least, questionable.

The evidence presented here also highlights the methodological limitations of alternative approaches to studying the role of the subways in the propagation of SARS-CoV-2. The test conducted in Figs. 2b–2e demonstrates the importance of studying changes in subway volume during the course of the COVID-19 outbreak. Less informative would be a study relating COVID-19 rates to static survey data on the proportion of individuals in each ZCTA regularly riding public transit prior to the epidemic. Our results also point to the importance of conducting tests of causation when baseline subway volume and COVID-19 incidence are high. A finding that coronavirus cases no longer relate to subway volume once subway use has plummeted to below 10 percent of baseline reveals little if anything about what happened back in March. The map of the Flushing Local line in Fig. 2b further highlights the pitfalls of studies that assign the entire volume of turnstile entries into a subway station to its enclosing ZCTA.7,8 Such a procedure, which effectively assumes that only people who live in the same ZCTA take the local subway, would erroneously discard the high-incidence ZCTAs 11369 and 11370, which have no subway within their boundaries.

If we are to successfully control future pandemic threats – and, for that matter, future outbreaks of COVID-19 – we need to understand in exhaustive detail how SARS-CoV-2 first took hold and then established hot spots in major urban epicenters throughout the world. Considerable effort has been made to understand exactly what happened in Wuhan.29,48 A study of Los Angeles County has tracked the initial seeding of imported infections in affluent areas as it spread radially to high-density neighborhoods, where the virus percolated through multi-generational households.14 While the outbreak in Italy has been traced phylogenetically to the Lombardy region, it remains unclear how exactly it started out and spread.49 Numerous investigators have relied upon compartmental models to understand the early dynamics of SARS-CoV-2 outbreaks.14,29,50,51 The evidence presented here for New York City points instead to a model of network-wide transmission followed by local percolation of infections.52-55

If the subway system indeed played a critical role in the early propagation of SARS-CoV-2, as supported by the evidence assembled here, we need to understand that the conventional methods of personal contact tracing are less likely to be useful in halting future outbreaks. That means more sophisticated contact tracing through the pings of mobile devices and records of electronic transactions will be necessary.56,57 To that end, the MTA will need to adopt a new system of digital passes, already in use in many cities worldwide, which would permit investigators to find out more than just the crude number of turnstile-clicks at each station.

In advance of the next outbreak, we will need to know whether the subways served principally as a rapid spatial disseminator of externally acquired infections, or as significant locus of in situ transmission. In the former case, social distancing and mandatory face coverings would not alone stop the rapid, widespread seeding of infections throughout the five boroughs that we observed in February and March of 2020. In the latter case, we will need to study now whether a policy of running only express lines with limited density might be a feasible alternative to the complete cordon sanitaire adopted in Wuhan more than a year ago.58

Methods

Data on Confirmed COVID-19 Cases

The NYC health department’s open data archive20 was our source of data on: confirmed COVID-19 cases and hospitalizations by borough and date of diagnosis (boroughs-case-hosp-death, used to construct Figs. 1a, 1b and 2b), aggregate, city-wide data on cases and hospitalizations by date of diagnosis (case-hosp-death, used in part to construct Fig. 2a); and cumulative cases by zip code tabulation area (ZCTA) (tests-by-zcta, used in part to construct Figs. 2c and 3d). Incidence per 10,000 population was based population counts described below.

Population Data

Data on the total populations of zip code tabulation areas (ZCTAs) were derived from the Census Bureau’s American Community Survey 5-year estimates for 2015–2019, accessed from the data server at the Missouri Census Data Center.59 Data on the total populations of census block groups (CBGs) were likewise derived from the Census Bureau’s American Community Survey 5-year estimates for 2015–2019, accessed from the Census Bureau’s website.60

Geography

The Metropolitan Transportation Authority (MTA) website for developers 61 was our source for the geocoordinates (longitude and latitude) of each of the subway stations, including the 22 stations on the Flushing Local (Number 7) line, as depicted in Figs. 1d, 2d, 3b, and 4b.

We downloaded the polygon shapes of all census block groups (CBGs) in New York City from the Census Bureau’s website.62 We relied on the Stata program geoinpoly,63 which uses a ray-casting algorithm to determine whether a point is contained in a polygon, to identify the unique CBG containing each subway station (as illustrated by the 82nd St–Jackson Heights station in Figs. 2d and 3a).

To map CBGs into ZCTAs, we proceeded in four steps. First, we used Stata mapping software to verify that most CBGs were uniquely contained in a given ZCTA (Supplement Fig. A). Second, we employed QGIS software to compute centroids of each CBG in New York City based upon the Census Bureau’s polygon shape files. Third, we downloaded the polygon shape files of all ZCTAs from the New York City health department’s data archive.64 Finally, we employed geoinpoly once again to determine the ZCTA shape polygon that contained the centroid of each CBG.

As discussed in detail below, our analysis of the prevalence of at-risk multi-generational households relied upon the Census Bureau’s American Community Survey Public Use Microdata Sample (PUMS) for the 5-year period 2015–2019 65. The data records for the PUMS are identified at the level of the Public Use Microdata Area (PUMA) 66, which is an aggregate of census tracts, which are in turn aggregates of CBGs. To map PUMAs into ZCTAs, we downloaded the scheme for aggregating New York City census tracts into PUMAs from the ESRI’s ArcGIS Hub 67, which in turn gave us a mapping from PUMAs to CBGs. We then relied on our prior mapping of CBGs into ZCTAs to go from PUMAs directly to CBGs, as seen in Supplement Fig. A.

Data on Phylogenetic Analysis of Viral Isolates

To construct Fig. 1c, we relied upon two data sources: (a) the tab entitled Clade A2a GISAID IDs within in the spreadsheet Data File S2, posted in the Supplementary Materials of Gonzalez-Reiche et al.;15 and (b) the spreadsheet Supplementary Table 2, posted in the Supplementary Materials of a later study of COVID-19 patients treated within the New York University Langone Hospital system.68 We merged the two files on the unique common identifier variable gisaid_epi_isi (where GISAID stands for Global Initiative on Sharing All Influenza Data). This gave a total of 78 MSHS viral samples authored by Gonzalez-Reiche et al. within the A2a clade, including date, location and strain identifier. These 78 samples formed the database for the vertical bars in the figure.

Next, we used the variable strain in the merged file to identify the 17 virus samples specifically highlighted as sharing the ORF1b:A1844V mutation in the New York Cluster 1 in Figure 2C of Gonzalez-Reiche et al.15 These 17 samples are indicated as the pink bubbles in Fig. 1c. This mutation resulted from a single amino acid substitution from alanine (A) to valine (V) at position #1844 in the stretch of the virus’ RNA coding for its ORF1b protein, which is one of the two replicase proteins common to SARS coronaviruses. In terms of the virus’ underlying genetic code, the mutation corresponded to a single base substitution in the virus’ positive-sense mRNA codon from GUX to GCX, where G = guanine, U = uracil, C = cytosine, A = adenine, and X = any of these four bases. This single RNA base substitution (or missense mutation) was shared by samples of infected persons residing in Manhattan, Queens, Brooklyn and Westchester County, collected during the space of only five days (March 14–18).

Data on Subway Turnstile Entries

The data on turnstile entries were similarly derived from the MTA’s website for developers.21 Since stations typically have multiple turnstiles, and since the turnstile counters are updated at intervals during each day, computation of entries by station and by date involved the aggregation of data points across large data sets with millions of individual observations. Accurate coding required us to take account of the fact that some turnstiles ran backwards, while others were reset when they reached their numerical limit. Still, the city-wide temporal patterns seen in Fig. 2a are consistent with other independent estimates.69

Classification of Flushing Line (Local 7) Stations

In Figs. 2d and 2e, we classified subway stations along the 7 (Flushing) Line into three groups: The six key stations within the Queens-Elmhurst hot spot, indicated in yellow from west to east, were: 74th St – Broadway; 82nd St – Jackson Hts; 90th St – Elmhurst Av; Junction Blvd; 103rd St – Corona Plaza; and 111 St. The stations within Manhattan, indicated in sky blue from west to east, were: 34th St – Hudson Yards, Times Sq – 42nd St, 5th Ave – Bryant Pk, and Grand Central – 42nd St. The remaining stations within the borough of Queens are indicated in pink.

Data on Smartphone Device Movements

Our data on smartphone device movements come from the Social Distancing database maintained by SafeGraph.32 Every device movement (or visit) had an origin and a destination. Each device’s unique origin was the CBG where it regularly spent the night. Every CBG in which the device stopped for more than 1 minute during a 24-hour period was the destination of a visit. For each calendar day and each CBG of origin, the database recorded the number of devices that visited each destination CBG. A destination CBG can be the same as the origin CBG.

We tested whether smartphone device movements whose destination CBG contained a subway station could serve as a proxy for subway turnstile entries. For each station, we compared two time series: the number of visits to the destination CBG containing that subway station, which we’ll call subway CBG visits, and the number of turnstile entries into that station. This comparison is illustrated for a particular subway station in Fig. 3a.

We further investigated the origins of those smartphone devices whose destination CBGs contained one of the six key stations within the Queens-Elmhurst hot spot. For each CBG, we determined two visit counts. The first count, which we denote n1, accumulated the total number of visits originating in that CBG with a destination at any one of the six stations during the months of January and February 2020. (The 74th Street–Broadway station on the 7 (Flushing) line shared the same CBG as the Jackson Heights–Roosevelt Ave. station on the intersecting 6th Avenue Local (M) line.) The second count, which we denote n2, accumulated the total number of visits originating in the same CBG during the same interval with a destination at either the Queensboro Plaza or Court Square stops, two of the principal destinations within the Queens portion of the 7 (Flushing) line. We then ranked each origin CBG by the statistic nmin = min{n1, n2}, which captured trips to and from the Queens-Elmhurst yellow stations and the Queensboro Plaza–Court Square complex. In Fig. 3b, the lighter-shaded CBGs correspond to 100 > nmin ≥ 50, while the darker shaded CBGs correspond to nmin ≥ 100.

To estimate subway visits by ZCTA, we aggregated the number of device visits to all destination CBGs containing a subway station, and then further aggregated these CBG-specific counts of subway visits at the ZCTA level. Supplement Fig. A illustrates the congruence between CBGs and ZCTAs. Supplement Fig. B illustrates the temporal evolution of visits to all subway station CBGs originating from four specific ZCTAs: 10003 (Manhattan), 11201 (Brooklyn), 11205 (Brooklyn), and 11368 (Queens).

Our reconstruction of the origins of subway visits from smartphone mobility data is to be distinguished from prior studies relying instead upon the SafeGraph Patterns Schema, a separate database which classifies visits by their destination points of interest.51,70 The latter database did not categorize subway stations as a point of interest

Prevalence of At-Risk Multi-Generational Households

We relied upon the five-year (2015–2019) public use microsample of the U.S. Census Bureau’s American Community (ACS)65 to estimate the proportion of households in New York City that were at risk for multi-generational transmission of SARS-CoV-2. Following an earlier study of intra-household transmission in Los Angeles County,14 we defined an at-risk household as having at least 4 persons, of whom at least one person was 18–34 years of age and at least one other person was at least 50 years of age. Based upon a subsample of 148,686 New York City households in the 5-year ACS database, we found that 18.3 percent of households satisfied this criterion. Across 55 public use microdata areas (PUMAs), the median proportion of at-risk households was 22.0 percent, with the 25th and 75th percentiles at 15.6 and 25.4 percent, respectively. As described above, we then mapped the PUMA-specific estimates into ZCTAs. Across 176 ZCTAs, the median proportion of households at risk was 22.4 percent, with the 25th and 75th percentiles at 13.7 and 24.8 percent, respectively. The minimum proportion was 3.2 percent (ZCTA 10017 in Manhattan), while the maximum proportion was 35.8 percent (11414 and 11420 in Queens).

Contiguity in Geographic and Subway Space

Our concepts of geographic and subway contiguity, including an accompanying formal matrix algebra, are developed in detail in the Supplement. Briefly, the map of ZCTAs in New York City can be regarded as a finite set of M > 0 compact polygons in a two-dimensional plane, indexed by i = 1, …, M. No two ZCTAs share any interior points in common, but they can share boundary points. When ZCTAs i and j do share boundary points, we say that they are geographically contiguous, or g-contiguous. By contrast, when ZCTA j is the next stop after ZCTA i on some subway line in some direction, we say that ZCTAs i and j are contiguous in subway space, or s-contiguous. G-contiguity does not imply s-contiguity, nor does s-contiguity imply g-contiguity.

As further elaborated in detail in the Supplement, we formulated compound relationships based on the elemental notions of g- and s-contiguity. Fig. 4a shows all ZCTAs that are (g+g2)-contiguous with ZCTA 11415. Equivalently, the figure displays all ZCTAs within a geographic contiguity radius of 2. Fig. 4b displays all ZCTAs that are (g+s+s2+s3+s4+s5)-contiguous with ZCTA 11415, that is all ZCTAs that are either g-contiguous with that ZCTA or within a subway radius of 5 stops along the same or a connecting line.

Non-Spatial Regressions

Let y denote a M × 1 column vector of ZCTA-specific observations of incremental COVID-19 incidence during April 1–8, 2020 (mapped in Fig. 3d). Let X0 denote the corresponding ZCTA-specific column vector of observations on the cumulative incidence of COVID-19 as of March 31 (Fig. 2c). Let X1 denote the corresponding vector of observations on relative subway volume as of March 16, 2020 (Fig. 3c), and let X2 denote the prevalence of at-risk multigenerational households (Fig. 3f). As detailed in the Supplement, we estimated non-spatial models of the form log y = α + β0 log X0 + β1 log X1 + β2 log X2, where the logarithm is assumed to operate separately on each vector coordinate.

Spatial Regressions

We then considered spatial regression models of the form log y = α + β0 log X0 + β1 log X1 + β2 log X2 + γ0 log WX0 + γ1 log WX1 + γ2 log WX2, where W is an M × M spatial weighting matrix. Each contiguity criterion necessarily had its own weighting matrix W. As detailed in the Supplement, pre-multiplication of each vector X0, X1, and X2 by W computed its respective population-weighted mean value among all ZCTAs satisfying the particular contiguity criterion.

Data Availability

This study relies on several large databases that are publicly accessible via the Internet. The Reference section contains the URLs for access to the underlying data on COVID-19 cases and hospitalizations, subway turnstile entries, and phylogenetic analyses of virus samples. The underlying data on smartphone movements, provided by SafeGraph, are available to researchers only by permission. We have posted our data analyses at the Open Science Framework (OSF) in a project entitled New York City COVID-19 Epidemic.

https://osf.io/v7k23/

Competing Interests Declaration

The author has no competing interests to declare.

Funding Declaration

The author has no funding sources to declare.

Human Subjects Declaration

This study relies exclusively on anonymized, publicly available data that contain no individual identifiers.

Data Availability Statement

This study relies on several large databases that are publicly accessible via the Internet. The Reference section contains the URLs for access to the underlying data on COVID-19 cases and hospitalizations, subway turnstile entries, and phylogenetic analyses of virus samples. The underlying data on smartphone movements, provided by SafeGraph, are available to researchers only by permission. We have posted our data analyses at the Open Science Framework (OSF) in a project entitled New York City COVID-19 Epidemic (https://osf.io/v7k23/).

Acknowledgments

This article represents the sole opinion of its author and does not necessarily represent the opinions of the Massachusetts Institute of Technology, Eisner Health, or any other organization or individual.

References

  1. 1.↵
    Metropolitan Transportation Authority. About Us. The MTA Network: Public Transportation for the New York Region. http://web.mta.info/mta/network.htm: Last accessed December 27, 2020 2018.
  2. 2.↵
    Metropolitan Transportation Authority. Subway and bus ridership for 2019. https://new.mta.info/agency/new-york-city-transit/subway-bus-ridership-2019: Updated April 14, 2020 2020.
  3. 3.↵
    Washington Metropolitan Area Transit Authority. Rail Ridership Data Viewer. https://www.wmata.com/initiatives/ridership-portal/Rail-Data-Portal.cfm: Last accessed December 27 2020.
  4. 4.↵
    Hu M, Lin H, Wang J, et al. The risk of COVID-19 transmission in train passengers: an epidemiological and modelling study. Clin Infect Dis 2020 doi: 10.1093/cid/ciaa1057 [published Online First: 2020/07/30]
    OpenUrlCrossRef
  5. 5.
    Du Z, Wang L, Cauchemez S, et al. Risk for Transportation of Coronavirus Disease from Wuhan to Other Cities in China. Emerg Infect Dis 2020;26(5):1049–52. doi: 10.3201/eid2605.200146 [published Online First: 2020/02/14]
    OpenUrlCrossRefPubMed
  6. 6.↵
    Zheng R, Xu Y, Wang W, et al. Spatial transmission of COVID-19 via public and private transportation in China. Travel Med Infect Dis 2020;34:101626. doi: 10.1016/j.tmaid.2020.101626 [published Online First: 2020/03/19]
    OpenUrlCrossRef
  7. 7.↵
    Sy KTL, Martinez ME, Rader B, et al. Socioeconomic Disparities in Subway Use and COVID-19 Outcomes in New York City. Am J Epidemiol 2020 doi: 10.1093/aje/kwaa277 [published Online First: 2020/12/30]
    OpenUrlCrossRef
  8. 8.↵
    Glaeser EL, Gorback C, Redding SJ. JUE insight: How much does COVID-19 increase with mobility? Evidence from New York and four other U.S. cities. https://www.sciencedirect.com/science/article/abs/pii/S0094119020300632?via%3Dihub: Journal of Urban Economics, ePub October 21 2020.
  9. 9.↵
    Fathi-Kazerooni S, Rojas-Cessa R, Dong Z, et al. Correlation of subway turnstile entries and COVID-19 incidence and deaths in New York City. Infect Dis Model 2021;6:183–94. doi: 10.1016/j.idm.2020.11.006 [published Online First: 2020/12/10]
    OpenUrlCrossRef
  10. 10.↵
    Centers for Disease Control and Prevention. Updated Guidance on Evaluating and Testing Persons for Coronavirus Disease 2019 (COVID-19). https://emergency.cdc.gov/han/2020/han00429.asp: Distributed via the CDC Health Alert Network, March 08, 2020, 8:20 PM ET, CDCHAN-00429 2020.
  11. 11.↵
    Lauer SA, Grantz KH, Bi Q, et al. The Incubation Period of Coronavirus Disease 2019 (COVID-19) From Publicly Reported Confirmed Cases: Estimation and Application. Ann Intern Med 2020;172(9):577–82. doi: 10.7326/M20-0504 [published Online First: 2020/03/10]
    OpenUrlCrossRefPubMed
  12. 12.↵
    Qin J, You C, Lin Q, et al. Estimation of incubation period distribution of COVID-19 using disease onset forward time: A novel cross-sectional and forward follow-up study. Sci Adv 2020;6(33):eabc1202. doi: 10.1126/sciadv.abc1202 [published Online First: 2020/08/28]
    OpenUrlFREE Full Text
  13. 13.↵
    Cohen PA, Hall LE, John JN, et al. The Early Natural History of SARS-CoV-2 Infection: Clinical Observations From an Urban, Ambulatory COVID-19 Clinic. Mayo Clin Proc 2020;95(6):1124–26. doi: 10.1016/j.mayocp.2020.04.010 [published Online First: 2020/05/27]
    OpenUrlCrossRefPubMed
  14. 14.↵
    Harris JE. Los Angeles County SARS-CoV-2 Epidemic: Critical Role of Multi-generational Intra-household Transmission. Journal of Bioeconomics 2021;23:55–83. doi: https://doi.org/10.1007/s10818-021-09310-2
    OpenUrl
  15. 15.↵
    Gonzalez-Reiche AS, Hernandez MM, Sullivan MJ, et al. Introductions and early spread of SARS-CoV-2 in the New York City area. Science 2020;369(6501):297–301. doi: 10.1126/science.abc1917 [published Online First: 2020/05/31]
    OpenUrlAbstract/FREE Full Text
  16. 16.↵
    Lemieux JE, Siddle KJ, Shaw BM, et al. Phylogenetic analysis of SARS-CoV-2 in Boston highlights the impact of superspreading events. Science 2020 doi: 10.1126/science.abe3261 [published Online First: 2020/12/12]
    OpenUrlAbstract/FREE Full Text
  17. 17.↵
    Metropolitan Transportation Authority (MTA). Letter from Chairman Patrick J. Foye to Senator Charles E. Schumer and Congresswoman Nita Lowey. April 16, 20202020.
  18. 18.↵
    Martinez J. NYC Subway Crews Hit Hardest by Coronavirus, MTA Numbers Show. https://www.thecity.nyc/2020/6/1/21277407/nyc-subway-crews-hit-hardest-by-coronavirus-pandemic-mta-numbers-show: The City, June 1 2020.
  19. 19.↵
    Metropolitan Transportation Authority (MTA). Queens Bus Map. http://web.mta.info/nyct/maps/busqns.pdf: Last Visited April 12, 2020 2018.
  20. 20.↵
    New York Department of Health and Mental Hygiene. NYC Coronavirus Disease 2019 (COVID-19) Data. https://github.com/nychealth/coronavirus-data: Last accessed, Jan. 2, 2021 2020.
  21. 21.↵
    Metropolitan Transportation Authority (MTA). Turnstile Data. http://web.mta.info/developers/turnstile.html: Last accessed April 4, 2020 2020.
  22. 22.↵
    Oran DP, Topol EJ. Prevalence of Asymptomatic SARS-CoV-2 Infection : A Narrative Review. Ann Intern Med 2020;173(5):362–67. doi: 10.7326/M20-3012 [published Online First: 2020/06/04]
    OpenUrlCrossRefPubMed
  23. 23.↵
    Havers FP, Reed C, Lim TW, et al. Seroprevalence of Antibodies to SARS-CoV-2 in Six Sites in the United States, March 23-May 3, 2020. https://www.medrxiv.org/content/10.1101/2020.06.25.20140384v1: MedRxiv June 26, 2020 2020.
  24. 24.↵
    Office of the Mayor. Statement From Mayor de Blasio on Bars, Restaurants, and Entertainment Venues. https://www1.nyc.gov/office-of-the-mayor/news/152-20/statement-mayor-de-blasio-bars-restaurants-entertainment-venues: City of New York, March 15 2020.
  25. 25.↵
    Vielkind J, Grayce M. Amid Coronavirus Fears, Officials Say New York City Subway Is Safe. https://www.wsj.com/amp/articles/amid-coronavirus-fears-officials-say-new-york-city-subway-is-safe-11583272664: Wall Street Journal, March 3 2020.
  26. 26.↵
    Parnell W, Shahrigian S. Mayor de Blasio says coronavirus fears shouldn’t keep New Yorkers off subways. https://www.nydailynews.com/coronavirus/ny-coronavirus-bill-de-blasio-coronavirus-subway-20200305-vmjdxjudbndlrjekashqs3hfou-story.html: New York Daily News, March 5 2020.
  27. 27.↵
    Metropolitan Transportation Authority. MTA New York City Subway: Large Print Edition with Railway Connections. https://new.mta.info/sites/default/files/2019-03/large_print_map_2019.pdf: Published March 2019. Last Accessed, January 2, 2021 2019.
  28. 28.↵
    Griffin JM, Collins AB, Hunt K, et al. A rapid review of available evidence on the serial interval and generation time of COVID-19. https://www.medrxiv.org/content/10.1101/2020.05.08.20095075v1: MedRxiv May 11, 2020 2020.
  29. 29.↵
    Hao X, Cheng S, Wu D, et al. Reconstruction of the full transmission dynamics of COVID-19 in Wuhan. Nature 2020;584(7821):420–24. doi: 10.1038/s41586-020-2554-8 [published Online First: 2020/07/17]
    OpenUrlCrossRefPubMed
  30. 30.↵
    Xinhua Net. China’s Wuhan reopens subway, railway station http://www.xinhuanet.com/english/2020-03/28/c_138926565.htm: March 28 2020.
  31. 31.↵
    D’Arienzo M, Coniglio A. Assessment of the SARS-CoV-2 basic reproduction number, R 0, based on the early phase of COVID-19 outbreak in Italy. Biosaf Health 2020;2(2):57–59. doi: 10.1016/j.bsheal.2020.03.004 [published Online First: 2020/08/25]
    OpenUrlCrossRefPubMed
  32. 32.↵
    SafeGraph Inc. Social Distancing Metrics. https://docs.safegraph.com/docs/social-distancing-metrics: Accessed September 24, 2020 2020.
  33. 33.↵
    Harris JE. Los Angeles County SARS-CoV-2 Epidemic: Critical Role of Multi-generational Intra-household Transmission. https://www.medrxiv.org/content/10.1101/2020.10.11.20211045v3: medRxIv, updated February 9, in press, Journal of Bioeconomics 2021.
  34. 34.
    Esteve A, Permanyer I, Boertien D, et al. National age and coresidence patterns shape COVID-19 vulnerability. Proc Natl Acad Sci U S A 2020;117(28):16118–20. doi: 10.1073/pnas.2008764117 [published Online First: 2020/06/25]
    OpenUrlAbstract/FREE Full Text
  35. 35.↵
    Aparicio Fenoll A, Grossbard S. Intergenerational residence patterns and Covid-19 fatalities in the EU and the US. Econ Hum Biol 2020;39:100934. doi: 10.1016/j.ehb.2020.100934 [published Online First: 2020/11/08]
    OpenUrlCrossRef
  36. 36.↵
    Franch-Pardo I, Napoletano BM, Rosete-Verges F, et al. Spatial analysis and GIS in the study of COVID-19. A review. Sci Total Environ 2020;739:140033. doi: 10.1016/j.scitotenv.2020.140033 [published Online First: 2020/06/14]
    OpenUrlCrossRefPubMed
  37. 37.↵
    World Health Organization. Transmission of SARS-CoV-2: implications for infection prevention precautions. https://www.who.int/news-room/commentaries/detail/transmission-of-sars-cov-2-implications-for-infection-prevention-precautions: July 9 2020.
  38. 38.↵
    Santarpia JL, Rivera DN, Herrera VL, et al. Aerosol and surface contamination of SARS-CoV-2 observed in quarantine and isolation care. Sci Rep 2020;10(1):12732. doi: 10.1038/s41598-020-69286-3 [published Online First: 2020/07/31]
    OpenUrlCrossRefPubMed
  39. 39.↵
    Azuma K, Yanagi U, Kagi N, et al. Environmental factors involved in SARS-CoV-2 transmission: effect and role of indoor environmental quality in the strategy for COVID-19 infection control. Environ Health Prev Med 2020;25(1):66. doi: 10.1186/s12199-020-00904-2 [published Online First: 2020/11/05]
    OpenUrlCrossRef
  40. 40.↵
    Centers for Disease Control and Prevention. Wear Face Masks on Public Transportation Conveyances and at Transportation Hubs. https://www.cdc.gov/coronavirus/2019-ncov/travelers/face-masks-public-transportation.html: November 10 2020.
  41. 41.↵
    Qian H, Miao T, Liu L, et al. Indoor transmission of SARS-CoV-2. Indoor Air 2020 doi: 10.1111/ina.12766 [published Online First: 2020/11/02]
    OpenUrlCrossRef
  42. 42.↵
    Gosce L, Barton DA, Johansson A. Analytical modelling of the spread of disease in confined and crowded spaces. Sci Rep 2014;4:4856. doi: 10.1038/srep04856 [published Online First: 2014/05/07]
    OpenUrlCrossRef
  43. 43.↵
    Gosce L, Johansson A. Analysing the link between public transport use and airborne transmission: mobility and contagion in the London underground. Environ Health 2018;17(1):84. doi: 10.1186/s12940-018-0427-5 [published Online First: 2018/12/06]
    OpenUrlCrossRef
  44. 44.↵
    Cooley P, Brown S, Cajka J, et al. The role of subway travel in an influenza epidemic: a New York City simulation. J Urban Health 2011;88(5):982–95. doi: 10.1007/s11524-011-9603-4 [published Online First: 2011/08/10]
    OpenUrlCrossRefPubMed
  45. 45.↵
    Howland RE, Cowan NR, Wang SS, et al. Public transportation and transmission of viral respiratory disease: Evidence from influenza deaths in 121 cities in the United States. PLoS One 2020;15(12):e0242990. doi: 10.1371/journal.pone.0242990 [published Online First: 2020/12/02]
    OpenUrlCrossRef
  46. 46.↵
    Biggerstaff M, Cauchemez S, Reed C, et al. Estimates of the reproduction number for seasonal, pandemic, and zoonotic influenza: a systematic review of the literature. BMC Infect Dis 2014;14:480. doi: 10.1186/1471-2334-14-480 [published Online First: 2014/09/05]
    OpenUrlCrossRefPubMed
  47. 47.↵
    Centers for Disease Control and Prevention. FluView Summary ending on January 23, 2021. https://www.cdc.gov/flu/weekly/weeklyarchives2020-2021/week03.htm2021.
  48. 48.↵
    Pekar J, Worobey M, Moshiri N, et al. Timing the SARS-CoV-2 index case in Hubei province. Science 2021;372(6540):412–17. doi: 10.1126/science.abf8003 [published Online First: 2021/03/20]
    OpenUrlAbstract/FREE Full Text
  49. 49.↵
    Micheli V, Rimoldi SG, Romeri F, et al. Geographical reconstruction of the SARS-CoV-2 outbreak in Lombardy (Italy) during the early phase. J Med Virol 2021;93(3):1752–57. doi: 10.1002/jmv.26447 [published Online First: 2020/08/21]
    OpenUrlCrossRef
  50. 50.↵
    Harris JE. Data from the COVID-19 epidemic in Florida suggest that younger cohorts have been transmitting their infections to less socially mobile older adults. Review of Economics of the Household 2020(August 11, 2020) doi: doi.org/10.1007/s11150-020-09496-w
    OpenUrlCrossRef
  51. 51.↵
    Chang S, Pierson E, Koh PW, et al. Mobility network models of COVID-19 explain inequities and inform reopening. Nature 2020 doi: 10.1038/s41586-020-2923-3 [published Online First: 2020/11/11]
    OpenUrlCrossRefPubMed
  52. 52.↵
    Sander LM, Warren CP, Sokolov IM, et al. Percolation on heterogeneous networks as a model for epidemics. Math Biosci 2002;180:293–305. doi: 10.1016/s0025-5564(02)00117-7 [published Online First: 2002/10/22]
    OpenUrlCrossRefPubMedWeb of Science
  53. 53.
    Keeling MJ, Eames KT. Networks and epidemic models. J R Soc Interface 2005;2(4):295–307. doi: 10.1098/rsif.2005.0051 [published Online First: 2006/07/20]
    OpenUrlCrossRefPubMedWeb of Science
  54. 54.
    Kenah E, Miller JC. Epidemic percolation networks, epidemic outcomes, and interventions. Interdiscip Perspect Infect Dis 2011;2011:543520. doi: 10.1155/2011/543520 [published Online First: 2011/03/26]
    OpenUrlCrossRefPubMed
  55. 55.↵
    Stegehuis C, van der Hofstad R, van Leeuwaarden JS. Epidemic spreading on complex networks with community structures. Sci Rep 2016;6:29748. doi: 10.1038/srep29748 [published Online First: 2016/07/22]
    OpenUrlCrossRef
  56. 56.↵
    Kang CR, Lee JY, Park Y, et al. Coronavirus Disease Exposure and Spread from Nightclubs, South Korea. Emerg Infect Dis 2020;26(10):2499–501. doi: 10.3201/eid2610.202573 [published Online First: 2020/07/08]
    OpenUrlCrossRef
  57. 57.↵
    Harris JE. Geospatial Analysis of the September 2020 Coronavirus Outbreak at the University of Wisconsin – Madison: Did a Cluster of Local Bars Play a Critical Role? https://www.nber.org/papers/w28132: National Bureau of Economic Research Working Paper No. 28132, November 18 2020.
  58. 58.↵
    BBC News. Coronavirus: Wuhan shuts public transport over outbreak. https://www.bbc.com/news/world-asia-china-51215348: January 23 2020.
  59. 59.↵
    Missouri Census Data Center. ACS Profiles Menu. https://mcdc.missouri.edu/applications/acs/profiles/: Accessed Jan 1, 2021 2020.
  60. 60.↵
    U.S Census Bureau. ACS Summary File Data. https://www2.census.gov/programs-surveys/acs/summary_file/2019/data/: November 28 2020.
  61. 61.↵
    Metropolitan Transportation Authority (MTA). Stations.csv (text file, comma-separated). http://web.mta.info/developers/data/nyct/subway/Stations.csv: Last accessed January 2, 2021 2020.
  62. 62.↵
    U.S Census Bureau. TIGER2017 Shape Files for Census Block Groups by State. https://www2.census.gov/geo/tiger/TIGER2017/BG/: September 27 2017.
  63. 63.↵
    Picard R. GEOINPOLY: Stata module to match geographic locations to shapefile polygons. https://ideas.repec.org/c/boc/bocode/s458016.html: August 16 2015.
  64. 64.↵
    New York City Department of Health and Mental Hygiene. nychealth/coronavirus-data/Geography-resources/. https://github.com/nychealth/coronavirus-data/tree/master/Geography-resources: Updated December 13 2020.
  65. 65.↵
    U.S Census Bureau. Accessing PUMS Data. https://www.census.gov/programs-surveys/acs/microdata/access.html: January 14 2021.
  66. 66.↵
    U.S Census Bureau. Public Use Microdata Areas (PUMAs). https://www.census.gov/programs-surveys/geography/guidance/geo-areas/pumas.html: August 18, 2020 2020.
  67. 67.↵
  68. 68.↵
    Maurano MT, Ramaswami S, Zappile P, et al. Sequencing identifies multiple early introductions of SARS-CoV-2 to the New York City Region. https://www.medrxiv.org/content/10.1101/2020.04.15.20064931v4: medRxiv August 19 2020.
  69. 69.↵
    Schneider TW. New York City Subway Usage. https://toddwschneider.com/dashboards/nyc-subway-turnstiles/: Last visited April 4, 2020 2020.
  70. 70.↵
    SafeGraph Inc. Places Schema. https://docs.safegraph.com/docs/places-schema: Data accessed July 30-31, September 24-26, 2020 2020.
Back to top
Posted July 07, 2021.
Download PDF

Supplementary Material

Data/Code
Email
Critical Role of the Subways in the Initial Spread of SARS-CoV-2 in New York City
Jeffrey E. Harris
medRxiv 2021.07.03.21259973; doi: https://doi.org/10.1101/2021.07.03.21259973
Twitter logo Facebook logo LinkedIn logo Mendeley logo
Citation Tools

Subject Area