1. Introduction
Since the Industrial Revolution, carbon dioxide emissions have increased every year, and so far, humans have put more than 1.5 trillion tons of carbon dioxide into the atmosphere. Most greenhouse gases come from heavy industry, vehicle exhaust, and deforestation. Once greenhouse gases exceed atmospheric limits, the greenhouse effect will occur. This will lead to problems such as melting ice caps, extreme weather and rising sea levels, which in turn threaten human existence. As a result, it is imperative for the whole world to effectively reduce greenhouse gas emissions and curb global warming. The entry into force of the Kyoto Protocol has greatly promoted the development of the international carbon market, which has gradually become an important tool to curb global carbon emissions. In order to better manage companies’ carbon dioxide emission reduction, the European Union has established the world’s first carbon emissions trading market on a large scale, the European Union Emissions Trading Scheme (EU ETS). By the end of 2018, 20 carbon markets had been established worldwide, covering eight percent of global greenhouse gas emissions and located in countries and regions that account for 37 percent of the global economy. Emissions trading schemes use the market economy to promote global greenhouse gas reduction and achieve emission peak and carbon neutrality [
1]. Carbon emission trading allows carbon emission rights to be exchanged as a commodity. Therefore, people began to focus on the link between the carbon market and the international financial market, energy, and other commodity markets. Current theoretical studies show that the carbon market is connected with other capital markets mainly through two channels. One is return spillovers, which are linked through the price discovery process. The other is “volatility spillover”, in which a shock in one market can adversely affect the willingness of market participants to hold risk in any market. As well, some studies show that energy assets and non-energy assets affect the carbon market in different ways [
2].
In recent years, most scholars have tended to focus on the carbon-energy markets links. Mansanet Bataller et al. [
3] first revealed the connection between the energy market and carbon price from an economic view. Alberola et al. [
4] verified that all the energy drivers emphasized by Bataller et al. did have an impact on EUAs’ price changes during the whole data period. Zachman [
5] confirmed the existence of nonlinear interaction between fuel, carbon, and electricity price, in which electricity price can be expressed as a linear combination of changes in fuel and carbon price. Hammoudeh et al. [
6] studied the influence of fluctuation between EUA and crude oil price through the nonlinear autoregressive distributed lag (NARDL) model, and then discovered the continuous negative impact of crude oil price on EUA price. Ortas et al. [
7] verified that carbon assets and energy commodities embody different leading/lagging motions at different time frequencies through a method based on wavelet coherence analysis. Zhang et al. [
8] adopted the threshold dynamic conditional correlation GARCH model and the complete BEKK GARCH model to conclude that both the coal-carbon market and the carbon-natural gas market had obvious one-way volatility spillover, and that the carbon market is observably positively correlated with the fossil energy market. Wang et al. [
9] used the DY spillover index to reveal the existence of an asymmetric spillover effect between the carbon market and energy market, and proved that such spillover effect changes over time. Ji et al. [
10] showed that brent crude oil price had a non-negligible impact on EUA price volatility and risk. Using VAR-GARCH model, Lee et al. [
11] found that the spillover effect between EUA and Brent crude oil market was stronger than the volatility spillover effect among EUA, biofuels and Brent crude oil market. Xu et al. [
12] made some studies on the information connection of the carbon-energy market system by introducing a multilayer recurrence network. EI Amri et al. [
1] found that the price determinants of EUA in the third stage of EU ETS included changes in brent crude oil price and coal price, and its price was affected by economic activities and Novel Coronavirus disease, etc.
With the integration and financialization of the world economy, more and more scholars tend to study the carbon market from a financial perspective in view of the obvious property of the carbon market as a financial asset [
13]. Most of the current research on the “carbon-financial” market system focus on the relationship between the carbon market and stock market. Specifically, Oberndorfer [
14] studied the development of the energy market from the perspective of the stock market. Chevallier [
15] then demonstrated the weak correlation between carbon prices and stocks and bonds prices. Subsequently, Chevallier [
16] confirmed the connection between macro economy and carbon price by fitting the FAVAR model with macroeconomic, financial, and commodity data. Tian et al. [
17] showed that strong market shocks to a large extent drove the relationship between the EUA market and stock returns of power companies, and the fluctuations of the two markets were remarkably consistent. Oestreich et al. [
18] demonstrated the explanation of carbon risk for cross-sectional changes in stock returns and the apparent carbon premium in stock returns.
At present, it has become a hot issue to integrate carbon market, energy market and financial market into one system and discuss the relationship between the three. The following is a brief review of relevant literature on the carbon-energy-financial system. Creti et al. [
19] studied the decisive factor of carbon price in the two phases of EU ETS and found that energy price was a key factor of carbon price, and the increase in oil price and conversion price would lead to the increase of carbon price. The growing role of carbon markets in the economy is also confirmed by studying the influence of carbon prices on stock prices in this paper. Koch [
20] based on the MGARCH model of smooth transition, studied the price connection between carbon allowance and market fundamentals in EU ETS and found that there are connections between the carbon market and financial market. Song et al. [
21] explored the spillover relationship between the price of EUA futures and the price of 19 financial, energy, and other asset categories through the ARMA (1,1)-CGARCH model. The research results show that there are volatility risks and information transmission among the carbon futures market and the financial, energy and other asset categories, in which the relationship with the financial market is slightly stronger, and this transmission relationship is dynamic with the development of the EU carbon market. Tan and Wang [
22] revealed that the carbon price was significantly affected by energy prices and financial risk factors in the three stages of EU ETS through quantile regression. By studying the price determination mechanism of the carbon market, Ji et al. [
23] found that the energy market influences enterprise emissions and then carbon market through changing energy prices, and financial market influences carbon market through price signals. Tan et al. [
2] investigated the system-wide correlation of returns and volatility among EUA and selected commodities and the two financial market assets in the “carbon-energy-financial” market system by the modified Lanne-Nyberg DY spillover index model. It is found that EUA has a higher correlation with energy markets than with financial markets, and the carbon market has a closer correlation with the stock market among the connection between carbon market and each financial asset. Liu et al. [
24] systematically investigated the return spillover effect among the carbon market, renewable energy markets, financial markets and energy markets from 2008 to 2021, and found that the carbon market was increasingly linked to the world capital markets. Yao et al. [
25] constructed a carbon-energy-stock system to compare information spillovers among the three subsystems under a unified framework, and adopted a connectivity network to identify the roles and positions of carbon market, energy market, and stock market, and discussed the dynamic evolution process of information spillovers.
As shown above, these literatures have explored the interdependence between carbon, conventional energy, and financial assets, but there are still some issues that need to be explored. Most of these studies have focused on two of these markets, and the systematic study of the linkages between “carbon-energy-finance” markets is incomplete. We can see that the relationship between markets is currently mainly studied through econometric models, such as the VAR model, GARCH model, wavelet method, Granger causality test, etc. Most of these econometric methods are based on smooth, linear data and are easily affected by the type and distribution of data. With the construction and development of carbon market, carbon price fluctuation is not only affected by the internal market mechanism, but also affected by the heterogeneous environment. Carbon price fluctuation presents a complex nonlinear dynamic evolution process. How to scientifically evaluate the role of carbon market in the whole social and economic system? How to accurately measure the signal transmission between carbon market price information and other markets? An in-depth analysis is needed from a new perspective.
For the past few years, recurrence quantitative analysis and recurrence networks have shown certain advantages in dealing with nonlinear time series [
26,
27,
28,
29], which have been widely used in finance, energy, engineering, life science and other fields to reveal nonlinear, autocorrelation, fractal and other characteristics of the system. In order to systematically analyze the linkage relationship of carbon-energy-financial market from the perspective of network, we made the following contributions in this study. Firstly, we combine multilayer networks with recurrence networks, then propose a method for constructing time-varying multilayer recurrence networks (TMRN) and delay time-varying multilayer recurrence networks (DTMRN), and further propose a new correlation index. Compared with traditional econometric methods, this method does not need to consider the data distribution and whether the data is stable. It also applies to nonlinear data, which is obviously more suitable for the complex data types in the current era of big data. Secondly, we use the construction of multilayer recurrence networks to integrate carbon market, energy market and financial market into a system to analyze the relationship between them. Using the two indexes we propose, we examine the correlation and lead relationship among the three markets. By comparing the results of the three stages of the EUcarbon market construction, we reveal the evolution process of the increasingly close connection between the three markets, and also find that the carbon market is playing an increasingly important role in the carbon-energy-financial system. Finally, we quantitatively describe the evolution of the price signal of the EU carbon market gradually coming into play, and reveal the information spillover effect among the carbon market, energy market, and financial market. Our results may provide some new ideas for policymakers to formulate policies in the process of carbon market construction and carbon neutrality. As the carbon market matures, it also provides an alternative asset investment channel for investors. Given the increasingly close links between the carbon market and other markets, it is reasonable to suggest that market participants should consider the development of the carbon market when investing in order to maximize profits. At the same time, based on the results of the leading relationship, that is, the carbon market leads the energy market and financial market in the current phase, it can also provide some basis for investors to make decisions. Investors can pick up on signals from the carbon market to optimize corporate portfolios [
23,
25,
30,
31].
The rest of the article is laid out as follows. In the second section, we present the research method used in this paper. In the third section, we conduct a numerical simulation to verify the effectiveness of the network-based indicators. In the fourth section, we inquire into the correlation and guiding relationship between carbon market, energy market and financial market by using the multilayer recurrence network. In the fifth section, we give the conclusion.
2. Materials and Methods
Numerous studies have indicated that carbon price, energy price and stock price of financial market all have chaotic characteristics to varying degrees [
32,
33,
34]. Phase space reconstruction is an effective method to recover and characterize the prime dynamical system from time series data with chaotic characteristics [
35]. In this paper, the proposed method is based on phase space reconstruction theory, the basic idea is as follows: Firstly, the phase space of carbon price, energy price and stock price data of financial market is reconstructed. Then, each phase point in the phase space is treated as a node of the network, and the connection between nodes is established under a certain threshold value, so as to construct the recurrence network of carbon price, energy price and stock price. Finally, the correlation analysis and guidance relationship analysis models are constructed respectively by using the microscopic topological structure and degree distribution characteristics of the network. In this section, we first review the definition of recurrence network and multilayer recurrence network, then introduce the construction steps of time-varying multilayer recurrence network and correlation indicators, and finally explain the construction steps of delay time-varying multilayer recurrence network and leading relationship indicators.
2.1. The Definition of RN and MRN
The recurrence network (RN) is defined as follows [
36]: Let
be a set of time series with
real numbers, then it can be reconstructed as a locus
with time delay embedding in its phase space, where
represents the embedding dimension and
represents the embedding delay. For the trajectory
, define the adjacency matrix
of RN, where, if
falls in the neighborhood of
(
-radius sphere), then
and
are connected, that is,
, otherwise
.
In order to study multivariable time series, the multi-layer recurrence network (MRN) method is introduced by combining recurrence network with multilayer network method [
37], which is defined as follows: Let
be a m-dimensional multivariate time series, where
. Then, we can construct a multilayer recurrence network by creating the RN of the
component of
, which corresponds to the
layer of the multilayer network. Its adjacent matrix is denoted as
, and if nodes
and
are connected in the
layer, then
otherwise,
.
2.2. The Correlation Measurement Based on TMRN
Carbon price, energy price and stock price of financial market show complex fluctu-ation characteristics over time. In order to describe the fluctuation characteristics of the several market prices in different periods, time parameter is introduced to construct time-varying multilayer recurrence network (TMRN).
Let be a multivariate time series data, where represents the number of time series, represents the length of each time series. The steps for constructing the time-varying multilayer recurrence network are as follows:
Step1: Build sliding window. Firstly, we determine the sliding window size and setting the sliding step , then we can obtain time window, where is Gauss function.
Step2: Construct time-varying single-layer recurrence network. The phase space is reconstructed in each time window and the corresponding recurrence network is obtained. The adjacency matrix on the layer in the sliding window is defined as . If the phase space distance between node and node on layer in the th sliding window is less than the threshold value, a connecting edge is established between node and node , i.e., , otherwise, .
Step3: Construct time-varying multilayer recurrence network. The recurrence network obtained by mapping each time series is regarded as a layer, and the layers are connected according to the corresponding time, and then the time-varying multilayer recursive network is constructed. The adjacency matrix
of time-varying multilayer recurrence network in the
sliding window can be expressed as:
where
is the identity matrix.
From the constructed TMRN, referring to indexes presented in the literature [
12], we propose the following indexes to measure the time-varying correlation between time series, and the calculation formula is
where
and
,
. If at time
, the network connection structure of each layer is completely the same, then
. If the connection structure of each layer is completely different at time
, then
.
Then, based on the degree distribution of each layer network nodes and the information entropy theory, we can define the mutual information between layers at time
, and the equation is as follows:
where,
,
and
respectively represent
layer and
layer in the multi-layer recurrence network,
represents the degree distribution of
layer network at time
, and
represents the joint degree distribution of network at
layer and
layer at time
. The equation for calculating the average mutual information of the whole system is as follows:
Therefore, we can measure the dynamic correlation between time series by using the two indicators defined from the micro and macro levels respectively. As well, compared with the traditional correlation measurement method between time series, the measurement indexes based on Equations (1)–(3) reflect the information of original time series through the network topology, including both linear and nonlinear information.
2.3. The Causality Measurement Based on DTMRN
The relationship between systems is very complex, so in addition to studying the interaction between systems, sometimes we also need to determine whether there is a causal relationship or mutual guiding relationship between systems. Therefore, based on the construction of time-varying multilayer recurrence complex network, we introduce the time-delay parameter , and then give the construction steps of delay time-varying multilayer recurrence network as follows:
Step1: Determine the time delay and build sliding window, and then we can obtain time window.
Step2: Reference time series and contrast time series are selected, and then recurrence network of reference time series at time and contrast time series at time or are constructed respectively.
Step3: We take the recurrence network mapped from the reference time series at time and the recurrence network mapped from the contrast time series at time as one layer of the multi-layer network respectively, and then we build the edge between layers on the basis of the corresponding time. Thus, the delay time-varying multilayer RN is constructed.
In this way, we obtained the construction schematic diagram of delay time-varying recurrence network mapped from multivariate time series, as shown in
Figure 1.
To measure the time-delay and time-varying correlation between time series, we take layer
as the reference layer and layer
as the contrast layer, and calculate the time-delay and time-varying correlation index between layer
and layer
as follows:
where
; and
and
are the adjacency matrix elements of the net- work of layer
at time
and layer
at time
respectively. Obviously
. If the network connection structure of layer
at time
and layer
at time
are exactly the same, then
; if the network connection structure of layer
at time
is completely different from that of layer
at time
, then
. According to Equation (4),
Therefore, the time delay correlation of the whole system can be expressed in the following matrix form:
where the main diagonal element represents the time-delay correlation coefficient of each time series.
Then, to measure the leading relationship among various time series, we define the following indexes:
Apparently, , and . If the two series are in complete synchronization at time , then ; if time series is ahead of the time series, then , otherwise, . Therefore, using the magnitude of and the positive value or negative value of , we can measure the mutual guiding relationship between the two sets of time series in different time periods.
Similarly, based on the degree distribution of network nodes of each layer and the information entropy theory, we can define the mutual information between layers with time-delay
at time
, and the calculation equation is as follows:
where,
,
,
and
represent layer
and layer
in the multi-layer recurrence network respectively, and
is the degree distribution of
layer network at time
, and
is the joint degree distribution of the
layer at time
and the
layer with delay
. Using Equation (9), similarly, we can define the measure index of mutual guiding relationship of multivariate time series based on mutual information between layers.
Similarly, if they are in synchronous state, then gets the maximum value; If the time series is ahead of the time series, then , otherwise, . Therefore, with the help of the value of and the positive or negative value of , we can also measure the mutual guiding relationship between the two sets of time series in different time periods.
To sum up, we can obtain the measurement flow chart of correlation (including correlation relation and leading relation) between systems based on multi-level recurrence network, as shown in
Figure 2.
5. Discussion and Conclusions
In the above analysis, we constructed the recurrence network of the EU carbon market, energy market and financial market respectively in different periods. The recurrence network structure of the three markets showed different topological characteristics in different periods, and the linkage relationship among the three markets also showed different characteristics, as shown in
Figure 15. We find that the relationship between carbon, energy and financial markets is getting closer and closer with the construction of the EU ETS. Specifically, the linkage index of carbon-energy increases from 0.4161 in Phase-I to 0.6265 in Phase-III, an increase of 51.79 percentage points. The linkage index of carbon-financial increased from 0.4330 in Phase-I stage to 0.7455 in Phase-III stage, an increase of 72.17 percentage points. The connection between energy market and financial market is also getting closer, rising from 0.5711 in Phase-I Phase to 0.5983 in Phase-III Phase, an increase of 4.76 percentage points. Carbon markets are playing an increasingly important role in the carbon-energy-financial system. From that energy market and financial market led the carbon market in Phase I, to that carbon market led energy market in Phase II, and then to that carbon market led energy market and financial market in Phase III, we can find that the carbon price can gradually effectively reflect immediate information of carbon market and then affect energy market and financial market through the change of the price with gradual improvement of various working mechanisms of EU carbon market.
Based on the construction methods of TMRN and DTMRN of multivariate time series data, the indexes of inter-system correlation measurement and inter-system guidance measurement are given in our work, and the validity of them is verified by numerical simulation. In the empirical analysis part, we use the price data of the European carbon market, brent crude oil in energy market, and S&P 500 in financial market from 7 July 2008 to 30 November 2021 as sample data. In order to facilitate comparison, we divided selected data into three parts according to different stages of EU carbon market construction, and constructed a multilayer recurrence network at each stage, and then used the proposed indicators to study the correlation of the three markets in different periods. We draw the following conclusions: (1) as for the correlation, the correlation between the three markets becomes closer with the development of carbon market. (2) For the leading relationship, energy market and financial market lead carbon market while financial market leads energy market in Phase I; Carbon market leads energy market while financial market leads carbon market and energy market in phase-II; The carbon market leads the energy market and the financial market while the financial market leads the energy market in Phase-III. As time goes by, the carbon market gradually integrates into the world capital market and plays a dominant role in the market system.
Here we compare the results with those of other relevant studies. Comparing the three phases, we find that the correlation between markets obviously increases with the development of EU carbon market, which is basically consistent with research results under the market system of “carbon-renewable energy-finance-energy” of literature [
24]. Our research found that there is information spillover between carbon market, energy market and financial market, and the direction of information spillover changes dynamically with the development of EU carbon market, which is obviously consistent with the results of literature [
21]. Tan et al. found in the “carbon-energy-finance” system that the correlation between EUA and energy assets was often higher than that of financial assets. In our research results, the correlation between carbon market and energy market was greater than that of carbon market and financial market in the early phase of Phase-I and Phase II, while it was opposite in the late phase of Phase-I and Phase III. For the former, we consider that the transition period of the carbon market leads to the instability of the connection, while for the latter, we consider that the difference in data selection leads to the instability. But in general, it does not affect the consistency of the core results. Overall, our results are in good agreement with those in the existing literature. Further, compared with previous studies [
20,
25], based on network topology indicators, we quantify the linkage relationship between the EU carbon market at different phases of development, energy market and financial market, quantitatively measure the position of the three markets in different phases in the whole system, and dynamically evaluate the evolution process of interaction between the three markets.
Based on the above research results, we give the following suggestions:
First, for market decision-makers, our results can provide some ideas. As a market-oriented means to promote energy conservation, emission reduction and green development, carbon trading aims to maximize the input and output of carbon reduction in the whole society by using the price discovery function of the market mechanism. Therefore, government departments should pay attention to information spillover of carbon market to other markets, and promote other markets to fulfill emission reduction requirements by utilizing the association between markets. As the core of carbon market, carbon pricing is also an effective policy tool. Therefore, government departments should consider the connection between markets when setting carbon prices based on the market, which can not only stimulate market productivity and innovation, but also achieve low-cost emission reduction. The gradual maturity of the EU carbon market also provides an important reference for countries to establish carbon trading market. In the process of national carbon market construction, government departments need to consider the operation of domestic carbon market and EU carbon market at the same time, so as to formulate reasonable business plans and regulatory policies. In addition, in the process of building carbon markets, the government should also explore whether there is a risk transmission mechanism between carbon markets and other markets by using the correlation between markets.
Second, we also have some advice for companies and investors. As the price fluctuation of the carbon market is affected by the prices of other markets, enterprises can receive signals from the carbon market through the connection between markets, and then make appropriate production plans to ensure that carbon emissions are reduced while production technology and energy efficiency are improved. In addition to completing the trading of carbon emission rights, a more important function of carbon market is to guide investment and control risks through carbon price signals. Our results indicate that there is information spillover from carbon market to energy market and financial market in the current stage. Therefore, investors who participate in other markets should also pay attention to the price fluctuation of carbon market and its signals, and timely adjust their investment strategies to maximize their interests. In addition, investors should maintain risk awareness, pay attention to the signals of carbon price fluctuations, and timely avoid risks.
Compared with the traditional indexes that measure the correlation and guiding relationship between systems, the indexes in this paper do not need to consider the distribution state and formal representation of data. Different from previous literature that focused on the relationship between the two markets, this paper systematically studied the relationship between the carbon market, energy market and financial market, including both correlation and guiding relationship. Of course, there are some limitations to our work. Although numerical simulation and empirical analysis have demonstrated the superiority of our proposed index, we have not strictly proved theoretically that it meets several properties of big data association metrics, such as equitability, monotonicity, etc., which will be the work we will focus on later. We also need to further verify whether the leading relationship measured by our indexes is economic causality. In addition, for the empirical part of studying the connection between carbon market, energy market and financial market, the data we selected are not sufficient, especially in Phase III, the fourth phase of the construction of EU carbon market. In the future work, we will further optimize the proposed index to adapt it to more types of data, and ulteriorly promote and apply the method in this paper to other fields, thereby enrich and complete the research system of carbon market and other economic fields based on complex network theory.