The empirical analysis of time series data conventionally begins with an examination of the order of integration of the series. Towards this end we apply the ADF and KPSS tests. As shown in the Table 1, all the series are integrated of order one as per the ADF test. The results of KPSS test reveal that while unemployment rate for men (UNEMP_M) and trade openness (TO) are stationary at levels, the rest of the series are I(1).
After ensuring that none of the variables is I(2), we conduct the cointegration analysis in the ARDL Bounds framework.
Table 1
Results of ADF and KPSS Tests
Variable
|
ADF Test
|
KPSS Test
|
|
Lags
|
t-statistics
|
Bandwidth
|
LM statistics
|
GCO2
|
0
|
3.751
|
4
|
0.189
|
∆GCO2
|
0
|
-5.655***
|
2
|
0.079
|
UNEMP_M
|
1
|
-2.863
|
3
|
0.073
|
∆UNEMP_M
|
0
|
-4.378***
|
1
|
0.072
|
UNEMP_F
|
2
|
-0.819
|
4
|
0.528
|
∆UNEMP_F
|
1
|
-8.173***
|
8
|
0.183
|
GDP_GR
|
0
|
-1.668
|
0
|
0.159
|
∆GDP_GR
|
0
|
-4.747***
|
0
|
0.085
|
TO
|
0
|
-0.380
|
4
|
0.141
|
∆TO
|
0
|
-4.892***
|
0
|
0.134
|
Note : Δ represents first difference. *** represents significance at 1%.For ADF test, the null hypothesis is of unit root. For KPSS test, the null hypothesis is of stationarity of series. |
Since we have annual time series data, we consider two as the maximum order of lags in the model (Narayan ,2004). The Akaike information criterion is used to determine the optimal number of lags. The results contained in Table 2 reveal that the F-statistic exceeds the upper bound critical value at all the conventional levels of significance. Therefore, the null hypothesis of no cointegration is rejected and we conclude that the variables are cointegrated.
Table 2
Results of ARDL bound test for Cointegration
|
Critical Values
|
F-statistic = 16.022
|
α = 0.01
|
α = 0.05
|
α = 0.10
|
I(0)
|
I(1)
|
I(0)
|
I(1)
|
I(0)
|
I(1)
|
4.590
|
6.670
|
3.354
|
4.774
|
2.752
|
3.994
|
Note : The lag order of the model is based on AIC |
Tables 3 and 4 report the estimated long run and short run coefficients of the model. As evident from Table 3, there is a significant negative impact of male unemployment rate (UNEMP_M) on the growth in CO2 emission in the long run. In particular, a one percentage point increase in male unemployment rate results in 0.0073 percentage point reduction in CO2 emission growth. The results thus confirm the presence of trade-off between male employment and pollution control. In other words, a decrease in CO2 emission growth may be achieved at the expense of employment of male workers.
Table 3
Estimated long-run coefficients
Variable
|
Coefficient
|
t-statistics
|
p value
|
UNEMP_M
|
-0.0073
|
-2.964
|
0.009
|
UNEMP_F
|
0.0114
|
2.964
|
0.009
|
GDP_GR
|
0.199
|
6.088
|
0.000
|
TO
|
0.060
|
1.228
|
0.230
|
R2
|
0.88
|
Breusch-Godfrey LM Test for Serial Correlation
|
F = 0.413 (p val = 0.669)
|
Adj R2
|
0.78
|
Breusch-Pagan-Godfrey Test for Heteroskedasticity
|
F = 0.247 (p val = 0.987)
|
|
|
Ramsey Reset Test
|
F = 1.424 (p val = 0.251)
|
Table 4
Estimated short-run coefficients
Variable
|
Coefficient
|
t-statistics
|
p value
|
UNEMP_M
|
-0.014
|
-2.092
|
0.016
|
UNEMP_F
|
0.005
|
3.2014
|
0.005
|
GDP_GR
|
0.183
|
2.083
|
0.050
|
TO
|
0.292
|
1.516
|
0.153
|
ECT
|
-1.568
|
-10.548
|
0.000
|
Interestingly, the coefficient associated with female unemployment rate is significant and positive in the long run. A one percentage point increase in unemployment of women in India leads to increase in the environmental pollution by 0.011 percentage points. In other words, a reduction in female unemployment results in lower environmental degradation and there is no trade-off relationship between female employment and environmental quality in the long run. The results also indicate that in the long run, while the GDP growth tends to increase the CO2 emission growth, the trade openness does not exert any impact on environmental degradation.
All the variables retain their sign and significance in the short run as well. As shown in Table 4, the coefficient associated with male unemployment rate is negative indicating a short run trade-off between male employment and environmental quality. In short run too, the regression coefficient on female unemployment rate is positive suggesting that a reduction in female unemployment is associated with lower CO2 emission in the short run.
The negative and significant error correction term in Table 4 confirms the cointegration amongst the variables. The ECT suggests that the system corrects its previous period’s disequilibrium at a speed of 156.8% annually. This speed may appear overadjusted as ECT conventionally lies between one and zero. According to Narayan and Smyth (2006), an ECT lying between one and two suggests that instead of a monotonic convergence to the equilibrium path, the error correction process fluctuates around the long-run value in a dampening manner. Once the path is complete, there is a rapid convergence to the equilibrium path (Singh and Shastri,2019).
Various diagnostic tests were applied to assess the appropriateness of the specification of the model. The Breusch Godfrey LM test for autocorrelation and Breusch-Pagan-Godfrey Test for Heteroskedasticity fail to reject the null hypothesis of no autocorrelation and homoskedasticity respectively. Further, the Ramsay reset test suggests that the functional form of the model is appropriate. Finally, Fig. 1 presents plots of CUSUM and CUSUMSQ test statistics that fall inside the critical bounds of 5% significance suggesting that the parameters in the model are stable.
In the interest of checking the robustness of the above evidence, we also apply Fully Modified Least Square (FMOLS) and Dynamic Ordinary Least Squares (DOLS) estimators. The results of FMOLS and DOLS provide the same evidence for a negative, significant impact of male unemployment rate on CO2 emission growth and a positive impact of female unemployment rate on CO2 emission growth.
Table 5
Long Run Coefficients under Alternative Estimators
Variable
|
FMOLS
|
DOLS
|
|
Coefficient
|
t-statistics
|
Coefficient
|
t-statistics
|
UNEMP_M
|
-0.0060
|
-2.208
(0.037)
|
-0.0059
(0.025)
|
-2.309
|
UNEMP_F
|
0.0115
|
2.244
(0.034)
|
0.016
|
2.458
(0.022)
|
GDP_GR
|
0.179
|
2.908
(0.050)
|
0.193
|
5.666
(0.000)
|
TO
|
0.061
|
1.545
(0.134)
|
0.075
|
1.695
(0.102)
|
R2
|
0.855
|
|
0.871
|
|
Adj R2
|
0.753
|
|
0.742
|
|
Note : Figures in parenthesis are p values |
All three methods, namely , ARDL,FMOLS and DOLS provide approximately the same outcomes, which validate existence of Environmental Philips Curve for male unemployment rate while the Environmental Philips Curve does not exist for women unemployment.
A few reasons explain this discrepancy and the absence of trade-off between female employment and environmental quality. First is the occupational distribution/segregation between men and women in India. The NSSO data (1970–2018) exhibits that women have largely been undertaking labor-intensive and informal work, concentrated in low-productivity sectors. While nearly 77% of rural women are employed in agriculture sector, for the urban women, service sector is the key employment provider with its share in employment rising from 35% in 1977-78 to 60% in 2017-18. In this sector, women are largely concentrated in professions such as teaching and nursing considered to be the cleaner sectors of the economy. Hence, women’s employment does not directly result in environmental degradation.
Not only that we find the absence of any trade-off between environmental quality and employment, we find a positive effect of women’s employment on the pollution control. As the female unemployment decreases, environmental degradation also reduces. The favourable effect of female employment on the environment may emerge from two channels. First, with the reduction in unemployment or access to work, women tend to go for family planning and birth control and wish to have children later in life. This may slow down population growth and relieve stress on the environment (Jameel et al., 2022).
Second, the employment of women increases women’s autonomy towards economic resources and results in women empowerment (Guinée, 2014; Peinado and Serrano 2018) giving them a greater role in decision making at various levels. Researches show that in capacity of decision makers, women are more sensitive towards environmental issues. For instance, illustrating the case of agriculture sector, O’Connor (2019) observes that giving women leadership roles leads to increased innovations to make farming more sustainable. Women often possess valuable knowledge and potential for sustainable innovation particularly with regard to agriculture technology. Desrochers et al.( 2019) opine that since women are more conscientious (i.e. goal directed, organized) than men, they are more likely to support environmental protection and indulge in pro environment decisions. Hence, an increase in women’s economic activity is conducive to more environment friendly decision making at various levels.
We next employ the block exogeneity Wald test (based on VECM) to examine causal relationship amongst the series. From the results presented in Table 5 it is evident that there is a unidirectional causality from male unemployment rate and GDP growth rate to environmental degradation. Interestingly, there is a bidirectional causality between female unemployment rate and environmental degradation. This implies that in the Indian context, it is not only the female employment that may reduce the environmental degradation, an improvement in environment also facilitates reduction in female unemployment[1]. This effect may come from two channels. First, an improvement in environmental quality is associated with lesser care work (as children and elders have lower probability of falling sick). This results in a greater availability of time to be invested in educational and vocational skills which in turn favourably affect women’s employment prospects. Also, environmental degradation often leads to curtail women’s employment opportunity owing to the reduction in natural assets caused by air pollution. (Unsworth and Ormrod, 1982).An improvement in environmental quality may in turn increase availability of work in agriculture sector where women are overrepresented.
Table 6
Results of Cuasality Test
Excluded
|
Chi Square
|
p value
|
Dependent Variable : GCO2
|
|
|
UNEMP_M
|
3.750
|
0.050
|
UNEMP_F
|
6.331
|
0.011
|
GDP_GR
|
7.437
|
0.006
|
TO
|
0.011
|
0.913
|
Dependent Variable : UNEMP_M
|
|
|
GCO2
|
1.115
|
0.290
|
UNEMP_F
|
0.0484
|
0.486
|
GDP_GR
|
14.569
|
0.000
|
TO
|
15.006
|
0.000
|
Dependent Variable : UNEMP_F
|
|
|
GCO2
|
8.511
|
0.014
|
UNEMP_M
|
7.001
|
0.008
|
GDP_GR
|
6.370
|
0.041
|
TO
|
6.751
|
0.009
|
Dependent Variable : GDP_GR
|
|
|
GCO2
|
0.040
|
0.841
|
UNEMP_M
|
13.060
|
0.000
|
UNEMP_F
|
9.020
|
0.011
|
TO
|
3.628
|
0.053
|
Dependent Variable : TO
|
|
|
GCO2
|
0.459
|
0.497
|
UNEMP_M
|
18.877
|
0.000
|
UNEMP_F
|
21.246
|
0.000
|
GDP_GR
|
13.476
|
0.000
|
Furthermore, we find that all three variables, namely, trade openness, unemployment rate for males and females cause CO2 emission growth via the GDP growth channel. Also, we find a bidirectional causality between trade openness and male unemployment rate as well as between trade openness and female unemployment rate. Bidirectional causality also exists between GDP growth and both male and female unemployment rates. Additionally, there is a unidirectional causality running from male unemployment rate to female unemployment rate.