This paper highlights a number of important gaps in the UK evidence base on the employment impacts of immigration, namely: (1) the lack of research on the local impacts of immigration – existing studies only estimate the impact for the country as a whole; (2) the absence of long-term estimates – research has focused on relatively short time spans – there are no estimates of the impact over several decades, for example; (3) the tendency to ignore spatial dependence of employment which can bias the results and distort inference – there are no robust spatial econometric estimates we are aware of. We aim to address these shortcomings by creating a unique data set of linked Census geographies spanning five Censuses since 1971. These yield a large enough sample to estimate the local impacts of immigration using a novel spatial panel model which controls for endogenous selection effects arising from migrants being attracted to high-employment areas. We illustrate our approach with an application to London and find that no migrant group has a statistically significant long-term negative effect on employment. EU migrants, however, are found to have a significant positive impact, which may have important implications for the Brexit debate. Our approach opens up a new avenue of inquiry into subnational variations in the impacts of immigration on employment.

Introduction A steady flow of articles from the UK populist press over the past decade have claimed or implied that migrants are taking the jobs of UK-born workers.1 This claim, and the debates surrounding it, have shaped the political agenda on immigration making it one of the defining issues in the Brexit2 referendum. Similar debates have been prominent in other countries that have experienced large inflows of migrants, particularly North America (Borjas, 2017) and Western Europe (Geddes and Scholten, 2016). The usual counter from economists is that such claims tend to fall prey to the ‘lump of labour fallacy’ (Schloss, 1981): the fallacious assumption that there is a fixed amount of work, and hence a fixed number of jobs, in the economy. Under this assumption, a job offered to a migrant worker is necessarily a job opportunity taken away from UK-born workers. The lump of labour assumption is dubious for a number of reasons. First, migrants are also consumers and so a rise in immigration potentially boosts aggregate demand for goods and services, which in turn creates more employment as firms hire more workers to meet the additional demand. Second, economic migrants are often more entrepreneurial than native workers, setting up new businesses and generating new employment opportunities (Levie, 2007). Third, skilled migrants make a disproportionate contribution to innovation (Kerr and Lincoln, 2010) which is likely to improve UK competitiveness, increasing long-run wages and employment (Devlin et al., 2014). Fourth, migrants often fill jobs that UK workers are unable or reluctant to accept, so without those migrants, much of the work would either not be done at all or be done by machines. Fifth, an increase in the share of migrants increases the probability that natives stay in school longer (Hunt, 2017), potentially boosting their long-term employability and productivity. Sixth, migrants increase cultural diversity, which in turn has the potential to boost innovation, social capital, tolerance, overseas trade links and growth (Elias and Paradies, 2016). Seventh, because they tend to be highly mobile and responsive to wage differentials, migrants help ‘grease the wheels of the labour market’ (Borjas, 2001) by responding to higher wages produced by regional labour shortages, improving labour market efficiency which in turn helps foster productivity and growth. Finally, because migrants are typically young and mobile, they can help rebalance the demographic profile of an ageing workforce (Bijak et al., 2007), reducing the dependency ratio, again boosting productivity, competitiveness and long-term employment growth. The extent to which these positive impacts offset the number of jobs taken by migrants is not something that can be predicted by theory alone, as the overall outcome depends on various contextual factors including the mix of skills among migrant and native workers and the types of jobs generated. So, what does the evidence to date tell us? UK empirical studies have consistently shown the impacts of migration on employment to be negligible or zero. For example, after reviewing the evidence to date, the most recent report of the Migration Advisory Committee (MAC) (2018), billed as the ‘most comprehensive-ever analysis of migration to Britain’ (Economist, 20183), concluded that migrants have negligible impact on the employment and unemployment outcomes of the UK-born workers (MAC, 2018: 2). Our contention, however, is that the broad consensus in the empirical literature belies a number of significant weaknesses in the methods used and in the scope of estimates. In particular, we argue that the existing literature has so far failed to provide robust evidence on the local and long-term impacts of immigration and has overlooked spatial spill-over effects between localities. The aims of this paper are: (1) to propose a way of linking data over a much longer time span (half a century) that would facilitate a new generation of research in the UK providing localised longer-term estimates of the impacts of immigration based on large samples; and (2) to develop a way of incorporating both spatial autocorrelation and endogeneity in a spatial dynamic framework. The paper is structured as follows. In the next section we provide a brief review of the literature with a view to identifying key data/methodological deficiencies. The following section describes our approach to creating a linked Census database that has the temporal and spatial attributes needed for robust long-term, large sample and local modelling. We then proceed to set out our strategy for econometric estimation, which we illustrate in the penultimate section with an application to London for the period 1971–2011. The final section concludes with a brief summary of the findings and limitations.

