Recovery Migration After Hurricanes Katrina and Rita: Spatial Concentration and Intensification in the Migration System

Abstract

Introduction

Hurricanes, typhoons, and other extreme environmental events place the growing proportion of the world’s population living in coastal cities at risk of displacement (Goodess 2013; McGranahan et al. 2007). The focus on the displacing potential of these events in existing research on disaster-related migration neglects the longer-run population change shaped by recovery migration (Black et al. 2013). Furthermore, the focus on potential displacement leads to exaggerated estimates of the number of “climate refugees,” which should be tempered by more theoretically informed and empirically based research (Gemenne 2011). Disaster-affected places rarely experience permanent population loss (Laczko and Aghazarm 2009; McLeman 2011; Wisner et al. 2004). Instead, they recover their populations through return migration and new in-migration (Fussell and Elliott 2009; McLeman and Smit 2006) and, we propose, adaptations of the migration system.

For scholars concerned with global climate change, the effects of Hurricanes Katrina and Rita on New Orleans and the Gulf Coast region provide an example of what could happen to coastal and nearby cities affected by hurricanes and coastal flooding (e.g., Adamo 2010). These hurricanes struck the Gulf of Mexico between Texas and Florida within weeks of each other in 2005. Although hurricanes and other damaging environmental events are not rare for this region, Katrina was the sixth most powerful and costly hurricane thus far recorded; and Rita ranked fourth most powerful, although it struck a less-populated region, so damage estimates were lower (Knabb et al. 2006). These devastating events left residents, politicians, planners, and scholars wondering whether and how the region would recover its built environment and population (Kates et al. 2006) as well as what it meant for areas that might confront similar disasters in the future.

Our study focuses on population recovery through migration to a disaster-affected region. We work from a migration systems framework to ask a foundational question about resiliency: Does an environmental shock alter the preexisting migration system of the affected region? Considering whether environmental events affect existing patterns of migration over time is vital to understanding large-scale and long-run impacts of environmental events on human populations (e.g., Hsiang et al. 2013). Most research on the demographic effects of Hurricane Katrina on New Orleans focuses on the unequal vulnerability of residents to displacement by sociodemographic and place-based characteristics (Cutter and Emrich 2006; Fussell et al. 2010; Groen and Polivka 2010; Myers et al. 2008). However, to our knowledge, no research has considered how the disaster impacted the broader system of migration flows to and from disaster-affected coastline counties over the more prolonged recovery period. Moreover, whereas migration systems research chiefly concerns identifying the system and factors that perpetuate it (e.g., Fawcett 1989; Kritz et al. 1992; Mabogunje 1970; Massey et al. 1998), our study concerns the dynamic response of the system to an exogenous shock. Indeed, our study is distinct from previous research in two ways. First, we consider places directly affected by an environmental event and places to which they are connected through migration. Second, we assess longer-run impacts by analyzing migration in the recovery period, as opposed to the immediate post-event period. In doing so, we speak to the broader spatial and temporal impact of environmental events on human settlement patterns. Our results show that population recovery occurred through a spatial concentration and intensification of the migration system in the years following the disaster.

A Systems Model of Recovery Migration

We take a systems approach to investigate recovery migration (Fussell et al. 2014). In this approach, dating to Ravenstein’s (1885) study of migrant streams and counter-streams in nineteenth century United Kingdom, the entire migration system is the object of study as opposed to individual migrants or their places of origin or destination (Fawcett 1989; Lee 1966). The central proposition is that when one place within the system experiences a change—that is, an environmental shock—its effects are felt throughout the entire system (Andrienko and Guriev 2004; Bakewell 2013b; Mabogunje 1970). A migration system is defined by both structure and process. The structural element of a migration system is the ties connecting places, which are the basis for measuring the size and attributes of flows of migrants between them (Mabogunje 1970). The process element of a migration system is the dynamics governing the ties—that is, the “rules of the game that govern . . . elements in the system” (Bakewell 2013a:15). The ties, flows, and their attributes and relationships interact to perpetuate and reinforce the system by encouraging migration and other types of exchanges (e.g., capital, commodities, information) along certain pathways and discouraging it along others (Mabogunje 1970:12; see also Kritz et al. 1992 and McHugh 1987). Although stability in the system over time and across space is often emphasized in work adopting a systems perspective (DeWaard et al. 2012; Massey et al. 1998), some scholars underscore the dynamic element of the system, at least conceptually, and focus on factors altering system elements—and, it follows, the migration system itself (Bakewell 2013a; Bennett and Haining 1985; Bennett et al. 1985; de Haas 2010; Fawcett 1989; Plane 1987; Plane and Rogerson 1986). There are several different approaches to viewing change in a migration system, including a focus on the speed of change (Bennett and Haining 1985) and the role of historical factors in shaping the type of change (DeWaard et al. 2012; Mueser 1989; see also Bell et al. 2002 on measuring connectivity).

