Next Article in Journal
Chinese Electric Power Development Coordination Analysis on Resource, Production and Consumption: A Provincial Case Study
Previous Article in Journal
Diamondback Moth, Plutella xylostella (L.) in Southern Africa: Research Trends, Challenges and Insights on Sustainable Management Options
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Stochastic Differential Equation Models for the Price of European CO2 Emissions Allowances

1
School of Economics and Management, Fuzhou University, No. 2, Xueyuan Road, Daxue New District, Fuzhou District, Fuzhou 350108, China
2
Department of Mathematics and Statistics, University of Strathclyde, Glasgow G1 1XH, UK
*
Author to whom correspondence should be addressed.
Sustainability 2017, 9(2), 207; https://doi.org/10.3390/su9020207
Submission received: 8 October 2016 / Revised: 23 January 2017 / Accepted: 24 January 2017 / Published: 4 February 2017
(This article belongs to the Section Economic and Business Aspects of Sustainability)

Abstract

:
Understanding the stochastic nature of emissions allowances is crucial for risk management in emissions trading markets. In this study, we discuss the emissions allowances spot price within the European Union Emissions Trading Scheme: Powernext and European Climate Exchange. To compare the fitness of five stochastic differential equations (SDEs) to the European Union allowances spot price, we apply regression theory to obtain the point and interval estimations for the parameters of the SDEs. An empirical evaluation demonstrates that the mean reverting square root process (MRSRP) has the best fitness of five SDEs to forecast the spot price. To reduce the degree of smog, we develop a new trading scheme in which firms have to hand many more allowances to the government when they emit one unit of air pollution on heavy pollution days, versus one allowance on clean days. Thus, we set up the SDE MRSRP model with Markovian switching to analyse the evolution of the spot price in such a scheme. The analysis shows that the allowances spot price will not jump too much in the new scheme. The findings of this study could contribute to developing a new type of emissions trading.

1. Introduction

In 2005, an emissions trading scheme (ETS) was adopted by the European Union (EU ETS) to reduce CO2 emissions by firms. According to the EU ETS, emissions-intensive firms have the right to release a certain amount of CO2 into the atmosphere. This amount of permitted CO2 emissions is allocated or auctioned to firms in the form of European Union Allowances (EUAs). The firm can emit 1 ton of CO2 via one EUA, which is the tradable commodity in the EU ETS.
The participating firms are forced to hold adequate emissions allowances for their outputs. The carbon market provides new business opportunities for market intermediaries, such as emissions-related firms, brokers, traders, and risk management consultants. For these groups, a valid spot price model becomes increasingly important.
Recently, to bridge the gap between theory and actual price behaviour, numerous empirical studies have investigated the time series of the emissions allowance price. Because the EU ETS is by far the largest, most developed market, empirical research has mainly focused on it. Paolella and Taschini [1] advocated using econometric frameworks to explain the heteroscedastic dynamics of emissions allowance prices under the EU ETS. The authors summarized the poor performance of previous methods, such as forecasting analysis based on demand/supply fundamentals, and supported the employment of a well-suited generalized autoregressive conditional heteroscedasticity (GARCH)-type model, which is suitable for depicting the stylized facts of daily returns [1]. By employing econometric testing procedures and trading strategies, Daskalakis and Markellos demonstrated that the market for European CO2 emissions allowances is inefficient [2]. Furthermore, within the framework of a GARCH approach, Oberndorfer claimed that the stock performance of electricity firms has a positive effect on the EUA price [3]. On the other hand, some research has used the stochastic method to understand the stochastic properties of spot volatility. In the study of Daskalakis et al. [4], several different diffusion and jump-diffusion processes were fitted to the European CO2 futures and options time series. This study suggested that the proscription of banking, which implies that emissions allowances cannot be used in the next phase, plays a significant role in determining the prices of EUA futures. Benz and Trück [5], employing a Markov switching model and AR-GARCH models, described the price dynamics and changes of the short-term allowance spot price. Similarly, Seifert et al. [6] built a stochastic equilibrium model to analyse the dynamics of CO2 allowances spot prices. The authors proposed that the CO2 emissions allowance spot price process should have a time- and price-dependent volatility structure. Recently, Kim et al. [7] estimated the dynamics of EUA futures prices based on Heston’s stochastic volatility model with or without jumps, and their empirical results revealed three important features of EUA futures prices: significant stochastic volatility, noticeable leverage effect, and inclusion of jumps.
There is extensive literature concerning the economic and policy aspects of the EU ETS, but there is little explicit study of the dynamic emissions allowance price in the presence of market uncertainty. Moreover, rare previous studies have covered the stochastic nature of emissions allowances in Phase III. There are three periods is EU-ETS protocol. Previous studies focused on Phases I and II, of which the main mechanism was free allocation based on past emissions. However, in Phase III, auctioning is expected to become the dominant allocation mechanism. It is necessary to examine adequate pricing models for EUA prices in this phase. Motivated by shortcomings of previous studies on carbon derivatives, we apply a new simple parameter estimation method for stochastic differential equations (SDEs) and discuss price forecast models based on different market data. We show that the mean reverting square root process (MRSRP) is the best of five SDEs to predict the EUA spot price.
In addition, the current ETS cannot reduce air pollution congestion (e.g., smog). This study is the first existing work to develop a new trading scheme to reduce air pollution congestion and explore the effect of air pollution on the EUA spot price by using the SDE MRSRP model with Markovian switching. The findings reveal that the spot price will not jump much under the new trading scheme when firms have to hand in more allowances on heavy pollution days.
The remainder of the paper is organized as follows. The following Section 2 presents the econometric analysis of EUA spot prices. Section 3 introduces a new parameter estimation methodology for SDEs and use empirical analysis to examine which SDE model is the best for reflecting the EUA spot price. Section 4 sets up an SDE model with Markovian switching to depict the EUA spot price considering the effect of air pollution. Finally, Section 5 concludes.