Data linkage Our definition of migrants is based on the country of birth variable from UK Census data. We define a migrant as someone born outside of the UK. While digital UK Census data exist going back to 1971, no two decades have the same definitions for country of birth. They also never use exactly the same geographical boundaries between decades; boundaries used in 1971, 1981 and 1991 in particular are very different from those used in 2001 and 2011. Any time-based analysis requires variables and geographies to both be harmonised: country of birth categories must be consistent and geographical zones must not change between Censuses. A contribution of this paper is to present such a harmonised data set over a five Census period from 1971 to 2011. We focus in this paper on London, as this represents a large labour market area with high population density yielding a large number of aerial units with large samples and relatively high numbers of migrants. Although this paper looks only at London, the harmonised data set is now freely available7 for the whole of Great Britain. It should now be possible, therefore, to apply the model proposed below to other parts of the country. Country of birth data have been harmonised at the lowest level that maximises the number of categories. For example, while later Censuses have many European countries listed, the earliest (1971) has only a single category for Europe. This single category imposes itself on all other decades when matching. Note, however, that while we have linked the data back to 1971, in order to include lagged employment, all the other variables in the model only go as far back as 1981. We use an altered version of 1991 wards as our common geographical zone. This choice was determined by the nature of data in the 1991 Census, where data are presented in two forms: ‘Small Area Statistics’ (SAS) tables are at small geographies but do not contain enough information owing to disclosure restrictions. ‘Local Base Statistics’ (LBS) have more information for country of birth but only at 1991 ward geography level. Choosing this geography as the common basis for the whole data set allows us to maximise country of birth categories across all five Censuses. However, LBS tables also have their own disclosure restrictions where some wards have values set to zero if counts are lower than 1000 people or 320 households. This is solved by creating a new variant of the 1991 ward geography. This takes advantage of the fact that zero-count LBS wards have their populations assigned to neighbouring wards. It is possible to work out which wards these are by comparing with population counts in the SAS tables. SAS geographical zones can be aggregated to wards and their counts subtracted from surrounding wards to detect which contain the re-assigned LBS counts. Once those wards are identified, neighbours are combined into a single new ‘ward’ containing the correct population count. This is only done for a small minority of wards overall but is a necessary step to avoid missing values. Census variables can then be assigned to this new geography. For 1991, borders match precisely, accounting for the new aggregated zones. For the other four Censuses, much smaller geographies are used as the source and so the majority are entirely contained within wards. Others that overlap ward boundaries have their values split according to zone area. The same process is also used for Census employment data, though this is easier than country of birth as there is rather less difficulty in harmonising employment proportions over time. Proof of concept application to London wards It is beyond the scope of the current paper to develop local estimates of the local employment impacts of immigration for the whole of the UK. Rather, we seek to demonstrate proof of concept by applying our proposed method to a single region. We have selected London because it is the pre-eminent destination of migrants in the UK and as such is of interest in its own right: The case of London is worth further study. Immigrant concentration in London as a whole far exceeds that elsewhere in any other city of the UK. Concentration and inflows of immigrants into London also differ widely according to area. (Dustmann et al., 2003: 51) London also has a large number of wards, the basic areal unit of analysis used in our longitudinal linkage of five Censuses, so it guarantees large samples for estimation. Nevertheless, our illustrative application to London should be extendable to other regions of the UK provided they have a sufficient number of wards and sufficient variation in migrant proportions across those wards. This may mean that some regions will need to be clustered in order to achieve sufficiently large samples and variation, but such applications of our method will nevertheless offer for the first time the opportunity to study subnational variation in the impacts of immigration on employment. Descriptive statistics on the data used in the model are given in the supplementary material (available online). Variable selection The advantage of our longitudinally linked ward-level Census data is that they offer both long time spans and the potential for comprehensive geographical coverage. However, they also bring significant limitations, most notably with respect to the choice of explanatory variables. In the modelling strategy described below we seek to explain the role of migration in determining the level of employment in each ward. Our selection of explanatory variables is limited to those we can extract or derive from the Census, namely: migrants born in Ireland, India, Pakistan, Europe and the Rest of the World, the number of UK-born residents and the unemployment rate location quotient (LQ) (explained in the ‘Econometric strategy’ section below).