Underlying the structure and processes defining a migration system are the decisions and behaviors of individuals and households. At the aggregate level, internal migration is often framed as labor migration, with populations redistributing from places with diminishing or less economic opportunities to areas with expanding or more opportunities (e.g., Greenwood 1997). Correspondingly, at the individual and household level, internal migration is viewed as an economic decision in which the costs and benefits of a potential move are weighed, and migration occurs when perceived gains exceed anticipated costs (e.g., Sjaastad 1962; Todaro 1976). Factors beyond labor ties and economic pulls also shape migration decisions. Personal networks lower barriers to migrating and inform destination selection in two ways: (1) by providing information about and access to housing, employment, and other resources at the destination; and (2) by reducing the psychic costs of migration (Greenwood 1969; Lansing and Mueller 1967; Nelson 1959). Such information, access, and costs are correlated with distances between origin and destination (Schwartz 1973; Sjaastad 1962). An exogenous shock to the migration system, such as a devastating hurricane, modifies the economic and social processes shaping decisions to move and destination choices and, we expect, the migration system itself.

As is widely noted, environmental migration is rarely distinguishable from other types of migration—that is, labor migration or network migration—because environmental forces operate indirectly through economic, political, and social structures (Black et al. 2011; see also Hunter 2005; McLeman and Smit 2006; and Perch-Nielsen et al. 2008). However, the spatial nature of environmental migration may alter the patterns of migration with a system. Drawing from a growing empirical literature on environmental drivers of migration, Findlay (2011:S51–S52) extracted several principles regarding environmental migrants’ destination choices that are consistent with general migration theory. First, most potential migrants prefer not to move. However, after a decision to move has been made, the second principle states migrants will move relatively short distances. Finally, to summarize the third through sixth principles, migrants prefer to go to places where they already have ties, allowing them to more easily and profitably exchange their human, social, and cultural capital. For most migrants, these places are nearby, although a few more advantaged migrants follow historical, cultural, or economic ties to more distant destinations. McLeman and Hunter (2010) arrived at similar generalizations by culling evidence from four cases of “climate migration.” They added that such migration is rarely permanent. Barring an environmental change destroying housing and livelihoods, few places are ever completely abandoned (McLeman 2011). From these propositions, we argue that return migration and new in-migration to an area after an environmental shock are to be expected, and the principles governing destination choice summarized by Findlay (2011) help identify the types of places that become the likely origins of recovery migration.

To understand how an environmental shock changes a migration system through recovery migration, the unit of analysis must shift from the household to the flows of households into affected places, and the time frame must include a longer time horizon than is typical of disaster research. In the pre-disaster period, the structure of the migration system results from migrant households’ destination choices under a general migration regime. In the period immediately before and after a rapid-onset disaster, households are likely to follow established ties to nearby places within the pre-disaster migration system. To a large extent, this response to an environmental shock defines the system structure in the recovery period for displaced households. In the recovery period, the processes driving migration to the disaster-affected region are economic and social, and are only indirectly environmental.¹ Displaced households attempt to minimize losses, and households may even search for opportunities in the recovering economy. Thus, a migration system affected by an environmental shock exhibits three distinct periods: a pre-disaster period, the rapid-onset disaster period, and the recovery period.

An additional consideration for models of recovery migration is that not all environmental events or changes are the same; therefore, any expectations of their impact on a migration system—short-term and, in our case, longer-term effects—must consider the character of the environmental shock in relation to migration decision processes (Black et al. 2011, 2013; Hunter 2005; McLeman and Smit 2006). Environmental events are often distinguished by the speed of onset. Recovery migration is more likely after a rapid-onset and short-duration event than a slow-onset and long-duration one because of the contrasting influence on migrants’ decisions. Potential migrants facing a rapid-onset event have comparably less time to make migration decisions and tend to make temporary moves (i.e., evacuations or short-term displacements), returning to their communities and livelihoods when the built environment is restored. These return migrations might persist for multiple years following the event, depending on the pace of rebuilding efforts (Kates et al. 2006). In contrast, potential migrants confronting slow-onset events face ongoing decisions about whether to engage in cyclical or permanent out-migration as a means of coping with losses affecting their livelihoods (Laczko and Aghazarm 2009). With storm warnings reported only days in advance, Hurricanes Katrina and Rita were rapid-onset events and extremely destructive for coastline counties affected by the storm surge. In these places, recovery of homes and livelihoods was prolonged or impossible. Thus, we expect recovery migration to concentrate in the harder hit areas in the years following the hurricanes.

Our conceptual model of recovery migration is rooted in research on historical environmental events and accumulating empirical evidence on migratory flows from and to Gulf of Mexico counties affected by the 2005 hurricane season. In their historical analysis, McLeman and Smit (2006) focused on Oklahoma and the 1930s Dustbowl Migration to describe several types of migration flows affecting the population size and composition of Oklahoma counties in subsequent years. Best known is the out-migration of Dustbowl refugees, mostly displaced tenant farmers, to rural California. Less noted are the return migration of some Dustbowl refugees and the migration of rural residents to nearby cities and towns. The key points are that (1) out-migration is not the only dynamic characterizing a migration system after an environmental event, and (2) return migration and new in-migration are important components of population recovery (Fussell and Elliott 2009).