2. Econometric Analysis of EUA Spot Prices

2.1. Discovery of the Spot Market

Under the EU ETS, there are two active EUA markets: the French Powernext and European Climate Exchange (ECX). In 2007, almost 79% of the EUA spot transactions were handled in the Powernext market, which plays a leading role in EUA price formation [4]. In both markets, spot contracts involve one EUA and are settled the day after the transaction. In our study, we use the datasets of daily settlement prices covering the period from 1 November 2012 to 31 October 2014 for Powernext and ECX, respectively. By comparing price levels from different trading platforms, we are able to study the potential impact of market conditions. Table 1 presents the descriptive statistics of the price levels and returns on EUA spot prices in both markets.
The second (fourth) column of Table 1 shows the time-series mean, median, maximum, minimum, standard deviation (Std. Dev.), skewness, kurtosis, Jarque–Bera test results, and autocorrelation coefficients of EUA prices. The third (fifth) column shows the time-series statistics of the returns of EUA prices. The coefficients of skewness and kurtosis indicate that the distribution of prices is leptokurtic in two markets, which is confirmed by the Jarque–Bera test results. The test leads to rejection of the null hypothesis that the spot price levels and returns come from normal distribution. Additionally, autocorrelation coefficients decrease slowly in price levels, which is consistent with possibly non-stationary, variation. The stationarity of the prices is confirmed through three unit root tests (for a detailed description of these, see [8]). The test results are presented in Table 2. The results of the augmented Dickey–Fuller (ADF) test and the Philips–Peron (PP) test are not statistically significant, which means we cannot reject the null hypothesis that the logarithmic EUA prices have a unit root. The results of the Kwiatkowski–Phillips–Schmidt–Shin (KPSS) test reveal that we can reject the null hypothesis of stationarity at the 1% level. All test results suggest that the logarithmic EUA prices in both markets are not stationary at statistically significant levels.
In addition, the descriptive analyses find that the EUA spot price behaviour in the two trading markets is very similar. This is as expected, because any temporal differences between the two markets will be wiped out quickly by rational arbitrage. Actually, the EUA spot prices in Powernext and ECX moved closely, for the absolute difference of the average mean is tiny. Furthermore, the weekly returns between the two markets have a strong correlation coefficient, at almost 95%.

2.2. Dynamics of Emissions Allowance Spot Prices