Econometric strategy Our approach is based on a dynamic spatial panel model developed by Baltagi et al. (2019) to estimate the relationship between the number of people from different countries of birth and the level of employment, controlling for a number of effects. The approach adopted is designed with a view to being able to use the model to simulate different employment outcomes on the basis of different totals of migrants in the future. The estimates below are for a time–space dynamic panel model for i = 1 , … , N where N is the number of districts, in this case N=760, which are the wards of Greater London. Also t = 1 , … , T where T = 5, corresponding to the Census years 1971, 1981, 1991, 2001 and 2011. y it = γ y it − 1 + ρ 1 w i y t + x it β + θ w i y t − 1 + ε it t = 2 , … , T ; i = 1 , … , N (1) In equation (1) y it = ln E it which is the log of the level of employment in ward i at time t. E is defined as the total economically active minus the number unemployed; x it is a (1×K = 7) vector and containing, for ward i at time t, the logs of the levels of migrants born in Ireland, India, Pakistan, Europe and the Rest of the World, together with the log of the number of UK-born residents and the log of the unemployment rate location quotient. The location quotient is defined as the share of the economically active that are unemployed in ward i at time t divided by the share at time t in Greater London as a whole. We have included the spatial lag of the temporal lag w i y t − 1 which helps eliminate bias in the estimation of γ , β and ρ 1 . Baltagi et al. (2019) give more detail of the rationale for its inclusion, based on equilibrium arguments, showing that we would expect to obtain a negative parameter θ relating to this variable. w i is a (1 ×N) vector which corresponds to the ith row of the (N×N ) matrix W N . W N is based on a first order contiguity matrix, so that prior to standardisation w ij = 1 if districts i and j share a boundary and w ij = 0 otherwise. This is subsequently row-standardised so that rows sum to 1. γ is the autoregressive time dependence parameter, ρ 1 is the spatial lag parameter and θ is the time–space diffusion parameter. In (2) we assume that there is a spatial moving average error process, so that: ε it = u it − ρ 2 w i u t (2) which implies that the errors in contiguous districts are interdependent. This local spill-over of unobserved variables and shocks captured by the errors mitigates the impact of omitting spatially lagged regressors ( w i x it ) from equation (1), which typically would be advocated to control for local spill-overs. As pointed out by Pace et al. (2012), Baltagi et al. (2019) and Fingleton et al. (2018), adopting the established convention (Kelejian and Prucha, 1998, 1999) which advises that optimal instruments should include spatial lags of regressors ( x it ) such as ( w i x it ) , the presence of spatially lagged regressors in (1) would require the use of ( w i 2 x it , w i 3 x it ) as instruments, but this appears to result in a weak instrument problem. The innovation u it is a compound process thus u it = μ i + ν it (3) in which the component μ i is a ward-specific time-invariant effect assumed to be iid ( 0 , σ μ 2 ) and ν it is the remainder effect assumed to be iid ( 0 , σ ν 2 ) . μ i and ν it are independent of each other and among themselves. The μ i control for unobserved heterogeneity across wards and the ν it account for random shocks across time and location. Given B N = ( I N − ρ 1 W N ) and C N = ( γ I N + θ W N ) in which I N is an identity matrix of dimension N, we can rewrite equation (1) as y t = B N − 1 C N y t − 1 + B N − 1 x t β + B N − 1 ε t (4) Under this specification, the short-run matrix of partial derivatives is: [ d y d x 1 k … d y d x Nk ] t = β k B N − 1 (5) Equation (5) is a matrix of partial derivatives of y t at time t with respect to the kth explanatory variable, giving the percentage change in employment due to a 1% change in the kth explanatory variable (for example, the number of migrants born in Ireland, etc.). Note that in conventional econometrics this elasticity would be simply the scalar β k , but here we are taking account of spill-over effects, resulting in the (N×N) matrix β k B N − 1 in which the derivative varies according to the ward incurring the change in the kth explanatory variable and the ward in which we measure the response. A simplified average measure of the total effect of a 1% change in the kth explanatory variable in all wards at time t is the total short-run elasticity (tse), which is the mean column sum of β k B N − 1 , thus ts e k = ∑ i = 1 N β k B Nij − 1 N (6) As shown in Fingleton and Szumilo (2019), this is exactly equal to the mean difference between the predicted log employment given by ts e k = ∑ i = 1 N ( y ^ it B − y ^ it A ) / N in which y ^ t A = B ^ N − 1 [ C ^ N y ^ t − 1 + x t β + H ^ N μ ¯ ] y ^ t B = B ^ N − 1 [ C ^ N y ^ t − 1 + ( x t + Δ x kt ) β + H ^ N μ ¯ ] H ^ N = ( I N − ρ ^ 2 W N ) Δ x kt = 1 (7) Matrix H N is defined so that it is non-singular and the time-invariant district heterogeneity effect μ ¯ is based on averaging simulated outcomes of μ = H N − 1 ( B N y t − C N y t − 1 − x t β ) − ν t taken over different realisations of ν ~ N ( 0 , σ ^ ν 2 ) . The total short-run elasticity ts e k gives the percentage change in employment given a temporary, one period, 1% change in variable k across N wards. In contrast the total long-run elasticity tl e k is the percentage change given a permanent 1% change in variable k across N wards. In this case the matrix of derivatives becomes: [ d y d x 1 k … d y d x Nk ] = [ − C + B ] − 1 β k I N (8) Again, the corresponding tl e k is given by the mean difference between the predicted log employment, in this case after iterating y ^ τ A = B ^ N − 1 [ C ^ N y ^ τ − 1 + x t β + H ^ N μ ¯ ] y ^ τ B = B ^ N − 1 [ C ^ N y ^ τ − 1 + ( x t + Δ x k τ ) β + H ^ N μ ¯ ] (9) over τ = 1 , … , T where T is a large number, with Δ x k τ = 1 for all τ for migrant group k. Observe that x t doesn’t change, so the log levels of the number of migrants from each origin is held constant as τ varies, and thus the total long-run elasticity of employment with respect to migrant group k is tl e k = ∑ i = 1 N ( y ^ it B − y ^ it A ) / N (10) Figure 1 illustrates the simulated paths of employment for two arbitrary London wards. Thus we see the paths of ward i (‘01ABFF’), with no change in European migrant numbers, as given by y ^ τ i A , τ = 1 , … , T and with a permanent 1% increase in European migrant numbers, given by y ^ τ i B , τ = 1 , … , T . Also shown are the paths for ward j (‘01ABFR’), given by y ^ τ j A , τ = 1 , … , T and y ^ τ j B , τ = 1 , … , T . We see convergence well before T = 50 and, because of row standardisation, each ward has the same long-run elasticity (equal to the mean of 0.28 given in Table 2) as given by the path differences. Download Open in new tab Download in PowerPoint Below we give the outcome of testing for dynamic stability and stationarity of the model. The rules are: e = vector of eigenvalues of W γ + ( ρ + θ ) e max < 1 if ρ + θ ≥ 0 γ + ( ρ + θ ) e min < 1 if ρ + θ < 0 γ − ( ρ − θ ) e max > − 1 if ρ − θ ≥ 0 γ − ( ρ − θ ) e min > − 1 if ρ − θ < 0 (11) Equivalently, dynamic stability and stationarity requires that the largest characteristic root of B N − 1 C N is <1. Given that these rules are adhered to, the paths of the dependent variable for each ward become stable, converging to levels as given by the prediction equation. Thus, the rules need to be satisfied to allow a long-run elasticity to exist. Further technical details on the rationale for the structural model specification, inference and estimation are presented in the supplementary material (available online).

