1 Introduction

Despite significant improvements in the last few decades, minority segregation is still a pervasive aspect of our society today. Indeed, while segregation at the national level is falling, pockets of highly segregated neighborhoods continue to persist (Glaeser and Vigdor 2001; Logan et al. 2004). Segregation manifests in several ways ranging from labor market discrimination to residential segregation. From a historical perspective, residential segregation has received growing attention since (Schelling’s 1969) seminal paper. Several works document the extent to which urban segregation occurs and negatively affects the wellbeing of the segregated groups by limiting their access to education, employment opportunities and health outcomes, along with favoring poverty and criminal behavior (Galster 1987; Orfield and Eaton 1996; Shihadeh and Flynn 1996; Cutler et al. 1999; Williams and Collins 2001; Card and Rothstein 2007; Shertzer et al. 2016; Rothstein 2017; Trounstine 2018; Aaronson et al. 2019; Troesken and Walsh 2019).

Theoretically, segregation outcomes are attributed to economic and social factors. Classical economic arguments going back to Tiebout (1956) and Rosen (1974) state that neighborhood sorting is a driver of residential segregation in the presence of heterogeneity in households’ incomes and preferences, which ultimately determine their willingness to pay for location characteristics. According to this view, the provision of public goods in specific neighborhoods can generate segregation through its effects on prices and housing demand (McGuire 1974; Card et al. 2008; Kollmann et al. 2018). Social arguments originating in Schelling’s (1969, 1971, 1978) works state that social interactions at microeconomic level are a cause of the segregation observed at the aggregate level. According to this view, even moderate individuals’ preferences for ethnic isolation can yield segregation in the entire neighborhood (Zhang 2004a, b, 2011; Pancs and Vriend 2007; Grauwin et al. 2012).

Clearly, these alternative theories emphasize how different policies can be used to mitigate segregation: on the one hand, economic theory suggests that urban planning policy, taking the form of public investments in specific neighborhoods, is a natural answer; on the other hand, social interactions theory identifies in social policy, aiming to promote integration between different ethnic groups, as another solution.Footnote 1 Few attempts have also been made to bridge these two theories by developing a mixed approach in which some social interaction mechanism (i.e., taste for ethnic isolation) is introduced in a classical economic setup of neighborhood sorting (Becker and Murphy 2000; Sethi and Somanathan 2004; Banzhaf and Walsh 2013). By following this last branch of literature, the goal of this paper consists of developing a framework in which economic and social factors jointly contribute to determining single households’ decisions which ultimately drive eventual segregation outcomes. Different to extant works, rather than extending an economic setting to account for some social elements, we introduce economic factors in a simple social interactions framework. Indeed, as pointed out by Arrow (1998), among the possible explanations of segregation, social interactions appear to provide the most convincing argument.

We build on the social interactions theory as outlined in Schelling’s works (1969, 1971, 1978) by introducing some economic mechanisms which determine the extent to which economic and social factors may interact to determine segregation. Schelling’s main conclusion is that even if the majority and minority groups have some degree of preference for integration, individual decisions could lead to aggregate outcomes of predominantly segregated neighborhoods. In particular, following in-migration of a minority group, provided that the minority share exceeds a critical threshold the neighborhood will “tip” to being composed entirely of the minority population.Footnote 2 The existence of such non-linearities and tipping behaviors has been firstly demonstrated empirically by Card et al. (2008) who, by applying a regression discontinuity approach to analyze tract-level data within US cities between 1970–2000, show that the majority population flees when the minority share of a neighborhood exceeds 5–20%. A number of later studies applying the same methodology confirm the existence of tipping over different time periods and in other countries as wellFootnote 3 (Alden et al. 2015; Shertzer and Walsh 2019; Kollmann et al. 2018). However, very little is known about the determinants of such dynamic effects and only a few attempts have been made to characterize tipping. Most of the theories proposed thus far focus on simple static approaches based upon traditional economic arguments (Card et al. 2008; Heal and Kunreuther 2010; Banzhaf and Walsh 2013), while attempts to discuss the role of social interactions in dynamic settings are limited (see Zhang 2011). However, as suggested by Schelling (1971), “the analysis of ‘tipping’ phenomena [...] requires explicit attention to the dynamic relationship between individual behavior and collective results”. Therefore, we aim to contribute to this latter literature in order to provide a dynamic explanation of tipping. Unlike previous works which present some modifications of Schelling’s original model and are particularly complicated to the extent that they preclude the understanding of the main mechanisms in place (see Zhang (2011), who analyzes Schelling’s setup from an evolutionary game point of view), we develop an alternative population dynamics model of heterogenous agents which maintains the basic features of Schelling’s model but gives rise to simpler, intuitive and analytically-tractable solutions.

To this end, we highlight the two main features of Schelling’s model that need to be preserved. First, Schelling is often considered as one the fathers of heterogeneous agents modeling, as his analysis represents one of the earliest agent-based model examples (Epstein and Axtell 1996; Zhang 2004a). Second, in Schelling’s setup there exists some degree of preference for ethnic isolation, which is essential to yield segregation. These two features are maintained in our framework which relies upon social interactions models based upon a random utility setting (Brock and Durlauf 2001; Blume and Durlauf 2003). The presence of heterogeneous agents interacting with one another through a social externality gives rise to preferences for ethnic isolation (due to homophilious considerations), allowing us to understand whether and how social interactions may result in majority and minority segregation driven by individuals’ decisions. In contrast to Schelling’s model and its refinements (Zhang 2011), our framework can be fully and analytically solved to give rise to closed-form results determining not only the long run equilibrium outcomes but also the tipping point, which to the best of our knowledge represents an important novelty in this literature. Specifically, our approach describes the dynamic evolution of the majority share in a neighborhood in which individuals are heterogeneous. This permits us to characterize the conditions under which either segregation or integration will occur, and in the former case, whenever segregation results from tipping, we identify three possible channels through which tipping may occur. Out of these three channels, two of them are associated with equilibrium outcomes (migration flows and changes in tipping points) and one with off equilibrium dynamics (endogenous polarization). To the best of our knowledge, the role of off equilibrium dynamics as a source of tipping has not been discussed in extant literature. Our neat characterization of integration and segregation conditions yields interesting policy implications for social and urban planning policies that mitigate segregation.

Note that our population dynamics setup with heterogenous agents substantially differs from the agent-based frameworks typically employed in extant literature which discuss how the dynamics of complex evolving systems driven by local micro-level interactions in a spatial context may result in the emergence of a macro-level systemic outcome, such as segregation (see Tesfatsion 2006; Namatame and Chen 2016). Indeed, in our baseline model we rely on an a-spatial characterization of a neighborhood which implies that all individuals are located into a single point in space. This allows the network structure of individuals’ interactions to be fully connected such that movements into or out of the neighborhood affect all individuals by the same degree, independent of the number of individuals considered. Such an abstraction from space allows us to derive closed-form solutions which clearly identify how the equilibrium value of the majority share depends on key parameters, which pinpoint mechanisms through which policy may affect the equilibrium outcomes. Nevertheless, we also present extensions of our baseline model showing that even by allowing for a spatial dimension (either in the form of a linear city or a random network) our main conclusions will continue to hold true. Therefore, we believe our modeling choices which permit analytical tractability represent mere simplifying assumptions and do not compromise the model’s ability to characterize real world phenomena.

This paper proceeds as follows. Section 2 briefly reviews the sociological studies on racism discussing how our paper relates to this literature by finding in structural racism the ultimate cause of residential segregation outcomes. Section 3 introduces our a-spatial social interactions framework to describe how single households decide whether to continue residing in a certain neighborhood by taking into account both economic factors (i.e., individual benefits and costs), along with social factors (i.e., the behavior of other own-group households through a social externality). As a matter of expositional simplicity and in line with previous literature, we shall refer to the major ethnic group as “whites” and to the minority as “blacks”. The aggregation of single households’ decisions allows us to characterize the neighborhood outcome in terms of the white share. Section 4 analyzes the asymptotic-population dynamics, which is entirely summarized by a differential equation describing the evolution of the white share. In this setting we can explicitly identify the long run equilibrium outcomes, which may represent either segregation or integration. In the case of minority segregation we determine under which conditions this may result from tipping and we identify two possible mechanisms related to exogenous shocks (migration flows and changes in tipping points) for tipping to occur. Section 5 analyzes the finite population dynamics, in which white flight is stochastic and characterized by some transition probabilities. In such a probabilistic setting, we show that our asymptotic long run predictions may not always come to pass and that a further potential tipping channel, related to endogenous polarization, may emerge. Section 6 presents some extensions of our baseline model to show that its qualitative results apply even in more complicated and general setups, confirming thus that our simple framework represents a good benchmark to characterize the dynamics of residential segregation. Section 7 presents some further extensions of our baseline setup aiming at introducing a spatial characterization, showing that despite the lack of analytical solutions from a qualitative point of view our main results continue to hold true even in a spatial context. Section 8 concludes and suggests directions for future research.

2 Residential segregation and racism

While our paper does not formally model racismFootnote 4 and its conclusions are consistent even in the absence of racism, it is nonetheless critically important to relate our model to the racism literature as residential segregation of ethnic groups may be linked to racial discrimination.