On average, stock prices rise because the investor is rewarded for his or her money’s time value over a long period. There will be a substantial increase in EUA prices when specific changes occur in a short period, such as policy or weather conditions. However, EUA prices tend to revert to a normal level in the long run. The properties of spot prices are the result of general mean-reversion behaviour and the spikes in prices caused by the supply and demand shocks.
As discussed previously, the EUA spot prices might fit the validity of the standard Brownian motion process. In line with the research of Chan et al. [9] and Dotsis et al. [10], we study the ability of different popular diffusion continuous-time models to analyse the dynamics of the EUA spot prices. It is necessary for the regulator to explore the dynamics to choose an appropriate pricing model and design a new trading scheme.
Under historical probability P, the following SDE represents the underlying stochastic properties of the EUA spot prices:
d s t = μ ( s t , t ) d t + σ ( s t , t ) d W t
where s t is the EUA spot price in time t, W t is a standard Wiener process, μ ( s t , t ) is the drift, and σ ( s t , t ) is the diffusion coefficient. The drift and diffusion are assumed a general function of the EUA spot price and time. We can obtain several different models by combining many assumptions for the components of μ ( s t , t ) and σ ( s t , t ) . There are five configurations, as follows:
Geometric Brownian motion process (GBMP)
d s t = μ s t d t + σ s t d W t
Square root process (SRP)
d s t = μ s t d t + σ s t d W t
Mean reverting process (MRP)
d s t = k ( θ s t ) d t + σ s t d W t
Mean reverting square root process (MRSRP)
d s t = k ( θ s t ) d t + σ s t d W t
Mean reverting logarithmic process (MRLP)
d l n s t = k ( θ l n ( s t ) ) d t + σ d W t
It is well known that Equations (1) and (2) have been used widely to depict the evolution of stock pricing, options, and commodity price indexes (e.g., [11,12,13]). In Equations (1) and (2), µ is the expected return of the asset per unit of time and σ is the volatility. However, the other three processes have mean reverting drifts. In Equations (3)–(5), k is the speed of mean reversion, θ denotes the unconditional mean, and σ measures the asset price volatility. Bierbrauer et al. [14] used Equation (3) to test electricity spot prices based on one-factor and two-factor models using data from the German EEX market. Equations (4) and (5) have been very popular for describing the interest rate and volatility in the literature (e.g., [15,16]).

3. EUA Spot Price Estimation

3.1. Parameter Estimation for the SDE

In this section, we describe general parameter estimation methodology for the five SDE models, Equations (1)–(5), based on least squares techniques. This method has the same accuracy and efficiency as the more complicated maximum likelihood estimation and is easier to apply [17]. We start with the general formulas from least squares theory and then, we develop formulas for the point and interval estimations for the MRSRP model. The process to develop the estimation for the other four models is similar, and thus, is omitted in this section.
We recall the MRSRP in the form of an Itô SDE:
d s t = k ( θ s t ) d t + σ s t d W t ,
where k, θ, and σ are positive constants and W t is a scalar Brownian motion. We assume that the initial condition S ( 0 ) 0 . To apply a numerical method to SDE (Equation (6)), it needs to be replaced by the equivalent problem
d s t = k ( θ s t ) d t + σ | S t | d W t .
Since the method could break down if negative values were supplied to the square root function, given stepsize Δ t , by applying the Euler–Maruyama method to Equation (7) and setting s 0 = S ( 0 ) , we obtain approximations s n S ( t n ) , where t n = n Δ t can be computed by
s n + 1 = s n ( 1 k Δ t ) + k θ Δ t + σ | s n | Δ W n ,
where Δ W n = W n + 1 W n .
Equation (8) can be rewritten as
y n + 1 = α v n + 1 + k u n + 1 + σ Z n + 1 ,
where y n + 1 =   s n + 1 s n Δ t | s n | , α = k θ , v n + 1 = Δ t | s n | , u n + 1 = Δ t | s n | and Z n + 1 ~ N ( 0 , 1 ) .
This is a multiple linear regression model and since data points s n and stepsize Δ t are provided, we can use a regression theorem to estimate the parameters k, α , and σ .
Rawlings et al. [18] discussed multiple linear regression in general matrix form:
Y = X β + ε ,
where
Y = ( y 1 y 2 y n ) , X = ( 1 x 11 x 12 x 1 p 1 x 21 x 22 x 2 p 1 x n 1 x n 2 x n p ) ,   β = ( β 0 β 1 β p ) ,   ε = ( ε 1 ε 2 ε n ) .

3.1.1. Point Estimations

The calculations work equally well for Equation (9), which can be written in matrix form Equation (10), where Y and ε remain the same while X and β become
X = ( v 1 u 1 v 2 u 2 v n u n ) ,   β = ( α k ) .
Then,
( α ^ k ^ ) =   β ^ =   ( X T X ) 1 ( X T Y ) = 1 v k 2 u k 2 ( v k u k ) 2 ( u k 2 v k y k u k v k u k y k v k 2 u k y k u k v k v k y k ) .
Then, we can obtain point estimations
θ ^ = α ^ k ^ = u k 2 v k y k u k v k u k y k v k 2 u k y k u k v k v k y k ,
k ^ = v k 2 u k y k u k v k v k y k v k 2 u k 2 ( u k v k ) 2 .
We obtain θ ^ by dividing the two point estimations here, which is a sensible thing to do. The results are in accordance with the estimations taken from S S k = S S θ = 0 , where
S S = ( y k k θ v k k u k ) 2

