Compound risks of hurricane evacuation amid the COVID-19 pandemic in the United States

Sen Pei, Kristina A. Dahl, Teresa K. Yamana, Rachel Licker, Jeffrey Shaman
doi: https://doi.org/10.1101/2020.08.07.20170555

Abstract

Current projections and unprecedented storm activity to date suggest the 2020 Atlantic hurricane season will be extremely active and that a major hurricane could make landfall during the global COVID-19 pandemic. Such an event would necessitate a large-scale evacuation, with implications for the trajectory of the pandemic. Here we model how a hypothetical hurricane evacuation from four counties in southeast Florida would affect COVID-19 case levels. We find that hurricane evacuation increases the total number of COVID-19 cases in both origin and destination locations; however, if transmission rates in destination counties can be kept from rising during evacuation, excess evacuation-induced case numbers can be minimized by directing evacuees to counties experiencing lower COVID-19 transmission rates. Ultimately, the number of excess COVID-19 cases produced by the evacuation depends on the ability of destination counties to meet evacuee needs while minimizing virus exposure through public health directives.

Introduction

The combination of the COVID-19 pandemic, existing racial and socioeconomic inequalities, and environmental stressors exacerbated by climate change is exposing the many ways in which “compound risks” threaten human lives and wellbeing while straining the ability of governments at all scales to limit the damage from any one threat on its own1. Intersections of climate extremes with the pandemic—recent widespread flooding in South Asia at a time of rapidly increasing COVID-19 caseloads, for example— have made clear that the consequences of such compound risk events can be lethal, though the underreporting of cases around the world2 and widely varying testing capabilities3 make it difficult to accurately quantify their magnitude.

With the anticipated peak of the 2020 Atlantic hurricane season approaching and COVID-19 cases widespread and abundant in many hurricane-prone areas of the United States, the nation is poised to experience the collision of two major disasters. This study therefore addresses how decision-making around one key aspect of hurricane response—evacuation—could influence the trajectory of the pandemic in the US and be optimized to limit excess COVID-19 cases. With future global warming expected to continue the observed trend toward increasingly intense Atlantic hurricanes4,5, understanding how to manage and minimize the impact of the combined risks associated with a major hurricane and a global pandemic could prove critical both later this year and in the long-term as the risks of such simultaneous disasters increase around the world.

Efficient, effective evacuations—whether voluntary or mandatory—are a critical component of ensuring public safety during natural disasters. The scale of recent evacuations from US Southeast and Gulf Coast states has been large: During Hurricanes Matthew (2016), Irma (2017), and Dorian (2019), for example, roughly 2.5, 6.5, and 1.1 million people, respectively, were under evacuation orders 6-8. By changing the distribution of people for days or weeks, a large-scale evacuation amid a pandemic has the potential to alter the trajectory and geographic distribution of infections. And by temporarily relocating from their own homes into potentially shared living arrangements where levels of social contact are higher, evacuees may experience greater transmission risk both during and after an evacuation.

This analysis evaluates how a large-scale evacuation of the southeast Florida coast from a hypothetical Category 3 hurricane would affect the total number of COVID-19 cases and their spatial distribution in evacuees’ origin and destination counties. This is accomplished by first building a hypothetical hurricane evacuation scenario from previously published studies of evacuation behavior in the southeast US region. We then use a simple, two-county infectious disease model to identify the most relevant evacuation and epidemiological characteristics influencing COVID-19 case counts. The findings from the these simulations are used to inform experiments with a larger metapopulation model representing SARS-CoV-2 transmission in all 3,142 US counties9, as well as various evacuation scenarios.

RESULTS

Identifying key parameters using a two-county model

We used a simplified, two-county metapopulation model, representing a generic pair of origin and destination counties, to determine the factors that have greatest influence on COVID-19 case numbers (Methods; Figure 1a). Specifically, we evaluated the effects of six evacuation and epidemiological characteristics on COVID-19 case numbers: transmission rates in origin and destination counties (quantified by the effective reproductive number, Re), the fraction of the origin county population that evacuates (peva), the duration of the evacuation period (Teva), and daily case numbers in the origin and destination counties (caseori and casedest). We simulated an evacuation by moving a fraction of the population from the origin to the destination county. Evacuees then mixed with the population of the destination county, before returning home. This simulation was repeated using different combinations of each of the six characteristics in order to determine the effects of each on COVID-19 case numbers during and following the evacuation.

Figure 1:
Figure 1:

Results from the two-county model showing that origin and destination transmission rates have the greatest influence on final case numbers. (a) A schematic diagram for the two-county model. Blue and orange boxes represent the origin and destination populations. Red dots within boxes represent infected individuals. (b) The marginal distribution of six parameters for the top 10% of combinations that lead to the lowest percentage increase (or highest percentage reduction) of reported cases in the origin county (solid red lines), the destination county (solid orange lines) and both counties combined (solid blue lines). Here Embedded Image and Embedded Image represent the transmission rates in the origin and destination; caseori and casedest represent the daily cases in the origin and destination; Teva is the duration of evacuation; and peva is the fraction of the origin population evacuating.

We found that transmission rates in the origin and destination counties were the primary determinant of case numbers (Methods; Figure 1b): evacuating individuals from a high-Re origin to a low-Re destination produced fewer additional cases in the origin county and in the origin and destination county combined. For the destination county alone, it was preferable to accept evacuees from a low-Re origin. However, in a real hurricane landing, the counties that require evacuation are determined by the path of the hurricane, i.e. a low-Re origin cannot be stipulated. The length of evacuation and the number of people evacuating also influenced case numbers; however, these two characteristics are also expected to be shaped more by the specific circumstances necessitated by a particular hurricane rather than by public health directives.

Full model simulations of hurricane evacuation scenarios

Next we used a national county-scale metapopulation model and a suite of scenarios to further explore how transmission rates and hurricane evacuation affect COVID-19 incidence in origin and destination counties (Methods). All scenarios assume that a Category 3 hurricane is approaching southeast Florida and that people living in Palm Beach, Broward, Miami-Dade, and Monroe Counties are ordered to evacuate. Based on studies of evacuation compliance and behavior in this region for Category 3+ hurricanes, we estimate that 48% of each county’s population would evacuate7,10-17 (Methods). Assuming that 19% of evacuees relocate elsewhere within their respective counties, this leads to a total of 2.3 million evacuees leaving the four affected counties7,16-20.

