# 1. Introduction

Changes and co-movement in gold and oil prices have always been the subject of worldwide attention. According to Shahbaz et al. (2017), there has been a positive price correlation between gold and oil prices more than 80% of the time in the past 50 years. The economic recession caused by the COVID-19 pandemic and unsuccessful trade negotiations between Russia and the Organization of the Petroleum Exporting Countries have led to a dramatic drop in West Texas Intermediate (WTI) crude oil prices (Atri et al., 2021). Uncertainty has always driven investors to seek safer assets, particularly gold. Many studies have demonstrated the safe haven property of gold (Baur & McDermott, 2010; Dutta et al., 2020; Reboredo, 2013; Salisu et al., 2021). As seen in Figure 1, the price of gold has witnessed a mild rise since March 2020, despite the decline in mid-March.

Figure 1.Daily oil and gold price between 2018 and 2021

This figure shows the time-series daily prices of gold and West Texas Intermediate (WTI) crude oil prices for the period April 3, 2018 to June 25, 2021. The data are sourced from the US Energy Information Administration and world gold council websites, see https://www.gold.org/cn/goldhub/data/gold-prices and https://www.eia.gov/.

Only a few studies have investigated the nexus between gold and oil prices, with the majority examining the mean and volatility spillovers between the returns of these two assets (Bampinas & Panagiotidis, 2015; Chen & Qu, 2019; Gharib et al., 2021). The spillover effect refers to the phenomenon where a large information shock increases the return correlation not just in the market in question, but also in other markets. Mean spillovers are also shown to be time varying. According to Yaya et al. (2016), the spillover effects between gold and oil prices disappeared after 2008, with only one-way spillover from gold to oil prices remaining. Using daily exchange-traded fund (ETF) prices of gold, silver, oil and global equity, Lau et al. (2017) show that the previous week’s information contributes to the transmission of oil and gold return spillovers between the gold and oil markets. Studies have also found that the volatility spillover effect differs over time. For instance, Kang et al. (2017) find more pronounced volatility spillovers across commodity futures markets, including oil and gold, in post-crisis periods, including the global financial crisis (from November 23, 2007, to May 22, 2009) and two European debt crises (from November 5, 2009, to April 22, 2010, and from May 2, 2011, to May 30, 2012).

Despite the above-mentioned works, relatively little is known about the asymmetry and direction of spillover effects between the oil and gold markets. Ross (1989) argues that asset returns would induce different volatility patterns and spillover effects, since the information innovation and processing times vary across markets. Asymmetric reactions in information shocks between pairs of assets are more common today, and understanding their possible implications will be useful to both investors and policy makers (Salisu & Oloko, 2015). Motivated by this research, in this paper we hypothesize that the oil and gold markets are asymmetrically related by their shocks on each other. Specifically, this study examines the asymmetric spillover effects in these two important commodity markets during the COVID-19 pandemic.

We first employ a vector autoregressive (VAR) model to test the mean spillover effects. The results suggest that the mean spillover effect is one way and there is no immediate return spillover effect from the gold market to the oil market. We use the Baba–Engle–Kraft–Kroner asymmetric generalized autoregressive conditional heteroskedasticity (BEKK-AGARCH) model to estimate the volatility spillover effects and find that a negative information shock in the oil market has a large impact on gold return volatility compared to a positive shock, suggesting that the effect is asymmetric. We show that the spillover effect from the oil to the gold market is stronger during the COVID-19 pandemic. Our study therefore contributes to the literature on gold–oil spillover effects by exploring not only the asymmetric nature of the relation between the two markets, but also how the pandemic has intensified the relation between the markets. This work thus adds to the broad and growing number of studies on COVID-19 (Narayan, 2021; Padhan & Prabheesh, 2021).

# 2. Methodology

First, a conditional $$\text{VAR}(p)$$ model is proposed to examine the mean spillover effects between gold and oil price returns $$(r_{i,t})$$, with the following specification:

$R_{t} = \alpha + \sum_{k = 1}^{p}{\beta_{k}R_{t - k}} + \varepsilon_{t}\ \tag{1}$

