Fourier quantile unit root test for the integrational properties of clean energy consumption in emerging economies

Fourier quantile unit root test for the integrational properties of clean energy consumption in emerging economies

Accepted Manuscript Fourier Quantile Unit Root Test for the Integrational Properties of Clean Energy Consumption in Emerging Economies Yifei Cai, Ang...

1MB Sizes 0 Downloads 48 Views

Accepted Manuscript Fourier Quantile Unit Root Test for the Integrational Properties of Clean Energy Consumption in Emerging Economies

Yifei Cai, Angeliki N. Menegaki PII: DOI: Reference:

S0140-9883(18)30452-3 https://doi.org/10.1016/j.eneco.2018.11.012 ENEECO 4221

To appear in:

Energy Economics

Received date: Revised date: Accepted date:

7 May 2018 15 September 2018 18 November 2018

Please cite this article as: Yifei Cai, Angeliki N. Menegaki , Fourier Quantile Unit Root Test for the Integrational Properties of Clean Energy Consumption in Emerging Economies. Eneeco (2018), https://doi.org/10.1016/j.eneco.2018.11.012

This is a PDF file of an unedited manuscript that has been accepted for publication. As a service to our customers we are providing this early version of the manuscript. The manuscript will undergo copyediting, typesetting, and review of the resulting proof before it is published in its final form. Please note that during the production process errors may be discovered which could affect the content, and all legal disclaimers that apply to the journal pertain.

ACCEPTED MANUSCRIPT

Fourier Quantile Unit Root Test for the Integrational Properties of Clean

T

Energy Consumption in Emerging Economies

IP

Yifei Cai1

Economics, Business School, University of Western Australia,

CR

Perth, Australia

[email protected]

US

Angeliki N.Menegaki

TEI STEREAS ELLADAS, UNIVERSITY OF APPLIED SCIENCES Economics and Management of Tourist and Culture Units

AN

Amfissa, Nea Poli, Greece

AC

CE

PT

ED

M

[email protected]

1

Corresponding author, Yifei Cai acknowledges that this paper was carried out while the author

was in receipt of a “University Postgraduate Award (UPA) and an Australian Government Research Training Program Scholarship (RTP)” at The University of Western Australia. [1]

ACCEPTED MANUSCRIPT

Fourier Quantile Unit Root Test for the Integrational Properties of Clean Energy Consumption in Emerging Economies

T

Abstract

IP

Since the adoption of the Kyoto protocol in 1997 and its entry into force in 2005, as well its aftermath such as the

CR

Doha amendment and Paris agreement, national policies have become more conscious of the usage of clean energy, mostly the different forms of renewable energy and nuclear energy. Ratifying countries and signatories had

US

commited themselves to binding targets for the reduction of greenhouse emissions by 8% with respect to 1990 levels until 2012, also based on the particular contribution to global emissions from each country. This paper

AN

examines the integrational properties of clean energy consumption from eight emerging economies which are also high greenhouse gas emitters. The empirical results show that the clean energy consumption is stationarity for

M

Brazil and Philippines by using a quantile unit root test without smooth breaks (Koenker and Xiao, 2004).

ED

However, after capturing the smooth breaks (Bahmani-Oskooee et al., 2014), we find the clean energy consumption of China, Pakistan and Thailand are stationary. The time-varying deterministic trend with smooth

PT

breaks are more fitted to the path of clean energy consumption in comparison to the deterministic trend without smooth breaks. The paper suggests economic insights useful for policy making.

AC

CE

Keywords: Clean energy; Emerging economies; Fourier quantile unit root; integration;

[2]

Introduction

M

1.

AN

US

CR

IP

T

ACCEPTED MANUSCRIPT

Within a framework of soaring global temperature, developing clean energy, namely without carbon

ED

emissions, has been viewed as the main policy to curb the CO2 emissions in the long-run. Previous studies have contributed much to investigate the role of foreign direct investment on clean energy use

PT

(Lee, 2013 and Paramati et al., 2016) and the nexus between clean energy consumption and economic

CE

growth (Cai et al., 2018). However, no studies focus on the integrational properties of clean energy use, although numerous papers test the integration of energy consumption (including primary energy

AC

consumption, coal consumption, natural gas consumption and renewable energy consumption) and CO2 emissions (Narayan and Smyth, 2007; Chen and Lee, 2007; Liddle, 2009; Hasanov and Telatar, 2011; Meng et al., 2013; Gozgor, 2016; Payne et al, 2017). The changes of industrial structure always affect the energy consumption structure. In other words, the heavy industries emit more carbon emissions in comparison to agriculture and service sectors because these industries consume a large amount of fossil fuel to support their long-run growth. With the rising of emerging market countries, these countries are experiencing a significant industrial structure change. The energy consumption structure should, without doubt, be altered with the polluting [3]

ACCEPTED MANUSCRIPT sectors becoming reduced or eliminated. In other words, clean energy will be paid increasingly more attention by countries, worldwide, due to their lower carbon emissions. However, the clean energy consumption always contains many structural breaks which significantly affect their stochastic convergence. This occurs because of the investment and the sunk costs taking place for renewable energy investments at some time and the subsidies they are usually accompanied with. This in turn,

T

entails changes (usually increases) in clean energy production which follows the new investments and

IP

integration into the electricity systems. In fact, many scholars (Lee and Chang, 2008 and Payne et al.,

CR

2014) have shown that the structural breaks in the energy consumption are the main threat against robust estimation. Many empirical papers use dummies to approximate the sharp structural breakpoints

US

in the series. However, Enders and Lee (2012) suggest smooth breaks which could be approximated by the Fourier function. Moreover, the stochastic convergence might perform differently at different

AN

quantiles.

Thus, this paper utilizes a newly proposed Quantile Unit Root test with Fourier Function of

M

Bahmani Oskooee et al. (2018) to examine the integrational properties of clean energy consumption for

ED

8 Emerging Market (EM) countries. This paper provides not only integration over the whole quantiles, but also the integration at each quantile. The merits of the Fourier quantile unit root test could be

PT

summarized as follows. First, the Fourier quantile unit root test could provide asymmetric dynamics by allowing different mean-reverting speed at different quantiles of clean energy consumption. Second, the

CE

asymmetries revealed by the quantile regression do not rely upon a special assumption about the functional form of nonlinearities and asymmetries. Third, the conventional unit root tests and even

AC

some which consider structural breaks, all focus on testing on the conditional central tendency. However, through the quantile regression, we could test on clean energy consumption at different quantiles, otherwise we would investigate the whole distribution of clean energy consumption. Specifically, the central tendency is broken down based on the quantile regression. Fourth, many shocks such as wars and political unrest or instability and regulatory policies would make the clean energy consumptions of many countries experience structural breaks. Fifth, because of the structural breaks occurring to clean energy consumption, it is reasonable to test the integrational properties in a framework taking structural breaks into account. Using the Fourier function in the trend function, those [4]

ACCEPTED MANUSCRIPT structural breaks could be approximated as a smooth process. Sixth, the clean energy consumption for some countries such as China, India, Indonesia, Malaysia and Thailand is not subject to normal distribution. As Koenker and Xiao (2004) noted, the non-normality of the series would result to the loss of power by traditional unit root tests2. Overall the originality of our paper is three-fold: First the clean energy (per se)

T

integration properties have not been examined so far. Literature has provided abundant

IP

examples on energy use in general and very few examples of renewable energy, but the clean

CR

energy has not been investigated as such. Second, the sample of emerging economies has not been investigated in relation to this concept either. The investigation of emerging economies is very

US

important, because they have a high economic growth which is fueled mostly by conventional energy, while clean energy is still under a development process and has not stabilized to a certain

AN

amount. Last, the employed econometric method is novel and recently developed. Thus, the paper

M

combines a unique contribution with all these three novel aspects: topic, sample, method. The remainder of the paper is organized as follows. Section 2 reviews the existing studies

ED

related to the test of convergence of energy consumption. Section 3 presents datasets and descriptive statistics. Section 4 introduces the econometric models. Section 5 shows the empirical results. The last

2.

CE

PT

section concludes the paper.

Literature Review

AC

The literature of energy consumption convergence has started more than twenty years ago with the work by Nilsson (1993) and Goldemberg (Goldemberg 1996). These studies found evidence that countries tended to converge to a common pattern of energy use. Afterwards, a number of studies has proliferated both in energy economics and in the broader field of environmental and resource economics with the study of carbon emissions convergence.

2 Meng et al. (2013) utilize the new LM tests based on residual augmented least square regression to solve the problems of non-normal errors when implementing nonlinear analysis.

[5]

ACCEPTED MANUSCRIPT Although convergence per se has been studied in many ways and with different methods, for example with the distribution approaches such as beta convergence (catch up process from high to low convergence), gamma convergence (intra generation mobility) or sigma convergence (spread of distribution). Another strand of energy convergence studies which we will not handle in this paper is the papers

T

which study convergence by means of decomposition analysis. Theil decomposition analysis, Divisia

Decomposition methods separate among various effects such as within and between group

CR

reported.

IP

Decomposition analysis, Laspeyers or Fisher Index decomposition analysis are the most frequently

formation and an account of within group inequalities.

US

decomposition, contribution of energy transformation indexes, contribution of energy to GDP

The interest on whether energy consumption per capita is stationary is motivated by the need to

AN

know about the effects shocks have on it. If a unit root is present, a permanent effect is due after the

M

shock occurence. There are numerous papers that examine the stationarity or non-stationarity of energy consumption and other variables typically perused in the energy-growth nexus. Actually, the

