Conditional Dependence between Oil Prices and Exchange Rates in BRICS Countries: An Application of the Copula-GARCH Model

: We studied the dependence structure between West Texas Intermediate (WTI) oil prices and the exchange rates of BRICS 1 countries, using copula models. We used the Normal, Plackett, rotated-Gumbel, and Student’s t copulas to measure the constant dependence, and we captured the dynamic dependence using the Generalized Autoregressive Score with the Student’s t copula. We found that negative dependence and signiﬁcant tail dependence exist in all pairs considered. The Russian Ruble (RUB)–WTI pair has the strongest dependence. Moreover, we treated ﬁve exchange rate–oil pairs as portfolios and evaluated the Value at Risk and Expected Shortfall from the time-varying copula models. We found that both reach low values when the oil price falls sharply.


Introduction
Crude oil, as a non-renewable resource and a crucial commodity, plays a significant role in the world's economy. Over the past few decades, the crude oil price has experienced great fluctuations compared to the period from the Second World War to the early 1970s. For example, in June 2014, the West Texas Intermediate (WTI) crude oil price reached a peak of 105 US dollars per barrel, and then fell sharply to 59 US dollars per barrel in December 2014. In the early 1980s, many studies pointed out that oil price dynamics influenced economic activity (Hamilton 1983;Gisser and Goodwin 1986;Mork 1989). On the other hand, currency markets have also suffered multiple crises in recent years, such as the Latin American currency crisis of 1994, the Asian financial crisis of 1997, and most recently the fluctuation in the exchange rate between the US dollar and Chinese Yuan after a trade war, as each country continued to dispute tariffs in 2018. With these notable events, every abnormal foreign exchange movement had a huge impact on the economy. Analyzing the relationship between crude oil prices and exchange rates yields considerable information about the increasing importance of oil and currency markets in economics for market operators, investors, and economists.
On the theoretical side, the extant literature (Bénassy-Quéré et al. 2007;Bodenstein et al. 2011;Habib et al. 2016;Beckmann et al. 2017) has pointed out that there are three channels between oil price shocks and exchange rates: the terms of trade channel, wealth effects channel, and portfolio reallocation channel. The terms of trade channel considers a simple two-country open economy static model, which has two goods sectors-a traded goods sector and a non-traded goods sector. When there is a positive term of trade shock, the price of a basket of traded goods or non-traded goods in the domestic economy is relatively higher than it is in a foreign economy, leading to the appreciation of 1 BRICS is the acronym of an association of five major emerging economies: Brazil, Russia, India, China and South Africa.
J. Risk Financial Manag. 2019, 12, 99 2 of 25 domestic currency. Hence, when oil prices rise, the currency of an energy-dependent country will appreciate, and vice versa. The wealth effects channel reflects the short-term effects while the portfolio reallocation channel reflects the medium-and long-term effects. When there is a positive oil price shock, the wealth transfers from the oil-importing country to the oil-exporting country, leading to large shifts in current account balances and portfolio reallocation. For this purpose, in order to adjust to clear the trade balance and asset markets, the currency of the oil-exporting country appreciates, and the currency of oil-importing country depreciates. Chen and Chen (2007) provide a simple model to illustrate the above theory: where p (p * ) is the log-linear approximation of the domestic (foreign) consumer price index, and p T (p T* ) and p N (p N* ) are the prices of the traded and non-traded sectors in the home country (foreign country), respectively. The α (α * ) weights correspond to expenditure shares on traded goods near the point of approximation for the domestic (foreign) country. The log of the real exchange rate q is defined as where s is the log of the nominal exchange rate. Therefore, the real exchange rate could be written as If α α * , when the oil price rises, the price of the traded sector in the country which is more dependent on oil importation will increase, and thereby cause a depreciation. 2 On the empirical side, there are also many studies that investigate the relationship between oil prices and exchange rates. Some of this literature finds a negative relationship. Akram (2009) uses quarterly data from 1990Q1 to 2007Q4, and a structural value at risk (VaR) model to analyze whether a decline in the US dollar contributes to higher commodity prices. Their results suggest that a weaker dollar leads to higher commodity prices, including oil prices. Basher et al. (2012) study the dynamic relationship between oil prices, exchange rates, and emerging market stock prices by modeling daily data from 1988 to 2008 using a structural VaR model. They find that positive shocks to oil prices lead to an immediate drop in the US dollar exchange rate, which has a statistically significant impact in the short-run. Reboredo and Rivera-Castro (2013) use a wavelet decomposition approach to investigate the relationship between oil prices (WTI) and US dollar exchange rates for a large set of currencies, including developed and emerging economies, net oil-exporting and oil-importing economies, and inflation-targeting countries. They find that oil price changes had no effect on exchange rates or vice versa in the pre-crisis period, but after the 2008 global financial crisis, they find a negative interdependence. By contrast, other research finds a positive relationship between oil prices and exchange rates. Bénassy-Quéré et al. (2007) study co-integration and causality between the real price of oil and the real value of the dollar from 1974 to 2005 using monthly data. Their results suggest that a 10% rise in the oil price coincides with a 4.3% appreciation of the dollar in the long run, and the causality runs from oil to the dollar. Ding and Vo (2012) use the multivariate stochastic volatility and the multivariate generalized autoregressive conditional heteroscedasticity (GARCH) models to investigate the volatility interactions between the oil market and the foreign exchange market. They divided the data into two parts to capture the structural breaks in the economic crisis, and they show that before the 2008 crisis, the two markets responded to shocks simultaneously and no interaction is found. However, during the financial crisis, they find a positive bidirectional volatility interaction between the two variables.
In addition to research in developed markets, other studies focus on this relationship in emerging markets. For example, Narayan et al. (2008) estimate the impact of oil prices on the nominal exchange rate on Fiji Island with GARCH and EGARCH 3 models. They find that a rise in oil prices leads to an appreciation of the Fijian dollar. Narayan (2013) uses exchange rate data from 14 Asian countries and demonstrates that a higher oil price leads to a future depreciation of the Vietnamese Dong, but a future appreciation of the local currencies of Bangladesh, Cambodia, and Hong Kong.
Analysis of the correlation between the currency and oil markets also provides abundant information about whether the oil price could act as a hedge or as a safe haven against exchange rate risk. Much literature investigates the relationship between different assets to examine whether these variables act as hedges against each other. The methods used in these papers are varied, such as linear regression, quantile regression, etc. Considering that the characteristics of skew and leptokurtic kurtosis often have asymmetric behaviors, and the oil price and exchange rates seem unsuitable for modeling with linear models, in this study, we propose a copula model to investigate the dependency between oil prices and exchange rates. Compared to traditional methods, a copula model has several advantages. First, it can provide a valid joint distribution in combination with any univariate distribution, which makes it more flexible in modeling multivariate distributions. Second, it can capture a wide range of dependence structures, such as asymmetric, nonlinear, and tail dependence in extreme events. Prior studies apply copula models to examine the relationship between oil prices and exchanges rates (Patton 2006;Wu et al. 2012;Aloui et al. 2013;Ji et al. 2019). However, most of them focus on developed markets.
In this study, we investigate the dependence structure between the oil and currency markets. Our contribution to the literature is twofold. First, we propose both constant and time-varying copula-GARCH-based models to measure the dependence between the two markets elastically. By using a copula model, we can study the co-movement between oil prices and exchange rates during bearish and bullish markets, and we can also examine whether oil acts as a hedge against exchange rate risk in BRICS countries. Second, we focus on exchange rate data from the BRICS countries, which is a group of five developing countries (Brazil, Russia, India, China, and South Africa), to determine whether the dependence structure in emerging markets is the same as that in developed countries.
The rest of this paper proceeds as follows. In Section 2, we introduce the empirical methodology and estimation strategy. Section 3 presents the data and discusses the empirical results. In the final section, we present our conclusions. The results of ARMA-GARCH model and the robustness analysis are presented in Appendices A and B, respectively. Sklar (1959) introduced the notion of the copula, in which an n-dimensional distribution function which can be decomposed into two parts, namely, the marginal distribution and the copula.

Empirical Methodology
We denote X ≡ [X 1 , . . . , X n ]' as the variable set of interest, and X i follows the marginal distribution F i for each i∈[0, n].
Following Sklar's theorem, there exists a copula C : [0, 1] n → [0, 1] with uniform marginals that can map the univariate marginal distribution F i to the multivariate distribution function F.
The density of the multivariate distribution function F is 3 EGARCH is an abbreviation for Exponential Generalized Autoregressive Conditional Heteroscedastic. EGARCH model was introduced by Nelson (1991).
We first obtain standardized residuals using an autoregressive moving-average GARCH (ARMA-GARCH) model for the daily returns of crude oil prices and exchange rates. Then, we transform the standardized residuals into a uniform distribution. Finally, we input the transformed uniform variates into the copula function to model the dependence structure.

ARMA-GARCH Model
To model the daily returns of crude oil prices and exchange rates, we use autoregressivemoving-average (ARMA) models with two kinds of GARCH models--the standard GARCH (sGARCH) model and the Glosten-Jagannathan-Runkle GARCH (GJR-GARCH) model. In both GARCH models, we assume that the standard errors from the GARCH models follow a skewed t distribution.
We denote r t as the daily returns, which follow where φ 0 is a constant and ε t is a white noise error item. Table 1 presents the functions for the two GARCH models.
Note: this table presents the functions of the GARCH model and the Glosten-Jagannathan-Runkle GARCH (GJR-GARCH) model. ε t is a zero-mean white noise item, which is assumed to follow a skewed t distribution, µ is the expected return, and z t is an i.i.d (independent and identically distributed). random variable with zero mean and a variance of one.