3.1.2. Variance of Estimated Parameters

To obtain the interval estimations for the parameters α and k, we need to calculate the variance of β ^ using the formula
v a r ( β ^ ) = ( X T X ) 1 σ 2 ,
where σ 2 can be estimated using
σ ^ 2 = ( Y X β ^ ) T ( Y X β ^ ) n p ,
where p is the number of parameters, and thus, p is 2 in the MRSRP model. Equation (15) can be simplified as
σ ^ 2 = Y T Y Y T X β ^ n 2 .
Equation (16) can be written as
σ ^ 2 = 1 n 2 ( y k 2 ( y k v k ) α ^ ( y k u k ) k ^ ) .
Substituting Equation (11) into Equation (16) yields
σ ^ 2 = y k 2 v k 2 u k 2 y k 2 ( v k u k ) 2 u k 2 ( v k y k ) 2 v k 2 ( u k y k ) 2 + 2 u k y k v k y k u k v k ( n 2 ) ( v k 2 u k 2 ( u k v k ) 2 )
σ ^ 2 is an asymptotically unbiased estimator to σ 2 (i.e., σ ^ n 2 σ 2   as   n ).
Thus, we obtain
v a r ( β ^ ) = v a r ( α ^ k ^ ) = 1 v k 2 u k 2 ( u k v k ) 2 ( u k 2 u k v k u k v k v k 2 ) σ ^

3.1.3. Interval Estimation for k

The 503 observations in our study is a sufficiently large number for the 95% confident interval (CI) for k to be
k ^ ± 1.96 v a r ( k ) = v k 2 u k y k u k v k v k y k v k 2 u k 2 ( u k v k ) 2 ± 1.96 v k 2 σ ^ v k 2 u k 2 ( u k v k ) 2
Note that as n , this 95% CI tends to
0 T 1 | s t | d t ( s T s 0 ) T 0 T 1 | s | d S 0 T 1 | s t | d t 0 T | s t | d t T 2 ± 1.96 σ ^ n 2 0 T 1 | s t | d t 0 T 1 | s t | d t 0 T | s t | d t T 2
We obtain the daily data of EUA spot prices and can work out the CI for k using Equation (19). We use the same methodology to estimate the parameter of other SDE models. Table 3 displays the estimated parameters, t-statistics results (the latter in brackets) and the Bayesian information criterion (BIC) for the five models under scrutiny for the full data of the study.
The results lead to several interesting insights. First, it is obvious that the mean reverting model is better than the non-mean reverting model regarding parameter significance and the BIC in both markets. The findings indicate that the model’s goodness of fit is increased by the addition of mean reversion. The results are consistent with the previous descriptive results that the EUA spot prices are the non-normality of returns and non-stationarity. Second, for mean-reversion models, the BIC of the MRSRP model is smaller than that of the other two mean-reversion models, and thus, the MRSRP model is slightly better than the MRP and MRLP models. Third, it should be noted that parameter µ in SRP is not statistically significant. Therefore, we further investigate the fitness of the other four models in Section 3.2.

3.2. Spot Price Estimation

To compare the four models further, we perform the following empirical analysis. First, we calculate theoretical prices using software R for the EUA spot prices in Powernext and ECX based on the four SDE models. The period analysed is extended from 1 November 2012 to 31 October 2014. Subsequently, we assess the accuracy of the price models by calculating the mean absolute per cent error (MAPE) between theoretical and actual spot prices. The MAPE expressed as a percentage is defined as
MAPE = 100 N t = 1 N F ( T ) t f F ( T ) t a F ( T ) t a ,
where N is the number of observations, F ( T ) t f is the theoretical spot price, while F ( T ) t a is the actual spot price. In addition, for comparison purposes, we compute the mean squared pricing error (MSE) for the whole period.
As shown in Table 4, the results present substantial pricing errors. In particular, the MAPE for MRSRP is 0.1798% and 0.188% in Powernext and ECX, respectively, and the MSE is 0.032% and 0.035%, respectively. These errors for MRSRP are far smaller than are those for other models. The results suggest that the pricing errors of MRSRP for the spot prices in the two markets are well below those of the other three models. In summary, it is clear that the MRSRP has the best fitness of the five models to forecast EUA spot prices.

4. MRSRP with Markovian Switching

4.1. Effect of Air Pollution on the EUA Spot Price

