I. Introduction

Climate change, primarily driven by fossil fuels, contributes around 70% of annual greenhouse gas emissions and 90% of CO2 emissions (United Nations, 2022) and is being aggravated by rapid industrialization in both OECD and developing countries (Lin & Omoju, 2017). The IPCC (2022) emphasizes that the transition towards clean energy is essential to limit global warming and unlock advantages such as employment generation and enhanced energy security. This “Energy transition” – substituting fossil fuels with cleaner renewable sources like wind, solar, biofuel, hydrogen, etc. (Harichandan et al., 2022) – also influences human mobility. Hence, net migration is increasingly affected by environmental factors besides traditional socio-economic determinants (World Bank, 2018).

The interaction between energy transition and migration patterns operates through multiple pathways. Global warming driven by burning fossil fuels creates environmental pressures that can trigger migration (IPCC, 2022; IOM, 2022). Similarly, transitioning to renewable energy acts as a pull factor for migration by creating new economic opportunities and improving the standard of living (Ram et al., 2019; United Nations, 2022). Renewable energy-fed electrification also improves services, infrastructure and rural prosperity, further attracting migrants while reducing rural push factors (Bhattacharya et al., 2015; Fried & Lagakos, 2021; Melkior et al., 2018) and shaping broader migration patterns (Łukaniszyn-Domaszewska et al., 2025).

Environmental quality also affects migration decisions, as higher renewable energy consumption reduces CO2 emissions and improves livability (Koengkan et al., 2020). In climate-vulnerable regions, environmental factors are increasingly shaping migration decisions (World Bank, 2018). Nevertheless, economic factors largely drive migration, with per capita GDP positively and unemployment negatively affecting net migration (Urbański, 2022). However, the interaction between these traditional economic factors and emerging environmental determinants of migration needs further investigation.

Although scholarly attention to the interlink between renewable energy and economic development has been intensified, its impact on migration outcomes remains unexplored. This study examines how renewable energy transitions affect net migration while accounting for traditional socio-economic indicators, using the Panel Autoregressive Distributed Lag (PARDL) methodology across 61 countries from 2000-2023. The aim is to assess both short-run and long-run relationships by hypothesizing that the renewable energy transition effect on net migration is mediated by both environmental and economic factors.

The remainder of the paper is organized as follows. Section II presents the data and methodology, followed by a discussion of the main findings in Section III. The final section concludes the paper.

II. Data and Methodology

A. Data

The current study employs an annual panel dataset for 61 countries from 2000-2023, selected according to data availability. These countries include: Algeria, Argentina, Australia, Austria, Azerbaijan, Bangladesh, Belarus, Belgium, Brazil, Canada, Chile, China, Colombia, Croatia, Czechia, Denmark, Ecuador, Egypt, Estonia, Finland, France, Germany, Greece, Hungary, Iceland, India, Indonesia, Iran (Islamic Republic), Iraq, Ireland, Italy, Japan, Kazakhstan, Latvia, Lithuania, Luxembourg, Malaysia, Mexico, Morocco, Netherlands, New Zealand, North Macedonia, Norway, Pakistan, Peru, Philippines, Poland, Portugal, Romania, Slovak Republic, South Africa, Spain, Sri Lanka, Sweden, Switzerland, Thailand, Türkiye, Ukraine, United Kingdom, United States, Vietnam. The outcome variable for this study is net migration (NM), and the main explanatory variables comprise renewable energy transition indicators such as share of renewable energy consumption (SRN), renewable electricity per capita (REPC), and energy intensity (EI). Additionally, GDP per capita (GDPPC), urbanization rate (UR), unemployment rate (UE), and carbon dioxide emissions (CO2) are utilized as control variables (see Table 1). All variables are transformed into natural logarithms for analysis.

Table 1.Variable description and summary statistics
Indicators Description Mean St. Dev Minimum Maximum Source
Net Migration (NM) Net migration = Immigrants – Emigrants (including both citizens non-citizens) 15.717 0.070 13.815 15.966 WDI
Share of consumption from renewable energy (SRN) Proportion of total energy consumption, obtained from renewable source 2.559 1.166 -2.302 4.418 WDI
Renewable Electricity Per Capita (REPC) Total electricity from renewables (kWh per capita) 6.468 1.721 0.485 10.934 EI- Statistical Review (2024)
Energy Intensity (EI) Primary energy use per dollar of GDP (kWh/USD) 1.411 0.432 -0.162 2.783 U.S. EIA (2023)
GDPPC (constant 2021 international $) GDPPC based on PPP 10.143 0.803 7.906 11.840 WDI
UR (% of total population) People living in urban areas 4.174 0.319 2.901 4.587 WDI
UE (% of total labor force) Percentage share of total labor 1.866 0.653 -1.390 3.619 WDI
CO2 Annual emissions of carbon dioxide (CO2), excluding LULUCF. 4.754 1.514 1.048 9.492 WDI

