1. Introduction
The environment, health and reasonable labor mobility are currently hot social topics. Since 1949, the scale of labor mobility has been unprecedented, and it has become an important factor for promoting social and economic development [
1]. From 1982 to 2015, China’s labor mobility scale first increased and then decreased, with the total population migration across counties enlarging 4.1 times from 1982 to 2010, and decreasing 9.59% from 2010 to 2015 [
2], becoming the main labor force for China’s economic development, especially in developed regions such as Guangdong, Shanghai and Beijing. Indeed, the characteristics of labor mobility in China can be divided into three stages. Firstly, the early stage of labor mobility (1949–1978). In this stage, only a small number of migrant workers can be migrated to another place, especially to the cities, and labor mobility is mainly based on the government’s planned policies to serve social stability and economic development. Secondly, the middle stage of labor mobility (1979–2010). In this stage, laborers begin to flee their regulated workplace for a better income, and the labor mobility scale across regions is growing. Both China’s labor mobility index and the annual average inflow of the mobile population are continuously increasing before 2010, as shown in
Figure 1. Thirdly, the transition stage of labor mobility (after 2011). In this stage, the characteristics of health-type mobility appeared, and health demands come to be a more important factor influencing labor’s mobility.
In the third stage, the deaths caused by air pollution increased. According to the Lancet’s “Global Burden of Disease Report 2019,” China’s PM2.5 caused 1.2 million early deaths and 25 million disabled-adjusted life-year losses in 2019. As the threat of air pollution and health damage continues to increase, the average annual labor mobility index (net labor inflow rate) in each province of China has continued to decline, and this indicates increased labor fleeing as shown in
Figure 2(the graph is consistent with Zhang and Wang (2020) [
3]). Based on this phenomenon, scholars believe that air pollution has a significant negative impact on labor mobility [
4,
5,
6].
The negative relationship between air pollution and labor mobility refutes the view of surplus labor depletion under the “Lewis Turning Point” from the perspective of exogenous influence. Especially in the developed and eastern coastal regions in China, people decide their own labor supply on facts as health and social investment, the labor income, the leisure replacement and healthy life, and we can also see the counter-urbanization phenomenon known as “fleeing Beijing, Shanghai and Guangzhou” [
6]. However, there are few empirical studies on the relationship between air pollution, health and labor mobility. Most studies do not provide direct evidence to confirm air pollution’s influence on labor mobility through damaging health, and focus especially on the healthy migration effect and the salmon bias effect after the self-selection of Latin immigrants (also known as the Hispanic paradox) and are limited to using the pollution amount or concentration and other related variables to construct indicators for indirect analysis [
5,
7]. The health paradox of Latino immigrants (the Hispanic paradox) refers to the fact that Latino immigrants with lower socioeconomic status in the United States who have limited access to healthcare generally have a higher health level than non-Hispanic whites [
8,
9]. Indeed, facing the health shocks caused by air pollution, the probability of labor migration across provinces is higher than that across counties, and health shocks on labor mobility are more sensitive in larger geographical regions [
10]. Therefore, this article focuses on providing direct evidence of air pollution’s impact on labor mobility, anchoring the health shocks perspective in provincial samples. The potential contributions of this article are as follows: Firstly, using the principal component-entropy weight combination method, this article calculated the air pollution index and the health shocks effect index, and provided a direct provincial sample for empirical analysis. Second, using the mediation effect model, we found that labor will be impacted by health shocks when it suffers from air pollution. Labor mobility will be promoted to change migration decisions, and the inflow rate of labor mobility across provinces will decrease. These findings will provide additional evidence for the debate on the Lewis turning point and the exploration of the health paradox of Latin immigrants.
The remaining sections of this article are as follows:
Section 2 is the literature review, summarizing the methods and conclusions about the relationship between air pollution, health shocks and labor mobility, and providing related hypotheses.
Section 3 is the methodology, introducing the mediation effect method and data resources and description.
Section 4 reveals the result, showing the impact of air pollution on labor mobility, and the mediation effect of health shocks and the threshold effect.
Section 5 is the discussion, focusing on the study’s results and their limitations. Finally, in
Section 6, the conclusions and policy implications are provided.
4. Results
4.1. Characterization of Main Variables
Based on the measurement results of the above methods, from the annual average of the past three years (see
Figure 5), the air pollution index (
pol), health shock index (
d) and the health shocks composite index of air pollution (
h) vary greatly among provinces, and the index of labor mobility (
l) is generally low.
Furthermore, the annual average growth rate of the labor mobility index (
lz,
lzt = (
lt/
lt−1) − 1), the annual average value of air pollution index (
pol), the annual average value of the health shock index (
d) and the annual average value of the health shocks composite index of air pollution (
h) of each province are shown in
Figure 6. We can see that the
lz line shows a positive tendency, and other three lines are showing a negative trend, indicating that the growth in the health shocks of air pollution may be closely related to the decline of the net inflow rate of labor mobility. We will confirm this in the next section.
4.2. The Baseline Results
A fixed or random model is adopted according to the Hausman test. A dynamic system GMM model was adopted to estimate the first-order hysteresis effect and to reduce the endogenous problem, such as the potential autocorrelation effect of the net inflow rate of labor mobility, and the potential impact of the agglomeration effect formed by labor mobility on environmental pollution [
49], and
Table 2 shows the estimation results.
In column (1), the coefficient of the air pollution index is significantly negative at the 5% confidence level, indicating that air pollution will reduce the net inflow of labor mobility, and the negative impact is further confirmed in column (4). As shown in column (4), for each 1% increase in the air pollution index, the net inflow of labor mobility will drop by 47.6%, significant at a 1% confidence level. Using a sample of 285 cities in China, Zhang (2019) [
50] suggested that the number of employed workers will be reduced 33.45% with the SO
2 emissions increasing 1%, and our results further confirm the negative impact of air pollution on labor mobility. Using the comprehensive calculated air pollution, not a single pollutant, as the air pollution index, the negative impact in our article is larger than others. In column (2), the coefficient of the interaction term between air pollution and health shocks is significantly negative at the 5% confidence level. This implies that the net inflow of labor mobility will be reduced by the interaction of air pollution and health shocks, and the endogenous problem underestimates the negative impact shown in column (5). In column (5), when controlling other conditions, for each 1% increase in the interaction item of air pollution and health shocks, the net labor inflow of labor mobility will drop by 15.7%, significant at the 1% confidence level. In column (3), the coefficient of the health shock index of air pollution is significantly negative at the 5% confidence level, indicating that the net inflow of labor mobility will be decreased 24.9%, with a 1% increase in the health shock index of air pollution. In column (6), a dynamic system GMM model is adopted that confirms the negative impact of the health shock index of air pollution on the net inflow of labor mobility. Therefore, we can see that the health shocks of air pollution are a significant hindrance to the net inflow of labor mobility, and confirm the related hypothesis in H1.
4.3. The Mediation Effect Results
Based on Equations (3) to (5), the mechanism effect results are shown in
Table 3. The variable pol has a significant negative impact on the net inflow of labor mobility at the 5% confidence level in column (7). In column (8), the coefficient of the air pollution index is significantly positive at the 10% confidence level, indicating that there is a negative impact of the air pollution index on health shocks, and health shocks will be increased 16.5% with an air pollution index increase of 1%. The coefficient of health shocks is significantly negative at the 1% confidence level in column (9). However, the coefficient of the air pollution index is not significant in column (9), which means that it exerted a full mediation effect. Therefore, it can be concluded that the impact of air pollution on the labor inflow is all caused by health shocks, and that rectifying health shocks is also an important factor influencing the decrease in the labor supply and migrating individuals with different health levels under different air pollution conditions.
In order to further ensure the reliability of mediating effect results, the Sobel test is used in this paper. In accordance with Sobel statistics [
51], the Sobel values and the
p values were directly calculated using the coefficients of each major effect, and the Sobel 95% confidence interval was [1.616, 1.738], with a significant value of −1.67 < −0.97 at the 5% confidence level, confirming the mediating effect. As mentioned above, the health shocks of air pollution restrict the net inflow of labor mobility, and the hypothesis H1 has been verified.
4.4. The Robustness Test
Based on the current study, we adopt two main ways to test the robustness of the results above, and
Table 4 shows the result. Firstly, we use one lagging term based on the main explanatory variables, and the results are displayed in column (10) to column (13). We can see that the coefficients of the variables pol*d and h are significantly negative at the 10% (or less) confidence level, and confirm the stability of the results above. Secondly, according to the idea of hedonic price, the housing price can measure the labor inflow in one region [
52]. However, the housing price does not reflect the real market demand in China influenced by the housing market bubble. Therefore, Xi and Liang (2015) [
53] used the idea of choosing a housing sales area as a proxy variable to analyze the effect of environmental migration. Meanwhile, to weaken the scale effect of the provincial population and the influence of natural growth, again referring to Zhang and Wang (2020) [
3], we use the natural logarithm of the commodity housing sales area per (lniv6) as the proxy variable of the labor mobility index, and the results are shown in column (14) to column (17). We can find that most coefficients of the variables h, pol and d are significantly negative at the 10% (or less) confidence level, indicating that the conclusion above is robust.
In addition, we use the lagged term of the air pollution index and the health shock index to test the medication effect, and
Table 5 shows the results. In
Table 5, we can see that the impact of air pollution on the net inflow of labor mobility is mainly caused by the mediating effect of health shocks. Thus, the reliability of the mediating effect above is rectified (MacKinnon et al. (1995) [
54] set 0.97 as the boundary value at the 5% significance level for testing distribution).
4.5. Results of Regional Conditions
Furthermore, from the regional perspective, we analyze the impact of the health shocks of air pollution (h) on the index of labor mobility (l). Since the Seventh Five-Year Plan, China has been strategically divided into four economic development zones in the east, center, west and northeast. Now, these have been redivided into five major economic belts of “two horizontal and three vertical.” In order to find the regional differentiates of the samples, we choose the eastern, central and western economic belts and the Yangtze River economic belt as the heterogenous conditions. Scholars believe that the difference in air pollution on the Huai River between the north and south is large, and the health-damaging effects are also accordingly different [
6,
11]. Therefore, it is necessary to analyze the regional differentiates according to the boundary of the Huai River. The results are shown in
Table 6.
The health shocks of air pollution in the east have significantly hindered the index of labor mobility. When other conditions remain unchanged, for every unit of the health shocks composite index of air pollution, the net inflow rate of labor mobility drops significantly by 33.1%, while the impact in the central and western regions is not significant. This result indicates that the labor force in the eastern region has fled, facing severe health shocks from air pollution. According to the sample, compared with 2007, in 2015, the net inflow rate of labor mobility in the eastern region declined by 86.35%, while that in the central and western regions increased by 85.24% and 196.44%, correspondently.
The impact of the health shocks of air pollution on labor mobility in the Yangtze River economic belt is significantly negative. When other conditions remain unchanged, for each unit increase in the health shocks of air pollution, the net inflow rate of labor mobility drops significantly by 34.6%. This is in line with the fact that there were 12,158 chemical enterprises above the designated size in the Yangtze River economic belt in 2016, accounting for 46% of those in China, with 41% heavy industry among them, forming a “chemical industry surrounding the river” with high pollution emissions. This result implies that in the process of promoting high-quality development, the negative impact of health shocks from air pollution on labor mobility cannot be ignored in the Yangtze River economic belt. Indeed, the rectification of chemical enterprises in the Yangtze River economic belt also took place after 2015, and this effectively avoided the exogenous effects of labor “flight.”
The impact of the health shocks of air pollution on labor mobility in the northern part of the Huai River is significantly negative, while not being significant in the southern part. When other conditions remain unchanged, for each unit increase in the health shocks of air pollution, the net inflow rate of labor mobility in the northern Huai River drops significantly by 31.8 %. Because the heating supply has increased more than air pollution only in the northern part of the Huai River, the impact of health shocks of air pollution on labor mobility in north is accordingly obvious [
6]. Therefore, the hypothesis H2 has been verified.
4.6. The Threshold Effect Results
According to the current study, the health shocks caused by air pollution changed greatly after 2010 [
39]. We chose 2010 as the time threshold variable, to confirm the impact of the health shock of air pollution on the Lewis inflection point theory. We also chose the health shocks of air pollution as a threshold variable to explore the level of the health shocks of air pollution’s influence on labor mobility, and the results are shown in
Table 6. In
Table 6, whether time or the health shocks of air pollution are the threshold variable, the null hypothesis cannot be accepted without a threshold in the single threshold situation, and the alternative hypothesis test with two or three thresholds fails, indicating that there is one threshold.
Table 7 and
Table 8 show the estimated results with a single threshold variable, and
Table 9 shows the impact effect estimation of threshold points.
In
Table 8, the time threshold value is 2010 at the 5% confidence level. In column (27), the coefficient of the time threshold before 2010 is not significant, and the coefficient of the time threshold after 2010 is significantly negative at the 1% confidence level, indicating that the negative effect of the health shocks of air pollution on labor mobility is not significantly obvious before 2010, and the negative effect is obvious after 2010. This result shows that the health shocks of air pollution in China have intensified since 2010.
In
Table 8, the threshold value of the health shocks of air pollution is 1.9873 at the 5% confidence level. In column (28), the coefficient of the variable h is −0.523 and −0.389 at the 1% confidence level, respectively. This indicates that the health shocks composite index of air pollution has one threshold effect on the net inflow of labor mobility, and the impact will decrease when the health shocks of air pollution cross the threshold value. In order to explore this interesting phenomenon,
Figure 7 shows the sample distribution of real income level under the two risk types of the health shocks of air pollution (there are two reasons to choose the variable of real income. First, since the reform and opening up, one of the main purposes of labor mobility between urban and rural areas in China is to increase family income, and improve family economic conditions. Second, the persistence of severe air pollution in China has been accompanied with extensive economic growth over a long period, and this growth can be reflected in real income factors. Therefore, we chose the real income as the representative factor of the gravity effect for analysis.). The average per capita real income level in type one and two regions is 6.656 and 6.339, respectively, and the percentage of the value is more than 7.106 (the threshold value of 7.106 is obtained by the clustering per capita real income (gp) according to the K-means). The proportion obtained by contingency table analysis is 31.2% and 31.1%, and the standard deviation value is 1.860 and 2.106, respectively, among them. We can find that the average per capita real income, and its percentage, with a value of more than 7.106 in type two is larger than in type one, while the standard deviation value corresponding to type two is less than type one, indicating that the stable and higher real income plays a strong attracting role on the net inflow of labor mobility when the health shocks of air pollution crosses the threshold value, and the attracting effect partially offsets the repulsion effect of the health shocks of air pollution on the inflow of labor mobility.
In addition,
Figure 8 shows the LR statistics value trend estimated for the threshold variable. We can find that the LR values in 4(a) and 4(b) are less than the critical value 7.35 at the 5% confidence level (the dashed line), further confirming the validity of the threshold effect (due to space limitations, other threshold LR test results are omitted). Therefore, hypothesis H3 has been verified, and labor mobility in China is currently simultaneously affected by income and the health shocks of air pollution.
5. Discussion
This article used the principal component entropy weighted method to calculate the air pollution composite index and health shock index. Most studies have chosen the AQI or the PM2.5, alternative variables of air pollution, to analyze the impact of air pollution on labor mobility [
5,
6], but one pollutant, such as in the PM2.5, cannot fully represent the total air pollution, and the calculation process of the AQI index does not eliminate the correlation between different pollutants. The principal component entropy weighted method can identify the most important information of air pollution and health shocks, and avoid the defect of the AQI index and the single air pollutant, to make the conclusion more reliable.
Based on the air pollution composite index and the health shock index above, this article analyzes the impact of air pollution on labor mobility, and the results show that labor mobility is affected by the health shocks of air pollution. Similar results were found by Cropper (1981) [
55] and Sun et al. (2019) [
39]. However, most studies focused on the impact factors such as the housing price or environmental regulation [
56,
57] and ignored the impact of health shocks. Adopting the mediation effect model, our results further imply that the impact of air pollution on labor mobility is all caused by health shocks, and health shock is a very important factor of air pollution’s impact on labor mobility at the macro level. We can explain these results from two angles: firstly, although individuals have the motivation to invest in their health and human capital, they will display mobility behavior to avoid air pollution, and show more as a moderating effect to avoid health shock [
58]. Secondly, when one region faces a serious health shock caused by environmental pollution, it will firstly exclude migrant labor in a poor health condition, and eventually extract the migrant population in a good health condition, resulting in the health level of the migrated population being higher than that of the locals, as described in the Latino health paradox. This phenomenon is similar to the health paradox of Latino immigrants with a health self-screening mechanism [
59].
Using the threshold model, our results further show that the health shocks of air pollution on labor mobility across provinces became more significant after 2010, and we confirm the view that China’s Lewis turning point was 2010 [
60]). However, compared with most current views [
61], our results more strongly emphasize health shocks in China’s Lewis turning process and its avoidance effect, and argue that the health shocks of air pollution are also an important reason for the decline of labor mobility supply across provinces. It has an obvious difference with the impact of different wages caused by the urban–rural dual structure. Indeed, the impact of the health shocks of air pollution on labor mobility will increase the labor’s burden of medical expenses, and decrease the labor’s expected income and life expectancy [
62]. In addition, when the labor anticipates that there may be a large health shock in one region with severe environmental pollution, the labor will reduce its mobility intention, and the regional labor inflow rate will decline.
Finally, based on the push–pull theory, we conclude that the attraction effect of regional stability and high real income will partially offset the repulsion effect of the health shocks of air pollution on labor mobility, when the health shock index of air pollution exceeds the threshold value of 1.9873, similar to the view of Li et al. (2020) [
63]. However, most studies focus only on the negative impact of air pollution on labor mobility. Environmental pollution has null jointness. It is not only a byproduct of production, but also an important input of production [
64], especially in developing countries or regions, and the regional economy can quickly achieve extensive growth with a corresponding improvement in pollution levels.
Compared with the empirical conclusions of previous studies, we found similar conclusions, with some differences. We confirm the negative impact of air pollution on labor mobility that Chen et al. found (2017) [
5,
6]. However, we find the impact of air pollution on labor mobility is all caused by health shocks, and this result is obviously different from most studies, such as that of Kahn and Mansur (2013) [
56]. In addition, we find the threshold effect of time and health shocks of air pollution, and confirm China’s Lewis turning point was 2010, in a manner similar to Kwan et al. (2018) [
60]. However, our study has several limitations. Firstly, the study period ranged from 2008 to 2015. For policy effect interference and data limitation, we did not extend the data to 2020 or 2021, and the statistical bias may underestimate the negative impact of air pollution on labor mobility. Second, the air pollution index and health shock index were calculated by the principal component entropy weight method, and do not consider the weight valued by related experts. This might have resulted in the underestimation of the mediation effect of health shocks. Thirdly, this article did not use the geographical method to display the spatial labor mobility and spatiotemporal impact of air pollution on labor mobility. Therefore, we may have to empirically analyze these in a further study.