ED

investigation of stationarity of energy consumption is the first step in every energy-growth nexus study (e.g. Magazzino 2014; 2016; 2017). Hence, practically, stationarity studies would be at least as many

PT

in number, as the number of energy-growth studies.

CE

However, the current literature review will not include pure energy-growth nexus studies per se. It includes only studies whose main and exclusive focus is to investigate stochastic convergence of

AC

energy consumption (through unit root investigation) and does not include studies that fall within the energy-growth nexus field whose focus is mainly on the causality analysis. A succinct literature review for this type of studies up to 2012, is provided by Smyth (2013). To make our search yardstick in the current literature review more concrete, we have employed Scopus bibliographic database (www.scopus.com) with keywords: “energy” plus “structural breaks” or “nonlinearity”, “unit root” and “fractional integration” in their titles, keywords and abstracts and have identified 27 studies on the investigation of the convergence of energy variables which appear to be mostly energy consumption per capita, total energy consumption (in aggregate or disaggregate form). [6]

ACCEPTED MANUSCRIPT None of the studies is concerned with the study of the convergence of clean energy consumption as such. The collected studies are shown in Tables 1 and 2. We provide separately studies with aggregate magnitudes (Table 1) and studies with disaggregate magnitudes. (This is also advised in Menegaki

AC

CE

PT

ED

M

AN

US

CR

IP

T

and Tsani (2018); Tsani and Menegaki (2018))

[7]

ACCEPTED MANUSCRIPT

Table 1. Energy consumption convergence studies Study

Period

Variable

Test

Countries

[1]

[2]

[3]

[4]

[5]

Narayan & Smith (2005)

1979–2000

Energy consumption per capita

T P

182 countries

ADF

Result

Agg/Disagg

[6]

[7]

I R

1

Univariate tests (with one or 2 breaks): No stationarity

1

Univariate test with break(s): Rejects H0 of unit root for 56 countries namely 31% of the sample

C S

Panel data unit root test: Stationarity (yes)

Ozcan (2013)

Hsu et (2008)

al.

1980-2009

Energy consumption per capita

1971-2003

Energy use

17 Middle East countries

LM (LS and IPS)

Panel SURADF cross-sectional and identifies countries contain root)

D E

M

A

U N

(hosts effects which a unit

LM tests: Stationarity (yes)

84 countries

Conventional unit root tests reveal stationarity even if only one series in the panel is strongly stationary.

1

Chen & Lee (2007)

1971-2002

Energy consumption per capita

CBL panel unit root

104 countries (7 regional panels)

Non stationarity (even with proper consideration of cross sectional correlations)

1

Mishra et al. (2009)

1980-2005

Energy consumption per capita

CBL panel unit root

13 Pacific island countries

60% of individual countries: Stationarity

1

Ozturk & Aslan (2011)

1970-2006

Energy consumption per capita

LM (IPS)

Turkey

LM univariate test without break Non-stationarity except for the residential sector

1

Narayan et al. (2008)

1973-2007

Energy consumption (sectoral)

LM (LS)

Stationarity except for South Australia

1

T P

C A

E C

(7 sectors) Australia states in industries)

[8]

(6 9

Whole sample: Stationarity

ACCEPTED MANUSCRIPT

Study

Period

Variable

Test

Countries

Result

[1]

[2]

[3]

[4]

[5]

[6]

Aslan & Kum (2011)

1970-2006

Energy consumption

LM (LS)

Turkey sectors)

(7

T P

Linear Unit root is rejected in 4/7 sectors (transport, non-energy use, final energy consumption, cycle & energy).

I R

Agg/Disagg [7] 1

LM test reveals stationarity with two structural breaks

Hasanov & Telatar (2011)

1980-2006

Primary consumption capita

energy per

Kum (2012)

1971-2007

Energy consumption per capita

LM

LM (FLM)

Yilanci & Tunali (2014)

1960-2011

Energy consumption per capita

Zhu & Guo (2016)

1993-2013

Nuclear consumption capita

Lean &Smyth (2014)

1978-2010

Kula al.(2012)

1960-2005

Kula (2014)

et

1923-2008

D E

T P

energy per

E C

C A

ADF, LM (KSS, S)

Simultaneous use of structural breaks and nonlinear unit root tests contribute to the higher probability of evidencing stationarity

1

15 East Asia & Pacific countries

No break

1

109 countries

FLM Stationarity for 24 countries after acknowledging the unknown nature of breaks. For the remaining countries, the standard LM test is use

1

27 countries

If no breaks Non-stationarity(low, upper-middle,developing countries)

0

U N

A

M

C S

178 countries

CBL panel unit root, Longitudinal clustering

One break

Nonstationarity Stationarity

If breaks Stationarity(total, lower-middle,high, developed countries)

Final energy demand (10 fuel types in 4 sectors)

LM (SP, LS)

Malaysia

Stationarity for 50%-70% of energy types in 25%-50% of sectors (no break, one break and two breaks)

0

Electricity consumption capita

LM (LS)

23 high income countries

Stationarity for 21/23 countries

0

LM (LS)

Turkey

Stationarity

0

Electricity

per

[9]

ACCEPTED MANUSCRIPT

Study

Period

[1]

[2]

Lean & Smyth (2009)

Variable

1973m1- 2008m7

[3] consumption capita

Test

Countries

Result

[4]

[5]

[6]

Petroleum consumption (24 types in 5 sectors)

LM tests

(FI)

USA sectors)

U N

1982-2007

Coal consumption

LM (IPS, CBL, DBC, W)

Apergis al.(2010b)

1980-2007

Total Natural gas consumption in billions of cubic feet (per state)

Panel unit roots (various)

Apergis & Payne (2010)

1960-2007

Petroleum consumption

Aslan (2011)

1960-2008

Natural consumption

C A

E C

Pestana Barros et al. (2012)

1981m1- 2010m10

Pestana Barros et al. (2013)

1973m1- 2011m10

gas

I R

Univariate unit roots: Non-stationarity for less than 50% of sectors and products for 9/24 series

0

Multivariate unit roots: Non-stationarity for the commercial and industrial sectors, Stationarity for the residential sector. Evidence for fractorial integration (FI)

USA states)

(50

Stationarity with structural breaks

0

USA states)

(50

Panel unit roots:

0

LM (LS, NP)

USA states)

(50

Stationarity; Endogenously structural breaks

LM (KSS)

USA states)