where
$$R_{t} = {(r_{oil,t},r_{gold,t})}^{'}$$, with $$r_{oil,t}$$ and $$r_{gold,t}$$ being the returns on the crude oil and gold markets, respectively, at time $$t$$;
$$\beta_{k}$$ is a $$(2 \times 2)$$ matrix of coefficients of the form $$\beta_{k} = \begin{pmatrix} \beta_{11} & \beta_{12} \\ \beta_{21} & \beta_{22} \\ \end{pmatrix}_{k};$$ and
$$\varepsilon_{t} = {(\varepsilon_{oil,t},\varepsilon_{gold,t})}^{'}$$, with $$\varepsilon_{oil,t}$$ and $$\varepsilon_{gold,t}$$ being the error terms from the mean equations for the oil and gold market returns, respectively, which has a conditional covariance matrix $$H_{t}$$, given the available information set $$I_{t - 1}$$.

When $$\beta_{\text{ij}}, i \neq j$$, is significantly different from zero, we can expect market $$i$$ to have a mean spillover effect on market $$j$$. The best lagged term $$p$$ is chosen based on the Akaike information, final prediction error, and Bayesian information criteria.

In terms of the conditional variance–covariance equations, we utilize the unrestricted BEKK-AGARCH(1,1) model, following Engle and Kroner (1995) and Salisu and Oloko (2015), to explore the volatility spillover effect. This AGARCH model allows for the interaction among the conditional variances and covariances across several time series. In addition, it incorporates the asymmetric effect coefficients to help us identify the different effects of positive and negative shocks. The conditional covariance–variance equation is expressed as follows:

\begin{align} H_{t} =&\ M^{'}M + A^{'}{\varepsilon_{t - 1}\varepsilon}_{t - 1}^{'}A + B^{'}H_{t - 1}B \\ &+ C^{'}I_{t - 1}{\varepsilon_{t - 1}\varepsilon}_{t - 1}^{'}C \end{align} \tag{2}

where M, A, B, and C are square matrices; $$H_{t}$$ is the conditional variance–covariance matrix, which represents market volatility; and $$I_{t - 1}$$ is a function of the error term. Specifically, we have
$$H_{t} = \begin{bmatrix} h_{11,t} & h_{12,t} \\ h_{21,t} & h_{22,t} \\ \end{bmatrix}$$, where $$h_{ii,t}$$ is the conditional variance of market $$i$$ at time $$t$$ and the off-diagonal elements $$h_{ij,t}$$ ($$i \neq j)$$ represent the conditional covariance between market $$i$$ and $$j$$ at time $$t$$;
$$M = \begin{bmatrix} m_{11} & 0 \\ m_{21} & m_{22} \\ \end{bmatrix}$$, which is a lower triangular matrix to ensure the positive definiteness of $$H_{t}$$ (Bekiros, 2014);
$$A = \begin{bmatrix} a_{11} & a_{12} \\ a_{21} & a_{22} \\ \end{bmatrix}$$, whose elements are the coefficients of the ARCH term, showing the effects of the volatility in the market in question and the volatility spillover from the other market on the conditional volatility of the specified market;
$$B = \begin{bmatrix} b_{11} & b_{12} \\ b_{21} & b_{22} \\ \end{bmatrix}$$, whose elements represent the coefficients of the GARCH term, showing the effect of past volatility in the market in question and past volatility from the other market on the conditional volatility of the specified market;
$$C = \begin{bmatrix} c_{11} & c_{12} \\ c_{21} & c_{22} \\ \end{bmatrix}$$, whose elements indicate the coefficients of the asymmetric effect, where we do not expect the positive and negative shocks to have the same effects on the conditional volatility; and $$I_{t-1}=\left[\begin{array}{cc} I_{\text {oil }, t-1} & 0 \\ 0 & I_{\text {gold,t-1 }} \end{array}\right]$$ where $$I_{\text {oil, } t-1}$$ and $$I_{i, t-1}=\left\{\begin{array}{l} 0, \varepsilon_{i, t}>0 \\ 1, \varepsilon_{i, t} \leq 0 \end{array}\right.$$ for the $$i$$ th market at time $$t$$.

Note that the interpretation of the estimated coefficient based on the variance and covariance matrix is not easy. Therefore, the conditional variance for each equation can be expanded as follows:

\begin{align} h_{11, t}=&\ m_{11}^{2}+a_{11}^{2} \varepsilon_{\text {oil }, t-1}^{2}\\ &+a_{21}^{2} \varepsilon_{\text {gold }, t-1}^{2}+2 a_{11} a_{21} \varepsilon_{o i l, t-1} \varepsilon_{\text {gold }, t-1}\\ &+c_{11}^{2} \varepsilon_{\text {oil }, t-1}^{2} I_{\text {oil }, t-1} +c_{21}^{2} \varepsilon_{g o l d, t-1}^{2} I_{o i l, t-1}\\ &+2 c_{11} c_{21} \varepsilon_{o i l, t-1} \varepsilon_{g o l d, t-1} I_{o i l, t-1}+b_{11}^{2} h_{11, t-1}\\ &+b_{21}^{2} h_{22, t-1}+2 b_{11} a_{21} h_{21, t-1} \end{align} \tag{3}

\begin{align} h_{22, \mathrm{t}}=&\ m_{22}^{2}+a_{12}^{2} \varepsilon_{\text {oil }, \mathrm{t}-1}^{2}\\ &+a_{22}^{2} \varepsilon_{\text {gold }, \mathrm{t}-1}^{2}+2 a_{12} a_{22} \varepsilon_{\text {oil }, \mathrm{t}-1} \varepsilon_{\text {gold }, \mathrm{t}-1}\\ &+c_{12}^{2} \varepsilon_{\text {oil }, \mathrm{t}-1}^{2} I_{\text {oil }, \mathrm{t}-1} +c_{22}^{2} \varepsilon_{\text {gold }, t-1}^{2} I_{\text {oil }, t-1}\\ &+2 c_{12} c_{22} \varepsilon_{\text {oil } t-1} \varepsilon_{\text {gold }, t-1} I_{o i l, t-1}+b_{12}^{2} h_{11, t-1}\\ &+b_{22}^{2} h_{22, t-1}+2 b_{12} a_{22} h_{21, t-1}\\ \end{align} \tag{4}

\begin{align} h_{21, t}=&\ m_{21} m_{22}+a_{11} a_{22} \varepsilon_{\text {oil }, t-1}^{2} +a_{21} a_{22} \varepsilon_{\text {gold }, t-1}^{2}\\ &+\left(a_{21} a_{12}+a_{11} a_{22}\right) \varepsilon_{\text {oil }, t-1} \varepsilon_{\text {gold }, t-1}\\ &+c_{11} c_{22} \varepsilon_{\text {oil }, t-1}^{2} I_{\text {oil } t-1}\\ &+c_{21} c_{22} \varepsilon_{\text {gold }, t-1}^{2} I_{\text {oil }, t-1}\\ &+\left(c_{21} c_{12}+c_{11} c_{22}\right) \varepsilon_{\text {oil } t-1} \varepsilon_{\text {gold }, t-1} I_{\text {oil, } t-1}\\ &+b_{11} b_{22} h_{11, t-1} +b_{21} b_{22} h_{22, t-1}\\ &+\left(b_{21} b_{12}+b_{11} b_{22}\right) h_{21, t-1} \end{align} \tag{5}

To be specific, Eqs. (3) and (4) indicate how shocks and volatility are transmitted within the market in question over time, while Eq. (5) indicates the volatility spillover across different markets.

# 3. Data

The dataset in our study covers two variables, crude oil prices and gold prices. NYMEX WTI is used to proxy for international crude oil prices, because of its high liquidity and pricing transparency. We employ the price for XAUUSD contracts to represent gold prices, since its trading volume is the largest in the world. The time period is from June 27, 2018, to June 25, 2021. The full sample is further divided into two subsamples, with the subperiod prior to the COVID-19 pandemic (from June 27, 2018, to December 30, 2019) and the COVID-19 subperiod (from December 30, 2019, to June 25, 2021) subperiod. All daily trading data are collected from the websites of the US Energy Information Administration[1] and the World Gold Council.[2] Both the oil and gold prices are quoted in US dollars and the return series $$(r_{i,t})$$ are calculated as follows:

\begin{align} r_{i,t} = lnP_{i,t} - lnP_{i,t - 1} \tag{6} \end{align}

where $$P_{i,t}$$ is the ith (gold or oil) market price at time $$t$$, and $$r_{i,t}$$ is the $$i$$th market return at time $$t$$. After eliminating missing values, we obtain a final full sample of 751 observations.

# 4. Empirical results

Table 1 shows the summary statistics for the two variables. The oil returns are more volatile than the gold returns, both during the full sample period and in the subperiods. According to the kurtosis statistic, the oil returns have fatter tails than the gold returns, and oil is thus risker. The Jarque–Bera results suggest that neither return follows a normal distribution, which is consistent with the kurtosis and skewness statistics. The Ljung–Box Q-statistics at both the 10th and 20th lags show that both return series display serial correlation at the 5% level of significance. Therefore, both returns exhibit conditional heteroscedasticity and ARCH effects. Additionally, augmented Dickey–Fuller and Phillips–Perron statistics suggest that all the series are stationary at the 1% level of significance. In conclusion, the VAR-BEKK-AGARCH(1,1) model can be applied in our study.

Table 1.Descriptive statistics of oil and gold price returns
Variables N Mean SD Min Max Skewness Kurtosis Jarque-Bera ADF PP test
A: full-sample
07/27/2018
06/25/2021
r oil 751 0.00191 0.029 -0.281 0.426 0.995 36.897 190000*** -65.32344*** -65.270***
r gold 751 0.00232 0.012 -0.096 0.068 -0.379 8.151 4394*** -62.749*** -62.764***
B: prior-COVID-19
07/27/2018
12/30/2019
r oil 375 0.00223 0.093 -1.000 0.999 0.204 74.046 79000*** -22.887*** -23.962***
r gold 375 0.00143 0.012 -0.053 0.051 -0.613 6.83 252*** -19.854*** -19.857***
C: COVID-19
12/30/2019
06/25/2021
r oil 376 -0.00029 0.024 -0.128 0.164 0.135 7.685 3229*** -61.185*** -61.224**
r gold 376 0.00112 0.012 -0.096 0.068 -0.353 8.312 4207** -59.471*** -59.474***

This table presents descriptive statistics for the gold and oil returns during the full period (Panel A) and two subperiods (Panel B and Panel C). $$r_{\text{oil}}$$ and $$r_{\text{gold}}$$ represent the market logarithm return for oil and gold prices respectively. ***, ** and * represent significance at the 1%, 5% and 10% levels, respectively.

## 4.1 Mean equation analysis

The estimation for the VAR(p) model requires setting the lag length p, for which we use the Akaike information and Bayesian information criteria. The optimal lag turns out to be two. Table 2 presents the results for the mean equation, which illustrate the shock and return spillover effects. It is evident that both the gold and oil returns are negatively affected by their own past shocks at the 5% level of significance, based on the negative and significant estimates of $$\beta_{11,1}$$, $$\beta_{11,2}$$, $$\beta_{22,1}$$, and $$\beta_{22,2}$$. As for the cross-market spillover effect, the oil market has a positive and significant return spillover effect on the gold market for the three periods, implied by the coefficients of $$\beta_{12,1}$$ and $$\beta_{12,2}$$. This spillover effect is intensified during the COVID-19 pandemic, since the coefficients of $$\beta_{12,1}$$ and $$\beta_{12,2}$$ are larger during the pandemic period. However, our study reveals that the oil market returns are not affected by past shocks to gold market returns, as evidenced by the nonsignificant coefficients of $$\beta_{21,1}$$ and $$\beta_{21,2}$$. In other words, there is no immediate return spillover effect from the gold market to the oil market. Our result is consistent with that of Salisu et al. (2021).

Table 2.Estimation result for VAP(p) model
A: full-sample
07/27/2018
06/25/2021
B: prior-COVID-19
07/27/2018
12/30/2019
C: COVID-19
12/30/2019
06/25/2021
α -0.00023 -0.00041 0.0234
β 11,1 -0.00435** -0.00231*** -0.00591***
β 12,1 0.00543*** 0.00031** 0.04398***
β 21,1 -0.03214 -0.00087 -0.06983
β 22,1 -0.28319** -0.01248*** -0.87745***
β 11,2 -0.15323 -0.24103 -0.09248
β 12,2 0.00658*** 0.00194*** 0.07815***
β 21,2 -0.00983 -0.02047 -0.02941
β 22,2 -0.00865** -0.00915** -0.01302***

This table presents the result for VAR (p) model. Based on criteria, the best lagged term is 2. Therefore, there are two estimation matrices for $$\beta_{k}$$. ***, ** and * represent significance at the 1%, 5% and 10% levels, respectively.

## 4.2 Conditional variance–covariance equation analysis

Table 3 presents the BEKK-AGARCH(1,1) estimation results, evaluating the shock and volatility spillover effects between the markets. As mentioned above, the diagonal coefficients in the matrices A and B capture the market’s own ARCH and GARCH effects. The significant coefficients of $$a_{\text{ii}}$$ and $$b_{\text{ii}}$$ across all three periods indicate that the volatilities of both the gold and oil markets are affected by their own past shocks and volatility. Notably, the ARCH and GARCH effects are magnified during the COVID-19 period for both markets.

Table 3.Estimation result for BEKK-AGARCH (1,1) model
A: full-sample
07/27/2018
06/25/2021
B: prior-COVID-19
07/27/2018
12/30/2019
C: COVID-19
12/30/2019
06/25/2021
m 11 0.1571** 0.1968 0.1980
m 21 0.0132 0.1791** 0.3241**
m 22 0.1821 -1.20E-07 3.90E-07
a 11 -0.0023*** -0.7152*** 0.1352***
a 12 0.0018** -0.02991*** 0.9017***
a 21 -0.1459 0.2549 -0.3819
a 22 -0.1483*** 0.1990*** -0.3196***
b 11 -0.0042*** 0.0049** 0.2415***
b 12 -0.0037*** 0.0049*** 0.2373***
b 21 -0.7831 -0.7329 -2.3272
b 22 -0.7851*** -0.8341** -2.255***
c 11 0.5910*** 0.0732*** 0.9871***
c 12 -0.0061*** -0.0022** 0.1133***
c 21 0.3013 0.5762 0.4539
c 22 0.9130*** 0.7282*** 0.5243***
Model diagnostics
AIC 6.982 7.193 4.948
Log-L -62734.139 -10731.193 -43612.573

The table represents the estimation result for BEKK-AGARCH (1,1) model, ***, ** and * represent significance at the 1%, 5% and 10% levels, respectively.

The off-diagonal elements of the matrices $$A$$ and $$\text{B\ }$$measure the shock and volatility spillover effects across the markets. The current conditional volatility of the gold market will be affected by both the conditional variance $$(b_{12}^{2})$$ and lagged shocks $$(a_{12}^{2})$$ of the oil market before COVID-19. This result is consistent with the work of (Ewing & Malik, 2013). One of the possible reasons is that the volatility spillover effect from the oil market to the gold market increases during the COVID-19 pandemic, given the larger positive and significant coefficient of$$\ b_{12}^{2}$$. However, no volatility or shock spillover effects from the gold market to the oil market are observed, given the nonsignificant coefficients of $$b_{21}^{2}$$ and $$a_{21}^{2}$$. Therefore, the volatility spillover effect is from the oil to the gold market and not from the gold to the oil market, which means that shocks to the oil market are more important for the gold market than shocks to the gold market are for the oil market. This finding could be due to the crude oil market being the largest commodity market in the world, and its volatility thus has a greater impact on the gold market.

Finally, the elements in matrix C represent asymmetric spillover effects in the market. A negative shock in the market results in larger volatility in the market in question compared to positive shocks, implying that the shocks are asymmetric. This result is widely observed in many financial asset markets. These effects are more pronounced during the COVID-19 period. In addition, a negative shock in the oil market leads to greater volatility in the gold market than a positive shock, since $$c_{12}^{2}$$ is significantly different from zero.

Our model diagnostic test results in Table 3 suggest the robustness of the results, based on the Akaike information criterion, and the log-likelihood statistics suggest that our estimated BEKK-AGARCH(1,1) model is well specified, reflecting the robustness of our model.

# 5. Conclusion

This study investigates the spillover effects between the international gold and oil markets. Using daily return data from June 27, 2018, to June 25, 2021, we employ a VAR(p)-BEKK-AGARCH(1,1) model to test whether there is an asymmetric spillover effect between the gold and oil markets and explore how this effect changes during the COVID-19 pandemic. Our results show no bidirectional return or volatility spillover effects between the gold and oil markets during the COVID-19 pandemic. Only a shock or volatility in the oil market leads to conditional volatility in the gold market, whereas the gold market’s volatility has no significant effects on the oil market. In addition, the spillover effect from the oil market to the gold market is amplified during the COVID-19 pandemic.

Second, both the oil and gold markets exhibit asymmetric spillover effects, such that a negative information shock has a large impact on market volatility. Moreover, a negative information shock in the oil market leads to greater volatility in the gold market than a positive shock. Our research analyzes how the COVID-19 pandemic influences the spillover effects between the gold and oil markets. Our paper also examines the asymmetric spillover effects between the gold and oil markets under the BEKK model. Further research could be conducted in terms of the asymmetric spillover effects among more asset markets, such as the exchange rate and commodity markets.