Air pollution is a major environmental and social problem in almost all large cities worldwide and “smog” became a serious problem in China in 2013. Smog occurs when emissions reach a certain concentration that impairs property, public health, and ecosystems. A recent US study found that about 4000 Chinese people are killed by air pollution per day. Therefore, finding an efficient policy instrument is increasingly prominent in the country’s overall development plans. The EU ETS is one such attempt to alleviate air pollution from emissions via a financial market mechanism. However, the current EU ETS cannot effectively solve the problem of pollution congestion (e.g., smog).
To reduce the degree of smog, we suggest applying a new scheme to update the current ETS. In the current ETS, a firm hands in one allowance to the regulator when the firm emits 1 ton of air pollution. Meanwhile, in the new trading scheme, a firm has to hand in two or more emissions allowances to the regulator, when the firm emits 1 ton of air pollution on heavy pollution days (e.g., daily mean concentration exceeds 58 µg/m3), compared to only one allowance on clean days. In other words, under the new scheme, when pollution congestion (e.g., smog) occurs, the EUA spot price is higher than it is on clean days. Rationally, a firm will reduce emissions in pollution days and ease the degree of smog in the environment. Then, the spot price changes according to the air conditions, which can help to control excessive air pollution from emissions through a self-regulating mechanism. Therefore, it is necessary to examine the pricing models for EUA derivatives in this new environment. To depict how the EUA spot price evolves under the new trading scheme, we apply the SDE MRSRP with Markovian switching.
Air conditions affect the EUA spot price under the new scheme. To describe air pollution levels simply, based on a regulation of the UK Department for Environment, Food & Rural Affairs, a clean day is defined as having a daily mean concentration of less than 58 µg/m3, otherwise, it is a pollution day. The frequency of “clean days” and “pollution days” in the UK in 2013 is shown in the following Figure 1.
As Figure 1 shows, “1” denotes clean days and “2” denotes pollution days. The probability distribution of the duration (in days) from clean to pollution (or from pollution to clean) follows an exponential distribution. Thus, we can apply the Markov chain to describe the switching.

4.2. Markovian Switching

The Markovian switching system was first introduced by Krasovskii and Lidskii [20]. This kind of stochastic model describes different types of dynamic systems that might experience abrupt changes in their parameters and structures. The advantages in modelling have been reported in the literature (see Mao [21], Boukas [22], and references therein).
We let ( Ω , , { t } t 0 , ) be a complete probability space with filtration { 𝓕 t } t 0 , which satisfies the usual constraints (i.e., it is increasing right-continuous while 𝓕 0 contains all null sets). Let r ( t ) , t 0 be a right continuous Markov chain on the probability space, which takes values in finite state space 𝕤 = { c l e a n , p o l l u t i o n } with the following generator
Γ = ( v c p v c p v p c v p c )
Here, v c p > 0 is the transition rate from the state of “clean” to that of “pollution”, while v p c > 0 is the transition rate from the state of “pollution” to that of “clean”, and thus,
{ r ( t + δ ) = p o l l u t i o n | r ( t ) = c l e a n } = v c p δ + o ( δ ) ,
and
{ r ( t + δ ) = c l e a n | r ( t ) = p o l l u t i o n } = v p c δ + o ( δ ) ,
where δ > 0 .
In our case, the Markov chain r ( · ) is independent of Brownian motion B ( · ) .
As is well known, almost every sample path of Markov chain r ( · ) is a right-continuous step function [23]. More precisely, there is a sequence { τ k } k 0 of finite-valued and 𝓕 t -stopping times, such that 0 = τ 0 < τ 1 < < τ k almost certainly and r ( t ) can be expressed as
r ( t ) = k = 0 r ( τ k ) I [ τ k , τ k + 1 ] ( t ) .
Moreover, given that r ( τ k ) = c l e a n , random variable τ k + 1 τ k follows exponential distribution with parameter v c p , namely,
( τ k + 1 τ k T | r ( τ k ) = c l e a n ) = e v c p T , T 0 ,
Meanwhile, given that r ( τ k ) = p o l l u t i o n , random variable τ k + 1 τ k follows exponential distribution with parameter v p c , namely,
( τ k + 1 τ k T | r ( τ k ) = p o l l u t i o n ) = e v p c T , T 0 .
We can easily simulate the sample paths of the Markov chain using the exponential distributions.

4.3. Simulation of the Spot Price in the New Trading Scheme