For each scenario, evacuees were assigned to a different set of destination counties based on a list of 165 possible destinations across 26 states identified during post-Hurricane Irma surveys7,17. In the baseline scenario, evacuees were assigned to all 165 destination counties in proportion to observed evacuation choices during Hurricane Irma. The number of evacuees assigned to each destination county for this baseline scenario is shown in Figure 2a. In order to explore the effects of destination county transmission levels on the number of evacuation-associated COVID-19 cases, we further proposed two hypothetical scenarios in which evacuees were assigned to locations with high Re or low Re. In the high (low) Re scenario, 90% of the evacuees assigned to each county in the baseline scenario were instead diverted to the 82 counties with the highest (lowest) Re, weighted by the proportion of evacuees sent to each of these counties in the baseline scenario (Methods).

Figure 2:
Figure 2:

Simulations for evacuation using the national county-level transmission model. (a) The number of evacuees accepted by 165 destination counties in the baseline scenario. (b) The estimated effective reproductive numbers Re for both origin and destination counties on July 23rd, 2020. (c) Comparison of excess cases in origin and destination counties combined (left column), only origin counties (middle column) and only destination counties (right column) for the baseline, low and high evacuation scenarios. Simulations were performed for three settings: no increase (top row), 10% increase (middle row) and 20% increase (bottom row) of the transmission rates in destination counties. Box plots show the median and interquartile and whiskers show the 95% CIs. Asterisks indicate that excess cases are significantly lower or higher than the baseline scenario (Wilcoxon signed rank test, p < 10−5). Results are obtained from 100 model simulations.

Scenario projections were initiated from the model state calibrated to observed county-level COVID-19 case and death data from February 21st through July 23rd, 2020 (Methods). The estimated effective reproductive numbers (Re) in origin and destination counties at the start of simulations are shown in Figure 2b. Evacuees tend to stay with friends or family, in hotels/motels, or in public shelters, each of which would likely increase transmission opportunities relative to simply staying home7. To reflect this, we assume the COVID-19 transmission rate in destination counties increases during the evacuation period by either 0%, 10%, or 20%. These settings implicitly differentiate the levels of control effected in destinations during evacuation, as well as differences in transmission potential associated with different types of accommodation (e.g., staying with friends/families, hotels or shelters). In addition, to reflect periods of hurricane preparation and recovery21-23, we elevated the transmission rate in the origin counties by 20% beginning 3 days prior to evacuation and ending 3 days after the return of evacuees (more detailed simulation settings are provided in Methods). For comparison, we also generated simulations for the same period but without evacuation.

In all scenarios, combined cases in origins and destinations are primarily driven by ongoing local COVID-19 transmission dynamics (Figure S1; Table S1); however, evacuation does alter disease outcomes. In the baseline scenario, total COVID-19 cases in the origin and destination counties increase significantly relative to the no-evacuation scenario (Figure 2c; Table 1; Wilcoxon signed rank test), indicating evacuation in and of itself can cause a statistically significant increase of COVID-19 cases. As indicated by the low and high scenarios, the number of COVID-19 cases resulting from evacuation is significantly lower (higher) than the baseline scenario if evacuees are directed to counties with lower (higher) transmission rates. However, as the transmission rate in the destination counties increases, the differences between the low and high scenarios become less pronounced. This result indicates that the benefits of a directed evacuation would be amplified by more stringent control efforts in destination counties.

View this table:
Table 1:

Full metapopulation model simulation of the median number of excess cases in origin and destination counties for different evacuation scenarios (baseline, low, high, and optimized) and different increases of transmission rates (Re) in destination counties (no change, 10% increase, 20% increase). Note that the high and low Re scenarios are not subject to the constraint of destination capacity, whereas the optimized scenario takes into account a hypothetical capacity for each destination county.

Greedy optimization method to minimize COVID-19 cases

The model simulations for our hypothetical evacuation scenarios indicate that a strategic evacuation plan could reduce excess COVID-19 infections. However, these scenarios neither accounted for the accommodation capacities of destination counties nor provided a framework for optimizing evacuation plans.

To address these issues, we developed a greedy search optimization algorithm aimed at minimizing total excess COVID-19 cases by strategically assigning evacuees to optimal destination counties. As indicated by the two-county model simulations, evacuating individuals to destinations with low Re reduces COVID-19 transmission. However, given the varying prevalence of infection in origin counties and the nonlinear transmission dynamics, it is not straightforward to determine the optimal number of evacuees from each origin county that should be prioritized and redirected to each of the lowest-Re destinations. In conducting the optimization, we imposed the following constraints on human movement: 1) we assumed that a fraction of evacuees cannot be redirected from their baseline destination county, representing individuals whose choice of destination will not be influenced by evacuation directives, perhaps due to financial constraints or preferences to stay with family; and 2) we prescribed a capacity limit on the number of evacuees received by each destination county. The greedy search starts from an evacuation matrix representing the evacuees who cannot be redirected from their destination and then iteratively directs the remaining evacuees to destination counties with lowest Re. In each iteration, the algorithm selects which origin counties will be assigned the evacuee slots available in a destination county. This search is repeated for each successive destination county until all evacuees are assigned a destination (Methods, Supplementary Information).

We repeated this evacuation optimization with three different settings: no increase, 10% increase, and 20% increase of transmission rates in destination counties, again to reflect differences in control efforts and accommodation type. We assumed that 10% of evacuees will maintain their original destination and are thus not redirected, and that each destination has a capacity of 120% of the evacuees accepted in the baseline scenario. We then generated the optimized evacuation plan for each setting. During the optimization process, assigning more evacuees to low-Re counties led to a reduced number of total infections compared to the baseline scenario (Figure S2). The optimized top 20 destinations for the four origin counties are reported in Tables S2-S4. In Figure 3a, we show the change in the number of evacuees to each destination for the optimized evacuation scenario with a 10% increased transmission rate in destination counties. In general, evacuees who traveled to high-Re destinations in the baseline scenario were redirected to low-Re destinations. For all three transmission rate scenarios, the optimized evacuation effectively reduces the number of excess cases in both origin and destination counties compared to the baseline scenario (Figure 3b). The reduction is greatest (up to 30%) for the scenario in which there is no increase in the transmission rate in destination counties, which highlights the crucial role of effective intervention during evacuation. In this optimization example, the fraction of non-allocable evacuees and destination capacity are hypothetical, as these quantities are unknown. If such information were available or could be estimated using socioeconomic and geographic characteristics, the optimization could be tailored to reflect more realistic constraints on evacuation.