Research specific to the 2005 hurricane season shows that the majority of residents in the most threatened Gulf of Mexico coastal counties evacuated in anticipation of Hurricane Katrina’s landfall (Elliott and Pais 2006; Groen and Polivka 2010; Haney et al. 2010). From 2006 through 2010, in-migration made several of the most devastated counties the fastest growing counties in the United States (U.S. Census 2008; U.S. Census Bureau 2011). Migratory in-flows included returning residents as well as newcomers, such as “hurricane chasers” seeking construction employment, young professionals pursuing urban development and entrepreneurial opportunities, energy sector workers repairing the damaged industrial infrastructure, and displaced residents from more severely affected counties within the recovering disaster-affected region (Ehrenfeucht and Nelson 2013; Fussell 2009).

In our study, we take as the unit of analyses the flows of migrant households among disaster-affected places and the places tied to them in order to understand how the migration system of Hurricanes Katrina- and Rita-affected coastline counties in the Gulf of Mexico changed between the pre-disaster and recovery periods. We examine recovery migration that brought population back to the disaster-affected areas several years after the disaster, not shorter-term moves better characterized as temporary evacuations. We test several hypotheses elaborated in the following section; if supported, these hypotheses will contribute to a general understanding of how environmental shocks change the drivers of migration systems. In this case, the migration system adapted to generate recovery migration.

Hypotheses About a Recovery Migration System

We investigate three qualities of the migration system of coastline counties most severely affected by Hurricanes Katrina and Rita to assess three corresponding hypotheses. First, we examine whether the migration system remains stable or changes after facing an environmental shock. Stability in the system would result if disaster-affected coastline counties’ populations migrated only to places in the pre-disaster migration system and if these were the only places sending in-migrants to disaster-affected counties in the recovery period. In this case, we would observe the same ties between counties in both periods. Alternatively, we might expect some change in the system given the large-scale and involuntary nature of population displacement from disaster-affected coastline counties in the immediate aftermath of the hurricanes. Some residents of these counties might have voluntarily or involuntarily relocated to counties that were not part of the pre-disaster migration system, thereby introducing new ties to the recovery migration system. New ties may also be introduced if people from counties outside the pre-disaster migration system relocated to disaster-affected coastline counties in search of work or other recovery-led opportunities. We examine the extent of stability and change in the migration system through an analysis of ties between specific county pairs occurring only in the pre-disaster or the recovery periods, but not both: the smaller the number of these unique ties in a given period, the more stable the migration system between the two periods.

Second, we further examine stability and change in the migration system by analyzing the size of in-flows among all ties within the pre-disaster and recovery periods. Environmental and disaster-driven migrations are shaped by the nature of the environmental change in the origin community. Given that hurricanes are rapid-onset, short-duration events—and, in this case, resources were available for recovery—we expect that displaced residents will return as the recovery progresses, and new in-migrants will arrive to pursue emerging opportunities. If recovery is underway, we would expect to see larger in-flows to disaster-affected coastline counties in the recovery period than in the pre-disaster period, and comparatively smaller out-flows in the recovery period. However, if recovery is faltering, there would be smaller in-flows of displaced residents and opportunity seekers into these counties, and relatively larger out-flows.

Third, extending the previous analysis, we examine change in the average size of in-flows to disaster-affected coastline counties from proximate and urban counties between the pre-disaster and recovery periods. If the principles of migrant destination choice are correct, disaster-affected migrants are most likely to have relocated to proximate counties and urban counties where they were best able to benefit from their human, social, and cultural capital. Therefore, as displaced residents return from these places in the recovery period, we expect the average size of these flows to be larger than they were in the pre-disaster period. New migrants originating in counties connected in the pre-disaster migration system are also likely to contribute to in-flows as they pursue new opportunities in recovering areas.

Data

We define three geographic regions: (1) disaster-affected coastline counties, (2) nearby Gulf of Mexico coastal counties, and (3) distant counties in the continental United States. We also identify urban counties within each of these regions because our third hypothesis includes consideration of urban counties. One could define regions in a number of ways. We focus on counties’ relationship to water because the storm surge accompanying a hurricane is extremely destructive. Storm surges bring powerful waves into coastline areas and push water through rivers and other water ways, destroying and damaging buildings and infrastructure. We identify these regions by combining the definitions used in a U.S. Census Bureau report on coastline population trends (Wilson and Fischetti 2010) with Federal Emergency Management Agency (FEMA) disaster declarations. In the U.S. Census Bureau report, a county adjacent to coastal waters or territorial seas is designated a coastline county and is part of a subset of coastal counties. A coastal county has at least 15 % of its land within the nation’s coastal watershed or a coastal cataloging unit (National Oceanic and Atmospheric Administration n.d.). After Hurricanes Katrina and Rita, 36 coastline counties were declared federal disaster areas by FEMA. We label these “disaster-affected coastline counties,” and they are the focal region in our analysis. Our second region includes 124 counties that are either coastline counties that were not declared federal disaster areas or coastal counties that may or may not have been declared federal disaster areas. The third region includes the 2,951 remaining distant counties in the continental United States. Slightly less than one-half (1,297) of all counties in the continental United States are urban. We refer to these three regions as “disaster-affected coastline counties,” “nearby counties,” and “distant counties,” as shown in Fig 1.

