El Niño, La Niña, and the Forecastability of the Realized Variance of Heating Oil Price Movements

: We use the heterogenous autoregressive (HAR) model to compute out-of-sample forecasts of the monthly realized variance (RV) of movements of the spot and futures price of heating oil. We extend the HAR–RV model to include the role of El Niño and La Niña episodes, as captured by the Equatorial Southern Oscillation Index (EQSOI). Using data from June 1986 to April 2021, we show evidence for several model conﬁgurations that both El Niño and La Niña phases contain information useful for forecasting subsequent to the realized variance of price movements beyond the predictive value already captured by the HAR–RV model. The predictive value of La Niña phases, however, seems to be somewhat stronger than the predictive value of El Niño phases. Our results have important implications for investors, as well as from the perspective of sustainable decisions involving the environment.


Introduction
As pointed out by Trenberth et al. [1], the El Niño Southern Oscillation (ENSO) is an irregularly periodic variation in winds and sea surface temperatures over the tropical eastern Pacific Ocean, which, in turn, impacts the climate of much of the tropical and subtropical zones. In this regard, the warming phase of the sea temperature is known as El Niño, while the cooling phase is called La Niña, with the two phases lasting several months each and typically occurring every few years with varying intensities. Because El Niño phases tend to be warmer and drier overall, people will require less heating and more cooling, while the opposite is true for La Niña phases. As a result, El Niño and La Niña phases are likely to cause fluctuations in the price of heating oil, given its crucial role as a source of heat in the winter months. Not surprisingly, while analyzing the role of the ENSO on the prices of oil, gasoline, heating oil, diesel fuel, jet fuel and propane, [2] detect the strongest impact for heating oil.
Against the backdrop, we build on this line of research by verifying whether this evidence of the impact of ENSO on the first moment of the price of heating oil can be extended to out-of-sample forecasting of its second moment (that is, the variance of price movements), which is likely to be affected via variable demand for heating oil during El Niño and La Niña phases, associated with extreme warm and cold temperatures during a year. In this regard, we can further hypothesize that with the La Niña events associated with extreme cold weather, they are likely to have a relatively stronger influence on the volatility of heating oil than the corresponding El Niño phases. Moreover, the ENSO has been held (partially) responsible for some of the world's greatest humanitarian disasters, with [3,4] documenting the global El Niño droughts of 1876-1878, 1888-1891, and [1896][1897][1898][1899][1900][1901][1902], which contributed to the deaths of between 30 and 60 million people in India, China and Brazil over the 26 year period. Moreover, these extreme weather conditions can constrain the supply of rain-driven agricultural commodities and induce food price and generalized inflation [5][6][7][8], which may trigger social unrest as shown empirically by [9] and both theoretically and empirically in [10]. In this regard, [11][12][13][14] show that El Niño and La Niña shocks may even have played a role in a substantial number of civil conflicts. In this sense, the ENSO cycle can serve as a proxy for rare disaster risks [15], which in turn contribute to jump risks [16,17] that account for a large part of the volatility process of asset and commodity prices [18][19][20]. Naturally, one can expect an impact on the variance of price movements of heating oil via the volatility jumps channel as well, over and above the heating oil demand fluctuations route due to El Niño and La Niña events.
Yet, from a statistical perspective, it is well-established that in-sample predictability does not necessarily guarantee out-of-sample predictability gains emerging from a specific predictor, and it is in fact the latter approach that often provides a more robust test of the appropriateness of an econometric model [21]. Hence, rather than focusing on in-sample predictions, we use an out-of-sample forecasting experiment to test the predictive value of El Niño and La Niña for the variance of movements of the price of heating oil. At the same time, heating oil is second only to gasoline in terms of products derived from crude oil, and given that commodity market participants need real-time volatility forecasts for effective investment, hedging risk and arbitrage strategies, the economic importance of the question we ask in this paper cannot be overemphasized. This is more so given the close linkage between energy-based commodities and the financial markets in the wake of the process of financialization of the former, and with investors looking for alternative safe investment assets following the Global Financial Crisis (see, for example, [22,23] for detailed discussions in this regard).
Given this, we forecast the monthly realized variance (RV) of movements of the spot and future heating oil price using an extended version of the widely-studied heterogeneous autoregressive (HAR)-RV model [24]. We extend the model to incorporate the role of El Niño and La Niña phases that occurred during the period from June 1986 to April 2021. In this regard, it is important to note that measuring variance in terms of RV, defined as the sum of squared daily returns of prices observed over a given month (see [25]), provides an observable and unconditional metric of variance, which is otherwise a latent process. Accordingly, we differ from the existing literature (see for example, [26][27][28][29][30][31][32]) on modeling and forecasting heating oil volatility based on various types of univariate GARCH models, under which the conditional variance is a deterministic function of model parameters and historical data, and, hence, is not model-free as in the case of RV. Moreover, as discussed in [33,34], the benchmark HAR-RV model captures the generally observed long-memory and multi-scaling properties of volatility, despite having a simplistic structure. Notably, one important feature of the HAR-RV model is that it employs volatilities from different time resolutions to forecast the realized heating oil price volatility. Accordingly, the model captures the main idea motivating the heterogeneous market hypothesis (EMH, [35]) in a simple and straightforward way. The EMH indicates that different groups of market participants populate financial markets (that is, in our case, the market for heating oil), and that these groups differ in their sensitivity to information flows at various time horizons. For example, in the energy market, speculators and traders are concerned with shortterm investment horizons, whereas institutional investors are concerned with long-term investment horizons.
At this stage, to emphasize our contribution, it is important that we provide some further details on the literature cited above dealing with heating oil volatility forecasting, with these studies also simultaneously often dealing with crude oil, natural gas and gasoline. While in what follows we present a brief outline of the existing literature on forecasting heating oil price volatility, the reader is referred to [36] for a more detailed discussion of the literature on volatility forecasting of multiple energy commodities.
In this regard, Sadorsky cite26 showed that within the GARCH family models, the threshold GARCH (TGARCH) framework proved to be the most appropriate for the volatility of heating oil futures returns over the daily periods from 1988 to 2003. Nomikos and Pouliasis [27] compared the forecasting ability of a Markov-Switching (MS)-GARCH model to the basic GARCH on a dataset containing information about daily futures prices from 1991 to 2008 for heating oil, to depict the dominance of the former due to its ability to account for the special volatility characteristics. Arouri [28], while analyzing both daily spot and futures prices from 1986 to 2011, highlighted the role of a Fractionally Integrated GARCH (FIGARCH) model in its ability to better describe the time-varying volatility process, as it can capture long-memory (which decreases after incorporating structural breaks), a phenomenon we also account for via the HAR-RV model. Given this, the FIGARCH model (suitably modified to account for structural breaks) is shown to perform better than any other related model in forecasting heating oil volatility. Along similar lines, ref. [29] showed that the Autoregressive Fractionally Integrated Moving Average (ARFIMA)-FIGARCH model dominated various other GARCH-class models in forecasting volatility for heating oil derived from daily futures data covering 1995-2012. Different to the existing studies, ref. [31] pointed out that the unbiased extreme value volatility estimator, namely the Add RS, when incorporated into an ARFIMA model of heating oil futures volatility, outperforms the alternative models from the GARCH family over the daily period from 1996 to 2013. In a similar vein, in terms of methodological innovation, ref. [32] compared support vector regression (SVR) with the GARCH models, and showed the superiority of the former approach in forecasting the volatility of heating oil futures based on daily data from 2015 to 2019. However, this study also indicated that, within the GARCH models, asymmetric versions tend to perform better-an observation drawn earlier by [26]. Interestingly, while all the above studies relied on univariate models, ref. [30] stressed the need to look at multivariate (the Full BEKK [37], the Diagonal BEKK and Scalar BEKK) models while forecasting weekly spot prices of heating oil, as these frameworks capture the time-varying correlation between the returns processes of the energy-based commodities.
Given the brief description of the academic literature, it appears that our current paper is the first paper that evaluates the out-of-sample forecasting power of the El Niño and La Niña phases for the realized variance of heating oil price movements, which are likely to be affected via demand and volatility jumps channels, based on an extended HAR-RV model. Our focus is not necessarily to contribute to the technicality of volatility forecasting models, but rather to consider the usage of innovative predictors, namely the El Niño and La Niña phases.
Section 2 provides the dataset. Section 3 describes the HAR-RV model. Section 4 presents the results from our forecasting experiment. Section 5 concludes the paper.