Figure 3:
Figure 3:

Optimization of evacuation plans. (a) The change in the number of evacuees to destination counties in the optimized evacuation plan compared with the baseline evacuation scenario. Evacuation was optimized for the setting in which transmission rates in destination counties increase by 10%. (b) Excess cases for the baseline and optimized evacuation scenarios are compared for the origin and destination counties combined (left column), only origin counties (middle column) and only destination counties (right column). Simulations were performed for three settings: no increase (top row), 10% increase (middle row) and 20% increase (bottom row) of the transmission rates in destination counties. Boxes and whiskers show the median, interquartile and 95% CIs. Asterisks indicate that excess cases are significantly lower than the baseline scenario (Wilcoxon signed rank test, p < 10−5). Results are obtained from 100 model simulations.

DISCUSSION

The results of this study have far-reaching consequences not only for hurricane evacuation this season but also for long-term US hurricane preparedness and evacuation planning.

Research suggests that people rely on past experiences when choosing their evacuation routes and destinations20,24. This study shows that excess COVID-19 cases could be minimized by instead directing evacuees to either counties with lower COVID-19 transmission rates or an optimized set of counties. While decisions about whether to evacuate and where to go ultimately fall to individual households, emergency communications from federal agencies and broadcast meteorologists can influence residents’ perceptions of hurricane threats and are seen as trusted sources of information in emergency situations25. Because a majority of US residents are concerned about COVID-1926, if the need for a large-scale evacuation arises, evacuees may turn to these same trusted sources for information on how best to stay safe while evacuating. Local, state, and federal officials who develop evacuation orders and communicate them to the general public may therefore want to consider whether their evacuation-related communications should include assessments of the relative safety of potential destination counties with respect to COVID-19 risk rather than allowing default evacuation patterns based on past storms to prevail.

This research shows that the magnitude of the impact of evacuation on COVID-19 caseloads is highly dependent on conditions in destination counties. The degree to which counties are prepared to host, isolate, and meet the needs of evacuees while also minimizing virus exposure through public health directives such as social distancing and mask wearing will be a key determinant of the impact of evacuation on COVID-19 case numbers. Preparedness within destination counties is particularly important because, as this analysis shows, destination counties will bear the brunt of the excess COVID-19 cases that result from an evacuation event. Destination counties must be aware of the influx of evacuees should it occur and must be allocated the financial and human resources needed to ensure the safety of both their residents and the evacuees they are sheltering.

The US response to the COVID-19 pandemic has varied widely from state to state and from county to county. As a result of policy, communication and ideological differences, compliance with mask wearing, for example, has varied substantially even within a given state27. This variability could extend into county-level hurricane preparedness measures, particularly given that guidance may be issued at the state level while implementation of specific measures is left to counties, as is the case for Florida’s current co-response guidance on hurricane evacuation and COVID-1928.

Centuries of systemic racism in the US have left Black, Native American, Latinx, and other non-White people with both higher exposure to and fewer resources to cope with environmental or health-related stressors compared with White populations29,30. For example, recent research suggests that federal financial aid after natural disasters is not equitably distributed among communities and may even exacerbate income inequality31,32. Low-income communities and communities of color consequently struggle to prepare in advance of and recover in the wake of natural disasters33,34.

Due to systemic health inequities, including higher exposure to air pollution29 and higher rates of underlying health conditions35, COVID-19 has also disproportionately affected Black, Native American, and Latinx people in the US36. These groups have experienced higher infection rates, poorer health outcomes, and deeper declines in employment during the pandemic36-38. The additional risks faced because of COVID-19 and the financial costs associated with evacuation39 could present additional challenges for these segments of the population during hurricane evacuation or discourage them from evacuating altogether.

This study has used idealized scenarios to model hurricane evacuation patterns. These scenarios cannot fully capture many household-level choices that could alter levels of social contact—and therefore potential COVID-19 exposure—during the evacuation period. For example, previous studies of hurricane evacuations along the US Southeast and Gulf Coasts show that evacuees strongly and consistently prefer to stay with friends and family over going to hotels/motels or public shelters7,14 15,20,39,40. Levels of social contact and potential virus transmission would likely differ across accommodation types, which implies a level of complexity and spatial heterogeneity that is not possible to incorporate within the model used in this study. Similarly, this study does not consider variable levels of exposure to COVID-19 based on evacuation transportation mode. While evacuees strongly prefer to travel in their own vehicles17, shared modes of transportation such as buses or carpools would increase potential virus transmission and exposure. Recent research suggests that for both transportation and shelter, the sharing economy— Internet-based transactions via companies like Airbnb that allow for peer-to-peer sharing of goods and services—could play a role in providing free or affordable resources to evacuees that would enable evacuees to maintain social distancing41.

As a result of the ongoing economic and physical toll of the pandemic, household-level decision making regarding evacuation may differ from that of past years. There are many sociodemographic factors associated with decision making around evacuation including experience with past hurricanes, length of residence, home ownership, age, income, race, employment status, level of social connectivity, social cues, perceived levels of self-efficacy and risk, and storm conditions42-46. While some of these factors (e.g. gender and race) are unchanged from last year, others (e.g. employment status and income) may be either changed or very much influenced by the current COVID-19 pandemic. In contrast, the idealized scenarios adopted for this study assume that people will choose destination counties and accommodation types that match past choices.

