1 Introduction

The COVID-19, which began in Wuhan, China, in December 2019, has rapidly evolved from a provincial health scare to a global meltdown, throwing the world to fear. According to the European Commission for Disease Control (ECDC 2020), this pandemic has now spread to more than 200 countries around the globe causing a severe global economic recession that began affecting the world on 20 February 2020. It is significantly worsened with the exceptional increase in the number of cases and deaths caused by this pandemic, as by the end of July 2020 more than 17 million cases of COVID-19, as well as 600.000 deaths, have been detected worldwide. Opinions about the COVID-19 pandemic are currently divided. Some research organizations believe that it could trigger another global financial crisis, while others believe that the impact of the pandemic, if not quickly stopped, could be worse than the outbreak of the SARS in 2003 in China, the 2008 global financial crisis, and the Second World War combined.

The COVID-19 pandemic generated a human and health crisis. The required measures to contain the virus have resulted in an economic downturn. The International Monetary Fund’s rapport indicated that, on 14 April 2020, all the G7 countries have already reached a “deep recession”. The IMF has announced that global growth, which is estimated at −3% in 2020, is “much worse” than it was during the 2009 Great Recession. According to Salisu et al. (2020) and Ghorbel and Jeribi (2021a, 2021b), the exclusion of the oil-stock nexus study from other variables during the pandemic is intentional. As the COVID-19 pandemic rages, oil prices have fallen by 30 percent, the biggest decline since the Gulf War of 1991 (Schneider & Domonoske, 2020). Also, the collapse of oil prices caused by the 2020 Russia–Saudi Arabia oil price war opened a new battlefront for pandemic-battered economies. In the absence of a major Middle East conflict, Jefferson (2020) foresees that crude oil prices may be expected to recover by the third quarter of 2020.

The negative economic consequences of this pandemic are significant for the G7 countries. There is a disruption in global supply chains. There is also a decline in demand for imported goods and services and a decrease in international tourism and business travel. We note that the global economy has become much more integrated if we compare this situation with the SARS outbreak in 2003. China, as the root of this pandemic, now plays a far more pronounced role in the industry, commerce, tourism, and FDI. China is now a manufacturing and a top power trade giant. The economic consequences of China’s negative shock would have a considerable impact on the global economy. Investors tend to look for new investment options that can deliver diversification and/or hedge benefits following periods of financial volatility experienced during the last decade. Financial markets have become more volatile and unstable due to the great uncertainty of the COVID-19 pandemic and its related economic losses (Zhang et al., 2020). At the beginning of the pandemic outbreak, Rizwan et al. (2020) and Gorbel and Jeribi (2021b) revealed a significant rise in systemic risk among the G7 countries. However, compared to the pre-COVID-19 period, China is showing some recovery.

For decades and through various crises, the yellow metal has historically been viewed as a safe haven asset in periods of financial uncertainty (Bouri et al., 2020; Ghorbel & Jeribi, 2021a, 2021b; Jareno et al., 2020; Selmi et al., 2018; Vardar et al., 2018). Gold values have risen significantly since the global economic and financial crises, while other commodities have sustained declines (Beckmann et al., 2015; Fakhfekh et al., 2021). During the COVID-19 pandemic, Conlon and McGee (2020), Conlon et al. (2020), Ghorbel and Jeribi (2021a), and Jeribi and Fakhfekh (2021) found that, contrary to gold, Bitcoin does not serve as a safe haven for the S&P500. In the same line of results, Corbet et al. (2020), Ghorbel and Jeribi (2021a), and Jeribi and Fakhfekh (2021) found that both Bitcoin and Ethereum do not serve as a safe haven for international stock markets.

Studying the relation between stock markets and oil prices and vice versa began with Hamilton’s (1983) work. Current pieces of Literature in this concern have expanded this research to the financial market and similar results have been identified (Arouri & Rault, 2010; Aydogan, 2017; Ghorbel & Jeribi, 2021b; Narayan & Gupta, 2015; Salisu et al., 2019). The latest instability in the global economy as a result of the COVID-19 outbreak has gained a great deal of interest, including the interaction among movements in oil prices, the economy, and stock markets. Ali et al. (2020) found that, unlike gold, crude oil has the highest volatility due to the COVID-19 pandemic and the regional disputes between oil-producing nations. However, Bakas and Triantafyllou (2020) and Albulescu (2020) suggested that the uncertainty surrounding the 2020 global pandemic had a significant negative effect on the volatility of commodity markets and, in particular, on the crude oil market, with a positive but less significant impact on the yellow metal market.

After the announcement of the COVID-19 pandemic, Sharif et al. (2020) studied the connections between the shock of oil price volatility, the US financial market, geopolitical risk, and economic policy instability. They demonstrate that oil prices were influenced by the COVID-19 pandemic, which can be justified by imposed travel bans and OPEC + agreements. Salisu et al. (2020) developed a panel Vector Autoregressive (pVAR) model to evaluate the behavior of the oil-stock nexus during the COVID-19 pandemic. Their findings show a unidirectional causality from oil price returns to stock returns in the pre-COVID-19 period, while during the 2020 global pandemic, a bi-directional causality between crude oil price and stock returns is reported. Chang et al. (2020) examined the risk spillover effects of the energy sectors in the USA, Europe, and Asia stock markets during three types of global crises. Their results indicated that for low extreme oil returns during the coronavirus crises, investors are more likely to display herding in the stock market. In addition, during the 2020 global pandemic, investors panic so they may unwisely sell their financial assets. In the same line of research, Mazur et al. (2020) found that natural gas, healthcare, software, and food stocks earn high positive returns, whereas, petroleum, entertainment, real estate, and hospitality sectors stocks fall dramatically during the COVID-19 pandemic. Adekoya and Oliyide (2021) and Gorbel and Jeribi (2021b) show that the latter cited pandemic is responsible for risk transmission across commodity and financial markets. Contrary to gold, Ghorbel and Jeribi (2021b) found that Bitcoin cannot be considered as a safe haven during the global pandemic when investing in crude oil.

Based on the above-mentioned studies, this paper extends the related literature and makes an original contribution in identifying the contagion between energy, traditional, and digital financial assets and indicates that Gold is a safe haven when investing in energy assets during the COVID-19 crisis. Bitcoin is considered a safe haven only when investing in the energy index. In this study, firstly, we apply the combined Markov Switching model and GARCH model on the G7 stock indices and the energy index (NYSE energy index), oil (WTI), gas (Henry Hub Natural Gas Spot Price), Bitcoin, and Gold prices. The choice of these countries is explained by the fact that they hold more than 2/3 of the world’s net wealth. In addition, they are the most affected by the COVID-19 (Amar et al., 2021). Secondly, we estimate the MSBEKK-GARCH model. Third, by using the MSDCC-GARCH model, we analyze the dynamic correlations between the studied assets to investigate the contagion effects of the COVID-19 pandemic.

The layout of this paper is as follows. Our econometric methodology is discussed in Sect. 2. Section 3 is devoted to highlighting the relevant data and empirical findings. Finally, Sect. 4 concludes.

2 Empirical methodology

First, referring to Ardia et al. (2018) and Sayed and Auret (2020), we estimate the Markov-Switching GARCH model for all series. Let yt be the daily series log-returns at time t. Following Ardia et al. (2018), we identify the MSGARCH model as:

$$y_{t} ~~/~\left( {s_{t} = k,I_{{t - 1}} } \right)\sim D\left( {0,h_{{k,t}} ,\xi _{k} } \right)$$
(1)

where \(D\left( {0,h_{{k,t}} ,\xi _{k} } \right)\) is a continuous distribution with a zero mean, a time-varying conditional variance\(h_{{k,t}}\), in regime k, and a vector \(\xi _{k}\) of additional shape parameters. The state variable \(s_{t}\) evolves according to a first-order homogeneous Markov chain with a finite number of states k. \(I_{{t - 1}}\) Denotes the information set available up to t − 1. For low and high volume regimes (\(s_{t} = \,1\) and \(s_{t} = 2),~\) we estimate the MSGARCH(1,1) model of Bollerslev (1986) for energy and financial assets:

$$h_{{k,t}}^{{}} = \omega _{k} + \alpha _{k} y_{{t - 1}}^{2} + \beta _{k} h_{{k,t - 1}}$$
(2)

Second, we extend the standard multivariate BEKK-GARCH model of Engle and Kroner (1995) to allow for the presence of regime shifts. We used the bivariate BEKK-MSGARCH for studying the transmission mechanism of shocks (volatility) originating from energy assets (energy index, oil, and gas) to financial assets (G7 stock indices, gold, and cryptocurrency).

According to Haas et al. (2009), we derive the multivariate BEKK-MSGARCH process:

$$\user2{r}_{\user2{t}} = \emptyset + \user2{e}_{\user2{t}}$$
(3)
$$\user2{e}_{{\user2{t},\user2{st}}} = \user2{H}_{{\Delta _{\user2{t}} \user2{t}}}^{{{\raise0.7ex\hbox{$1$} \!\mathord{\left/ {\vphantom {1 2}}\right.\kern-\nulldelimiterspace} \!\lower0.7ex\hbox{$2$}}}} \user2{E}_{\user2{t}} E_{t} /\Omega _{{t - 1}} \to N(0_{{M \times 1}} ,I_{M} )$$

Both the return r and the variance H are made regime-dependent. Let rt be the return matrix at time t, modeled as a constant plus a disturbance term. \(\Phi\) constitutes the constant vector, IM denotes the identity matrix of dimension M. The transition between the successive states is governed by a first-order Markov process {Δt} with finite state space S = {1, 2,..., k} and a primitive (i.e., irreducible and aperiodic) fixed k × k transition probability matrix P.

Where the transition probabilities are given by pij = p \((\Delta _{t} = j/\Delta _{{t - 1}} = i),\) ij = 1, …, k. The regime-dependent covariance matrix H is assumed to follow a Multivariate Markov Switching GARCH (p,q,k) introduced by Bollerslev (1986);

$$h_{{it}} \, = \,\omega _{i} \, + \,\,\sum\limits_{{i = 1}}^{q} {\alpha_{ij}\varepsilon _{{t - i}}^{2} \, + \,\sum\limits_{{j = 1}}^{p} {\beta _{{ij}} h_{{j,t - j}} } ,}\quad j\, = \,1, \ldots ,k$$
(4)

where \(\alpha _{i} = [\alpha _{{i1}} ,...,\alpha _{{ik}} ],\) i = 1,…,q and \(\beta _{i} = [\beta _{{i1}} ,...,\beta _{{ik}} ],\) i = 1,…,p are parameter matrices of appropriate dimension. The number of the independent element of the regime-dependent conditional covariance matrices Hjt is N = M/(M + 1)/2.

To make the application trustworthy, parameter constraints are required. Such a parameterization is provided by the Baba et al. (1990) (BEKK) which specifies the conditional volatility as:

$$H_{{jt}} = \omega _{{0j}}^{*} \omega _{{0j}}^{{*'}} + \sum\limits_{{l = 1}}^{L} {\sum\limits_{{i = 1}}^{q} {\alpha _{{ij,l}}^{*} e_{{t - i}} e_{{t - i}}^{'} \alpha _{{ij,l}}^{{*'}} + } } \sum\limits_{{l = 1}}^{L} {\sum\limits_{{i = 1}}^{q} {\beta _{{ij,l}}^{*} H_{{t - i}} \beta _{{ij,l}}^{{*'}} } } ,\,\,j\, = \,\left\{ {1, \ldots ,k} \right\}$$
(5)

Where \(\omega _{{0j}}^{*}\) are k × k lower triangular matrices of state-dependent coefficients, L is the lag operator. \(\omega _{{0j}}^{*}\),\(\alpha _{{ij}}^{*}\) and \(\beta _{{ij}}^{*}\) are state-dependent matrices. In variance decomposition, first, we estimate \(\alpha _{{ij}}^{*}\) = \(\alpha _{{energy,asset;S_{t} = i}}^{*}\) and \(\beta _{{ij}}^{*}\) = \(\beta _{{energy,asset;S_{t} = i}}^{*}\) from energy index (NY energy index) to financial assets (i = energy index and j = asset) for high volume regime (St = 1) and low volume regime (St = 2). Second, we estimate \(\alpha _{{ij}}^{*}\) = \(\alpha _{{oil,asset;S_{t} = i}}^{*}\) and \(\beta _{{ij}}^{*}\) = \(\beta _{{oil,asset;S_{t} = i}}^{*}\) from Oil (SSE) to financial assets (i = oil and j = asset) for high volume regime (St = 1) and low volume regime (St = 2). Finally, we estimate \(\alpha _{{ij}}^{*}\) = \(\alpha _{{gas,asset;S_{t} = i}}^{*}\) and \(\beta _{{ij}}^{*}\) = \(\beta _{{gas,asset;S_{t} = i}}^{*}\) from the gas spot to financial assets (i = gas and j = asset) for high volume regime (St = 1) and low volume regime (St = 2).

Third, referring to Billio and Caporin (2005), we consider the MSDCC-GARCH model in the following equation:

$$Q_{t} = \left[ {1 - \alpha \left( {s_{t} } \right) - \beta \left( {s_{t} } \right)} \right]\bar{Q}\left( {s_{t} } \right) + \alpha \left( {s_{t} } \right)\eta _{{t - 1}} \eta _{{t - 1}} + \beta \left( {s_{t} } \right)Q_{{t - 1}}$$
(6)
$$\bar{Q}\left( {s_{t} } \right) = \frac{1}{T}\mathop \sum \limits_{{i = 1}}^{T} \eta _{{t - 1}} \eta _{{t - 1}}$$

The Markov chain is governed by the following transition matrix: \(P = \left\{ {p_{{ji}} ~i,j = 1, \ldots S} \right\}\). Where S is the number of regimes, S = 1 for the low regime volatilities and S = 2 for the high regime volatilities.Footnote 1 For a high volume regime (S = 2), we estimate the dynamic conditional correlation between energy assets (energy index, oil, and gas) and financial assets.

3 Data and results

3.1 Data and descriptive statistics

Adjusted closing-price data are relevant to the G7 stock indices (S&P500; FTSE; Nikkei; CAC40; Dax30; S&P/TSX and FTSEMIB), energy index (NYSE energy index), oil (WTI), and gas (Henry Hub Natural Gas Spot Price). Also, we add the Bitcoin and gold price, regarding the period ranging from January 1, 2016, to July 23, 2020, on a daily frequency basis.Footnote 2 We used Oxmetrics software.Footnote 3 Daily returns are defined by \(y_{t} = {\text{ln}}\left( {{\raise0.7ex\hbox{${p_{t} }$} \!\mathord{\left/ {\vphantom {{p_{t} } {p_{{t - 1}} }}}\right.\kern-\nulldelimiterspace} \!\lower0.7ex\hbox{${p_{{t - 1}} }$}}} \right)\), with pt standing for the series respective closing price on day t.

Table 1 presents the statistics summary of the energy index, oil, gas, stock indices, Bitcoin, and gold returns. All assets recorded mean positive returns during this period, except the energy index, gas, and Italian stock index (FTSEMIB). Whereas Gold presents the lowest risk, and Bitcoin presents the highest risk. All asset returns have kurtosis values higher than three and the distribution of returns is negatively and positively skewed for all-time assets, which indicates that all return markets are far from the normal distribution. Therefore, the assumption of Gaussian returns is rejected by the Jarque–Bera test for all assets. The empirical statistics of the Engle (1982) test for conditional heteroscedasticity are significant for all cases suggesting the presence of ARCH effects in returns which justifies our choice of GARCH models.

Table 1 Summary statistics of financials assets

3.2 Empirical results

Figure 1 illustrates the volatilities of the high (S = 2) and low (S = 1) regimes for each asset extracted by the MSGARCH(1,1) model of Eq. 2. Figure 1 plots the time-varying smoothed probabilities for the low- and high-volatility regimes. The stock market, as well as Bitcoin and Gold, are assessed in regime i if the associated smoothed probability P(St = i) is higher than 0.5. Concerning the plots, we note that there are similar regime shift patterns in English, Japanese, Italian, and French stock markets before the 2020 global pandemic. More specifically, all these markets experienced a low volatility regime. However, German and Canadian stock markets registered a high volatility regime. The low– and high–volatility regimes switch among them for the US stock market as well as the energy index, Gas, oil, Gold, and Bitcoin until the end of 2019. During the COVID-19 outbreak, all assets experienced a high volatility regime except for Gold, which can be considered a safe-haven asset for investors during a crisis.

Fig. 1
figure 1

Volatilities of low and high regimes for bitcoin, gold and stock indices

Table 2 presents the estimation results of the Markov Switching-BEKK-GARCH(1,1,2) model, from the energy index to financial assets. The estimation results of the MSGARCH-BEKK(1,1,2) model support the presence of significant volatility spillover from the energy index to the stock indexes, in two regimes (high and low volatility). In other words, the current conditional volatility of each stock index depends not only on its past volatility but also on the past volatility of the energy index. On the contrary, the estimation results of the Markov Switching -BEKK-GARCH(1,1,2) model, from the energy index to Bitcoin show that the GARCH parameters for the two regimes are not significant. So, the current conditional volatility of Bitcoin depends only on its past volatility. In Table 2, the results of estimating ARCH (\(\alpha _{{energy,Bitcoin,s_{t} = 1}}^{{}}\) and \(\alpha _{{enegy,Bitcoin,s_{t} = 2}}^{{}}\)) parameters show evidence of shock transmission effect from energy index to cryptocurrency (Bitcoin). It should be noticed, though, that for the pair energy index—cryptocurrency, (\(\alpha _{{energy,Bitcoin,s_{t} }}^{{}}\)) is significantly positive, in the two regimes. Consequently, the past news about shocks in the energy index positively affects the current conditional volatility of cryptocurrency.

Table 2 Volatility spillover between energy index and financials assets

In Table 2, the results of estimating ARCH (\(\alpha _{{energy,Gold,s_{t} = 1}}^{{}}\) and \(\alpha _{{energy,Gold,s_{t} = 2}}^{{}}\)) parameters show evidence of shock transmission effect from the energy index to Gold returns. It should be noticed, though, that for the pair energy index-Gold, (\(\alpha _{{energy,Gold,s_{t} = 1}}^{{}}\); t = 1 and t = 2) is significantly negative in the two regimes. Consequently, the past news about shocks in the energy index negatively affects the current conditional volatility of Gold. The estimation of GARCH (\(\beta _{{energy,Gold,s_{{t = 1}} }}^{{}}\) and \(\beta _{{energy,Gold,s_{{t = 2}} }}^{{}}\)) parameters is significantly positive between energy and Gold. Consequently, the current conditional volatility of Gold depends not only on its past volatility but also on the past volatility of the energy index.

Table 3 presents the estimation results of the Markov Switching-BEKK-GARCH (1,1,2) model from oil to financial assets. The empirical results from the Markov Switching GARCH-BEKK(1,1,2) model support the shock transmission volatilities from oil to G7 stock indexes in a regime of lower and higher volatilities. In essence, unanticipated events in the oil returns can stimulate higher volatility in the stock indexes. The estimation of the GARCH parameters supports the presence of significant volatility spillover from the Oil to the stock indexes, in two regimes (high and low volatility). In other words, the current conditional volatility of the stock index depends not only on its past volatility but also on the past Oil volatility.

Table 3 Volatility spillover between Oil and financials assets

In Table 3, the results of estimating ARCH (\(\alpha _{{oil,Bitcoin,s_{t} }}^{{}}\); t = 1 and t = 2) parameters show evidence of shock transmission effect between Oil and cryptocurrency (Bitcoin). It should be noticed, though, that for the pair Oil-cryptocurrency, (\(\alpha _{{oil,Bitcoin,s_{t} }}^{{}}\); t = 1 and t = 2) is significantly positive in the two regimes. Consequently, the past news about shocks in Oil positively affects the current conditional volatility of cryptocurrency. The estimation of GARCH (\(\beta _{{oil,Bitcoin,s_{t} }}^{{}}\); t = 1 and t = 2) parameters is significantly positive between Oil and cryptocurrency. Consequently, the current conditional volatility of cryptocurrency depends not only on its past volatility but also on the past volatility of the Oil.

In the two regimes (lower and higher volatilities), the ARCH coefficients (\(\alpha _{{oil,Gold,s_{t} }}^{{}}\); t = 1 and t = 2) document a negative effect of shock transmission volatilities from Oil to Gold. In essence, unanticipated events in the Oil price can stimulate higher volatility in the Gold returns.

By focusing on the GARCH parameter, we found that there is strong evidence of a significant negative transmission of volatility (\(\beta _{{oil,Gold,s_{t} }}^{{}}\); t = 1 and t = 2) from Oil to Gold returns in the two regimes. Consequently, the current conditional volatility of Gold depends not only on its past volatility but also on the past volatility of Oil.

Table 4 presents the estimation results of the Markov Switching-BEKK-GARCH (1,1,2) model from Gas price to financial assets. The estimation results of the MSGARCH-BEKK(1,1,2) model support the presence of significant volatility spillover from the Gas to the stock indexes in the two regimes (high and low volatility). In other words, the current conditional volatility of the stock index depends not only on its past volatility but also on the past volatility of the Gas price. Also, the results of estimating ARCH (\(\alpha _{{gas,asset,s_{t} }}^{{}}\); t = 1 and t = 2) parameters show evidence of shock transmission effect between gas and the stock indexes. It should be noticed, though, that for the pair Gas-stock index, (\(\alpha _{{gas,asset,s_{t} }}^{{}}\); t = 1 and t = 2) is significantly positive in the two regimes. Consequently, the past news about shocks in Gas prices positively affects the current conditional volatility of the stock index.

Table 4 Volatility spillover between gas and financials assets

In Table 4, the estimation results of the Markov Switching-BEKK-GARCH (1,1,2) model, from Gas price to Bitcoin show that the GARCH parameters for the two regimes are significant. The results of estimating ARCH (\(\alpha _{{gas,Bitcoin,s_{t} = 1}}^{{}}\) and \(\alpha _{{gas,Bitcoin,s_{t} = 2}}^{{}}\)) parameters show evidence of shock transmission effect from Gas price to cryptocurrency (Bitcoin). It should be noticed, though, that for the pair Gas-cryptocurrency, (\(\alpha _{{gas,Bitcoin,s_{t} }}^{{}}\)) is significantly positive in the two regimes. Consequently, the past news about shocks in Gas prices positively affects the current conditional volatility of cryptocurrency. The estimation of GARCH (\(\beta _{{gas,Bitcoin,s_{{t = 1}} }}^{{}}\) and \(\beta _{{gas,Bitcoin,s_{{t = 2}} }}^{{}}\)) parameters is significantly positive between Gas and cryptocurrency. Consequently, the current conditional volatility of cryptocurrency depends not only on its past volatility but also on the past volatility of the Gas prices.

In Table 4, the results of estimating ARCH (\(\alpha _{{gas,Gold,s_{t} = 1}}^{{}}\) and \(\alpha _{{gas,Gold,s_{t} = 2}}^{{}}\)) parameters show evidence of shock transmission effect from Gas to Gold returns. It should be noticed, though, that for the pair Gas-Gold, (\(\alpha _{{gas,Gold,s_{t} }}^{{}}\); t = 1 and t = 2) is significantly negative in the two regimes. Consequently, the past news about shocks in Gas prices negatively affects the current conditional volatility of Gold. The estimation of GARCH (\(\beta _{{gas,Gold,s_{{t = 1}} }}^{{}}\) and \(\beta _{{gas,Gold,s_{{t = 2}} }}^{{}}\)) parameters is significantly positive between Gas and Gold. Consequently, the current conditional volatility of Gold depends not only on its past volatility but also on the past volatility of the Gas prices.

Figures 2, 3, and 4 present the dynamic correlation between the G7 stock market indices and the NYSE energy index, WTI, and Natural Gas Spot, respectively, for the high regime. Figure 2 shows high and positive dynamic relationships between the NYSE energy index and the S&P500 as well as CAC 40. However, we observe a low and positive correlation between DAX30 and the NYSE energy index. The correlation became negative for the British, Canadian, and Italian stock market indices during the 2020 global pandemic. This can be explained by the fact that investors’ and policy uncertainties have significantly raised. In addition, the global financial cycle is immensely altered. The dynamic relationship between the NYSE energy index and the yellow metal is always negative. As for Bitcoin, the dynamic correlation with the NYSE energy index switch between negative and positive values and reaches its lowest level during the COVID-19 pandemic. Gold and Bitcoin can be considered as safe havens during the global pandemic when investing in the NYSE energy stocks. This result is consistent with that of Gaspareniene et al. (2018) and Robiyanto et al. (2019). However, it is inconsistent with that of Conlon and McGee (2020) and Corbet et al. (2020). Figures 3 and 4 show that, the dynamic relationships between Oil and Gas prices on the one hand, and G7 stock market prices, on the other hand, switch between low and high levels with different levels for each stock market. It proves the volatility spillover from energy assets to financial assets. Contrary to the Japanese stock market, the US and Canadian stock markets are the most correlated with the WTI. In addition, the US, Italian, and Canadian stock markets are the most correlated with Gas prices. These results are consistent with Chang et al. (2020) and Adekoya and Oliyide (2021). Few important findings are evidenced in Table 5. In fact, our findings indicate also that the correlations for the couples WTI-Gold, WTI-Bitcoin, and Gas-Gold as well as Gas-Bitcoin are symmetric. During the COVID-19 pandemic, the correction between crude oil and Gas on the one hand, and the yellow metal, on the other hand, reached its lowest level. However, the correction between crude oil as well as Gas and Bitcoin reached its highest level. Contrary to Gold, Bitcoin cannot be considered as a safe haven during the COVID-19 pandemic for investors that purchase crude oil and Gas as well as energy producers. This result confirms that of Conlon and McGee (2020) and Corbet et al. (2020).