Data
We use daily spot and futures prices of heating oil (U.S. dollars and cents per gallon) as represented by the New York Harbor No. 2 Heating Oil Spot Price FOB and the continuous NYMEX New York HARBOR ULSD futures contract, respectively, with the price data extracted from DataStream. After having computed the continuously compounded (i.e., log) daily returns of the spot and futures prices, we compute the monthly realized variance by taking the sum of daily squared returns over a month, where the overall sample period runs from June 1986 to April 2021. Figure 1 shows the daily price data that we used to construct the monthly realized variances. Table 1 depicts the key summary statistics of daily returns. Both returns series are slightly negatively skewed, are leptokurtic and they are slightly positively (negatively) skewed in the case of spot (futures) returns. The two top panels of Figure 2 show the monthly realized variances of the spot and futures returns of the heating oil price. While the overall pattern of RV is similar across the spot and futures data, there are also noticeable differences visible in the figure, for example, in the early 1990s, in 2000 and in 2020 at the end of the sample period.   In order to measure the ENSO cycle, we rely on the Southern Oscillation Index (SOI) as compiled and published by the Bureau of Meteorology, Government of Australia (http: //www.bom.gov.au/climate/current/soihtm1.shtml, accessed on 21 May 2021). The SOI summarizes the development and intensity of El Niño or La Niña events in the Pacific Ocean. From a methodological point of view, the SOI is computed by taking the pressure differences between Tahiti and Darwin. Such pressure differences are important because low atmospheric pressure tends to occur over warm water, while high pressure tends to occur over cold water, in part because of deep convection over the warm water. El Niño phases are, by definition, periods of sustained warming of the central and eastern tropical Pacific Ocean. La Niña phases are periods of sustained cooling of the central and eastern tropical Pacific Ocean. As such, the result is a decrease and an increase in the strength of the Pacific trade winds, respectively. Sustained negative (positive) realizations of the SOI below (above) −7 (+7) usually indicate El Niño (La Niña) phases.
It should be mentioned, however, that a drawback of the SOI data is because both Darwin and Tahiti are located well south of the Equator, implying that the surface air pressure at both locations is less directly linked to ENSO. In order to address this drawback, a new index called the Equatorial Southern Oscillation Index (EQSOI) was constructed. For further details, see https://www.climate.gov/news-features/blogs/enso/why-arethere-so-many-enso-indexes-instead-just-one (accessed on 21 May 2021). The data for this index are collected at two new regions, both centered on the Equator, where the western one is located over Indonesia while the eastern one is located over the equatorial Pacific, close to the South American coast. We collect data on the EQSOI from the Climate Prediction Center of the National Oceanic and Atmospheric Administration (US Department of Commerce). See https://www.cpc.ncep.noaa.gov/data/indices/ (accessed on 21 May 2021).
Like [38], we first compute the SOI and EQSOI "anomalies", defined as the deviation of the SOI and EQSOI in any given month from its historical average, normalized (divided) by its historical standard deviation. The two middle panels of Figure 2 show the EQSOI data and the EQSOI anomaly. Equipped with the anomalies, we identify sustained negative SOI and EQSOI anomaly values below −1 (above +1) as indicators of El Niño (La Niña) phases. In order to capture the two phases of the ENSO cycle, we create a dummy variable which assumes the value of one whenever the SOI and EQSOI anomalies are negative (positive), and zero otherwise. Finally, we multiply the dummy variable with the SOI and EQSOI anomalies to obtain a metric for the intensity of the El Niño (La Niña) phases. In our empirical analysis, we use the El Niño and La Niña phases, depicted in the two lower panels of Figure 2, associated with EQSOI as a measure of the ENSO, given its stronger reliability. We use the El Niño and La Niña phases based on the SOI data for a robustness check.

The HAR-RV Model
The HAR-RV model proposed by [24] has become the workhorse model in empirical finance for the analysis of the dynamics of the realized variance (RV) of asset-price returns. While the HAR-RV model is typically applied to studying high-frequency data, we use a variant of the HAR-RV model to study the monthly realized variance (RV) of spot and futures movements of the price of heating oil. Our benchmark HAR-RV model is given by: where h = the forecast horizon, RV t+h = the average (spot or futures) realized variance observed during the forecast horizon, RV q,t = the average quarterly realized variance from month t − 3 to t − 1, and RV y,t = the average yearly realized variance from month t − 12 to t − 1). We estimate the following coefficients β j , j = 0, 1, 2, 3, and t = a disturbance term. The predictive value of the benchmark HAR-RV model is compared to that of each of the following extended HAR-RV models: In our empirical analysis, we focus on the out-of-sample predictive performance of our various HAR-RV models. We compute out-of-sample forecasts by estimating the models on rolling-estimation windows of different lengths. Equipped with the out-ofsample forecasts, we then evaluate the forecasts by means of the Clark and West [39] test. The null hypothesis of this test is that the two models that a researcher compares have an equal out-of-sample mean-squared prediction error (MSPE). The Clark-West test is easy to implement because it requires regressing the quantityf t+h = ( A hat over a variable marks a forecast of RV and the subindices A and B denote the two models being compared (B denotes the larger model). We reject the null hypothesis if the t-statistic of the constant in this regression model is significantly positive (one-sided test; we use Newey-West robust standard errors to compute the t-statistic). In addition, we use the root mean squared forecasting error (RMSFE) and the mean absolute forecasting error (MAFE) statistics to compare the forecasting performance of the various HAR-RV models. We carry out our forecasting experiments using the R language for statistical computing [40].
As an extension and a robustness check (and given the results of the causality-inquantiles test reported in Table A1 at the end of the paper), we consider a quantileregression version of the HAR-RV model (see [41,42], among others). The quantileregression variant of the HAR-RV model is tailored to study whether the predictive value of the El Niño and La Niña phases differs across the quantiles of the conditional distribution of the realized variance. The quantile-regression HAQR-RV model is estimated by minimizingb where ρ α = t+h (α − 1( t+h < 0) denotes the standard check function, α denotes a quantile, 1 denotes the indicator function, and b α denotes a vector of quantile-specific parameters to be estimated (a hat denotes an estimated parameter). In Equation (5), we simplify the notation by letting the matrix X t denote the predictors of the HAQR-RV model. In order to evaluate forecasts for the various quantiles of the conditional distribution of RV, we use an out-of-sample relative loss criterion [43,44]. The out-of-sample relative-loss criterion, R o , is defined as: whereR V t denotes the forecast derived from the benchmark model given in Equation (1), τ is the first period for which a forecast is available. It should be noted that the benchmark forecast is also indexed by the quantile parameter, α, showing that the benchmark model is also estimated as an HARQ-RV model. The parameter, p, governs the functional form used to compute the out-of-sample relative loss criterion. Specifically, in case we set p = 1, a forecaster evaluates forecast errors under a lin-lin loss function (L1; this case corresponds to the check function of the quantile regression framework), while in the case of p = 2, a forecaster uses a quad-quad loss function (L2) to evaluate forecast errors [45].
It should be noted that the quantile parameter, α ∈ (0, 1), governs the asymmetry of the loss function. The loss function is symmetric for α = 0.5 and is symmetric and quadratic if, in addition, p = 2. Setting α = 0.5 and p = 1, in turn, implies that the the loss from a forecast error increases in the absolute magnitude of the forecast error. If α > 0.5 (α < 0.5), then the loss from underpredicting (overpredicting) RV exceeds the loss from an overprediction (underprediction) of the same size, that is, for α = 0.5 the loss function is asymmetric.
For a given quantile parameter and a given functional form of the loss function, the HARQ-RV model, extended to include the El Niño and/or La Niña phases, outperforms the benchmark HARQ-RV model when the out-of-sample relative loss criterion is positive, R o > 0. When the out-of-sample relative loss criterion assumes a negative value, R o < 0, the benchmark model performs better than the extended model. Table 2 summarizes our main forecasting results. We report, for the EQSOI data, results for the CW test (p-values) for different forecast horizons (h = 1, 3,6,12,24,36,48) and various lengths of the rolling estimation window (48, 60, 72, 96 and 120 months). In Panels A to C, we document the results for spot RV, while Panels D to F contain the results for futures RV. In the case of spot RV (Panel A), we observe that El Niño phases add to the predictive performance of the HAR-RV model when we use an intermediate window length and forecast horizons of up to two years. The test results become largely insignificant for the short and the long window length and for the longer forecast horizons. When we consider La Niña phases in Panel B, in contrast, most test results are significant. Only for a combination of a long estimation window and a short to intermediate forecast horizon are the test results not statistically significant. A similar result emerges when we combine El Niño and La Niña phases in a unified forecasting model (Panel C). It is interesting to examine whether the test results are robust to a change in the definition of realized variance. To this end, we report in Table 3 the test results for the square root of realized variance, √ RV, which researchers in empirical finance also often consider as a metric of realized "volatility". The results that we document in Panel A witness that the contribution of El Niño phases to the predictive performance of the HAR-RV model, in the case of spot √ RV, is concentrated at the shorter forecast horizons. The overall picture that emerges, however, is that La Niña phases (Panel B) have a more systematic statistically significant effect on the forecasting performance of the HAR-RV model than on El Niño phases. Only the test results for the long rolling estimation window are not significant. We conclude that La Niña phases are a major driver of the results, especially as far as the longer forecast horizons are concerned, for the combined model in Panel C. For the futures-based √ RV (Panels D to F), the results are not too different. We observe that the contribution of El Niño phases to the forecasting performance of the HAR-RV model, in terms of statistical significance, is concentrated at the two short rolling estimation windows and at the intermediate forecast horizons. The overall picture that emerges from a comparison with the results for La Niña phases in Panel D, however, is that, with the exception of the long rolling estimation window, the test results for La Niña phases are statistically significant for a broader array of combinations of window length and forecast horizons. As a result, it is not surprising that the test results in Panel F, for a model that features both El Niño and La Niña phases as predictors, are largely significant.   Our empirical findings show that both El Niño and La Niña have predictive value in an out-of-sample forecasting experiment for the realized variance (and volatility) of heating oil price movements, as measured by means of spot and futures price data. On balance, however, the predictive value of La Niña appears to be somewhat more robust in our out-of-sample forecasting experiment across the various model configurations that we have studied than the predictive value of El Niño. The results that we report in Table 5 illustrate this finding. Specifically, we report in Table 5 (for the EQSOI data in Panel A and the SOI data in Panel B) the average RMSFE statistics and the average ratios of the MAFE statistics as computed, for every model, across the five different lengths of the rolling-estimation window and the six different forecast horizons. The ratios are defined by dividing the RMSFE (MAFE) obtained for the baseline HAR-RV model by the corresponding RMSFE (MAFE) obtained for the HAR-RV model extended to include the El Niño and/or La Niña phases. Hence, a ratio larger than one indicates that the extended model produces, on average, more accurate forecasts than the baseline HAR-RV model. The punchline of the results summarized in Table 5 is that extending the the HAR-RV model to include El Niño phases tends to improve forecast accuracy, but extending the HAR-RV model to incorporate La Niña phases leads to a larger improvement of forecast accuracy. In fact, including both El Niño and La Niña phases in the forecasting model yields RMSFE and MAFE ratios that are close to the ratios reported for the forecasting model that features the La Niña phases as the only additional predictor. We next turn to a quantiles-based analysis. As for in-sample predictability, we find, based on the multivariate version of the the nonparametric k-th order causality-in-quantiles test of [46], which allows us to test for causality in (spot and futures) returns and squared returns (i.e., variance) simultaneously over their respective conditional distributions (thus defining the states of the heating oil market) strong evidence of in-sample predictability from both El Niño and La Niña phases, when controlling for La Niña or El Niño in the model, respectively-see the results reported in Table A1 at the end of the paper. The results further show that the evidence of in-sample predictability is not concentrated at one or a few quantiles, justifying our decision to use the HARQ-RV model to compute quantilespecific forecasts of RV, and to use the out-of-sample relative-loss criterion to evaluate the forecasts relative to the forecasts we obtain from the benchmark HARQ-RV model. In order to save journal space, we compare the models that feature either El Niño or La Niña phases as an additional predictor. Figures 3 and 4 summarize the results for the spot RV under the lin-lin and quad-quad loss function. For the lin-lin loss function, the colored areas depicted in the figure show that, according to the out-of-sample relative-loss criterion, using El Niño phases contributes to the overall forecasting performance of the HARQ-RV model mainly at the short to intermediate forecast horizons, while the areas where the benchmark HARQ-RV model performs better than the extended model become larger as the forecast horizon increases. The results for the La Niña phases, in contrast, show that the contribution of this predictor to the forecasting performance of the model increases in terms of the out-of-sample relative-loss criterion when we study the intermediate to long forecast horizons. For the quad-quad loss function, in turn, the HARQ-RV model extended to include La Niña phases as an additional predictor outperforms the benchmark HARQ-RV model in terms of the out-of-sample relative-loss criterion for most combinations of the quantile parameter and the length of the rolling estimation window for all forecast horizons being studied. A key message conveyed by the results is that the overall contribution of La Niña phases to the forecasting performance of the model is stronger than that of El Niño phases. In this regard, it should also be noted that, as indicated by the legends of the figures, the maximum gains a forecaster can realize in terms of the out-of-sample relative-loss criterion are often larger for La Niña phases than for El Niño phases.

HAR-RV vs. HAR-RV-EL-NINO / L2 / h=48
Rolling-window length In this regard, we note that the higher values of the statistic for the k-th order nonparametric causality-in-quantiles test, as reported in Table A1, confirm the finding of a stronger impact of La Niña phases relative to El Niño ones, especially at higher conditional quantiles, which corresponds to higher returns and volatility. In fact, based on the average derivative (AD) using the conditional pivotal quantile, based on the approximation or the coupling approach of [47], both El Niño and La Niña phases are found to have a positive relationship with the entire conditional distributions of returns and the volatility of heating oil price. These results are available from the authors upon request. Figures A1 and A2 at the end of the paper (Appendix A) summarize the results for the futures RV. Again, we observe that the range of combinations of the quantile parameter and the length of the rolling estimation window that result in positive realizations of the out-of-sample relative-loss criterion tends to be wider for La Niña phases than for El Niño phases. Hence, to summarize, the message to take home from the results is that the out-of-sample relative loss criterion takes on positive values for many combinations of the quantile parameter, window length, and the the forecast horizon being studied, especially as far as the La Niña phases are concerned.
Finally, we present in Table 6 the results for upside ("good") and downside ("bad") RV, where the former is computed from positive daily returns, while the latter is computed from negative daily returns. This is an important issue, as [48] stresses that financial market participants care not only about the nature of volatility, but also about its level, with traders typically differentiating between good and bad volatilities, which, in turn, provides us the motivation to look at such disaggregation of RV. We present in Table 6 the results for the futures RV, but results for spot RV are qualitatively similar (and available from the authors upon request). The results show that El Niño phases add to the forecasting performance of the HAR-RV model estimated for upside and downside RV only for the two shortest rolling estimation windows. La Niña phases, in contrast, also carry information for several forecast horizons, which is useful for forecasting upside and downside RV when we study the longer rolling-estimation windows, in the case of downside RV mainly for the longer forecast horizons. As a consequence, it is not surprising that using both El Niño and La Niña phases as predictors of upside/downside RV results in significant test results for several combinations of the rolling estimation window and the forecast horizon.

Conclusions
Our empirical findings show that both El Niño and La Niña phases contain predictive value for the realized variance (and volatility) of heating oil price movements. Across the various model configurations that we have studied, however, the predictive value of La Niña phases appears to be stronger than the predictive value of El Niño phases in terms of both statistical significance of the test results and in terms of forecast accuracy, which is not surprising given that the latter is associated with the cooling of temperatures and, hence, greater demand for heating oil and a stronger impact on volatility via trading in the market.
Our results are useful for investors who need to obtain reliable information on the response of the future path of heating oil price fluctuations to changes in climate patterns for portfolio management and hedging decisions. Our results also have implications from the context of "sustainability"; in this case related to the environment. Note that La Niña phases lead to a greater demand for heating oil, and hence drive the volatility of heating oil price. Now, the rise in heating oil demand is expected to cause an increase in the emissions of carbon dioxide, which is known to result in global warming [49][50][51][52]. The implication of this is that the warming of the ocean's surface is faster than the subsurface, which is likely to induce occurrences of the El Niño phenomena [53]. In this regard, the National Oceanic and Atmospheric Management Climate Prediction Center states that global warming makes the extreme weather more severe, including an increase in the frequency of the super El Niño phenomenon. Therefore, there is a close feedback relationship between the El Niño phenomenon and the movements in prices in the heating oil market, resulting from its high demand due to La Niña phases. Naturally, policy authorities need to ensure stability in the price of the heating oil market, which would only be possible if the demand is kept from fluctuating wildly in the wake of La Niña events, and this would entail investing more in greener, more sustainable and renewable energy sources such as solar, wind, hydro, ocean, geothermal, biomass and hydrogen.
In future research, it will be interesting to build on our empirical research and examine how the returns and the realized variance of the prices of non-energy (for example, agricultural) commodities respond to El Niño and La Niña phases, and whether these phases help to forecast price movements. Such a forecasting exercise for the returns and the variance should be particularly useful for policymakers in developing countries, and in countries whose exports consist, to a large extent, of agricultural commodities.  Acknowledgments: We thank two anonymous reviewers for helpful comments. The usual disclaimer applies.

Conflicts of Interest:
The authors declare no conflict of interest.
Code Availability: Code is available from the authors upon request.

HAR-RV vs. HAR-RV-LA-NINA / L2 / h=48
Rolling-window length α Figure A2. Results for the HARQ-RV Model (Futures; Quad-Quad Loss Function). Note: This figure reports the out-ofsample relative loss criterion for a quad-quad loss function (L2) and an extended HARQ-RV models that features either El Niño or La Niña phases as an additional predictor. The baseline model is the benchmark HARQ-RV model. White spaces indicate areas where the out-of-sample relative loss criterion takes on negative values (the baseline model performs better than the extended model). The rolling-estimation window is varied from 48, 49, 50, ... to 120 observations. The results are based on the SOI data. The parameter h denotes the forecast horizon (in months).