The movement of people in and out of hurricane-affected counties does not simply cease after all evacuees have returned to their homes. For instance, communities affected by hurricanes often experience an influx of workers who assist with rebuilding and recovery efforts47-49, which could also influence infection rates in the affected counties long after the evacuation period. Post-evacuation movement patterns are beyond the scope of the present study.

Critically, hurricane evacuation is intended to save lives and prevent serious injuries to residents of hurricane-prone regions. While this study evaluates excess COVID-19 cases resulting from evacuation, it does not evaluate non-COVID-19 related risks to human health and lives in the event that people choose to remain in their homes despite receiving evacuation orders-risks that could increase if people are afraid to evacuate out of concern for contracting COVID-19. Nor does it address evacuations of hospitals, nursing homes, prisons, or other facilities. It will be critical for emergency managers to factor in these— and other—complexities when developing plans50.

Finally, the results presented here are based on scenarios that, while plausible, are strictly hypothetical. While the overall notion that distributing evacuees to destination counties with low transmission rates minimizes excess cases should theoretically apply to geographies outside Florida or the US, additional model simulations of such scenarios should be generated.

CONCLUSION

The data presented here show that while a large-scale hurricane evacuation would increase the total number of COVID-19 cases in the US, directing evacuees to plausible destination counties with low COVID-19 transmission rates would minimize the excess cases induced by the evacuation event. These results have far-reaching implications for immediate emergency management and communications practices, as well as long-term disaster preparedness.

Faced with the prospect of tens of thousands of additional cases arising from a hurricane evacuation, states and counties at both ends of evacuation routes must be allocated the necessary financial and human resources required to meet evacuees’ needs while also ensuring community safety and health through measures intended to reduce COVID-19 transmission rates. Further, resource distribution must prioritize the nation’s most vulnerable groups.

Data Availability

Data and code are available on github.

https://github.com/SenPei-CU/Hurricane-COVID

METHODS

Hurricane evacuation scenario development

To develop a hypothetical hurricane evacuation scenario for southeast Florida, we drew from previously published hurricane evacuation studies focused on Category 3+ hurricanes that have affected the Southeast or Gulf Coast regions of the US. In this scenario, we assumed a Category 3 hurricane approaching the southeast Florida coast along a track that would necessitate evacuations from Palm Beach51, Broward52, Miami-Dade53, and Monroe54 Counties.

We obtained the population under mandatory evacuation orders in each county from GIS shapefiles of the zones of mandatory evacuation from a Category 3+ hurricane for each county55. We then calculated the percent of the population living within mandatory evacuation zones that would actually comply with evacuation orders by averaging the compliance rate from eight regionally relevant studies of evacuation behavior7,10-16. We found that the average evacuation order compliance rate from these studies was 66%. Because many people living outside of the mandatory evacuation zones voluntarily choose to evacuate during hurricanes as well, we used the same approach to calculate the percent of the population living outside of mandatory evacuation zones but within the affected counties that would evacuate. Based on four studies of this “shadow evacuation” phenomenon, we determined that an average of 47% of county residents outside mandatory evacuation zones would also evacuate7,13-15,17. The mandatory and voluntary evacuees together represent 48% of each origin county’s population. We then determined that 19% of these evacuees would relocate within their respective counties based on the average from four evacuation behavior studies7,16-20.

Finally, to determine the destination counties of evacuees from each of the four origin counties, we obtained raw post-Hurricane Irma survey data7,17. These data allowed us to identify both the destination counties and the percent of evacuees choosing each destination county. We then apportioned evacuees leaving the four origin counties to each destination county.

Two-county model of COVID-19 transmission

In order to identify the most sensitive factors driving the COVID-19 transmission during evacuation, we first ran simulations using a simple two-county model. This model describes the transmission dynamics of COVID-19 during a hurricane evacuation from an origin (county 1) to a destination (county 2). Mathematically, the transmission dynamics in the origin and destination are depicted by a susceptible-exposed-infected-recovered (SEIR) model. We simulate the disease transmission as a stochastic Markov process using the following equations:

Embedded Image Embedded Image Embedded Image Embedded Image Embedded Image

Here Ni, Si(t), Ei(t), Embedded Image, Embedded Image and Ri(t) are the total, susceptible, exposed, reported infected, unreported infected and recovered population in county i on day t; βi is the transmission rate in county i; μ is the relative transmissibility for unreported infections; Z is the average latency period; D is the average duration of infectiousness; α is the fraction of reported infections. The effective reproductive number, which quantifies the local transmission rate, is computed as Re = βiD[α + μ(1 − α)]Si/Ni using the next generation matrix approach.

During evacuation, we assume a fraction (peva) of the population is evacuated from county 1 to county 2 and mixes with the local population for Teva days. Individuals within each compartment are randomly drawn from the population in the origin. We track the infections in the evacuated population in county 2, which then return to the origin after the evacuation and mix with the population therein. To account for the increased human interactions associated with evacuating to shared living spaces7, we additionally assume the transmission rates in the origin and destination are elevated during a period that spans the evacuation process. Specifically, the transmission rate in the origin is increased by 20% starting from 3 days prior to the evacuation until 3 days after the return of evacuees. The transmission rate in the destination is increased by 20% during the evacuation.

In model simulations, we fixed the following parameters in Eqs. (1)-(5): total population N1 = N2 = 106; reporting rate α = 0.1; relative transmissibility μ = 0.64; latency period Z = 4 days; infectious period D = 4 days. Denote the daily reported cases in the origin and destination as casei. To initiate model simulations, we set Embedded Image and Embedded Image. Model simulations were generated for the following stages after day 0: 14 days of local transmission, 3 days of pre-evacuation (with elevated transmission rate in the origin), Teva days of evacuation (with elevated transmission rates in both the origin and destination counties), 3 days of post-evacuation (with elevated transmission rate in the origin), and 28 days of additional post-evacuation simulation.

We varied six parameters to generate a large number of parameter combinations for use in model simulation: β1, β2 = 0.1, 0.2,…,0.8; Teva = 3, 4,…,10 days; peva = 0.1,0.2,…,0.8; case1, case2 = 50,100,…, 400. In total, 86 = 262,144 parameter combinations were simulated. For comparison, we also ran a simulation without evacuation for each parameter combination and computed the percentage change of total cases in the origin and destination counties attributed to evacuation. We selected the top 10% of combinations that lead to the lowest percentage increase (or highest percentage reduction) of reported cases in the origin county, the destination county and both counties combined, and inspected the marginal distributions of those six parameters.