Having reviewed the details of the Markov chain, we now study the SDE MRSRP model with Markov switching. Our aim is to compute the EUA spot prices under the new trading scheme if the firm has to hand many more emissions allowances to the regulator for pollution days. Then, we can compare the spot prices between the current and new trading schemes.
We assume the system parameter σ is constant. Given that, on clean days, the parameter k c = 0.021 , and on pollution days, the parameter k p = 0.037 . We estimate the transition rate from clean to pollution using v c p = 0.127 and the transition rate from pollution to clean using v p c = 0.347 .
In 2013, there were 250 business days in the EU ETS. Moreover, the starting level of the spot price on 2 January 2013 was 6.52, and we set it as the starting level for 2013, that is, p 1 = 6.52 . In this study, we design a function in R-software to perform simulations of the spot price for the SDE MRSRP model (Equation (4)) with Markovian switching with step time Δ t = 1 . Combined with the spot price of the current scheme without considering the effect of air pollution, we obtain Figure 2.
Figure 2 depicts the computer simulation of the EUA spot price in 2013. The black line is the simulation of the EUA spot price without switching, while the red line is the simulation of the EUA spot price with switching. The former is the EUA spot price under the current EU ETS scheme; the latter is the EUA spot price under the new trading scheme considering the effect of air pollution. As shown in Figure 2, the red line is more stable than the black one. The spot prices do not jump too much under the new ETS. Thus, it would help to improve the ETS and achieve the primary goal—reduction of air pollution by the government at least possible cost.

5. Discussions and Conclusions

This study examines the stochastic nature of emissions allowances under different schemes pertaining to the relevant models that previous studies failed to achieve. First, we analyse the two primary markets for CO2 emissions allowance spot price under the scheme of the EU ETS: Powernext and ECX. In line with previous studies focused on the EUA spot prices in Phases I and II (e.g., [4,5,6]), the empirical analysis provides evidence that EUA spot prices display non-stationary behaviour in Phase III. In addition, the EUA spot prices approximately follow geometric Brownian motion. An empirical evaluation using actual market data demonstrates that the MRSRP has the best fitness of the five SDEs for representing the dynamics of the EUA spot prices.
More importantly, to reduce smog, this study develops a new scheme in which a firm has to hand many more allowances to the regulator when it emits one unit of air pollution on heavy pollution days, compared to only one allowance on clean days. In Section 4, we show that under the new trading scheme, the air condition can be improved in the pollution days. Using real-life data, we find that the time gap between clean days and pollution days follows exponential distribution, and therefore we develop a SDE MRSRP model with Markovian switching to predict the spot prices in the new trading scheme. The model simulation has shown that the spot prices are expected to not jump too much under the new trading system.
We contribute to the literature in three ways. First, we explore the fitness of the SDEs mainly based on the data of the Phase III of ETS, which has rarely been covered in previous studies. Second, this study is one of the earliest works to develop a new trading scheme with two different states to reduce pollution congestion. The theoretical conception would serve as the basis for more in-depth studies in the future and as a tool for formulating policies aimed at reducing pollution disasters. Third, although a lot of research has been done on spot price modelling, none of the existing literature used SDE with Markovian switching to study the effect of the air condition, and this model is better in our study because there are two different air conditions under our new scheme.
There are several important implications arising from this study for the regulators and managers. First, the current CO2 emissions market behaves like the MRSRP model. One explanation is that the supply of EUA is fixed in ETS, when there is an increase in demand due to policy or weather conditions, pushing prices higher. When demand returns to normal levels, prices will fall. Besides the mean reverting characteristic, the spot price is proportion to its variance in the current CO2 emission market [7]. This important property improves the assessment of production costs incorporating CO2 costs since the introduction of emission trading system, or supports emissions-related investment decisions. The ability of managers to predict the EUA spot prices helps to maintain market efficiency and sustain a healthy trading volume. Second, the estimations in this study showed that the EUA price would not fluctuate too much under the new trading scheme, which greatly enhances the political acceptability of the scheme. Specifically, it implies that the new trading scheme does not lead to more complications in the pricing of emission allowances or to the adverse effect on market liquidity and efficiency. We suggest the government to apply our scheme to upgrade the current trading scheme.
The study has methodological limitations. Primarily, we compare only five SDEs to the EUA spot prices, neglecting other SDEs. In addition, for lack of the EUA price, we cannot compare the theoretical and actual spot price under the new scheme. When future data become available in years to come, it would be worthwhile to re-perform the work in greater detail than the current data allow. Finally, more advanced econometric techniques to estimate the model parameters, including the MCMC method, will be widely elaborated on in the future.

Acknowledgments

The work described in this paper is supported by the National Natural Science Foundation of China (71503049) and the Social Science Foundation of Fujian Province (2014C041). We would like to thank three anonymous referees for insightful comments which led to improvements of the paper.

Author Contributions