Note: The table presents statistical overview, namely mean, standard deviation, minimum and maximum values for logged values of NM, SRN, REPC, EI, GDPPC, UR, UE and CO2, based on 1,464 observations for each variable; LULUCF: land use, land-use change, and forestry; EI Statistical Review: Energy Institute - Statistical Review of World Energy (2024); U.S. EIA: U.S. Energy Information Administrations

B. Model specification

Figure 1 briefly illustrates how the transition to renewable energy influences net migration through three key pathways, namely economic, environmental, and societal—by enhancing opportunities, reducing carbon emissions, and improving infrastructure and employment outcomes. It highlights the interconnectedness of these multidimensional factors in shaping migration trends.

Figure 1
Figure 1.Conceptual Framework

Notes: This figure depicts the conceptual framework of the study based on literature review.

The study utilizes the following model to assess the connection between net migration and renewable energy transition:

lnNMit=β0+β1lnSRN+β2lnREPC+β3EI+β4lnGDPPC+β5lnUR+β6lnUE+β7lnCO2+εit

where t represents time and εit represents the error term. For evaluating the above model, we employ the panel ARDL approach, which is not only effective in analyzing short-term fluctuations and long-term relationships among variables, but is also particularly useful when there is a mixed order of integration and potential cross-sectional dependence.

The PARDL long-run equation is specified as:

ΔNMit=α1+pk=1βijNMitj+qk=0δijθitj+εit

where i and t represent cross-sectional unit and timeframe, respectively. α1 is the intercept term, pk=1βijNMitj is the autoregressive (AR) term, which captures the lagged effects of NM, qk=0δijθitj is distributed lag of explanatory variables like GDPPC, etc., βij and δij are the coefficients, εit is the error term, j is maximum lag, θitj shows the lag of independent variables. p and q are the lag order.

The PARDL short-run equation is specified as:

ΔNMit=φi+pk=1XijΔNMitj+qk=0YijΔθitj+ijECTti+εit

where ECT signifies error correction term. ij is the convergence from long-run to short-run equilibrium.

Further, the study uses bias-adjusted Langrange multiplier (LM) test and the Pesaran cross-sectional dependence (CD) test, to check for the cross-sectional dependence among panels. To assess the stationarity of the variables, the cross-sectionally augmented Im-Pesaran-Shin (CIPS) unit root test is applied. We also conducted the Durbin-Wu-Hausman test, which shows that the model is free from endogeneity concerns. The model incorporated one lag per variable, selected using a general-to-specific approach based on parsimony and statistical significance, consistent with the panel ARDL framework.

III. Results

Table 2 presents the CSD & CIPS test statistics along with diagnostic tests. Cross-sectional dependence is confirmed for all variables. The CIPS test shows that lnSRN, lnREPC, lnGDPPC, lnUR, and lnUE are nonstationary at levels (i.e., I(1)), yet stationary in first differences (i.e., I(0)), confirming the presence of unit roots; while lnNM, lnEI, and lnCO2 exhibit stationarity at levels (i.e., I(0)). This mixed order of stationarity justified our model selection and also indicates the likely presence of cointegration. The poolability test result indicates heterogeneity across countries and justifies the use of a panel method that accounts for heterogeneity rather than pooled OLS. The Wooldridge test result suggests that our model does not suffer from significant serial correlation issues. The modified Wald test result validates the existence of heteroscedasticity across panels. The Hausman test confirms the fixed effects specification. Lastly, the Durbin-Wu-Hausman test result indicates no endogeneity concern.

Table 2.Findings from CSD test, CIPS unit root tests, and diagnostic checks
Variables Cross-sectional Dependence test CIPS
Pesaran Cross-section Dependency Bias-corrected scaled LM Level First difference
lnNM 15.553*** 66.374 -2.598** -4.241***
lnSRN 24.282*** 324.419 -2.124 -4.325***
lnREPC 111.190*** 299.424 -2.394 -4.967***
lnEI 140.822*** 455.348 -2.615** -4.834***
lnGDPPC 164.518*** 480.768 -2.048 -3.392***
lnUR 145.322*** 523.530 -2.053 -3.690***
lnUE 23.620*** 126.341 -2.129 -3.859***
lnCO2 17.585*** 336.899 -2.566* -4.338***
Diagnostic Test
Poolability Test F( 60, 1396) = 9.02***
Wooldridge Test F( 1, 60) = 2.872
Modified Wald Test chi2 (61) = 1511203.04***
Hausman Test chi2(6) = 29.50***
Durbin-Wu-Hausman Test (‘REPC’ variable) chi2(1) = 0.141
F(1,1333) = 0.159

Note: The table provides CSD and CIPS unit root test results. *, **, *** denote statistical significance at 10%, 5%, and 1% levels, respectively.