Simulating COVID-19 transmission in the US

Following the two-county model analysis, we conducted in-depth COVID-19 simulations using a nation-wide metapopulation SEIR model representing all 3,142 US counties9. In this model, disease transmission in each county follows SEIR dynamics, but is also influenced by movement to and from other counties. We consider two types of movement: daily work commuting and random movement. To simulate movement prior to March 15, 2020, we use information on county-to-county work commuting that is publicly available from the US Census Bureau. After March 15, the census survey data are no longer representative due to changes of mobility behavior in response to COVID-19 control measures.

Therefore, in simulating movement after March 15, 2020, we use estimates of the reduction of inter-county visitors to points of interest (POI) (e.g., restaurants, stores, etc.) to inform the decline of inter-county movement on a county-by-county basis. We generated these estimates using data from SafeGraph56. We further assume the number of random visitors between two counties is proportional to the average number of commuters between them. As population present in each county is different during daytime and nighttime, we model the transmission dynamics of COVID-19 separately for these two time periods. The model equations are presented in the supplementary information. Similar models have been used to simulate COVID-19 transmission in China57 and influenza transmission in the United States58,59. The local effective reproductive number is derived as Re = βiD[α + μ(1 − α)]Si/Ni. To account for reporting delays in COVID-19 case and death observations, we mapped simulated documented infections to confirmed cases using a separate observational delay model fitted to the US case data59.

We calibrated the transmission model against county-level case and death data reported from February 21, 2020 through July 23, 202056,60, which produced an estimate of model parameters and state variables. We then ran simulations representing the following stages: 3 days of pre-evacuation, 7 days of evacuation, 3 days of post-evacuation and 14 days after post-evacuation. For comparison, we also ran simulations without evacuation and increase of transmission rates in origin and destination counties. An ensemble of 100 trajectories were generated to represent the uncertainty arising from different initial conditions and stochastic dynamics.

In the modeled hurricane evacuation, Vji evacuees travel from origin i to destination j and mix with the local population for Teva = 7 days, before returning to origin i. As for the two-county model, we increased the transmission rate in origin counties by 20% during the 3-day pre-evacuation, 7-day evacuation and 3-day post-evacuation. Three scenarios in which the transmission rate in destination counties is increased by 0%, 10% and 20% were simulated to compare different effects of hosting evacuees on local disease transmission.

The greedy algorithm to optimize evacuation

We developed a greedy optimization algorithm aimed at minimizing total excess COVID-19 cases by strategically assigning evacuees to optimal destination counties. In the evacuation optimization, we assume that a fraction p of evacuees from an origin to a destination won’t change their evacuation plans (for reasons such as personal connections at the destination or budgetary limitations) and the capacity of accepting evacuees for each destination j is Cj. Denote the baseline evacuation matrix as V, where Vji represents the number of evacuees from origin i to destination j. The optimization objective is to assign the rest of evacuees (i.e., (1 − p) × ∑j Vji from origin i) to destinations in an optimal way that minimizes the total infections in both origin and destination counties. Finding the exact solution to this combinatorial optimization problem is computationally challenging due to the large number of options.

In this study, we use a practically feasible greedy optimization approach that prioritizes moving the un-assigned evacuees to destinations with low Re. Specifically, we start from the evacuation matrix pV that represents evacuees assigned a destination. In each step of greedy search, we run a series of simulations, each one filling the available evacuee slots in the destination with the lowest Re from one of the origin counties. We select the origin county that generates the minimum number of reported cases and assign them to the destination county. We repeat this greedy search for each successive destination county until all evacuees are assigned a destination. The pseudo-code for this greedy algorithm is provided in the supplementary information. In this study, we assume 10% of evacuees from an origin to a destination in the baseline evacuation matrix V cannot be reallocated (i.e., p = 0.1), and the capacity of each destination is 120% of the evacuees in the baseline scenario V (i.e., Cj = 1.2 ∑i Vji).

Supplementary Information

Transmission model for 3,142 US counties

We formulate COVID-19 transmission as a discrete Markov process during both day and night. Daytime transmission lasts for dt1 days and the nighttime transmission dt2 days (dt1 + dt2 = 1). Here, we assume daytime transmission lasts for 8 hours and nighttime transmission lasts for 16 hours, i.e., dt1 = 1/3 day and dt2 = 2/3 day. The transmission dynamics are depicted by the following equations.

Daytime transmission:

Embedded Image Embedded Image Embedded Image Embedded Image Embedded Image

Nighttime transmission:

Embedded Image Embedded Image Embedded Image Embedded Image Embedded Image

Here, Sij, Eij, Embedded Image, Embedded Image Rij and Nij are the susceptible, exposed, reported infected, unreported infected, recovered and total populations in the subpopulation commuting from county j to county i (i → j), where Embedded Image; β is the transmission rate of reported infections; μ is the relative transmissibility of unreported infections; Z is the average latency period (from infection to contagiousness); D is the average duration of contagiousness; α is the fraction of documented infections; θ is a multiplicative factor adjusting random movement; Embedded Image is the average number of commuters between counties i and j; and Embedded Image and Embedded Image are the daytime and nighttime populations of county i.

The pseudo-code for the greedy optimization algorithm

Input:

Origin i = 1,2,…,n

Destination j = 1,2,….,m, where Re(1) ≥ Re(2) ≥…≥ Re(m)

Evacuation matrix V = {Vji}, Vji is the number of evacuees from origin i to destination j in the baseline scenario

Capacity of evacuees that can be accommodated by each destination: Cj

The fraction of evacuees that can’t be reallocated for each origin-destination pair: p

Variables:

ui:: the current number of evacuees in origin i that could be reallocated to different counties

ν: the currently available destination county with lowest Re.

M: the current evacuation matrix, Mνi is the number of evacuees assigned from origin i to destination ν.

Initial conditions:

ν = m

ui = (1−p)∑jVji

M = pV

Algorithm:

