Most existing literature focus on investigating long-term memory in stock markets returns, fixed income markets returns, or commodity markets returns. Mixed results were found on the presence of apparent temporal dependencies in financial market volatility. Several studies conclude that capital markets are characterized by long memory processes  -  . However, many studies do not find any significant and robust evidence of positive long-term persistence in the financial markets  -  . Simultaneously, other stream of literature finds temporary or little evidence of long-term memory in different stock markets  -  .
Moreover, studies that investigate the temporal dependencies in financial market volatility, employ various methodologies such as, classical rescaled-range (R/S) analysis   , modified rescaled-range (R/S) analysis  , the spectral regression method  , and different GARCH specifications. However, a systematic review of the literature indicates that FIGARCH models outperform many of the other conditional heteroscedastic models in predicting and modeling different classes of assets, such as stock returns   , exchange rate returns    , and futures returns   in different market settings. Again, as per the best of authors’ knowledge, not sufficient applications of volatility models are available in the existing literature.
We believe that the existing literature on investigating market efficiency and/or inspecting the temporal dependencies lead to very different or even contradicting conclusions, and hence still need further examination. And, surprisingly, the Energy Exchange-Traded Funds (ETFs) do not attract the attention of researchers despite its diverse nature and potential for future investors.
This study aims at investigating the market efficiency of ETFs. We estimate FIGARCH model proposed by Baillie et al.  to examine the long memory characteristics in the ETFs volatility of both renewable and unrenewable energy ETFs using daily returns calculated by Thomson Reuters Eikon. We use nine major indexes of renewable and unrenewable energy ETFs, namely, the Fidelity MSCI Energy Index ETF (FENY), First Trust Energy AlphaDEX Fund (FXN), iShares US Oil & Gas Exploration & Production ETF(IEO), iShares Global Energy ETF (IXC), iShares US Energy ETF (IYE), SPDR S&P North American Natural Resources ETF (NANR) and VanEck Vectors Oil Services ETF (OIH), Vanguard Energy Index Fund ETF (VDE), Energy Select Sector SPDR Fund (XLE), to inspect the temporal dependencies in depth.
We adopt FIGARCH model to investigate the long-range dependence in ETFs volatility of both renewable and unrenewable energy ETFs. We find evidence of long-term memory in the volatility of all the ETFs. This implies that all the ETFs under investigation are weak-form inefficient funds. The results also indicate that the volatility has a predictable structure in all the ETFs of both renewable and unrenewable energy ETFs, indicating the potential of diversification for the international investors. We argue that a better understanding of the long-range dependence in ETFs volatility (long memory processes), within the EFTs market, is inevitable for international investors, multinational corporations, and portfolio managers, to achieve superior diversification and manage their financial risk exposure.
The rest of the paper is organized as follows: next section describes the FIGARCH model used to study long-range dependence in EFTs returns. Section 3 provides a detailed outlook of the data. In Section 4, we present and discuss the results of the empirical analysis and finally Section 5 concludes the paper.
2. Model Specifications
The Autogressive Conditional Heteroscedasticity (ARCH) process proposed by Engle  and generalized ARCH (GARCH) by Bollerslev  are well known for volatility modeling and forecasting of stock returns. More precisely, to model the most prominent features of the time series data (also called stylized facts) such as volatility clustering, excess kurtosis, and fat-tailedness. However, to explain how persistent volatility is, GARCH process can easily be extended to identify the long memory process, a common observation in actual data, through a fractionally integrated procedure proposed by Baillie et al.  . Specifically, the Fractionally Integrated Generalized Autoregressive Conditional Heteroscedastic (FIGARCH) process.
We start our empirical specification with the GARCH (p, q) process introduced by Bollerslev  , we can write the conditional variance as:
The primary constraint of this model is that all the expounding variables must be positive i.e., this is known as the non-negativity restriction. Further, for stationarity we require that is less than unity. However, if this restriction violates, i.e., we conclude that the shocks are persistent. Hence, to account for the persistency of shocks an IGARCH (1, 1) model proposed by Engle and Bollerslev  can be written as
The IGARCH model implies infinite persistence of the conditional variance to a shock in squared returns. The IGARCH process can also be illustrated as an ARMA (m, p) process
The fractionally integrated GARCH or FIGARCH class of models is obtained by replacing the first difference operator (1 − L) in above model with the fractional differencing operator (1 − L)d. where d is a fraction 0 < d < 1. Thus, the FIGARCH class of models can be obtained by considering following equation.
The FIGARCH process identifies potential presence of long memory or the subsistence of dependencies in financial time series mainly due to the hyperbolically decaying autocorrelation function, or in other words, long memory process can be illustrated through a fractionally integrated procedure, Meanings, the level of integration is fewer than one however superior than zero, implying that the impact of a shock continue over an extensive period of time. The main advantage of FIGARCH process is that it allows for long memory in the conditional variance which is characterized by the fractional integration parameter d and the short-term dynamics can be modeled through the traditional GARCH parameters. Following Baillie et al.  we adopt the Quasi maximum likelihood estimation (QMLE) technique.
3. Data and Descriptive Statistics
The data comprise daily returns of both renewable and unrenewable energy ETFs. We use nine major indexes of renewable and unrenewable energy ETFs (namely, the Fidelity MSCI Energy Index ETF (FENY), First Trust Energy AlphaDEX Fund (FXN), iShares US Oil & Gas Exploration & Production ETF(IEO), iShares Global Energy ETF (IXC), iShares US Energy ETF (IYE), SPDR S&P North American Natural Resources ETF (NANR), VanEck Vectors Oil Services ETF (OIH), Vanguard Energy Index Fund ETF (VDE), Energy Select Sector SPDR Fund (XLE), to inspect the temporal dependencies in depth.
The beginning of the sample period is dictated by the availability of data for each index investigated. The end of the period is May 2017 for all indices. Figure 1 shows the development of different indexes of both renewable and unrenewable energy ETFs. All the data are retrieved from Thomson Reuters Eikon and daily returns are constructed as the first difference of logarithmic prices multiplied by 100.
Before formal investigation of long memory in energy ETFs, we inspect the time-series properties of our data set using primary techniques, for instance, Stationarity in the time series is checked by applying the Augmented Dickey Fuller (ADF) test. To check the null hypothesis of normal distribution we calculate Jarque-Bera test statistic. Finally, to investigate the null that autocorrelation coefficients up to 20 lags are zero, we compute Ljung and Box  test statistic, together with the ARCH LM-statistic (five lags) on each returns series. The results shown in Table 1, in general, support the findings of prior studies and allow us to reject the null hypothesis that returns have unit root in favor of alternate hypothesis of stationarity, P-Values are reported in the table (even at 1% MacKinnon critical value). The Jarque-Bera normality and Engle’s Lagrange Multiplier ARCH tests both revels that the return data of all the energy ETFs exhibit non-normality and ARCH effects, P-Values are reported in the table. These primary findings grant confirmation against the market efficiency hypothesis and allow us to use GARCH specification through LB statistics and ARCH LM-statistic.
4. Empirical Results
Using the Energy Exchange-Traded Funds, we re-examine the subject of whether or not actual returns reveal temporal dependence. We utilize both renewable and unrenewable energy ETFs. Our empirical investigation is based on the GARCH family models. First, to model the volatility dynamics of ETFs, we utilize the traditional GARCH (1, 1) model. Second, to check the volatility persistence, we adopt IGARCH (1, 1) technique. Finally, to investigate the long run dependence in ETFs, we employ FIGARCH (1, 1) framework. The results of estimated GARCH, IGARCH and FIGARCH models are reported in Table 2.
Panel A of Table 2 presents GARCH (1, 1) estimations of all the ETFs under investigation. The results show that both the ARCH and GARCH parameters
Figure 1. Development of different energy ETFs.
Table 1. Descriptive statistics of Energy Exchange-Traded Funds for the daily return indices.
Table 2. Long memory in ETFs market volatility, estimated from a univariate GARCH (1, 1), IGARCH (1, 1) and FIGARCH (1, 1) model of daily return indices.
(α and β) are statistically significant for all the ETFs under investigation, which confirm the existence of the time-varying conditional variance. It is also evident from Panel A of Table 2 that the parameters of the conditional variance equations are all positive and meet the positivity constraint for the GARCH (1, 1) specification. However, the sum of α and β parameters is very close to the unity, indicating the persistence of the volatility in all the indices. One shortcoming of the traditional GARCH model is its failure to capture long-range dependence or to account for persistence of volatility in the data. Hence, we utilize IGARCH process proposed by Engle and Bollerslev  . The IGARCH model implies infinite persistence of the conditional variance to a shock in squared returns. Panel B of Table 2 presents IGARCH (1, 1) estimations of all the ETFs under investigation. The results are very similar to the standard GARCH (1, 1) estimations presented in Panel A, confirming the temporal dependencies in all the ETFs.
Finally, FIGARCH (1, 1) model is employed in order to investigate the existence of possible temporal dependencies in the volatility of all the ETFs under investigation. The FIGARCH process identifies potential presence of long memory or the subsistence of dependencies in financial time series mainly due to the hyperbolically decaying autocorrelation function. Results from this model are shown in Panel C of Table 2. As per our results, the fractional differencing parameter, d, is found to be significantly different from zero and is within the theoretical value (i.e., 0 < d < 1). This indicates that the volatility of all the ETFs under investigation clearly exhibits a long memory process. It is our connotation that our findings show the importance of modelling long memory in volatility and suggests that future volatility depends on its past realizations and, as a result, is predictable. Our findings also support the findings of prior studies on both stock and commodity markets.
To conclude, we report the sample skewness and kurtosis for the standardized residuals, (denoted by b3 and b4 in Table 2), also Ljung-Box portmanteau tests for up to 20th-order serial correlation in the standardized and the squared standardized residuals (denoted by Q20 and Q220 in Table 2) as diagnostic tests for all three models. While comparing different GARCH family models based on diagnostic tests, we found FIGARCH model performs better than the other two models, which is again consistence with the most recent prior research on the topic.
In summary, we provide evidence of long memory in the volatility of all the ETFs, which suggest that, all the ETFs under examination are week form inefficient. This is an evidence of violation of efficient market hypothesis, which can lead to the arbitrage opportunities for international investors who are interested to invest in one of the fastest growing market―Energy Exchange-Traded Funds. Furthermore, our results confirm that the volatility has a predictable structure in all Energy Exchange-Traded Funds, indicating the need of regulatory and economic reforms within the Energy Exchange-Traded Funds system.
In this study, we attempt to re-examine the market efficiency of Energy Exchange-Traded Funds. Specifically, we look at the long memory in ETFs market volatility. To inspect the temporal dependencies in depth we utilize nine major ETFs, representing both renewable and unrenewable EFTs. To study the long memory we estimate FIGARCH model proposed by Baillie et al.  using daily returns calculated by Thomson Reuters Eikon. We find evidence of long memory in the volatility in all the ETFs. This implies that all the ETFs under investigation are week form inefficient. Our results show that the volatility in different ETFs has a predictable structure. Our results indicate the need of regulatory and economic reforms within the Energy Exchange-Traded Funds system. As per our empirical investigation FIGARCH model performs better than the tradition GARCH models.
 Cajueiro, D.O. and Tabak, B.M. (2004) Evidence of Long Range Dependence in Asian Equity Markets: The Role of Liquidity and Market. Physica A: Statistical Mechanics and its Applications, 342, 656-664.
 Panas, E. (2001) Estimating Fractal Dimension Using Stable Distributions and Exploring Long Memory through ARFIMA Models in Athens Stock Exchange. Applied Financial Economics, 11, 395-402.
 Chow, K.V., Pan, M.-S. and Sakano, R. (1996) On the Long-Term or Short-Term Dependence in Stock Prices: Evidence from International Stock Markets. Review of Quantitative Finance and Accounting, 6, 181-194.
 Cajueiro, D.O. and Tabak, B.M. (2005) Possible Causes of Long-Range Dependence in the Brazilian Stock Market. Physica A: Statistical Mechanics and Its Applications, 345, 635-645.
 Geweke, J. and Porter-Hudak, S. (1983) The Estimation and Application of Long Memory Time Series Models. Journal of Time Series Analysis, 4, 221-238.
 Beine, M., Laurent, S. and Lecourt, C. (2002) Accounting for Conditional Leptokurtosis and Closing Days Effects in FIGARCH Models of Daily Exchange Rates. Applied Financial Economics, 12, 589-600.
 Baillie, R.T., Bollerslev, T. and Mikkelsen, H.O. (1996) Fractionally Integrated Generalized Autoregressive Conditional Heteroskedasticity. Journal of Econometrics, 74, 3-30.
 Antonakakis, N. and Darby, J. (2013) Forecasting Volatility in Developing Countries’ Nominal Exchange Returns. Applied Financial Economics, 23, 1675-1691.
 Baillie, R.T., Han, Y.W., Myers, R.J. and Song, J. (2007) Long Memory Models for Daily and High Frequency Commodity Futures Returns. The Journal of Futures Markets, 27, 643-668.
 Jin, H.J. and Frechette, D. (2004) Fractional Integration in Agricultural Futures Price Volatilities. American Journal of Agricultural Economics, 86, 432-443.