Our study concerns the connections between places, which we identify as ties between origin-destination pairs of counties and the size of migrant flows across these ties. A tie refers to the presence of a flow of any size, whereas a flow refers to the number of households migrating from i to j. We focus on in-flows and out-flows between pairs—not net flows to or from sending or receiving counties—because the meaning of the flow depends on its directionality (Rogers 1990). In our case, out-flows from disaster-affected coastline counties to counties within each region capture the out-migration dimension of the system. The in-migration dimension of the system is measured by in-flows to disaster-affected coastline counties from counties within each region. Based on existing research, we assume that in-flows to disaster-affected coastline counties consist of returning residents (e.g., Fussell et al. 2010; Groen and Polivka 2010; Myers et al. 2008), newcomers attracted by opportunities related to the recovery (Fussell 2009), and disaster-affected residents migrating between counties within the region. Together, these migrants comprise recovery migration. We examine the size of in-flows and the number of ties from other counties to disaster-affected coastline counties in the recovery period to test our hypotheses about environment-induced change in the migration system from a longer-run perspective.

We measure migration flows and their attributes with the Internal Revenue Service (IRS) Statistics of Income Division (SOI) County-to-County Migration Data files. The data include all U.S. federal income taxpayers, thereby underrepresenting the very poor and older populations, who are less likely to file income tax returns or be included as dependents on others’ tax returns, as well as the small percentage of tax returns filed after late September of the filing year (Gross n.d.). The data lack information about migrants other than their household income and broad age groups. Despite these limitations, researchers agree that the IRS migration data are the best available source for tracking changes in internal migration in the United States (Engels and Healy 1981; Isserman et al. 1982; Molloy et al. 2011). The Current Population Survey (CPS) indicates that in each year between 1992 and 2009, approximately 87 % of household heads filed tax returns, making the IRS data reliable for identifying population-level trends (Molloy et al. 2011). Although some researchers have used adjustment procedures to improve coverage of the IRS data (i.e., Plane 1999), there is no clean way of making such adjustments; moreover, the issue of coverage is just one limitation common to migration estimates (Raymer et al. 2013). These data are ideal for our study because they capture annual intercounty migration flows pre-dating and following the 2005 hurricane season in the style of a natural experiment. Evacuation behavior is not our interest, and it is better measured by the American Community Survey (ACS) (e.g., Johnson et al. 2008; Koerber 2006), CPS (e.g., Groen and Polivka 2010), or by specialized surveys (e.g., Sastry 2009). These data sources are inadequate for our analysis given their limited geographic representation or time frame.

For both conceptual and practical reasons, we use the years 1999–2004 to measure the period before the 2005 hurricane season (the pre-disaster migration system) and 2007–2009 to measure the recovery period (the recovery migration system). We do not examine the disaster period (2004–2006) because migration during this time is conceptually difficult to distinguish and because the data are of poorer quality. Johnson et al. (2008) found a general decline in match rates (i.e., coverage) between the tax filing years 2004–2005 and 2005–2006, which was greatest in areas affected by Hurricanes Katrina and Rita. By comparing the pre-disaster and recovery periods, we compare periods in which migration systems were relatively stable and data are of comparable quality. Our approach of averaging across individual tax filing years within the pre-disaster and recovery periods produces annualized estimates, thereby negating possible problems associated with imbalanced samples, as well as those associated with using just one tax filing year.

Methods

Our methodological approach moves beyond description to test hypotheses concerning recovery migration patterns in a natural experiment framework. As such, we address two current problems in research on population-environment interactions. First, we focus on population-level patterns to examine change over time using a cross-sectional data series. Second, we use smaller geographic units than many studies examining local-level responses to environmental change, which typically focus on the region as a whole (e.g., Grübler et al. 2007; Lutz et al. 2007). Moreover, by using all counties in the contiguous 48 U.S. states, we more completely represent the migration system; our study thus complements Curtis and Schneider’s (2011) approach by connecting areas directly and indirectly affected by an environmental shock through recovery migration.²

Modeling Migration Systems

Our analysis adapts two methods for studying geographic mobility: (1) descriptive estimates and maps modeling migration systems, and (2) confirmatory tests of hypotheses based on regression models. For our models of migration systems, we use the IRS data to develop estimates and maps of changes in the Gulf of Mexico migration system taking place between the pre-disaster (1999–2004) and recovery (2007–2009) periods. Characterizing these changes requires modeling migration systems in a way that simultaneously considers the population of households that were “at risk” of migrating in each sending county (Rogers 1975, 1990, 1995) from the vantage points of both sending and receiving counties (DeWaard 2013; DeWaard and Raymer 2012).