Illustrative application to London The estimates given in Tables 1 and 2 are for two estimators with corresponding long-run elasticities and indications that we have dynamic stability and stationarity. Two alternative assumptions are made for the moments conditions underpinning the parameter estimates. One is that the regressors are exogenous. This means that the whole temporal sequence of the regressors is independent of the (differenced) errors and hence the dependent (endogenous) variable, log employment level, so that the matrix of instruments includes x t , W N x t , W N 2 x t , t = 1 , … , T . Table 1. Parameter estimates and elasticities assuming exogenous regressors. View larger version Table 2. Parameter estimates and elasticities controlling for selection effects. View larger version In contrast, the endogenous variables y t , W N y t and W N y t − 1 are lagged by two decadal Census periods (i.e. 20 years) to retain zero covariance with the difference errors. Assuming variables are endogenous, it is standard to use only observations that are lagged by two time periods in order to satisfy moments conditions. For example: E ( y i l Δ ν i t ) = 0 , ∀ i , l = 1 , … , T − 2 ; t = 3 , … , T E ( w i y l Δ ν i t ) = 0 , ∀ i , l = 1 , … , T − 2 ; t = 3 , … , T (12) For these to hold, following Arellano and Bond (1991), we require that ν it is serially uncorrelated so that E ( Δ ν it , Δ ν it − 2 ) = 0 , but unfortunately the test statistic m 2 = cov ( Δ ν it , Δ ν it − 2 ) / s . e . is not defined with so few periods. We simply assume that the moments conditions hold by virtue of the length of time between t and t − 2 . The second assumption is that, alternatively, the regressors are themselves endogenous. This seems reasonable in the context, for as Bond (2002) observes, ‘strict exogeneity rules out any feedback from current or past shocks to current values of the variable, which is often not a natural restriction in the context of economic models relating several jointly determined variables’. Accordingly, we prefer to assume that our regressors are endogenous, in other words variation in the regressors both causes, and is caused by, variation in the level of employment. For example, a reasonable proposition is that the number of resident migrants born in Ireland will partly depend on the employment level of the ward. Consequently, we assume feedback from the dependent variable, and hence shocks embodied within the dependent variable, to the regressors and assuming that this is not the case tends to magnify the causal impact of the regressors, as we show subsequently. In order to allow for endogeneity, the regressors are also lagged by two periods, hoping to retain zero covariance as required by the moments conditions. Therefore, the set of instruments only includes x t , W N x t , W N 2 x t , t = 1 , … , T − 2 , and this has the beneficial advantage of reducing the number of instruments from 121 in the case of assuming exogeneity to 51, thus helping to minimise weak instrument problems that tend to occur with a surfeit of instruments. We see the effects of the different estimation techniques in Tables 1 and 2. Note first that we are controlling for temporal and spatial spill-overs. In other words, employment levels tend to have some kind of memory, regardless of the other factors affecting them. The level of employment in a ward is significantly related to the level observed in the previous Census. They also are spatially organised, tending to occur in clumps across space as employment in one district may cause, or be caused by, employment in a nearby, contiguous ward. These are more or less autonomous processes, which we have attempted to isolate so as to obtain the real effect of different country of birth concentrations. Also, some of the heterogeneity across wards, which is assumed to be constant over time, is represented by the term σ μ 2 which denotes the variance of μ i . So, with this error component we pick up the net effect of unobserved factors that make each ward distinctive and which also influence each ward’s employment level. In addition, this is spatially dependent, according to a spatial moving average error process, with the negative coefficient indicating positive local error interdependence, recognising that proximate wards tend to have similar socio-economic and environmental attributes that are omitted as explicit regressors and therefore present in the errors. Additionally, we have controlled for the level of unemployment, or rather the log of the location quotient for unemployment in each Census year. Higher levels of unemployment may be a characteristic of different ethnic groups, so the idea here is to isolate the unemployment effect on the level of employment so as to get a sharper focus on each country of birth group per se, rather than its higher or lower unemployment level. By introducing the different country of birth population levels, one can see if they carry any additional information about the level of employment, over and above that carried by the other variables in the model. So, for example, does knowing the level of Irish-born migrants in a district provide any additional information about the employment level given knowledge of the other variables (unemployment rate location quotient, Indian-born residents etc)? Table 1 gives the parameter estimates and elasticities assuming that the regressors are exogenous. Evidently there are some significant causal impacts, though, as we show below, some of these are illusory. Controlling for the temporal and spatial spill-over effects due to y t , W N y t and W N y t − 1 , evidently the long-run elasticity indicates that a 1% increase in migrants from Ireland leads to a 0.079% fall in the level of employment. The elasticities for Indian-, Pakistani-, European-, Rest of the World- and UK-born residents are all positive. We next consider the outcomes under an assumption that the regressors are endogenous. For example, the statistically significant effects obtained assuming exogeneity may be the results of reverse causation, where an increase in the level of employment causes country of birth numbers to increase, maybe attracted by employment opportunities. For example, Indian-born residents may be sorted into areas with a high level of employment rather than causing a high level of employment. Table 2 gives the details, indicating that allowing for reverse causation, or bidirectional effects, there are no significant changes in local employment levels as a result of change in the levels of Irish-, Indian- and Pakistani-born residents. In other words, the significant negative relationship between Irish-born migrants and employment level, and the positive relation between Indian-born migrants and employment level, given in Table 1, appears to be the outcome of sorting, with Irish migrants attracted to lower employment wards, and Indian migrants attracted to higher employment wards. These different outcomes may be the consequence of social segregation processes and differences in the housing markets as they impact the distribution of these migrant groups. Once we control for sorting or selection effects, as in Table 2, the links between Irish, Pakistani and Indian migrant numbers and employment levels become insignificant, suggesting that the number of Irish or Indian migrants does not cause variation in employment levels. On the other hand, the significant relations between European, UK and Rest of the World residents and employment evident in Table 1 do not disappear after controlling for endogeneity. From Table 2 it appears that there are causal effects whereby a 1% increase in the number of residents born in Europe, the Rest of the World or in the UK leads to rising employment levels. A permanent 1% increase in European-born migrants causes the level of employment to rise by 0.28%. For the UK-born, the impact is a 0.67% increase in employment, and for migrants from the Rest of the World, a 1% increase causes employment to increase by 0.11%. In order to highlight the scope of the methodology, the not insubstantial causal effect of a 1% change in the number of European migrants, and to illustrate possible Brexit-induced impacts, we compare the equilibrium level of employment with the anticipated level if the number of European migrants became 1% lower than the 2011 level in each London ward. Figure 2(a) shows the outcome, which is a variegated pattern of job reduction. The anticipated job loss is about 500 in the financial district of Canary Wharf, with Figure 2(b) illustrating that more than 200 of the 760 wards are predicted to have a job loss of at least 130. Summing over the 760 wards gives an overall total job loss of 117,410 from predicted a total of 4,185,100 London-wide jobs. Of course, this preliminary analysis could be extended to explore the impact of changes in migrant populations in individual or groups of wards and allow different assumptions about other drivers of employment levels. Download Open in new tab Download in PowerPoint