All authors contributed equally to the writing of this paper. All authors read and approved the final manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Paolella, M.S.; Taschini, L. An econometric analysis of emission allowances prices. J. Bank. Financ. 2008, 32, 2022–2032. [Google Scholar] [CrossRef] [Green Version]
  2. Daskalakis, G.; Markellos, R.N. Are the European carbon markets efficient? Rev. Futures Mark. 2008, 17, 103–128. [Google Scholar]
  3. Oberndorfer, U. EU emission allowances and the stock market: Evidence from the electricity industry. Ecol. Econ. 2009, 68, 1116–1126. [Google Scholar] [CrossRef]
  4. Daskalakis, G.; Psychoyios, D.; Markellos, R.N. Modeling CO2 emission allowance prices and derivatives: Evidence from the European trading scheme. J. Bank. Financ. 2009, 33, 1230–1241. [Google Scholar] [CrossRef]
  5. Benz, E.; Trück, S. Modeling the price dynamics of CO2 emission allowances. Energy Econ. 2009, 31, 4–15. [Google Scholar] [CrossRef]
  6. Seifert, J.; Uhrig-Homburg, M.; Wagner, M. Dynamic behavior of CO2 spot prices. J. Environ. Econ. Manag. 2008, 56, 180–194. [Google Scholar] [CrossRef]
  7. Kim, J.; Park, Y.J.; Ryu, D. Stochastic volatility of the futures prices of emission allowances: A Bayesian approach. Phys. A 2017, 465, 714–724. [Google Scholar] [CrossRef]
  8. Mills, T.C.; Markellos, R.N. The Econometric Modelling of Financial Time Series, 3rd ed.; Cambridge University Press: Cambridge, UK, 2008. [Google Scholar]
  9. Chan, K.C.; Karolyi, G.A.; Longstaff, F.A.; Sanders, A.B. An empirical comparison of alternative models for the short-term interest rate. J. Financ. 1992, 47, 1209–1227. [Google Scholar] [CrossRef]
  10. Dotsis, G.; Psychoyios, D.; Skiadopoulos, G. An empirical comparison of continuous-time models of implied volatility indices. J. Bank. Financ. 2007, 31, 3584–3603. [Google Scholar] [CrossRef]
  11. Merton, R.C. Option pricing when underlying stock returns are discontinuous. J. Financ. Econ. 1976, 3, 125–144. [Google Scholar] [CrossRef]
  12. Dassios, A.; Nagaradjasarma, J. The square root process and Asian options. Quant. Financ. 2006, 6, 337–347. [Google Scholar] [CrossRef] [Green Version]
  13. Kourtis, A.; Markellos, R.N.; Psychoyios, D. Wine price risk management: International diversification and derivative instruments. Int. Rev. Financ. Anal. 2012, 22, 30–37. [Google Scholar] [CrossRef]
  14. Bierbrauer, M.; Menn, C.; Rachev, S.T.; Trück, S. Spot and derivative pricing in the EEX power market. J. Bank. Financ. 2007, 31, 3462–3485. [Google Scholar] [CrossRef]
  15. Detemple, J.; Osakwe, C. The valuation of volatility options. Eur. Financ. Rev. 2000, 4, 21–50. [Google Scholar] [CrossRef]
  16. Bahmani-Oskooee, M.; Kutan, A.M.; Zhou, S. Do real exchange rates follow a nonlinear mean reverting process in developing countries? Sou. Econ. J. 2008, 74, 1049–1062. [Google Scholar]
  17. Pan, J.; Gray, A.; Greenhalgh, D.; Mao, X. Parameter estimation for the stochastic SIS epidemic model. Stat. Inference Stoch. Process. 2014, 17, 75–98. [Google Scholar] [CrossRef] [Green Version]
  18. Rawlings, J.O.; Pantula, S.G.; Dickey, D.A. Applied Regression Analysis: A Research Tool; Springer: New York, NY, USA, 1998. [Google Scholar]
  19. European Environment Agency. Air Quality e-Reporting (AQ e-Reporting). Available online: http://www.eea.europa.eu/data-and-maps/data/aqereporting#tab-data-by-country (accessed on 25 January 2017).
  20. Krasovskii, N.N.; Lidskii, E.A. Analysis design of controller in systems with random attributes—Part 1, Automat. Rem. Control 1961, 22, 1021–1025. [Google Scholar]
  21. Mao, X. Exponential stability of stochastic delay interval systems with Markovian switching. IEEE Trans. Autom. Control 2002, 47, 1604–1612. [Google Scholar]
  22. Boukas, E.K. Stabilization of stochastic nonlinear hybrid systems. Int. J. Innov. Comput. Inf. Control 2005, 1, 131–141. [Google Scholar] [CrossRef]
  23. Anderson, W.J. Continuous-Time Markov Chains; Springer: Berlin/Heidelberg, Germany, 1991. [Google Scholar]
