
Figure 1
Foreign-born population in Colombia.
Sources: Colombian National Integrated Household Survey (GEIH) (2013–2019); Population Census (1993, 2005).

Figure 2
Venezuelan migrants in Colombia (2019).
Sources: Colombian National Integrated Household Survey (GEIH) (2019).
Table 1
Labor market effects of immigration
| (1) | (2) | (3) | |
|---|---|---|---|
| OLS | 2SLS | Test (1) = (2) (p-value) | |
| ln(hourly wage) | −0.73* (0.42) | −1.05*** (0.22) | 0.209 |
| ln(hours/week) | 0.08 (0.26) | 0.27 (0.27) | 0.399 |
| Unemployment | 0.01 (0.13) | −0.08 (0.07) | 0.255 |
| Labor force participation | −0.10 (0.12) | −0.21* (0.11) | 0.355 |
| K-P Wald stat | 23.35 | ||
| N | 474 | 474 | |
| Year FE, City FE | X | X |
Notes: Outcomes are residualized and multiplied by 100. Observations are weighted by city–year population. Column 3 presents a Hausman test for endogenous regressors with robust standard errors. Cluster-robust standard errors are in parentheses.2SLS, two-stage least squares; FE, fixed effect; K-P Wald stat., Kleibergen-Paap Wald statistic; OLS, ordinary least squares.
Table 2
2SLS estimates by demographic groups
| (1) | (2) | (3) | (4) | (5) | |
|---|---|---|---|---|---|
| ln(hourly wage) | ln(hours/week) | Unemployment | LFP | K-P Wald stat. | |
| All | −1.05*** (0.22) | 0.27 (0.27) | −0.08 (0.07) | −0.21* (0.11) | 23.35 |
| Male | −1.16*** (0.22) | 0.37** (0.17) | −0.05 (0.08) | −0.18** (0.08) | 24.34 |
| Female | −0.91*** (0.25) | 0.19 (0.39) | −0.12 (0.08) | −0.22* (0.13) | 22.22 |
| Age 15–24 years | −0.94*** (0.20) | 0.14 (0.46) | −0.01 (0.12) | −0.54*** (0.18) | 21.34 |
| Age 25–34 years | −0.97*** (0.22) | 0.31 (0.21) | −0.13 (0.10) | −0.08 (0.07) | 22.11 |
| Age 35–44 years | −1.10*** (0.23) | 0.29 (0.20) | −0.06 (0.07) | −0.07 (0.06) | 23.39 |
| Age 45–54 years | −1.01*** (0.26) | 0.34 (0.25) | −0.10** (0.05) | −0.11 (0.10) | 26.11 |
| Age 55–64 years | −1.21*** (0.30) | 0.24 (0.38) | −0.02 (0.05) | −0.10 (0.11) | 23.71 |
| Less than secondary | −1.42*** (0.23) | 0.50 (0.34) | −0.05 (0.06) | −0.23* (0.12) | 28.72 |
| Secondary | −0.86*** (0.22) | 0.32 (0.21) | −0.02 (0.09) | −0.12** (0.06) | 25.64 |
| Postsecondary | −0.75* (0.38) | −0.02 (0.23) | −0.18* (0.10) | −0.24 (0.16) | 17.22 |
Notes: Outcomes are residualized and multiplied by 100. All models include Year FE and City FE. Observations are weighted by city–year–group population. Cluster-robust standard errors are in parentheses. FE, fixed effect; K-P Wald stat., Kleibergen-Paap Wald statistic; LFP, labor force participation.