The sum of the last column vector in Eq. (3) is a count of the number of households from our starting hypothetical population of U.S. households in Eq. (1), which, in fact, migrated to disaster-affected coastline county j. Dividing this quantity by the size of the hypothetical population of households at risk of migrating to j—that is, the trace of the matrix in Eq. (1)—gives the proportion, $p$⁠, of households that were at risk of migrating to j and actually did so: that is, as governed by the probabilities in Eq. (2).⁴ Subtracting this quantity for the pre-disaster period from that for the recovery period, we derive an estimate of how the system of migration flows from all sending counties in the contiguous United States to disaster-affected coastline county j changed over time while accounting for the risk of migration in each sending county. We repeat these steps for each disaster-affected coastline county and present the combined results in tables and maps.

In addition to modeling migration to disaster-affected coastline counties from the vantage point of receiving areas, we likewise consider migration to disaster-affected coastline counties from the vantage point of sending counties, and from disaster-affected coastline counties to all counties in the contiguous United States. With respect to the former, for each row in Eq. (2), we sum those elements where receiving county j is a disaster-affected coastline county and subtract this quantity for the pre-disaster period from that for the recovery period. We map these results to show how migration to disaster-affected coastline counties, from the vantage point of sending counties, changed over time. To further model migration from disaster-affected coastline counties to each county in the contiguous United States, we revise the approach in Eqs. (1)–(3) for each U.S. county j, with the matrix in Eq. (1) respecified so that the starting population in each U.S. county is set to 0, excluding disaster-affected coastline counties in which the starting population is arbitrarily set to 1,000 households. For each U.S. county j, the sum of the last column vector in Eq. (3) gives an estimate of the number of households from our hypothetical population in Eq. (1) that, in fact, migrated to j from disaster-affected coastline counties. We then compare the resulting figure for the pre-disaster period to the corresponding figure for the recovery period. As before, we repeat the steps for each U.S. county j, and report the combined results.

Hypothesis Tests

In the second part of our analysis, we seek to determine whether the change in the number of ties and the size of migration flows in the disaster-affected coastline counties’ migration system between the pre-disaster and recovery periods is in the predicted direction and statistically significant. Our data allow us to construct an experimental framework offering counterfactuals used to distinguish a secular time trend from changes due to the treatment of interest (e.g., exposure to storm surge from Hurricanes Katrina and Rita). This requires examining the experiences of a comparison group. In our study, we define the comparison group as the nearby Gulf of Mexico coastal counties and the experimental group as the disaster-affected coastline counties. Although the nearby counties are unlikely to be identical to the disaster-affected coastline counties, a comparison of the two groups is the closest approximation of a natural experiment.⁵

For our first hypothesis concerning change in ties, we assess whether the number of all ties and the number of unique ties—ties that exist in only one of the two periods—differ between the pre-disaster and recovery periods. We focus on unique ties as opposed to all ties because the former, by definition, are indicative of a change in the migration system. To observe change, we compare the proportion of all possible ties in the pre-disaster period with those in the recovery period, where the number of all possible ties corresponds with the number of sending counties multiplied by the number of receiving counties in the specific region (less one because a county cannot be “tied” to itself). For example, for flows between disaster-affected coastline counties, there are 1,260 (36 × 35) possible ties; comparatively, for flows to disaster-affected coastline counties from nearby counties, there are 4,464 (36 × 124) possible ties. We test whether the difference in the proportion of unique ties is statistically significant for both in-flows and out-flows using a two-sample difference in proportions test for a comprehensive analysis of system contraction or expansion.

For our second and third hypotheses relating to the size of migration flows, our outcome of interest is the percentage change in the size of migration in-flows to disaster-affected coastline counties (experimental group) relative to nearby counties (comparison group) between the pre-disaster and recovery periods. In these analyses, we a use regression approach to build on the descriptive results and formally test for changes (i.e., contraction and expansion) in the system, focusing on changes in the size of flows across ties. The advantage of our approach is that it produces estimates of percentage change in the size of flows across ties in the system between the pre-disaster and recovery periods.

Results

Changes in the Migration System of the Disaster-Affected Coastline Counties

Our first task is to identify which counties were connected to disaster-affected coastline counties before Hurricanes Katrina and Rita so that we can describe the pre-disaster and recovery migration systems. We examine the migration system from the perspective of disaster-affected coastline counties as the ties through which migrants flow to or from those counties and the attributes of those ties—specifically, their size, both total and average, and the types of counties they connect.

To determine whether there is stability in the migration system, our first hypothesis, we compare the number of unique ties of disaster-affected coastline counties in the pre-disaster and recovery periods. If the system is perfectly stable, there will be no unique ties in either period because the ties will be common to both the pre-disaster and recovery periods. If the system is expanding, there will be more unique ties in the recovery period than the pre-disaster period; if it is contracting, there will be fewer. We find a 57.8 % decrease in the number of unique out-ties from disaster-affected coastline counties to all types of counties (Table 1, out-ties), with no significant change (–3.3 %) in the number of unique in-ties to disaster-affected coastline counties (Table 1, in-ties). This suggests that the recovery migration system of disaster-affected coastline counties contracted with respect to out-ties but not in-ties. The pattern is the same when we limit the data set to only urban counties.

The spatial concentration of the out-flow side of the system is evident by distinguishing the ties by the proximity of the counties they connect. Change was largest and statistically significant for distant counties (–63.1 %), followed by nearby counties (–43.3 %). The decrease in out-ties to other disaster-affected coastline counties was also large (–34.8 %) but not statistically significant. Change in the in-flow side of the system is also evident when considered this way. Unique in-ties to nearby counties grew by 33.3 %, and in-ties among disaster-affected coastline counties (–34.8 %) and, to a lesser extent, in-ties from distant counties (–6.8 %) were eliminated. Although this change was not significant overall, when we narrow the data to only urban origin counties, which make up the majority of ties in the migration system, we observe a statistically significant increase of 41.8 % in in-ties from nearby urban counties.

The significant changes in ties between types of places allow us to reject the hypothesis of stability in the migration system. The in-flow side of the system increased unique in-ties with nearby urban counties, consistent with assumptions that these were the counties most likely to have sheltered long-term displaced residents and to have provided recovery period workers. Further, we see the out-flow side of the system withdrew out-ties to nearby and distant counties, which is likely due to population losses suffered by disaster-affected coastline counties that in turn reduced the likelihood of out-migration in the recovery period. These results describe a recovery migration system that was more spatially concentrated than the pre-disaster system.

Changes in the Size of In-migration Flows to the Disaster-Affected Coastline Counties

Our second set of hypotheses concerns the size of the in-flows to disaster-affected coastline counties in the recovery period. If recovery migration is strong, we expect to see larger in-flows in the recovery period than in the pre-disaster period, and smaller out-flows than in-flows in the recovery period. If recovery migration is weak, we should see smaller or no different in-flows in the recovery period than in the pre-disaster period, and larger out-flows than in-flows in the recovery period. The descriptive evidence shows the total flow size into disaster-affected coastline counties grew by 19.4 % overall, and was larger than out-flows from these counties (144,854 vs. 137,424) (Table 2).

We conclude that recovery migration was strong. In-flows from nearby counties increased the most, by 30.1 %, although they were followed closely by in-flows from distant counties, which grew by 25.9 % (Table 2). These increases are somewhat larger for in-flows from urban counties (32.7 % and 26.4 %, respectively). In contrast, and as a point of comparison, out-flows from disaster-affected coastline counties increased relatively little, by 4.6 %, with the largest flows going to other disaster-affected coastline counties (8.2 %) and nearby counties (9.2 %), while flows to distant counties actually diminished (–1.3 %). The patterns are similar for urban counties. These results are consistent with the second hypothesis that in-migration to disaster-affected coastline counties would be higher in the recovery period than in the pre-disaster period. Furthermore, the spatial concentration of the migration system, evident in the test of our first hypothesis, is accompanied by the intensification of flows, especially in-flows. Such churning of migrants suggests that out-migration during the immediate disaster period was mostly temporary and intensified in-migration best characterizes the recovery period.

We illustrate these changes in the in-flow side of the migration system geographically in Fig. 2, panel A, which identifies tied counties for which the number of migrants changed the most between the pre-disaster and recovery periods. Change estimates produced by the multiregional migration model reflect an increase or a decrease in the number of migrants between periods. Counties highlighted in the darkest shade are among the top 5 % of counties that sent increased numbers of migrants to the disaster-affected coastline counties in the recovery period compared with the pre-disaster period. Counties shaded in medium gray are the bottom 5 %, which sent comparatively fewer migrants.⁶ The patterns in the maps are consistent with the increase in in-ties and in-flows from nearby counties and the decrease in in-ties and smaller in-flows from distant counties (Tables 1 and 2). Only a small number of distant counties, largely southern cities, were among the top 5 % of senders in the recovery period. The percentage change was also high in a few distant rural counties, such as places with a strong energy industry (e.g., Sweetwater County, Wyoming, and La Plata County, Colorado) from which migrants pursuing employment in the Gulf Coast’s damaged oil industry may have originated. The majority of tied distant counties—for example, counties composing the metropolitan areas of Boston, Chicago, Denver, New York, and Washington, DC—sent comparatively fewer migrants to disaster-affected coastline counties in the recovery period than in the pre-disaster period. Most top-sending counties are clustered around the Gulf of Mexico in nearby counties or disaster-affected coastline counties. The spatial concentration of top-sending counties is what we would expect if recovery migration were strong and composed of pre-disaster residents of the disaster-affected coastline counties who relocated to nearby counties and new migrants who were connected to employment opportunities in the recovering region.

As a point of comparison, counties receiving the largest increases in out-flows from the disaster-affected coastline counties’ migration system between the two periods are identified in Fig. 2, panel B. Recovery-period out-flows from disaster-affected coastline counties to nearly all tied counties outside the Gulf of Mexico were lower than in the pre-disaster period (medium or light gray). Instead, out-flows from disaster-affected coastline counties concentrated in other disaster-affected coastline counties and nearby counties (dark gray). There are a few exceptions, however, with larger recovery-period in-flows to counties composing the southern metropolitan areas of Miami, Nashville, Oklahoma City, and Shreveport, and a few more distant metropolitan areas such as Boston, Chicago, Denver, Philadelphia, San Diego, and Seattle. Exceptions aside, we see spatial concentration and intensification of out-flows among disaster-affected coastline counties, as predicted, with only a few distant and mostly urban counties becoming important recovery period destinations.

Local spatial concentration is illustrated in Fig. 3, which shows changes in the number of in-flows between disaster-affected coastline counties only. Recovery-period in-migration to disaster-affected metropolitan counties grew, specifically to counties forming the metropolitan areas of Corpus Christi, Houston, New Orleans, and Gulfport. Rural counties hit hardest by Hurricane Rita also received larger in-flows, presumably by returning residents (Jefferson County, Texas, and Cameron and Vermilion Parishes, Louisiana). In contrast, in-migration diminished to rural Texas and Louisiana coastline counties. Although there are fewer in-ties (Table 1) and only very small increases in in-flows (Table 2) among disaster-affected coastline counties, this map demonstrates that in-flows within the recovering region were directed toward urban areas.

Results from the modified gravity regression model buttress our descriptive findings and provide more rigorous support for our contention that the migration system of the disaster-affected coastline counties became more geographically concentrated and movement-intensified in the recovery period by documenting change over time in the average size of migration flows between pairs of sending and receiving counties with controls for population size in sending counties. Positive λ coefficients for the experimental group in all eight models show a statistically significant increase in in-flows to disaster-affected coastline counties between the pre-disaster and recovery periods, whereas negative λ coefficients for the comparison group show a statistically significant decrease in in-flows to nearby counties (Table 3). Model 1 compares the growth in in-flows with disaster-affected and nearby counties from all counties and confirms the second hypothesis by showing that in-flows to disaster-affected coastline counties in the recovery period increased by 0.2 %, while in-flows to nearby counties decreased by 1.1 %. Models 2–4 show that in-flows to disaster-affected coastline counties for all three regions grew in the recovery period, whereas in-flows to nearby coastal counties consistently declined.

The size of the coefficients decreases as the distance from disaster-affected counties increases, thereby supporting the third hypothesis that in-flows will be greater from proximate counties. In-flows from disaster-affected coastline counties and nearby counties grew by 4.0 % and 2.6 %, respectively, between the pre-disaster and recovery periods, and by only 0.2 % for distant counties. The third hypothesis poses that in-flows will be larger from urban counties, for which we find support when each of these models is duplicated for urban counties (Models 5–8). In the urban counties model, the coefficients are larger than in the all counties model, with the exception of the λ coefficient in Model 8, which is equal to its all-counties counterpart in Model 2. Evidence of spatial concentration and intensification of the migration system is demonstrated most clearly in the contrast between in-flows from disaster-affected coastline counties and nearby urban counties, which grew by 5.3 % and nearly 4.9 %, respectively, and in-flow from distant urban counties, which grew by only 0.2 %.

To isolate the disaster recovery effect, we fit a difference-in-difference model. The positive δ coefficients for each difference-in-difference model (Table 4) show that the average change in in-flows to disaster-affected counties was larger than in nearby counties. Moreover, results are consistent with our second and third hypotheses that proximate and urban counties disproportionately attract environmental migrants, and will therefore be a strong source of recovery migration. The coefficient is larger among disaster-affected coastline counties (Model 2) and nearby counties (Model 3) compared with distant counties (Model 4), and is even larger for urban counties (Models 6, 7, and 8).

Our analysis has shown that the migration system of disaster-affected coastline counties became more spatially concentrated in the recovery period by subtracting out-ties with all regions except other disaster-affected coastline counties. Where ties were added to the system, they were mostly in-ties to nearby counties. This spatial concentration was accompanied by growth in the size of in-flows to disaster-affected coastline counties from all regions, especially from nearby and urban counties. Thus, in addition to geographic concentration and intensification, we also see an urbanization of the recovery migration system.

Conclusion

Coastal populations are expected to experience more intense and frequent coastal weather events and inundation resulting from climate change. Rooted in a concern for the human impacts of such environmental events, our study investigated changes in migration systems resulting from Hurricanes Katrina and Rita, two of the most severe hurricanes on record. We examine the effects of these events on U.S. migration systems to gain insights into the migratory consequences of disasters caused by extreme coastal storms, and introduce the concept of recovery migration. Recovery migration occurs when displaced residents of the disaster-affected area return and new in-migrants arrive. We propose that the migration system acts as a conduit for population recovery of disaster-affected regions, thereby moving beyond the immediate post-disaster period in which mobility is best characterized as evacuation and short-term displacement.

Before summarizing the contributions of our research we consider its limitations. The IRS flow data measures only the mobility of taxpayers and their dependents, which excludes the very poor and older populations. This bias may exaggerate mobility rates because excluded groups tend to be less mobile than employed and working-age populations. Additionally, these data do not permit us to analyze the composition of the flows, which may be equally important for recovery as the size of flows. Although we have augmented the IRS flow data by adding measures of county geography and urbanity, refined geographic measures or measures from additional sources could be added to further test hypotheses related to migrants’ destination choice. Finally, more a caveat than a limitation, disaster responses differ across countries with and without private and governmental disaster insurance programs. In the United States, disaster-stricken areas have a comparatively greater capacity to recover than those in nations lacking such programs. Hence, our findings regarding recovery migration are more generalizable to areas with private and public disaster recovery programs. Despite these limitations, we feel confident that our analysis describes the dominant changes in the migration system of the Hurricanes Katrina- and Rita-affected coastline counties between the pre-disaster and recovery periods.

In our investigation, we focus on recovery migration and ask whether an environmental event alters the preexisting migration system of the affected region. In doing so, we make methodological, substantive, and theoretical contributions to research on environment-related migration. Methodologically, instead of focusing on individual and household out-migration from the area affected by an environmental crisis, we use flow data between counties to model the disaster’s broad impact on the complete migration system of the most severely disaster-affected coastline counties. Moreover, by leveraging the unique experimental quality and fine temporal and geographic scale offered by the IRS flow data, we are able to test confirmatory hypotheses about the human impacts of environmental events. Our findings inform the emergent literature on environment-migration relationships by lengthening the time frame and extending the geographic scope of our understanding of this mobility, shifting focus to recovery migration, and moving the study of population-environment interactions into a new and fertile domain.

Substantively, we contribute evidence to research on environmental migration by showing how pre-disaster migration systems channel migration flows in predictable ways after an environmental event. We find that the recovery migration system of disaster-affected coastline counties became more spatially concentrated, including mostly nearby counties, especially urban counties, and only a few urban counties outside the Gulf of Mexico and the South more generally. At the same time, the size of in-flows to disaster-affected coastline counties from these counties grew. This spatial concentration and intensification of in-flows was predicted by the principles of environmental migration destination choice. However, these principles did not anticipate our finding of increased mobility within and between disaster-affected coastline counties and nearby counties. This heightened mobility suggests that migratory churning is part of the recovery as the population adjusts to changed social, economic, political, and environmental structures in a disaster-affected region.

Theoretically, we build on the general principles of environmental migration destination choice (Findlay 2011), principles with obvious roots in general migration theory (Black et al. 2011). These principles propose that out-migration from areas experiencing environmental crises tends to be short-distance and reliant on connections in the pre-crisis migration system, especially ties to urban areas. Our contribution is to consider what occurs after the environmental crisis has subsided and recovery is underway, thereby shifting the focus to in-migration of former residents as well as newcomers—what we call “recovery migration.” By examining recovery migration, we are able to address a pressing question facing places impacted by environmental events: From where will they recover their population? The answer permits residents, planners, politicians, and scholars to understand where households who are likely to migrate to disaster-affected areas are located, information that is useful for a planful recovery. We confirm the general principles of migration destination choice by demonstrating that nearby and urban counties will become the origins of in-migrants to the crisis-affected areas in the recovery period, and extend them by showing that these in-migration streams will be larger in the recovery period if the disaster-affected area is reconstructed.

Regardless of whether the destructive power of Hurricanes Katrina and Rita was compounded by global warming, coastal erosion, or technological failures, their effects on the population of the Gulf of Mexico provide an analog to the potential impacts of climate change (Adamo 2010). By considering analogs for climate change, we can develop more realistic and comprehensive scenarios of how climate change will affect human populations and settlements. Overall, our findings dampen alarmist concerns that climate change will produce drastic population redistribution (e.g., large numbers of poor migrants from the global South flooding the global North) while also providing a set of testable hypotheses to guide empirical research.

Acknowledgments

This research was supported by center Grant #R24 HD047873 and training Grant #T32HD07014 awarded to the Center for Demography and Ecology at the University of Wisconsin at Madison, center Grant #R24 HD041023 awarded to the Minnesota Population Center at the University of Minnesota–Twin Cities by the Eunice Kennedy Shriver National Institute of Child Health and Human Development, and by funds to Curtis from the Wisconsin Agricultural Experimental Station and the Wisconsin Alumni Research Foundation. We thank Lori Hunter and Patricia Romero-Lankao for organizing the CUPC-NCAR Workshop on Migration, Urbanization, and Climate Change held May 7–8, 2012, in Boulder, Colorado, where Curtis and Fussell began this project; the Department of Global Health Systems and Development at the Tulane University School of Public Health, where Fussell was hosted from 2012 to 2013; Lori Hunter and the University of Colorado Population Center for organizing the summer short course, Environmental Demography, that DeWaard attended on June 13–14, 2013; Lori Hunter and the anonymous reviewers at Demography for their insightful feedback; and William R. Buckingham at the Applied Population Laboratory at the University of Wisconsin–Madison for his cartographic expertise.

Notes

References