Conclusion This paper has highlighted important deficiencies in the UK evidence base on the employment impacts of immigration. Perhaps most problematic of these is the dearth of robust estimations of the local impacts of immigration – existing studies only estimate the impact for the country as a whole. While the impact of migration on employment and the economy as a whole may be positive, it is possible that the local impacts vary considerably. This potentially raises questions of social justice and whether there is a political imperative for regions that have gained from immigration to compensate areas that have lost out. We also noted that existing studies tend to focus on short- and medium-term effects – we are not able to find any UK studies that provide robust estimates of the employment impact after several decades, for example. This is important as some of the impacts of immigration may take many years to affect employment outcomes. Existing studies also tend to ignore spatial dependence of employment, which can bias the results and distort inference. Our goal has been to address these shortcomings by creating a unique data set of linked Census geographies spanning five Censuses since 1971. These linked data sets yield a large enough sample to estimate the local impacts of immigration using a novel spatial panel model which controls for endogenous selection effects arising from migrants being attracted to high-employment areas. We illustrated our approach with an application to London and found that no migrant group had a statistically significant long-term negative effect on employment. European migrants and those born in the Rest of the World were found to have a significant positive impact. It would be of interest to see whether these findings are replicated in other city regions of the UK. Our approach is not without limitations. Because our focus has very much been on the employment outcomes of immigration, there are a number of important effects we do not consider including hours worked, wages, productivity and the wider economic and social impacts of immigration. Our approach does have the scope to introduce additional covariates, including a more disaggregated breakdown of migrant groups, were data available, and this could challenge the conclusions of our analysis. However, we are aware of no source of data on these variables at the local level over the time span of our study period. There is perhaps an unavoidable trade-off, therefore, between having a richer model (with wages, etc.) for a shorter time period for the UK as a whole and having a more parsimonious model that provides large sample estimates at the local level over a longer time horizon. We argue that in demonstrating how the latter can be achieved we provide an important complementary perspective on migration research, and one that opens up a new avenue of inquiry into subnational variations in the impacts of immigration on employment. Another limitation of our study is that, despite uniquely spanning five Censuses, the number of periods at our disposal is insufficient to formally test the assumptions made regarding the viability of the moments equations used in model estimation. This might be possible given additional periods, but the data set at our disposal currently is at the cutting edge of the data technology: it is probably not feasible to add locally geocoded Census data on the variables in our model before 1971. However, when the 2021 Census data come online, it should be possible to add this extra wave of data to the model, which may make it possible to formally test the moments equations.