Figure 1. UK air condition in 2013 [19].
Figure 1. UK air condition in 2013 [19].
Sustainability 09 00207 g001
Figure 2. Comparison of EUA spot prices under different schemes.
Figure 2. Comparison of EUA spot prices under different schemes.
Sustainability 09 00207 g002
Table 1. Descriptive statistics of EUA price levels (P) and returns (R).
Table 1. Descriptive statistics of EUA price levels (P) and returns (R).
PowernextECX
PRPR
# Obs.503502503502
Mean5.235−0.00335.229−0.0034
Median5.1200.0105.1200.000
Maximum9.0700.9009.0600.760
Minimum2.700−1.6502.720−1.590
Std. Dev.1.0530.2181.0460.216
Skewness0.306−1.0190.320−1.020
Kurtosis3.14511.0733.15510.385
Jarque–Bera8.310 *1450.438 **9.087 *1227.364 **
ρ ( 1 ) 0.9920.0000.9930.000
ρ ( 2 ) 0.9700.0480.9690.077
ρ ( 3 ) 0.938−0.1010.936−0.130
ρ ( 4 ) 0.910−0.1420.907−0.127
Note: * (**) denotes (denote) significance at the 5% (1%) level; ρ ( t ) are autocorrelation coefficients at lag t.
Table 2. Unit root test on logarithmic EUA price levels.
Table 2. Unit root test on logarithmic EUA price levels.
TestNull HypothesisPowernext ECX
CTCCTC
ADFUnit root−2.780−3.129−2.773−3.441
PPUnit root−2.858−3.134−2.870−3.415
KPSSStationarity0.577 **0.379 **0.748 **0.313 **
Note: C (TC) refers to a constant (a time trend and constant) in the test equation. * (**) denotes (denote) significance at the 5% (1%) level. ADF is the augmented Dickey–Fuller test, PP refers to the Philips–Peron test, and KPSS refers to the Kwiatkowski–Phillips–Schmidt–Shin test.
Table 3. Estimation results of five SDE models: 1 November 2012–31 October 2014.
Table 3. Estimation results of five SDE models: 1 November 2012–31 October 2014.
ParameterPowernextECX
GBMPSRPMRPMRSRPMRLPGBMPSRPMRPMRSRPMRLP
µ0.00058 *
(0.287)
−0.0006
(−0.361)
0.00056 *
(0.278)
−0.0006
(−0.349)
k 0.032 **
(3.287)
0.029 **
(3.126)
0.029 **
(2.952)
0.03 **
(3.372)
0.029 **
(3.148)
0.029 **
(2.940)
0.163 **
(3.421)
0.153 **
(3.122)
0.048 **
(2.900)
0.165 **
(3.506)
0.151 **
(3.141)
0.047 **
(2.888)
σ0.0460.0990.0450.0980.0460.0450.0980.0450.0970.046
BIC−1666−881−1671−1884−1647−1670−895−1676−1899−1649
Note: * (**) denotes (denote) significance at the 5% (1%) level. GBMP is geometric Brownian motion process, SRP is square root process, MRP is mean reverting process, MRSRP is mean reverting square root process, and MRLP refers to mean reverting logarithmic process.
Table 4. Comparison of different models for EUA spot prices.
Table 4. Comparison of different models for EUA spot prices.
PowernextECX
GBMPMRPMRSRPMRLPGBMPMRPMRSRPMRLP
MAPE (%)0.44330.18890.17980.25170.3340.2110.1880.297
MSE (%)0.19650.03560.0320.0630.1160.0450.0350.088

Share and Cite

MDPI and ACS Style

Cai, W.; Pan, J. Stochastic Differential Equation Models for the Price of European CO2 Emissions Allowances. Sustainability 2017, 9, 207. https://doi.org/10.3390/su9020207

AMA Style

Cai W, Pan J. Stochastic Differential Equation Models for the Price of European CO2 Emissions Allowances. Sustainability. 2017; 9(2):207. https://doi.org/10.3390/su9020207

Chicago/Turabian Style

Cai, Wugan, and Jiafeng Pan. 2017. "Stochastic Differential Equation Models for the Price of European CO2 Emissions Allowances" Sustainability 9, no. 2: 207. https://doi.org/10.3390/su9020207

APA Style

Cai, W., & Pan, J. (2017). Stochastic Differential Equation Models for the Price of European CO2 Emissions Allowances. Sustainability, 9(2), 207. https://doi.org/10.3390/su9020207

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop