Explosive Episodes and Time-Varying Volatility: A New MARMA–GARCH Model Applied to Cryptocurrencies

Financial assets often exhibit explosive price surges followed by abrupt collapses, alongside persistent volatility clustering. Motivated by these features, we introduce a mixed causal–noncausal invertible–noninvertible autoregressive moving average generalized autoregressive conditional heteroskedasticity (MARMA–GARCH) model. Unlike standard ARMA processes, our model admits roots inside the unit disk, capturing bubble-like episodes and speculative feedback, while the GARCH component explains time-varying volatility. We propose two estimation approaches: (i) Whittle-based frequency-domain methods, which are asymptotically equivalent to Gaussian likelihood under stationarity and finite variance, and (ii) time-domain maximum likelihood, which proves to be more robust to heavy tails and skewness—common in financial returns. To identify causal vs. noncausal structures, we develop a higher-order diagnostics procedure using spectral densities and residual-based tests. Simulation results reveal that overlooking noncausality biases GARCH parameters, downplaying short-run volatility reactions to news (α) while overstating volatility persistence (β). Our empirical application to Bitcoin and Ethereum enhances these insights: we find significant noncausal dynamics in the mean, paired with pronounced GARCH effects in the variance. Imposing a purely causal ARMA specification leads to systematically misspecified volatility estimates, potentially underestimating market risks. Our results emphasize the importance of relaxing the usual causality and invertibility assumption for assets prone to extreme price movements, ultimately improving risk metrics and expanding our understanding of financial market dynamics.