| Literature DB >> 35694370 |
Shaobin Zhang1, Hao Ji1, Maoxi Tian1, Binyao Wang1.
Abstract
In July 2021, China began its national emissions trading scheme, marking a new stage of development for the country's carbon market. This study analyzes the multidimensional correlation between carbon prices in the Guangdong pilot market and eight influencing factors from three perspectives (the international carbon market, energy prices, and China's economic situation), using the ARMA-GARCH-vine copula model. The CoVaR between the carbon price and each factor is then calculated using copula-CoVaR. The results show that the crude oil market plays the primary role in the vine structure, and that the carbon market is not strongly correlated with other markets. China's carbon market is still a regional market driven by government policy, and the international carbon and energy markets (especially the crude oil market) have upward risk spillover effects upon it. This indicates an asymmetric risk spillover between influencing factors and the carbon market. The findings of this study will help market participants prepare risk management strategies and make related investment decisions, and provide a reference for policy makers to formulate national emission trading scheme policies.Entities:
Keywords: China’s carbon market; High-dimensional nonlinear dependence; Major influence factors; Risk spillovers
Year: 2022 PMID: 35694370 PMCID: PMC9167038 DOI: 10.1007/s10479-022-04770-9
Source DB: PubMed Journal: Ann Oper Res ISSN: 0254-5330 Impact factor: 4.820
Fig. 1Five-dimensional vine copula structures
Fig. 2Monthly price and monthly trade volume of Guangdong pilot carbon market
Descriptive statistics of the assets returns
| Cpr | WTI | BRE | CER | EUA | Gas | Coal | CSI | Ind | |
|---|---|---|---|---|---|---|---|---|---|
| Mean | − 0.30 | − 0.13 | − 0.22 | − 0.04 | 0.54 | − 0.03 | 0.01 | 0.23 | 0.16 |
| Median | 0.00 | 0.37 | 0.12 | − 0.45 | 0.66 | − 0.24 | − 0.07 | 0.28 | 0.42 |
| Max | 64.26 | 139.01 | 19.49 | 279.63 | 25.46 | 24.57 | 5.22 | 10.66 | 8.55 |
| Min | − 54.67 | − 163.08 | − 36.27 | − 52.33 | − 32.12 | − 30.82 | − 2.24 | − 14.02 | − 14.77 |
| S. D | 11.78 | 12.73 | 4.92 | 17.35 | 5.82 | 4.82 | 0.85 | 3.22 | 3.33 |
| Skewness | 0.25 | − 2.38 | − 1.51 | 12.48 | − 0.23 | 0.32 | 2.31 | − 0.71 | − 1.08 |
| Kurtosis | 5.36 | 125.27 | 14.33 | 205.21 | 7.55 | 16.33 | 13.53 | 5.58 | 6.36 |
| J.B.(p) | 0.00 | 0.00 | 0.00 | 0.00 | 0.00 | 0.00 | 0.00 | 0.00 | 0.00 |
| LB-Q (12) | 0.00 | 0.01 | 0.00 | 0.03 | 0.00 | 0.00 | 0.00 | 0.78 | 0.65 |
| ARCH-LM (12) | 0.00 | 0.00 | 0.00 | 1 | 0.73 | 0.00 | 0.00 | 0.00 | 0.00 |
| ADF(p) | 0.01 | 0.01 | 0.01 | 0.01 | 0.01 | 0.01 | 0.01 | 0.01 | 0.01 |
Cpr represents the carbon price of Guangdong ETS, WTI represents WTI crude oil futures prices, Bre. represents the brent oil futures prices, CER represents CER futures price, EUA represents EUA futures price, Gas represents Chinese natural gas prices, Coal represents Chinese coal price index, CSI represents CSI300, Ind represents Industry Index of China. J.B. is the p-value of Jarque–Bera statistics, LB-Q is the p-value of Ljung-Box test with lagged order 12, ARCH-LM is the p-value of ARCH test with lagged order 12, ADF is the p-value of the augmented Dicky-Fuller unit root test
Fig. 3Vine structures of Guangdong pilot carbon market. Note: each node denotes the marginal distribution of the corresponding variable, each edge represents a bivariate copula for two connected nodes, and the label is the family of corresponding bivariate copula. 1 = carbon price of Guangdong ETS, 2 = WTI, 3 = BRENT, 4 = CER, 5 = EUA, 6 = Natural Gas, 7 = Coal, 8 = CSI300, 9 = Industry Index
The estimation result of vine copula model of Guangdong pilot carbon market
| Tree | Edge | Copula | Par.1 | Par.2 | τ | ||
|---|---|---|---|---|---|---|---|
| 1 | 5,6 | t-copula | − 0.11 (0.06) | 5.28 (1.58) | − 0.07 | 0.03 | 0.03 |
| 3,2 | t-copula | 0.88 (0.02) | 7.19 (2.99) | 0.68 | 0.42 | 0.42 | |
| 9,8 | t-copula | 0.91 (0.01) | 6.51 (2.33) | 0.73 | 0.57 | 0.57 | |
| 5,4 | Gaussian copula | 0.27 (0.06) | – | 0.18 | – | – | |
| 3,5 | Gumbel copula | 1.24 (0.06) | – | 0.20 | 0.25 | – | |
| 1,7 | rot270 Joe copula | − 1.16 (0.06) | – | − 0.08 | – | – | |
| 9,1 | survival Joe copula | 1.08 (0.04) | – | 0.08 | – | 0.10 | |
| 9,3 | Clayton copula | 0.31 (0.10) | – | 0.13 | – | 0.11 | |
| 2 | 3,8;9 | Frank copula | − 0.12 (0.33) | – | − 0.01 | – | – |
| 5,2;3 | Frank copula | − 0.31 (0.34) | – | − 0.03 | – | – | |
| 4,6;5 | survival Gumbel copula | 1.09 (0.05) | – | 0.08 | – | 0.11 | |
| 3,4;5 | rot270 Clayton copula | − 0.06 (0.08) | – | − 0.03 | – | – | |
| 9,5;3 | Gumbel copula | 1.06 (0.04) | – | 0.05 | 0.07 | – | |
| 9,7;1 | rot270 Joe copula | − 1.06 (0.06) | – | − 0.03 | – | – | |
| 3,1;9 | Gaussian copula | 0.05 (0.06) | – | 0.03 | – | – | |
| 3 | 5,8;3,9 | Frank copula | − 0.14 (0.37) | – | − 0.01 | – | – |
| 4,2;5,3 | Joe copula | 1.07 (0.06) | – | 0.04 | 0.08 | – | |
| 3,6;4,5 | Gaussian copula | − 0.01 (0.06) | – | − 0.01 | – | – | |
| 9,4;3,5 | Clayton copula | 0.05 (0.08) | – | 0.03 | – | 0.00 | |
| 1,5;9,3 | survival Clayton copula | 0.10 (0.09) | – | − 0.05 | 0.00 | – | |
| 3,7;9,1 | t-copula | 0.01 (0.06) | 7.83 (4.09) | 0.00 | 0.02 | 0.02 | |
| 4 | 4,8;5,3,9 | rot270 Joe copula | − 1.02 (0.05) | – | − 0.01 | 0.00 | – |
| 6,2;4,5,3 | survival Clayton copula | 0.12 (0.09) | – | 0.06 | 0.00 | – | |
| 9,6;3,4,5 | Frank copula | 0.05 (0.35) | – | 0.01 | – | – | |
| 1,4;9,3,5 | Joe copula | 1.05 (0.05) | – | 0.03 | 0.06 | – | |
| 7,5;1,9,3 | rot90 Joe copula | − 1.03 (0.06) | – | − 0.02 | – | – | |
| 5 | 6,8;4,5,3,9 | Frank copula | 0.56 (0.34) | – | 0.06 | – | – |
| 9,2;6,4,5,3 | survival Clayton copula | 0.01 (0.08) | – | 0.00 | 0.00 | – | |
| 1,6;9,3,4,5 | Gaussian copula | 1.02 (0.04) | – | 0.02 | 0.03 | – | |
| 7,4;1,9,3,5 | rot90 Joe copula | − 1.05 (0.06) | – | − 0.03 | – | – | |
| 6 | 2,8;6,4,5,3,9 | rot270 Joe copula | − 1.01 (0.05) | – | − 0.01 | – | – |
| 1,2;9,6,4,5,3 | Clayton copula | 0.03 (0.08) | – | 0.01 | – | 0.00 | |
| 7,6;1,9,3,4,5 | Frank copula | 0.45 (0.35) | – | 0.05 | – | – | |
| 7 | 1,8;2,6,4,5,3,9 | Gaussian copula | 0.01 (0.06) | – | 0.00 | – | – |
| 7,2;1,9,6,4,5,3 | Clayton copula | 0.07 (0.08) | – | 0.03 | – | 0.00 | |
| 8 | 7,8;1,2,6,4,5,3,9 | rot90 Clayton copula | − 0.02 (0.08) | – | − 0.01 | – | – |
1 = cprice, 2 = WTI,3 = BRENT, 4 = CER, 5 = EUA, 6 = Natural Gas, 7 = Coal, 8 = CSI300, 9 = Industry Index. The table reports parameter estimates for vine copula models and their standard errors (in brackets). For t-copula, Par.1 denotes correlation and Par.2 represents degree of freedom , respectively. And τ, and denote the Kendall coefficients, upper tail dependence coefficient and lower tail dependence coefficient, respectively
The upside VaR (95 percent quantile) and the downside VaR (5 percent quantile) of returns
| Cpr | WTI | BRE | CER | EUA | Gas | Coal | CSI | Ind | |
|---|---|---|---|---|---|---|---|---|---|
| VaRup | 15.26 | − 5.88 | 6.23 | 8.98 | 8.94 | 4.34 | 0.81 | 5.13 | 4.93 |
| VaRdown | − 16.68 | − 7.98 | − 7.41 | − 9.97 | − 8.15 | − 4.24 | − 0.87 | − 4.88 | − 4.95 |
The bivariate copula and corresponding parameters
| type | AIC | par1 | par2 | τ | |||
|---|---|---|---|---|---|---|---|
| cprice-WTI | Gumbel copula | 0.65 | 1.05*** (0.04) | – | 0.04 | 0.06 | – |
| cprice–BRENT | Gumbel copula | 0.18 | 1.05*** (0.04) | – | 0.05 | 0.07 | – |
| cprice-CER | Gumbel copula | − 2.03 | 1.04 *** (0.04) | – | 0.04 | 0.05 | – |
| cprice-EUA | Gumbel copula | − 1.95 | 1.08*** (0.04) | – | 0.07 | 0.10 | – |
| cprice-Gas | Gumbel copula | 1.13 | 1.03*** (0.04) | – | 0.03 | 0.04 | – |
| cprice-Coal | t-copula | 0.51 | − 0.13** (0.06) | 10.00** (4.31) | − 0.08 | 0.00 | 0.00 |
| cprice-CSI300 | Clayton copula | − 1.09 | 0.15* (0.08) | – | 0.07 | – | 0.01 |
| cprice-Industry Index | Clayton copula | − 2.50 | 0.17* (0.09) | – | 0.08 | – | 0.02 |
cprice represents carbon price of Guangdong ETS, ⁎⁎⁎ Statistical significance at the 1% level, ⁎⁎ Statistical significance at the 5% level, ⁎ Statistical significance at the 10% level
Fig. 4Dynamic VaR of carbon price and CoVaR conditional on other markets. Note: cprice represents the VaR of carbon price in Guangdong ETS, and subfigures A, B, C, D, E, F, G, H show the VaR of carbon price conditional on the fact that WTI, Brent, CER, EUA, Natural Gas, Coal, CSI300, Industry Index markets experienced an extreme movement
Risk spillover from other markets to the Guangdong pilot carbon market
| Cprice-WTI | Cprice-BRENT | Cprice-CER | Cprice-EUA | Cprice-Gas | Cprice-Coal | Cprice-CSI300 | Cprice-Industry Index | |
|---|---|---|---|---|---|---|---|---|
| CoVaR | − 4.560 | − 5.029 | − 5.465 | − 4.948 | − 2.189 | − 0.579 | − 3.254 | − 3.551 |
| ∆ CoVaR_D | 0.034 | 0.043 | 0.024 | 0.052 | 0.005 | 0.000 | 0.122 | 0.225 |
| %∆CoVaR_D | 0.75 | 0.85 | 0.45 | 1.07 | 0.25 | 0 | 3.88 | 6.76 |
| Rank | 5 | 4 | 6 | 3 | 7 | 8 | 2 | 1 |
| CoVaR_U | 12.705 | 11.649 | 17.379 | 18.987 | 6.592 | 0.663 | 6.514 | 6.397 |
| ∆CoVaR_U | 4.309 | 3.228 | 6.520 | 7.238 | 2.000 | − 0.084 | 0.079 | 0.076 |
| % ∆CoVaR_U | 51.32 | 38.32 | 60.47 | 61.60 | 43.55 | − 11.21 | 1.23 | 1.23 |
| Rank | 3 | 5 | 2 | 1 | 4 | 8 | 6 | 7 |
The order is based on the value of %∆CoVaR, cprice represents carbon price of Guangdong ETS
The risk spillover from other markets to the national average carbon price
| Cprice-WTI | Cprice-BRENT | Cprice-CER | Cprice-EUA | Cprice-Gas | Cprice-Coal | Cprice-CSI300 | Cprice-industry index | |
|---|---|---|---|---|---|---|---|---|
| CoVaR-D | − 5.03 | − 5.05 | − 5.84 | − 4.80 | − 2.36 | − 0.58 | − 3.10 | − 3.22 |
| ∆ CoVaR-D | 0.00 | 0.00 | − 0.01 | 0.00 | − 0.03 | − 0.00 | − 0.02 | − 0.01 |
| %∆CoVaR-D | 0.00 | 0.00 | 0.12 | 0.00 | 1.14 | 0.51 | 0.53 | 0.30 |
| ranked | 6 | 7 | 5 | 8 | 1 | 3 | 2 | 4 |
| CoVaR-U | 13.61 | 10.28 | 12.25 | 10.65 | 4.54 | 0.96 | 7.99 | 6.72 |
| ∆CoVaR-U | 5.34 | 2.43 | 1.58 | 0.02 | 0.11 | 0.01 | 0.00 | 0.60 |
| % ∆CoVaR-U | 64.53 | 30.87 | 14.77 | 0.19 | 2.56 | 0.71 | 0.00 | 9.82 |
| Ranked | 1 | 2 | 3 | 7 | 5 | 6 | 8 | 4 |
Cprice is the national average carbon price