To the best of our knowledge, there have been few attempts in the economics literature (Phelps 1972) to address racism’s contentious issue but none of the works analyze the determinants of racism and its implications on socio-economic outcomes. In contrast, racism and its relation with racial discrimination and racial inequality have been widely discussed in sociology (see Clair and Denis (2015), for a concise and precise literature review).

It is not possible to talk about racism without first defining “race”. The concept of race as applied to humans has no biological meaning, it is merely a social construct. Extensive sociological studies have attributed this social construction to two factors: individual attitudes at the micro-level and an aggregate social process at the macro-level (Clair and Denis 2015). Micro-level racism originates from prejudicial beliefs and attitudes among single individuals. Such prejudices may be either conscious or unconscious (due to an implicit bias, unaware negative beliefs and feelings about racial groups) and affect single individuals’ expectations and actions (Lane et al. 2007). Macro-level racism results from a social process in which individual-level homophily and social norms reinforce each other, influencing aggregate-level attitudes. Such attitudes can give rise, even if without overtly racist intent, to broad-scale effects determining rules, opportunities and policies impacting different ethnic groups in a heterogeneous manner – a phenomenon referred to as structural (or institutional or systemic) racism (Stokely and Hamilton 1967; Denis 2015). These two starkly different views suggest that racism may emanate differently, and thus the strategies to address the problem may require different approaches. For the former, it is an individual-based problem that could be alleviated by favoring cross-racial interactions (Allport 1954). For the latter, this macro-structural process is sensitive to and could be aggravated by increasing cross-racial interactions (Nagel 1995).

Following a significant decline in overt racist attitudes witnessed since World War II, the structural racism theory has become the most widely spread view in contemporary sociological research, even if the individual racism theory remains appropriate for normative analysis.

Indeed, the goal of the individual-level approach to racism is “not to blame the innocent but to improve understanding so that the policies, practices, and ideas that perpetuate racial inequality can be identified and dismantled. Conversely, institutional racism’s focus on the extraindividual might obscure and absolve the role of individual actors in maintaining racism” (Clair and Denis 2015).

Against this sociological debate’s backdrop, our paper bridges the two sociological theories. We endogenize the structural nature of racism by providing micro-foundations through the homophilous behavior of single individuals.Footnote 5 In our model setup, residential segregation may result from single individuals’ decisions regarding staying in or leaving a neighbourhood according to their homophilous preferences for ethnic isolation.Footnote 6 Such preferences affect social norms and determine the desirability of conforming to the same group individuals’ behavior, which in turn affects others’ choices.

In so doing, homophilous considerations and social norms drive social interactions and determine residential outcomes, giving rise to a social process with heterogeneous effects on different ethnic groups at the aggregate level. Residential segregation can thus be attributed primarily to structural racism as individual-level homophily affects residential outcomes only to the extent that it gives rise to social effects at the macro level through the social interaction channel. Such a characterization of residential segregation follows Schelling’s (1969) view according to which micro motives translate into macro outcomes, and homophily is an important driver of individual behavior (Paolillo and Lorenz 2018).

Despite the critiques that it has received over the years (see Jost et al. (2004), for a reviews of such criticisms), homophily is a very well established fact of human nature (even among children), and the primary driver of homophilius relations is ethnicity as “Homophily in race and ethnicity creates the strongest divides in our personal environments” (McPherson and Cook 2001). Even in the housing market’s specific context, some studies show that homophily drives residential choices (Bayer et al. 2014; Krysan and Farley 2002; Bakens and Pryce 2019).Footnote 7 Building on this evidence, our approach recognizes that individual-level homophily is an important driver of human behavior, determining the structure of aggregate and social outcomes. As such, to mitigate segregation successfully, it is essential to focus on the determinants of both individual behavior and social interactions. Our normative analysis, as we are going to see later, confirms this prediction by showing that policies aiming to affect individuals’ incentives and social phenomena can play both an important role in reducing ethnic inequalities.

3 The model

Consistent with previous literature, we consider a setting in which the housing supply is fixed. Specifically, we focus on a neighborhood D (for district) endowed with a large number N of dwellings, each of which can be occupied by either one white or one black household. We assume that both white and black households may leave the neighborhood and when they do so they are replaced by one household of either the same or the opposite ethnicity. In the former case the ethnic composition of the neighborhood does not change and this has no effects on single individuals’ decisions, while in the latter the relative size of the different ethnic groups changes and this impacts individuals’ residential choices. Since we are interested in understanding how population dynamics in the neighborhood may depend on the taste for ethnic isolation, which is clearly relevant only in situations in which the ethnic composition within the neighborhood changes, in the following we will focus only on the implications of the latter scenario in which the departure of a household is replaced by another of a different ethnicity.

For the sake of simplicity, while we model white household’s behavior we assume that blacks’ is exogenously given, and in particular with a given probability 0 < p < 1 they will decide to leave the neighborhood.Footnote 8 Let us assume that initially most dwellings are occupied by whites and each white household i = 1,...,M with M < N attempts to maximize the utility associated with its residential choice. White households are heterogeneous in their degree of preference towards own-group neighbors. The utility function of any white household i is associated with the choice ωi = {0,1} such that ωi = 1 (ωi = 0) denotes that the household leaves (stays in) the neighborhood. The decision to leave or stay is determined by the utility associated with residing in that neighborhood. This utility depends on three elements: a private and a social component, which are common to all households, and an idiosyncratic component, which is household-specific. The private component is given by the net benefit of leaving the neighborhood \(b-c\in \mathbb {R}\), determined by the difference between benefits b > 0 and costs c > 0. The social component is associated with the characteristics of the ethnic composition of the neighborhood. In particular, white households prefer to live in a white neighborhood, thus their propensity to leave increases with the prevalence of black individuals. This is captured by the term \(J \tilde {b_{i}^{e}}\), where \(\tilde {b_{i}^{e}}\) is the expectations of (white) household i about the proportion of blacks in the neighborhood and J > 0 quantifies the importance of the social component relative to the private one in each household’s utility.Footnote 9 Since the share of blacks is complementary to the share of whites, the social component can be rewritten as \(J(1-\tilde {w_{i}^{e}})\) where \(\tilde {w_{i}^{e}}\) is the expectations of household i about the proportion of whites in the neighborhood. Since blacks may move in the neighborhood only when some white dwelling is vacated (i.e., when some white household leaves the neighborhood), the social component represents the fact that whites tend to conform to the behavior of other whites: when one white household expects many whites to leave, that household will more likely decide to leave as well. To properly define \(\tilde {w_{i}^{e}}\), we introduce a vector of state variables \(s_{j}, \ j=1,\dots , N\) counting the number of white dwellings, by taking either the value 0 if the site j is occupied by a black or 1 if it is occupied by a white. This allows us to define the expectations about the proportion of whites as follows: \(\tilde {w_{i}^{e}}=\frac {1}{N-1}\mathbb {E}[{\sum }_{j\neq i}s_{j}]\). Therefore, the decision of the single (white) household depends on the expectations of other (white) households’ choices. The idiosyncratic component is given by εi, which is household-specific random term independently and identically distributed (i.i.d.) across households with zero mean (i.e., E(εi) = 0), determining the specific type of each household. The utility function of each white household is, therefore, given by the following expression:

$$ u_{i}(\omega_{i})=\omega_{i}\left( b-c+\varepsilon_{i} +J(1-\tilde{w}_{i}^{e})\right). $$
(1)

The utility associated with the decision to leave is \(u_{i}(1)=b-c+\varepsilon _{i} +J(1-\tilde {w}_{i}^{e})\) while the utility to stay is ui(0) = 0. Clearly, as long as the utility associated with staying in the neighborhood is larger than the utility associated with leaving, the household will continue to stay in the neighborhood. In other words, whenever ui(1) is larger (smaller) than zero, the white household will leave (stay).