The Panel ARDL results given in Table 3 indicate that SRN and REPC are both positively and significantly associated with higher NM, likely because cleaner and more sustainable regions attract more migrants (Łukaniszyn-Domaszewska et al., 2025). Similarly, EI also has a significant positive influence on NM, suggesting that countries with high EI, reflecting greater industrialization and manufacturing activities as well as higher employment opportunities, are more likely to attract migrants. This is in alignment with Environmental Migration Theory and the New Economics of Labour Migration theory (Stark & Bloom, 1985). Expectedly, the UE is found to have a significant negative long-run effect, indicating that higher unemployment in a host country discourages immigration, which aligns with the push-pull theory of migration (Lee, 1966). GDPPC is found to be positively related to NM (Jennissen, 2003), albeit insignificant for the sample countries. The UR does not exhibit statistically significant effects in the long run, implying that other structural or environmental factors may play a more decisive role in influencing migration patterns.

Table 3.PARDL test results: (DV- Net Migration)
Variables PARDL Results Robustness Check (Sub-sample Analysis)
Long-Run Dynamics
lnSRN 0.001*
(0.000)
0.006
(0.006)
lnREPC 7.54e-07**
(3.35e-07)
0.019***
(0.005)
lnEI 6.057e-4**
(2.768e-4)
0.104***
(0.021)
lnGDPPC 0.001
(0.001)
0.006
(0.052)
lnUR -0.005
(0.007)
0.106
(0.119)
lnUE -0.002***
(5.024e-4)
-0.007
(0.009)
lnCO2 -0.002*
(0.001)
-0.022
(0.022)
Short-Run Dynamics
D(SRN) 0.015
(0.028)
0.094
(0.105)
D(REPC) -2.98e-5
(3.41e-5)
0.034
(.021)
D(EI) -0.002
(0.009)
-0.016
(.216)
D(GDPPC) -0.016
(0.094)
-0.294
(0.329)
D(UR) -3.322
(2.215)
-6.133
(8.143)
D(UE) 0.029
(0.031)
0.029
(0.051)
D (CO2) 0.080*
(0.047)
0.227
(0.153)
ECT -0.557**
(0.049)
-0.492***
(0.084)

Note: *, **, *** imply statistical significance at 10%, 5%, and 1% level, respectively. DV means dependent variable. Estimates are based on a PARDL (1,1,1,1,1,1,1,1) specification, with one lag for each variable, for simplicity, given T=24.

In the short run, except for CO2 emissions, all other explanatory variables are found to be insignificant in their interaction with migration. CO2 is found to be positively related to NM, possibly suggesting that short-run economic booms might temporarily attract migrants. Furthermore, the negative error correction term is found to be significant at the 5% level, indicating a stable long-run relationship among variables.

Additionally, we’ve conducted a robustness check, employing sub-sample analysis, in which we considered 20% of countries from each income group category (high income, lower- and upper-middle income), which further confirms the direction of our main findings.

To substantiate the long-run relationships established through the PARDL framework, we employed a suite of complementary panel cointegration tests, i.e., Pedroni (2004), Kao (1999), and Westerlund (2007). As presented in Table 5, the uniform refutation of the null hypothesis indicates strong evidence for a stable long-run equilibrium among the variables. These findings enhance the credibility and robustness of the primary PARDL estimates.

Table 4.Panel cointegration test results
Name of Testing Category Test Statistic
Pedroni Modified Philips-Perron t-statistics 6.959***
Philips-Peron t-test -10.104***
Augmented Dickey-Fuller (ADF) t-test -11.044***
Kao Modify Dickey-Fuller (DF) t-statistics -5.492***
Dicky-Fuller t-test -12.242***
Augmented Dickey-Fuller t-test 6.229***
Unadjusted modified Dickey-Fuller t-test -41.710***
Unadjusted Dickey-Fuller t-test -26.809***
Westerlund Variance ratio test -3.188***

Note: Here, *** implies 1% level of significance.

IV. Conclusion

The findings confirm that the share of renewable energy, renewable electricity per capita, and energy intensity positively impact net migration. Among the socio-economic variables, GDPPC is found to positively influence NM, while factors like CO2 emissions, the unemployment rate, and the urbanization rate negatively affect migration. To address regional disparities in migration, policymakers and stakeholders should leverage energy transitions and promote sustainable development. In developing nations, subsidies and credit schemes for solar and wind facilities can help create employment opportunities, and expanding vocational training in collaboration with renewable energy firms will equip workers with installation and maintenance skills, promoting inclusive green jobs. Developed economies should improve their energy infrastructure and invest in innovation clusters to attract skilled migrants and support long-term urban growth. Carbon pricing and emissions trading can curb CO2 emissions and climate-induced displacement. More importantly, region-specific renewable investments can bridge regional disparities and alleviate migratory pressures.