(50

60% of states follow non-linear behavior

A

D E

T P

(5

[7]

C S

Apergis et al. (2010a)a et

T P

per

Agg/Disagg

M

Nonstationarity

Panel unit roots with endogenously determined structural breaks: Stationarity determined

0 0

For 27/50 states: Nonstationarity For 23/50 states: Stationarity

Renewable consumption

energy

FI

USA

Nonstationarity

0

Nuclear consumption

energy

FI

USA

Nonstationarity

0

[10]

ACCEPTED MANUSCRIPT

Study

Period

Variable

Test

Countries

Result

[1]

[2]

[3]

[4]

[5]

[6]

Gil- Alana et al. (2010)

1973m1- 2009m5

Electricity consumption

Apergis Tsoumas (2012)

&

1989m1-2009m12

Fossil, coal electricity consumption

Apergis Tsoumas (2011)

&

Narayan al.(2008)

et

Borozan Borozan (2015)

&

T P

FI

USA

Nonstationarity

and

FI

USA (different sectors)

Stationarity

1989m1-2009m12

Solar, geothermal and biomass energy consumption

FI

USA (different sectors)

Stationarity

1971-2003

Crude oil production

LM

60 countries

No breaks

2001-2013

Electricity consumption capita

Wang et al. (2016)

1965-2011 (nuclear)

Herrerias al. (2017)

et

1995-2011

Mishra and Smyth (2017)

1973-2014

1990-2011 (RES)

D E

Non-fossil consumption

energy

Electricity, consumption

Coal

T P

E C

C A

LM per

Energy consumption per capita

I R

C S

U N

A

M

LM (LS, IPS)

Breaks

[7] 0 0

0

Nonstationarity

0

Stationarity

Croatia

Nonstationarity

0

Japan

Nonstationarity (nuclear)

0

FLM

Club modeling

Agg/Disagg

Stationarity (RES) convergence

LM and RALS-LM

Chinese regions

Different convergence between rural and urban areas

0

Australia (sector level)

Convergence is confirmed for 6 out of 7 sectors

1

Note: Column [4], ADF: Augmented Dickey-Fuller, LS: Lee and Strazicich 2003-one break: SP:Schmidt and Phillips (1992)-no break, IPS: Im, Pesaran et al. 2003, SURADF: Panel Seemingly unrelated regressions (Breuer et al. 2001), CBL: Carrion-i-Silvestre et al. (2005), KSS: Kapetanios and Shin(2003) and Sollis (2004), F.I. fractorial integration ((Nielsen 2005), (DBC: Del Barrio-Castro et al. 2005; W: Westerlund 2005); NP: Narayan and Popp (2010), FLM: Fourier type Lagrange multiplier (Enders and Lee 2012), S: Sollis (2004) with simultaneous structural change and nonlinearity

[11]

ACCEPTED MANUSCRIPT

1: Aggregate, 0: Disaggregate

T P

I R

C S

A

U N

D E

M

T P

E C

C A

[12]

ACCEPTED MANUSCRIPT Studies either focus on whole economies or group of economies or even sectors of a particular economy such as is done in Ozturk and Aslan (2011) or Aslan and Kum (2011). Typically studies employ tests that have not been used before in a particular country or set of countries and thus they corroborate their results and prove their robustness to make sound policy recommendations. They also demonstrate the difference of results when breaks are inserted (Ozturk and Aslan, 2011) and linear tests vs non-linear unit root tests are employed (Aslan and Kum, 2011). Apparently only four three studies in Table 1, focus on the convergence of clean energy

T

consumption. These are Zhu & Guo (2016), Pestana Barros et al. (2012) and Pestana Barros et al. (2013) for nuclear energy and renewable energy consumption and Apergis & Tsoumas (2011) for solar,

IP

geothermal and biomass energy consumption. Particular focus on these eight emerging economies has

CR

not taken place so far. Thus our study is novel both in the topic and the employed sample. Apparently, none of the studies in Table 1 deals with the investigation of clean energy convergence. Most of the studies deal with total energy or various disaggregate forms of energy

US

consumption. Most of the studies find evidence of convergence which is improved after the employment of structural breaks. Breaks are some crucial time moments after which energy consumption behavior stabilizes towards a certain amount or pattern. These events may be new

AN

policy generation occurrences or economic and political events.

Furthermore, when one observes the results of previous studies on the convergence of renewable energy, or nuclear energy or natural gas or other clean form, then one sees that

M

non-convergence is the main result. This reveals that the penetration of renewable energy sources is still under full development in most economies around the world and that the relevant studies

ED

are still so few that do not allow us to reach a firm conclusion. This justifies the need for the research carried out in the current paper.

PT

Throughout the time span of the current analysis, namely from 1965 to 2016, all the emerging economies have experienced a high average GDP annual growth rate. Brazil has grown at a yearly 4.04%, China at 6.50%, India at 7.17%, Indonesia at 4.43%, Malaysia at 5.89%, Philippines at 4.63%

CE

and Thailand at 5.09% (WDI, 2017). The corresponding figure for high income countries has been 2.76% and for low middle income countries 4.56%. This shows the high GDP growth experienced in the

AC

emerging economies. Particularly Brazil, since 2015 has had a negative growth, but the rest of the sampled countries have pertained a positive economic growth despite the world economic crisis. That been said, it is more than likely that for this economic growth to be maintained, a considerable amount of energy consumption is necessary. This generates greenhouse gas emissions. Thus, emerging economies are facing the issue of energy security (Sadorsky, 2009), given that most of fossil energy proven reserves are found in other countries such as the Middle East, Russia etc. Energy security combined with environmental considerations, brings forward the development of clean energy sources in emerging countries. Two of them, namely China and India are among the top five renewable energy producers. According to Financial Times (2017), emerging economies are close to overtaking the highly developed countries in the renewable energy production with an installed capacity of wind and solar to 307GW and 272GW respectively, fulfilling a 51% and 53% respectively. Evidence for convergence is an issue of huge importance for international environmental organizations since the [13]

ACCEPTED MANUSCRIPT increase of renewable and generally of clean energy is a priority in environmental policy agendas and a crucial result for emerging economies which are to blame for their significant amount of emissions.

3.

Datasets and Descriptive Statistics

To test the stochastic integrational properties of clean energy consumption, we should first define the concept of clean energy. Based on the previous studies of Lee (2013), Paramati et al. (2016) and Cai et

T

al. (2018), the term clean energy is regarded as the energy which does not produce CO 2 emissions. That been said, clean energy is regarded as the sum of the energy consumption from nuclear energy and

IP

renewable energy3. All the data used in this paper has been selected from BP Statistical Review of

CR

World Energy covering the period from 1965 to 2016 expressed in million tonnes of oil equivalent. The countries selected for this paper include Brazil, China, India, Indonesia, Malaysia, Pakistan, Philippines and Thailand, based upon data availability. The descriptive statistics are listed in Table 2. Moreover,

US

these countries have signed and ratified the Kyoto protocol. As reported in Table 2, the maximum of clean energy consumption is in China with 483.5916 million tonnes of oil equivalent. However, the minimum of clean energy consumption is in Malaysia with 4.8633 million tonnes of oil equivalent.

AN

Besides, the standard deviation of China is 112.6459 which is the highest among the countries. For skewness, only Philippines negatively skewed. In terms of the kurtosis, China, India, Malaysia and Thailand is over 3. For the normality test, we find that null is rejected for China, India, Indonesia,

M

Malaysia and Thailand suggesting non-normality of the series.

Econometric Models

PT

4.

ED



Bahmani-Oskooee et al. (2018) propose a quantile unit root test with smooth breaks. We assume a data

CE

generation process y which could be determined by a deterministic trend 𝑑(𝑡) and a stationary error

AC

term with variance δ and zero mean, as follows, y = 𝑑(𝑡) + 𝜀

(1)

To capture the time-varying intercepts with a smooth process, we follow the method proposed by Enders and Lee (2012) and Bahmani-Oskooee et al. (2018) and directly use the terms of 𝛼 sin and 𝛽 cos

, 𝑑(𝑡) = 𝑐 + 𝑎𝑡 + 𝛼 sin

2𝜋𝑘𝑡 2𝜋𝑘𝑡 + 𝛽 cos 𝑇 𝑇

(2)

3 Based on BP Statistical Review of World Energy, the clean energy consumption includes hydroelectricity consumption, nuclear energy consumption, solar energy consumption, wind energy consumption, geothermal, biomass energy consumption and other renewable energy consumption.

[14]

ACCEPTED MANUSCRIPT to be noted, c, k, t and T represent the constant, frequency of the Fourier function, time trend and sample size, respectively. As usual, we set π = 3.1416. Then, the equation (1) can be simplified as, y = 𝑐 + 𝑎𝑡 + 𝛼 sin

2𝜋𝑘𝑡 2𝜋𝑘𝑡 + 𝛽 cos +𝜀 𝑇 𝑇

(3)

where, 𝛼 and 𝛽 measure the amplitude and displacement of the frequency component. In particular, 𝛼 = 𝛽 = 0 is a special case of standard linear specification. Becker et al. (2006) create a more powerful test to detect structural breaks under an unknown form. We set the maximum of 𝐾 = 5 when

T

we determine an optimal 𝑘. For any 𝐾 = 𝑘, we estimate equation (3) by employing the ordinary least

IP

squares (OLS) method and save the sum of squared residuals (SSR). Frequency 𝑘 ∗ is setting as optimum frequency at the minimum of SSR. After searching the optimal frequency 𝑘 ∗ ,

= y − 𝑐̂ − 𝑎𝑡 − 𝛼 sin

2𝜋𝑘𝑡 2𝜋𝑘𝑡 − 𝛽 cos 𝑇 𝑇

(4)

US

y

through the following equation,

CR

Bahmani-Oskooee et al. (2018) obtain the adjusted y

where, 𝑐̂ , 𝑎, 𝛼 and 𝛽 could be estimated through the OLS after searching the optimal frequency 𝑘 ∗ . Next, Bahmani-Oskooee et al. (2018) use the quantile unit root test proposed by Koenker and Xiao

could be presented as follows,

M

y

AN

(2004) to investigate the stationarity of y . The ADF regression model (Diceky and Fuller, 1979) on

ED

y

=𝛼 y

where, 𝑝 is the lag order. In equation (1), 𝛼

+

𝛼 ∆y

+𝜀

(5)

is used to measure the persistency of the data

generation process y . Commonly, if 𝛼 = 1, y

contains a unit root with persistency, and if

PT

|𝛼 | < 1, y is stationary with mean-reverting properties. Based on equation (1), Koenker and Xiao

𝜏y

AC

where, 𝑄

CE

(2004) extend the model to a quantile regression model, which could be expressed as follows,

information set (y

,…,y

,…,y

𝑄

𝜏y

,…,y

= 𝑄 (𝜏) + 𝜃(𝜏)y

+

denotes the 𝜏𝑡ℎ conditional quantile of 𝑦

𝜑 ∆y

(6)

conditional on the

) . 𝑄 (𝜏) is the 𝜏𝑡ℎ conditional quantile of 𝜀 . 𝜃(𝜏) is used to

capture the mean-reverting speed of y at different quantiles. Here, the quantiles is set to be 𝜏 ∈ (0.1, 0.3, 0.5, 0.7, 0.9) . To obtain the coefficient 𝜃(𝜏) and ∑

𝜑 , we could minimize the

following equation,

min



𝜏−𝐼

y

< 𝑄 (𝜏) + 𝜃(𝜏)y

𝜑 ∆y

+

𝜑 ∆y

y

(7)

[15]

− 𝑄 (𝜏) − 𝜃(𝜏)y

ACCEPTED MANUSCRIPT here, I (∙) = 1 if y

< 𝑄 (𝜏) + 𝜃(𝜏)y

+∑

, otherwise I (∙) = 0. Koenker and Xiao

𝜑 ∆y

(2004) further propose t–ratio statistic with the null non-stationary hypothesis α(τ) = 1 against different alternative hypothesis: α(τ) < 1, α(τ) > 1 and α(τ) ≠ 1 to check the unit root hypothesis at specific quantiles, which could be expressed as,

t (𝜏 ) =

𝑓 𝐹

(𝜏 )

𝜏 (1 − 𝜏 )

(𝑌 P

,∆

,…,∆

𝑌 ) 𝜃(𝜏) − 1

(8)

T

where 𝑓(∙) is probability functions of y , and 𝐹(∙) is cumulative density function of series y .

orthogonal to 𝑋 = 1, ∆y

, … , ∆y

. 𝑓 𝐹

) and 𝑃 is the projection matrix onto the space (𝜏 )

IP

is the vector of lagged dependent variables (y

𝑌

is a consistent estimator of 𝑓 𝐹

(𝜏 ) =

(𝜏 − 𝜏 ) 𝐺 (𝜔(𝜏 ) − 𝜔(𝜏

))

(9)

US

𝑓 𝐹

CR

indicated by Koenker and Xiao (2004), which can be expressed as,

(𝜏 )

here 𝜔(𝜏 ) = (𝑐(𝜏 ), 𝜃(𝜏 ), 𝜑 (𝜏 ), … , 𝜑 (𝜏 ) ) and τ ∈ [λ, λ] . We set λ = 0.1 and λ = 0.9 .

AN

Obviously, we test the unit root hypothesis at different quantiles in comparison with traditional ADF test only emphasizing on the conditional central tendency.

To assess the performance of Quantile Unit Root test, Koenker and Xiao (2004) suggested a quantile

𝑄𝐾𝑆 = 𝑆𝑢𝑝

∈[ , ]

|t (𝜏 )|

(10)

ED

M

auto-regression based Kolmogorov-Smirnov (𝑄𝐾𝑆) test which could be presented as,

In this paper, we select the maximum of t (𝜏 ) to build the QKS-Fourier statistics over the

PT

quantiles𝜏 ∈ (0.1, 0.3, 0.5, 0.7, 0.9) . Although the limiting distributions of both 𝑡 (𝜏 ) and 𝑄𝐾𝑆 tests are nonstandard, Koenker and Xiao (2004) propose re-sampling procedures to derive critical

Empirical Findings and Economic Implications

AC

5.

CE

values. In this paper, we make the bootstrap iterations to 10000 times to accurate the empirical results. 4

To avoid the basis of heterogeneity, we first take logarithm to the clean energy consumption of 8 countries before we implement the unit root tests. For comparative purposes, we first use univariate unit root tests including ADF test, PP test, KPSS test and DF-GLS test. For the results of ADF test, the null stationary hypothesis is rejected for Brazil and Philippines with only considering constant, and rejected for Indonesia and Thailand when taking both constant and trend into account. When implementing PP test, the stationary null is rejected for Brazil and Pakistan at 10% significant level with constant only and for Thailand with constant and trend. The KPSS test is rejected for nearly all cases except for Malaysia when considering the constant only. The DF-GLS unit root test is also

4

For detailed bootstrap iterations, please see Koenker and Xiao (2004)

[16]

ACCEPTED MANUSCRIPT utilized to examine the unit root hypothesis. Based on the empirical results, the unit root hypothesis is rejected only for China. No evidence to support stationarity for the rest of the cases.
These conventional unit root tests are notorious for leading to lower testing power if no structural breaks are considered. Perron (1989) has suggested that the unit root tests without considering structural breaks would result to the bias in the empirical results. Given that suggestion, many studies further extend the conventional unit root tests by pouring dummies into

T

the regression model to approximate the sharp breaks leaded by financial crisis, wars and policy

IP

changes. For comparative purposes, we employ the unit root tests with structural breaks including Lumsdaine and Papell (1997); Zivot and Andrews (2002) and Lee and Strazicich (2003).

CR

Table 4 reports the results from the Zivot and Andrews (2002) test with one break. When considering the break in the intercept only, we find the unit root hypothesis is rejected for all countries by determining the optimal frequency by SIC. By taking both intercept and trend into

US

account, we find the null hypothesis cannot by rejected for India, Pakistan and Thailand. For the rest of the countries, we believe the clean energy consumption of the rest of the cases is stationary.

AN

Table 5 reports the results from unit root test of Lumsdaine and Papell (1997). When considering one break, the unit root hypothesis is rejected for Brazil and China. By taking two breaks into account, no evidence supports stationarity of clean energy consumption. Table 6 reports the

M

results from unit root test proposed by Lee and Strazicich (2003). By considering two breaks, the clean energy consumption of China, India and Thailand is stationary. In obvious, the unit root

ED

test of Zivot and Andrews (2002) tends to reject the unit root test. For Lumsdaine and Papell (1997) and Lee and Strazicich (2003), the test is more likely to fail rejections of the unit root

PT

hypothesis.



CE

Since Koenker and Xiao (2004) first proposed the quantile unit root test, the method is widely used to test the stationarity of economic indicators, such as testing on exchange rate including Bahmani-Oskooee and Ranjbar (2016) and Bahmani-Oskooee et al. (2017). The quantile

AC

regression could divide the whole sample into many sub-samples at different selected quantiles. Furthermore, the unit root test built upon the quantile regression could provide specific stationarity at different quantiles. In other words, the shocks to series would make different impacts (permanent and transitory) on the variables to be tested. However, the method proposed by Koenker and Xiao (2004) does not take structural breaks into account. Perron (1989) suggests the failures of conventional unit root tests are always resulted from the failures in detecting the structural breaks. To the best of our knowledge, there are two types of the methods to approximate the structural breaks, i.e., dummy variables and the Fourier function. The dummy variable is more likely to be used to capture the sharp breaks. However, in this paper, we approximate the smooth breaks in the trend through a Fourier function. In fact, after capturing the smooth breaks, we find that the testing power has significantly increased. For comparative purposes, we first utilize the quantile unit root tests without smooth breaks

[17]

ACCEPTED MANUSCRIPT of Koenker and Xiao (2004). The results are listed in Table 7.
The 𝑄𝐾𝑆 statistics are lower than 90% significance level indicating non-stationarity for clean energy consumption for 6 countries, i.e., China, India, Indonesia, Malaysia, Pakistan and Thailand, over the whole quantiles. This means that the shocks on clean energy consumption would make permanent effects on the path for these 6 countries. In other words, through results from the quantile unit root test without smooth breaks, these countries should mandatorily implement energy encouragement policy

T

related to develop clean energy when testing on the overall quantiles. However, the QKS statistic of

IP

Brazil and Philippines is 4.687 and 9.145 indicating the rejection of null non-stationary hypothesis. The energy encouragement policy would only transitorily affect clean energy consumption of Brazil and

CR

Philippines. From the perspective of the unit root behavior at specific quantiles τ = 0.1,0.3,0.5,0.7,0.9, we need to observe the value of t (𝜏 ) with the p-value. Obviously, the clean energy consumption of

US

Brazil, China, Pakistan and Philippines is stationary at upper quantiles. However, for the rest of the countries, we find no evidence to support the stationarity at specific quantiles.

AN



Although the results of quantile unit root tests without smooth breaks incline to non-stationary conclusions for most of the countries, Perron (1989) proves that the ignorance of

M

structural breaks would lead to spurious rejection of the unit root tests. Enders and Lee (2012) propose a version of unit root tests with a Fourier function to approximate the smooth breaks. Moreover,

ED

Bahmani-Oskooee et al. (2018) propose a quantile unit root test with smooth breaks which could efficiently approximate the unknown structural breaks in the series. To implement the Fourier quantile unit root test, we should first specify some parameters in the model. As mentioned above, the optimal

PT

frequency in equation (2) is selected by an F test and the lag is selected by AIC criteria. Table 8 reports the results of the Fourier quantile unit root test. We first focus on the Fourier QKS statistics with the null of unit root over the range of quantiles, namely [0.1,0.3,0.5,0.7,0.9]. The unit root hypothesis is

CE

rejected for 3 out of 8 countries, i.e., China, Pakistan and Thailand. That is, after approximating the structural breaks in the deterministic trend, we find that more are countries in favor of the stationary

AC

hypothesis than the quantile unit root test without smooth breaks proposed by Koenker and Xiao (2004). As noted by Perron (1989), structural breaks in the series would significantly lower the testing power of the unit root tests than without considering structural breaks. In other words, the shocks to clean energy consumption for China, Pakistan and Thailand would make transitory effects on the path. The encouragement policy for the clean energy would transitorily affect the path. In other words, no matter what the shock (wars, regulatory policies, financial crisis and oil crisis) to the series is, the clean energy consumption would finally return to the mean due to the mean-reverting properties. For the rest of the countries, we find no supports to the divergence for the clean energy consumption. The energy encouragement policy for those countries would only make permanent impacts on the path. The policy authorities in these countries should implement related policies to develop clean energy to achieve green growth.

[18]

ACCEPTED MANUSCRIPT To further examine the unit root behavior in each quantile of the clean energy consumption for selected countries, we should see the results of t (𝜏 ) with the p-values. Some interesting findings could be summed up. First, the t (𝜏 ) of the clean energy consumption for the countries like China, India and Philippines is significant at 10% level at the upper quantiles (τ = 0.7,0.9). The shocks to clean energy consumption of these countries would permanently affect the path at upper quantiles for countries like China, India and Philippines. Second, the t (𝜏 ) of the clean energy consumption for Indonesia is only significant at the extreme and lowest quantile τ = 0.1. Shocks to clean energy consumption would permanently affect the path in the long-run at the lower quantiles of the series.

T

Third, the t (𝜏 ) of the clean energy consumption for Thailand is nearly significant at overall

IP

quantiles except extremely the low quantile (τ = 0.1). Except for the lower quantiles, the shocks would affect the path in a transitory way. Fourth, we find unit root behavior for Brazil, Philippines and

CR

Malaysia at all quantiles τ = 0.1,0.3,0.5,0.7,0.9. In other words, the unit root performance shows asymmetric properties at different quantiles even when we focus on a given series. These findings are

US

the first time to be revealed in the field of testing the convergence of clean energy consumption. Besides, we also calculate the half-life of a shock for these emerging market countries. Based on the half-life for specific countries, we could find that even the half-lives for a country are

AN

different at various quantiles too. The shocks to the clean energy consumption would affect the mean-reverting process at a different magnitude. The half-lives of a shock to a given series are also

M

varied over the different quantiles, which again reconfirm the asymmetric properties in half-lives. Finally, we plot the path of clean energy consumption with fitted intercepts with and without

ED

smooth breaks in Figure 1. The clean energy consumption is colored in blue. The fitted intercepts without smooth breaks are colored in black. The fitted intercepts with smooth breaks are colored in red. Obviously, the fitted intercepts with the Fourier function could be more fitted to the path of clean

PT

energy consumption, especially for the countries of Brazil, Indonesia, Pakistan and Philippines. The path of clean energy consumption for these countries appears to be more cupped than in the rest of the countries. The Fourier function in the intercepts could efficiently approximate the overall trend. In

CE

other words, the estimations from quantile unit root test with smooth breaks are more reliable. For the countries like China, India and Thailand, the fitted intercepts, with and without smooth breaks, are

AC

closer due to their less concavity. In other words, the Fourier function poured into the intercepts would automatically fit the trend of clean energy consumption for these countries. After comparing the fitted intercepts with and without smooth breaks, we support that the Fourier function would efficiently approximate the structural breaks in the series. Therefore, future studies may consider using Fourier function to test the structural breaks in energy consumption. Overall, the asymmetric performances of the integrational properties of clean energy consumption are first been suggested through the quantile unit root test with and without smooth breaks. This means that each one of them is following a different energy consumption path and a different environmental strategy and that the production of clean energies is still under increasing evolution in these countries. Besides, even within the same series, the integrational properties of clean energy consumption show asymmetric properties too. Furthermore, the asymmetric half-lives also suggest the

[19]

ACCEPTED MANUSCRIPT necessities to introduce the quantile regression into testing on the integrational properties of energy consumption. At this point, we need to underline the fact that these countries do not belong to a broader coalition or economic union as for example the European Union does. Thus, the development of clean energies is relied on the countries’ individual planning and their general ratification of Kyoto protocol demands and the agreements that have followed, but it is not nevertheless dictated by obligatory directives and guidelines, as is the case within an official economic union. Up-to-date the results produced in this paper have not been produced elsewhere. The few studies on renewable energy in general, or some types of renewable energy, are not about

T

emerging economies but for the USA. Therefore, we cannot compare the results of our study with

IP

the results of previous studies, because they are incomparable cases. This fact shows the

Conclusion

US

6.

CR

originality of our research and the need for further research in this field.

The current paper has used an innovative method of unit root testing for the investigation of integrational properties of clean energy consumption in eight emerging economies. The reason we have

AN

focused on these countries is that they have experienced high GDP growth so far, they are still under development with high industrialization and urbanization and they are bound to consume even more energy in the future than they currently do. Consequently, they are major greenhouse gas emitters. At

M

the same time, some of those emerging countries are at the forefront of clean energy production, e.g.

ED

China ranks among the five largest producers of renewable energy worldwide. The innovative method employed in the paper uses a Fourier function approach in the unit root testing, also in different quantiles of the distribution of the investigated variable. This approach

PT

provides individual information for every quantile, it tolerates nonlinearities and provides adequate smoothing for structural breaks.

CE

Based on the results specific to this paper, the clean energy consumption of 2 countries, i.e., Brazil and Philippines, is stationary. However, after approximating the smooth breaks, we find that 3 of 8 countries, i.e., China, Pakistan and Thailand, are stationary. Besides, based on the integrational

AC

properties at each quantile, we find asymmetric performance of stationarity of clean energy consumption for a given country. The estimations from quantile unit root test with smooth breaks are more reliable because the time-varying intercepts with smooth breaks are more fitted to the path. This has various interpretations. First the emerging economies follow a different energy mix from each other and do not conform to common energy strategies, since they are heterogeneous countries both geographically and politically. Thus, some countries are leaders in clean energy production and others are not. Besides, the integrational properties at each quantile are also provided by quantile unit root test. Specifically, energy policy regulators should realize the asymmetric behavior of the integrational properties at different quantiles. Although the clean energy consumption of India and Indonesia is stationary when testing on the whole quantile, the stationarities are asymmetric at specific quantiles. [20]

ACCEPTED MANUSCRIPT Regulators of these two countries should pay attention to these asymmetric effects. Moreover, the clean energy consumption for countries like Brazil and Malaysia and Philippines is non-stationary even from the perspective of quantile-level. Lastly, the clean energy consumption of China, Pakistan and Thailand is not stationary not only at the overall quantiles, but also at some specific quantiles. Integrating properties of clean energy consumption would entail that countries are resilient to shocks and that the latter will not leave permanent effects. Convergence also means that countries have reached a stable point of clean energy production and that there is a kind of dynamic equilibrium at that

T

point. Conversely, non-convergence entails that countries are vulnerable to shocks and thus more support should be provided by governments and energy authorities for the establishment and the

IP

promotion of clean energies. With respect to policy implications, our results clearly imply that the

CR

green energy encouragement policy would have a transitory effect for countries which exhibit non- convergence. Conversely, the policy would have a permanent effect for countries with evidence of convergence. Thus, emerging economies in pursuit of cleaner production through the

US

employment of clean energies will manage to have clean energy penetration in cases of convergence.

AN

References

M

Aslan, A. (2011). Does natural gas consumption follow a nonlinear path over time? Evidence from 50 US States. [Article]. Renewable and Sustainable Energy Reviews, 15(9), 4466-4469, doi:10.1016/j.rser.2011.07.105.

ED

Aslan, A., & Kum, H. (2011). The stationary of energy consumption for Turkish disaggregate data by employing linear and nonlinear unit root tests. [Article]. Energy, 36(7), 4256-4258, doi:10.1016/j.energy.2011.04.018.

PT

Apergis, N., & Payne, J. E. (2010). Structural breaks and petroleum consumption in US states: Are shocks transitory or permanent? [Article]. Energy Policy, 38(10), 6375-6378, doi:10.1016/j.enpol.2010.06.015.

CE

Apergis, N., Loomis, D., & Payne, J. E. (2010a). Are fluctuations in coal consumption transitory or permanent? Evidence from a panel of US states. [Article]. Applied Energy, 87(7), 2424-2426, doi:10.1016/j.apenergy.2009.12.005.

AC

Apergis, N., Loomis, D., & Payne, J. E. (2010b). Are shocks to natural gas consumption temporary or permanent? Evidence from a panel of U.S. states. [Article]. Energy Policy, 38(8), 4734-4736, doi:10.1016/j.enpol.2010.03.016. Apergis, N., & Tsoumas, C. (2012). Long memory and disaggregated energy consumption: Evidence from fossils, coal and electricity retail in the U.S. [Article]. Energy Economics, 34(4), 1082-1087, doi:10.1016/j.eneco.2011.09.002. Apergis, N., & Tsoumas, C. (2011). Integration properties of disaggregated solar, geothermal and biomass energy consumption in the U.S. [Article]. Energy Policy, 39(9), 5474-5479, doi:10.1016/j.enpol.2011.05.015. Bahmani-Oskooee, M., & Ranjbar, O. (2016). Quantile unit root test and PPP: evidence from 23 OECD countries. Applied Economics, 48(31), 2899-2911. Bahmani-Oskooee, M., Chang, T., Chen, T. H., & Tzeng, H. W. (2017). Quantile unit root test and the PPP in Africa. Applied Economics, 49(19), 1913-1921. Bahmani-Oskooee, M., Chang, T., Elmi, Z., & Ranjbar, O. (2018). Re-testing Prebisch–Singer hypothesis: new evidence using Fourier quantile unit root test. Applied Economics, 50(4), 441-454.

[21]

ACCEPTED MANUSCRIPT Bai, J., & Perron, P. (2003). Computation and analysis of multiple structural change models. Journal of applied econometrics, 18(1), 1-22. Becker, R., Enders, W., & Lee, J. (2006). A stationarity test in the presence of an unknown number of smooth breaks. Journal of Time Series Analysis, 27(3), 381-409. Borozan, D., & Borozan, L. The stationarity of per capita electricity consumption in Croatia allowing for Structural Break(s). In Proceedings of the 13th International Symposium on Operational Research, SOR 2015, 2015 (pp. 337-342) Cai, Y., Sam, C. Y., & Chang, T. (2018). Nexus between clean energy consumption, economic growth and CO2 emissions. Journal of Cleaner Production, 182, 1001-1011.

IP

T

Chen, P. F., & Lee, C. C. (2007). Is energy consumption per capita broken stationary? New evidence from regional-based panels. [Article]. Energy Policy, 35(6), 3526-3540, doi:10.1016/j.enpol.2006.12.027.

CR

Christopoulos, D. K., & León-Ledesma, M. A. (2010). Smooth breaks and non-linear mean reversion: Post-Bretton Woods real exchange rates. Journal of International Money and Finance, 29(6), 1076-1093. Enders, W., & Lee, J. (2012). A unit root test using a Fourier series to approximate smooth breaks. Oxford Bulletin of Economics and Statistics, 74(4), 574-599.

US

Financial Times (2017).Emerging markets poised to lead pack on renewable energy, available from: [https://www.ft.com/content/b0a5f72a-9fa4-11e7-9a86-4d5a475ba4c5], accessed on 12/04/2018

AN

Gil-Alana, L. A., Loomis, D., & Payne, J. E. (2010). Does energy consumption by the US electric power sector exhibit long memory behavior? [Article]. Energy Policy, 38(11), 7512-7518, doi:10.1016/j.enpol.2010.07.018.

M

Goldemberg, J. (1996). Communication: A note on the energy intensity of developing countries. [Article]. Energy Policy, 24(8), 759-761, doi:10.1016/0301-4215(96)00045-6.

ED

Gozgor, G. (2016). Are shocks to renewable energy consumption permanent or transitory? An empirical investigation for Brazil, China, and India. Renewable and Sustainable Energy Reviews, 66, 913-919.

PT

Hasanov, M., & Telatar, E. (2011). A re-examination of stationarity of energy consumption: Evidence from new unit root tests. [Article]. Energy Policy, 39(12), 7726-7738, doi:10.1016/j.enpol.2011.09.017.

CE

Herrerias, M. J., Aller, C., & Ordóñez, J. (2017). Residential energy consumption: A convergence analysis across Chinese regions. [Article]. Energy Economics, 62, 371-381, doi:10.1016/j.eneco.2016.06.006. Hsu, Y. C., & Lee, C. C. (2008). Revisited: Are shocks to energy consumption permanent or temporary? New evidence from a panel SURADF approach. [Article]. Energy Economics, 30(5), 2314-2330, doi:10.1016/j.eneco.2007.09.007.

AC

Koenker, R., & Xiao, Z. (2004). Unit root quantile autoregression inference. Journal of the American Statistical Association, 99(467), 775-787. Kula, F. (2014). Is per capita electricity consumption non-stationary? A long-span study for Turkey. [Article]. Energy Sources, Part B: Economics, Planning and Policy, 9(2), 161-164, doi:10.1080/15567249.2010.483447. Kula, F., Aslan, A., & Ozturk, I. (2012). Is per capita electricity consumption stationary? Time series evidence from OECD countries. [Article]. Renewable and Sustainable Energy Reviews, 16(1), 501-503, doi:10.1016/j.rser.2011.08.015. Kum, H. (2012). Are fluctuations in energy consumption transitory or permanent? Evidence from a panel of East Asia & Pacific countries. [Article]. International Journal of Energy Economics and Policy, 2(3), 92-96. Lean, H. H., & Smyth, R. (2009). Long memory in US disaggregated petroleum consumption: Evidence from univariate and multivariate LM tests for fractional integration. [Article]. Energy Policy, 37(8), 3205-3211, doi:10.1016/j.enpol.2009.04.017. [22]

ACCEPTED MANUSCRIPT Lean, H. H., & Smyth, R. (2014). Are shocks to disaggregated energy consumption in Malaysia permanent or temporary? Evidence from LM unit root tests with structural breaks. [Article]. Renewable and Sustainable Energy Reviews, 31, 319-328, doi:10.1016/j.rser.2013.10.040. Lee, C. C., & Chang, C. P. (2008). New evidence on the convergence of per capita carbon dioxide emissions from panel seemingly unrelated regressions augmented Dickey–Fuller tests. Energy, 33(9), 1468-1475. Lee, J. W. (2013). The contribution of foreign direct investment to clean energy use, carbon emissions and economic growth. Energy Policy, 55, 483-489. Lee, J., & Strazicich, M. C. (2003). Minimum Lagrange multiplier unit root test with two structural breaks. Review of economics and statistics, 85(4), 1082-1089.

IP

T

Magazzino, C., (2014), A Panel VAR Approach of the Relationship Among Economic Growth, CO2 Emissions, and Energy Use in the ASEAN-6 Countries, International Journal of Energy Economics and Policy, 4, 4, 546-553

CR

Magazzino, C., (2016), CO2 emissions, economic growth and energy use in the Middle East countries: A panel VAR approach, Energy Sources, Part B: Economics, Planning, and Policy 11, 10, 960-968

US

Magazzino, C., (2017), The relationship among economic growth, CO2 emissions, and energy use in the APEC countries: a panel VAR approach, Environment Systems and Decisions, 37, 3, 353-366

AN

Menegaki, A. N. and S. Tsani (2018). Chapter 5 - Critical Issues to Be Answered in the Energy-Growth Nexus (EGN) Research Field. The Economics and Econometrics of the Energy-Growth Nexus, Academic Press: 141-184.

M

Tsani, S. and A. N. Menegaki (2018). Chapter 11 - The Energy-Growth Nexus (EGN) Checklist for Authors. The Economics and Econometrics of the Energy-Growth Nexus, Academic Press: 347-376. Liddle, B. (2009). Electricity intensity convergence in IEA/OECD countries: Aggregate and sectoral analysis. Energy Policy, 37(4), 1470-1478.

ED

Lumsdaine, R. L., & Papell, D. H. (1997). Multiple trend breaks and the unit-root hypothesis. Review of economics and Statistics, 79(2), 212-218.

PT

Meng, M., Payne, J. E., & Lee, J. (2013). Convergence in per capita energy use among OECD countries. Energy Economics, 36, 536-545.

CE

Mishra, V., Sharma, S., & Smyth, R. (2009). Are fluctuations in energy consumption per capita transitory? Evidence from a panel of Pacific Island countries. [Article]. Energy Policy, 37(6), 2318-2326, doi:10.1016/j.enpol.2009.02.022. Mishra, V., & Smyth, R. (2017). Conditional convergence in Australia's energy consumption at the sector level. [Article]. Energy Economics, 62, 396-403, doi:10.1016/j.eneco.2016.07.009.

AC

Narayan, P. K., & Smyth, R. (2005). The residential demand for electricity in Australia: An application of the bounds testing approach to cointegration. [Article]. Energy Policy, 33(4), 467-474, doi:10.1016/j.enpol.2003.08.011. Narayan, P. K., & Smyth, R. (2007). Are shocks to energy consumption permanent or temporary? Evidence from 182 countries. Energy policy, 35(1), 333-341. Narayan, P. K., Narayan, S., & Smyth, R. (2008). Are oil shocks permanent or temporary? Panel data evidence from crude oil and NGL production in 60 countries. [Article]. Energy Economics, 30(3), 919-936, doi:10.1016/j.eneco.2007.07.005. Nilsson, L. J. (1993). Energy intensity trends in 31 industrial and developing countries 1950-1988. [Article]. Energy, 18(4), 309-322, doi:10.1016/0360-5442(93)90066-m. Sadorsky, P. (2009). "Renewable energy consumption and income in emerging economies." Energy Policy 37(10): 4021-4028.

[23]

ACCEPTED MANUSCRIPT Smyth, R. (2013). Are fluctuations in energy variables permanent or transitory? A survey of the literature on the integration properties of energy consumption and production. [Review]. Applied Energy, 104, 371-378, doi:10.1016/j.apenergy.2012.10.069. Ozcan, B. (2013). Are shocks to energy consumption permanent or temporary? the case of 17 middle east countries. [Article]. Energy Exploration and Exploitation, 31(4), 589-605, doi:10.1260/0144-5987.31.4.589. Ozturk, I., & Aslan, A. (2011). Are fluctuations in energy consumption per capita transitory? Evidence from Turkey. [Article]. Energy Exploration and Exploitation, 29(2), 161-167, doi:10.1260/0144-5987.29.2.161.

T

Payne, J. E., Miller, S., Lee, J., & Cho, M. H. (2014). Convergence of per capita sulphur dioxide emissions across US states. Applied Economics, 46(11), 1202-1211.

IP

Payne, J. E., Vizek, M., & Lee, J. (2017). Stochastic convergence in per capita fossil fuel consumption in US states. Energy Economics, 62, 382-395.

CR

Paramati, S. R., Ummalla, M., & Apergis, N. (2016). The effect of foreign direct investment and stock market growth on clean energy use across a panel of emerging market economies. Energy Economics, 56, 29-41.

US

Perron, P. (1989). The great crash, the oil price shock, and the unit root hypothesis. Econometrica: Journal of the Econometric Society, 1361-1401.

AN

Pestana Barros, C., Gil-Alana, L. A., & Payne, J. E. (2012). Evidence of long memory behavior in U.S. renewable energy consumption. [Article]. Energy Policy, 41, 822-826, doi:10.1016/j.enpol.2011.11.052. Pestana Barros, C., L. A. Gil-Alana, et al. (2012). "Evidence of long memory behavior in U.S. renewable energy consumption." Energy Policy 41: 822-826.

M

Solarin, S.A., Gil-Alana, L.A. and Al-Mulali, U. (2018). Stochastic convergence of renewable energy consumption in OECD countries: a fractional integration approach, Environmental Science and Pollution Research 25: 17289-17299.

ED

Wang, Y., Li, L., Kubota, J., Zhu, X., & Lu, G. (2016). Are fluctuations in Japan's consumption of non-fossil energy permanent or transitory? [Article]. Applied Energy, 169, 187-196, doi:10.1016/j.apenergy.2016.02.009. Indicators

(WDI),

available

from:

PT

Worldbank (2017). World Development [https://data.worldbank.org/products/wdi]

CE

Wang, Y., Li, L., Kubota, J., Zhu, X., & Lu, G. (2016). Are fluctuations in Japan's consumption of non-fossil energy permanent or transitory? [Article]. Applied Energy, 169, 187-196, doi:10.1016/j.apenergy.2016.02.009. Worldbank (2017). World Development [https://data.worldbank.org/products/wdi]

Indicators

(WDI),

available

from:

AC

Yilanci, V., & Tunali, Ç. B. (2014). Are fluctuations in energy consumption transitory or permanent? Evidence from a Fourier LM unit root test. [Article]. Renewable and Sustainable Energy Reviews, 36, 20-25, doi:10.1016/j.rser.2014.04.002. Zhu, H., & Guo, P. (2016). Are shocks to nuclear energy consumption per capita permanent or temporary? A global perspective. [Article]. Progress in Nuclear Energy, 88, 156-164, doi:10.1016/j.pnucene.2015.12.013. Zivot, E., & Andrews, D. W. K. (2002). Further evidence on the great crash, the oil-price shock, and the unit-root hypothesis. Journal of business & economic statistics, 20(1), 25-44.

[24]

ACCEPTED MANUSCRIPT Table 2 Descriptive Statistics Mean

Median

Maximum

Minimum

Std. Dev.

Kurtosis

Jarque-Bera

Probability

BRAZIL

54.7381

50.2202

128.5310

5.4250

36.1708

Skewness 0.3526

2.0204

3.1569

0.2063

CHINA

78.4559

28.4789

483.5916

4.3860

112.6459

2.0914

6.6742

67.1567***

0.0000

INDIA

16.9489

70.6508

4.3370

17.2319

1.3737

3.8865

18.0560***

0.0001

2.8903

1.7769

8.3836

0.2400

2.7654

0.7171

2.0530

6.3998**

0.0408

MALAYSIA

1.1737

1.0068

4.8633

0.1350

1.0379

1.6107

5.6543

37.7482***

0.0000

PAKISTAN

4.1388

4.0805

9.8448

0.4970

2.6699

0.2889

1.9606

3.0641

0.2161

PHILIPPINES

3.8248

3.9152

8.2793

0.1590

2.6295

-0.1108

1.5993

4.3571

0.1132

THAILAND

1.5945

0.9967

6.4360

0.2480

1.4386

1.6502

5.1963

34.0531***

0.0000

AC

CE

PT

ED

M

AN

US

CR

IP

Note: ***, ** and * denote 1%, 5% and 10% significant levels.

T

22.2428

INDONESIA

[25]

ACCEPTED MANUSCRIPT

Table 3 Univariate Unit Root Test ADF(k) Constant BRAZIL CHINA INDIA INDONESIA MALAYSIA PAKISTAN PHILIPPINES THAILAND

-4.8471(0)*** 2.3843(0) -0.2574(0) -0.1788(1) -0.6154(0) -2.0889(0) -3.0096(7)** -0.4099(1)

Constant and Trend -1.7598(0) -0.8580(0) -1.7900(0) -4.3939(5)*** -2.2645(0) -1.8452(0) -1.6641(7) -4.5138(0)***

PP(k) Constant

Constant and Trend

KPSS(k) Constant

Constant and Trend

-4.6848(1)*** 2.5057(3) -0.2692(1) -0.4664(0) -0.5863(6) -2.8641(6)* -2.1156(2) -0.3761(5)

-1.7400(1) -0.8915(2) -1.9150(2) -1.9835(1) -2.4009(2) -1.6889(2) -1.1854(2) -4.5379(2)***

0.9147(5)*** 0.9624(5)*** 0.9517(5)*** 0.8789(5)*** 0.9367(5)*** 0.9223(5)*** 0.8473(5)*** 0.9635(5)***

0.2290(5)*** 0.2024(5)** 0.1338(5)* 0.1244(5)* 0.1020(5) 0.2352(5)*** 0.2245(5)*** 0.1371(3)*

U N

C S

I R

T P

DF-GLS(k) Constant Constant and Trend 0.2707(7) -0.8836 2.0879(1)** -0.7704 2.3301(0) -1.7120 0.1819(5) -2.3332 1.1138(0) -2.2078 0.9402(0) -1.3783 -0.1437(9) -1.4372 1.5463(4) -4.4233

Note: ***, ** and * denote 1%, 5% and 10% significant levels. K in the bracket of ADF test is determined by SIC with the maximum 10 lags. The k in the brackets of PP test and KPSS test are Newey-West automatic bandwidth using Bartlett Kernel.

A

D E

M

T P

E C

C A

[26]

ACCEPTED MANUSCRIPT

Table 4 Zivot-Andrews Unit Root test

BRAZIL CHINA INDIA INDONESIA MALAYSIA PAKISTAN PHILIPPINES THAILAND

Intercept p-value 0.0177 0.0234 0.0023 0.0043 0.0112 0.0099 0.0022 0.0174

t-statistic(K) -2.4590(0)** -2.8566(0)** -3.3995(0)** -6.0769(10)*** -3.1404(0)** -3.2809(7)*** -6.5653(7)*** -5.3829(0)**

breaking point 1979 2008 1982 1998 1996 1997 1980 1978

t-statistic -4.4605(1)*** -4.4229(0)*** -3.4319(0) -8.1589(10)*** -2.9829(0)** -4.5327(0) -6.4176(0)*** -5.7791(0)

A

T P

I R

C S

U N

Intercept and trend p-value 0.0000 0.0003 0.2176 0.0014 0.0198 0.4557 0.0000 0.1217

Note: ***, ** and * denote 1%, 5% and 10% significant levels. K in the bracket of ADF test is determined by SIC with the maximum 10 lags.

D E

M

T P

E C

C A

[27]

breaking point 1973 2003 2000 1998 2001 1987 1980 1998

ACCEPTED MANUSCRIPT

Table 5 Lumsdaine and Papell Unit Root test

BRAZIL CHINA INDIA INDONESIA MALAYSIA PAKISTAN PHILIPPINES THAILAND

One break t-statistic(K) -4.6695(2)** -4.9469(2)*** -4.0717(2) -2.4358(2) -2.1783(2) -3.7726(2) -4.8281(2)*** -3.9891(2)

T P

Two breaks breaking date 1979 2003 2002 1975 2009 1987 1984 2001

t-statistic(K) -5.6231(2) -4.9533(2) -4.9365(2) -5.5250(2) -3.0398(2) -3.9032(2) -5.8997(2) -4.1802(2)

C S

I R

A

U N

1978 1983 1978 1980 1988 1988 1977 1975

breaking date 1985 2003 2002 1992 2007 2001 1984 1999

Note: ***, ** and * denote 1%, 5% and 10% significant levels. K in the bracket of ADF test is determined by SIC with the maximum 5 lags. The critical values at 1%, 5% and 10% level are -4.93, -4.42 and -4.11 when only considering one break in the trend and -7.19, -6.62 and -6.37 when considering two breaks in the trend, respectively.

D E

M

T P

E C

C A

[28]

ACCEPTED MANUSCRIPT

Table 6 Lee and Strazicich Unit Root test

BRAZIL CHINA INDIA INDONESIA MALAYSIA PAKISTAN PHILIPPINES THAILAND

1% C.V. -6.6910 -6.8210 -6.6910 -6.9320 -6.6910 -6.9780 -7.0320 -7.1960

t-statistic(K) -5.4200(4) -10.2518(2)*** -5.8825(1)* -5.5416(0) -5.5011(5) -5.7476(4) -5.3259(5) -6.0646(5)*

5% C.V. -6.1520 -6.1660 -6.1520 -6.1750 -6.1520 -6.2880 -6.3750 -6.3120

10% C.V. -5.7980 -5.8320 -5.7980 -5.8250 -5.7980 -5.9980 -6.0110 -5.8930

I R

C S

U N

T P

breaking points

1974 1981 1980 1982 1993 1985 1978 1978

Note: ***, ** and * denote 1%, 5% and 10% significant levels. K in the bracket of ADF test is determined by SIC with the maximum 5 lags.

A

D E

M

T P

E C

C A

[29]

1985 2003 2005 1992 1998 1997 1989 2005

ACCEPTED MANUSCRIPT

Table 7 Results for quantile unit root test without smooth breaks Brazil

CHINA

INDIA

INDONESIA

MALAYSIA

PAKISTAN

PHILIPPINES

THAILAND

θ(τ ) t (τ ) QKS θ(τ ) t (τ ) QKS θ(τ ) t (τ ) QKS θ(τ ) t (τ ) QKS θ(τ ) t (τ ) QKS θ(τ ) t (τ ) QKS θ(τ ) t (τ ) QKS θ(τ ) t (τ ) QKS

0.1 0.973(0.172) -0.935(0.474) 4.687** 1.016(0.350) 1.052(0.922) 2.524 1.019(0.391) 0.614(0.912) 1.040 1.103(0.293) 0.603(0.977) 1.921 0.950(0.253) -1.257(0.349) 1.509 0.994(0.480) -0.072(0.827) 3.598 0.971(0.294) -0.797(0.446) 9.145*** 1.128(0.293) 1.373(0.986) 1.778

D E

T P

E C

0.3 0.949(0.006)*** -3.470(0.005)*** 3.871(90% C.V.) 1.029(0.186) 1.980(0.994) 4.077(90% C.V. 1.023(0.347) 1.040(0.980) 4.035(90% C.V. 0.976(0.342) -0.638(0.476) 2.997(90% C.V. 0.958(0.295) -1.053(0.487) 3.304(90% C.V.) 0.967(0.327) -0.908(0.488) 3.976(90% C.V. 0.945(0.021)** -2.536(0.034)** 3.709(90% C.V. 1.069(0.326) 0.735(0.985) 3.624(90% C.V.

0.5 0.957(0.003)*** -3.074(0.009)*** 4.670(95% C.V.) 1.014(0.329) 0.919(0.967) 4.855(95% C.V.) 0.983(0.385) -0.617(0.676) 4.787(95% C.V.) 0.983(0.377) -0.676(0.527) 3.460(95% C.V.) 1.020(0.404) 0.571(0.965) 3.793(95% C.V.) 0.964(0.265) -1.771(0.162) 4.687(95% C.V.) 0.943(0.015)** -2.392(0.044)** 4.437(95% C.V.) 1.008(0.477) 0.145(0.923) 4.143(95% C.V.)

U N

M

A

Note: ***, ** and * denote 1%, 5% and 10% significant level. Bootstrap procedures with 20000 replications are used to generate the p-value.

C A

[30]

T P

I R

C S

0.7 0.946(0.000)*** -3.402(0.005)*** 7.127(99%.C.V.) 1.011(0.397) 0.762(0.935) 7.116(99%.C.V.) 1.002(0.489) 0.077(0.902) 6.918(99%.C.V.) 0.955(0.238) -1.405(0.184) 4.610(99%.C.V.) 0.973(0.354) -0.840(0.501) 5.008(99%.C.V.) 0.953(0.200) -2.333(0.065)* 6.943(99%.C.V.) 0.945(0.052)* -1.978(0.062)* 6.654(99%.C.V.) 0.980(0.440) -0.444(0.750) 5.652(99%.C.V.)

0.9 0.917(0.006)*** -4.687(0.021)** 1.026(0.299) 2.524(0.987)*** 1.021(0.399) 0.650(0.918) 0.889(0.158) -1.921(0.110) 0.999(0.495) -0.014(0.828) 0.938(0.207) -3.598(0.053)* 0.848(0.010)*** -9.145(0.001)*** 0.934(0.313) -1.778(0.251)

ACCEPTED MANUSCRIPT

Table 8 Results for Quantile Unit Root Test with Smooth Breaks

BRAZIL

CHINA

INDIA

INDONESIA

Frequency

F stat.

0.1

1653.677

0.1

0.1

0.7

2872.182

569.403

450.752

0.1

0.3

0.5

θ(τ )

0.9290(0.330)

1.0400(0.3540)

0.9870(0.4430)

PAKISTAN

PHILIPPINES

0.1

0.1

0.1

271.47

1289.429

543.883

-0.699(0.624)

-0.417(0.435)

14.3985

7.8030

4.639 (99% C.V.) 0.6510(0.0020)***

0.6880(0.0120)**

-1.599(0.329)

-3.291(0.006)***

-2.773(0.098)*

2.8484

1.6148

1.8535

4.187 (95% C.V.) 0.8090(0.0740)*

5.745 (99% C.V.) 0.7490(0.0010)***

0.7540(0.0530)**

t (τ )

-0.424(0.617)

0.434(0.954)

-0.205(0.778)

9.4118

NA

52.9717

Fourier QKS θ(τ )

0.765 0.7890(0.0860)*

3.12 (90% C.V.) 0.9080(0.2790)

3.511(95% C.V.) 0.7840(0.0700)*

t (τ )

-1.556(0.3400)

-0.592(0.707)

half-life

2.9248

7.1821

Fourier QKS θ(τ )

3.876* 0.9100(0.2420)

3.624 (90% C.V.) 0.8750(0.1390)

t (τ )

-0.718(0.311)

-0.887(0.569)

-1.931(0.234)

-2.839(0.017)**

-1.599(0.210)

half-life

7.3496

5.1909

3.2702

2.3983

2.4548

Fourier QKS θ(τ )

2.839 0.2880(0.000)***

3.159 (90% C.V.) 0.7140(0.0100)***

3.498 (95% C.V.) 0.8000(0.0400)**

4.449 (99% C.V.) 0.8320(0.1230)

0.8740(0.2520)

t (τ )

-2.334(0.050)**

-1.417(0.194)

-1.608(0.167)

-0.978(0.394)

-0.789(0.475)

half-life

0.5568

2.0576

3.1063

3.7687

5.1468

2.334 0.7750(0.0540)*

3.064 (90% C.V.) 0.8430(0.0840)*

3.486 (95% C.V.) 0.8460(0.0990)*

4.613 (99% C.V.) 0.8850(0.1760)

0.8920(0.2970)

t (τ )

-1.499(0.1980)

-1.449(0.268)

-1.32(0.424)

-1.021(0.398)

-0.321(0.608)

half-life

2.7194

4.0585

4.1447

5.6737

6.0649

Fourier QKS θ(τ )

1.597 0.5620(0.1100)

3.121 (90% C.V.) 0.7580(0.0210)**

3.521 (95% C.V.) 0.6620(0.000)***

4.544 (99% C.V.) 0.6060(0.000)***

0.4570(0.001)***

D E

T P

θ(τ )

E C

C A

T P

0.9 0.9150(0.2680)

half-life

Fourier QKS MALAYSIA

0.7 0.9530(0.3050)

M

C S

U N

A

I R

t (τ )

-0.844(0.526)

-1.018(0.296)

-3.114(0.012)**

-3.423(0.004)***

-2.417(0.043)**

half-life

1.2029

2.5017

1.6804

1.3839

0.8852

Fourier QKS θ(τ )

3.745** 0.7590(0.107)

3.107 (90% C.V.) 0.9510(0.320)

3.515 (95% C.V.) 0.8100(0.049)**

4.559 (99% C.V.) 0.7990(0.075)**

0.5780(0.042)**

t (τ )

-0.781(0.434)

-0.476(0.569)

-1.737(0.217)

-1.565(0.131)

-2.179(0.172)

half-life

2.5136

13.7964

3.2894

3.0890

1.2644

[31]

ACCEPTED MANUSCRIPT

Frequency

F stat.

0.1

0.3

0.5

0.7

0.9

θ(τ )

2.976 0.4960(0.063)*

3.627 (90% C.V.) 0.4350(0.000)***

4.325 (95% C.V.) 0.2900(0.000)***

6.186 (99% C.V.) 0.1950(0.000)***

0.3210(0.000)***

t (τ )

-0.887(0.365)

-2.397(0.033)

-4.598(0.000)

-6.074(0.000)

-3.309(0.013)**

half-life

0.9885

0.8327

0.5599

0.4240

0.6100

Fourier QKS

6.074***

3.125 (90% C.V.)

3.492 (95% C.V.)

Fourier QKS THAILAND

0.1

309.436

C S

I R

T P

4.529 (99% C.V.)

Note: ***, ** and * denote 1%, 5% and 10% significant level. Bootstrap procedures with 20000 replications are used to generate the p-value. The number in the brackets are the p-value.

A

U N

D E

M

T P

E C

C A

[32]

ACCEPTED MANUSCRIPT

Figure 1 The path of clean energy consumption with fitted intercepts with and without smooth breaks 2.4

2.8

2.2 2.4 2.0 1.8

2.0

1.6 1.6 1.4 1.2

1.2

1.0 0.8

0.6

0.4 65

70

75

80

85

90

95

00

05

10

15

65

70

75

80

85

90

95

1.2

1.8 0.8 1.6

1.2

0.0

1.0 -0.4 0.8

75

80

85

90

95

00

05

INDIA INDIA Fitted Intercepts without smooth breaks INDIA Fitted Intercepts with smooth breaks

10

15

15

65

70

75

80

85

90

95

00

05

10

15

INDONESIA INDONESIA Fitted Intercepts without smooth breaks INDONESIA Fitted Intercepts with smooth breaks

1.2

M

0.8

AN

-0.8

70

10

US

0.4

1.4

65

05

CR

2.0

0.6

00

CHINA CHINA Fitted Intercepts without smooth breaks CHINA Fitted Intercepts with smooth breaks

IP

BRAZIL BRAZIL Fitted Intercepts without smooth breaks BRAZIL Fitted Intercepts with smooth breaks

T

0.8

0.6 0.4 0.2

1.0 0.8

ED

0.6

0.0 -0.2 -0.4

-0.8 -1.0 65

PT

-0.6

70

75

80

85

90

0.4 0.2 0.0 -0.2 -0.4

95

00

05

10

15

65

70

75

80

85

90

95

00

05

10

15

10

15

PAKISTAN PAKISTAN Fitted Intercepts without smooth breaks PAKISTAN Fitted Intercepts with smooth breaks

CE

MALAYSIA MALAYSIA Fitted Intercepts without smooth breaks MALAYSIA Fitted Intercepts with smooth breaks

1.6

1.0 0.8

1.2

0.6

0.8

AC

0.4

0.4

0.2 0.0

0.0

-0.2

-0.4

-0.4

-0.8

-0.6

-1.2

-0.8 65

70

75

80

85

90

95

00

05

10

15

PHILIPPINES PHILIPPINES Fitted Intercepts without smooth breaks PHILIPPINES Fitted Intercepts with smooth breaks

65

70

75

80

85

90

95

00

05

THAILAND Thailand Fitted Intercepts without smooth breaks Thailand Fitted Intercepts with smooth breaks

Note: The blue line is the clean energy consumption, the black line is fitted intercepts value without smooth breaks, the red line is fitted intercepts with smooth break.

[33]

ACCEPTED MANUSCRIPT Highlights The unit root hypothesis of clean energy consumption is first visited.

2.

Quantile unit root test presents asymmetric performance.

3.

Smooth breaks are considered.

4.

The deterministic trend is closely fitted to clean energy consumption.

AC

CE

PT

ED

M

AN

US

CR

IP

T

1.

[34]