Note that the utility level depends on a number of factors. (i) The private component bc measures the net private utility associated with leaving the neighborhood D. This may be thought of as the net utility obtained from a wide range of sources, including the differential benefits between amenities in other neighborhoods and those in neighborhood D, and the costs related to relocation decisions. The differential benefits may take into account housing prices, population density, degree of safety, location, availability of public transport, school, parks or leisure facilities. The costs may include accommodation search, relocation fees (related to selling their property for homeowners, buying a new property, depositing bonds, or interrupting a current lease for tenants). We do not restrict a priori how benefits and costs compare such that net benefits can be positive, negative or null. Note that the private component term captures economic factors related to the availability of private and public goods, consistent with what is discussed by extant economics literature (McGuire 1974; Card et al. 2008; Kollmann et al. 2018). (ii) The social component \(J(1-\tilde {w}_{i}^{e})\) captures the (negative) externality generated by the presence of individuals of a different ethnicity in the neighborhood on the residential decision of (white) households where J determines the magnitude of such an externality. Specifically, J quantifies the importance of a social norm measuring the extent to which conforming to the behavior of other whites is desirable from the point of view of the single (white) household. Indeed, the larger \((1-\tilde {w_{i}^{e}})\) the higher the expected number of whites leaving by household i, and thus the higher the number of incoming blacks. We can think of the utility arising from this social component as the support and benefits of informal networks in the local community shared by people of the same ethnic group: a larger share of whites leaving will reduce the strength of such informal networks providing each white household with a stronger incentive to leave the neighborhood as well. Therefore, the social component captures social factors similar to those discussed in the social interactions literature (Zhang 2004a, b, 2011; Pancs and Vriend 2007; Grauwin et al. 2012), along with factors capturing how institutions affect individual decisions (Rothstein 2017; Trounstine 2018; Troesken and Walsh 2019),indeed, such a social externality by incentivizing white households to mimic the behavior of other whites determines to some extent their preference for ethnic isolation. Thus, such a social externality determining preference for ethnic isolation may be viewed as a homophily-driven factor leading the majority population to reduce its interactions with minority groups (Durrett and Levin 2005)).Footnote 10 Indeed, the primary driver of homophilius relations is ethnicity, and ethnic considerations have played an important role in the formation of social networks in modern societies (Wimmer and Lewis 2010). (iii) The idiosyncratic component εi captures any form of heterogeneity within the white population, which may be related to either the private or the social components. In particular, it may be useful to think of it as the individual degree of preference towards own-group neighbors due to ethnic prejudice. Similar to what is discussed in the mixed economic-social interactions literature (Becker and Murphy 2000; Sethi and Somanathan 2004; Banzhaf and Walsh 2013), we assume that such a degree of preference towards own-group neighbors is different from household to household. This suggests that, everything else equal, according to their specific type (the specific value of εi) different households may prefer leaving the neighborhood, while others staying. The type of each household, εi, is determined by the realization of i.i.d. random shocks drawn from a common distribution \(\eta (z)=\mathbb P(\epsilon _{i}\leq z)\).

As in Brock and Durlauf (2001), we can show that the utility function (1) can be cast into a probabilistic choice model such that:

$$ \mathbb P(\omega_{i}=1| s_{i}=1,\ \tilde{w}_{i}^{e})=\eta\left( b-c+\varepsilon_{i} +J(1-\tilde{w}_{i}^{e})\right). $$
(2)

By following the same approach as in Blume and Durlauf (2003), we can recast the above model into a dynamic probabilistic choice model. Let us define \({w^{N}_{t}}=\frac {1}{N} {\sum }_{j=1}^{N} s_{j,t}\) as the proportion of white households in the neighborhood at time t, which we will refer to as the “white share” for the sake of expositional simplicity. Under this dynamic version of the model, therefore, the number of whites in the neighborhood changes over time as follows: \(M_{t}=N {w^{N}_{t}}\). Also, note that the white share determines the stock of whites within the neighborhood and needs to be distinguished from the share of white households that leave the neighborhood at time t, which represents what can be referred to as the “white flight share” (see Banzhaf and Walsh, 2016). Then, as in the static model, white households decide whether to leave or stay in the neighborhood at any given time t with the following probability:

$$ \mathbb P(\omega_{i,t+{\Delta} t}=1| \omega_{i,t}=0, s_{i,t}=1, {{w^{N}_{t}}})= \eta\big(b-c+\varepsilon_{i} +J(1- {{w^{N}_{t}}})\big). $$
(3)

Recall that si,t indicates whether the site i is occupied by a black (s = 0) or by a white (s = 1) individual, and when it is occupied by a black the household may decide to leave the neighborhood with probability p ∈ [0,1] as follows:

$$ \mathbb P(\omega_{i,t+{\Delta} t}=1| \omega_{i,t}=0, s_{i,t}=0, {{w^{N}_{t}}})= p. $$
(4)

The Markovian dynamics of the stochastic process \(w^{N}\equiv ({w^{N}_{t}})_{t\geq 0}\) in Eq. 3, representing the evolution of the white share, cannot be explicitly analyzed in a finite dimensional population framework. Therefore, in order to derive some analytical results in the next section we will focus on an infinite dimensional population version of the model, and we will get back to its finite population version in the following section. This will help understand the nature of neighborhood dynamics by clearly identifying when either integration or segregation due to tipping will occur, and in the latter case to pinpoint three alternative tipping channels, distinguishing between those associated with equilibrium and off equilibrium dynamics.

Note that our model is substantially different from Schelling’s (1969, 1971, 1978) original setup but at the same time it maintains its idiosyncrasies, even if in a simplified fashion. In particular Schelling has developed two alternative frameworks to discuss residential segregation, the spatial proximity model and the bounded–neighborhood model (see Zhang (2011), for a concise summary of their characteristics), both of which define the neighborhood from a spatial perspective. Such a spatial structure makes the dynamic analysis of segregation outcomes particularly cumbersome and precludes the possibility of obtaining intuitive closed-form results. In order to increase the degree of tractability, in our model the neighborhood is completely a-spatial, namely it is a point in space where the entire (white and black) population resides.Footnote 11 By analyzing the evolution of the white share we can understand the population dynamics within the neighborhood.

4 Asymptotic dynamics

As mentioned above, even though analyzing the (stochastic) finite population dynamics is not possible, it is possible to analyze the (deterministic) dynamics associated with its asymptotic version in which the number of dwellings N is infinitely large. By following the same argument as in Blume and Durlauf (2003) and Marsiglio and Tolotti (2018), it is possible to derive the following convergence result.

Proposition 1

Suppose that η(⋅) is absolutely continuous. Then, the sequence of stochastic processes wN converges almost surely to w ≡ (wt)t≥ 0, where wt is the solution of the following differential equation:

$$ \dot w_{t}= p (1-w_{t}) - \eta\left( b-c+J (1-w_{t})\right) w_{t}, $$
(5)

with fixed and given initial condition w0 ∈ (0,1).

Proposition 1 states that in such an asymptotic framework we can describe the evolution of the white share, which is entirely determined by two terms: the probability of blacks leaving the neighborhood, p(1 − wt), and the probability of whites leaving it, \(\eta \left (b-c+J (1-w_{t})\right ) w_{t}\), which quantifies the white flight and depends on whites’ net benefits of leaving and the social externality parameter. Clearly, while the former term increases the share of whites (since each black leaving is replaced by a white individual), the latter decreases it.Footnote 12 Equation 5 implicitly determines also the evolution of the share of black households in the neighborhood, bt = 1 − wt, which since complementary to the white share is given by the following equation:

$$ \dot b_{t}= \eta\left( b-c+J b_{t}\right) (1-b_{t})- p b_{t}. $$
(6)

Moreover note that (5) describes the aggregate behavior of whites, and through aggregation we lose track of the behavior of each individual household (and their type). Since the aggregate outcome in the neighborhood is entirely described by a differential equation, this can be explicitly analyzed in order to characterize the long run equilibrium outcomes.

For ease of illustration, in the following we will focus on household types uniformly distributed over the unit interval. This allows us to explicitly derive a number of conclusions, different from previous literature which fails to provide a neat characterization of neighborhood outcomes. However, apart from the eventual loss of closed-form expressions for the equilibrium outcomes, the same qualitative results will hold true also whenever the distribution of types will be single-pickedFootnote 13. Under the uniform distribution assumption, η(z) is given by the following expression:

$$ \eta(z)= \left\{\begin{array}{lll} 0 & \text{if} & z<0\\ z & \text{if} & 0 \leq z < 1\\ 1 & \text{if} & z \geq 1 \end{array}\right. $$
(7)

In such a setting it is straightforward to show that the equilibrium outcome depends on the size of blacks’ leaving probability (p), the size of whites’ social externality parameter (J) and net benefits (bc). If J is small, there always exists a unique stable equilibrium, \(\overline {w}^{*}\), while as soon as J is large enough there might exist either a unique stable equilibrium \(\overline {w}^{*}\) or a multiplicity of equilibria \(\overline {w}_{L}<\overline {w}_{M}<\overline {w}_{H}\) with the middle being unstable and the extremes stable. In order to characterize the possible outcomes, it may be convenient to define the values w+ and w as follows:

$$ \begin{array}{@{}rcl@{}} w^{+}&=&\frac{(b-c+J+p)+\sqrt{(b-c+J+p)^{2}-4Jp}}{2J} \end{array} $$
(8)
$$ \begin{array}{@{}rcl@{}} w^{-}&=&\max\left\{\frac{(b-c+J+p)-\sqrt{(b-c+J+p)^{2}-4Jp}}{2J}; \frac{p}{1+p}\right\}. \end{array} $$
(9)

The above expressions deserve some explanations. Equation 8 represents the larger solution to the quadratic equation resulting from setting \(\dot w=0\) in Eq. 5, in the case in which the realization of w lies within the support of the uniform distribution η. Concerning (9), the first term represents the smaller solution of the same quadratic equation, while the second the lower bound for the share of whites driven by the presence of p in Eq. 5. Indeed, in the extreme case in which \(b-c\to \infty \) the share of whites will naturally converge to \(\frac {p}{1+p}\): even if all whites are keen to (and actually do) leave the neighborhood, a fraction p(1 − wt) of blacks is always exogenously leaving the neighborhood to be replaced by some white entrants. In this respect in the following we will interpret the outcome \(w^{-}=\frac {p}{1+p}\equiv w^{-}_{S}\) as a situation of “minority segregation”, which may be either complete segregation (in the case in which p = 0 such that \(w^{-}_{S}=0\)) or partial segregation (in the case in which p > 0 such that \(w^{-}_{S}>0\)).Footnote 14 Since it is also possible that not all white households choose to leave the neighborhood accepting thus to co-reside with blacks, and this happens whenever \(w^{-}=\frac {(b-c+J+p)-\sqrt {(b-c+J+p)^{2}-4Jp}}{2J}\equiv w^{-}_{I}\), we will interpret such an outcome as a situation of “ethnic integration”, even if whites and blacks may be co-residing in the neighborhood in different (proportionate or disproportionate) amounts. Equation 9 states that minority segregation and ethnic integration are mutually exclusive, and in particular segregation may be possible whenever \({\Delta }\equiv (1+p)\sqrt {(b-c+J+p)^{2}-4Jp}-(1+p)(b-c+J+p)-2Jp>0\) while integration whenever Δ < 0. Therefore, in the case of convergence to the equilibrium w (both in the case of a unique equilibrium or multiple equilibria), we might observe either integration or minority segregation according to which of the above parameter conditions holds true.

The next two propositions summarize all the possible outcomes.

Proposition 2

Assume that the social externality is small (J ≤|bc + p|). If the net benefit is positive (bc > 0) either integration (\(\overline {w}^{*}=w^{-}_{I}\) provided that Δ < 0 holds true) or minority segregation (\(\overline {w}^{*}=w^{-}_{S}\) provided that Δ > 0 holds true) will occur, while if the net benefit is negative (bc < 0) majority segregation will occur (\(\overline {w}^{*}=1\)).

Proposition 3

Assume that the social externality is large (J > |bc + p|). If the net benefit is positive (bc > 0) either integration (\(\overline {w}^{*}=w^{-}_{I}\) provided that Δ < 0 holds true) or minority segregation (\(\overline {w}^{*}=w^{-}_{S}\) provided that Δ > 0 holds true) will occur, while if the net benefit is negative (bc < 0) there exists a tipping point (\(\overline {w}_{M}=w^{+}\)) determining whether majority segregation (\(\overline {w}_{H}=1\)) or either integration (\(\overline {w}^{*}=w^{-}_{I}\) provided that Δ < 0 holds true) or minority segregation (\(\overline {w}^{*}=w^{-}_{S}\) provided that Δ > 0 holds true) will occur.

Proposition 2 and Proposition 3 show that according to the specific parameter configurations the neighborhood might end up in a situation of either segregation (of either the majority of the minority) or integration, and which equilibrium outcome effectively occurs is determined by economic (i.e., net benefit), social (i.e., social externality) and behavioral (i.e., blacks’ leaving probability) factors.

Intuitively, independent of the size of the social externality, if the net benefit from leaving is positive whites will be keen to leave the neighborhood which becomes populated (either entirely or partially) by blacks. Alternatively, if the net benefit is negative, whites will generally not be keen to leave and the neighborhood will remain completely populated by whites unless the initial share of blacks is particularly large, since also in this case the neighborhood will become populated (either entirely or partially) by blacks. If the equilibrium is unique, which outcome is effectively achieved is entirely determined by how economic, social and behavioral factors compare, and the end results are to a large extent intuitive. If the equilibria are multiple, the possible results are more complicated (and interesting) since the outcome depends, other than the relative size of the different parameters, also on the initial share of whites (or, equivalently, of blacks) in the neighborhood. If the initial share of whites is sufficiently large (w0 > w+), then white segregation will always occur. If the initial share of whites is low (w0 < w+), then either integration will occur if Δ < 0 or black segregation will occur if Δ > 0. In the case of integration, the equilibrium share of blacks will increase with whites’ net benefit from leaving and the social externality parameter and decrease with the blacks’ probability of leaving, while in the case of black segregation the equilibrium share of blacks depends inversely on the blacks’ leaving probability.Footnote 15

These different outcomes can be illustrated graphically using a specific parametrization. Specifically, we set blacks’ leaving probability at p = 0.1, we let the social externality parameter vary as J ∈ [0,2] and we consider two different values for the net benefit from leaving, one positive (bc = 0.25 > 0) and one negative (bc = − 0.25 < 0). Figure 1 shows how the equilibrium white share changes with the social externality parameter in the case in which the net benefit from leaving is either positive (panel A) or negative (panel B). We can observe that when the net economic benefit is positive, the equilibrium is unique and decreasing in J, and in particular it is associated with integration until Δ is negative while it is associated with minority segregation as soon as Δ becomes positive. When the net benefit is negative, for small J the equilibrium is unique and corresponds to majority segregation while for large J there coexist two stable equilibria (black curves), one associated with majority segregation and one with either minority segregation or integration, separated by the unstable equilibrium (red curve). Interestingly, there exists a discontinuity in the equilibrium value of the white share, which occurs at the point in which the majority segregation outcome collapses to either the minority segregation or integration outcome.

Fig. 1
figure 1

Equilibrium values of the white share for different values of the social externality parameter, in the case in which the net benefit is either positive (panel A) or negative (panel B)

From the above analysis, we can conclude that tipping can eventually occur only in a situation in which the social externality is large (J > |bc + p|), whites’ net benefit from leaving is negative (bc < 0) and Δ > 0. Specifically, tipping refers to a situation in which the neighborhood starts from a white status-quo (\(w_{0}> \overline {w}_{M}\)) and then it eventually tips to become predominantly black (\(\overline {w}_{L}=w^{-}_{S}\)). This requires that the white share changes basin of attraction by crossing the tipping point (see Zhang, 2011, for a similar discussion of the tipping process). The next proposition summarizes the situations under which this can occur.

Proposition 4

Tipping may occur as a result of either (i) in-migration of the minority, or (ii) an increase in the tipping point.

Proposition 4 outlines two channels through which tipping may occur, and in both cases this can only be due to exogenous shocks. (i) Clearly, an in-migration of blacks is accompanied by an out-migration of whites and thus, by reducing the overall share of whites in the neighborhood, it also increases the white flight. This is consistent with the interwar US experience where migration flows have been a major determinant of residential segregation (Wright 1980; Shertzer and Walsh 2019; Kollmann et al. 2018). (ii) Further, a rise in the tipping point, caused by either an increase in whites’ net benefit from leaving or an increase in the social externality parameter or a decrease in the blacks’ probability of leaving, can directly induce tipping. This is consistent with the post-war US experience where housing price considerations have played a major role in black segregation (Boustan 2010; Brooks 2011; Card et al. 2008). Note that these two channels for tipping are similar to those presented in the economics literature and in the original Schelling’s (1971) work, even if they have been derived in more informal ways or in purely static frameworks. For example, Card et al. (2008) analyzing a local housing market model show that exogenous variations in the growth of blacks and thus in their demand for dwellings may induce whites to flee the neighborhood giving rise to tipping. Similarly, Banzhaf and Walsh (2013) show that in a neighborhood sorting model tipping may arise from exogenous changes in the level of provision of a local public good since it modifies the amenities differential between neighborhoods. Both mechanisms are consistent with our results in Proposition 4 suggesting that exogenous shocks may lead a neighborhood to tip.

The existence of these two channels (i.e., in-migration of blacks and increase in the tipping point) confirms that our theory is supported by empirical evidence which explains real world experiences. However, to the best of our knowledge, no other work has thus far explicitly characterized the tipping point, and its closed-form expression clearly shows that it depends on economic, social and behavioral factors which in turn explain the large variation in tipping point estimates across space and over time (Card et al. 2008; Kollmann et al. 2018) as driven by some parameter differences. Moreover, in the empirical literature on tipping (Card et al. 2008; Alden et al. 2015; Shertzer and Walsh 2019; Kollmann et al. 2018), a common assumption underlying the regression discontinuity approach is that tipping points are exactly the same in all neighborhoods within a city (Banzhaf and Walsh 2013). However, our model suggests this is unlikely to be the case due to differential amenities across neighborhoods, and therefore in certain neighborhoods white flight might not take place till the black share exceeds particularly large values. A similar conclusion is presented in different contexts by Banzhaf and Walsh (2013), and Caetano and Maheshri (2017).

This is the first study to identify the conditions under which integration or tipping may occur in an intuitive and tractable approach. What our analysis suggests is that even a small difference in economic net benefit, social externality or behavioral factors like blacks’ leaving probability can yield diametrically opposite segregation outcomes, which has important policy implications. For example, suppose that we observe a neighborhood in which growing black segregation occurs as a result of tipping, which is the case whenever J > |bc + p|, along with bc < 0 and Δ > 0. Clearly, in this context understanding how to intervene with specific policies in order to mitigate segregation is a concern for local policymakers, and integration may be achieved by simply decreasing the value of Δ enough to turn it negative Δ < 0. In order to effectively achieve integration two different types of policy intervention may be used, aiming at decreasing either the whites’ net benefit or the social externality, respectively. A decrease in the net benefit from leaving for the majority can be implemented through specific urban planning policies aimed at increasing the provision of public goods (like safety, public transport, school, parks or leisure facilities) in the neighborhood, which in turn increases the attractiveness of a dwelling located in the neighborhood D disincentivizing households to leave. A decrease in the social externality can be implemented through social policies aimed at favoring multiculturalism and tolerance with respect to diversity in the local community, in order to lower the concerns of the majority for the presence of minor ethnic groups in the neighborhood. Provided that these policies affect the size of the net benefit and the social externality sufficiently to ensure that the inequality Δ < 0 holds, then it will be possible for the neighborhood to experience integration. Our conclusions are consistent with extant literature which suggests that urban planning and social policies are natural solutions of residential segregation issues (Pancs and Vriend 2007; Banzhaf and Walsh 2013).

Some additional comments may be worthwhile at this point. Our theoretical model suggests that, despite the existence of preferences for ethnic isolation within a neighborhood, the dynamic evolution and equilibrium of its ethnic composition are not obvious at all. As we might expect the taste for own-group neighbors can naturally gives rise to segregation (either of the majority or the minority) but this is only one of the possible outcomes. Indeed, another equally likely possibility is integration suggesting thus that the nature of preferences at individual level per se does not yield a specific residential pattern, but this is rather determined by how economic, social and behavioral factors combine and relate one another. Since integration can be considered as a desirable result which does not require the implementation of specific mitigation policies, in the above (and following) discussion we have emphasized segregation which appears more interesting from a normative perspective as it clearly requires careful consideration to understand how specific policy tools can be employed to improve the residential outcome favoring integration. Nevertheless, this does not mean that integration is not appealing or less likely than segregation, and in fact an important contribution of our paper involves showing that integration and segregation (eventually due to tipping) may alternatively arise within the same framework from the same household preferences structure, and identifying clear cut conditions under which these two alternative scenarios may occur.

5 Finite dynamics

Thus far we have explicitly characterized the possible equilibrium outcomes within a neighborhood, but when studying social dynamic processes it is essential not to rely exclusively on equilibrium analysis but consider also off equilibrium dynamics (Pancs and Vriend 2007). Recall that we have focused only on the deterministic asymptotic version of the stochastic finite population framework characterized by transition probabilities as in Eq. 3, however such an approach may lead to the loss of important information in the presence of multiple equilibria as discussed in Marsiglio and Tolotti (2018).Footnote 16 Indeed, in this case (which occurs whenever J > |bc + p| and bc < 0), due to the probabilistic nature of the model and the presence of two locally stable equilibria, namely \(\overline {w}_{L}\) and \(\overline {w}_{H}\), it may well be possible that, despite our deterministic theory suggests that as long as \(w_{0}> \overline {w}_{M}\) the long run equilibrium will be \(\overline {w}_{H}=1\), the equilibrium achieved will be \(\overline {w}_{L}=w^{-}\) meaning that the neighborhood will endogenously tip (provided that Δ > 0 holds true) giving rise to minority segregation (without any migration flow or change in the tipping point). This feature is related to the metastability properties of the two stable equilibria, whereby metastability is defined as in Benaim and Weibull (2003).Footnote 17

We illustrate this outcome in Fig. 2 which represents the evolution of the white share as predicted by the asymptotic deterministic model (solid curve) and resulting from a simulation of the stochastic finite population model (dotted curve). We arbitrarily set the parameter values as follows: J = 1.3, bc = − 0.25, x0 = 0.806, N = 100, implying that the tipping point is \(\overline {w}_{M}\approx 0.8\), which is consistent with the estimate of Card et al. (2008). However, it is possible to show that the results will be qualitatively identical in any other parametrization giving rise to tipping (see Proposition 3). Note that since \(w_{0}>\overline {w}_{M}\) we would expect, consistent with our asymptotic theory, that white flight will not occur and the neighborhood will become entirely populated by whites, that is \(\overline {w}_{H}=1\). In panel A, we show the typical situation in which the white share converges straight towards the expected equilibrium of majority segregation. In panel B, in contrast, we show that the white share converges towards the unexpected equilibrium of minority segregation, \(\overline {w}_{L}=w^{-}_{S}=0.099\), and such a shift from the basin of attraction of one equilibrium to the other occurs endogenously leading to polarization in agents’ decisions: each agent decides whether to leave by forming his own expectations about the behavior of other agents, and when many agents agree on their expectations of a large number of agents leaving (even if the share of whites is higher than the tipping point) this will give rise to a process of white flight resulting in minority segregation (in which the actual white share falls short of the tipping point). By recalling that in our setup the white share does not necessarily coincide with the white flight share, to gain better understanding on the dynamics underlying the above figure, we report in Fig. 3 the fraction of whites arriving and leaving the neighborhood per unit of timeFootnote 18 associated with the two panels in Fig. 2. In panel A we can observe that convergence to the expected equilibrium \(\overline {w}_{H}=1\) occurs when the share of incoming whites (the black curve) exceeds the white flight share (the red curve). In contrast, in panel B we have the opposite situation in which the red curve is always above the black one, thus resulting in converge to the unexpected equilibrium of minority segregation. We can thus state the following result.

Fig. 2
figure 2

White share dynamics in the asymptotic (in red) and finite population (in black, dashed) model. The simulated white share \({w^{N}_{t}}\) shows different types of behavior: converging straight towards the expected equilibrium (A), or converging towards segregation (B)

Fig. 3
figure 3

Fraction of whites leaving (red, dashed) and arriving (black, dash-dot) per unit of time under the same situations in Fig. 2: converging straight towards the expected equilibrium (A), or converging towards the unexpected equilibrium with segregation (B)

Proposition 5

Tipping may also occur as a result of endogenous polarization in the majority’s decisions.

Proposition 5 outlines a third channel through which tipping may occur, but differently from the two identified in Proposition 4, this arises endogenously from off equilibrium dynamics. Note that endogenous polarization in whites’ decisions is a peculiarity of our framework and deals with its probabilistic properties. As long as the trajectory of the stochastic process is still in its transient phase (i.e., out of equilibrium), it can leave the basin of attraction of the theoretically predicted equilibrium and reach an unexpected equilibrium. In other words, due to the randomness induced by the households-specific random terms \((\varepsilon _{i})_{i=1,\dots , N}\), the dynamics of the model with a finite number N of dwellings may behave differently from that predicted by the asymptotic model. In order to observe this phenomenon, N must be neither too large (otherwise a suitable law of large numbers would push the dynamics quickly towards its deterministic limit) or too small (otherwise pure noisy trajectories would emerge). In our simulations, we consider N in the order of hundreds as a suitable intermediate value to observe such a phenomenon.Footnote 19

To the best of our knowledge, there in no other work in which tipping occurs endogenously except for Zhang’s (2011), in which tipping is driven by a probabilistic framework similar to ours. However, in his model black segregation is the only possible outcome, while in our setting white segregation or integration of blacks and whites may also occur. Moreover, following Schelling’s (1969, 1971, 1978) model his analysis is performed in terms of the pair of different-ethnicity neighbors, while ours is in terms of the white share which is the variable commonly employed in the empirical analysis of segregation (Card et al. 2008). However, note that unlike the other two channels, tipping due to endogenous polarization in whites’ decisions cannot be easily captured in empirical works, and this raises concerns about the robustness of the regression discontinuity approach given that it precludes the possibility of tipping that results from endogenous polarization.

The possibility that a neighborhood will endogenously tip even in the absence of exogenous shocks further stresses the relevance of our previous policy discussion. The need to implement specific policies in order to mitigate segregation does not arise only when growing minority segregation occurs but also whenever there is multiplicity of equilibria such that tipping may occur. In fact, even if the white share is large enough to exceed the tipping point, black segregation may potentially occur via endogenous coordination. In this case relying on urban planning and social policies might be the best strategy in order to effectively mitigate segregation.

6 Extensions

We now discuss a few extensions of our framework in order to show that, despite its simplicity, the conclusions of our baseline model are robust to a number of generalizations, and thus it represents a good benchmark to analyze issues related to residential segregation in an intuitive setup.

6.1 Private component depending on the minority share

We first consider a situation in which the private component in white households’ utility is partly affected by the size of the minority share. Indeed, the private component may capture housing and rental prices which are likely to change when the demographic composition of the neighborhood changes. One of the main driver of the white flight during the postwar period in certain US cities was the white population’s concern about a potential fall in real estate value as a large number of blacks moved from the southern to the northern states (Wright 1980; Shertzer and Walsh 2019; Kollmann et al. 2018). This may be captured in our setup by modifying the private utility component as follows: \(b-c=(b-c)(1+\alpha \tilde {b_{i}^{e}})\), stating that the net benefit from leaving increases with the expected share of blacks in the neighborhood, and the parameter α ≥ 0 quantifies the extent to which the net benefit depends on the demographic composition of the neighborhood. By recalling that \(\tilde {b_{i}^{e}}=1-\tilde {w_{i}^{e}}\), the utility function of white household i can be written as follows:

$$ u_{i}(\omega_{i})=\omega_{i}\left( b-c+\varepsilon_{i} +[J+\alpha(b-c)](1-\tilde{w}_{i}^{e})\right). $$
(10)

By defining \(\overline {J}=J+\alpha (b-c)\), it straightforward to note that the model boils down to our benchmark model in which the social externality parameter is now given by \(\overline {J}\). Therefore, all our conclusions from Proposition 1 to Proposition 5 still hold true (by replacing J with \(\overline {J}\) eventually).

6.2 Endogenous minority dynamics

We now consider a situation in which blacks’ behavior is no longer exogenously given, but it is endogenous and mimics whites’. Blacks, rather than leaving the neighborhood with a given probability p, decide whether to leave or stay in the neighborhood by following a utility maximization process. The utility function of black household k with k = 1,...,MN is given by the following expression:

$$ u_{k}(\omega_{k})=\omega_{k}\left( b^{b}-c^{b}+{\varepsilon_{k}^{b}} +J^{b}(1- \tilde{b}_{k}^{e})\right), $$
(11)

where the parameters with the superscript b denote the fact that blacks’ benefit, cost and social externality parameter may differ from whites’.

By assuming that \({\varepsilon _{k}^{b}}\) are i.i.d. with absolutely continuous distribution function ηb, it is not difficult to show that, under this more general setting, Eq. 5 reads as follows:

$$ \dot w_{t}= \eta^{b}(b^{b}-c^{b}+J^{b} w_{t}) (1-w_{t}) - \eta\left( b-c+J (1-w_{t})\right) w_{t}. $$
(12)

The above equation states that, exactly as in our baseline model, the share of whites in the neighborhood depends on the probability of blacks leaving the neighborhood and the probability of whites leaving it. The only difference is related to the fact that blacks’ leaving probability is given by ηb(⋅) and is determined by blacks’ choices. By assuming an uniform distribution for η, it is clear that we can expect that the behavior of Eq. 12 will closely resemble that of Eq. 5, even if it is computationally more involved, and so our main conclusions from Proposition 4 and Proposition 5 will still hold true. One peculiarity of this more general setting is related to the possibility that complete segregation of the minority may occur as a result of an endogenous optimization problem: if bbcb is negative and very low, then ηb will be zero (corresponding to the p = 0 case in Eq. 5) and thus one of the possible equilibrium outcomes (both in the case of uniqueness and multiplicity) will be \(w^{-}_{S}=0\). In order to get a better sense of the possible outcomes, we present a specific example in which the parameters are set as follows: bc = 0.25, J = 0.5 = Jb, with bbcb free to vary. It is possible to show that for bbcb ≤ 0 the only equilibrium for the white share is \(w^{-}_{S}=0\). For 0 < bbcb < 0.25, the equilibrium \(w^{-}_{S}\) increases with the net benefit in the range [0,0.5], and by symmetry when bbcb = 0.25, the equilibrium is exactly 0.5. When bbcb > bc, the situation is completely reversed and the equilibrium white share is larger than 0.5. In order to allow for multiple equilibria (and hence tipping, eventually), either J or Jb needs to increase. For example, with bbcb = bc = − 0.25 and Jb = J = 1, the two extreme equilibria, \(\overline {w}_{H}=1\) and \(w^{-}_{S}=0\), coexist and are both stable, with their basin of attraction separated by the tipping point.Footnote 20

6.3 Two white subpopulations

We now consider a situation in which the white population is heterogeneous in its taste for ethnic isolation. Some studies show that individuals may differ in their willingness to interact with individuals of a different ethnicity, simply because they may have been differently exposed to diversity (Christ et al. 2014). Individuals interacting frequently with people from diverse ethnic groups tend to develop more tolerance towards diversity than individuals who seldom interact with them. In our setting, this suggests that it may be possible that the majority may be divided into multiple subgroups in which individuals are characterized by a different value of the social externality parameter. In order to keep the analysis as simple as possible, similar to Barucci and Tolotti (2012), let us assume that there exist two subpopulations of whites, W1 and W2, in which the social externality parameter is given by J1 and J2, respectively. Without loss of generality, we assume that J1J2, that is subpopulation 1 is less concerned than subpopulation 2 about the eventual presence of blacks in the neighborhood. We also assume that at time 0, α ∈ (0,1) is the proportion of whites of type 1 and 1 − α the proportion of type 2, such that \(W_{0}=\alpha {W^{1}_{0}}+(1-\alpha ){W^{2}_{0}}\). We generalize the site indicator sj to account for three different types of residents: sj ∈{0,1,2}, in which sj = 0 denotes that site i is occupied by a black, sj = 1 that it is occupied by a white of type 1, and sj = 2 that it is occupied by a white of type 2. When a black household leaves the site (again with exogenous probability p), now it may be replaced by a white of type 1 (with probability 0 < q < 1) or by a white of type 2 (with probability 1 − q). By defining \(w^{1,N}_{t}\) and \(w^{2,N}_{t}\) as the share of whites of type 1 and type 2, respectively, the probabilistic choice model can be stated as follows:

$$ \begin{array}{@{}rcl@{}} \mathbb P\left( \omega_{i,t+{\Delta} t}\vphantom{w^{2,N}_{t}}\right.\! &=&\! 1|\left. \omega_{i,t}=0, s_{i,t}=1, w^{1,N}_{t}, \ w^{2,N}_{t}\right)\\ &=&\! \eta\left( b\! -\! c+\varepsilon_{i} +J^{1}\left( 1\! -\! w^{1,N}_{t}\! -\! \ w^{2,N}_{t}\right)\right) \end{array} $$
(13)
$$ \begin{array}{@{}rcl@{}} \mathbb P\left( \omega_{i,t+{\Delta} t}\vphantom{w^{2,N}_{t}}\right.&=&\! \left. 1| \omega_{i,t}=0, s_{i,t}=2, w^{1,N}_{t}, \ w^{2,N}_{t}\right)\\ &=&\! \eta\left( b-c+\varepsilon_{i} +J^{2} \left( -1w^{1,N}_{t}- w^{2,N}_{t}\right)\right) \end{array} $$
(14)
$$ \begin{array}{@{}rcl@{}} \mathbb{P}\left( \omega_{i,t+{\Delta} t}=1| \omega_{i,t}=0, s_{i,t}=0\right) \! & = &\! p. \end{array} $$
(15)

Similar to what we have discussed in our baseline model, it is possible to prove the following result.

Proposition 6

Suppose that η(⋅) is absolutely continuous. Then, the sequence of stochastic processes \(w^{1,N}_{t}\) and \(w^{2,N}_{t}\) converges almost surely to \(w\equiv ({w^{1}_{t}},{w^{2}_{t}})_{t\geq 0}\), where

$$ \begin{array}{@{}rcl@{}} \dot {w^{1}_{t}}&=& p q (1-{w^{1}_{t}}-{w^{2}_{t}})-\eta\left( b-c+J^{1} (1-{w^{1}_{t}}-{w^{2}_{t}})\right) {w^{1}_{t}}, \end{array} $$
(16)
$$ \begin{array}{@{}rcl@{}} \dot {w^{2}_{t}}&=& p (1-q) (1-{w^{1}_{t}}-{w^{2}_{t}})- \eta\left( b-c+J^{2} (1-{w^{1}_{t}}-{w^{2}_{t}})\right) {w^{2}_{t}}, \end{array} $$
(17)

with fixed and given initial conditions \({w_{0}^{1}}=\alpha w_{0}\) and \({w_{0}^{2}}=(1-\alpha ) w_{0}\), with w0 ∈ (0,1).

Proposition 6 shows that the determinants of the evolution of the share of type 1 and type 2 whites are the same as the determinants of the evolution of white share in the case of a unique white population (economic, social and behavioral factors). The dynamics of the type 1 and type 2 shares is characterized by the difference between inflows in and outflows from the neighborhood of the relevant group. It may thus be possible to observe whites of type 1 fleeing, whites of type 2 fleeing, both white groups fleeing, or none fleeing out of the neighborhood. Again, relying on the assumption of an uniform distribution for η, the behavior of the system (16) and (17) will closely resemble that of Eq. 5, and so our main conclusions from Proposition 4 and Proposition 5 will still hold true.

In order to visualize this, we present a numerical simulation in the case of tipping in Fig. 4. In panels A and B we show the time evolution of three different populations: the (total) white share (in black), the share of type 1 whites (in red) and the share of type 2 (in blue). Parameters are chosen arbitrarily in order to make the figure as clear as possible. In particular, we set bc = − 0.1, J1 = 0.5, J2 = 1, p = 0.22, α = q = 0.5 such that wM ≈ 0.735 and \(w_{L}=w^{-}_{S}\approx 0.51\), and the initial condition, w0 = 0.74, is above the tipping point such that the expected equilibrium will be of majority segregation, wH = 1. In panel A we can observe that the expected equilibrium is effectively achieved, after an initial transient period in which the dynamics seem to depart from the expected outcome. In contrast, in panel B the finite dimensional model converges to the unexpected equilibrium of minority segregation, \(w_{L}=w^{-}_{S}\). These results are qualitatively identical to those discussed for our baseline model, confirming the robustness of our conclusions. Note, however, that the presence of two whites subpopulations allows for more complicated dynamics: by looking at the total share of whites (black curve) in panel A, we can see that the dynamics start decreasing because of a large flight of type 2 whites, but the persistence of type 1 whites makes the dynamics flip behavior and at some point also the type 2 whites return; the end result of this process is represented by the equilibrium of majority segregation in which the neighborhood becomes populated entirely by whites.

Fig. 4
figure 4

White share dynamics for the two subpopulations model both for the asymptotic and finite dimensional model, with N = 200. The simulated white shares show different types of behavior: converging straight towards the expected equilibrium (A), or converging towards segregation (B)

6.4 Vacant dwellings

Finally, we extend our baseline model to allow for the possibility that dwellings remain vacant. Specifically, when an occupant (white or black) decides to leave, the dwelling may be occupied by an occupant of opposite ethnicity (with probability 0 ≤ q ≤ 1) or it remains empty (with probability 1 − q). In turn, a vacant site may become occupied by a newcomer with probability 0 ≤β≤ 1. For the sake of simplicity, we assume that q and β are homogeneous across the ethnic groups and thus take the same values for both whites and blacks. Similar to previous extensions, Eq. 5 can be easily modified to account for this new generalization. For this purpose, we introduce a new variable, ξt = 1 − wtbt, counting the proportion of vacant dwellings such that 0 ≤ ξt ≤ 1. Under the same technical assumption employed in our previous settings, it is possible to prove that the evolution of the share of whites, blacks and vacant dwellings is given by the following expressions, respectively:

$$ \begin{array}{@{}rcl@{}} \dot w_{t}&=& p q b_{t}-\eta\left( b-c+J b_{t}\right) w_{t}+\upbeta \xi_{t}, \end{array} $$
(18)
$$ \begin{array}{@{}rcl@{}} \dot b_{t}&=& -p b_{t} + \eta\left( b-c+J b_{t} \right) q w_{t} + \upbeta \xi_{t}, \end{array} $$
(19)
$$ \begin{array}{@{}rcl@{}} \dot \xi_{t}&=& p (1-q) b_{t} + \eta\left( b-c+J b_{t} \right) (1-q) w_{t} - 2\upbeta \xi_{t}, \end{array} $$
(20)

with fixed and given initial conditions ξ0 ∈ [0,1], w0 = α(1 − ξ0) and b0 = (1 − α)(1 − ξ0), where 0 < α < 1 represents the initial proportion of whites in occupied dwellings. Note that the third equation is redundant since, at each point in time, ξt = 1 − wtbt, thus it is possible to analyze the dynamic of our neighborhood by focusing only on the evolution of the white and black shares, which will in turn determine the evolution of the vacant dwellings. Note also that the model above is a mere generalization of our baseline setup, since whenever q = 1, β = 0, and ξ0 = 0, it boils down to our original model summarized by Eq. 5, thus it should come as no surprise that our main conclusions will also apply in this extended framework.

We now present a numerical simulation to illustrate the implications of introducing the possibility of vacant dwellings in our analysis. In Fig. 5, we show that, similar to our baseline model, multiple equilibria still emerge in the asymptotic population model and thus tipping may occur in its finite population version. In order to push the new vacant dwellings feature as far as possible we set q = 0, meaning that when an occupant leaves a dwelling, this will always remain vacant before getting reoccupied by a newcomer with probability β. We set all parameters as follows: bc = − 0.25, J = 1.3, p = 0.2, q = 0, β = 0.5 and ξ0 = 0.1, and consider two different values for α, quantifying the initial proportion of whites: α = 0.6 (panel A) and α = 0.8 (panel B). In this scenario initially 10% of dwellings are vacant, and the proportion of whites occupying non-vacant dwelling is either 60% or 80%. We can observe that in panel A the initial values of the white, black and vacant dwelling shares are (w0,b0,ξ0) = (0.54,0.36,0.1) and their final values are (wT,bT,ξT) = (0.24,0.54,0.22) giving rise to an integration outcome in which the shares of whites and blacks are both positive. In panel B, instead, their initial and final values are (w0,b0,ξ0) = (0.72,0.17,0.1) and (wT,bT,ξT) = (1,0,0) respectively giving rise to a majority segregation outcome in which the share of blacks is completely null. Given the same parameter values, the convergence to different equilibria when the initial conditions change confirms that multiple locally stable equilibria exist, which is the essential feature (also in our baseline model) for tipping through endogenous polarization to occur. Indeed, minority segregation may occur if the social externality parameter takes a higher value (set at J = 2) while all other parameters remain constant (panel C); we can see that the white share converges to a very low value and the neighborhood becomes predominantly populated by black households, characterizing thus a minority segregation outcome.

Fig. 5
figure 5

Dynamics of the white share (red curve), black share (black dash-dotted curve) and vacant dwellings (blue dashed curve) in the asymptotic model with vacant dwellings. The panels show different types of behavior: integration (panel A), majority segregation (panel B) and minority segregation (panel C)

7 Spatial extensions

One of the most important features of Schelling’s (1969, 1971, 1978) works, which is missing in our baseline setup, is some spatial characterization of neighborhoods. We thus return to our baseline setup summarized by Eq. 5 and discuss further extensions of such an a-spatial model in order to show that, even by introducing some spatial dimension, from a qualitative point of view the main conclusions of our a-spatial analysis still apply. In particular, we model an urban agglomerate by considering two different characterizations of space, one in terms of neighborhoods interconnected within a linear city and one in terms of neighborhoods interconnected within a city through a random network, showing that both setups give rise to results consistent with those outlined in our previous analysis.

7.1 A linear City

We start by extending our previous analysis to allow for the existence of multiple interconnected neighborhoods within a city. Specifically, we now consider an urban agglomerate, which we refer to as a “city”, that is composed of a multiplicity of interconnected neighborhoods. We assume that the city is not a single point in space where all dwellings are located, but it develops along a line (Hotelling 1929) and thus white and black households may be heterogeneously distributed along the city. This implies that in different neighborhoods the initial distribution of the ethnic groups may be different, but individuals may move from one neighborhood to another and thus over time the ethnic distribution in the city might change. Specifically, we denote with Wt,x the number of white households in the position x at date t, in a compact interval \([x_{a},x_{b}]\subset \mathbb {R}\), where, for simplicity, we assume xa = 0 and xb = 1. We assume that there are no migration flows through the borders of [xa,xb], namely the directional derivatives are null, \(\frac {\partial W_{t,x}}{\partial x}=0\), at x = xa and x = xb. In this framework, any position x may be interpreted as a specific neighborhood while a set of adjacent locations is interpreted as a region in the spatial city. Since dwellings are spatially distributed, by defining \(N={\int \limits }_{x_{a}}^{x_{b}}N_{x}\) where Nxdx are the dwellings in neighborhood x, we can compute the share of whites in location x as \(w_{t,x}=\frac {W_{t,x}}{N}\). The behavior of white households in each single neighborhood is exactly as discussed in our a-spatial model, thus we can describe the spatio-temporal evolution of the share of whites through the following partial differential equation:Footnote 21

$$ \begin{array}{@{}rcl@{}} \frac{\partial w_{t,x}}{\partial t}= d \frac{\partial^{2} w_{t,x}}{\partial x^{2}} + p (1-w_{t,x}) - \eta\left( b-c+J (1-w_{t,x})\right) w_{t,x}, \end{array} $$
(21)

where d ≥ 0 is the diffusion coefficient measuring the extent to which individuals tend to move from one to another neighborhood. The term \(d \frac {\partial ^{2} w_{t,x}}{\partial x^{2}}\) captures some form of preference for diversity: white households tend to move from neighborhoods highly populated by whites to neighborhoods less populated by whites. This term tends to counterbalance whites’ preference for residing in neighborhoods predominantly populated by white households. The initial spatial distribution of whites households w0,x captures the eventual initial heterogeneity in the ethnic composition of neighborhoods, and our goal in this section is to understand whether such heterogeneity may tend to persist or die out over time.

Analytically investigating the model’s outcomes would require discussing the solution and properties of partial differential equations, which is quite complicated and goes beyond the scope of this section. We therefore proceed via numerical simulations to illustrate the spatio-temporal dynamics of the white share in our linear city. Figure 6 presents a significant example in which the main parameters are set as in Section 5 (J = 1.3, p = 0.1, bc = − 0.25) while the initial distribution of whites is uniform along the city (w0,x = x) and the diffusion parameter d takes three different values: d = 0 (panel A), d = 0.0001 (panel B) and d = 0.001 (panel C). We can observe that in the absence of spatial diffusion each neighborhood is independent from the others and thus its white share follows exactly the same dynamic pattern we have characterized in our a-spatial theory: neighborhoods in which the initial share is lower than the tipping point (\(\overline {w}_{M}=0.8\)) converge toward a minority segregation outcome, \({w}_{S}^{-} =0.1\), while those in which the initial share is larger than the tipping point the convergence is towards a situation of majority segregation, \(\overline {w}_{H}=1\) (panel A).

Fig. 6
figure 6

Spatio-temporal dynamics of the white share within the linear city for different values of the diffusion coefficient d: d = 0 (panel A), d = 0.0001 (panel B), d = 0.001 (panel C)

In the presence of diffusion, all neighborhoods are interdependent and minority segregation becomes more prevalent, but the final outcome depends on whether spatial diffusion is fast or slow. We can see that there are some neighborhoods which even in the absence of diffusion are meant to converge to a majority segregation outcome, but because of diffusion these neighborhoods end up with minority segregation. In other words, these neighborhoods endogenously tip to become populated by the minority. However, if the diffusion parameter is low, there is still a small region that is predominantly populated by whites in which the equilibrium white share is \(\overline {w}_{H}=1\), while most of the city becomes populated by blacks where the equilibrium is \({w}_{S}^{-} =0.1\) (panel B). If the diffusion parameter is high, such a region disappears and the entire city becomes predominantly populated by blacks giving rise to pervasive minority segregation and the equilibrium white share in every neighborhood is \({w}_{S}^{-} =0.1\) (panel C).

These results can be explained as follows: by capturing preference for ethnic diversity, diffusion tends to increase the concentration of white households in neighborhoods which are initially and relatively less white-populated. If diffusion is fast most whites move rapidly to neighborhoods which are predominantly black-populated, but the evolution of the white share in those neighborhoods follows exactly the same pattern discussed in our a-spatial theory and, since there are only few whites in those neighborhoods, white flight will occur giving rise to an urban outcome of minority segregation. If diffusion is slow a few whites move to neighborhoods that are predominantly black-populated, thus there will still be some predominantly white-population neighborhoods remaining. Therefore, in our setup the preference for ethnic diversity tends to smooth out initial spatial differences across neighborhoods which in the long run leads to more pronounced minority segregation outcomes at urban level. These results suggest that, exactly as in our a-spatial theory, a neighborhood initially populated by whites may tip to become predominantly populated by blacks and identifies preference for ethnic diversity as another possible driver of such an outcome.

7.2 A City as a random network

We now consider a spatial framework where a city forms a network. Each neighborhood (a node on the graph) is characterized by a degree which measures the number of interconnected neighborhoods for a referenced neighborhood. Specifically, different neighborhoods are located on a graph characterized by a certain degree distribution pk, whose average is λ > 0 which measures the average number of connections between neighborhoods. Accordingly, this parameter captures the intensity of the interconnections between individuals across neighborhoods and the extent to which social interactions at urban level matter in driving households’ residential choices. In particular, the probability that a white household will leave a neighborhood of degree k (and thus will be replaced by a black household) depends on the average probability that the city is predominantly black-populated, measured by \(v_{t}={\sum }_{k=0}^{\infty } q^{k} (1-{w^{k}_{t}})\), where qk represents the excess degree distribution related to the distribution pk.Footnote 22 The idea is that all households located in a neighborhood of degree k have the same likelihood of leaving, independent of their exact location. Specifically, the probability that a white household leaves is determined by the degree of the neighborhood, say k, multiplied by the probability that a directly connected peer does, and this latter probability is represented exactly by vt.Footnote 23. Therefore, the higher is k (or vt), the higher is the probability of leaving. Accounting for such urban interconnections, the evolution of the share of whites residing in a neighborhood with degree k is given by the following expression:Footnote 24

$$ \dot {w^{k}_{t}}=p(1-{w^{k}_{t}})-\eta\left( b-c+ J \frac{k}{\lambda}v_{t}\right){w^{k}_{t}}, $$
(22)

where the term \(\frac {k}{\lambda }v_{t}\) multiplied to the social externality parameter generalizes the term (1 − wt) = bt in our a-spatial model to account for the influence of the city ethnic composition in the white households’ decision process. The only difference with respect to our baseline a-spatial model is the “city composition effect”, which depends on k and λ. The higher the degree, k, the higher the probability that individuals’ decisions are influenced by the near urban context. The higher the average degree of the network, λ, the smaller the effect of a single neighborhood on the individual household’s residential decision. Indeed, if all neighborhoods have the same degree (kλ), our extended model in Eq. 22 reduces to our a-spatial model in Eq. 5; in fact, in this case, \(J\frac {k}{\lambda } =J\) and vt = 1 − wt, such that the city is an exact replica of all the identical neighborhoods that compose it. The proportion of neighborhoods populated by whites can be obtained as the weighted average of all wk, that is \(w_{t}={\sum }_{k=0}^{\infty } p^{k} {w^{k}_{t}}\). Note that we now have one equation as Eq. 22 for each k, thus with potentially infinitely many equations it is not possible to explicitly analyze our model’s dynamics. Therefore, since there is no closed-form solution we will proceed with numerical simulations to illustrate the model’s possible outcomes.

Consider the case of a specific distribution function, namely the Poisson distribution with parameter λ.Footnote 25 Figure 7 presents a significant example in which the main parameters are set out as in Section 5 (J = 1.3, p = 0.1, bc = − 0.25) with initial condition \({w^{k}_{0}}=0.81\) above the tipping point \(\overline {w}_{M}=0.8\), and the parameter λ taking two different values: λ = 3 (panel A) and λ = 30 (panel B). The black curves represent the dynamics of \({w^{k}_{t}}\), that is the white share in neighborhoods characterized by different values of k, and for the purpose of exposition, we have only depicted the curves for k ∈{0, 1, 3, 30,100, 300}; the dashed red curve represents the dynamics of wt, the white share in the city, whereas the two dotted blue curves show the dynamics of our a-spatial model with two initial conditions leading either to minority or majority segregation . We can observe that because of the city composition effect some neighborhoods endogenously tip to become predominantly black-populated, while others do not. In particular, neighborhoods characterized by a high value of k (k = 100,k = 300) tend to converge towards a situation of minority segregation in which the white share is low (in our example, \({w}_{S}^{-} =0.1\)), while those characterized by a low degree (k = 0,k = 1) converge towards a situation of majority segregation in which the white share is close to one (remember that \(\overline {w}_{H}=1\)) and those characterized by a degree close to λ show a less extreme behavior and their dynamics is similar to the white share dynamics in the city. Although the dynamics of the white share at urban level is similar in the two simulations, we stress that the proportion of neighborhoods characterized by a low k (and thus they do not experience tipping) have a different impact and outcome in the two simulations: for λ = 3, the proportion of neighborhoods whose white share ends up close to \(\overline {w}_{H}=1\) is about 6% (panel A), while in the case of λ = 30, this proportion is close to zero (panel B). This means that, if the city is less interconnected (small λ), a majority segregation outcome is more likely to occur in some regions in the city, at least in regions which are characterized by poorly interconnected neighborhoods. If the city is more interconnected (high λ), the equilibrium is characterized by pervasive minority segregation in which almost all neighborhoods are predominantly populated by black households.

Fig. 7
figure 7

Dynamics of the white share within the random network for different values of the average degree λ: λ = 3 (panel A), λ = 30 (panel B)

These results can be explained as follows. The degree of interconnections between neighborhoods determines the extent to which the decision of a single household located in a specific neighborhood influences the residential outcome in that neighborhood as well as the other neighborhoods that it is connected to. If the degree of interconnections amongst neighborhoods is high there exist numerous interrelated neighborhoods, and the white flight in one neighborhood will induce a similar white flight in other neighborhoods leading to minority segregation in a substantial part of the city, which could expand to the entire city. If the degree of interconnections is low there exist some neighborhoods which might be almost independent from others (as it could be in the case of elite neighborhoods) and so the white flight in another neighborhood might not have any impact on the residential choices of whites residing in those neighborhoods, which might remain predominantly white-populated. Thus, our setup suggests that tight interconnections between neighborhoods might lead in the long run to more pronounced minority segregation outcomes at urban level. These results suggest that, exactly as in our a-spatial theory, a neighborhood initially populated by whites may tip to become predominantly populated by blacks and identifies interconnections between neighborhoods as another possible driver of such an outcome.

8 Conclusion

There is an established literature arguing that social interactions may be an essential determinant of residential segregation. We build on Schelling’s works (1969, 1971, 1978) by introducing some economic mechanisms in a social interactions framework to determine the extent to which economic and social factors jointly determine segregation. We develop a simple analytically tractable population dynamics model of heterogeneous agents to explicitly characterize and identify the conditions under which integration or minority segregation via tipping might occur as a result of optimizing households’ behavior when choosing whether to leave or stay in a neighborhood. Their decision crucially depends on the relative size of economic (net benefits from leaving), social (social externality) and behavioral (other group’s leaving probability) factors. We explicitly analyze the asymptotic population model to characterize the long run equilibrium and obtain a closed-form expression for the tipping point, which depends on economic, social and behavioral factors. We also analyze the finite population model to identify the possibility of crossing the basin of attraction of the (locally) stable equilibria. This allows us to determine three possible channels, two related to exogenous shocks (migration flows and changes in tipping points) and one to the model’s endogenous probabilistic features (endogenous polarization in the majority’s decisions), through which tipping may occur. The characterization of integration and segregation conditions point to the need for intervention through urban planning or social policies to mitigate segregation. We also present few extensions of our baseline model showing that from a qualitative point of view its main conclusions apply also in more sophisticated and realistic setups. Despite our baseline model abstracts completely from space, we present some further extensions that introduce some characterizations of space (in terms of a linear city and of a random network) which demonstrate that our conclusions are robust even to the incorporation of a spatial dimension.

Despite its simplicity we believe that our model is striking and intuitive in capturing the possible determinants of residential segregation, but further research is needed to clarify the mechanisms underlying tipping. Our closed-form expression for the tipping point has the potential to inform future empirical research, since it provides some testable predictions which permit testing of our model’s implications without relying on the regression discontinuity approach. Moreover, as discussed in Schelling (1971) and in the later social interactions literature, residential decisions often give rise to segregation clusters, meaning that there are some spatial patterns which need to be further examined. Bringing theory to the data in order to validate our conclusions and better explaining the spatial dynamics of residential segregation are left for future research.