Marginal Distribution
Following Patton (2013), we consider both parametric and non-parametric models in modeling the marginal distribution.
For the parametric estimation, we choose the skewed t distribution as in Hansen (1994). The density function is The skewed t distribution has two parameters to control its shape: a skewness parameter λ∈(−1, 1), which determines whether the mode of the density is skewed to the left or to the right, and a degrees of freedom parameter ν∈(2, ∞), which controls the kurtosis of the distribution. When λ = 0, we obtain the Student's t distribution; when ν→∞, we obtain the skewed normal distribution; and when λ = 0 with ν→∞, we obtain the standard normal distribution. For this property, the skewed t distribution has more flexibility in fitting the data than do the normal distribution or the Student's t distribution.
For the non-parametric estimation, we use the empirical distribution function (EDF): where T is the length of the data andε it indicates the standardized residuals from the GARCH model.

Constant Copula Model
There are three families of copula model: elliptical copulas (Normal and Student's t); Archimedean copulas (Gumbel, Clayton, and Frank); and quadratic copulas (Plackett). In this study, we chose the Normal, Plackett, Rotated Gumbel, and Student's t copulas to investigate the dependence structure between oil prices and exchange rates. 4 The normal and Student's t copulas are based on an elliptical distribution, such as the normal or Student's t distribution. The normal copula is symmetric and has no tail dependence. θ > 0 and θ < 0 lead to positive and negative dependence, respectively. The Student's t copula is also symmetric, but can capture tail dependence in extreme events. When θ 1 = 0 with θ 2 →∞, then it becomes independent between the two variables. Like the Normal copula, the Plackett copula is also symmetric and cannot capture either lower or upper tail dependence. 0 < θ < 1 implies negative dependence and θ > 1 implies positive dependence. The Gumbel copula can only capture extreme lower dependence, θ = 1 and θ→∞ indicate independence and perfectly negative dependence, respectively. Table 2 summarizes the properties of the abovementioned bivariate copula models.

Time-Varying Copula Model
In this study, we examined the time-varying Student's t copula using a Generalized Autoregressive Score (GAS) model. Creal et al. (2013) proposed the GAS model, which provides a general framework for modeling time variation. We consider δ t as the parameter vector of the time-varying copula models, which we can express as the following function: where f t is the time-varying parameter vector in the GAS model. Following Creal et al. (2013), we can assume that f t is given by the familiar autoregressive updating equation: where ω is a constant, S t−1 is a time-dependent scaling matrix, c (·) is the density function of the copula model, and u t is the vector of the probability integral transforms using the univariate marginals. According to Creal et al. (2013), we can set the scaling matrix S t−1 to be equal to the pseudo-inverse information matrix: For the time-varying Student's t copula, we use the function to ensure that the parameter lies within (−1, 1).
Note: this table presents the basic properties of the Normal, Plackett, rotated-Gumbel, and Student's t copula models. u 1 and u 2 are uniform variates. In the rotated-Gumbel copula, denotes the inverse cumulative distribution function of the univariate standard normal distribution, and t −1 (·) indicates the inverse cumulative distribution function of the univariate Student's t distribution. F t (·) denotes the probability density function of the standard Student's t distribution. This table is based on Aloui et al. (2013), Patton (2013) and Nelsen (2006).

Data and Summary Statistics
Our dataset contained the WTI oil prices in dollars and the nominal dollar-denominated exchange rates for the Brazilian Real (BRL), Russian Ruble (RUB), Indian Rupee (INR), Chinese Yuan (CNY), and South African Rand (ZAR). After 2000, China reformed its exchange rate regime twice. The first time was in 2005, when the authorities adopted a managed floating exchange rate system with reference to a basket of currencies, and the second time was in 2010, when the CNY shifted from a bilateral to a multilateral reference, with greater emphasis on the currency basket. Before each reform, CNY fluctuated with no change or within a narrow range. We thus chose the offshore exchange rate for CNY for our analysis. Our sample period was from 4 October 2010 to 11 December 2018, to match the data availability for the offshore CNY exchange rate. We downloaded all data from the DataStream database.
We obtained the stationary return series using the following function: Since many events occurred during our sample period, such as China's central bank devaluation of its currency in 2015, the United Kingdom's vote to leave the European Union in 2016, etc., which may be directly or indirectly relevant to oil prices and exchange rates, we applied a methodology from Perron (1998, 2003) for testing the structural breaks, to avert unexpected changes in our dataset. The p-values of structural break tests on the return series are summarized in Table 3. The results showed an acceptance of the null hypothesis that there are no breaks, therefore we were able to use this sample period for further analysis.  Perron (1998, 2003). Figure 1 shows the plot of the time-path of the return series, and Table 4 reports the descriptive statistics for daily returns. The mean of CNY is close to zero, and its standard deviation is the lowest of all return series. Thus, in some way, CNY remained stable under the government regulations. RUB has the most fluctuations among the five exchange rates, especially at the end of 2014, when the price of crude oil declined by nearly 50%, which was a major cause of the financial crisis in Russia. All return series are skewed and have high kurtosis, indicating that the distributions of the return series show obvious non-normality characteristics. The Jarque-Bera test performed on the return series also verifies our analysis.
Before beginning the copula analysis, we provide the results of the rank correlation and linear correlation in Table 5 for comparison. All coefficients are significant, indicating negative relationships between all currency-oil pairs. The RUB-WTI pair has the highest coefficient for both the rank and linear correlations, possibly because Russia is the second largest oil-exporting country in the world, and Russia exports more oil than the other countries in our sample.
Note: this table represents the Kendall's τ correlation coefficient, which measures the strength of dependence between two variables; the Spearman's ρ correlation coefficient, which measures the degree of association between two variables; and the Pearson correlation coefficient, which measures the degree of the relationship between linearly related variables. *** indicates significance at 1% level.

Constant Copula Result
Tables 6-10 report the estimations of the four constant copulas (Normal, Plackett, rotated-Gumbel, and Student's t) for five currency-oil price pairs (BRL-WTI, RUB-WTI, INR-WTI, CNY-WTI, and ZAR-WTI), respectively. The parameters from the copula model are estimated based on the maximum likelihood method. In all pairs, the Student's t copula has the highest value of log-likelihood, while the rotated-Gumbel copula has the lowest. Almost all estimated copula parameters are significant at the 1% level. The results of the copula models reveal a significant dependence between the oil price and exchange rates. All copula models suggest a negative dependence among all considered pairs, for both the parametric and semi-parametric cases, which could be because oil is considered as a hedge against inflation in BRICS countries. For oil-importing countries, like China, India, and South Africa, the negative dependence is in line with the theoretical method, which illustrates that the currency of an oil-importing country will depreciate when oil price goes up. This result is also in accordance with other studies, such as that of Ghosh (2011), who assesses the relationship between the Indian Rupee and oil price and finds the existence of a negative correlation. Kin and Courage (2014) find that an increase in oil prices leads to a depreciation of the South African rank exchange rate. However, for oil-exporting countries such as Russia and Brazil, the result is the opposite with the theoretical model. We still found some studies which have the same results as those found in our analysis, such as Blokhina et al. (2016), who find a negative relationship by applying a linear regression to Russia. In each copula model, the RUB-WTI pair has the highest dependence coefficient, followed by ZAR-WTI, BRL-WTI, CNY-WTI, and INR-WTI, which is consistent with the rank and linear correlation results.
According to Genest et al. (2009), the Cramer-von Mises (CvM) test is the most powerful test to check the goodness-of-fit of copula models. Here, we used both the CvM test and the Kolmogorov-Smirnov (KS) test to check the quality of the data fit. We defined the CvM and KS tests for the bivariate case as:Ĉ Table 11 presents the p-values of the goodness-of-fit test. For all pairs, the rotated-Gumbel copula is rejected at the 1% level of significance, except for in the parametric case of INR-WTI. The Student's t copula generates the highest p-value in most cases, indicating that it provides the best fit for the data.
Tail dependence describes co-movements between two variables during extreme events. In our case, lower tail dependence measures the probability of currency depreciation when the oil price falls. In contrast, upper tail dependence measures the probability of currency appreciation when the oil price rises. In Table 12, we present the tail dependence coefficients from the constant Student's t copula model, which can capture symmetrical tail dependence. The results show significant extreme market co-movements between the exchange rate and oil price. INR-WTI has the highest tail dependence, while RUB-WTI has the lowest for both the parametric and semi-parametric cases.        Table 13 presents the estimated time-varying Student's t copula model using the GAS model. For all pairs, theβ estimate is highly significant and has a large value, except for the INR-WTI pair, indicating that the time-varying copula model parameters in this period have a great impact on the copula parameters of the following period. Figure 2 presents the plot of the dynamic estimated parameters from time-varying Student's t copula in the semi-parametric case. The red line indicates the average of the estimated time-varying parameters. As in the results from the constant copula models, all pairs reveal negative dependence structures in the following order: RUB-WTI, ZAR-WTI, BRL-WTI, CNY-WTI, and INR-WTI. This result is in line with our results from constant copula models. Figure 3 reports the plot of the conditional tail dependence with time deviation in the semi-parametric case. For clarity, we used different scales for the y axis for the different pairs. In general, for all pairs, the tail dependence from the time-varying copula, which is less than 0.0025, is much lower than that of the constant copula models (about 0.4). parameters. As in the results from the constant copula models, all pairs reveal negative dependence structures in the following order: RUB-WTI, ZAR-WTI, BRL-WTI, CNY-WTI, and INR-WTI. This result is in line with our results from constant copula models. Figure 3 reports the plot of the conditional tail dependence with time deviation in the semiparametric case. For clarity, we used different scales for the y axis for the different pairs. In general, for all pairs, the tail dependence from the time-varying copula, which is less than 0.0025, is much lower than that of the constant copula models (about 0.4).  Following Patton (2013), we tested the goodness-of-fit of the time-varying copula model by transforming the data through the Rosenblatt method, which is a multivariate version of the Following Patton (2013), we tested the goodness-of-fit of the time-varying copula model by transforming the data through the Rosenblatt method, which is a multivariate version of the Following Patton (2013), we tested the goodness-of-fit of the time-varying copula model by transforming the data through the Rosenblatt method, which is a multivariate version of the probability integral transformation proposed in Diebold et al. (1999) and developed further by Rémillard (2017). The transform uses the function: V 1t = U 1t ∀t,

Time-Varying Copula Result
Table 14 reports the p-values of the goodness-of-fit test. In the parametric case, the time-varying Student's t copula model is accepted with high p-values, except for INR-WTI pair with the KS test. In the semi-parametric case, the time-varying Student's t copula model is rejected for the BRL-WTI, INR-WTI, and ZAR-WTI pairs with the CvM test. Based on the results from the goodness-of-fit test, we consider that the constant copula model provides a better fit than the time-varying copula model for the INR-WTI pair.

Value at Risk and Excepted Shortfall
VaR and expected shortfall (ES) are statistics that quantify the risk of a portfolio. VaR measures the maximum expected loss with a confidence level, and ES measures the expected loss when the portfolio return is greater than the value of the VaR calculated with that confidence level. In this section, we treat the five exchange rate-oil price pairs (BRL-WTI, RUB-WTI, INR-WTI, CNY-WTI, and ZAR-WTI) as equally weighted portfolios, and then we calculate dynamic VaR and ES using a Monte Carlo simulation. The procedure has the following steps.
Step 1: at each time t (t is from 1 to T, where T is the length of the return series), generate 5000 random trials of uniform variables using the parameters estimated from the time-varying Student's t copula model. Step 2: transform the random uniform variables into standardized residuals, then conduct a series of new asset returns using the estimated parameters from the GARCH models.
Step 3: use the portfolio weights (1% in this study) and the new asset returns to calculate the portfolio returns and then calculate the VaR and ES.
Step 4: repeat steps 1-3 T times. Figure 4 plots the results. We can see that the VaR begins at 3%, decreases to 6% at the beginning of 2015 (9% for RUB-WTI), and from the middle of 2015 it decreases to reach approximately 8% by the beginning of 2016 for all pairs, due to the falling crude oil prices. Similarly, the ES also reaches a low value at that period.

Conclusions
This study examined both the constant and dynamic dependencies between crude oil prices (WTI) and exchange rates for the BRICS countries, which includes oil-exporting and oil-importing economies, from 4 October 2010 to 11 December 2018. We employed constant and time-varying copula-GARCH-based models to capture the nonlinear, asymmetric, and tail dependence between the two variables.
The empirical results showed, first, that a significant dependence existed in all exchange rateoil price pairs, and the dependence is negative, which means that oil serves as a hedge against rising inflation in BRICS countries. For example, when one currency appreciates, consumers from that country will find the oil price less expensive and then increase their demand for oil, which will lead to a rise in the oil price. Second, in both the constant and time-varying copula models, the RUB-WTI pair has the strongest dependence, possibly because Russia is a large oil producer and has a high share of global oil exports, making Russia itself more dependent on oil. Third, from the goodness-offit test results, the Student's t copula model has the best performance among all the constant copula models considered, while for the time-varying copula models, it was rejected for the BRL-WTI, INR-WTI, and ZAR-WTI pairs. We also measured risk management by treating the five exchange rate-oil price pairs as portfolios and calculating the dynamic values of VaR and ES. Both reached very low values for all pairs when the oil price experienced a sharp decrease.
The results of this paper offer us several suggestions. First, the significant correlation between oil price and exchange rate indicates that enterprises and governments in BRICS countries, when the oil price fluctuates, should pay attention to the exchange rate risk. Second, the negative relationship between oil price and exchange rate shows that the oil could be a hedge against inflation in BRICS countries, which is useful for foreign investors to manage exchange rate risk.
However, in this paper, we only considered bivariate copula models. For the future extension of this research, we will focus on a high-dimensional copula model, which could provide a more flexible analysis.

Conclusions
This study examined both the constant and dynamic dependencies between crude oil prices (WTI) and exchange rates for the BRICS countries, which includes oil-exporting and oil-importing economies, from 4 October 2010 to 11 December 2018. We employed constant and time-varying copula-GARCH-based models to capture the nonlinear, asymmetric, and tail dependence between the two variables.
The empirical results showed, first, that a significant dependence existed in all exchange rate-oil price pairs, and the dependence is negative, which means that oil serves as a hedge against rising inflation in BRICS countries. For example, when one currency appreciates, consumers from that country will find the oil price less expensive and then increase their demand for oil, which will lead to a rise in the oil price. Second, in both the constant and time-varying copula models, the RUB-WTI pair has the strongest dependence, possibly because Russia is a large oil producer and has a high share of global oil exports, making Russia itself more dependent on oil. Third, from the goodness-of-fit test results, the Student's t copula model has the best performance among all the constant copula models considered, while for the time-varying copula models, it was rejected for the BRL-WTI, INR-WTI, and ZAR-WTI pairs. We also measured risk management by treating the five exchange rate-oil price pairs as portfolios and calculating the dynamic values of VaR and ES. Both reached very low values for all pairs when the oil price experienced a sharp decrease.
The results of this paper offer us several suggestions. First, the significant correlation between oil price and exchange rate indicates that enterprises and governments in BRICS countries, when the oil price fluctuates, should pay attention to the exchange rate risk. Second, the negative relationship between oil price and exchange rate shows that the oil could be a hedge against inflation in BRICS countries, which is useful for foreign investors to manage exchange rate risk.
However, in this paper, we only considered bivariate copula models. For the future extension of this research, we will focus on a high-dimensional copula model, which could provide a more flexible analysis.

Acknowledgments:
We are grateful to the editor and three anonymous referees for many helpful comments and suggestions.

Conflicts of Interest:
The authors declare no conflict of interest. Table A1 summarizes the results of ARMA-GARCH models. We first use ARMA models to estimate the return series, and then use GARCH model to estimate the residuals from the ARMA models to obtain the standard residuals. We chose the ARMA model lags using the AIC criterion. All returns series are modeled by standard GARCH model at first, but the standard residuals in RUB and ZAR returns, could not pass the Ljung-Box test and Lagrange Multiplier test, which means the GARCH model is not fitted well in these returns. Hence, we chose GJR-GARCH model in fitting RUB and ZAR return data. As mentioned in the main text, each GARCH model tests the hypothesis that the error item follows a skewed t distribution. Almost all estimated parameters from the GARCH models are significant at the 1% level. Persistence coefficients (α + β in GARCH model, α + β + γ/2 in GJR-GARCH model) in all models are very close to 1, which reveal the existence of volatility clustering. This appears show that when the volatility is high, it is likely to remain high and vice versa when the volatility is low. This consequence is also in line with previous research. Table A2 provides the results of the Ljung-Box test of the standardized residuals and squared standard residuals (up to 20 lags), as well as the Lagrange Multiplier test of the standardized residuals. All values are insignificant, meaning that the standardized residuals obtained from the GARCH model have no serial correlation and autoregressive conditional heteroscedasticity. Note: this table presents the estimated parameters from the marginal distribution. ϕ, AR (p), and MA (q) are the estimated parameters from the ARMA (p, q) model. ω, α, β (, γ) are the estimated parameters from the GARCH (GJR-GARCH) model. λ, ν are the estimated parameters from skewed t model for the distribution of the error term. The persistence is calculated by α + β in GARCH model and by α + β + γ/2 in GJR-GARCH model. The values in parenthesis represent the standard errors of each parameter. *, **, and *** indicate significance at the 10%, 5%, and 1% levels, respectively.

Appendix B. Robustness Analysis
To examine the robustness of our results, we used the Brent crude oil price to analyze five currency-oil price pairs: BRL-BRENT, RUB-BRENT, INR-BRENT, CNY-BRENT, and ZAR-BRENT. We modeled the dependence structure for exchange rates and oil prices.
We summarize the rank correlation and linear correlation results, and the tail dependence from the constant Student's t copula model in Tables A3 and A4, respectively. We plot the estimated dynamic parameters and tail dependences in Figures A1 and A2, respectively, and illustrate the results for the dynamic VaR and ES Figure A3.
All results are consistent with those for the exchange rate-WTI pairs, which confirm the suitability of our proposed approach to capture the dependence structure between exchange rates and crude oil prices. Table A3. Rank correlation and linear correlation.