Fig. 2
figure 2

Dynamic correlation between NYSE energy index and financial assets for high volatility regime

Fig. 3
figure 3

Dynamic correlation between Oil and financial assets for high volatility regime

Fig. 4
figure 4

Dynamic correlation between Gas and financial assets for high volatility regime

Table 5 Summary of the safe haven and spillover analysis

4 Conclusion

This paper performs a general Markov–Switching GARCH model for modeling structural breaks for the energy sector index, Oil, Gas, cryptocurrency, G7 stock indices, and Gold. We show that the energy assets, Bitcoin, gold, and stock indices exhibit regime changes in their volatilities dynamics. The estimation results of the MSGARCH-BEKK(1,1,2) model support the presence of significant volatility spillover from the energy index to the stock indexes in two regimes, especially in high volatility. In fact, the estimation results of the Markov Switching-BEKK-GARCH(1,1,2) model, show that the past news about shocks in the energy index positively affects the current conditional volatility of stock indices and cryptocurrency, however, it negatively affects the current conditional volatility of Gold. In other words, the current conditional volatility of each stock index and Gold depends not only on its past volatility but also on the past volatility of the energy index.

The estimation results of the MSBEKK-GARCH prove a significant effect of choc from Oil to financial assets (G7 stock indexes, Gold, and Bitcoin), on the one hand, and from Gas to financial assets, on the other hand. In other words, the stock indices and Gold depend not only on their past volatility but also on past volatilities of Oil and Gas. When focusing on the time-varying correlation between the G7 stock market indices and the energy index, WTI and Natural Gas Spot, respectively, we found high and positive dynamic relationships between the energy index and the S&P 500 as well as CAC 40. The correlation between the energy index and other G7 indices switches between low and high levels. In addition, we argue that Gold and Bitcoin can be considered as safe-haven assets during the global pandemic when investing in the stock indices. We show also that the dynamic relationships between Oil and G7 stock markets, on the one hand, and Gas price and G7 stock markets, on the other hand, switch between low and high levels with different levels for each stock index. Contrary to Gold, Bitcoin cannot be considered as a safe haven during the COVID-19 pandemic when investing in energy assets (crude oil and Gas).

The empirical results have implications for both institutional and retail investors, risk management, as well as for economic and financial policies, highlighting the prevalence of tail behavior and persistent asymmetric nature of risk spillovers. Thus, we draw two main managerial implications. First, the G7 financial market regulators should have a solution to immediately support the market when fear is overwhelming. In this case, market crash risk can be managed in extreme cases. Second, both investors and regulators should implement hedging or safe haven strategies that have been largely ignored by many analysts and policymakers especially during the SARS crisis, not least by those who have been involved with energy issues. From an investing point of view, Oil prices can affect stock prices through their effect on potential corporate cash flows and indirectly through the interest rate used to discount future cash flows. However, macroeconomic variables adjust based on the oil-stock relationship aids policy-making by helping to understand the economy’s sensitivity to both internal (due-to-stock) and external (due-to-oil) shocks. Thus, the behavior of the Oil and stock markets, especially during the pandemic, is therefore essential to policy-making and to achieving better macroeconomic outcomes.