While max(ui) > 0

For i = 1 to n

Mi = M: reallocating from origin i to ν

Embedded Image

Run projection using Mi

Infi: total infection in all origin and destination counties

End

Embedded Image

M = Mk

Embedded Image

If ∑j Mνj == Cν

End

End

Output M as the optimized evacuation matrix.

Fig. S1.
Fig. S1.

Comparison of total cases in the origin and destination counties combined (left column), the origin counties only (middle column) and destination counties only (right column) for the no-evacuation, baseline, low and high evacuation scenarios. Simulations were performed for three settings: no increase (top row), 10% increase (middle row) and 20% increase (bottom row) of transmission rates in destination counties. Box plots show the median and interquartile and whiskers show the 95% CIs. Asterisks indicate that excess cases are significantly higher than the no-evacuation scenario (Wilcoxon signed rank test, p < 10−5).

Fig. S2.
Fig. S2.

Evolution of the total cases (blue lines) and the number of assigned evacuees (red lines) in the greedy algorithm. Results are shown for the settings with no increase, 10% increase and 20% increase of transmission rates in destination counties. The optimization starts from an evacuation matrix 0.1 × V, where V represents the evacuation matrix in the baseline scenario.

View this table:
Table S1.

The median number of total cases in origin and destination counties for different evacuation scenarios (no evacuation, baseline, low and high) and levels of elevated transmission rates (Re) in destination counties (no change, 10% increase, 20% increase). In the no evacuation scenario (Re) in destination counties is not increased. Note that the median excess cases shown in Table 1 is not necessarily the difference of the median total cases between the evacuation and non-evacuation scenarios shown in Table S1. That is, median{totalcaseeva(i) − totalcasenoeva(i)) ≠ median{totalcaseeva(i)) - median{totalcasenoeva(i)), where totalcaseeva(i) and totalcasenoeva(i) are the total numbers of cases for evacuation and non-evacuation scenarios in the ith simulation.

View this table:
Table S2.

The optimized evacuation plan (no increase of transmission rates in destination counties). We show the top 20 destinations for each origin county.

View this table:
Table S3.

The optimized evacuation plan (10% increase of transmission rates in destination counties). We show the top 20 destinations for each origin county.

View this table:
Table S4.

The optimized evacuation plan (20% increase of transmission rates in destination counties). We show the top 20 destinations for each origin county.

ACKNOWLEDGEMENTS

The authors would like to acknowledge Stephen Wong, who graciously provided data that informed this analysis. Astrid Caldas, Rachel Cleetus, Juan Declet-Barreto, Adrienne Hollis, and Erika Spanger-Siegfried provided inspiration and insight during the course of the research. Work on this study was supported by NSF grant DMS-2027369, a gift from the Morris-Singer Foundation, and the generous support of the Barr Foundation, the Common Sense Fund, the Energy Foundation, Farvue Foundation, the Fresh Sound Foundation, the MacArthur Foundation, New York Community Trust, the Rauch Foundation, Sand County Charitable Trust, The Scherman Foundation, three anonymous funders, and members of the Union of Concerned Scientists.