Funding

The author(s) disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: This research was funded by the Economic and Social Research Council (ESRC) through the Urban Big Data Centre (Grant Reference: ES/L011921/1) and the Understanding Inequalities (Grant Reference: ES/P009301/1) projects. ORCID iD

Gwilym Pryce https://orcid.org/0000-0002-4380-0388

Notes 1.

For example ‘Immigration is reducing jobs for British workers and David Cameron must act now’, Daily Mail, 11 January 2012; ‘Job hopes of 4 million Brits hit by an “unlimited” pool of EU migrants who are willing to work for low wages’, Daily Mail, 17 May 2018. 2.

‘Brexit’ is the shorthand term used to denote Britain’s exit from the European Union. Analysis of polls has found that one of the main reasons people voted for Brexit was to restrict immigration–see https://blogs.lse.ac.uk/brexit/2018/05/04/leavers-have-a-better-understanding-of-remainers-motivations-than-vice-versa/ (accessed 7 November 2019). 3.

Available at: https://www.economist.com/britain/2018/09/20/what-immigration-system-should-britain-adopt-after-brexit. 4.

Migration Advisory Committee (2012)‘Analysis of the Impacts of Migration’, cited in Devlin et al. (2014). 5.

Devlin et al. (2014). 6.

We define ‘native workers’ as those born in the UK, irrespective of race and ethnicity. 7.

Available at: https://github.com/SheffieldMethodsInstitute/HarmonisedCountryOfBirthDatasets.