Policy Research Working Paper 9500 Do Immigrants Shield the Locals? Exposure to COVID-Related Risks in the European Union Laurent Bossavie Daniel Garrote Sanchez Mattia Makovec Caglar Ozden Social Protection and Jobs Global Practice & Development Research Group December 2020 Policy Research Working Paper 9500 Abstract This paper investigates the relationship between immigra- potential bias and unobservable factors that would lead tion and the exposure of native workers to the health and migrants to self-select into more vulnerable occupations labor-market risks arising from the COVID-19 pandemic. across regions. The results show that the presence of immi- Using various measures of occupational risks based on Euro- grant workers has a causal effect in reducing the exposure of pean Union labor force survey data, the paper finds that native workers to various risks by enabling the native-born immigrant workers, especially those from lower-income workers to move into jobs that could be undertaken from member countries in Eastern Europe or from outside the the safety of their homes or with lower face-to-face interac- EU, face greater exposure than their native-born peers to tions. The effects on the native-born population are more both income and health-shocks related to COVID-19. The pronounced for high-skilled workers than for low-skilled paper also shows that native workers living in regions with workers, and for women than for men. The paper does a higher concentration of immigrants are less exposed to not find a significant effect of immigration on wages and some of the income and health risks associated with the employment—indicating that the effects are mostly driven pandemic. To assess whether this relationship is causal, a by a reallocation of natives from less safe jobs to safer jobs. . Bartik-type shift-share instrument is used to control for This paper is a product of the Social Protection and Jobs Global Practice and the Development Research Group. It is part of a larger effort by the World Bank to provide open access to its research and make a contribution to development policy discussions around the world. Policy Research Working Papers are also posted on the Web at http://www.worldbank.org/prwp. The authors may be contacted at lbossavie@worldbank.org, dgarrotesanchez@worldbank.org, mmakovec@worldbank.org, and cozden@worldbank.org. The Policy Research Working Paper Series disseminates the findings of work in progress to encourage the exchange of ideas about development issues. An objective of the series is to get the findings out quickly, even if the presentations are less than fully polished. The papers carry the names of the authors and should be cited accordingly. The findings, interpretations, and conclusions expressed in this paper are entirely those of the authors. They do not necessarily represent the views of the International Bank for Reconstruction and Development/World Bank and its affiliated organizations, or those of the Executive Directors of the World Bank or the governments they represent. Produced by the Research Support Team Do Immigrants Shield the Locals? Exposure to COVID-Related Risks in the European Union Laurent Bossavie Daniel Garrote Sanchez Mattia Makovec Caglar Ozden World Bank Originally published in the Policy Research Working Paper Series on December 2020. This version is updated on March 2021. To obtain the originally published version, please email prwp@worldbank.org. JEL codes: F22, J28, J61 Keywords: Immigration, COVID-19 pandemic, labor market vulnerability, occupational choice Acknowledgments: The findings, interpretations, and conclusions expressed in this paper are entirely those of the authors. They do not necessarily represent the views of the International Bank of Reconstruction and Development/World Bank and its affiliated organizations, or those of the Executive Directors of the World Bank or the countries they represent. We are also thankful to the Knowledge for Change Program, the World Bank Research Support Budget and the Multi-Donor Trust Fund on Trade for financial support. All errors are our responsibility. 1. Introduction COVID-19 pandemic continues to expose almost everyone in the world to unforeseen health and economic risks. This paper explores the ways in which the prevalence of immigration in a labor market affects different types of workers’ exposure to these risks. We provide evidence that not only the immigrant workers were more exposed to the economic and health-related shocks of the pandemic, but they also served as a protective shield to the native workers. By selecting into higher risk occupations, immigrants enable the native workers to move into jobs that can be undertaken from the safety of their homes or with lower face-to-face interaction with customers and co- workers. The findings come in contrast to the political narrative prevailing in the main high-income destination countries at the beginning of 2020. The United Kingdom completed its withdrawal from the EU on January 31, 2020, in part due to the desire to have more control over its immigration policies and limit migrants’ flows. Intense political debates and polarization on immigration partly fueled the rise of the right-wing parties in Europe and the political controversies over the border wall or the Dream Act in the US. Despite these high-profile examples of the popular and political backlash against immigration, the academic literature provides ample evidence on how immigrant workers often fill the difficult and dangerous jobs that locals are not willing to undertake (Orrenius and Zavodny 2009, 2012; Sparber and Zavodny 2020). As this paper shows, COVID-19 provides another example where immigrant workers are, in effect, “protecting” the native worker by virtue of taking on the jobs that turn out to be the riskiest during the pandemic. Different types of workers sort into different types of jobs according to their abilities, preferences, educational attainment levels, age, and other personal characteristics. These distinctions are especially pronounced between native and immigrant workers. Immigrants face additional sorting mechanisms that define and constrain their occupational choices: they are affected by legal constraints, immigration networks, educational downgrading, limited language skills, and 2 limitations on the recognition of their foreign qualifications. 1 In response to the arrival of immigrants, native-born workers may leave certain geographic areas, occupations or sectors, and enter other markets to better exploit their comparative advantages (World Bank 2018). Such endogenous (and pre-pandemic) location and occupation decisions of both immigrants and natives shape their relative exposure to and the impacts of the labor market and health shocks created by both the pandemic itself and government measures to control its spread. The empirical analysis in this paper accounts for these endogenous decisions while quantifying different types of workers’ health and economic exposure to COVID-related risks. A critical feature of the COVID-19 pandemic is that it generated simultaneous supply and demand shocks that rapidly increased unemployment levels, health risks, and income losses. On the supply side, many sectors faced various disruptions, resulting mainly from mobility restrictions and closures that governments continue to impose to control the spread of the epidemic. Many sectors experienced severe contractions that were most pronounced for occupations that were not deemed essential, and for tasks that could not be performed from home. On the labor-demand side, unemployment and earnings losses have been widespread, due to sudden declines in overall economic activity, especially in sectors related to transportation, tourism, restaurant and food services, direct retail and personal services. Demand growth in sectors such as health care or internet retail have failed to compensate for these losses in other sectors. In addition, workers in many occupations face health risks from direct exposure to potentially infectious people with whom they interact in the course of their work. This is a vital concern for essential occupations that cannot be carried out from home and require extensive and hard-to-avoid face-to-face interactions with customers or coworkers. Many jobs in retail, health and education fall into this category. To assess the exposure of immigrant and native workers to the economic and health risks posed by the pandemic, we construct various measures of vulnerability. We first show that immigrants are generally employed in occupations that are more vulnerable to these COVID-19-related risks. We then examine whether the presence of immigrants in local labor markets impacts the 1 Prominent examples in the literature include studies of network effects (Munshi 2003; Beaman 2012), educational downgrading (Chiswick, Lee and Miller 2005; Dustmann et al. 2013), limited language skills (Borjas 2015; Friedberg 2000), and barriers to recognition of foreign occupational qualifications (Brucker et al., forthcoming). 3 vulnerability of native workers in the same geographic areas. Our empirical analysis is motivated by a general equilibrium model of comparative advantages in task performance between immigrant and native workers (Peri and Sparber, 2009). In the model, native workers reallocate to other occupations in response to an influx of immigrant workers. We use a Bartik-type of instrument to account for the endogeneity of migrants’ location choices that respond to local job opportunities. We focus on destination countries in Western Europe, including the 15 countries that were the initial members of the European Union (prior to the 2004 enlargement), Norway, and Switzerland. This region is the destination for an estimated 60 million of some 272 million immigrants worldwide. The analysis is based on a harmonized labor force dataset that contains detailed information on personal characteristics (such as age, education, occupation and sector) of native workers and labor migrants in hundreds of local labor markets in subregions within European countries. 2 We find that native-born workers in European subregions that have a higher share of immigrants are significantly less likely to be exposed to various dimensions of occupational vulnerability associated with COVID-19. We show that the relationship is causal, and results get stronger once the endogeneity of immigrants’ location choices is taken into account. We also find heterogeneous effects depending on characteristics of native workers. The effects of immigration on job safety are stronger for highly (tertiary) educated native workers, who benefit from the presence of both high-skilled and low-skilled migrants. By contrast, the effects are smaller and statistically insignificant for less (non-tertiary) educated native workers. We also assess whether these compositional effects on employment of certain types of native workers are accompanied by effects on their total employment and wages. We find no evidence of wage or employment impacts among native workers, suggesting that the increase in job safety among native workers is driven by their reallocation from vulnerable jobs to safer jobs. In short, we show that immigration to Western Europe reduced the economic exposure of native-born workers to some of the COVID- 2 We perform the empirical analysis at the NUTS2 (Nomenclature of Territorial Units for Statistics) level. This is a geocode standard adopted by the EU to reference subdivisions of countries for statistical purposes. Each NUTS region includes between 800 thousand and 3 million people and is based on administrative boundaries of countries. Some smaller countries (such as Cyprus, Estonia or Malta) have 1 NUTS2 region while larger countries have a larger number of subregions (such as 39 in Germany and 27 in France). Our dataset has 124 regions in total. 4 19-related labor market shocks by pushing them into occupations that are more amenable to working from home. The literature on the labor market impact of the COVID-19 pandemic has grown exponentially over a very short period. As these papers all show, the impact of COVID-19 on labor markets has been very profound, with asymmetric effects across occupations and population groups (see for example, Adam-Prassl et al. 2020; Beland, Brodeur and Wright, 2020; or Coibilion, Gorodnichenko, and Weber, 2020). Several studies indicate that the employment risks associated with COVID-19 disproportionately affect workers who were already exposed to wage and job uncertainties prior to emergence of the pandemic (Casarico and Lattanzio 2020; Rio Chanona et al., 2020; Mongey, Pilossoph, and Weinberg, 2020; Papanikolaou and Schmidt, 2020; Yasenov, 2020). Similarly, Garrote-Sanchez et al. (2020) show that COVID-19 is likely to exacerbate inequality between workers – especially in richer countries, where better paid and better educated workers with secure jobs are more insulated from the shock. Using monthly data in the US for 2020, Montenovo et al. (2020), Liu and Mai (2020) and Cajner et al. (2020) show occupations that require greater physical proximity and less suitability for teleworking exhibited larger employment losses. By focusing on another category of vulnerable workers – immigrants – our paper contributes to the recent literature on the occupational exposure to COVID-19 related risks. The handful of papers specifically comparing the COVID-19 exposure of immigrant and native workers finds that immigrant workers are typically more exposed to both economic and health shocks associated with the pandemic (Basso et al., 2020; Fasani and Mazza, 2020a and 2020b; Yasenov, 2020). Among the reasons proposed by Borjas and Cassidy (2020) is the higher concentration of migrants in jobs with lower potential to be performed remotely. Yasenov (2020) finds that only 30 percent of the migrants in the United States have jobs suitable for telework, compared to 45 percent of natives. While these papers compare migrants and natives, the impact of immigrants on native workers’ economic and health exposure to the COVID-19-related shocks – the causal linkages – have not yet been fully explored to the best of our knowledge. Although our paper focuses on health and economic vulnerability in the context of the COVID-19 shock, it also contributes to the broad and older literature on the general relationships between 5 immigration and natives’ occupational choices, especially in jobs less prone to negative economic shocks. Studies published before the COVID-19 pandemic suggest that immigration tends to push natives towards better jobs in the US. Peri and Sparber (2009, 2011a) show that immigration pushes natives away from manual occupations into occupations that require more communication- intensive skills, since they have comparative advantage in these areas. Dillender and McInerney (2020) show that Mexican immigration leads US natives to move into safer jobs. Bond, Giuntella, and Lonsky (2020) find that immigration leads US natives to move away from shift work into daytime jobs. Orrenius and Zavodny (2009) and Sparber and Zavodny (2020) show that migration to the US leads less-educated native-born workers to move into jobs with less exposure to hazardous conditions, and that those effects are stronger for female native workers. D’Amuri and Peri (2014) and Foged and Peri (2016) extend the analysis for Europe, again presenting evidence that immigrants, by taking manual occupations, push natives towards more abstract and communication-intensive jobs. Alacevich and Nicodemo (2019) show that inflows of foreign-born workers drive reductions in the injury rate, paid sick leave, and severity of impairment for natives in Italy. In the United Kingdom, Giuntella et al. (2019), find that, in response to migration inflows from 2003 to 2013, medium-skilled, native workers reallocated toward occupations with a lower injury risk index score and lower physical burdens. In Spain, Bellés-Obrero et al. (forthcoming) show that Spanish-born workers are shifting away from manual occupations to those involving more interpersonal interactions. The remainder of the paper is organized as follows: Section 2 describes our conceptual framework. Section 3 discusses the data and methods used to measure economic exposure to COVID-19. Section 4 presents our identification strategy and empirical estimation. Section 5 reports and discusses our empirical results. Section 6 concludes. 2. Conceptual Framework Our analysis builds on the highly influential contribution of Peri and Sparber (2009) who provide a simple and intuitive conceptual framework to explore natives’ occupational choices in response to immigration. Their model is built around constant elasticity of substitution (CES) assumptions 6 under which native-born and immigrant workers are imperfect substitutes within each skill or education category. In turn, these skill categories are imperfect substitutes for each other in the production of goods and services. In this setting, immigrant workers have a comparative advantage in occupations that require manual tasks, as compared to native-born workers who have relatively better communication and language skills. Communication-intensive tasks are defined as activities such as directing, training, and organizing people; “non-communication / manual” tasks refer to physical skills such as dexterity or strength. As a result of the imperfect substitutability between native and foreign workers, increased immigration causes natives to specialize in language- intensive jobs. In addition to labor specialization, wage effects emerge, showing that compensation for manual tasks declines relative to communication tasks, reducing (and possibly fully compensating for) the adverse wage effects of immigration on native-born workers. Native-born workers’ task specialization (in response to immigration) prior to the pandemic can have implications for the exposure of natives and immigrants to labor market risks after the pandemic begins to spread. In other words, ex-ante (pre-COVID-19) labor market forces arising from immigration may have prompted natives to specialize in and move to occupations that also happen to be less economically exposed to ex-post (post-COVID-19) risks that we describe in the data section below. Occupational characteristics that became critical in determining unemployment risks during the COVID-19 pandemic could not have been foreseen prior to the outbreak. 3 However, pre-pandemic occupational choices of migrants are linked to high-risk occupational exposure to the pandemic (Bossavie et al., 2021). As migrants concentrate in high COVID-19-risk occupations, natives’ risk levels would be expected to fall. This would arise, for example, if natives were pushed (or encouraged) to move to higher-skilled occupations that also turned out to be more amenable to home-based work once the COVID-19 shocks hit. We use the relative task content of jobs to construct different types of risks that arise due to the pandemic (Section 3). We then identify the correlation between the relative intensity of manual and communication tasks that different occupations entail with the exposure risks of occupations. 3 As noted by Fasani and Mazza (2020b), while the economic vulnerability measures such as essentiality, temporariness, teleworkability, and industry resilience strongly predict employment outcomes after COVID-19, these characteristics did not have any predictive power prior to the pandemic. In addition, COVID-exposed occupations had a negative wage premium, contrary to the expected compensating wage differential that one would expect if workers in these occupations had realized such increased risks (Basso et al., 2020). 7 This approach allows us to link Peri and Sparber’s conceptual framework directly with measures of labor market vulnerabilities created by COVID-19 related risks. 3. Data and Descriptive Statistics 3.1. Data The primary data source used throughout this paper is the European Labour Force Survey (EU- LFS), which is a harmonized dataset composed of nationally representative labor force surveys conducted in all of the 28 EU countries, in addition to Switzerland and Norway. The sample size of the EU-LFS ranges from 4.5 million to 5.2 million observations depending on the year. Data are available on a yearly basis between 2001 and 2018, the latest available year prior to COVD- 19. The data provide information on whether the individuals in the sample were born in the same country in which they work. While the specific origin countries if migrant workers are not available in the data, broad world regions of origin are observed. These are as follows: the 15 original member states of the European Union (EU15) 4, the 13 new member states of the EU as of 2013 (NMS13) 5, non-EU European countries, Latin America, Asia, the Middle East and North Africa, Sub-Saharan Africa, and other non-European high—income countries that are members of the Organisation for Economic Co-operation and Development (OECD). 6 The dataset includes detailed information on various labor market outcomes of the individuals sampled, along with their demographic characteristics such as gender, age, and educational attainment. In terms of labor market status, observations are labelled as (i) out of the labor force, (ii) seeking employment, or (iii) currently employed. Detailed information on the characteristics of the jobs of employed workers are classified according to (1) the sector of activity at the one- 4 EU15 comprises the higher-income EU countries that were members before accession of the 10 new members in 2004. These are Austria, Belgium, Denmark, Finland, France, Germany, Greece, Ireland, Italy, Luxembourg, Netherlands, Portugal, Spain, Sweden and United Kingdom. 5 The New Member States (NMS13) and their accession dates are the following: Cyprus, Czech Republic, Estonia, Hungary, Latvia, Lithuania, Malta, Poland, Slovak Republic, Slovenia (2004); Romania, Bulgaria (2007); Croatia (2013). 6 These high-income OECD members are Australia, Canada, Japan, New Zealand and the United States. 8 digit disaggregated level of the Statistical Classification of Economic Activities in the European Union (NACE); and (2) by occupational categories up to the three-digit level of the International Standard Classification of Occupations (ISCO-08). Because the ISCO classification scheme changed significantly between 2010 and 2011, 7 we restrict our final sample to the period from 2011 to 2018. In addition, the dataset includes information on the wage decile of the worker in the national wage distribution for individuals who are wage employees. 8 Finally, the dataset provides information for each observation within destination countries at a subnational, regional level: the European Commission’s second level of nomenclature of territorial units for statistics (NUTS2). This allows us to measure immigrants’ presence and labor market outcomes at a rather disaggregated geographical level, and to conduct our estimation at the regional level within European countries. For our analysis, we restrict the sample to the 124 NUTS2 regions in Western Europe, which includes the EU15 countries, Switzerland, and Norway. The share of migrants in the New Member States (NMS13) countries is very small, so the sample captured in the EU Labour Force Survey (EU-LFS) can be noisy. 3.2. Measuring Economic and Health Exposure to COVID-19 and Correlation with the Task Content of Occupations Labor-force surveys are not yet available to provide post-COVID-19, harmonized information of labor-market outcomes for the European Union, at the sub-region level. We categorize economic exposure to COVID-19 using information on workers’ occupations and labor market outcomes prior to the pandemic. This approach reflects the method used in the literature (e.g., Dingel and Neiman, 2020; Garrote Sanchez et. al., 2020) to measure risk exposure before an event takes place. We observe the detailed occupational categories in the data, including the three-digit ISCO-08 categories, which we use to estimate the economic exposure of different types of workers to COVID-19. 7 The scheme changed from ISCO-88 to ISCO-08. 8 The dataset does not include information on the exact wage of individuals. 9 We use four different metrics of vulnerability throughout this paper to assess occupation-specific exposure to the COVID-19 risks. These metrics address the definition of essential workers, the ability to work at home, the degree of income safety, and the extent of health risk. The first three measures capture the effects of supply-side shocks while the fourth is a directly related to exposure to infection. The measures are defined as follows: 1- Essential workers: We define jobs as essential if they are exempt from mobility restrictions even in the presence of government mandated lockdowns. When categorizing jobs as “essential,” we follow Fasani and Mazza (2020a) who use the European Commission’s guidelines concerning the exercise of the free movement of workers during COVID-19 outbreak in the first half of 2020, 9 and supplement it with the Dutch government’s definition of “key workers.” Throughout this uncertain period, governments provided lists of occupations considered essential and thus exempt from mobility restrictions and closures during lockdowns. While the criteria for essential occupations vary somewhat from country to country, there is a significant overlap, such as core occupations delivering critical health care (doctors, nurses and other medical professionals) and providing basic goods and services (food, utilities, transport, security, ICT, or research and science). Based on these criteria, we categorize jobs at the three-digit ISCO-08-classification level (which includes 182 occupations) as either essential or non-essential. According to this definition, about 37.8% of the working-age population in the EU works in occupations that are deemed essential. 10 2- The ability to work from home: Our “amenability to home-based work” metric is based on the premise that a person faces lower employment and income risk if she can work at home, and higher risk if she is required to be physically present in the workplace to carry out the required tasks. Recent literature on the impact of COVID-19 on labor markets highlights the feasibility of home-based work (telework) as a key determinant to overcome governments’ mobility restrictions and associated income and/or unemployment risks. The ability to telework also limits the exposure to health risks among workers that face supply shocks. Dingel and Neiman (2020) use the Occupational Information Network (O*NET) surveys to calculate the amenability of different 9 https://eur-lex.europa.eu/legal-content/EN/TXT/?uri=CELEX:52020XC0330(03) 10 This result is slightly higher than Fasani and Mazza (2020a) who estimate that 33% of workers in the EU work in essential sectors. The disparities emerge as we include the occupations related to “Agriculture, Forestry and Fishery Laborers” (ISCO-08 code 921), as well as an approximation of essential occupations at the 2-digit level for Bulgaria, Poland and Slovenia as the EU-LFS does not include 3-digit occupations for those countries. 10 occupations to telework in the United States, including questions such as the importance of the use of email, the frequency of operating and handling big machinery or processes, or performing physical and outdoor activities. Dingel and Neiman 2020’s measure, however, is based on the Standard Occupational Classification (SOC) system used in the United States at a very granular (6-digit) level of disaggregation. We then concord and aggregate this classification to the International Standard Classification of Occupations (ISCO08) at the 3-digit level, which is widely used in the analysis of EU labor markets. The final measure used in the empirical analysis, a continuous variable between 0 and 1, is the weighted average of occupations in each ISCO 3-digit category that are amenable to telework. As Dingel and Neiman (2020) acknowledge, their Home- Based Work (HBW) index is likely to represent an “upper bound” on the number of jobs that could feasibly be performed entirely from home, as it “neglects” many characteristics that would make working from home difficult. 3- Income risk: This metric characterizes the degree of income safety of the job and is a combined index of our first two measures described above. It takes the value of one if the job is essential (shielded from all mobility restrictions) based on our first metric. If the job is not deemed essential by the government, then it is income-safe only to the extent that it can be performed from home (based on the second metric). So, for non-essential jobs, this measure is equal to the measure of amenability to telework for non-essential jobs. 4- Exposure to health risks: This measure characterizes the health risks involved in the performance of the tasks entailed by an occupation. It quantifies the extent to which a given essential job (one that is exempt from mobility restrictions) involves face-to-face interactions with clients, co-workers or suppliers. When assessing occupation-specific exposure to health risks, we quantify the extent of face-to-face (FTF) interactions with customers, co-workers and suppliers. 11 Our measure of health-safe jobs is a continuous variable between zero and one. It takes the value one if jobs are non-essential. Otherwise, the value varies depending on the amenability to telework (higher levels are associated with higher health safety) and the intensity of face-to-face interactions 11 Jobs with high face-to-face interactions involve tasks such as establishing and maintaining personal relationships, assisting and caring for others, performing for or working directly with the public, or selling to or influencing others through extensive personal interaction (Garrote-Sanchez el al., 2020). 11 (higher levels are associated with lower health safety). To quantify the later, we use the index based on O*NET surveys task descriptions created by Avdiu and Nayyar (2020). The relationships between these four measures and the share of workers in the labor force who fall into each category are shown in Figure 1 below. For example, 38% of the workers in our sample are considered essential (Measure #1 – blue circle). Similarly, 38% of the workers are able to work from home (Measure #2 – red circle). Some 12 % of these teleworkable and essential jobs overlap, which means 64% of the workers in our sample have “income-safe” jobs (Measure #3 – combined area of blue and red circles). We find that 48% of the jobs do not require face-to-face interaction (black circle). An estimated 15% of these jobs are also essential, and 20% are telework jobs. In total, 17% of all jobs are exposed to health risks (Measure #4 – area of the blue circle that does not overlap with either of the red of black circles). These are the essential jobs that require face-to- face interaction but cannot be performed from home via telework. In our sample, 6% of jobs can be considered the “best” jobs with the minimal economic and health risk exposure to COVID-19. That is, they are essential, can be done from home, and do not require face to face interaction. These jobs fall into the overlap of the three circles. By contrast, 17% of the jobs (those outside the three circles) can be considered to have the highest risk level. These jobs are non-essential, cannot be done from home, and require extensive face-to-face interactions. The measures of income safety by occupation are strongly correlated with the relative task content of those occupations, as shown by the literature (Autor, Levy and Murnane, 2003; Autor and Handel, 2013). To calculate the task content of a given occupation, we use the US Department of Labor O*NET abilities survey, which provides information on various characteristics of occupations. We follow the same “basic” definitions of manual and communication skills used in the literature. We thus define manual tasks as those that require “movement and strength,” and we define communication tasks as those that require oral and written expression and comprehension. We then determine the extent to which various jobs entail such tasks by merging these measures with the EU-LFS data and constructing concordances between SOC and ISCO categories of occupations. 12 Figure 1. Relative size of Telework, Essential and Non face-to-face jobs in EU-LFS Telework 6% Essential 12% 17% 6% 14% 9% 17% 19% No face-to-face Total employment Source: Own elaboration based on EU-LFS data (2018) and O*NET surveys (2019). There is a strong and positive correlation between the relative intensity of communication and manual tasks required in occupations and the share of teleworkable jobs (our Measure #2) at the three-digit ISCO level (Figure 2). The Pearson correlation coefficient is 0.87 and statistically significant at the 1% level. As the O*NET survey shows, jobs classified as teleworkable also tend to require a higher intensity of communication-oriented tasks. The correlation between income- safe jobs (teleworkable or essential) and the relative share of communication to manual tasks is also positive, although slightly weaker than for telework jobs alone. 12 This is due to the fact that certain essential jobs, such as low-skilled delivery jobs, require fewer communication tasks. Similarly, there is also a positive correlation 13 between health-safe jobs and the relative share of communication tasks. 12 The correlation is 0.45 but still statistically significant at the 1% level. 13 This correlation coefficient is 0.31 and again significant at the 1% level. 13 Figure 2: Correlation between the relative intensity of communication versus manual tasks and the share of teleworkable jobs in each of ISCO three-digit occupations in 2018 531 232 351 413 142 252 112 234411 215 262 133 334 251 332 121 233 212 412 243 Share of telework jobs in occupation group 111 235 134 231 241 261 242 1 122 216 143 331 431 264 .8 213 352 333 214 .6 342 341 211 263 422 311 265 421 .4 441 132 516 831 612 622 961 315 343 511 432 335 732 .2 226 613 933 754 731 325 221 962 541 522 621 753 832 524 514 834 833 611 711 532 921631 633 721 812 752 912 817 7 8 814 931 723 835 81122 632 15 932818 713 821 9 515 11 751 712 816512 941 813 741 324 742 513 521 322 313314 321 225 224 312 523 141 131 323222 223 952 0 1 1.5 2 2.5 3 3.5 Relative intensity of communication / manual tasks Source: Own elaboration based on EU-LFS data and O*NET surveys. Note: The size of the circles represents the relative size of each ISCO three-digit occupation in employment. The correlation uses the relative size of each occupation as weights. 3.3 Migrant-Native Gaps in Income and Health Risk Exposure to COVID-19 Before illustrating the econometric analysis, we present several descriptive statistics of our sample of working age individuals (ages 25-64) from the EU Labour Force Survey. As of 2018, there were 270 million working-age people residing in Western European countries. Almost one out of five (18.7%) of them were foreign-born (Table 1). These immigrants originally come from various high-income Western European countries (3.1%), new EU member states such as Poland and Romania (3.6%), other European countries that are not EU members (3.5%), and other locations outside of Europe (8.1%). Migrant workers are slightly younger than native workers on average (mean age is 43 versus 45.4) and have a slightly lower share of men than the native workforce (48.4% versus 50.2%). Migrants also have slightly lower education levels, with 31.2% having 14 tertiary education degrees compared to 33.5% of natives. There is, however, significant variation with respect to education levels depending on migrants’ region of origin. The share of workers with tertiary education is higher among immigrant workers from other EU15 countries (43%) compared to those from non-EU European countries (21.2%) and new EU member states (27.3%). Female migrants to the EU also tend to be more educated than male migrants (32.3% versus 30.0%), especially among migrants from non-EU European countries and new member states. Table 2 displays the employment outcomes of native and immigrant workers, by education level, including the type of occupations in which they are employed, based on the classifications described in the previous subsection. The employment rate of immigrant workers in the EU is significantly lower than that of natives, both for high-skilled and low-skilled individuals. Among immigrant workers, 70.7% are employed, 9.0% are unemployed, and the remaining 20.3% are out of the labor force. The corresponding figures for native-born workers are: 77.6% employed, 5.4% unemployed and 17% out of the labor force. Among the tertiary educated, employment levels are 79.2% for immigrants and 87.5% for native workers. Employment figures decline significantly for less-educated groups, to 67% for immigrants and 72.6% for natives. The main contributor for the employment gap is the lower labor force participation rate of immigrant women (see Table 3). Pre-COVID occupational decisions of migrants or natives were not shaped by concerns about economic vulnerability in a future pandemic since COVID-19 was a true shock in the econometric sense. Yet, we still observe significant and systematic differences in economic exposure to COVID-19 related risks. Among the employed, we observe significant differences in the types of jobs immigrant and native workers have. For example, 27.3% of immigrant workers are employed in teleworkable jobs compared to 40.9% of native-born workers. This difference is also observed when we compare natives and migrants within the same education group. Immigrant workers, however, are more likely to be employed in essential jobs (41.8%) compared to natives (37.6%). Although the share of essential jobs varies across destination countries, migrants are over- represented in essential jobs in almost every host country in the EU15. Essential jobs tend to be concentrated among highest-and lowest-skilled occupations (such as physicians and grocery store employees) both for natives and immigrants. As seen in Figure 3, immigrant workers in Western European countries represent a very large share of essential workers in the health sector, specifically among doctors and nurses, as well as ICT. However, the share of immigrant workers 15 is even larger in essential jobs that also have the lowest wages and the lowest levels of education. We find that 40.4% of migrants with low levels of education have essential jobs, compared to 34.1% of native workers. Finally, immigrants are less likely to be employed in income-safe jobs (53.9% versus 65.3% of native workers) and in health-safe jobs (77.1% versus 84.4) though the difference between the immigrants and the natives is narrower in the latter (84.4% versus 77.1%). Figure 3: Share of foreign workers in key essential occupations by destination country Doctors Nurses ICT UK Personal care Drivers Cleaners All occupations Doctors Nurses ICT Germany Personal care Drivers Cleaners All occupations Doctors Nurses ICT France Personal care Drivers Cleaners All occupations Doctors Nurses ICT Spain Personal care Drivers Cleaners All occupations Doctors Nurses ICT Italy Personal care Drivers Cleaners All occupations 0 .1 .2 .3 .4 .5 Share of foreign workers Source: Own calculations based on EU-LFS 2018 data. 16 Figure 4: Share of workers by region of origin and risk type Panel A. Telework jobs Panel B. Essential jobs 50% 50% 45% 40% 45% 35% 30% 40% 25% 20% 35% 15% 10% 30% 5% 0% 25% Natives Foreign EU15 NMS13 Non-EU28 Natives Foreign EU15 NMS13 Non-EU28 born migrant migrant migrant born migrant migrant migrant Panel C. Income risk jobs Panel D. Health risk jobs 50% 50% 48% 45% 46% 40% 44% 35% 42% 30% 40% 25% 38% 20% 36% 15% 34% 10% 32% 5% 30% 0% Natives Foreign EU15 NMS13 Non-EU28 Natives Foreign EU15 NMS13 Non-EU28 born migrant migrant migrant born migrant migrant migrant Source: Own calculation based on EU-LFS 2018 data, following EC directive (2020) and Fasani and Mazza (2020). Figure 4 highlights the heterogeneity in exposure to labor market risks between immigrants from different regions of the world, especially between migrants from high-income EU15 countries relative to the rest of immigrants. For telework jobs, the share of immigrant workers from EU15 countries is even higher than natives’ share, whereas it is significantly lower among migrants from 17 NMS13 and non-EU countries (Panel A). While the share of immigrant workers in essential jobs from NMS13 countries and the rest of the world is significantly higher than it is for native workers, the share of migrants from other EU15 countries in those jobs is very similar to the share for natives (Panel B). Likewise, the share of EU15 migrants exposed to income risks is lower than the share of natives, while it is significantly higher for migrants from other parts of Europe or the rest of the world (Panel C). Regarding health risks, migrants from EU15 countries are exposed to the same level of risk as natives, while migrants from other regions are significantly more exposed (Panel D). We observe significant differences in exposure to economic and health-related risks associated with COVID-19 between male and female migrants. As shown in Table 3, women, both native- born and immigrants, are more often employed in teleworkable jobs compared to men. This pattern is mainly driven by women with low levels of education who are more likely to be in teleworkable clerical or office-based occupations; by contrast men of similar education levels tend to be employed in manufacturing-oriented or menial jobs. The shares of male and female workers with high levels of education in teleworkable jobs are about the same, regardless of whether they are immigrants or native-born. The share of male and female native workers in essential jobs is about the same (37%); however, the share of female migrants in essential jobs (47.8%) is much higher than male migrants (36.9%). This pattern is again driven by less-educated workers. We see that 50.4% of female migrants with low levels of education are in essential jobs, compared to 32.8% of male migrants with low levels of education. Female migrants are also more likely to be in income-safe jobs (68.6%) compared to male migrants (51.6%). A similar pattern is observed among native workers, although the gap is less pronounced. These gender patterns are again driven by low-skilled migrants. Our data show that 63.4% of female migrants with low levels of education are employed in income-safe jobs compared to only 40.7% of male migrants with similarly low levels of education. On the other hand, female migrants are much less likely to be employed in health-safe jobs (66.7%) compared to male immigrants (85.7%) or female natives (79.7%). Income levels also show variation by origin and education levels. Unfortunately, the survey only reports the income decile of a worker (within a country) rather than an exact income level. Hence, 18 our analysis is less precise. Table 4 reports the share of workers in the top three deciles of the national wage distribution, by a worker’s migration status and the type of occupation. Migrant workers are less likely to be in the top three wage deciles (24.4%) compared to natives (35%). This income gap between natives and immigrants is observed for workers with both high and low levels of education. For example, 47.1% of immigrant workers with tertiary educated are in the top three deciles, as opposed to 57.4% of native workers with a similar education level. The highest wages are observed among workers in teleworkable jobs. We find that 47.7% of all workers with teleworkable jobs are in the top three deciles. This ratio is almost identical between migrant and native workers. For other types of jobs, however, the wage penalty of immigrant workers is much more salient. An estimated 37.2% of natives in essential jobs are in the top three deciles, compared to 23.2% of migrants in the same types of occupations. For income-safe jobs, the proportion is 40.3% for natives compared to 28.7% for immigrant workers. The pattern is similar for health-safe jobs. Significant wage gaps are also observed between male and female workers, regardless of education level and migration status. As shown in Table 5, 47% of male native workers are in the top three income deciles, compared to 22.4% of native women. Among migrants, 32.6% of men and only 14.8% of women are in the top three deciles. These significant gender wage gaps exist among workers (regardless of education or migration status) in teleworkable jobs, essential jobs, and income-safe jobs. For example, 55.8% of native males in income-safe jobs are in the top three deciles, compared to 26.6% of native women. Among immigrants, 41.1% of males in income-safe jobs are in the top three deciles, compared to 18% for immigrant women. 4. Empirical Strategy and Identification We now investigate whether immigration to high-income Western European countries had a causal impact on the extent of exposure of native workers to COVID-19-related economic and health risks via the jobs they hold. As discussed above, our empirical specification to study this question is derived from the conceptual framework of Peri and Sparber (2009). A critical issue to identify 19 causal effects in this context is that the employment of migrants in specific geographic regions or in specific occupations is likely to be endogenous, driven by market conditions, expectations, social networks, and historical precedents. This section discusses the empirical strategy used to address such selection biases. We first construct measures of different types of risks (at the NUTS2 level) to which native workers are exposed, based on their occupations. Since the analysis is at the regional level, we first must control for differences in terms of the distribution of workers by their education, gender, and age. This step is necessary because those personal characteristics affect the supply of skills required for different types of occupations, and, at the same time, they may be correlated with the size of immigration in the region. To control for this spurious correlation, we first run individual- level regressions to predict how occupation- and labor-market-related characteristics for each individual i depend on individual observable characteristics. We perform these regressions separately for each of the years of study (2011 to 2018) and by immigration status (native vs. migrant): = + ∗ + (4) Status stands for each of the main six variables capturing natives’ and migrants’ status in the labor market. The first four are occupation-level measures defined in the previous section: essential jobs (Measure # 1); amenability to telework (Measure #2); degree of income-safety (Measure #3) or health-safety (Measure #4). The last two are dummy variables based on the employment status (employed or unemployed) or the income level (in the top three income deciles) of the individual. Xi is a vector of individual-level characteristics such as gender, education, and dummies for age groups. Using this measure, we calculate the gap between the actual measure and the predicted � one ( ) for each individual based on the estimated parameters: � _ = − (5) � where � = [ | ] = � . This gap is a measure of each of main six labor- � + market outcomes of interest “cleaned” of differences in individual level characteristics. We then calculate the weighted mean of the variable g_Statusi across individuals in each NUTS2 region using their personal weights in the sample, and the hours worked. 20 The obtained variable, now denoted g_Statusst for each regions s and year t, becomes the dependent variable in our estimation framework, which we regress on the fraction of migrant workers in the working-age population in region s and year t (fst ), together with region and time fixed effects. The baseline equation we estimate is given therefore by the following: _ = + + + (6) In this specification, g_Statusst is the adjusted measure of risk exposure among the working-age population of natives in NUTS2 region s and year t; τt are the time fixed effects that account for common, time-varying technological parameters; αs are region fixed effects that account for variation due to time-invariant unobserved characteristics of the population in that region; and fst is the share of foreign workers in each region and year. Therefore, γ is the main parameter of interest in the estimation. From the EU Labour Force Survey, our sample consists of 123 NUTS regions in 14 EU countries observed with yearly frequency from 2011 to 2018. Estimating Equation 6 using a standard OLS approach can raise endogeneity concerns, even after including fixed effects to control for unobserved regional characteristics and time patterns. In particular, OLS estimates for γ may be inconsistent if the share of immigrants in EU regions is correlated with the unobservable term in equation 6 ( . . [εst |fst ] ≠ 0). This might arise, for instance, when migration flows to a given region are driven by unexpected contemporaneous region-specific shocks. Further, OLS estimates could suffer from reverse causality if migrants choose to locate in areas where natives are employed in occupations that complement migrants’ skills. These occupations, in turn, may differ in their degree of exposure to economic and health risks due to COVID-19. Similarly, unobserved technology and demand shocks, which may differ across EU regions due to sectoral variation, may simultaneously affect the productivity of (and demand for) certain occupations as well as the demand for immigrant workers. A large body of literature indeed shows that migrants’ location choices are endogenous with respect to local labor market conditions, including native workers’ characteristics and labor market outcomes (World Bank, 2018). Therefore, OLS estimates of the relationship between the vulnerability of the jobs that natives hold to COVID-19 and the degree of presence of immigrants in the region may not reflect true causal effects. 21 We use an instrumental variables (IV) approach to address the possible endogeneity of the share of migrants in a given destination region. We follow a consolidated approach in the literature, and construct a shift-share Bartik instrument to control for the current migrants’ share in a given region. 14 The instrument consists of both a “share” and a “shift” component. To build the instrument, we use the spatial distributions of immigrants in 2002 across European NUTS regions (the initial “share” prior to the period of analysis), interacted with the aggregate growth rate of the stock of immigrants at the EU level between 2002 and the year of study (the “shift”). 15 The “share” component provides predictive power to the instrument. It relies on the fact that new immigrants’ settlement decisions are partly determined by the presence of earlier immigrants, mainly through network externalities (Gross and Schmitt, 2003). Recent migrants learn of opportunities in those locations from the network of immigrants, and they may prefer the benefits of living near their co- nationals. 16 As a result, the past distribution of immigrants across European regions can help predict their current location. The “shift” component provides exogeneity to the instrument if the events influencing the flow of migrants to the EU as a whole are arguably exogenous to changes to the labor market conditions in each NUTS region. These criteria will be satisfied as long as each region is small relative to the EU and business cycles of the regions are not highly correlated. Formally, our first Bartik instrument– denoted PS, hereafter – for the share of immigrants in each region can be expressed as the following: �1+ � _ = ℎ_,2002 ∗ [ (1+ ,2002 ] (7) ),2002 where Share_migs,2002 is the share of migrants in NUTS2 destination region s in 2002; (1+g_mig)t,2002 is the growth rate of working-age migrants in EU15 countries between 2002 and year t; and (1+g_nat)t,2002 is the growth rate of working-age natives in all EU15 countries between 2002 and year t. 14 For prior studies using Bartik-type instruments to study the impact of migration shocks see Jaeger (2006); Edo et al. (2019), Mitaritonna, Orefice, and Peri (2017), Peri and Sparber (2009). Jaeger, Ruist, and Stuhler (2018) provide a detailed review of this literature and offer potential caveats about the validity of this instrument. 15 The earliest year for which information on workers’ country of origin is available in the EU-LFS is 2002. Our instrument, however, does not distinguish the lagged immigrant shares by migrants’ region of origin as some destination countries in the EU-LFS do not provide the region of origin breakdown. 16 A large body of literature documents the role played by migrants’ networks and diasporas is affecting the location choices of new migrants. See for example, Munshi and Rosenzweig (2006), Beine, Docquier and Ozden (2011), Beaman (2012), Patel and Vella (2013) and Kramarz and Nordström Skans (2014), among others. 22 One potential concern with this baseline version of the Bartik instrument is that it uses the contemporary number of native individuals in a given European region. As a result, the instrument may capture spurious effects due to the potentially endogenous native population growth in a destination region – which could be due, for example, to internal labor mobility of native workers in response to local demand shocks. To address this concern, we construct a second version of the Bartik instrument in which we fix the native population at the initial year (2002, in our case) as suggested in Mitaritonna, Orefice, and Peri (2017). We first define: = � ,2002 � ∗ , ,2002 where Mst is the number of migrants in time t in destination NUTS2 region s; Ms,2002 is the number of migrants in destination region s in 2002; Meu,2002 is the number of migrants in the EU15 in 2002; and Meu,t is the number of migrants in the EU15 in year t. Then the instrument for the share of migrants Pct_migst based on Mitaritonna, Orefice, and Peri (2017) – denoted MOP instrument, hereafter – is: _ = � � (8) +,2002 where Ns,2002 is the number of natives in destination region s in 2002. The Bartik instrument is valid under the condition that the past migration networks (which create differential exposure to common immigration shocks) do not have a direct effect on natives’ occupational choices – such as, in our case, the allocation between jobs that are economically exposed or not exposed to COVID-19 (Goldsmith-Pinkham, Sorkin and Swift, 2020). In other words, the distribution of immigrants across regions of Europe in 2002 should be orthogonal to subsequent labor market shocks in the region, conditional on region and year fixed effects. Therefore, our identification strategy posits that the reasons to co-locate where previous immigrants lived are driven by preference and information networks, and arguably not by region- specific labor market shocks. A standard concern in the literature is that the initial distribution of immigrants is correlated with persistent local factors that also influences natives’ allocation across occupations. As discussed 23 in Dustmann et al. (2005), Basso and Peri (2015) and Edo et al. (2019), a way to minimize the correlation between past immigration and current outcomes is to use a sufficiently long time lag to predict the actual number of immigrants. To mitigate this concern, we derive initial migration shares using data from 2002 – the first year for which data are available in the EU Labour Force Survey, and benchmark time that is nine years before the beginning of our period of analysis. In addition, we control for region fixed effects that take care of unobservable regional characteristics in the IV specification, to ensure that the exclusion restriction is met. Another possible concern with the shift-share instrument is that using an earlier distribution pattern of immigrants as an instrument can violate the exclusion restriction in the presence of general equilibrium adjustments and strong serial correlation of immigrant flows (Jaeger et al., 2018). If it takes time for markets to adjust to shocks, the error term can also include other influences that reflect ongoing general equilibrium adjustment effects of past migration shocks. This is a greater concern in contexts where the mix and share of immigrants does not vary much over time, like in the US, and where, as a result, the first stage of the IV estimation is stronger due to more stable migration inflows. As noted by Jaeger et al. (2018), the serial correlation in the number of immigrants is lower in Europe compared to the United States as immigration pattern to Europe have changed substantially over time. 17 In particular, the period between our initial measure of immigrant share in 2002 and our period of study for labor market outcomes, which starts in 2011, has been marked by the EU enlargement which triggered important changes in migration flows to the EU (Elsner, 2013). A final concern is that the total inflows of immigrants in the EU (the shift) could be correlated with local demand-pull factors, given that the inflows include the immigration growth rates in the NUTS2 region of reference. To mitigate this issue, prior studies have adjusted the measure of the exogenous supply-push factor by using migration flows to destination areas other than the one of reference (Ajzenman, Dominguez, and Undurraga, 2020). As a robustness check, we adjust the PS and MOP instruments including just the growth rate of immigration to all EU NUTS2 regions besides the one of analysis in the “shift” component. Our main results are robust to this alternative 17 In section 5, we evidence that the first stage of the IV estimation is weaker than what has been evidenced in the US context, once region fixed effects are accounted for. 24 instrument, although the instrument becomes slightly weaker (as measured by the F-stat of excluded instrument). 18 5. Results 5.1 Impact of Immigration on Occupational Exposure to COVID-19. We now report our results on the impact of the extent of immigration across EU regions on several dimensions of occupational exposure to the COVID-19 pandemic, based on the data and empirical approach detailed in Section 3. The first indicator of occupational safety with respect to COVID- 19 exposure is the share of natives employed in teleworkable jobs (i.e., those jobs that can be carried out from home – our Measure #1). Table 6 displays the naïve OLS estimates for the effects of the immigrant share on this measure. We observe a large, statistically significant, and positive association between the share of immigrants in EU regions and the share of natives employed in teleworkable jobs. The magnitude of the coefficient implies that an increase of 1 percentage point in the share of immigrant workers in a given region increases the share of natives employed in teleworkable jobs by around 0.25 percentage points. The magnitude of the coefficient is reduced to 0.13 when region fixed effects are also included (column 3), but the estimate remains highly statistically significant. To address the endogeneity concerns outlined in Section 4 regarding the dependent variable in the naïve OLS estimation, we present our IV estimates of the effect of migrants on natives’ share in teleworkable jobs (Table 7). The first stage results of the IV estimation are reported in the top panel. The F-statistic is larger than the Stock and Yogo (2005) cutoff of 10 in the specification without fixed effects (columns 1 and 4), and with both time effects and country fixed effects (columns 2 and 5). The instrument remains strong in the most restrictive specification that includes both time and NUTS2-level-region fixed effects (columns 3 and 6), and clusters the standard errors at the NUTS2 level, despite the relatively short time period of study. The results suggest that, 18 Results with these alternative instruments are available upon request. 25 according to the criteria of Stock and Yogo (2005), both versions of the instrument are powerful enough to explain variation of the share of migrants across regions, and to identify the treatment effect of interest. In all these specifications, the association between the share of immigrants and our network-based instrument is statistically significant at the 1 percent level, although the coefficient estimate is slightly reduced once we control for NUTS2 fixed effects (columns 3 and 6). The bottom panel of columns 1 through 6 of Table 7 report the second-stage results of our IV estimation, using the two alternative instruments. As in the OLS estimation, the share of immigrants in the region has a positive impact on the share of natives employed in jobs that are amenable to telework once the endogeneity of the migrant share is accounted for. The IV coefficient estimates for γ remain statistically significant at the 1 percent level and are slightly higher than the OLS estimates when we compare columns 1 and 4 in Table 7 to column 1 in Table 6 (or columns 2 and 5 in Table 7 to column 2 in Table 6). The comparison between IV and OLS estimates suggests that OLS estimates are likely have downward bias, underestimating the true effect of migrants’ presence on natives’ take-up of teleworkable jobs. One possible explanation for this downward bias is that immigrants locate in areas where there are more manual, blue-collar jobs rather than teleworkable, white-collar jobs for both native and immigrant workers. As a result, not accounting for the systematic selection of migrants into the regions where there is less supply of teleworkable jobs would lead the OLS result to underestimate the true effect of immigration on natives’ likelihood to be employed in such occupations. The difference between the IV and the OLS estimates increases significantly in our preferred specifications in columns 3 and 6 in Table 7 which include year and NUTS2-region fixed effects. Potential explanations might be the weaker power of the instrument when more granular geographical effects are accounted for, as well as the endogenous sorting of immigrant workers into the EU regions. Tables 8 and 9 examine the potentially heterogeneous impact of immigration on the share of natives in teleworkable jobs across gender and age groups. In both cases, we observe a positive impact of migrants’ presence across different subgroups of natives. In terms of age groups (Table 8), the effect is slightly larger for younger workers (aged 25-34) when compared to prime-age workers (aged 35-54). By gender, the coefficient is somewhat larger for female native-born 26 workers than for male native-born workers, although the differences are not statistically significant (Table 9). Heterogeneous effects are more prevalent across different educational groups within the labor force (Table 10). We split the first stage of the IV specification between the share of migrants with high levels of education in the total population of highly educated workers, and the share of migrants with lower levels of education in the total population with lower levels of education. The instruments are still strong with country fixed effects, although somewhat weaker when NUTS2- region fixed effects are included. The F-statistic for the excluded instrument falls below 10 (between 7 and 8.5), possibly signaling a weak instrument problem. The second-stage results show that the coefficient estimates for γ are large in the case of the impact on highly educated native workers from the presence of immigrants with both high and low levels of education. The magnitude of the coefficient indicates that a 1 percentage point increase in the share of immigrants with tertiary education increases the share of tertiary-educated natives in teleworkable jobs by around 1.5 percentage points. An increase in the number of immigrants with lower levels of education causes highly educated migrants to move to teleworkable occupations at an even higher rate – around 2 percentage points. In contrast, immigration has no significant impact on the native workers with lower levels of education, and we do not see any change in their share of telework jobs. These heterogeneous effects by education level point towards a higher degree of complementarity between immigrant workers and highly educated natives – possibly due to a greater ability of well-educated native workers to switch occupations, and to stronger complementarities between natives and immigrants at higher skill levels. The possibility of working from home shields workers from the negative economic outcomes arising from mobility restrictions and lockdowns during the COVID-19 pandemic. As discussed in Section 3, certain occupations are deemed essential by governments. These occupations do not face the same mobility restrictions because they are important for the functioning of the economy and the delivery of critical services such as health care, transportation, and safety. A job that has at least one these features – either essential or teleworkable – protects the job holder from the supply-side threats, making such occupation “income-safe.” Table 11 reports the results of a similar IV estimation for the different categories of jobs, where the dependent variable is the share of natives employed in income-safe jobs which is our third measure of risk exposure. We include 27 the OLS and the IV results for teleworkable jobs (originally reported in column 1 in Table 6, and columns 3 and 6 in Table 7) in columns 1 to 3 for ease of comparison. Interestingly, while there is a strong positive impact of migration on the expansion of telework jobs, the effect is not significant for essential occupations. This is likely to be related to the different nature of these jobs. While telework jobs tend to be concentrated among the highly educated and highly paid information technology and service jobs, a large share of essential jobs is found in low-paid occupations in agriculture, transport, and domestic work. While there are high-skilled occupations such as medical professionals among the essential jobs, they do not necessarily dominate the list. Thus, we find there is a positive but mostly insignificant impact of immigration on income-safe jobs. The impact of immigration on the share of native workers in essential jobs does not show much variation across education categories. As seen in Table 12, the impact of immigrants of high or low education levels on natives with the same education levels is not significant in terms of essential jobs. However, we see that highly educated natives tend to move out of essential jobs, in the presence of more immigrants with low levels of education (columns 3 and 4, Table 12). In terms of combined effects, as it was the case with the previous table, there is no significant impact on movement toward or away from income-safe jobs for natives in either education category. Another source of exposure to COVID-19 risk that has become salient with the onset of the pandemic is the potential exposure to health risks. These tend to arise from the necessity to have face-to-face interactions with other people as workers perform the duties required by their jobs. As previously mentioned, we define jobs that place one’s health at risk as those that are essential (continuing even during the peak of the pandemic) but cannot be performed from home; these jobs instead require frequent in-person contact with customers, suppliers or other colleagues. Around 17% of the jobs in our sample fall into the high health-risk category. As Table 13 shows, the impact of migration on the share of health-safe jobs among natives is positive and highly significant. In other words, migrants encourage and enable natives to move into safer jobs both in terms of income and health exposure. A 1 percentage point increase in the share of immigrants in the labor force leads to an almost 0.5 percentage point increase in the share of health-safe jobs among natives. By education group, this positive impact is again concentrated among tertiary-educated natives in locations where the impact of immigration of workers with low levels of education is especially 28 strong. In contrast, natives with low levels of education experience no significant change in the share of health-safe jobs due to migration. 5.2 Wage and Employment Effects of Immigration Results reported in the previous section show strong compositional effects where native-born workers move to certain occupations in response to labor supply shocks derived from the unequal composition of migration waves across Western European regions. In this section, we examine whether immigration impacted natives’ overall employment levels. Our goal is to shed light on whether the increase in the share of native workers employed in jobs protected from COVID-19 was purely driven by a reallocation from less safe to safer jobs, or also by changes in the overall employment levels of native workers. Table 14 reports the OLS and IV estimates for the employment effects of immigration on natives, which also control for year and NUTS2-region fixed effects. The OLS estimates (column 1, Table 14) show a slightly positive impact of migration on employment rates of natives aged 25-64. However, once the endogeneity issues are addressed, the coefficients for the IV regressions become smaller, and the overall effect becomes statistically insignificant. We extend the analysis to different age groups as shown in Table 15. Similarly, we find again a “close to zero” impact on employment for the prime working age population groups between 25- 34 and 35-54 years old. Interestingly, there is a significant reduction in employment rates for the youngest group (between the ages of 15 to 24). This result is not associated with a significant drop in labor force participation, but rather with higher unemployment rates. In contrast, older working- age natives (ages 55-64) see an increase in their employment rate due to higher immigration rates, indicating that older workers are able to stay in the labor force longer. The overall null effect of migration on employment however masks important gender differences (Table 16). While there is a non-significant reduction in employment rates of male natives (columns 1 and 2), there is an increase in employment rates of females due to the arrival of migrants (Columns 3 and 4 of Table 16). These findings are in line with earlier studies in the US (Cortes and Tessada, 2011), and suggest that migrants – especially those working in personal care, domestic care or household 29 services – can alleviate the pressure on women due to activities (such as housekeeping, and caring for children and the elderly) that traditionally keep them detached from the labor market. Finally, when we explore the impact across different education groups, we again fail to find any statistically significant effect. Even though the employment effect is positive for the highly educated and negative for those with lower levels of education, the effects are not precisely estimated. Finally, we examine the impact of immigration on natives’ wages, the focus of many papers in the literature (see World Bank 2018 for a review). Results in Table 17 show that we again fail to find any significant impact of migration on natives’ wages. As a caveat to these findings, we underscore that the EU Labour Force Survey only provides information on the wage decile of workers, not their actual wages. In order to circumvent this problem, we construct a variable that measures the gap between the actual and expected probability of being in the top three income deciles given individual-level characteristics (age, gender, education levels), and that varies by country and year. The impact of migrants’ presence on the wage-gap variable, however, remains not statistically different from zero, a result that holds when looking at the share of top three wage deciles of teleworkable, essential, income-safe and health-safe jobs. These results are robust to alternative versions of the dependent variable, such as the share of employees in the bottom three deciles, or the income decile estimated using ordered probit regressions in the initial stage. Therefore, the adjustment in the labor market after the arrival of migrants seems to happen through the composition of employment but does not affect the relative or absolute wages of the different types of occupations in line with a more competitive and mobile labor market. 6. Conclusion This paper provides evidence that exposure to negative economic shocks and health risks associated with the COVID-19 pandemic is greater among immigrant workers than among native- born people in high-income Western European countries. Using data from the European Union Labour Force Survey for 2011-2018, we show that this greater risk exposure stems from systematic differences in occupations of immigrant and native workers. Immigrant workers have significantly lower presence in occupations that can be performed from home (telework), and higher presence 30 in occupations listed as “essential.” Moreover, immigrants are more likely to be employed in jobs that combine both features – taking on work that is deemed to be essential and requires face-to- face contact. Thus, immigrant workers are more likely to bear a greater brunt of the economic and health risks of the pandemic than native-born workers. The arrival of immigrants to a regional labor market in the EU, in itself, had a causal impact on the occupational shift of natives, contributing to lower levels of exposure to the pandemic related risks. Our results show that native-born workers – especially those with tertiary education living in regions with a higher share of migrant workers – are more often employed in occupations that can be carried out from home, which we refer to as the ability to telework. These effects are large in magnitude and hold when Bartik-type instruments are used to account for the endogeneity of immigrants’ locations and occupation choices. However, the positive impact of immigration – of migrants with both high and low levels of education – is concentrated among native workers with high levels of education. Those native workers with less schooling, in contrast, do not seem to be affected much from immigration. In addition to the ability to telework, we explore several other occupational characteristics in relation to COVID-19 related risks. The first one is whether an occupation is considered as essential. Based on governmental regulations, essential jobs cover a wide skill range from low skilled ones (in certain retail sectors) to high-skilled ones in medical professions. Next, we classify certain jobs as income-safe if they are either essential or could be performed via telework. Our results show that the presence of immigrants in a certain subregion does not have a significant effect on the share of native workers employed in essential or income safe jobs. This is different from the case of teleworkable jobs and it is likely to be due to the high share of low-skilled jobs in essential occupations (such as retail and transportation) that attract large numbers of migrants. However, these low-skilled essential jobs tend to entail larger health risks as they often require physical interactions with other people. Based on our last measure of job exposure to health-risks we also find that the presence of immigrants causes native workers to move towards occupations that are less exposed to health risks. These effects are again stronger for the highly educated natives. 31 We finally analyze the potential impact of immigration on changes in the overall employment levels and wages of native workers as it impacts occupational composition. In contrast with the large occupational shifts, we find no absolute employment or wage effects in the EU. These results are consistent across education levels. In many countries, many migrants perform the dirty, dangerous and demeaning jobs – sometimes referred to as “3D jobs” – that natives are unwilling to do. The evidence from the COVID-19 pandemic is pointing in the same direction. As more evidence emerges, especially on labor market outcomes during the pandemic, we expect the immigration policies of the post-COVID era will be shaped accordingly. 32 References Adams, Richard H., and John Page. 2005. “Do International Migration and Remittances Reduce Poverty in Developing Countries?” World Development, 33 (10): 1645–69. Adams-Prassl, Abi, Teodora Boneva, Marta Golin and Christopher Rauh. 2020. “Inequality in the impact of the coronavirus shock: Evidence from real time surveys,” Journal of Public Economics, Volume 189, September 2020, 104245 Alacevich, Caterina, and Catia Nicodemo. 2019. “Immigration and work-related injuries: Evidence from Italian administrative data.” IZA Discussion Paper No. 12510. Bonn: Germany: Institute of Labor Economics. Avdiu, Besart, and Gaurav Nayyar. 2020. “When face-to-face interactions become an occupational hazard: Jobs in the time of COVID-19,” World Bank Policy Research Working Paper, (9240). Basso, Gaetano and Giovanni Peri. 2015 “The association between immigration and labor market outcomes in the United States” IZA Discussion Paper No. 9436. Bonn: Germany: Institute of Labor Economics. Basso, Gaetano, Tito Boeri, Alessandro Caiumi and Marco Paccagnella, 2020. “The new hazardous jobs and worker reallocation,” OECD Social, Employment and Migration Working Papers 247, OECD Publishing. Beaman, Lori A. 2012. “Social Networks and the Dynamics of Labor Market Outcomes: Evidence from Refugees Resettled in the U.S.” Review of Economic Studies, 79 (1): 128-161. Beine, Michel, Frederic Docquier and Caglar Ozden. 2011. “Diasporas,” Journal of Development Economics 95 (1): 30-41 Béland, Louis-Philippe, Abel Brodeur and Taylor Wright. 2020. “The Short-Term Economic Consequences of COVID-19: Exposure to Disease, Remote Work and Government Response.” IZA Discussion Paper No. 13159. Bellés-Obrero, Cristina, Nicolau M. Bassols, and Judit V. Castello. [forthcoming]. “Safety at work and immigration.” Journal of Population Economics. Bond, Timothy, Osea Giuntella and Jakub Lonsky. 2020. “Immigration and Work Schedules: Theory and Evidence,” IZA Discussion Paper No. 13236. Bonn: Germany: Institute of Labor Economics. Borjas, George J. 2015. “The Slowdown in the Economic Assimilation of Immigrants: Aging and Cohort Effects Revisited Again.” Journal of Human Capital 9 (4): 483–517. Borjas, George J. and Hugh Cassidy. 2020. “The Adverse Effect of the COVID-19 Labor Market Shock on Immigrant Employment,” NBER Working Paper 27243 (June), National Bureau of Economic Research, Cambridge, MA. Brücker, Herbert, Albrecht Glitz, Adrian Lerche and Agnese Romiti. [forthcoming]. “Occupational Recognition and Immigrant Labor Market Outcomes,” Journal of Labor Economics. Cajner, Tomaz, Andrew Figura, Brendan M. Price, David Ratner, and Alison Weingarden. 2020. “Reconciling Unemployment Claims with Job Losses in the First Months of the COVID-19 Crisis.” FEDS Working Paper 2020-055. 33 Casarico, Alessandra, and Salvatore Lattanzio. 2020. “The Heterogenous Effects of COVID-19 on Labor Market Flows: Evidence from Administrative Data.” Covid Economic, 52. Chiswick, Barry R., Yew Liang Lee, and Paul W. Miller. 2005. “A Longitudinal Analysis of Immigrant Occupational Mobility: A Test of the Immigrant Assimilation Hypothesis.” International Migration Review 39 (2): 332–53. Coibion, Olivier, Yuriy Gorodnichenko and Michael Weber. 2020. “Labor Markets During the COVID-19 Crisis: A Preliminary View,” NBER Working Paper 27017. Congressional Budget Office. 2006. “Potential influenza pandemic: Possible macroeconomic effects and policy issues”. https://www.cbo.gov/sites/default/files/109th-congress-2005-2006/reports/12-08- birdflu.pdf Cortes, Patricia and Jose Tessada. 2011. “Low-Skilled Immigration and the Labor Supply of Highly Skilled Women,” American Economic Journal: Applied Economics 3 (3): 88-123. del Rio-Chanona, R. Maria, Penny Mealy, Anton Pichler, Francois Lafond and Doyne Farmer. 2020. “Supply and demand shocks in the COVID-19 pandemic: An industry and occupation perspective,” arXiv preprint arXiv:2004.06759. D’Amuri, Francesco, and Giovanni Peri. 2014. “Immigration, jobs, and employment protection: evidence from Europe before and during the great recession.” Journal of the European Economic Association, 12(2), pp.432-464. Dillender, Marcus, and Melissa McInerney. 2020. “The role of Mexican immigration to the United States in improved workplace safety for natives from 1980 to 2015,” Journal of health economics, 70, p.102280. Dingel, Jonathan I, and Brent Neiman. 2020. “How many jobs can be done at home?” Journal of Public Economics, Volume 189, 104235. Dustmann, Christian, Francesca Fabbri and Ian Preston. 2005 “The impact of immigration on the British labour market” The Economic Journal, 115 (507): F324-F341. Dustmann, Christian, Tommaso Frattini and Ian P. Preston, 2013. “The Effect of Immigration along the Distribution of Wages,” Review of Economic Studies, 80 (1): 145-173. Edo, Anthony, Yvonne Giesing, Jonhathan Öztunc, J. and Panu Poutvaara. 2019. “Immigration and electoral support for the far-left and the far-right.” European Economic Review, 115: 99-143. Elsner, Benjamin, 2013. “Emigration and wages: The EU enlargement experiment.” Journal of International Economics, 91(1): 154-163. European Commission. 2020. “Communication from the Commission Guidelines concerning the exercise of the free movement of workers during COVID-19 outbreak,” 2020/C 102 I/03. https://eur- lex.europa.eu/legal-content/EN/TXT/PDF/?uri=CELEX:52020XC0330(03)&from=EN Fasani, Francesco and Jacopo Mazza. 2020a. “Immigrant Key Workers: Their Contribution to Europe's COVID-19 Response,” IZA Policy Papers No. 155 (April), Institute of Labor Economics, Bonn, Germany. Fasani, Francesco, and Jacopo Mazza. 2020b. “Being on the Frontline? Immigrant Workers in Europe and the COVID-19 Pandemic.” IZA Discussion Paper 13963 (December), Institute of Labor Economics, Bonn, Germany. 34 Friedberg, Rachel M. 2000. “You Can’t Take It with You? Immigrant Assimilation and the Portability of Human Capital.” Journal of Labor Economics 18 (2): 221–51. Foged, Mette and Giovanni Peri. 2016. “Immigrants' Effect on Native Workers: New Analysis on Longitudinal Data,” American Economic Journal: Applied Economics, 8 (2): 1-34. Garrote Sanchez, Daniel, Nicolas Gomez Parra, Caglar Ozden, and Bob Rijkers. 2020. “Which jobs are most vulnerable to COVID-19? Analysis of the European Union,” World Bank Research & Policy Briefs No. 34. Washington, D.C.: World Bank Group. Garrote Sanchez, Daniel, Nicolas Gomez Parra, Caglar Ozden, Bob Rijkers, Mariana Viollaz and Hernan Winkler. 2021. “Who on Earth Can Work from Home?, World Bank Research Observer. Giuntella, Osea, Fabrizio Mazzonna, Catia Nicodemo, and Carlos Vargas-Silva. 2019. “Immigration and the reallocation of work health risks.” Journal of Population Economics, 32: 1009–42. Goldsmith-Pinkham, Paul, Isaac Sorkin, and Henry Swift. 2020. “Bartik instruments: What, when, why, and how.” American Economic Review 110, 2586–2624. Gross, Dominique M. and Nicholas Schmitt. 2003. “The role of cultural clustering in attracting new immigrants.” Journal of Regional Science, 43 (2): 295-318. Jaeger, David A. 2006. “Green Cards and the Location Choices of Immigrants in the United States.” IZA Discussion Paper No 2145. Jaeger, David A., Joakim Ruist, and Jan Stuhler. 2018. “Shift-share instruments and the impact of immigration.” NBER Working Papers 24285, National Bureau of Economic Research, Inc. Kramarz, Francis and Oskar Nordström Skans. 2014. “When Strong Ties are Strong: Networks and Youth Labour Market Entry,” Review of Economic Studies, 81 (3): 1164–1200. Liu, Ou, and Tam Mai. 2020. “Employment during the COVID-19 Pandemic: Collapse and Early Recovery.” Mimeo. Mitaritonna, Cristina, Gianluca Orefice, and Giovanni Peri. 2017. “Immigrants and firms’ outcomes: Evidence from France,” European Economic Review, 96, issue C, 62-82. Mongey, Simon, Laura Pilossoph and Alex Weinberg. 2020. “Which Workers Bear the Burden of Social Distancing Policies?”, University of Chicago, Becker Friedman Institute for Economics Working Paper No. 2020-51. Montenovo, Laura, Xuan Jiang, Felipe Lozano Rojas, Ian M. Schmutte, Kosali I. Simon, Bruce A. Weinberg and Coady Wing. 2020. “Determinants of disparities in covid-19 job losses.” NBER Working Paper 27132. Munshi, Kaivan. 2003. “Networks in the Modern Economy: Mexican Migrants in the U. S. Labor Market,” The Quarterly Journal of Economics, 118 (2): 549-599. Munshi, Kaivan and Mark Rosenzweig. 2006. “Traditional Institutions Meet the Modern World: Caste, Gender, and Schooling Choice in a Globalizing Economy,” American Economic Review, 96(4): 1225-1252. Orrenius Pia M., and Madeline Zavodny. 2013. "Immigrants in risky occupations," in: Amelie F. Constant and Klaus F. Zimmermann (eds.), International Handbook on the Economics of Migration, chapter 11, 214- 226, Edward Elgar Publishing. 35 Orrenius, Pia M, and Madeline Zavodny. 2009. “Do immigrants work in riskier jobs?” Demography, Volume 46-Number 3, August 2009: 535–551. Papanikolaou, Dimitris and Lawrence D.W. Schmidt. 2020. “Working Remotely and the Supply-side Impact of Covid-19.” NBER Working Paper 27330. Patel, Krishna, and Francis Vella. 2013. “Immigrant Networks and Their Implications for Occupational Choice and Wages”, Review of Economics and Statistics, 95 (4): 1249-1277. Peri, Giovanni, and Chad Sparber. 2009. Task specialization, immigration, and wages. American Economic Journal: Applied Economics, 1(3), pp.135-69. Peri, Giovanni, and Chad Sparber. 2011. Highly educated immigrants and native occupational choice. Industrial Relations: a journal of economy and society, 50(3): 385-411. Sparber, Chad and Madeline Zavodny. 2020. "Immigration, Working Conditions, and Compensating Differentials," IZA Discussion Papers 13663, Institute of Labor Economics (IZA). Stock, James H. and Motohiro Yogo. 2005. “Testing for Weak Instruments in Linear IV Regression,” Ch. 5 in J.H. Stock and D.W.K. Andrews (eds), Identification and Inference for Econometric Models: Essays in Honor of Thomas J. Rothenberg, Cambridge University Press. World Bank. 2018. Moving for Prosperity: Global Migration and Labor Markets. World Bank, Washington DC. Yasenov, Vasil. 2020. “Who Can Work from Home?” IZA Discussion Paper No. 13197. 36 Table 1: Descriptive statistics of the EU15+EFTA population in 2018 All population 25-64 Males 25-64 Females 25-64 Size % % (% Mean % tertiary Mean % tertiary Mean tertiary males total) age education age education age education Native-born 81.3 50.2 45.4 33.5 45.2 32.8 45.6 34.3 Migrants 18.7 48.4 43.0 31.2 43.0 30.0 43.0 32.3 Migrants from EU15+EFTA 3.1 51.8 45.3 43.0 45.1 41.9 45.4 44.2 Migrants from NMS13 3.6 46.1 41.9 27.3 41.7 23.7 42.0 30.3 Migrants from non-EU Europe 3.5 47.6 43.5 21.2 43.6 18.0 43.5 24.0 Migrants from Rest of World 8.1 48.5 42.5 31.4 42.4 31.4 42.5 31.3 Source: Own elaboration based on EU-LFS. Table 2: Labor market statistics of the EU15 population (25-64) in 2018 Total High education Low education Total Natives Migrants Natives Migrants Natives Migrants Population (1,000) 270,026 219,628 50,398 73,571 15,602 145,775 34,440 Share 100.0 81.3 18.7 82.5 17.5 80.9 19.1 Unemployment rate 6.1 5.4 9.0 3.4 6.6 6.6 10.2 Employment rate 76.3 77.6 70.7 87.5 79.2 72.6 67.0 % telework jobs 38.4 40.9 27.3 62.1 52.1 28.0 14.1 % essential jobs 38.4 37.6 41.8 43.4 44.5 34.1 40.4 % income-safe jobs 64.2 65.3 59.3 81.4 75.6 55.6 50.6 % health-safe jobs 83.1 84.4 77.1 88.2 83.1 82.1 73.9 Source: Own elaboration based on EU-LFS. 37 Table 3: Labor market statistics of the EU15 population (25-64) in 2018 by gender Total High education Low education Natives Migrants Natives Migrants Natives Migrants Males Females Males Females Males Females Males Females Males Females Males Females Population (1,000) 110,159 109,469 24,418 25,980 36,083 37,487 7,267 8,335 73,928 71,848 16,956 17,484 Share 81.9 80.8 18.1 19.2 83.2 81.8 16.8 18.2 81.3 80.4 18.7 19.6 Unemployment rate 5.2 5.5 8.3 9.8 3.0 3.7 5.7 7.6 6.4 6.7 9.6 11.1 Employment rate 82.3 72.8 79.9 62.1 90.2 84.8 86.5 72.8 78.5 66.6 77.1 57.2 % telework jobs 36.4 45.9 24.9 30.2 61.2 63.1 52.7 51.5 22.6 34.5 11.7 17.3 % essential jobs 37.7 37.6 36.9 47.8 44.4 42.4 45.4 43.7 33.9 34.4 32.8 50.4 % income-safe jobs 61.0 70.2 51.6 68.6 81.0 81.7 74.1 77.1 49.8 62.6 40.7 63.4 % health-safe jobs 88.6 79.7 85.7 66.7 92.2 84.1 88.0 78.0 86.6 76.8 84.7 59.8 Source: Own elaboration based on EU-LFS. Table 4: Share of workers in the top 3 deciles of wages in 2018 Total High education Low education Total Natives Migrants Natives Migrants Natives Migrants Total employees 33.0 35.0 24.4 57.4 47.1 23.1 15.0 Telework employees 47.7 47.9 46.8 62.0 59.0 30.9 25.5 Essential employees 34.6 37.2 23.2 59.5 50.7 22.1 10.9 Income-safe employees 38.2 40.3 28.7 60.2 53.0 24.9 13.5 Health-safe employees 35.6 37.3 28.1 58.5 49.5 24.8 17.7 Source: Own elaboration based on EU-LFS. 38 Table 5: Share of workers in the top 3 deciles of wages in 2018 by gender Total High education Low education Natives Migrants Natives Migrants Natives Migrants Males Females Males Females Males Females Males Females Males Females Males Females Total employees 47.0 22.4 32.6 14.8 34.3 42.9 22.5 34.6 34.3 11.2 22.5 5.6 Telework employees 66.6 33.0 61.9 32.6 77.0 47.8 73.1 44.3 51.0 17.8 39.1 14.9 Essential employees 50.1 23.3 35.2 12.3 72.1 45.6 63.4 35.6 34.9 8.4 19.9 3.6 Income-safe employees 55.8 26.6 41.1 18.0 74.6 46.0 66.5 39.3 39.7 12.9 22.9 6.2 Health-safe employees 48.4 24.6 34.1 18.9 72.4 43.7 61.0 36.4 35.0 12.6 23.3 7.5 Source: Own elaboration based on EU-LFS. 39 Table 6: Impact of migration on the share of natives 25-64 in telework jobs (OLS estimates) Dependent variable: share of natives in telework jobs (% of the total native population 25-64) Explanatory variable (fborn): Share of migrants 25-64 (% of the total population 25-64) (1) (2) (3) VARIABLES OLS OLS OLS Fborn 0.276*** 0.243*** 0.126** (0.056) (0.052) (0.055) Constant -0.054*** -0.038*** -0.000 (0.009) (0.012) (0.016) Observations 982 982 982 R-squared 0.236 0.795 0.959 Year FE NO YES YES Country FE NO YES NO Nuts2 FE NO NO YES Clustered standard errors at the nuts2 level in parentheses; *** p<0.01, ** p<0.05, * p<0.1 40 Table 7: Impact of migration on the share of natives 25-64 in telework jobs (IV estimates) Dependent variable second stage: share of natives in telework jobs (% of the total native population 25- 64) Explanatory variable second stage (fborn): Share of migrants 25-64 (% of the total population 25-64) (1) (2) (3) (4) (5) (6) VARIABLES IV-PS IV-PS IV-PS IV-MOP IV-MOP IV-MOP First stage (dependent variable: fborn) Shiftshare (PS) 0.932*** 1.005*** 0.462*** (0.082) (0.0671) (0.124) Shiftshare (MOP) 0.879*** 0.952*** 0.701*** (0.078) (0.065) (0.181) Second stage (dependent variable: share of natives in telework jobs) fborn 0.286*** 0.247*** 0.686*** 0.290*** 0.250*** 0.661*** (0.056) (0.051) (0.225) (0.056) (0.050) (0.220) Constant -0.056*** -0.039*** -0.160** -0.056*** -0.040*** -0.153** (0.009) (0.011) (0.064) (0.009) (0.011) (0.063) Observations 982 982 982 982 982 982 F-stat first stage 130.4 224.6 13.9 125.5 212.2 15.0 R-squared (second 0.236 0.795 0.934 0.235 0.795 0.937 stage) Year FE NO YES YES NO YES YES Country FE NO YES NO NO YES NO Nuts2 FE NO NO YES NO NO YES Clustered standard errors at the nuts2 level in parentheses; *** p<0.01, ** p<0.05, * p<0.1 41 Table 8: Impact of migration on the share of natives 25-64 in telework jobs by age group Dependent variable: share of natives (either 25-34 or 35-54) in telework jobs Explanatory variable (fborn): Share of migrants 25-64 (% of the total population 25-64) (1) (2) (3) (4) (5) (6) (7) (8) Age 25-34 Age 35-54 Variable IV-PS IV-PS IV-MOP IV-MOP IV-PS IV-PS IV-MOP IV- MOP Fborn 0.242*** 0.789** 0.247*** 0.783** 0.266*** 0.719*** 0.270*** 0.628** (0.065) (0.313) (0.064) (0.311) (0.059) (0.279) (0.058) (0.290) Constant 0.010 -0.145 0.009 -0.143 -0.036*** -0.152* -0.037*** -0.126 (0.014) (0.089) (0.014) (0.089) (0.014) (0.080) (0.014) (0.083) Observations 982 982 982 982 982 982 982 982 R-squared 0.764 0.912 0.764 0.912 0.768 0.916 0.768 0.921 Year FE YES YES YES YES YES YES YES YES Country FE YES NO YES NO YES NO YES NO Nuts2 FE NO YES NO YES NO YES NO YES Clustered standard errors at the nuts2 level in parentheses; *** p<0.01, ** p<0.05, * p<0.1 42 Table 9: Impact of migration on the share of natives 25-64 in telework jobs by gender Dependent variable: share of natives 25-64 (either males or females) in telework jobs Explanatory variable (fborn): Share of migrants over the total working-age population (1) (2) (3) (4) (5) (6) (7) (8) Males Females VARIABLES IV-PS IV-PS IV-MOP IV-MOP IV-PS IV-PS IV-MOP IV-MOP Fborn 0.290*** 0.565* 0.293*** 0.526** 0.193*** 0.814*** 0.197*** 0.799*** (0.058) (0.290) (0.056) (0.261) (0.045) (0.214) (0.044) (0.236) Constant -0.068*** -0.150* -0.068*** -0.139* -0.008 -0.170*** -0.009 -0.166** (0.012) (0.083) (0.012) (0.075) (0.012) (0.061) (0.012) (0.067) Observations 982 982 982 982 982 982 982 982 R-squared 0.743 0.924 0.743 0.926 0.804 0.910 0.804 0.912 Year FE YES YES YES YES YES YES YES YES Country FE YES NO YES NO YES NO YES NO Nuts2 FE NO YES NO YES NO YES NO YES Clustered standard errors at the nuts2 level in parentheses; *** p<0.01, ** p<0.05, * p<0.1 43 Table 10: Impact of migration on the share of natives 25-64 in telework jobs by education levels Dependent variable second stage: share of high- educated or low- educated natives 25-64 in telework jobs Explanatory variable second stage (fborn): Share of high educated (hs) or low educated (ls) migrants (% of the high- educated or low educated population 25-64) (1) (2) (3) (4) (5) (6) (7) (8) (9) (10) (11) (12) High-educ migrants on high- educ natives Low- educ migrants on high- educ natives Low- educ migrants on low- educ natives VARIABLES IV-PS IV -PS IV -MOP IV-MOP IV-PS IV -PS IV -MOP IV-MOP IV-PS IV -PS IV -MOP IV-MOP First stage (dependent variable: fborn) shiftshare_ps2002 0.856*** 0.690*** 0.955*** 0.203 0.955*** 0.203 (0.099) (0.258) (0.073) (0.165) (0.073) (0.165) shiftshare_mop2002 0.681*** 0.576*** 0.941*** 0.516* 0.941*** 0.516* (0.080) (0.205) (0.072) (0.301) (0.072) (0.301) Second stage Fborn 0.289*** 1.316** 0.303*** 1.531** 0.214*** 2.027** 0.218*** 2.055*** 0.195*** 0.122 0.198*** -0.023 (0.084) (0.553) (0.083) (0.632) (0.058) (0.792) (0.056) (0.747) (0.042) (0.266) (0.041) (0.260) Constant -0.043** -0.330** -0.045*** -0.391** -0.029** -0.535** -0.030** -0.543** -0.034*** -0.006 -0.035*** 0.035 (0.017) (0.159) (0.016) (0.181) (0.013) (0.226) (0.012) (0.213) (0.011) (0.076) (0.011) (0.074) Observations 980 980 980 980 982 982 982 982 982 982 982 982 R-squared second 0.803 0.851 0.802 0.824 0.813 0.776 0.813 0.771 0.740 0.954 0.740 0.954 stage F-stat first stage 78.6 7.0 73.4 7.0 258.2 7.3 241.2 8.5 258.2 7.3 241.2 8.5 Year FE YES YES YES YES YES YES YES YES YES YES YES YES Country FE YES NO YES NO YES NO YES NO YES NO YES NO Nuts2 FE NO YES NO YES NO YES NO YES NO YES NO YES Clustered standard errors at the nuts2 level in parentheses; *** p<0.01, ** p<0.05, * p<0.1 44 Table 11: Impact of migration on the share of natives 25-64 in telework, essential and income-safe jobs Dependent variable: share of natives in each type of occupation (% of the total native population 25-64) Explanatory variable (fborn): Share of migrants 25-64 (% of the total population 25-64) (1) (2) (3) (4) (5) (6) (7) (8) (9) Telework jobs Essential jobs Income-safe jobs VARIABLES OLS IV-PS IV-MOP OLS IV-PS IV-MOP OLS IV-PS IV-MOP fborn 0.126** 0.686*** 0.661*** 0.080 -0.272 -0.248 0.178* 0.364 0.362 (0.055) (0.225) (0.220) (0.063) (0.194) (0.216) (0.091) (0.256) (0.275) Constant -0.000 -0.160** -0.153** -0.015 0.086 0.079 -0.018 -0.071 -0.070 (0.016) (0.064) (0.063) (0.018) (0.055) (0.062) (0.026) (0.073) (0.078) Observations 982 982 982 982 982 982 982 982 982 R-squared 0.959 0.934 0.937 0.909 0.892 0.895 0.941 0.939 0.939 Year FE YES YES YES YES YES YES YES YES YES Country FE NO NO NO NO NO NO NO NO NO Nuts2 FE YES YES YES YES YES YES YES YES YES Clustered standard errors at the nuts2 level in parentheses; *** p<0.01, ** p<0.05, * p<0.1 45 Table 12: Impact of migration on the share of natives 25-64 in essential and income-safe jobs by education levels Dependent variable: share of high- educated or low- educated natives 25-64 in essential or income-safe jobs Explanatory variable (fborn): Share of high educated (hs) or low educated (ls) migrants (% of the high- educated or low educated population 25- 64) (1) (2) (3) (4) (5) (6) (7) (8) (9) (10) (11) (12) Essential jobs Income-safe jobs High-educ mig on Low-educ mig on Low-educ mig on High-educ mig on Low-educ mig on Low-educ mig on high-educ natives high-educ natives low-educ natives high-educ natives high-educ natives low-educ natives VARIABLES IV -PS IV-MOP IV -PS IV -MOP IV -PS IV-MOP IV -PS IV-MOP IV -PS IV-MOP IV -PS IV-MOP Fborn -0.409 -0.337 -1.144** -1.069* -0.246 -0.177 0.745 0.922* 0.278 0.338 -0.043 -0.103 (0.319) (0.363) (0.518) (0.549) (0.328) (0.298) (0.475) (0.539) (0.492) (0.466) (0.436) (0.424) Constant 0.115 0.094 0.325** 0.304* 0.081 0.061 -0.189 -0.240 -0.058 -0.075 0.045 0.061 (0.090) (0.103) (0.148) (0.157) (0.094) (0.085) (0.136) (0.155) (0.141) (0.133) (0.125) (0.121) Observations 980 980 982 982 982 982 980 980 982 982 982 982 R-squared 0.860 0.864 0.775 0.786 0.852 0.860 0.909 0.899 0.927 0.926 0.910 0.908 Year FE YES YES YES YES YES YES YES YES YES YES YES YES Country FE NO NO NO NO NO NO NO NO NO NO NO NO Nuts2 FE YES YES YES YES YES YES YES YES YES YES YES YES Clustered standard errors at the nuts2 level in parentheses; *** p<0.01, ** p<0.05, * p<0.1 46 Table 13: Impact of migration on the share of natives 25-64 in health-safe jobs Dependent variable: share of natives in health-safe jobs (not essential | not F2F | teleworkable) (% of the total native population 25-64) Explanatory variable (fborn): (1) (2) (3) (4) (5) (6) (7) (8) Total High-educ migrants Low-educ migrants Low-educ migrants on high-educ natives on high-educ natives on low-educ natives VARIABLES IV-PS IV-MOP IV-PS IV-MOP IV-PS IV-MOP IV-PS IV-MOP Fborn 0.451*** 0.485*** 0.624** 0.748** 1.686** 1.714** 0.030 -0.019 (0.148) (0.175) (0.314) (0.361) (0.729) (0.728) (0.292) (0.276) Constant -0.111*** -0.121** -0.140 -0.175* -0.444** -0.452** -0.004 0.011 (0.042) (0.050) (0.089) (0.102) (0.207) (0.207) (0.084) (0.079) Observations 982 982 980 980 982 982 982 982 R-squared 0.786 0.778 0.530 0.487 0.598 0.599 0.835 0.835 Year FE YES YES YES YES YES YES YES YES Country FE NO NO NO NO NO NO NO NO Nuts2 FE YES YES YES YES YES YES YES YES Clustered standard errors at the nuts2 level in parentheses; *** p<0.01, ** p<0.05, * p<0.1 47 Table 14: Impact of migration on the share of natives 25-64 in employment Dependent variable: share of natives in employment (% of the total native population 25-64) Explanatory variable (fborn): (1) (2) (3) VARIABLES OLS IV -PS IV -MOP Fborn 0.149** 0.087 0.121 (0.068) (0.330) (0.357) Constant -0.028 -0.011 -0.020 (0.019) (0.094) (0.102) Observations 982 982 982 R-squared 0.939 0.939 0.939 Year FE YES YES YES Country FE NO NO NO Nuts2 FE YES YES YES Clustered standard errors at the nuts2 level in parentheses; *** p<0.01, ** p<0.05, * p<0.1 48 Table 15: Impact of migration on the share of natives in employment by age group Dependent variable: share of natives (either 25-34 or 35-54) in employment Explanatory variable: Share of migrants over the total working-age population (1) (2) (3) (4) (5) (6) (7) (8) Age 15-24 Age 25-34 Age 35-54 Age 55-64 VARIABLES IV-PS IV-MOP IV-PS IV-MOP IV-PS IV-MOP IV-PS IV-MOP Fborn -0.753** -0.802** 0.104 0.095 -0.279 -0.347 0.940* 1.074** (0.341) (0.337) (0.480) (0.497) (0.514) (0.544) (0.488) (0.454) Constant 0.276*** 0.290*** 0.221 0.193 0.128 0.147 -0.343** -0.382*** (0.097) (0.096) (0.141) (0.141) (0.147) (0.155) (0.139) (0.129) Observations 982 982 982 982 982 982 982 982 R-squared 0.930 0.927 0.835 0.838 0.895 0.893 0.898 0.891 Year FE YES YES YES YES YES YES YES YES Country FE NO NO NO NO NO NO NO NO Nuts2 FE YES YES YES YES YES YES YES YES Clustered standard errors at the nuts2 level in parentheses; *** p<0.01, ** p<0.05, * p<0.1 49 Table 16: Impact of migration on the share of natives in employment by gender and education Dependent variable: share of natives 25-64 in employment (males, females, high-educated and low- educated) Explanatory variable: Share of migrants over the total working-age population (1) (2) (3) (4) (5) (6) (7) (8) (9) (10) Males Females High-educ Low-educ Low-educ migrants on high- migrants on high- migrants on low- educ natives educ natives educ natives VARIABLES IV-PS IV-MOP IV-PS IV-MOP IV-PS IV-MOP IV-PS IV-MOP IV-PS IV-MOP Fborn -0.343 -0.296 0.568* 0.600* 0.579 0.738 -0.158 -0.009 -0.455 -0.448 (0.459) (0.480) (0.318) (0.339) (0.501) (0.541) (0.465) (0.468) (0.670) (0.613) Constant 0.073 0.060 -0.107 -0.116 -0.125 -0.170 0.084 0.041 0.132 0.130 (0.131) (0.137) (0.091) (0.097) (0.143) (0.155) (0.133) (0.134) (0.191) (0.175) Observations 982 982 982 982 980 980 982 982 982 982 R-squared 0.854 0.856 0.944 0.943 0.904 0.898 0.911 0.913 0.900 0.908 Year FE YES YES YES YES YES YES YES YES YES YES Country FE NO NO NO NO NO NO NO NO NO NO Nuts2 FE YES YES YES YES YES YES YES YES YES YES Clustered standard errors at the nuts2 level in parentheses; *** p<0.01, ** p<0.05, * p<0.1 50 Table 17: Impact of migration on the share of natives in top 3 deciles of income by type of occupation (1) (2) (3) (4) (5) (6) (7) (8) (9) (10) Total employees Income-safe Telework employees Essential employees Health-safe employees employees VARIABLES IV-PS IV-MOP IV-MOP IV-MOP IV-MOP IV-MOP IV-MOP IV-MOP IV-MOP IV-MOP Fborn -0.075 -0.085 -0.770 -0.880 -0.168 -0.322 -0.208 -0.314 -0.068 -0.025 (0.330) (0.378) (0.478) (0.552) (0.429) (0.487) (0.362) (0.413) (0.495) (0.570) Constant 0.036 0.039 0.308** 0.366** 0.071 0.114 0.117 0.146 0.039 0.028 (0.091) (0.105) (0.132) (0.153) (0.119) (0.135) (0.100) (0.114) (0.136) (0.157) Observations 781 781 781 781 781 781 781 781 781 781 R-squared 0.898 0.898 0.803 0.787 0.830 0.826 0.879 0.876 0.770 0.771 Year FE YES YES YES YES YES YES YES YES YES YES Country FE NO NO NO NO NO NO NO NO NO NO Nuts2 FE YES YES YES YES YES YES YES YES YES YES Clustered standard errors at the nuts2 level in parentheses; *** p<0.01, ** p<0.05, * p<0.1 51