REFERENCES

  1. 1.↵
    Phillips, C. A. et al. Compound climate risks in the COVID-19 pandemic. Nat. Clim. Change 10, 586–588 (2020).
    OpenUrl
  2. 2.↵
    Lau, H. et al. Evaluating the massive underreporting and undertesting of COVID-19 cases in multiple global epicenters. Pulmonology (2020) doi: 10.1016/j.pulmoe.2020.05.015.
    OpenUrlCrossRef
  3. 3.↵
    Kavanagh, M. M. et al. Access to lifesaving medical resources for African countries: COVID-19 testing and response, ethics, and politics. The Lancet 395, 1735–1738 (2020).
    OpenUrl
  4. 4.↵
    Seneviratne, S. I. et al. Changes in Climate Extremes and their Impacts on the Natural Physical Environment. in Managing the Risks of Extreme Events and Disasters to Advance Climate Change Adaptation (eds. Field, C. B., Barros, V., Stocker, T. F. & Dahe, Q.) 109-230 (Cambridge University Press, 2012). doi:10.1017/CBO9781139177245.006.
    OpenUrlCrossRef
  5. 5.↵
    Kossin, J., Knapp, K. R., Olander, T. L. & Velden, C. S. Global increase in major tropical cyclone exceedance probability over the past four decades. Proc. Natl. Acad. Sci. 117, 11975–11980 (2020).
    OpenUrlAbstract/FREE Full Text
  6. 6.↵
    Han, S. Y., Tsou, M.-H., Knaap, E., Rey, S. & Cao, G. How Do Cities Flow in an Emergency? Tracing Human Mobility Patterns during a Natural Disaster with Big Data and Geospatial Data Science. Urban Sci. 3, 51 (2019).
    OpenUrl
  7. 7.↵
    Wong, S., Shaheen, S. & Walker, J. Understanding Evacuee Behavior: A Case Study of Hurricane Irma. https://escholarship.org/uc/item/9370z127 (2018).
  8. 8.↵
    Roache, Madeline. Hurricane Dorian Evacuation Orders by County | Time. Time (2019).
  9. 9.↵
    Pei, S., Kandula, S. & Shaman, J. Differential Effects of Intervention Timing on COVID-19 Spread in the United States. medRxiv 2020.05.15.20103655 (2020) doi:10.1101/2020.05.15.20103655.
    OpenUrlAbstract/FREE Full Text
  10. 10.↵
    Baker, E. Warning and Response. in Hurricane Hugo 202-210 (National Academy Press, 1994).
  11. 11.
    Zhang, Y., Prater, C. S. & Lindell, M. K. Risk Area Accuracy and Evacuation from Hurricane Bret. Nat. Hazards Rev. 5, 115–120 (2004).
    OpenUrl
  12. 12.
    Stein, R. M., Duenas-Osorio, L. & Subramanian, D. Who Evacuates When Hurricanes Approach? The Role of Risk, Information, and Location. Soc. Sci. Q. 91, 816–834 (2010).
    OpenUrlCrossRefPubMed
  13. 13.↵
    Cutter, S. L., Emrich, C. T., Bowser, G., Angelo, D. & Mitchell, J. T. 2011 South Carolina Hurricane Evacuation Behavioral Study. http://artsandsciences.sc.edu/geog/hvri/sites/sc.edu.geog.hvri/files/attachments/HES_2011_Final_Report.pdf (2011).
  14. 14.
    Lindell, M. K., Kang, J. E. & Prater, C. S. The logistics of household hurricane evacuation. Nat. Hazards 58, 1093–1109 (2011).
    OpenUrl
  15. 15.↵
    Yin, W., Murray-Tuite, P., Ukkusuri, S. V. & Gladwin, H. An agent-based modeling system for travel demand simulation for hurricane evacuation. Transp. Res. Part C Emerg. Technol. 42, 44–59 (2014).
    OpenUrl
  16. 16.↵
    Martin, Y., Li, Z. & Cutter, S. L. Leveraging Twitter to gauge evacuation compliance: Spatiotemporal analysis of Hurricane Matthew. PLOS ONE 12, e0181701 (2017).
    OpenUrl
  17. 17.↵
    Wong, S. D., Pel, A. J., Shaheen, S. A. & Chorus, C. G. Fleeing from Hurricane Irma: Empirical Analysis of Evacuation Behavior Using Discrete Choice Theory. Transp. Res. Part Transp. Environ. 79, 102227 (2020).
    OpenUrl
  18. 18.
    Dow, K. & Cutter, C., Susan L. Emerging hurricane evacuation issues: Hurricane Floyd and South Carolina. Nat. Hazards Rev. 3, 12–18 (2002).
    OpenUrl
  19. 19.
    Dash, N. & Morrow, B. Return delays and evacuation order compliance: the case of Hurricane Georges and the Florida Keys. Glob. Environ. Change Part B Environ. Hazards 2, 119–128 (2000).
    OpenUrl
  20. 20.↵
    Wu, H.-C., Lindell, M. K. & Prater, C. S. Logistics of hurricane evacuation in Hurricanes Katrina and Rita. Transp. Res. Part F Traffic Psychol. Behav. 15, 445–461 (2012).
    OpenUrl
  21. 21.↵
    Yin, W. An Agent-based Travel Demand Model System for Hurricane Evacuation Simulation. (Virginia Polytechnic Institute and State University, 2013).
  22. 22.
    Noltenius, M. S. Capturing Pre-evacuation Trips and Associative Delays: A Case Study of the Evacuation of Key West, Florida for Hurricane Wilma. (University of Tennessee, Knoxville, 2008).
  23. 23.↵
    Lindell, M. K., Sorensen, J. H., Baker, E. J. & Lehman, W. P. Community response to hurricane threat: Estimates of household evacuation preparation time distributions. Transp. Res. Part Transp. Environ. 85, 102457 (2020).
    OpenUrl
  24. 24.↵
    Mesa-Arango, R., Hasan, S., Ukkusuri, S. V. & Murray-Tuite, P. Household-Level Model for Hurricane Evacuation Destination Type Choice Using Hurricane Ivan Data. Nat. Hazards Rev. 14, 11–20 (2013).
    OpenUrl
  25. 25.↵
    Lazrus, H., Morrow, B. H., Morss, R. E. & Lazo, J. K. Vulnerability beyond Stereotypes: Context and Agency in Hurricane Risk Communication. Weather Clim. Soc. 4, 103–109 (2012).
    OpenUrl
  26. 26.↵
    Dunn, A., Suite 800Washington & Inquiries, D. 20036USA202-419-4300 | M.-857-8562 | F.-419-4372 | M. As the U.S. copes with multiple crises, partisans disagree sharply on severity of problems facing the nation. Pew Research Center https://www.pewresearch.org/fact-tank/2020/07/14/as-the-u-s-copes-with-multiple-crises-partisans-disagree-sharply-on-severity-of-problems-facing-the-nation/ (2020).
  27. 27.↵
    Katz, J., Sanger-Katz, M. & Quealy, K. A Detailed Map of Who Is Wearing Masks in the U.S. The New York Times (2020).
  28. 28.↵
    Florida Division of Emergency Management. Florida Co-Response Pre-Landfall Tropical Weather Guidance (V2).pdf. (2020).
  29. 29.↵
    Bell, M. L. & Ebisu, K. Environmental Inequality in Exposures to Airborne Particulate Matter Components in the United States. Environ. Health Perspect. 120, 1699–1704 (2012).
    OpenUrlPubMedWeb of Science
  30. 30.↵
    Singh, G. K. et al. Social Determinants of Health in the United States: Addressing Major Health Inequality Trends for the Nation, 1935-2016. Int. J. MCH AIDS 6, 139–164 (2017).
    OpenUrl
  31. 31.↵
    Emrich, C. T., Tate, E., Larson, S. E. & Zhou, Y. Measuring social equity in flood recovery funding. Environ. Hazards 19, 228–250 (2020).
    OpenUrl
  32. 32.↵
    Howell, J. & Elliott, J. R. Damages Done: The Longitudinal Impacts of Natural Hazards on Wealth Inequality in the United States. Soc. Probl. 66, 448–467 (2019).
    OpenUrl
  33. 33.↵
    Cleetus, R., Bueno, R. & Dahl, K. A. Surviving and Thriving in the Face of Rising Seas. https://www.ucsusa.org/sites/default/files/attach/2015/11/surviving-and-thriving-full-report.pdf (2015).
  34. 34.↵
    Baker, E. J. Household preparedness for the Aftermath of Hurricanes in Florida. Appl. Geogr. 31, 46–52 (2011).
    OpenUrl
  35. 35.↵
    Colen, C. G., Ramey, D. M., Cooksey, E. C. & Williams, D. R. Racial Disparities in Health among Nonpoor African Americans and Hispanics: The Role of Acute and Chronic Discrimination. Soc. Sci. Med. 1982 199, 167–180 (2018).
    OpenUrl
  36. 36.↵
    Larsen, J., Wimberger, E., King, B. & Houser, T. A Just Green Recovery. https://rhg.com/research/a-just-green-recovery/ (2020).
  37. 37.
    Oppel, R. A., Gebeloff, R., Lai, K. K. R., Wright, W. & Smith, M. The Fullest Look Yet at the Racial Inequity of Coronavirus. The New York Times (2020).
  38. 38.↵
    Wu, X., Nethery, R. C., Sabath, B. M., Braun, D. & Dominici, F. Exposure to air pollution and COVID-19 mortality in the United States: A nationwide cross-sectional study. medRxiv 2020.04.05.20054502 (2020) doi:10.1101/2020.04.05.20054502.
    OpenUrlAbstract/FREE Full Text
  39. 39.↵
    Wu, H.-C., Lindell, M. K., Prater, C. S. & Huang, S.-K. Logistics of hurricane evacuation during Hurricane Ike. in Logistics: Perspectives, Approaches and Challenges (eds. Cheung, J. & Song, H.) 127-140 (Nova Science Publishers, 2013).
  40. 40.
    Bian, R., Wilmot, C. G., Gudishala, R. & Baker, E. J. Modeling household-level hurricane evacuation mode and destination type joint choice using data from multiple post-storm behavioral surveys. Transp. Res. Part C Emerg. Technol. 99, 130–143 (2019).
    OpenUrl
  41. 41.↵
    Wong, S. D., Walker, J. L. & Shaheen, S. A. Bridging the gap between evacuations and the sharing economy. Transportation (2020) doi: 10.1007/s11116-020-10101-3.
    OpenUrlCrossRef
  42. 42.↵
    Lazo, J. K., Bostrom, A., Morss, R. E., Demuth, J. L. & Lazrus, H. Factors Affecting Hurricane Evacuation Intentions: Factors Affecting Hurricane Evacuation Intentions. Risk Anal 35, 1837–1857 (2015).
    OpenUrl
  43. 43.
    Demuth, J. L., Morss, R. E., Lazo, J. K. & Trumbo, C. The Effects of Past Hurricane Experiences on Evacuation Intentions through Risk Perception and Efficacy Beliefs: A Mediation Analysis. Weather Clim. Soc. 8, 327–344 (2016).
    OpenUrl
  44. 44.
    Huang, S.-K., Lindell, M. K. & Prater, C. S. Who Leaves and Who Stays? A Review and Statistical Meta-Analysis of Hurricane Evacuation Studies. Environ. Behav. 48, 991–1029 (2016).
    OpenUrlCrossRef
  45. 45.
    Collins, J., Ersing, R., Polen, A., Saunders, M. & Senkbeil, J. The Effects of Social Connections on Evacuation Decision Making during Hurricane Irma. Weather Clim. Soc. 10, 459–469 (2018).
    OpenUrl
  46. 46.↵
    Metaxa-Kakavouli, D., Maas, P. & Aldrich, D. P. How Social Ties Influence Hurricane Evacuation Behavior. Proc. ACMHum.-Comput. Interact. 2, 1–16 (2018).
    OpenUrl
  47. 47.↵
    Fussell, E. Hurricane Chasers in New Orleans: Latino Immigrants as a Source of a Rapid Response Labor Force. Hisp. J. Behav. Sci. 31, 375–394 (2009).
    OpenUrlCrossRefWeb of Science
  48. 48.
    Theodore, N. After the Storm: Houston’s Day Labor Markets in the Aftermath of Hurricane Harvey. https://greatcities.uic.edu/wp-content/uploads/2017/11/After-the-Storm_Theodore_2017.pdf (2017).
  49. 49.↵
    Jordan, M. Hurricane Chasers: An Immigrant Work Force on the Trail of Extreme Weather - The New York Times. The New York Times (2019).
  50. 50.↵
    Wong, S. & Shaheen, S. Double the Trouble: Evacuations During COVID-19. Medium https://medium.com/@SusanShaheen1/double-the-trouble-evacuations-during-covid-19-916626c37cf4 (2020).
  51. 51.↵
    Palm Beach County EMS. Map A.9. Coastal evacuation zones & routes.
  52. 52.↵
    Evacuation zones, routes and shelters. (2016).
  53. 53.↵
    Miami Dade County. Storm Surge Planning Zones. https://www.miamidade.gov/global/emergency/hurricane/storm-surge-zones.page (2020).
  54. 54.↵
    Florida State Emergency Response Team & South Florida Regional Planning Council. Chapter VI: Regional evacuation transportation analysis. in Statewide Regional Evacuation Study Program (2019).
  55. 55.↵
    Florida Hurricane Evacuation Zones (FeatureServer). ArcGIS REST Services Directory https://services1.arcgis.com/CY1LXxl9zlJeBuRZ/arcgis/rest/services/Florida_Hurricane_Evacuation_Zones/FeatureServer (2019).
  56. 56.↵
  57. 57.↵
    Li, R. et al. Substantial undocumented infection facilitates the rapid dissemination of novel coronavirus (SARS-CoV-2). Science 368, 489–493 (2020).
    OpenUrlAbstract/FREE Full Text
  58. 58.↵
    Pei, S., Kandula, S., Yang, W. & Shaman, J. Forecasting the spatial transmission of influenza in the United States. Proc. Natl. Acad. Sci. 115, 2752–2757 (2018).
    OpenUrlAbstract/FREE Full Text
  59. 59.↵
    Pei, S. & Shaman, J. Initial Simulation of SARS-CoV2 Spread and Intervention Effects in the Continental US. medRxiv 2020.03.21.20040303 (2020) doi:10.1101/2020.03.21.20040303.
    OpenUrlAbstract/FREE Full Text
  60. 60.↵
    US Coronavirus Cases and Deaths. USAFacts.org/visualizations/coronavirus-covid-19-spread-map (2020).
Back to top
Posted August 11, 2020.
Download PDF
Data/Code
Email
Compound risks of hurricane evacuation amid the COVID-19 pandemic in the United States
Sen Pei, Kristina A. Dahl, Teresa K. Yamana, Rachel Licker, Jeffrey Shaman
medRxiv 2020.08.07.20170555; doi: https://doi.org/10.1101/2020.08.07.20170555
Twitter logo Facebook logo LinkedIn logo Mendeley logo
Citation Tools

Subject Area