Figure 3
The effect of migration on residual ln(hourly wage) is separately estimated via 2SLS within occupation skill groups.
Note: 95% confidence intervals are presented around each point estimate.
Table 3
2SLS estimates’ robustness
| (1) | (2) | (3) | (4) | (5) | (6) | |
|---|---|---|---|---|---|---|
| Original 2SLS | Drop <100 km from border | Control trade with Venezuela | Year trend × inverse distance to border | Year trend × region | Year trend × pre-trend | |
| ln(hourly wage) | ||||||
| All | −1.05*** (0.22) | −1.02 (0.68) | −1.07*** (0.21) | −0.53 (0.67) | −0.59** (0.28) | −1.06*** (0.22) |
| Less than secondary | −1.42*** (0.23) | −1.43** (0.64) | −1.45*** (0.22) | −0.45 (0.57) | −0.57 (0.40) | −1.42*** (0.23) |
| Secondary | −0.86*** (0.22) | −0.93 (0.68) | −0.86*** (0.22) | −0.54 (0.58) | −0.32 (0.23) | −0.87*** (0.22) |
| Postsecondary | −0.75* (0.38) | −0.65 (0.99) | −0.76** (0.35) | −0.33 (1.18) | −0.71** (0.30) | −0.82** (0.36) |
| ln(hours/week) | ||||||
| All | 0.27 (0.27) | −0.50* (0.28) | 0.27 (0.28) | −0.68*** (0.24) | 0.47* (0.24) | 0.28 (0.26) |
| Less than secondary | 0.50 (0.34) | −0.58 (0.37) | 0.51 (0.35) | −0.74** (0.31) | 0.67** (0.30) | 0.52 (0.34) |
| Secondary | 0.32 (0.21) | −0.29 (0.28) | 0.32 (0.21) | −0.37 (0.25) | 0.50** (0.23) | 0.32* (0.17) |
| Postsecondary | −0.02 (0.23) | −0.52* (0.28) | −0.04 (0.25) | −0.80*** (0.22) | 0.16 (0.22) | −0.02 (0.22) |
| Unemployment | ||||||
| All | −0.08 (0.07) | −0.15 (0.20) | −0.09 (0.07) | −0.38** (0.16) | −0.25* (0.13) | −0.12 (0.08) |
| Less than secondary | −0.05 (0.06) | −0.13 (0.17) | −0.06 (0.06) | −0.34** (0.13) | −0.17* (0.10) | −0.07 (0.07) |
| Secondary | −0.02 (0.09) | −0.21 (0.19) | −0.02 (0.08) | −0.41*** (0.16) | −0.21 (0.14) | −0.04 (0.08) |
| Postsecondary | −0.18* (0.10) | −0.21 (0.24) | −0.18** (0.09) | −0.47** (0.24) | −0.37** (0.18) | −0.19* (0.10) |
| LFP | ||||||
| All | −0.21* (0.11) | 0.07 (0.18) | −0.20* (0.11) | 0.05 (0.23) | −0.20** (0.09) | −0.13 (0.09) |
| Less than secondary | −0.23* (0.12) | 0.09 (0.23) | −0.23* (0.12) | 0.04 (0.31) | −0.20* (0.11) | −0.17 (0.11) |
| Secondary | −0.12** (0.06) | −0.02 (0.15) | −0.09 (0.07) | 0.00 (0.17) | −0.14*** (0.05) | −0.11 (0.07) |
| Postsecondary | −0.24 (0.16) | 0.13 (0.16) | −0.25 (0.16) | 0.17 (0.19) | −0.26** (0.13) | −0.21 (0.15) |
| K-P Wald stat. | 23.35 | 13.07 | 28.14 | 13.95 | 88.40 | 23.61 |
| Number of metro areas | 79 | 73 | 79 | 79 | 79 | 79 |
Notes: Outcomes are residualized and multiplied by 100. All models include Year FE and City FE. See text for description of robustness checks. Observations are weighted by city–year–group population. Cluster-robust standard errors are in parentheses. FE, fixed effect; K-P Wald stat., Kleibergen-Paap Wald statistic; LFP, labor force participation.
Table 4
2SLS effects on employment by occupation skill group
| (1) | (2) | (3) | (4) | (5) | (6) | (7) | |
|---|---|---|---|---|---|---|---|
| Occupation | Underemployed | K−P Wald stat. | |||||
| Group 1 | Group 2 | Group 3 | Group 4 | Group 5 | |||
| All | −0.06*** (0.02) | −0.03 (0.08) | 0.07 (0.04) | −0.01 (0.02) | −0.09*** (0.02) | 0.29** (0.13) | 25.69 |
| Female | 0.02 (0.03) | −0.05 (0.05) | 0.02 (0.06) | −0.03 (0.02) | −0.07*** (0.02) | 0.24* (0.12) | 25.35 |
| Male | −0.15*** (0.03) | −0.00 (0.10) | 0.12*** (0.02) | 0.02 (0.03) | −0.11*** (0.02) | 0.35** (0.14) | 26.10 |
| Age 15–24 years | −0.16*** (0.03) | −0.19** (0.07) | −0.01 (0.07) | −0.04 (0.05) | −0.01 (0.02) | 0.15 (0.11) | 24.77 |
| Age 25–34 years | −0.11*** (0.03) | 0.01 (0.13) | 0.11** (0.05) | 0.15*** (0.04) | −0.11** (0.05) | 0.49*** (0.16) | 23.15 |
| Age 35–44 years | 0.00 (0.03) | 0.08 (0.09) | 0.08** (0.04) | −0.03 (0.04) | −0.18*** (0.04) | 0.36** (0.14) | 25.02 |
| Age 45–54 years | −0.04 (0.05) | 0.07 (0.05) | 0.12*** (0.04) | −0.08* (0.04) | −0.04 (0.03) | 0.35** (0.15) | 29.18 |
| Age 55–64 years | 0.09** (0.04) | −0.12 (0.09) | 0.02 (0.08) | −0.05 (0.03) | −0.06*** (0.02) | 0.09 (0.08) | 29.02 |
| Less than secondary | −0.09*** (0.03) | −0.04 (0.07) | −0.01 (0.06) | −0.02* (0.01) | −0.01** (0.00) | 0.18 (0.11) | 31.28 |
| Secondary | −0.14*** (0.03) | −0.04 (0.06) | 0.15*** (0.05) | −0.00 (0.05) | 0.00 (0.01) | 0.28* (0.15) | 26.91 |
| Postsecondary | −0.02 (0.02) | −0.02 (0.09) | 0.12** (0.05) | 0.04 (0.04) | −0.20*** (0.07) | 0.49*** (0.17) | 19.22 |
Notes: All models are 2SLS linear probability models for probability of employment in a group (multiplied by 100), not conditional on working, with Year FE and City FE. Occupations are ranked according to mean education of natives pre-2015 and grouped into quintiles. Workers are underemployed if they say they would like to change jobs to improve use of skills or training. Observations are weighted by city–year–group population. Cluster-robust standard errors are in parentheses.2SLS, two-stage least squares; FE, fixed effect; K-P Wald stat., Kleibergen-Paap Wald statistic.
Table 5
2SLS effects on employment by type of work
| (1) | (2) | (3) | (4) | (5) | |
|---|---|---|---|---|---|
| Formal salaried | Informal salaried | Own account | Employer | K-P Wald stat. | |
| All | −0.09*** (0.02) | −0.02 (0.07) | −0.04 (0.09) | 0.01 (0.02) | 25.69 |
| Female | −0.09*** (0.03) | −0.00 (0.06) | −0.03 (0.10) | 0.01 (0.01) | 25.35 |
| Male | −0.08** (0.03) | −0.03 (0.10) | −0.05 (0.09) | 0.02 (0.03) | 26.10 |
| Age 15–24 years | −0.10*** (0.03) | −0.02 (0.08) | −0.31*** (0.11) | −0.01 (0.00) | 24.77 |
| Age 25–34 years | −0.02 (0.05) | −0.01 (0.08) | 0.07 (0.11) | −0.00 (0.02) | 23.15 |
| Age 35–44 years | −0.09** (0.04) | −0.04 (0.08) | 0.10 (0.09) | 0.04 (0.03) | 25.02 |
| Age 45–54 years | −0.10** (0.04) | −0.00 (0.06) | 0.02 (0.08) | 0.03 (0.03) | 29.18 |
| Age 55–64 years | −0.10** (0.05) | −0.03 (0.05) | −0.01 (0.11) | 0.05** (0.02) | 29.02 |
| Less than secondary | −0.10*** (0.02) | −0.05 (0.08) | −0.06 (0.10) | 0.02 (0.03) | 31.28 |
| Secondary | −0.05 (0.04) | −0.07 (0.07) | 0.06 (0.11) | −0.00 (0.02) | 26.91 |
| Postsecondary | −0.08* (0.05) | 0.01 (0.06) | −0.04 (0.10) | 0.03 (0.02) | 19.22 |
Notes: Models are 2SLS linear probability models for probability of employment in a group (multiplied by 100), not conditional on working, with Year FE and City FE. Wage-sector workers are split by formality, defined according to compliance with mandatory health and pension schemes. Self-employed workers are split into own account (with no employees, predominantly informal) and employer. Observations are weighted by city–year–group population. Cluster-robust standard errors are in parentheses. 2SLS, two-stage least squares; FE, fixed effect; K-P Wald stat., Kleibergen-Paap Wald statistic.
Table 6
Effects on native internal migration
| (1) | (2) | (3) | (4) | (5) | (6) | (7) | |
|---|---|---|---|---|---|---|---|
| Outmigration | Inmigration | ||||||
| OLS | 2SLS | Sample mean | OLS | 2SLS | Sample mean | K-P Wald stat. | |
| All | 0.06 (0.06) | 0.03 (0.03) | 6.58 | −0.04 (0.25) | −0.08 (0.11) | 10.17 | 17.78 |
| Male | 0.08 (0.06) | 0.06 (0.04) | 6.76 | 0.01 (0.26) | −0.04 (0.12) | 10.16 | 18.13 |
| Female | 0.05 (0.06) | 0.01 (0.03) | 6.42 | −0.07 (0.23) | −0.11 (0.11) | 10.17 | 17.48 |
| Age 15–24 years | 0.07 (0.08) | 0.05 (0.06) | 7.67 | −0.02 (0.32) | −0.06 (0.16) | 14.88 | 16.69 |
| Age 25–34 years | 0.05 (0.08) | −0.03 (0.05) | 9.49 | −0.09 (0.32) | −0.13 (0.15) | 13.02 | 16.74 |
| Age 35–44 years | 0.05 (0.09) | 0.03 (0.05) | 6.64 | −0.06 (0.25) | −0.10 (0.11) | 9.03 | 17.02 |
| Age 45–54 years | 0.07 (0.05) | 0.07** (0.03) | 3.97 | −0.04 (0.14) | −0.09 (0.06) | 5.62 | 20.18 |
| Age 55–64 years | 0.04 (0.03) | 0.05* (0.03) | 3.09 | −0.04 (0.13) | −0.09 (0.07) | 4.24 | 20.09 |
| Less than secondary | 0.10* (0.06) | 0.07** (0.04) | 5.54 | −0.01 (0.25) | −0.04 (0.12) | 8.87 | 21.90 |
| Secondary | 0.08 (0.08) | 0.10*** (0.03) | 6.62 | −0.08 (0.23) | −0.17 (0.11) | 10.04 | 17.97 |
| Postsecondary | 0.00 (0.08) | −0.07 (0.05) | 7.70 | −0.02 (0.28) | −0.01 (0.14) | 11.71 | 13.30 |
Notes: Sample includes the years 2014 and 2019. In- and outmigration rates are taken as shares of the 5-year lag population and multiplied by 100. All models include Year FE and City FE. Observations are weighted by city–year–group population. Cluster-robust standard errors are in parentheses.2SLS, two-stage least squares; FE, fixed effect; K-P Wald stat., Kleibergen-Paap Wald statistic; OLS, ordinary least squares.
Table 7
2SLS wage effects interacted with regional characteristics
| ln(hourly wage) | |||||||||
|---|---|---|---|---|---|---|---|---|---|
| (1) | (2) | (3) | (4) | (5) | (6) | (7) | (8) | (9) | |
| Migrant share | −0.97*** (0.26) | −1.14** (0.55) | −1.02*** (0.23) | −0.55 (0.39) | −1.28*** (0.23) | −0.94*** (0.30) | −1.23*** (0.20) | −1.01*** (0.37) | −1.01*** (0.30) |
| Migrant share interacted with: 2014 mean ln(hourly wage) | −0.36 (0.46) | ||||||||
| 2014 unemployment rate | 0.07 (0.29) | ||||||||
| 2014 informal rate | −0.62* (0.34) | ||||||||
| 2014 own-account rate | −0.75** (0.35) | ||||||||
| 2017 WB DB (starting a business) | 0.29* (0.16) | ||||||||
| 2017 WB DB (construction permits) | −0.31 (0.50) | ||||||||
| 2017 WB DB (registering property) | −0.36 (0.43) | ||||||||
| 2017 WB DB (paying taxes) | −0.06 (0.23) | ||||||||
| 2014 per capita GDP | −0.13 (0.32) | ||||||||
| K-P Wald stat. | 12.38 | 9.51 | 14.11 | 13.27 | 36.46 | 11.55 | 8.73 | 9.03 | 13.31 |
| N | 474 | 474 | 474 | 474 | 474 | 474 | 474 | 474 | 474 |
Notes: Both the endogenous variable and the instrument are interacted with the interaction variable indicated in each row. All interaction variables are centered around zero and divided by the standard error. 2014 metro-level economic indicators are measured using the GEIH. Doing Business measures at the department level were taken from the World Bank 2017 Doing Business report. The values of 2014 per capita GDP at the department level downloaded from DANE. All models include Year FE and City FE. Observations are weighted by city–year population. Cluster-robust standard errors are in parentheses. DANE, National Administrative Department of Statistics; FE, fixed effect; GDP, gross domestic product; GEIH, Colombian National Integrated Household Survey; K-P Wald stat., Kleibergen-Paap Wald statistic; WB DB, World Bank Doing Business.
Table A1
Characteristics of migrants and nonmigrants
| Nonmigrants | Migrants | |
|---|---|---|
| Male (%) | 48.5 (50) | 49.6 (50) |
| Age (years) | 36.5 (13.9) | 31.5 (11.3) |
| Labor force participation (%) | 71.7 (45) | 79.4 (40.5) |
| Unemployment (%) | 11.4 (31.7) | 14.8 (35.5) |
| Median hourly wage (2010 USD) | 2.3 (6.1) | 1.6 (4.8) |
| Hours/week | 45.2 (15.9) | 49.6 (17.4) |
| Own account (%) | 25.3 (43.5) | 32.1 (46.7) |
| Informal (%) | 56 (49.6) | 88.3 (32.1) |
| N | 447,264 | 21,730 |
[i] Notes: Means are presented, SDs are in parentheses; data are restricted to urban residents of age 15–64 years. Migrant is defined as anyone in 2019 who was living in Venezuela 5 years ago. Population weights are applied. GEIH, Colombian National Integrated Household Survey; SD, standard deviation.
[ii] Source: GEIH (2019).
Table A2
Preperiod correlates of 2005 Venezuelan share of population
| (1) | (2) | (3) | (4) | (5) | |
|---|---|---|---|---|---|
| Population (100,000) | −0.001 (0.005) | 0.005 (0.003) | 0.004 (0.003) | 0.004 (0.003) | 0.002 (0.003) |
| 100*ln(hourly wage) | −0.017** (0.008) | −0.018* (0.009) | −0.018* (0.009) | −0.024* (0.026) | |
| 100*ln(hours/week) | 0.021 (0.017) | 0.021 (0.017) | 0.026 (0.019) | ||
| 100* unemployment rate | 0.047 (0.125) | 0.049 (0.126) | |||
| 100* LPP rate | 0.041 (0.047) | ||||
| N | 79 | 79 | 79 | 79 | 79 |
| R2 | 1.6e−04 | 0.051 | 0.073 | 0.073 | 0.09 |
| 2014 metro characteristic | Mc,2005 | Mc,2005 | Mc,2005 | Mc,2005 | Mc,2005 |
Notes: Mc,2005 is the 2005 Venezuelan share of the metro area, which ranges from 0% to 0.9%, and has been normalized to a mean of zero and SD of one. Observations are weighted by 2014 population. Robust standard errors are in parentheses. LPP, Labor force participation; SD, standard deviation.
Table A3
Wage estimate sensitivity
| Panel A: including return migrants | ||||
|---|---|---|---|---|
| (1) | (2) | (3) | (4) | |
| OLS | IV: 2005 census | IV: 1993 census | IV: inverse distance | |
| Geographic unit: CZ metro areas | ||||
| Migrant share | −0.73* (0.42) | −1.05*** (0.22) | −1.13*** (0.21) | −1.16*** (0.17) |
| K-P Wald stat. | 23.35 | 21.78 | 28.54 | |
| Number of units | 79 | 79 | 79 | 79 |
| Geographic unit: administrative metro areas | ||||
| Migrant share | −0.67 (0.49) | −0.94*** (0.26) | −1.03*** (0.20) | −1.08*** (0.16) |
| K-P Wald stat. | 29.58 | 19.88 | 46.11 | |
| Number of units | 23 | 23 | 23 | 23 |
| Geographic unit: department | ||||
| Migrant share | −0.83 (0.51) | −1.19*** (0.28) | −1.21*** (0.30) | −1.27*** (0.18) |
| K-P Wald stat. | 19.07 | 20.90 | 35.20 | |
| Number of units | 24 | 24 | 24 | 24 |
| Panel B: excluding return migrants | ||||
|---|---|---|---|---|
| (5) | (6) | (7) | (8) | |
| OLS | IV: 2005 census | IV: 1993 census | IV: inverse distance | |
| Geographic unit: CZ metro areas | ||||
| Migrant share | −0.87** (0.42) | −1.24*** (0.26) | −1.32*** (0.25) | −1.34*** (0.21) |
| K-P Wald stat. | 26.42 | 21.52 | 29.50 | |
| Number of units | 79 | 79 | 79 | 79 |
| Geographic unit: administrative metro areas | ||||
| Migrant share | −0.79 (0.48) | −1.09*** (0.30) | −1.18*** (0.23) | −1.24*** (0.18) |
| K-P Wald stat. | 30.17 | 17.76 | 39.35 | |
| Number of units | 23 | 23 | 23 | 23 |
| Geographic unit: department | ||||
| Migrant share | −0.97* (0.52) | −1.41*** (0.34) | −1.44*** (0.35) | −1.48*** (0.21) |
| K-P Wald stat. | 19.55 | 19.38 | 31.37 | |
| Number of units | 24 | 24 | 24 | 24 |
Notes: The outcome, namely, log hourly wage, is residualized and multiplied by 100. All models include Year FE and Unit FE. Panel A includes Colombian-born migrants in the migrant share, while Panel B excludes them. Columns 2–4 use different share components of the shift-share IV, based respectively on the complete 2005 census (as in the main analysis), the 10% subsample of the 1993 census, and the inverse minimum driving distance from the metro area or department capital to the closest Venezuelan border crossing. Geographic units include the 79 CZ-defined metro areas (as in the main analysis), the 23 official metro areas of which the GEIH is representative, and the 24 departments of Colombia. Observations are weighted by unit–year population. Cluster-robust standard errors are in parentheses. CZ, commuting zone; FE, fixed effect; GEIH, Colombian National Integrated Household Survey; IV, instrumental variable; K-P Wald stat., Kleibergen-Paap Wald statistic; OLS, ordinary least squares.
Table A4
Wage estimate sensitivity: 2014–2018 difference model
| Panel A: including return migrants | ||||
|---|---|---|---|---|
| (1) | (2) | (3) | (4) | |
| OLS | IV: 2005 census | IV: 1993 census | IV: inverse distance | |
| Geographic unit: CZ metro areas | ||||
| Foreigner share (2018 GEIH) | −0.26 (0.60) | −1.07*** (0.39) | −1.05*** (0.35) | −1.24*** (0.24) |
| Foreigner share (2018 census) | −0.54 (0.49) | −1.26*** (0.46) | −1.21*** (0.42) | −1.45*** (0.28) |
| Number of units | 79 | 79 | 79 | 79 |
| Geographic unit: administrative metro areas | ||||
| Foreigner share (2018 GEIH) | −0.09 (0.63) | −0.87** (0.43) | −0.89*** (0.33) | −1.07*** (0.23) |
| Foreigner share (2018 census) | −0.39 (0.54) | −1.13** (0.57) | −1.09** (0.43) | −1.35*** (0.29) |
| Number of units | 23 | 23 | 23 | 23 |
| Geographic unit: department | ||||
| Foreigner share (2018 GEIH) | −0.25 (0.59) | −1.02*** (0.39) | −0.97** (0.41) | −1.19*** (0.23) |
| Foreigner share (2018 census) | −0.46 (0.79) | −1.36*** (0.50) | −1.32** (0.53) | −1.53*** (0.29) |
| Number of units | 24 | 24 | 24 | 24 |
| Panel B: excluding return migrants | ||||
|---|---|---|---|---|
| (5) | (6) | (7) | (8) | |
| OLS | IV: 2005 census | IV: 1993 census | IV: inverse distance | |
| Geographic unit: CZ metro areas | ||||
| Foreigner share (2018 GEIH) | −0.62 (0.75) | −1.58*** (0.58) | −1.53*** (0.52) | −1.80*** (0.34) |
| Foreigner share (2018 census) | −0.65 (0.55) | −1.55*** (0.56) | −1.46*** (0.51) | −1.77*** (0.34) |
| Number of units | 79 | 79 | 79 | 79 |
| Geographic unit: administrative metro areas | ||||
| Foreigner share (2018 GEIH) | −0.36 (0.76) | −1.27** (0.63) | −1.26*** (0.48) | −1.54*** (0.32) |
| Foreigner share (2018 census) | −0.48 (0.61) | −1.38** (0.70) | −1.31** (0.53) | −1.65*** (0.36) |
| Number of units | 23 | 23 | 23 | 23 |
| Geographic unit: department | ||||
| Foreigner share (2018 GEIH) | −0.55 (0.75) | −1.54*** (0.59) | −1.43** (0.61) | −1.77*** (0.33) |
| Foreigner share (2018 census) | −0.55 (0.99) | −1.81*** (0.67) | −1.72** (0.70) | −2.02*** (0.38) |
| Number of units | 24 | 24 | 24 | 24 |
Notes: See notes to Table A3. This model is a regression of the change in residual wages from 2014 to 2018 on the 2018 foreigner share, calculated using the GEIH or Census as a share of the 2018 population. The instrument is created analogously (interacted with the national migration rate from the GEIH or census). In no case does the K-P Wald statistic fall below 10. CZ, commuting zone; FE, fixed effect; GEIH, Colombian National Integrated Household Survey; IV, instrumental variable; K-P Wald stat., Kleibergen-Paap Wald statistic; OLS, ordinary least squares.
Table A5
Wage estimate sensitivity: 1-year migration measure
| Panel A: including return migrants | ||||
|---|---|---|---|---|
| (1) | (2) | (3) | (4) | |
| OLS | IV: 2005 census | IV: 1993 census | IV: inverse distance | |
| Geographic unit: CZ metro areas | ||||
| 1-year migrant share (data until 2019) | −1.45 (1.29) | −3.53*** (0.80) | −3.96*** (0.90) | −3.82*** (0.67) |
| 1-year migrant share (data until 2017) | −2.18 (2.35) | −4.80*** (1.36) | −5.15*** (1.25) | −5.74*** (0.91) |
| Number of units | 79 | 79 | 79 | 79 |
| Geographic unit: administrative metro areas | ||||
| 1-year migrant share (data until 2019) | −1.38 (1.67) | −3.05*** (0.89) | −3.60*** (0.82) | −3.58*** (0.59) |
| 1-year migrant share (data until 2017) | −2.76 (3.03) | −4.40** (1.74) | −5.03*** (1.35) | −5.65*** (0.96) |
| Number of units | 23 | 23 | 23 | 23 |
| Geographic unit: department | ||||
| 1-year migrant share (data until 2019) | −1.80 (1.65) | −4.01*** (1.01) | −4.31*** (1.19) | −4.17*** (0.63) |
| 1-year migrant share (data until 2017) | −3.58 (3.15) | −5.38*** (1.59) | −5.29*** (1.63) | −6.33*** (0.92) |
| Number of units | 24 | 24 | 24 | 24 |
| Panel B: excluding return migrants | ||||
|---|---|---|---|---|
| (5) | (6) | (7) | (8) | |
| OLS | IV: 2005 census | IV: 1993 census | IV: inverse distance | |
| Geographic unit: CZ metro areas | ||||
| 1-year migrant share (data until 2019) | −1.74 (1.29) | −3.75*** (0.85) | −4.28*** (1.06) | −3.94*** (0.75) |
| 1-year migrant share (data until 2017) | −3.64 (2.85) | −5.51*** (1.48) | −5.92*** (1.30) | −6.64*** (0.96) |
| Number of units | 79 | 79 | 79 | 79 |
| Geographic unit: administrative metro areas | ||||
| 1-year migrant share (data until 2019) | −1.58 (1.57) | −3.03*** (0.87) | −3.55*** (0.79) | −3.55*** (0.59) |
| 1-year migrant share (data until 2017) | −3.78 (3.23) | −4.85** (1.88) | −5.49*** (1.40) | −6.42*** (1.02) |
| Number of units | 23 | 23 | 23 | 23 |
| Geographic unit: department | ||||
| 1-year migrant share (data until 2019) | −1.97 (1.63) | −4.14*** (1.03) | −4.45*** (1.21) | −4.30*** (0.67) |
| 1-year migrant share (data until 2017) | −4.82 (3.41) | −6.68*** (1.93) | −6.50*** (1.93) | −8.07*** (1.08) |
| Number of units | 24 | 24 | 24 | 24 |
Notes: See notes to Table A3. This is the original two-way fixed-effects model, replacing the 5-year migrant share with the 1-year migrant share (the share of the population that arrived from Venezuela in the past 12 months) and varying the end of the period of analysis. In no case does the K-P Wald stat. fall below 10.
CZ, commuting zone; FE, fixed effect; GEIH, Colombian National Integrated Household Survey; IV, instrumental variable; K-P Wald stat., Kleibergen-Paap Wald statistic; OLS, ordinary least squares.
Table A6
2SLS estimates additional robustness
| (1) | (2) | (3) | (4) | (5) | (6) | |
|---|---|---|---|---|---|---|
| Original 2SLS | Drop Bogotá | Drop metro yearly sample <1,000 | 2014–2019 Difference model | Kronmal specification | Assign natives to past metro | |
| ln(hourly wage) | ||||||
| All | −1.05*** (0.22) | −1.03*** (0.22) | −1.04*** (0.23) | −1.09*** (0.27) | −0.66*** (0.15) | −1.04*** (0.20) |
| Less than secondary | −1.42*** (0.23) | −1.41*** (0.24) | −1.41*** (0.25) | −1.33*** (0.27) | −0.95*** (0.17) | −1.38*** (0.22) |
| Secondary | −0.86*** (0.22) | −0.87*** (0.22) | −0.85*** (0.23) | −0.79*** (0.28) | −0.54*** (0.15) | −0.80*** (0.21) |
| Postsecondary | −0.75* (0.38) | −0.69* (0.37) | −0.75* (0.39) | −0.96** (0.48) | −0.46** (0.23) | −0.80** (0.36) |
| ln(hours/week) | ||||||
| All | 0.27 (0.27) | 0.28 (0.27) | 0.27 (0.27) | 0.27 (0.32) | 0.17 (0.17) | 0.29 |
| Less than secondary | 0.50 (0.34) | 0.52 (0.34) | 0.51 (0.34) | 0.51 (0.39) | 0.34 (0.22) | 0.51 (0.34) |
| Secondary | 0.32 (0.21) | 0.32 (0.21) | 0.30 (0.21) | 0.26 (0.26) | 0.20 (0.12) | 0.32 (0.21) |
| Postsecondary | −0.02 (0.23) | −0.02 (0.22) | −0.00 (0.22) | 0.01 (0.25) | −0.01 (0.12) | 0.02 (0.22) |
| Unemployment | ||||||
| All | −0.08 (0.07) | −0.08 (0.07) | −0.09 (0.08) | −0.07 (0.07) | −0.05 (0.04) | −0.07 (0.07) |
| Less than secondary | −0.05 (0.06) | −0.05 (0.07) | −0.06 (0.07) | −0.03 (0.06) | −0.04 (0.04) | −0.06 (0.06) |
| Secondary | −0.02 (0.09) | −0.01 (0.08) | −0.03 (0.09) | 0.02 (0.09) | −0.01 (0.05) | −0.01 (0.08) |
| Postsecondary | −0.18* (0.10) | −0.18* (0.10) | −0.17* (0.10) | −0.20* (0.10) | −0.11* (0.06) | −0.15 (0.10) |
| LFP | ||||||
| All | −0.21* (0.11) | −0.21** (0.11) | −0.21** (0.10) | −0.17* (0.10) | −0.13** (0.07) | −0.20** (0.10) |
| Less than secondary | −0.23* (0.12) | −0.24* (0.12) | −0.24** (0.11) | −0.16 (0.10) | −0.16* (0.08) | −0.24** (0.12) |
| Secondary | −0.12** (0.06) | −0.12** (0.06) | −0.12** (0.06) | −0.10 (0.07) | −0.07** (0.04) | −0.12** (0.05) |
| Postsecondary | −0.24 (0.16) | −0.25 (0.16) | −0.23 (0.17) | −0.27 (0.18) | −0.15* (0.08) | −0.21 (0.15) |
| K-P Wald stat. | 23.35 | 23.41 | 19.28 | 30.89 | 24.17 | 24.47 |
| Number of metro areas | 79 | 78 | 27 | 79 | 79 | 79 |
Notes: Outcomes are residualized and multiplied by 100. All models include Year FE and City FE. Column 4 regresses the change in outcomes from 2014 to 2019 on the change in the migrant share. Column 5 replaces the migrant share with the number of migrants (divided by 100) and controls for baseline population interacted with year FEs. Column 6 assigns natives to their 5-year lagged metro area of residence. Observations are weighted by city–year–group population. Cluster-robust standard errors are in parentheses. 2SLS, two-stage least squares; FE, fixed effect; K-P Wald stat., Kleibergen-Paap Wald statistic; LFP, labor force participation.
Table A7
2SLS estimates with migrant share quadratic
| (1) | (2) | (3) | (4) | |
|---|---|---|---|---|
| ln(hourly wage) | ||||
| All | Less than secondary | Secondary | Postsecondary | |
| Migrant share | −1.40 (1.12) | −2.37** (1.11) | −1.14 (1.02) | −0.79 (1.67) |
| Migrant share2 | 0.02 (0.05) | 0.05 (0.05) | 0.02 (0.05) | 0.00 (0.08) |
| Joint F-statistic | 62.36 | 55.46 | 39.73 | 25.46 |
| K-P Wald stat. | 58.02 | 35.66 | 33.97 | 89.30 |
Notes: Linear and quadratic migrant shares are instrumented with the linear and quadratic forms of the IV. Wages are residualized and multiplied by 100. All models include Year FE and City FE. Observations are weighted by city–year–group population. Cluster-robust standard errors are in parentheses. 2SLS, two-stage least squares; FE, fixed effect; K-P Wald stat., Kleibergen-Paap Wald statistic.
Table A8
Robustness check of 2SLS effects on occupation skill group and type of work
| (1) | (2) | (3) | (4) | (5) | (6) | (7) | (8) | (9) | (10) | (11) | |
|---|---|---|---|---|---|---|---|---|---|---|---|
| Original 2SLS | Drop <100 km from border | Control trade with Venezuela | Year trend × inverse distance to border | Year trend × region | Year trend × pre-trend | Drop Bogotá | Drop metro yearly sample <1,000 | 2014–2019 difference model | Kronmal specification | Assign natives to past metro | |
| Occup skill Group 1 | −0.06*** (0.02) | −0.05 (0.08) | −0.05** (0.03) | 0.04 (0.07) | −0.00 (0.04) | −0.06** (0.02) | −0.07*** (0.02) | −0.07*** (0.02) | −0.06*** (0.02) | −0.04*** (0.02) | −0.06** (0.02) |
| Occup skill Group 2 | −0.03 (0.08) | 0.15 (0.10) | −0.05 (0.07) | 0.21** (0.10) | 0.03 (0.08) | −0.01 (0.08) | −0.03 (0.08) | −0.03 (0.08) | −0.04 (0.07) | −0.02 (0.04) | −0.03 (0.08) |
| Occup skill Group 3 | 0.07 (0.04) | 0.16** (0.08) | 0.08 (0.05) | 0.21*** (0.07) | 0.05 (0.05) | 0.09** (0.04) | 0.06 (0.04) | 0.06 (0.04) | 0.08* (0.05) | 0.04 (0.03) | 0.06 (0.04) |
| Occup skill Group 4 | −0.01 (0.02) | 0.03 (0.05) | 0.00 (0.03) | 0.04 (0.06) | 0.02 (0.02) | −0.01 (0.02) | −0.01 (0.02) | −0.00 (0.02) | −0.00 (0.02) | −0.01 (0.02) | −0.01 (0.02) |
| Occup skill Group 5 | −0.09*** (0.02) | −0.11* (0.06) | −0.09*** (0.02) | −0.14** (0.07) | −0.11*** (0.04) | −0.09*** (0.02) | −0.08*** (0.02) | −0.09*** (0.02) | −0.09*** (0.03) | −0.06*** (0.02) | −0.09*** (0.02) |
| Underemployed | 0.29** (0.13) | 0.52 (0.32) | 0.26** (0.12) | 0.57* (0.32) | 0.31** (0.12) | 0.29** (0.13) | 0.28** (0.13) | 0.30** (0.13) | 0.33** (0.13) | 0.18** (0.08) | 0.28** (0.12) |
| Formal salaried | −0.09*** (0.02) | −0.10** (0.05) | −0.08*** (0.03) | −0.07 (0.04) | −0.09* (0.05) | −0.09*** (0.02) | −0.08*** (0.02) | −0.09*** (0.02) | −0.07*** (0.03) | −0.06*** (0.01) | −0.06*** (0.02) |
| Informal salaried | −0.02 (0.07) | 0.07 (0.16) | 0.01 (0.07) | 0.18 (0.15) | 0.09 (0.06) | −0.04 (0.06) | −0.03 (0.07) | −0.01 (0.07) | −0.07 (0.08) | −0.01 (0.05) | 0.00 (0.06) |
| Own account | −0.04 (0.09) | 0.09 (0.21) | −0.06 (0.08) | 0.13 (0.22) | −0.08 (0.13) | −0.04 (0.09) | −0.04 (0.09) | −0.05 (0.09) | 0.01 (0.07 | −0.03 (0.05) | −0.07 (0.09) |
| Employer | 0.01 (0.02) | 0.04 (0.05) | 0.02 (0.02) | 0.07 (0.04) | 0.06*** (0.02) | 0.02 (0.02) | 0.01 (0.02) | 0.02 (0.02) | 0.02 (0.02 | 0.01 (0.01) | 0.01 (0.02) |
| Kleibergen-Paap Wald stat. | 25.69 | 12.80 | 30.99 | 14.50 | 102.73 | 24.59 | 25.77 | 20.69 | 33.17 | 26.38 | 26.37 |
| Number of metro areas | 79 | 73 | 79 | 79 | 79 | 79 | 78 | 27 | 79 | 79 | 79 |
Notes: Outcomes are multiplied by 100. All models are 2SLS linear probability models for probability of employment in a group, not conditional on working, with Year FE and City FE. Observations are weighted by city–year–group population. See paper for description of robustness checks. Cluster-robust standard errors are in parentheses. 2SLS, two-stage least squares; FE, fixed effect; Occup, occupational.
Table A9
Robustness of 2SLS effects on native internal migration
| (1) | (2) | (3) | (4) | (5) | (6) | (7) | (8) | |
|---|---|---|---|---|---|---|---|---|
| Original 2SLS | Drop <100 km from border | Control trade with Venezuela | Year trend × inverse distance to border | Year trend × region | Drop Bogotá | Drop metro yearly sample <1,000 | Kronmal specification | |
| Outmigration | ||||||||
| All | 0.03 (0.03) | 0.01 (0.08) | 0.03 (0.05) | −0.09 (0.12) | −0.11** (0.05) | 0.05* (0.03) | 0.03 (0.03) | 0.02 (0.02) |
| Less than secondary | 0.07** (0.04) | 0.07 (0.10) | 0.07 (0.06) | −0.04 (0.14) | −0.09 (0.06) | 0.10*** (0.03) | 0.04 (0.03) | 0.05** (0.03) |
| Secondary | 0.10*** (0.03) | 0.05 (0.10) | 0.13*** (0.04) | 0.04 (0.13) | 0.04 (0.05) | 0.10*** (0.03) | 0.07** (0.03) | 0.06*** (0.02) |
| Postsecondary | −0.07 (0.05) | −0.06 (0.08) | −0.13 (0.08) | −0.16 (0.15) | −0.24*** (0.08) | −0.04 (0.04) | −0.04 (0.05) | −0.04 (0.03) |
| Inmigration | ||||||||
| All | −0.08 (0.11) | 0.03 (0.28) | −0.08 (0.12) | 0.26 (0.32) | −0.03 (0.09) | −0.09 (0.12) | −0.07 (0.12) | −0.05 (0.07) |
| Less than secondary | −0.04 (0.12) | 0.20 (0.29) | −0.05 (0.12) | 0.44 (0.30) | −0.01 (0.10) | −0.05 (0.12) | −0.02 (0.13) | −0.02 (0.08) |
| Secondary | −0.17 (0.11) | −0.04 (0.30) | −0.19* (0.10) | 0.11 (0.36) | −0.08 (0.08) | −0.18* (0.11) | −0.16 (0.11) | −0.11* (0.06) |
| Postsecondary | −0.01 (0.14) | −0.06 (0.30) | 0.03 (0.15) | 0.26 (0.41) | 0.07 (0.14) | −0.02 (0.14) | −0.04 (0.13) | −0.00 (0.09) |
| Kleibergen-Paap Wald stat. | 17.78 | 8.88 | 38.19 | 6.30 | 61.57 | 17.79 | 14.00 | 19.49 |
| Number of metro areas | 79 | 73 | 79 | 79 | 79 | 78 | 27 | 79 |
Notes: Sample includes the years 2014 and 2019. Migration rates are taken as a share of the 5-year lag population. All models include Year FE and City FE. Observations are weighted by city–year–group population. See paper for description of robustness checks. Cluster-robust standard errors are in parentheses. 2SLS, two-stage least squares; FE, fixed effect.
Table A10
Robustness of 2SLS wage effects interacted with regional characteristics
| Interaction variable | (1) | (2) | (3) | (4) | (5) | (6) | (7) | (8) | (9) | (10) | (11) |
|---|---|---|---|---|---|---|---|---|---|---|---|
| Original 2SLS | Drop <100 km from border | Control trade with Venezuela | Year trend × inverse distance to border | Year trend × region | Year trend × pre-trend | Drop Bogotá | Drop metro yearly sample <1,000 | 2014–2019 difference model | Kronmal specification | Assign natives to past metro | |
| 2014 mean ln(hourly wage) | −0.36 (0.46) | −0.69 (0.60) | −0.36 (0.46) | −0.31 (0.45) | −0.20† (0.42) | −0.36 (0.46) | −0.47 (0.49) | −0.68 (0.62) | −0.48† (0.50) | −0.63** (0.29) | −0.35 (0.40) |
| 2014 unemployment rate | 0.07 (0.29) | 0.35 (0.82) | 0.11 (0.31) | 0.58** (0.25) | −0.29 (0.33) | 0.08 (0.29) | 0.08 (0.30) | 0.04 (0.33) | 0.08 (0.33) | 0.11 (0.18) | 0.07 (0.26) |
| 2014 informal rate | −0.62* (0.34) | −0.80** (0.38) | −0.63*, † (0.34) | −0.56* (0.34) | −0.76*** (0.27) | −0.60* (0.34) | −0.62* (0.37) | −0.64 (0.40) | −0.62† (0.42) | −0.34 (0.48) | −0.53* (0.31) |
| 2014 own-account rate | −0.75** (0.35) | −1.15** (0.45) | −0.81**, † (0.38) | −0.73** (0.33) | −0.86***, † (0.33) | −0.73** (0.35) | −0.75** (0.37) | −0.80* (0.42) | −0.71*, † (0.43) | −0.74* (0.40) | −0.61* (0.33) |
| WB DB (starting a business) | 0.29*, † (0.16) | 0.43 (0.35) | 0.34*, † (0.18) | 0.41* (0.24) | 0.28† (0.25) | 0.29*, † (0.16) | 0.28*, † (0.16) | 0.27*, † (0.15) | 0.26† (0.16) | 0.16† (0.11) | 0.26*, † (0.14) |
| WB DB (construction permits) | −0.31 (0.50) | −0.31 (0.51) | −0.31 (0.51) | −0.29 (0.52) | −0.08† (0.38) | −0.28 (0.49) | −0.36 (0.51) | −0.30 (0.53) | −0.35 (0.61) | −0.18 (0.23) | −0.26 (0.47) |
| WB DB (registering property) | −0.36 (0.43) | −0.37 (0.43) | −0.38 (0.50) | −0.34 (0.39) | −0.38† (0.37) | −0.37 (0.43) | −0.50 (0.45) | −0.33 (0.45) | −0.34 (0.54) | −0.20 (0.16) | −0.30 (0.39) |
| WB DB (paying taxes) | −0.06 (0.23) | −0.28 (0.53) | −0.05 (0.25) | 0.30 (0.25) | −0.49* (0.27) | −0.05 (0.24) | −0.06 (0.23) | −0.07 (0.24) | −0.11 (0.27) | 0.01 (0.12) | −0.05 (0.22) |
| 2014 per capita GDP | −0.13 (0.32) | −0.15 (0.33) | −0.14† (0.31) | −0.20 (0.33) | 0.09† (0.37) | −0.11 (0.31) | −0.18 (0.30) | −0.12 (0.33) | −0.18† (0.35 | −0.14 (0.15) | −0.14 (0.30) |
| N | 474 | 438 | 474 | 474 | 474 | 474 | 468 | 162 | 79 | 474 | 474 |
Notes: Each row presents the interaction coefficient from a separate regression using the indicated interaction variable.
See notes to Table 7. Cluster-robust standard errors are in parentheses. 2SLS, two-stage least squares; GDP, gross domestic product; K-P Wald stat., Kleibergen-Paap Wald statistic; WB DB, World Bank Doing Business.
Table A11
School attendance 2SLS linear probability models
| In school | In school | In school | |
|---|---|---|---|
| In labor force | Out of labor force | ||
| All ages (including border locations) | |||
| Migrant share | 0.07** (0.03) | −0.03 (0.08) | 0.11* (0.06) |
| Sample mean | 18.76 | 7.14 | 11.62 |
| Age ≤24 years (including border locations) | |||
| Migrant share | 0.20* (0.10) | −0.19 (0.24) | 0.39* (0.23) |
| Sample mean | 54.60 | 14.43 | 40.17 |
| All ages (excluding border locations) | |||
| Migrant share | 0.13 (0.11) | 0.20* (0.10) | −0.07 (0.09) |
| Sample mean | 18.82 | 7.17 | 11.65 |
| Age ≤24 years (excluding border locations) | |||
| Migrant share | 0.28 (0.34) | 0.49* (0.27) | −0.21 (0.33) |
| Sample mean | 54.69 | 14.45 | 40.24 |
Notes: All outcomes are multiplied by 100. All models include Year FE and City FE. Observations are weighted by city–year–group population. Excluding border locations drops cities within 100 km driving distance of the Venezuelan border. Cluster-robust standard errors are in parentheses. 2SLS, two-stage least squares; FE, fixed effect.

Figure A1
Education by migrant status.
Notes: Data restricted to urban residents of age 15–64 years. GEIH, Colombian National Integrated Household Survey.
Source: GEIH (2019).

Figure A2
Concentration of migrants and natives across occupations.
Notes: Data restricted to urban residents of age 15–64 years. Occupations ranked according to mean years of completed schooling for natives in the GEIH between 2010 and 2015. GEIH, Colombian National Integrated Household Survey.
Source: GEIH (2019).

Figure A3
Relationship between 2019 and 2005 migrant shares.
Note: RMSE, root mean square error.

Figure A4
Pre-trends in labor market outcomes.
Note: Pre-trends estimated according to Eq. (5) in Section 5. They are 95% confidence intervals.

Figure A5
2SLS coefficient plot – by work sector.
Note: 2SLS, two-stage least squares. They are 95% confidence intervals.

Figure A6
LOWESS plot.
Note: LOWESS, locally weighted scatterplot smoothing. The red dotted line displays the least square line of best fit.