key: cord-1022126-l1as6oeg authors: Park, Beum-Jo title: The COVID-19 pandemic, volatility, and trading behavior in the bitcoin futures market() date: 2021-08-26 journal: Res Int Bus Finance DOI: 10.1016/j.ribaf.2021.101519 sha: 59288d803f26478476e438b84403abc65cc26aca doc_id: 1022126 cord_uid: l1as6oeg This paper contributes to the literature on the coronavirus (COVID-19) pandemic impacts on the Bitcoin futures (BTCF) market and to the ongoing consideration of the dynamic relationship between volatility (or returns) and trading behavior variables, such as volume and open interest as a proxy for belief dispersion. This paper focuses on the role of the unprecedented market stress induced by the COVID-19 pandemic in the interrelations among the variables. Accordingly, this paper proposes a structural change (SC)-VAR-MGARCH model and finds the COVID-19 pandemic has initiated a significant regime change. Furthermore, the relationship between the variables in the pre-pandemic regime is notably unclear, whereas an increase in belief dispersion in the pandemic regime due to market stress reduces BTCF returns but raises trading volume and volatility evidently. The outcomes in the pandemic regime are remarkably consistent with the difference of opinions model, though existing evidence on the dynamic relations is ambiguous. Moreover, the outcomes support our hypothesis that, in addition to information flows, market stress causing traders’ behavioral biases should be considered as one of the crucial factors of tremendous price variability. The current coronavirus pandemic forces us to live in a socio-economic situation that is rarely experienced in the past. 1 Economic collapse and widespread financial market stress will be featured in the literature in the following years. Market stress is apparent in other financial markets 2 , including the crypto market, though sharp price decline and tremendous equity market volatility have garnered considerable attention. COVID-19 has spread to countries outside of Mainland China. As a result, many investors switched from holding risky assets to perceived safe-haven assets due to mounting worries about the COVID-19 pandemic. Thus, worldwide equity markets that trade representative risky assets reported their largest single-week declines in the final week of February 2020 since the 2008 financial crisis. Moreover, the US equity market lost 12 percent of its value under enormous pressure. ☆ I would like to express my appreciation to seminar participants at the Financial Economic Meeting for their valuable comments. Certain COVID-19 pandemic impacts on financial markets (e.g., Ashraf, 2020; Shehzad et al., 2020; Harjoto et al., 2021) are concerning from the perspective of investment and risk management. Insights into futures price volatility associated with trading behavior under market stress can adversely impact financial market participants, such as speculators, hedgers, and arbitrageurs, using index futures contracts. Thus, this paper focuses on analyzing how the pandemic-induced market stress has affected the dynamic relationship between volatility and trading behavior variables. These variables include trading volume and open interest as a proxy for dispersion of traders' beliefs on asset prices following Bessembinder et al. (1996) . In the regime of market stress, traders' biases increase with the belief that updating overreactions to news and traders with relatively different beliefs are more dominant (Basak and Atmaz, 2018) . Therefore, market stress leading to the increase in belief dispersion significantly impacts the relationship between volatility and trading behavior. In financial literature, their relationship has received significant attention from economists and market participants. Among the theoretical studies, Clark (1973) explains the relationship between volatility and trading volume through the mixture of distribution hypothesis (MDH). His findings suggest that volatility and trading volume are positively related due to its dependence on a common latent mixing variable, namely, rate of information arrival (e.g., Epps and Epps, 1976; Tauchen and Pitts, 1983; Fleming et al., 2006) . The sequential information arrival model (SIAM) is another theory describing the relationship and generally considered (e.g., Copeland, 1976; Jennings et al., 1981) . SIAM assumes that the sequential arrival of new information to the market generates trading volume and price movements. As a result, volatility and volume measured over the period of full response to information arrival exhibit a positive relationship. The primary implication of this model is that volatility can be potentially forecasted for the future price based on trading volume. In addition to trading volume 3 , the information role of open interest has gained much scholarly attention because new information inflow to the market may expand the dispersion of traders' beliefs, increasing trading volume (Kim and Verrecchia, 1991) . The theory of information-flow paradigm represented by MDH and SIAM illustrates a positive relationship between volatility, trading volume, and open interest. However, the empirical evidence is ambiguous or conflicting. For instance, various studies reveal ambiguous and negative relationships between these variables. Galati (2000) ; Liesenfeld (2001) ; Amatyakul (2010) , and Park (2011) among others present the evidence of a negative relationship between volatility and trading volume. Bessembinder and Seguin (1993) ; Diether et al. (2002) , and Yu (2011) among others also identify that the relationship of other variables can be negative. Other researchers argue that the relationship changes over time and relies on strategies and certain characteristics (e.g., Admati and Pfleiderer, 1988; Doukas et al., 2006; Avramov et al., 2009; Czudaj, 2019) . Admati and Pfleiderer (1988) assert that the relationship between volume and volatility depends on a strategic trading behavior of informed and liquidity traders. Czudaj (2019) emphasizes that the relationship between trading volume, volatility, and open interest in agricultural futures markets varies substantially over time. Volatility is driven by previous periods' trading volume and open interest, whereas the reversed relationship from lagged volatility to trading volume and open interest is limited to certain periods. In sum, the mixed empirical evidence indicates that information flows are not the only potential factor for price variability but also other crucial factors that should not be overlooked. This study aims to gain further insights into the feasible reasons for the drastic increase in volatility in financial markets during the COVID-19 pandemic. This paper concentrates on pandemic-induced market stress, its impact on trading behavior variables, such as volume and open interest, and the relationships of these variables with return volatility. Accordingly, this paper considers the Bitcoin futures (BTCF) market for empirical analysis 4 , proposes a SC-VAR-MGARCH model incorporating a structural changing vector autoregressive (Sims, 1980) with multivariate GARCH effects (Bollerslev, 1986) and introduces a two-step estimation method. The multivariate model is set up to consider certain properties of the BTCF market. Risks of the COVID-19 pandemic prevail in this market, i.e., volatility clustering in trading behavior variables and returns, dynamic relations, and expected structural break. Hence, this paper contributes to the literature in several ways. First, to the best of the author's knowledge, this study is the first to explore the impacts of pandemic-induced market stress on the dynamic relationship between the variables in futures markets. This study is significant in the literature as the ongoing COVID-19 pandemic raises concerns about future market stress under greater economic costs and uncertainty. Second, this study discovers volatility dynamics in BTCF returns relative to trading behavior under extreme market stress. Expanding the understanding of volatility sources in the BTCF market could lead to superior trading and portfolio management strategies. Third, previous studies (e.g., Thies and Molnár, 2018; Ardia et al., 2019; Koutmos, 2020) have found that Bitcoin exhibits regime-switching behavior and pointed out that ignoring regime changes in volatility dynamics can lead to biased estimation results in volatility models. Similarly, this study investigates whether the COVID-19 pandemic resulted in further regime changes. For this investigation, this paper suggests a suitable test method, proposed by Ross (2015) , as it considers the model 3 Open interest is sometimes confused with trading volume, but the two terms refer to different measures. That is, trading volume is widely known to represent the number of contracts traded in a specific period, whereas open interest is the number of outstanding traded contracts at a time point. 4 For our empirical analysis on the effects of the COVID-19, the BTCF market is well suited because of several reasons. First, despite the short COVID-19 period, it has a relatively sufficient accumulation of data compared to other asset markets due to the nature of trading worldwide around the clock (24 hours/7 days) in multiple online exchanges. Second, trading behavior such as volume and open interest are closely linked to demand for speculation and hedging of futures contracts (Bessembinder and Seguin, 1993; Lucia and Pardo, 2010) that are well represented in the BTCF market (Jena et al., 2018) . Third, empirical studies (e.g., Fassas et al., 2020) show that the BTCF market has a greater influence on the value of the Bitcoin than the BTC market, implying the futures dominance in the case of Bitcoin. Fourth, crypto markets including the BTCF market do not involve financial institution intermediaries and are linked to the peer-to-peer transactions which reduce transaction costs (Fang et al., 2021) . This means that crypto markets are liable to attract speculative interest and reveal traders' behavioral biases (e.g., Kallinterakis and Wang, 2019) . Fifth, crypto markets have gained both market value and academic attention. specification. The rest of this paper is outlined as follows. Section 2 describes data set and preliminary statistics. Section 3 introduces an empirical methodology, including the SC-VAR-MGARCH model and a two-step estimation method. Section 4 runs a test for detecting a structural change caused by the COVID-19 pandemic and provides empirical evidence of the interrelations among volatility, trading volume, and belief dispersion. Then, the feasibility of the econometric model is derived. Finally, Section 5 provides conclusions and suggestions for further work. This paper collected data from the Bitfinex exchange using the exchange API from October 29, 2019 to July 5, 2020. 5 The Bitfinex exchange records individual trades for BTCF contracts in real time. Thus, trade-level data match our research purpose to analyze the impact of pandemic-induced market stress on trading behavior and dynamic volatility. In addition, the exchange has the largest volumes of BTC/BTCF-USD exchanges during the selected period 6 and provides data of continuous futures price and trade series. The variables in the daily data are closing BTCF price, trading volume, and open interest. To avoid problems arising from the non-stationary state usually observed in asset prices, we take the natural logarithmic differences between two successive trading days. Let p t be the price at time t, then the log differences of prices are defined as BTCF returns, r t = ln(p t /p t− 1 ) · 100. Following Schwert (1990) and Cao and Tsay (1993) , the volatility of the BTCF returns (ϑ t ) is estimated by the absolute value of returns minus its mean (μ r ), pre-multiplied by As in the previous literature, the trading volume and open interest series are converted into logarithms, i.e., log-volume (LV t ) and log-open interest (LO t ), to stabilize their variability and to reduce the non-normality of their distribution. Table1 sets out some summary statistics and preliminary test for the BTCF returns, volatility, trading volumes, and open interest. In particular, both return and volatility series exhibit highly dynamic movements and include the extreme values that may be created from the COVID-19 pandemic in terms of their major descriptive statistics such as minimum and maximum values, means and standard deviations. Sample distributions of return and volatility series are quite skewed and highly leptokurtic, and for all the series the Jarque-Bera test statistics firmly reject the null hypothesis of normality. The Ljung-Box Q test statistics for five and ten lags are also significant suggesting the presence of significant serial correlation in all the series, implying that they tend to be clustered together over time. In particular, the high correlation in the volume and open interest series indicates that the rate of information arrival is serially correlated and supports the validity of the MGARCH-type model proposed in this paper. The first panel of Fig. 1 shows the dynamic movements of the returns and the second panel exhibits turbulent volatility during the pandemic period. This turbulence was likely caused by market stress that probably formed under high uncertainty in the BTCF market, leading to the more dispersion of traders' belief and their overreaction to news, when the pandemic crisis introduced a shock to the market. The most conspicuous feature in Fig. 1 is that all series including the trading activity series do show up clear distinction in behavior between the pre-pandemic and the pandemic periods. This distinction implies primarily that we should make inference of a regime change in the BTCF market due to the COVID-19 pandemic. 29.18*** (0.000) 18.74*** (0.002) 516.84*** (0.000) 1174.9*** (0.000) Q(10) 32.48*** (0.000) 28.16*** (0.001) 908.62*** (0.000) 2258.3*** (0.000) Notes: Skewness and Kurtosis are coefficients of skewness and kurtosis, respectively. JB is the Jarque-Bera normality test, Q(M) is the Ljung-Box Q statistics at lag M. Significance levels (p-values) are in parentheses, and asterisks *** denote statistical significance at the 1% level. After Sims' (1980) pioneering work, the VAR model has become the most popular among multivariate time-series models for correlated series, e.g., returns, volatility, and trading behavior series in the BTCF market. However, the VAR model is useful only for modeling the first-order moment of the series and presupposes the implausible assumption of a constant variance. Thus, the VAR model is extended by employing a class of stochastic process called MGARCH for the inherent conditional heteroscedasticity and correlation structures of its residuals, which may change over time (Ku, 2008; Carnero and Eratalay, 2013) . Hence, the dynamics of the relationship between future returns, open interest, and trading volume can be modeled using the following SC-VAR-MGARCH model. This model combines the VAR model and MGARCH-DCC model (Engle and Sheppard, 2001; Engle, 2002) with structural changes. The concept of the MGARCH-DCC models implies that the conditional covariance matrix, H t , can be decomposed into conditional standard deviations, D t , and correlation matrix, R t , which are time-varying. Conditionally on regime s t = k and information set ϕ t− 1 , we specify the SC-VAR-MGARCH model as: is a n × 1 vector at time t, n is the number of series, s t is a state variable with a number of states K, C k is a n × 1 vector of constants, ε k t is a n × 1 vector of error terms, and B k p is a n × n matrix of autoregressive coefficients as follows: In regime k, the error terms ε k t is defined as: is a n × n positive definite matrix at time t and may be obtained by a Cholesky factorization of H k t . Note that H k t can be decomposed as follows: Where D k t is a n × n diagonal matrix of conditional standard deviation and R k t is a positive definite correlation matrix that allows to be time-varying with dynamic movement (Engle, 2002) The elements in the diagonal matrix D k t are conditional standard deviations: The conditional variance h k i,t is specified as the following GARCH(1,1) model in a vector form: Where , and d k is a n×1 vector and α k , β k are n × n matrices of coefficients. In case of DCC, the conditional correlation matrix R k t is time varying and is denoted as: Where R k t is decomposed into: Where a k and b k are scalars including 0 and must be satisfied of a k + b k < 1 to bound stationarity and Q k is the unconditional This paper presents SC-VAR-MGARCH models to investigate structural changes and trading behavior in the BTCF market. Under general conditions (Engle and Sheppard, 2001) , we conduct estimation methods to allow for two-step estimation. In the first step, a structural change caused by the COVID-19 pandemic is detected. In the second step, the parameters of the SC-VAR-MGARCH model are estimated using the structural break point estimate from the first step. This paper uses the Ross's (2015) test method to detect a structural change in the return series caused by the COVID-19 pandemic. The test method accounts for unknown structural changes and thus has important advantages to our analysis. Moreover, the method can effectively detect structural changes without prior distribution parameter information in various segments and structural changes in not only the first-order moment but also in the second-order moment of non-Gaussian sequences. The test method generally assumes a fixed sequence of T observations x 1 , … , x T from the random variables X 1 , … , X T . These random variables may contain at least one state changes in the probability distribution at the unknown change points δ 1 , δ 2 , . . ., δ K− 1 . Then, K is a finite number of states. If structural change points do not exist, then the observation sequence is independent and identically distributed according to a distribution, denoted by F 0 . In contrast, if structural changes are present, then the observation sequence after the structural change point is distributed as F 1 , where F 1 ∕ = F 0 . Hence, the distribution of the sequence can be written as: Where the F i stands for the distribution in each segment. The following hypotheses can be established to test whether a structural change occurs immediately after a particular observation ω: Where ψ i is the potentially unknown parameters of each distribution. This set-up becomes a standard problem that can be solved by a two-sample hypothesis test with test statistic S ω,T . In other words, when S ω,T is greater than the value of an appropriately chosen threshold π ω,T , 7 we reject the null hypothesis that the two samples have identical distributions. Thus, we can conclude that a structural change point has occurred immediately after observation x ω . Therefore, for all points ω (1 < ω < T), the test statistic S ω,T is evaluated, and a point ω with a maximum value is considered as the structural break point. The test statistic is defined as: Where S ω,T are obtained from standardization of S ω,T that have the means μS ω,T and the standard deviations σS ω,T . Since it is assumed that the change point is unknown in advance, Ross (2015) suggests immediately following the value of ω that maximizes S T as the best estimate of the change point location: Given the estimates of break points in the first step, the SC-VAR-MGARCH model can be estimated using maximum likelihood. In practice, since there are too many parameters in the model, for compactness let's define the following vector of the parameters in the model: ( Where 7 The threshold is chosen to limit the Type 1 error rate α, like the standard statistical hypothesis test. That is, α is the probability of incorrectly deciding that a change has occurred if no change has actually occurred. In this context, π ω,T must be chosen as the upper α quantum of the S ω,T distribution in the null hypothesis. where Λ k is the vector of parameters in the mean equation (Eq. 1), Π k is the vector of parameters in the conditional variance equation (Eq. 7), P k is the vector of parameters in the conditional correlation equation (Eq. 9), the vec operator is the vectorization of the indicated matrix, and the vech operator is the vectorization of the lower triangular part of the indicated matrix. Under the assumption of a normal distribution of error terms in the models, the conditional MLE (maximum likelihood estimator) of Ψ k on s t = k should be derived by maximizing the multivariate Gaussian log-likelihood function: Where T k is the number of observations in regime k (t = 1, 2, … , T k ). From Eq. 4, the equation can be rewritten as Where w k t = D k t − 1 ε k t are the standardized errors. The COVID-19-induced market stress has created unprecedented volatility that may be closely linked to trading behavior variables, such as open interest as a proxy for belief dispersion. Traders with relatively different beliefs are dominant in the regime of market stress because their biases increase with the belief updating to respond to news (Basak and Atmaz, 2018) . Therefore, in the context of the COVID-19-induced market stress, we investigate the dynamic relationships between returns, volatility, and trading behavior variables in the BTCF market using the SC-VAR-MGARCH model and Granger causality tests. This subsection focuses on whether the COVID-19 pandemic has initiated a structural change by applying the test method introduced in subsection 3.2.1. to detect a structural change in the returns series of the BTCF. As mentioned above, the test method does not require a known break point and prior distribution parameter information in various regimes. We implement the test and estimate the change point location. The test statistic S T is 17.2084, which is greater than the value of an appropriately chosen threshold. Thus, we reject the null hypothesis of no structural change. The value of ω that maximizes S T is 129, which is the best estimate of the change point location. In other words, the date of the estimated change point is March 5, 2020. This date is a few days ahead of the date when the COVID-19 pandemic was declared, which may be a potential date of structural break. Table 2 shows summary statistics for each regime. While Regime I can be regarded as a normal situation with volatility up to 11.0240, Regime II is a pandemic situation, which is a very stressful with volatility up to 62.8833. It is worthwhile to note that in We estimate the SC-VAR(1)-MGARCH(1,1) model given the detected break point as follows 9 : ' 10 , and DLO k t = detrended log-open interest, DLV k t = detrended log-volume 11 , the conditional variance of the error terms ε k t is defined as a GJR-GARCH model to consider asymmetric volatility: Where I t− 1 is an indicator function, i.e., I t− 1 = 1 if ε k i,t− 1 < 0, I t− 1 = 0 otherwise. Orders of the model are determined to minimize the value of SIC criterion and to avoid over-parameterization leading to converging problems. Multivariate ARCH tests for the residuals obtained from VAR models are performed first to determine the validity of this model. Table 3 shows that the test statistics for 5 and 10 lags are significantly large to indicate the presence of ARCH effects. These results explicitly indicate the need for further use of the MGARCH-type model to expand the understanding of the series relations. In addition, the strong ARCH effects in all specifications of each regime should be closely related to the time series dependence of the trading behavior variables shown in the previous data analysis. The empirical results for the SC-VAR(1)-MGARCH(1,1) models of in Tables 4 and 5 , respectively. This study mainly examines how the time-varying relationship between returns (or volatility), trading volume, and open interest reacts to the COVID-19 pandemic. Therefore, the results of regime I (pre-pandemic) and regime II (pandemic) must be compared. Table 4 indicates that returns and trading behavior in regime I have no significant relationship, Multivariate ARCH (5) Notes: DLO t = detrended log-open interest, DLV t = detrended log-volume, ARCH(M) is the LM statistics for ARCH test at lag M. Significance levels (pvalues) are in parentheses. asterisks *, **, and *** denote statistical significance at the 10 %, 5%, and 1% levels, respectively. 8 Clear identification of financial market stress is challenging due to the difficulty in measuring sentiment in markets. In that perspective, the use of OFR Financial Stress Index, which measures disruptions in the normal functioning of financial markets as a systemic financial stress, can be an alternative solution to the problem of identifying market stress. The OFR FSI is constructed from 33 financial market variables, such as yield spreads, valuation measures, and interest rates. 9 According to existing studies (e.g., Najand and Yung, 1991) , when contemporary relationships are considered in the GARCH-type models, the estimates might suffer from simultaneity bias because the variables are endogenous as in theoretical models based on information flow. Hence, this study does not reflect contemporary relationships. 10 Some empirical studies (e.g., Bessembinder and Seguin, 1993; Galati, 2000) distinguish the variables of volume and open interest into the expected and unexpected parts. As far as the author knows, however, there is no existing theoretical analysis to distinguish the variables into the expected and the unexpected parts. Furthermore, this study focuses on how behavioral biases due to COVID-19 pandemic affect the relationship between volatility and the variables that include all information on trading behavior such as serial correlation. Therefore, in order not to obscure the point, it is not necessary to divide the variables into the expected and the unexpected parts in this paper. 11 As in existing literature (Galati, 2000; Park, 2010; Czudaj, 2019, etc.) , it is necessary to filter out trends that are present in the time series of the trading volume and open interest. whereas one lagged trading volume negatively affects returns and one lagged return negatively affects open interest in regime II. As mentioned above, the relationships are theoretically controversial. When traders have rational expectations and use prices correctly, open interest acting as a proxy for the dispersion of traders' beliefs and trading volume are positively related to expected returns and volatility. By contrast, when traders do not use the price or have behavioral biases, as in a difference of opinions model (Bessembinder et al., 1996) , these relationships are reversed (Diether et al., 2002; Yu, 2011) . Therefore, the latter interprets that these relationships are negatively significant in regime II (pandemic) due to behavioral biases, such as overreaction to COVID-19-related news, in contrast to regime I (pre-pandemic) in which the relationships are not significant. The peculiarity that the relationship between volatility and trading behavior differs in the two regimes is also confirmed in Table 5 . Trading behavior variables have no significant relation to volatility in Regime I. However, one lagged open interest and trading volume positively influence volatility, and one lagged volatility positively influences open interest in Regime II. In line with the difference of opinions model, belief dispersion is a crucial factor and may reflect information asymmetry or information uncertainty. Therefore, an increase in belief dispersion reduces BTCF prices but contributes to higher volatility. Despite conflicting claims about the relationship between returns (or volatility) and trading behavior, the increase of belief dispersion is known to mostly raise trading volume (Goetzmann and Massa, 2005) , and the estimates of the relevant coefficients are also consistent with this. According to estimates of GARCH parameters (α, β, γ), all series for the sample period in Regime II tend to be more clustered over time and to have the greater persistence of shocks as compared with the sample period in Regime I. Overall, it can be assured that pandemic-induced market stress plays a significant role in trading behavior variables and their correlations to returns and volatility on the BTCF market. As a next step, to clarify the causal relationship between the variables, Granger causality tests are implemented and reported in Table 6 . In Regime I, returns do not Granger-cause trading volume/open interest nor does it cause a reverse relationship at conventional levels of significance, whereas in Regime II, returns Granger-causes trading volume/open interest and trading volume Grangercauses returns at conventional levels of significance. Using volatility instead of returns, their causal relationship becomes more significant. Further, a bi-directional causality relationship between volatility and trading volume/open interest in Regime II presumably support our intuition based on the difference of opinions model that severe market stress caused by the COVID-19 pandemic might generate the trader's behavioral biases, leading to a drastic change in the causal relations of the variables. The dynamic relations between returns (r t ) and trading variables such as DLO t and DLV t in the BTCF market: Parameter estimates of the SC-VAR- Notes: DLO t = detrended log-open interest, DLV t = detrended log-volume, Av. lnL represents the mean value of the maximized log-likelihood function, asterisks *, **, and *** denote statistical significance at the 10 %, 5%, and 1% levels, respectively. Turning to dynamic conditional correlations (DCC) between returns/volatility and trading behavior variables, we find that since the sum of estimates of DCC parameters a and b for each specification is less than 1 in Tables 4 and 5, the correlation itself is mean reverting and it fluctuates around the unconditional correlation. That is, a + b = 0.5843 for Regime I and 0.7403 for Regime II in the . On the other hand, a + b = 0.0592 for Regime I and 0.0298 for Regime II in the model of . Since a and b are the parameters of the news term and the decay term respectively, the closer the sum is to 1, especially the greater a is, the more severely the conditional correlation changes over time. Thus, the conditional correlations for the sample period in Regime II are stronger and their movements are more stable than Regime I, which are also visually evident by Figs. 2 and 3. From the figures, we can get insights on the dynamics between returns/volatility, trading volume, and open interest on the BTCF Notes: DLO t = detrended log-open interest, DLV t = detrended log-volume, Av. lnL represents the mean value of the maximized log-likelihood function, asterisks *, **, and *** denote statistical significance at the 10 %, 5%, and 1% levels, respectively. Notes: Chi-sq test statistics and their p values of the test are reported with the null hypothesis of unconditional Granger, and asterisks *, **, and *** denote statistical significance at the 10 %, 5%, and 1% levels, respectively. market. In the presence of belief dispersion, the BTCF price and trading volume are more sensitive to news but the variability of their conditional correlations become lower in comparison to relatively normal states. The empirical analyses highlight that it is important to account for regime change initiated by pandemic-induced market stress when modeling the relationship between returns/volatility, trading volume and open interest for the BTCF market. This paper explains the dynamic relationship between returns (or volatility) and trading behavior variables, such as trading volume and open interest, as a proxy for dispersion of traders' beliefs on asset prices. Moreover, this study investigates how the COVID-19 pandemic impacts their interrelations. Accordingly, this paper considers the BTCF market for empirical analysis given its certain advantages and proposes the SC-VAR-MGARCH model. This study contributes to literature in several ways. First, this study is significant in that it is the first to explore the impact of pandemic-induced market stress on the dynamic relationship between price and trading behavior variables in the futures market. Furthermore, this study also explains that the ongoing COVID-19 pandemic raises concerns about future market stress with greater economic risk. Second, this study finds that the COVID-19 pandemic has initiated a significant regime change. The estimated change point is on March 5, 2020, which is a few days ahead when the COVID-19 pandemic was declared. The declaration date may be a potential date of structural break. Third, this study can help investors, hedgers, and arbitrageurs expand their understanding of returns volatility sources and the dynamics of the relationships under extreme financial market stress beyond the BTCF market, thus establishing excellent trading and portfolio management strategies. The estimates of the SC-VAR-MGARCH models are remarkably consistent with the difference of opinions model during the pandemic-induced market stress, though existing evidence of the relationships are ambiguous. As a result, an increase in belief dispersion due to market stress reduces BTCF returns but significantly raises trading volume and volatility in the pandemic regime in contrast with the pre-pandemic regime. Therefore, our empirical results sufficiently support the SC-VAR-MGARCH model considering the pandemic-induced market stress and emphasize the intrinsic role of market stress causing traders' behavioral biases, such as overreaction to news. Our study suggests future research to consider the dynamic relationship between price and trading variables as a more reliable model for reflecting traders' herd behavior generated from market stress and uncertainty. In addition, when high-frequency time-series data of the COVID-19 pandemic period sufficiently accumulate in the future, the COVID-19 impacts on realized volatility, realized jumps, and trading behavior in asset markets, in addition to the BTCF market, must be closely analyzed. A theory of intraday patterns: volume and price variability The relationship between trading volume and jump processes in financial markets Regime changes in Bitcoin GARCH volatility dynamics Stock markets' reaction to COVID-19: cases or fatalities? Dispersion in analysts' earnings forecasts and credit rating Belief dispersion in the stock market Price volatility, trading volume, and market depth: evidence from futures markets An empirical examination of information, differences of opinion, and trading activity Nonlinear time-series analysis of stock volatilities Estimating VAR-MGARCH models in multiple steps A subordinated stochastic process model with finite variance for speculative process A model of asset trading under the assumption of sequential information arrival Dynamics between trading volume, volatility and open interest in agricultural futures markets: a Bayesian time-varying coefficient approach Differences of opinion and the cross section of stock returns Divergence of opinion and equity returns Dynamic conditional correlation: a simple class of multivariate Garch models Theoretical and Empirical Properties of Dynamic Conditional Correlation Multivariate Garch The stochastic dependence of security price changes and transaction volumes: implications for the mixture-of-distribution hypothesis Price discovery in bitcoin futures Stochastic volatility, trading volume, and the daily flow of information Trading volumes, volatility and spreads in foreign exchange markets: evidence from emerging market countries. BIS Working Papers 93 How do equity markets react to COVID-19? Evidence from emerging and developed countries Index futures volatility and trading activity: measuring causality at a multiple horizon An equilibrium model of asset trading with sequential information arrival Do investors herd in cryptocurrencies -and why? Trading volume and price reactions to public announcements Market risk and Bitcoin returns Student-t distribution based VAR-MGARCH: an application of the DCC model on international portfolio risk management A generalized bivariate mixture model for stock price volatility and trading volume On measuring speculative and hedging activities in futures markets from volume and open interest data A GARCH examination of the relationship between volume and price variability in futures markets Surprising information, the MDH, and the relationship between volatility and trading volume Parametric and nonparametric sequential change detection in R: the cpm package Stock volatility and the crash of' 87' COVID-19's disasters are perilous than Global Financial Crisis: a rumor of fact? Macroeconomics and reality The price variability-volume relationship on speculative markets Bayesian change point analysis of Bitcoin returns Disagreement and return predictability of stock portfolios