Skip to Main Content
Purpose

This study investigates systemic risk and risk spillovers among 21 commercial banks in Uzbekistan.

Design/methodology/approach

To this end, we adopt the Conditional Value-at-Risk framework and employ quantile regression, panel regression and a panel vector autoregression with exogenous variables (Panel VARX).

Findings

The empirical results show that banks with higher individual risk contribute disproportionately to system-wide risk, and that systemic risk surged sharply during the COVID-19 period before rapidly subsiding. Macroeconomic factors are found to play a more prominent role in driving systemic risk than bank-specific characteristics, highlighting the banking sector's vulnerability to external shocks. Dynamic analysis further indicates that global uncertainty exerts a stronger influence on risk spillovers than domestic real-sector shocks.

Practical implications

These findings underscore the importance of strengthening bank-level risk management, enhancing system-wide supervision and establishing a robust macroprudential policy framework to mitigate externally driven systemic risks in Uzbekistan's banking sector.

Originality/value

This paper contributes to the emerging-economy literature by providing the first comprehensive CoVaR-based assessment of systemic risk in the Uzbek banking sector.

This paper examines systemic risk in Uzbekistan's banking sector using the Conditional Value-at-Risk (CoVaR) framework proposed by Adrian and Brunnermeier (2016), one of the most widely used approaches for measuring systemic risk. Specifically, we analyze bank-level risk spillovers under extreme market conditions and empirically investigate the determinants of systemic risk.

The global financial crisis highlighted that the social costs of financial distress arise primarily from system-wide spillovers rather than from the failure of individual institutions. In response, a substantial body of literature has developed measures designed to capture the externalities imposed by financial institutions on the broader financial system. Among these, the CoVaR framework has emerged as a leading market-based measure of systemic risk.

CoVaR is defined as the value-at-risk (VaR) of the financial system conditional on a specific institution being in distress, while ΔCoVaR captures the marginal contribution of that institution to system-wide tail risk by comparing distress and median states (Adrian and Brunnermeier, 2016). The framework is typically implemented using quantile regression, which accommodates asymmetric and nonlinear dependence structures and avoids restrictive distributional assumptions. These features make CoVaR particularly well suited to financial data characterized by fat tails and state-dependent behavior.

Growing literature applies CoVaR-based measures to emerging and developing economies, where banking systems are often dominant, highly exposed to external shocks and subject to structural transformations. Across regions, studies consistently find that systemic risk is strongly time-varying and increases sharply during periods of financial stress (Girardi and Ergün, 2013).

Empirical evidence from emerging Asia demonstrates that CoVaR and ΔCoVaR are effective in identifying systemically important banks and tracking systemic risk dynamics over financial cycles. Studies on China and other Asian economies show that large, state-owned or highly leveraged banks contribute disproportionately to system-wide tail risk, particularly during crisis episodes (Jiang et al., 2020; Wen et al., 2020; Xu et al., 2021). Similar findings are reported for India and ASEAN economies, where quantile-based measures reveal strong interbank spillovers and heightened sensitivity to macro-financial conditions and external shocks (Verma et al., 2019; Khan et al., 2021; Pham et al., 2021). Evidence from other emerging economies further confirms that systemic risk is closely linked to bank characteristics such as size, leverage, capitalization and liquidity, as well as to global financial conditions (Rizwan and Ahmad, 2019; De Jonghe et al., 2015; Mendonça et al., 2019).

Overall, the emerging-economy literature highlights three stylized facts. First, systemic risk is strongly time-varying and rises nonlinearly during stress periods. Second, global financial conditions often play a dominant role in shaping domestic systemic risk, sometimes exceeding the explanatory power of bank-specific factors. Third, rankings of systemically important banks can vary across risk measures, underscoring the need for robustness analysis and complementary indicators.

Uzbekistan represents a distinctive case among emerging and transition economies. Its banking system has undergone rapid transformation in recent years, including financial liberalization, monetary policy reform and a strategic shift toward commercialization and privatization of state-owned banks. In this context, international institutions emphasize the need to strengthen system-wide risk monitoring and macroprudential surveillance (IMF, 2025). Consistent with this policy focus, the Central Bank of the Republic of Uzbekistan has highlighted systemic risk indicators conceptually aligned with the ΔCoVaR framework in its financial stability assessments (Central Bank of the Republic of Uzbekistan, 2022). However, despite extensive evidence for other emerging economies, bank-level empirical analysis of systemic risk in Uzbekistan remains scarce.

This paper fills this gap by providing the first comprehensive CoVaR-based assessment of systemic risk in Uzbekistan's banking sector. Using data for 21 commercial banks, we combine quantile regression, panel regressions and a PVARX framework to analyze risk spillovers and their determinants. The results show that banks with higher individual risk contribute disproportionately to systemic risk, that systemic risk surged during the COVID-19 period before subsiding and that macroeconomic and global financial factors dominate bank-specific characteristics. These findings have important implications for bank-level risk management and macroprudential policy in transition economies.

The remainder of the paper is organized as follows. Section 2 presents the static and dynamic CoVaR models and estimation results. Section 3 analyzes the determinants of systemic risk using panel regressions and PVARX models. Section 4 concludes with policy implications.

VaR refers to the minimum loss of a financial asset or portfolio at a given confidence level α, where X denotes an investment with stochastic returns, as defined in Equation (1).

(1)

VaR provides a concise and intuitive framework for assessing risk in complex investment portfolios composed of multiple financial instruments, and its theoretical foundations are relatively straightforward. However, extensive empirical research has identified several fundamental shortcomings of VaR. First, VaR may violate the principle of diversification as the risk of a portfolio can exceed the sum of risks calculated separately for individual assets. Second, VaR fails to capture information about extreme tail losses, which are often of greatest interest for risk management. Third, VaR typically relies on the assumption of normally distributed returns, which is unrealistic for financial time series that commonly exhibit volatility clustering and fat tails. Finally, VaR measures only institution-specific risk and does not account for spillover effects across financial institutions.

These limitations of VaR became particularly evident following the Global Financial Crisis (GFC) of 2008. To address these shortcomings, Adrian and Brunnermeier (2016) proposed the CoVaR, a risk measure that captures the tail risk of the financial system conditional on an institution being under distress. CoVaR is derived from the conditional distribution of system-wide losses beyond the VaR threshold. Compared with VaR, CoVaR provides a more informative framework for assessing systemic risk and portfolio optimization. If the loss of a financial asset exceeds its VaR at confidence level α, the conditional expectation is expressed in Equation (2).

(2)

When the stochastic loss of the portfolio is denoted by −X (−X<0), VaRα at the 1α quantile represents the loss level that is exceeded with probability 1α, as shown in Equation (3).

(3)

CoVaR optimization yields results consistent with VaR optimization under normal and elliptical distributions. However, unlike VaR, which may be discontinuous, CoVaR is a continuous and convex function with respect to portfolio positions, making it more suitable for systemic risk analysis. To examine spillover effects in a banking system consisting of multiple institutions, we specify the quantile regression model in Equation (4).

(4)

where Xˆqsystem|Xi denotes the estimated losses at q quantile when the loss of institution i is X under a certain confidence level. Based on the definition of VaR, CoVaR is obtained as the predicted system loss conditional on institution i’s loss reaching its VaR level, as shown in Equations (5) and (6). Once VaRqi is estimated from the return series of institution i, the risk spillover effect, measured by ΔCoVaR, can be computed as shown in Equation (7).

(5)
(6)
(7)

Given that financial time series typically follow leptokurtic rather than normal distributions, conventional linear regression methods are often inappropriate for risk estimation. Moreover, ordinary least squares (OLS) regression focuses on conditional means and thus provides limited insights into tail behavior. Quantile regression addresses these limitations by modeling the relationship between explanatory variables and specific quantiles of the dependent variable. As a result, quantile regression is particularly well suited for analyzing tail risk and extreme events. It is also more robust to outliers and does not impose restrictive assumptions such as homoscedasticity. These features make quantile regression widely used in financial economics and systemic risk analysis.

According to internal data from the Central Bank of Uzbekistan (CBU), 33 commercial banks are recorded. Among them, 12 banks (banks 4, 6, 10, 12, 15, 17, 19, 25, 30, 31, 32 and 33) are excluded due to data limitations and extreme outliers. Based on data availability and research feasibility, the final sample includes 21 banks over the period from March 2019 to March 2022, yielding 37 monthly observations.

Table 1 reports summary statistics for individual bank returns. All banks exhibit positive average returns. Although skewness varies substantially across banks, all return series display leptokurtic behavior. Except for Bank 20 and Bank 29, the Jarque–Bera test rejects the null hypothesis of normality at the 10% significance level for all banks, consistent with the typical characteristics of financial return series, such as fat tails and asymmetry. Due to the lack of stock price data, earnings before interest and taxes (EBIT) are used to compute individual bank returns. The return of the banking system is calculated as the weighted sum of individual bank returns, where weights are based on each bank's share of total equity. The book value of equity (BVE) is employed as a proxy for market value, and end-of-period data are used to maintain consistency.

Table 1

Summary statistics of individual bank returns

BankMeanMaxMinStd. devSkew-nessKurtosisJ-BADFN
10.181.09−0.890.42−0.794.360.040.0037
20.170.97−0.880.40−0.844.540.020.0037
30.396.66−0.951.184.1523.230.000.0037
50.278.35−2.811.514.0323.870.000.0037
70.303.77−0.980.752.4613.560.000.0037
80.262.91−0.930.651.679.290.000.0037
90.221.31−0.910.46−0.504.710.050.0037
110.161.15−0.840.39−0.865.230.000.0037
130.201.13−0.910.42−0.935.170.000.0037
140.172.22−0.900.560.886.400.000.0037
160.201.17−0.920.44−0.674.640.030.0037
180.231.14−0.930.44−0.954.780.010.0037
200.301.97−0.960.580.264.330.210.0037
210.4210.02−2.942.112.9713.760.000.0037
220.170.82−0.890.38−1.335.400.000.0037
230.301.93−0.960.610.514.560.070.0037
240.293.49−1.741.031.095.800.000.0037
260.201.34−0.940.43−0.775.450.000.0037
270.9141.18−17.728.193.2018.070.000.0037
280.1418.38−25.395.46−1.7317.300.000.0037
290.231.37−0.920.50−0.053.780.620.0037

Note(s): The p-values for J-B (Jarque–Bera) and ADFs (augmented Dickey–Fuller statistics) are reported

Source(s): Authors’ calculation

To quantify risk spillovers, we first estimate the quantile regression model specified in Equation (8). The estimated coefficients are then substituted into Equation (9) to compute CoVaR for each bank. Using these estimates, ΔCoVaR is obtained from Equation (10), which measures the marginal contribution of an institution to systemic risk by comparing distress and median states.

(8)
(9)
(10)

where VaR0.05Banki denotes Bank i's return value at 5% quantile, ΔCoVaR0.05SYS|Banki represents an institution's contribution to systemic risk (spillover effects) and CoVaR0.05SYS|Banki and CoVaR0.5SYS|Banki are CoVaR under distress and median (normal) circumstances, respectively.

The extent to which the banking system is affected by an individual bank is a reliable indicator of systemic risk. The impact is inversely correlated with the stability of the bank's own returns, which will affect the fluctuations in the aggregate yields of the banking sector. As seen in Table 2, Bank 2 is the largest, followed by Bank 11 and Bank 1, in terms of the absolute magnitude of risk spillovers from the specific bank to the banking system. There is little change in the ranking of unconditional VaR, which suggests that banks with higher risk will exert greater spillover effects on the banking sector, and vice versa. The positive association demonstrates that by enhancing and perfecting internal risk management, the whole banking sector can effectively withstand threats.

Table 2

Results of VaR, CoVaR and ΔCoVaR

BankVaR0.05SYS|BankiCoVaR0.05SYS|BankiΔCoVaR0.05SYS|Banki
2−3.58−4.29−2.95
11−3.25−3.98−2.88
1−3.52−4.19−2.65
26−3.02−3.80−2.44
22−2.96−3.66−2.39
13−2.78−3.54−2.30
9−3.05−3.71−2.22
20−2.81−3.50−2.14
8−2.89−3.62−2.08
18−2.68−3.37−1.97
23−2.45−3.20−1.88
29−2.83−3.51−1.85
14−2.66−3.42−1.80
16−2.76−3.47−1.77
21−2.38−3.05−1.65
3−2.24−2.93−1.60
5−1.95−2.65−1.51
24−2.12−2.85−1.50
7−2.18−2.78−1.41
28−1.96−2.72−1.22
27−1.91−2.69−1.09
Source(s): Authors’ calculation

For ease of interpretation, all values are multiplied by −1 and converted to positive numbers. Figure 1 illustrates that VaR and ΔCoVaR exhibit similar dynamics, indicating that both measures capture key aspects of systemic risk. Comparing these indicators allows for the identification of systemically important banks and provides valuable information for risk monitoring and early-warning systems.

Figure 1
A bar chart comparing V a R and Delta C o V a R across banks with dashed lines showing average V a R and Delta C o V a R.The horizontal axis lists banks labeled “Bank 2”, “Bank 11”, “Bank 1”, “Bank 26”, “Bank 22”, “Bank 13”, “Bank 9”, “Bank 20”, “Bank 8”, “Bank 18”, “Bank 23”, “Bank 29”, “Bank 14”, “Bank 16”, “Bank 21”, “Bank 3”, “Bank 5”, “Bank 24”, “Bank 7”, “Bank 28”, and “Bank 27”. The vertical axis ranges from 0 to 4 in increments of 1 unit. Each bank has two vertical bars representing “V a R” and “Delta C o V a R”. Two horizontal dashed lines represent the averages labeled “avg underscore V a R” and “avg underscore Delta C o V a R”. The line for “avg underscore V a R” is positioned near 2.65, and the line for “avg underscore Delta C o V a R” is positioned near 1.95. The data from each bar are as follows: For Bank 2: V a R 3.58; Delta C o V a R 2.95. For Bank 11: V a R 3.25; Delta C o V a R 2.88. For Bank 1: V a R 3.52; Delta C o V a R 2.65. For Bank 26: V a R 3.02; Delta C o V a R 2.44. For Bank 22: V a R 2.96; Delta C o V a R 2.39. For Bank 13: V a R 2.78; Delta C o V a R 2.30. For Bank 9: V a R 3.05; Delta C o V a R 2.22. For Bank 20: V a R 2.81; Delta C o V a R 2.14. For Bank 8: V a R 2.89; Delta C o V a R 2.08. For Bank 18: V a R 2.68; Delta C o V a R 1.97. For Bank 23: V a R 2.45; Delta C o V a R 1.88. For Bank 29: V a R 2.83; Delta C o V a R 1.85. For Bank 14: V a R 2.66; Delta C o V a R 1.80. For Bank 16: V a R 2.76; Delta C o V a R 1.77. For Bank 21: V a R 2.38; Delta C o V a R 1.65. For Bank 3: V a R 2.24; Delta C o V a R 1.60. For Bank 5: V a R 1.95; Delta C o V a R 1.51. For Bank 24: V a R 2.12; Delta C o V a R 1.50. For Bank 7: V a R 2.18; Delta C o V a R 1.41. For Bank 28: V a R 1.96; Delta C o V a R 1.22. For Bank 27: V a R 1.91; Delta C o V a R 1.09. Note: All numerical values are reported in absolute ones terms in Table 2.

Visualization of VaR, average VaR, ΔCoVaR and average ΔCoVaR. Source: Authors’ calculation

Figure 1
A bar chart comparing V a R and Delta C o V a R across banks with dashed lines showing average V a R and Delta C o V a R.The horizontal axis lists banks labeled “Bank 2”, “Bank 11”, “Bank 1”, “Bank 26”, “Bank 22”, “Bank 13”, “Bank 9”, “Bank 20”, “Bank 8”, “Bank 18”, “Bank 23”, “Bank 29”, “Bank 14”, “Bank 16”, “Bank 21”, “Bank 3”, “Bank 5”, “Bank 24”, “Bank 7”, “Bank 28”, and “Bank 27”. The vertical axis ranges from 0 to 4 in increments of 1 unit. Each bank has two vertical bars representing “V a R” and “Delta C o V a R”. Two horizontal dashed lines represent the averages labeled “avg underscore V a R” and “avg underscore Delta C o V a R”. The line for “avg underscore V a R” is positioned near 2.65, and the line for “avg underscore Delta C o V a R” is positioned near 1.95. The data from each bar are as follows: For Bank 2: V a R 3.58; Delta C o V a R 2.95. For Bank 11: V a R 3.25; Delta C o V a R 2.88. For Bank 1: V a R 3.52; Delta C o V a R 2.65. For Bank 26: V a R 3.02; Delta C o V a R 2.44. For Bank 22: V a R 2.96; Delta C o V a R 2.39. For Bank 13: V a R 2.78; Delta C o V a R 2.30. For Bank 9: V a R 3.05; Delta C o V a R 2.22. For Bank 20: V a R 2.81; Delta C o V a R 2.14. For Bank 8: V a R 2.89; Delta C o V a R 2.08. For Bank 18: V a R 2.68; Delta C o V a R 1.97. For Bank 23: V a R 2.45; Delta C o V a R 1.88. For Bank 29: V a R 2.83; Delta C o V a R 1.85. For Bank 14: V a R 2.66; Delta C o V a R 1.80. For Bank 16: V a R 2.76; Delta C o V a R 1.77. For Bank 21: V a R 2.38; Delta C o V a R 1.65. For Bank 3: V a R 2.24; Delta C o V a R 1.60. For Bank 5: V a R 1.95; Delta C o V a R 1.51. For Bank 24: V a R 2.12; Delta C o V a R 1.50. For Bank 7: V a R 2.18; Delta C o V a R 1.41. For Bank 28: V a R 1.96; Delta C o V a R 1.22. For Bank 27: V a R 1.91; Delta C o V a R 1.09. Note: All numerical values are reported in absolute ones terms in Table 2.

Visualization of VaR, average VaR, ΔCoVaR and average ΔCoVaR. Source: Authors’ calculation

Close modal

In practice, the risk faced by a financial institution is influenced not only by its own characteristics but also by shocks originating from other financial sectors. The transmission of such shocks is commonly referred to as the risk spillover effect. Given its ability to capture tail dependence, CoVaR is particularly suitable for analyzing dynamic systemic risk. To model time-varying risk spillovers, we estimate a quantile regression framework in which bank returns are regressed on lagged state variables, as shown in Equations (11) and (12).

(11)
(12)

where Xti represents the return rate of financial institution i at time t, q is the quantile and Mt1 denotes the state variable with a lag of one period. The parameters of the model can be estimated using quantile regression. The inclusion of lagged variables accounts for delayed transmission effects in financial markets. Based on the estimated coefficients, time-varying VaR and CoVaR series are computed using Equations (13) and (14). Subsequently, dynamic ΔCoVaR is calculated as shown in Equation (15).

(13)
(14)
(15)

Following the literature, extreme risk scenarios are captured using the 1 and 5% quantiles, while the 50% quantile represents normal market conditions. The same sample of 21 banks is used, and several banking-sector risk factors are included as state variables, as summarized in Table 3.

Table 3

List of state variables and description

NoState variableDescriptionTime span
1Liquidity spread (LS)First difference value of [weighted average interest rates on deposits (up to 30 days) – treasury bill rate (1 year)]2019M3∼ 2022M3
2Interest rate spread (IRS)First difference value of [weighted average interest rates on deposits (over 1 year) – weighted average interest rates on deposits (up to 30 days)]2019M3∼ 2022M3
3Credit spread (CS)First difference value of [1-year deposit rate – 1-year treasury bill rate]2019M3∼ 2022M3
4Stock market volatility index (VUCI)Conditional variance of Uzbekistan Composite index returns estimated by GARCH (1,1)2019M3∼ 2022M3
Source(s): Authors’ calculation using raw data

All state variables can measure the bank's systemic risk from different perspectives. We consider the liquidity spread (LS) as a key indicator of both the bank's credit policy and the market participants' capability to finance. The LS is the difference between the weighted average short-term interest rate and the 1-year treasury bill rate, which is known as the risk-free interest rate. The interest rate spread (IRS) is generally measured by the difference between long-term and short-term treasury bill rates, which are proxied by the weighted average long-term deposit rate and the weighted average short-term deposit rate, respectively. The credit spread (CS) is the difference in yields between government bills with the same maturity and bonds issued by commercial banks. In this study, the one-year deposit rate is used in place of corporate bond yields. IRS and CS reflect interest rate risk and credit risk in the market and exert a significant influence on banks' asset allocation and leverage ratio. During economic expansion (contraction), owing to decreasing (increasing) credit risk and credit spreads, overall risks are alleviated (intensified). To achieve a more efficient asset allocation, banks are expected to ramp up leverage and raise liabilities.

Uzbekistan's sole stock exchange, the Tashkent Stock Exchange (TSE), provides the basis for measuring domestic market volatility. Since an official volatility index is unavailable, we estimate the conditional variance of the Uzbekistan Composite Index (UCI) using a GARCH (1,1) model. As reported in Table 4, all state variables are stationary, with LS and IRS approximately normally distributed, while CS and VUCI exhibit leptokurtosis and asymmetry.

Table 4

Summary statistics of state variables

VariableMeanMaxMinStd. devSkewnessKurtosisJ-BADFN
LS−0.134.45−4.932.18−0.123.040.960.0037
IRS0.165.46−4.062.060.353.660.490.0037
CS0.003.58−4.341.33−0.096.410.000.0037
VUCI0.050.090.040.012.558.140.000.0437

Note(s): The p-values for J-B (Jarque–Bera) and ADFs (augmented Dickey–Fuller statistics) are reported

Source(s): Authors’ calculation

Figure 2 shows that systemic risk measured at the 1% quantile consistently exceeds that at the 5% quantile. Systemic risk rose sharply during 2020 and peaked in March 2021, reflecting the impact of the COVID-19 shock, before declining rapidly in subsequent months. More recently, systemic risk has shown an upward trend, approaching pre-crisis levels. Figure 3 indicates that while the magnitude of spillovers varies across banks, their dynamic patterns are largely synchronized, with some degree of heterogeneity.

Figure 2
A line graph showing Delta Co V a R at 1 percent and 5 percent from March 2019 to March 2022.The horizontal axis shows time labeled “Mar-19”, “Jun-19”, “Sep-19”, “Dec-19”, “Mar-20”, “Jun-20”, “Sep-20”, “Dec-20”, “Mar-21”, “Jun-21”, “Sep-21”, “Dec-21”, and “Mar-22” in increments of 3 months. The vertical axis ranges from 0.00 to 12.00 in increments of 4.00. A legend below the chart identifies two series: “Delta Conditional Value at Risk (1 percent)”, represented by a solid line, and “Delta Conditional Value at Risk (5 percent)”, represented by a dashed line. The solid line for “Delta Conditional Value at Risk (1 percent)” begins near (Mar-19, 4.6), rises to a peak around (Jun-19, 9.5), declines near (Sep-19, 3.5), then follows a zig zag pattern between Dec-19 and Sep-20 in the range of 5.00 to 8.50, increases again around (Mar-21, 12.0), drops sharply near (Jun-21, 2.7), rises again near (Sep-21, 9.5), and ends near (Mar-22, 7.6). The dashed line for “Delta Conditional Value at Risk (5 percent)” begins near (Mar-19, 4.0), rises to a peak around (Jun-19, 9.0), declines near (Sep-19, 3.0), then follows a zig zag pattern between Dec-19 and Sep-20 in the range of 5.50 to 7.80, increases again around (Mar-21, 11.3), drops sharply near (Jun-21, 2.6), rises again near (Sep-21, 9.0), and ends near (Mar-22, 7.2). Note: All numerical values are approximated.

Comparison of ΔCoVaR in different scenarios (q = 1% VS. q = 5%). Source: Authors’ calculation

Figure 2
A line graph showing Delta Co V a R at 1 percent and 5 percent from March 2019 to March 2022.The horizontal axis shows time labeled “Mar-19”, “Jun-19”, “Sep-19”, “Dec-19”, “Mar-20”, “Jun-20”, “Sep-20”, “Dec-20”, “Mar-21”, “Jun-21”, “Sep-21”, “Dec-21”, and “Mar-22” in increments of 3 months. The vertical axis ranges from 0.00 to 12.00 in increments of 4.00. A legend below the chart identifies two series: “Delta Conditional Value at Risk (1 percent)”, represented by a solid line, and “Delta Conditional Value at Risk (5 percent)”, represented by a dashed line. The solid line for “Delta Conditional Value at Risk (1 percent)” begins near (Mar-19, 4.6), rises to a peak around (Jun-19, 9.5), declines near (Sep-19, 3.5), then follows a zig zag pattern between Dec-19 and Sep-20 in the range of 5.00 to 8.50, increases again around (Mar-21, 12.0), drops sharply near (Jun-21, 2.7), rises again near (Sep-21, 9.5), and ends near (Mar-22, 7.6). The dashed line for “Delta Conditional Value at Risk (5 percent)” begins near (Mar-19, 4.0), rises to a peak around (Jun-19, 9.0), declines near (Sep-19, 3.0), then follows a zig zag pattern between Dec-19 and Sep-20 in the range of 5.50 to 7.80, increases again around (Mar-21, 11.3), drops sharply near (Jun-21, 2.6), rises again near (Sep-21, 9.0), and ends near (Mar-22, 7.2). Note: All numerical values are approximated.

Comparison of ΔCoVaR in different scenarios (q = 1% VS. q = 5%). Source: Authors’ calculation

Close modal
Figure 3
A multi-line graph showing Delta Co V a R series across quarterly periods from 2019 to 2022.The horizontal axis shows time from 2019 to 2022, divided into quarterly periods labeled 1st, 2nd, 3rd, and 4th for each year. The vertical axis is labeled Coefficient of variation and ranges from negative 2.5 to 0.5 in increments of 0.5. A legend below the chart lists multiple series labeled “Delta Co V a R 1”, “Delta Co V a R 14”, “Delta Co V a R 20”, “Delta Co V a R 23”, “Delta Co V a R 29”, “Delta Co V a R 5”, “Delta Co V a R 8”, “Delta Co V a R 28”, “Delta Co V a R 16”, “Delta Co V a R 21”, “Delta Co V a R 24”, “Delta Co V a R 26”, “Delta Co V a R 27”, “Delta Co V a R 9”, “Delta Co V a R 13”, “Delta Co V a R 18”, “Delta Co V a R 22”, “Delta Co V a R 2”, “Delta Co V a R 3”, “Delta Co V a R 7”, and “Delta Co V a R 11”, each represented by different solid lines. Across the graph, most series begin slightly below zero in the 1st quarter of 2019, followed by a sharp decline in the 2nd quarter of 2019, where several lines reach values between negative 1.5 and negative 2.0. The lines then rebound close to zero in the 3rd quarter of 2019. During 2020, the series fluctuate in a zigzag pattern, generally ranging between negative 1.5 and negative 0.2. A pronounced drop occurs again in the 1st quarter of 2021, where several lines reach approximately negative 2.3 to negative 2.4, followed by a rebound toward zero in the 2nd quarter of 2021. In later quarters of 2021 and into 2022, the lines fluctuate moderately, typically between negative 1.2 and negative 0.2. Note: All numerical values are approximated.

Comparison of banks' dynamic systemic risk. Source: Authors’ calculation

Figure 3
A multi-line graph showing Delta Co V a R series across quarterly periods from 2019 to 2022.The horizontal axis shows time from 2019 to 2022, divided into quarterly periods labeled 1st, 2nd, 3rd, and 4th for each year. The vertical axis is labeled Coefficient of variation and ranges from negative 2.5 to 0.5 in increments of 0.5. A legend below the chart lists multiple series labeled “Delta Co V a R 1”, “Delta Co V a R 14”, “Delta Co V a R 20”, “Delta Co V a R 23”, “Delta Co V a R 29”, “Delta Co V a R 5”, “Delta Co V a R 8”, “Delta Co V a R 28”, “Delta Co V a R 16”, “Delta Co V a R 21”, “Delta Co V a R 24”, “Delta Co V a R 26”, “Delta Co V a R 27”, “Delta Co V a R 9”, “Delta Co V a R 13”, “Delta Co V a R 18”, “Delta Co V a R 22”, “Delta Co V a R 2”, “Delta Co V a R 3”, “Delta Co V a R 7”, and “Delta Co V a R 11”, each represented by different solid lines. Across the graph, most series begin slightly below zero in the 1st quarter of 2019, followed by a sharp decline in the 2nd quarter of 2019, where several lines reach values between negative 1.5 and negative 2.0. The lines then rebound close to zero in the 3rd quarter of 2019. During 2020, the series fluctuate in a zigzag pattern, generally ranging between negative 1.5 and negative 0.2. A pronounced drop occurs again in the 1st quarter of 2021, where several lines reach approximately negative 2.3 to negative 2.4, followed by a rebound toward zero in the 2nd quarter of 2021. In later quarters of 2021 and into 2022, the lines fluctuate moderately, typically between negative 1.2 and negative 0.2. Note: All numerical values are approximated.

Comparison of banks' dynamic systemic risk. Source: Authors’ calculation

Close modal

To examine the determinants of systemic risk, we estimate the panel regression model specified in Equation (16). The lagged dependent variable is included to account for persistence and potential endogeneity in systemic risk measures.

(16)

Book value of equity (BVE), total assets (AST) and debt-to-equity ratio (DTE) are bank-specific financial variables represented by Banki,t1. Fint denotes financial variables including the CBOE Volatility Index (VIX), Federal Funds Rate (FED) and Uzbekistan Stock Market Volatility Index (VUCI). Besides, the BVE and AST are log-differenced, the VIX and FED are first differenced series and the DTE is original data that was obtained from the CBU.

The specific explanations on independent variables can be summarized as follows. First, the BVE is vital because it provides a snapshot of a bank's financial state at a specific point in time. It can be used to compare a bank's current equity position to its equity position in previous periods. Additionally, it can be used to evaluate the financial leverage of a bank. Second, the AST is defined as the assets owned by a bank that has an economic value whose benefits can be derived in the future. To improve market value and sustainability for the future, the bank should look healthy and a bank's health will be decided on various parameters, among which assets are the most crucial ones as it will help in predicting the range of profit a bank can earn on its current investment over the period. Third, the DTE shows the level of financial leverage being used, reflecting the ability of owners' equity to protect creditors' equity and the bank's long-term solvency. The three financial indicators listed above and the degree of risk are closely related. Fourth, VIX is a popular measure of the stock market's expectation of volatility based on S&P 500 index options, which is often referred to as the fear index. The influence of the US financial market on the global market, especially in emerging markets economies (EMEs), cannot be ignored. Fifth, FED is the target interest rate set by the Federal Open Market Committee (FOMC). FED is the representative overnight call rate, which sensitively reflects the surplus and shortage of excess reserves in the interbank market. It affects monetary and financial conditions, which in turn have a bearing on critical aspects of the broader economy including employment, growth and inflation. And finally, VUCI is an indicator capturing the uncertainty of domestic equity market returns.

Table 5 shows summary statistics for the variables. Applying the Levin et al. (2002) and Im et al. (2003) tests to check if panel datasets are stationary demonstrates that all series except for DTE reject the null hypothesis at the 1% level. Series VIX, FED and VUCI are also stationary.

Table 5

Summary statistics

VariableMeanStd. devSkewnessKurtosisJ-BLLCIPSN
ΔCoVaR−0.680.40−0.874.210.000.000.00756
BVE0.020.073.6574.570.000.000.00756
AST0.020.06−1.4823.790.000.000.00756
DTE4.712.55−0.252.650.000.360.13756
VIX0.010.300.152.680.03756
FED−0.060.26−4.7526.490.00756
VUCI0.050.012.408.520.00756

Note(s): The p-values for J-B (Jarque–Bera), LLC (Levin–Lin–Chu) and IPS (Im–Pesaran–Shin) are reported

Source(s): Authors’ calculation

Table 6 shows the lag terms of ΔCoVaR are significantly positive at a statistical level of 5% or 10%, which accords with the autocorrelation of most financial time series. In other words, the current level of risk spillovers is positively correlated with the level of risk spillovers in the past period.

Table 6

Empirical results of panel regression with fixed effects

Variable(1)(2)(3)(4)(5)(6)
CoVaR (−1)0.032** (2.121)0.030* (2.014)0.043** (2.708)0.041** (2.597)0.031** (2.539)0.029** (2.496)
BVE(−1)−0.231 (−1.636) −0.279 * (−1.843) −0.232 (−1.642) 
AST(−1) −0.080 (−0.382) −0.084 (−0.390) −0.136 (−0.779)
DTE(−1)0.007* (2.027)0.009** (2.671)0.008** (2.279)0.011*** (2.923)0.014** (2.250)0.017** (2.234)
VIX  0.201*** (4.688)0.197*** (4.741)0.218*** (4.878)0.215*** (4.916)
FED  −0.055*** (−2.945)−0.056*** (−3.009)−0.075*** (−3.681)−0.078*** (−3.785)
VUCI    9.960*** (7.551)10.080*** (7.706)
Cons−0.687*** (−35.092)−0.704*** (−34.075)−0.688*** (−32.992)−0.708*** (−32.939)−1.136*** (−20.009)−1.159*** (−21.022)
N756756756756756756

Note(s): (1) The Hausman test results choose fixed effects

(2) ***, ** and * denote significance level at 1%, 5% and 10%, respectively

(3) t-values are in the parentheses

For the banks with steady earnings and cash flows, a greater DTE is preferred since it uses leverage to boost return on equity (ROE) and lower weighted average cost of capital. However, if the DTE is excessively high, the bank may not be able to repay its debts or even go bankrupt, which drives up risk spillovers. The estimated coefficients consistently demonstrate that a larger DTE is linked to a higher systemic risk, despite their differing degrees of significance.

Based on baseline models (1) and (2), we add two global financial variables to models (3) and (4). Global risks are bound to be transmitted to the domestic financial market and commodity market, which is why the coefficient of VIX is positive and significant. The change in interest rates will have an impact on the bank's lending and investment activities, as well as risk spillovers. For instance, a rise (fall) in interest rate will boost (curb) the cost of borrowing for businesses, and from the standpoint of banks, it will increase (decrease) the yield, which could improve (deteriorate) risk-taking capacity. The significance and sign of the FED coefficients are as expected. Domestic equities market volatility has a positive impact on ΔCoVaR in models (5) and (6) as VIX does.

Regarding the short panel (known as the large-N and small-T case), it is difficult to determine whether the error term exhibits autocorrelation in that the information of each unit is insufficient. However, the error term is assumed to be independent and identically distributed (IID). For a long panel (large-T and N > T case), the assumption can be loosened. Meanwhile, given the possible in-group autocorrelation, between-group heteroscedasticity and contemporaneous correlation, we employ the panel generalized least squares (GLS) as a robustness check. Table 7 substantiates that the results are robust since most variables are slightly different in significance level and magnitude. While the sign of FED is as expected, it is no longer significant.

Table 7

Empirical results of panel GLS

Variable(1)(2)(3)(4)(5)(6)
CoVaR (−1)0.334*** (9.534)0.332*** (9.447)0.308*** (8.806)0.305*** (8.664)0.282*** (8.271)0.282*** (8.226)
BVE(−1)−0.224 (−1.557) −0.267* (−1.875) −0.237* (−1.751) 
AST(−1) −0.049 (−0.282) −0.075 (−0.451) −0.179 (−1.115)
DTE(−1)0.008* (1.650)0.008* (1.726)0.009* (1.838)0.009* (1.937)0.009** (2.000)0.010** (2.152)
VIX  0.151*** (4.723)0.145*** (4.605)0.167*** (5.394)0.161*** (5.295)
FED  −0.042 (−1.111)−0.046 (−1.214)−0.055 (−1.507)−0.059 (−1.623)
VUCI    9.960*** (7.551)8.878*** (6.664)
Cons−0.478*** (−12.939)−0.485*** (−13.053)−0.504*** (−13.398)−0.513*** (−13.586)−0.952*** (−12.540)−0.966*** (−12.769)
N756756756756756756

Note(s): (1) The Hausman test results choose fixed effects

(2) ***, ** and * denote significance level at 1%, 5% and 10%, respectively

(3) t-values are in the parentheses

Additionally, we estimate a three-variable panel vector autoregression (Panel VAR) model to examine the dynamic relationship among several key variables and pave the way for providing policy implications as ΔCoVaR is significantly affected by the VIX and FED, both of which are exogeneous variables within the system. Therefore, we can estimate an extended model by adding lagged endogenous and exogenous simultaneously, Panel VAR with exogenous variables, i.e. Panel VARX as in Equation (17).

(17)

As a nontheoretical model, the VAR framework does not require any prior constraints on variables. Frequently, we are more interested in the dynamics that are predicted by VAR models than the actual coefficients that are estimated. For this reason, it is most common that VAR studies report impulse response functions (IRFs) and forecast error variance decomposition (FEVD). The IRFs are conducted to trace the dynamic path of variables in the system in response to shocks to other variables, focusing more on the dynamic impact on the system when an error term changes or the model is affected by a certain shock, rather than analyzing the impact of a change in one variable on another variable. A common approach to identifying the shocks of a VAR model is to employ orthogonal impulse responses (OIRs). The basic idea is to decompose the variance–covariance matrix so that Σ = PP’, where P is a lower triangular matrix with positive diagonal elements, which is often obtained by a Cholesky decomposition. Note that the results of an OIR might be sensitive to the order of the variables. We specify Cholesky ordering of endogenous variables as follows: VUCI, DTE and ΔCoVaR, as the volatility of equity market could affect the equity-to-debt ratio. In addition, Granger causality tests show that DTE does help predict VUCI. Therefore, we subsequently adjust the order of endogenous variables as a robustness check. FEVD separates the forecast error variance into proportions attributed to each variable in the system. To be more specific, FEVD further evaluates the importance of different structural shocks by analyzing the contribution of each structural shock to the change of endogenous variables (usually measured by variance).

Figure 4 presents that the positive shock of DTE will be transmitted to ΔCoVaR, which will cause the systemic risk to rise in the initial stage. The magnitude reaches a peak in the first period and remains stable until the second period before it starts to fall. The shock only lasts for three months and then becomes insignificant. When an external shock hits the VUCI, ΔCoVaR will ramp up dramatically in the next month and subsequently moves towards the horizontal line. Albeit statistically significant, the shock does not last for long. The results are robust, regardless of the ordering of endogenous variables as illustrated in Figure 5.

Figure 4
A 3 by 3 grid of panels showing “I R F” with shaded “95 percent C I” for delta, d t e, and v u c i across steps 0 to 10..The figure contains 9 impulse response function panels arranged in a 3 by 3 grid. Each panel displays a line labeled “I R F” with a surrounding shaded region representing the “95 percent C I”. The horizontal axis in all panels is labeled “step” and ranges from 0 to 10. A legend below the panels shows a shaded box labeled “95 percent C I” and a line labeled “I R F”. The bottom left text reads “impulse : response”. First panel positioned at left in first row 1, titled “delta : delta”, the vertical axis ranges approximately from 0 to 1. The line begins near 1.0 at step 0, drops sharply to about 0.15 at step 1, approaches approximately 0.0 by step 2, and remains close to zero through step 10. A narrow shaded confidence interval surrounds the line and is slightly wider during the early steps. Second panel positioned at center in first row 1, titled “delta : d t e”, the vertical axis ranges approximately from minus 0.4 to 0.2. The line begins near 0.0 at step 0, decreases to around minus 0.10 at step 1, and gradually increases toward approximately 0.0 by step 10. A wide shaded confidence interval surrounds the line across most steps. Third panel positioned at right in first row 1, titled “delta : v u c i”, the vertical axis ranges approximately from minus 0.006 to 0. The line begins at 0.0, drops to about minus 0.004 near step 1, and gradually rises toward approximately 0.0 by step 10. A shaded confidence interval surrounds the line and is wider in the early steps. Fourth panel positioned at left in second row 2, titled “d t e : delta”, the vertical axis ranges approximately from minus 0.1 to 0.1. The line begins near 0.0 at step 0, rises to about 0.06 around step 2, and gradually declines toward approximately 0.02 by step 10. The shaded confidence interval widens in later steps. Fifth panel positioned at center in second row 2, titled “d t e : d t e”, the vertical axis ranges approximately from minus 1 to 2. The line begins near 1.0 at step 0 and gradually declines toward approximately 0.15 by step 10. A broad shaded confidence interval surrounds the line and becomes wider as the steps increase. Sixth panel positioned at right in second row 2, titled “d t e : v u c i”, the vertical axis ranges approximately from minus 0.005 to 0.005. The line begins near 0.0, increases slightly to about 0.003 near step 2, and gradually decreases toward approximately 0.001 by step 10. A relatively narrow shaded confidence interval surrounds the line. Seventh panel positioned at left in third row 3, titled “v u c i : delta”, the vertical axis ranges approximately from 0 to 4. The line begins at 0.0, rises sharply to about 2.2 near step 1, then steadily declines toward approximately 0.0 by step 10. A wide shaded confidence interval surrounds the line and narrows gradually in later steps. Eighth panel positioned at center in third row 3, titled “v u c i : d t e”, the vertical axis ranges approximately from minus 10 to 5. The line begins near 0.0 at step 0, decreases to about minus 1.5 around steps 3 to 4, and gradually increases toward approximately minus 0.5 by step 10. A wide shaded confidence interval surrounds the line across the panel. Ninth panel positioned at right in third row 3, titled “v u c i : v u c i”, the vertical axis ranges approximately from 0 to 1. The line begins near 1.0 at step 0, declines to about 0.6 at step 1, continues decreasing toward approximately 0.1 by step 5, and approaches 0.0 by step 10. A narrow shaded confidence interval surrounds the line and gradually narrows as the values approach zero. Note: All numerical values are approximated.

Impulse response functions (1). Source: Authors’ calculation

Figure 4
A 3 by 3 grid of panels showing “I R F” with shaded “95 percent C I” for delta, d t e, and v u c i across steps 0 to 10..The figure contains 9 impulse response function panels arranged in a 3 by 3 grid. Each panel displays a line labeled “I R F” with a surrounding shaded region representing the “95 percent C I”. The horizontal axis in all panels is labeled “step” and ranges from 0 to 10. A legend below the panels shows a shaded box labeled “95 percent C I” and a line labeled “I R F”. The bottom left text reads “impulse : response”. First panel positioned at left in first row 1, titled “delta : delta”, the vertical axis ranges approximately from 0 to 1. The line begins near 1.0 at step 0, drops sharply to about 0.15 at step 1, approaches approximately 0.0 by step 2, and remains close to zero through step 10. A narrow shaded confidence interval surrounds the line and is slightly wider during the early steps. Second panel positioned at center in first row 1, titled “delta : d t e”, the vertical axis ranges approximately from minus 0.4 to 0.2. The line begins near 0.0 at step 0, decreases to around minus 0.10 at step 1, and gradually increases toward approximately 0.0 by step 10. A wide shaded confidence interval surrounds the line across most steps. Third panel positioned at right in first row 1, titled “delta : v u c i”, the vertical axis ranges approximately from minus 0.006 to 0. The line begins at 0.0, drops to about minus 0.004 near step 1, and gradually rises toward approximately 0.0 by step 10. A shaded confidence interval surrounds the line and is wider in the early steps. Fourth panel positioned at left in second row 2, titled “d t e : delta”, the vertical axis ranges approximately from minus 0.1 to 0.1. The line begins near 0.0 at step 0, rises to about 0.06 around step 2, and gradually declines toward approximately 0.02 by step 10. The shaded confidence interval widens in later steps. Fifth panel positioned at center in second row 2, titled “d t e : d t e”, the vertical axis ranges approximately from minus 1 to 2. The line begins near 1.0 at step 0 and gradually declines toward approximately 0.15 by step 10. A broad shaded confidence interval surrounds the line and becomes wider as the steps increase. Sixth panel positioned at right in second row 2, titled “d t e : v u c i”, the vertical axis ranges approximately from minus 0.005 to 0.005. The line begins near 0.0, increases slightly to about 0.003 near step 2, and gradually decreases toward approximately 0.001 by step 10. A relatively narrow shaded confidence interval surrounds the line. Seventh panel positioned at left in third row 3, titled “v u c i : delta”, the vertical axis ranges approximately from 0 to 4. The line begins at 0.0, rises sharply to about 2.2 near step 1, then steadily declines toward approximately 0.0 by step 10. A wide shaded confidence interval surrounds the line and narrows gradually in later steps. Eighth panel positioned at center in third row 3, titled “v u c i : d t e”, the vertical axis ranges approximately from minus 10 to 5. The line begins near 0.0 at step 0, decreases to about minus 1.5 around steps 3 to 4, and gradually increases toward approximately minus 0.5 by step 10. A wide shaded confidence interval surrounds the line across the panel. Ninth panel positioned at right in third row 3, titled “v u c i : v u c i”, the vertical axis ranges approximately from 0 to 1. The line begins near 1.0 at step 0, declines to about 0.6 at step 1, continues decreasing toward approximately 0.1 by step 5, and approaches 0.0 by step 10. A narrow shaded confidence interval surrounds the line and gradually narrows as the values approach zero. Note: All numerical values are approximated.

Impulse response functions (1). Source: Authors’ calculation

Close modal
Figure 5
A 3 by 3 grid of impulse response panels with “I R F” and shaded “95 percent C I” for delta, d t e, and v u c i.The figure contains 9 impulse response function panels arranged in a 3 by 3 grid. Each panel displays a line labeled “I R F” with a surrounding shaded region representing the “95 percent C I”. The horizontal axis in all panels is labeled “step” and ranges from 0 to 10. A legend below the panels shows a shaded box labeled “95 percent C I” and a line labeled “I R F”. The bottom left text reads “impulse : response”. A legend below the panels shows a shaded box labeled “95 percent C I” and a line labeled “I R F”. The bottom left text reads “impulse : response”. First panel positioned at left in first row 1, titled “delta : delta”, the vertical axis ranges approximately from 0 to 1 in increments of 0.5. The line begins near 1.0 at step 0, drops sharply to about 0.15 at step 1, approaches approximately 0.0 by step 2, and remains close to zero through step 10. A narrow shaded confidence interval surrounds the line and is slightly wider during the early steps. Second panel positioned at center in first row 1, titled “delta : v u c i”, the vertical axis ranges approximately from minus 0.006 to 0.002 in increments of 0.002. The line begins at 0.0, drops to about minus 0.004 near step 1, and gradually rises toward approximately 0.0 by step 10. A shaded confidence interval surrounds the line and is wider in the early steps. Third panel positioned at right in first row 1, titled “delta : d t e”, the vertical axis ranges approximately from minus 0.5 to 0.5 in increments of 0.5. The line begins near 0.0 at step 0, decreases to around minus 0.10 at step 1, and gradually increases toward approximately 0.0 by step 10. A wide shaded confidence interval surrounds the line across most steps. Fourth panel positioned at left in second row 2, titled “v u c i : delta”, the vertical axis ranges approximately from minus 2 to 4 in increments of 2. The line begins at 0.0, rises sharply to about 2.2 near step 1, then steadily declines toward approximately 0.0 by step 10. A wide shaded confidence interval surrounds the line and narrows gradually in later steps. Fifth panel positioned at center in second row 2, titled “v u c i : v u c i”, the vertical axis ranges approximately from 0 to 1 in increments of 0.5. The line begins near 1.0 at step 0, declines to about 0.6 at step 1, continues decreasing toward approximately 0.1 by step 5, and approaches 0.0 by step 10. A narrow shaded confidence interval surrounds the line and gradually narrows as the values approach zero. Sixth panel positioned at right in second row 2, titled “v u c i : d t e”, the vertical axis ranges approximately from minus 20 to 10 in increments of 10. The line begins near 0.0 at step 0, decreases slightly below zero in the early steps, and gradually approaches approximately 0.0 by step 10. A wide shaded confidence interval surrounds the line across the panel. Seventh panel positioned at left in third row 3, titled “d t e : delta”, the vertical axis ranges approximately from minus 0.1 to 0.2 in increments of 0.1. The line begins near 0.0 at step 0, rises to about 0.06 around step 2, and gradually declines toward approximately 0.02 by step 10. The shaded confidence interval widens in later steps. Eighth panel positioned at center in third row 3, titled “d t e : v u c i”, the vertical axis ranges approximately from minus 0.005 to 0.01 in increments of 0.005. The line begins near 0.0, increases slightly to about 0.003 near step 2, and gradually decreases toward approximately 0.001 by step 10. A relatively narrow shaded confidence interval surrounds the line. Ninth panel positioned at right in third row 3, titled “d t e : d t e”, the vertical axis ranges approximately from minus 2 to 2 in increments of 1. The line begins near 1.0 at step 0 and gradually declines toward approximately 0.15 by step 10. A broad shaded confidence interval surrounds the line and becomes wider as the steps increase. Note: All numerical values are approximated.

Impulse response functions (2). Source: Authors’ calculation

Figure 5
A 3 by 3 grid of impulse response panels with “I R F” and shaded “95 percent C I” for delta, d t e, and v u c i.The figure contains 9 impulse response function panels arranged in a 3 by 3 grid. Each panel displays a line labeled “I R F” with a surrounding shaded region representing the “95 percent C I”. The horizontal axis in all panels is labeled “step” and ranges from 0 to 10. A legend below the panels shows a shaded box labeled “95 percent C I” and a line labeled “I R F”. The bottom left text reads “impulse : response”. A legend below the panels shows a shaded box labeled “95 percent C I” and a line labeled “I R F”. The bottom left text reads “impulse : response”. First panel positioned at left in first row 1, titled “delta : delta”, the vertical axis ranges approximately from 0 to 1 in increments of 0.5. The line begins near 1.0 at step 0, drops sharply to about 0.15 at step 1, approaches approximately 0.0 by step 2, and remains close to zero through step 10. A narrow shaded confidence interval surrounds the line and is slightly wider during the early steps. Second panel positioned at center in first row 1, titled “delta : v u c i”, the vertical axis ranges approximately from minus 0.006 to 0.002 in increments of 0.002. The line begins at 0.0, drops to about minus 0.004 near step 1, and gradually rises toward approximately 0.0 by step 10. A shaded confidence interval surrounds the line and is wider in the early steps. Third panel positioned at right in first row 1, titled “delta : d t e”, the vertical axis ranges approximately from minus 0.5 to 0.5 in increments of 0.5. The line begins near 0.0 at step 0, decreases to around minus 0.10 at step 1, and gradually increases toward approximately 0.0 by step 10. A wide shaded confidence interval surrounds the line across most steps. Fourth panel positioned at left in second row 2, titled “v u c i : delta”, the vertical axis ranges approximately from minus 2 to 4 in increments of 2. The line begins at 0.0, rises sharply to about 2.2 near step 1, then steadily declines toward approximately 0.0 by step 10. A wide shaded confidence interval surrounds the line and narrows gradually in later steps. Fifth panel positioned at center in second row 2, titled “v u c i : v u c i”, the vertical axis ranges approximately from 0 to 1 in increments of 0.5. The line begins near 1.0 at step 0, declines to about 0.6 at step 1, continues decreasing toward approximately 0.1 by step 5, and approaches 0.0 by step 10. A narrow shaded confidence interval surrounds the line and gradually narrows as the values approach zero. Sixth panel positioned at right in second row 2, titled “v u c i : d t e”, the vertical axis ranges approximately from minus 20 to 10 in increments of 10. The line begins near 0.0 at step 0, decreases slightly below zero in the early steps, and gradually approaches approximately 0.0 by step 10. A wide shaded confidence interval surrounds the line across the panel. Seventh panel positioned at left in third row 3, titled “d t e : delta”, the vertical axis ranges approximately from minus 0.1 to 0.2 in increments of 0.1. The line begins near 0.0 at step 0, rises to about 0.06 around step 2, and gradually declines toward approximately 0.02 by step 10. The shaded confidence interval widens in later steps. Eighth panel positioned at center in third row 3, titled “d t e : v u c i”, the vertical axis ranges approximately from minus 0.005 to 0.01 in increments of 0.005. The line begins near 0.0, increases slightly to about 0.003 near step 2, and gradually decreases toward approximately 0.001 by step 10. A relatively narrow shaded confidence interval surrounds the line. Ninth panel positioned at right in third row 3, titled “d t e : d t e”, the vertical axis ranges approximately from minus 2 to 2 in increments of 1. The line begins near 1.0 at step 0 and gradually declines toward approximately 0.15 by step 10. A broad shaded confidence interval surrounds the line and becomes wider as the steps increase. Note: All numerical values are approximated.

Impulse response functions (2). Source: Authors’ calculation

Close modal

Table 8 shows the contribution of each variable to ΔCoVaR variation. At the beginning, approximately 80% of the variation in the ΔCoVaR is from shocks to ΔCoVaR itself, and most of the remaining 20% is from VUCI. The contribution of DTE to the variation in the ΔCoVaR changes slowly over time and eventually seems to converge at around 5.4%. Table 9 reveals that the adjustment of endogenous variables has little impact on the results. The system becomes stable after one year, with the contribution of ΔCoVaR, DTE and VUCI to variation in ΔCoVaR leveling off at roughly 74%, 5% and 21%, respectively.

Table 8

Forecast error variance decomposition of ΔCoVaR (1)

Period (month)CoVaRDTEVUCI
179.840.3719.79
278.770.8920.34
377.741.5520.71
476.981.9621.06
576.422.2521.33
675.033.3321.64
774.753.5121.74
873.564.3622.08
972.434.8922.68
1072.344.9522.71
1171.285.4023.32
1270.245.4424.32
Source(s): Authors’ calculation
Table 9

Forecast error variance decomposition of ΔCoVaR (2)

Period (month)CoVaRDTEVUCI
179.790.4919.72
278.991.0120.00
378.011.6320.36
477.601.9820.42
577.042.4420.52
676.322.9520.73
775.753.4820.77
875.233.9120.86
974.734.3720.90
1074.074.8921.04
1173.895.0221.09
1273.665.2221.12
Source(s): Authors’ calculation

This study empirically analyzes risk spillovers among 21 commercial banks in Uzbekistan using the CoVaR framework combined with quantile regression, panel regression and PVARX models. The main findings can be summarized as follows. First, static CoVaR results show that banks with higher individual risk generate larger spillover effects on the banking system, as reflected in the close correspondence between VaR and ΔCoVaR rankings. This implies that strengthening internal risk management at the bank level can effectively enhance system-wide resilience. Dynamic CoVaR estimates further indicate that systemic risk surged sharply during the COVID-19 period, peaking in March 2021, before declining rapidly in subsequent months. Risk spillovers across banks also exhibit broadly synchronized dynamics over the sample period.

Second, panel regression results indicate that macroeconomic and global financial variables exert a stronger influence on systemic risk than bank-specific financial indicators. This finding reflects the heightened vulnerability of Uzbekistan's banking sector to external shocks. Finally, PVARX analysis using impulse response functions and forecast error variance decomposition reveals that domestic equity market volatility contributes more to variations in systemic risk than bank leverage, reinforcing the view that Uzbekistan's banking system is highly exposed to global financial conditions.

Based on these findings, several policy implications emerge. Individual banks should prioritize reducing their risk exposure through stronger internal governance and more comprehensive risk measurement systems. At the system level, risk management and supervisory frameworks should be further strengthened. Many Uzbek commercial banks exhibit volatile earnings and high leverage, while asset quality and profitability face increasing pressure from elevated inflation and slowing economic growth. Establishing dedicated risk management institutions and developing multiple risk indicators, rather than relying on a single metric, would improve resilience to diverse financial risks. In addition, a nationwide system for risk monitoring and financial supervision is advisable, particularly given the systemic importance of large banks.

Overall, the results indicate that Uzbekistan's banking sector remains vulnerable to external shocks. Accordingly, the timely establishment of an effective macroprudential policy framework and the proactive design of policy tools to address external systemic risks are essential for safeguarding financial stability.

Earlier version of this manuscript was prepared as part of the final report for the Bank of Korea’s 2022 Knowledge Partnership Program.

Adrian
,
T.
and
Brunnermeier
,
M.K.
(
2016
), “
CoVaR
”,
The American Economic Review
, Vol. 
106
No. 
7
, pp. 
1705
-
1741
, doi: .
Central Bank of the Republic of Uzbekistan
(
2022
),
Financial Stability Review
,
Central Bank of the Republic of Uzbekistan
,
Tashkent
.
De Jonghe
,
O.
,
Diepstraten
,
M.
and
Schepens
,
G.
(
2015
), “
Banks' size, scope and systemic risk: what role for CoVaR?
”,
Journal of Banking and Finance
, Vol. 
55
, pp. 
197
-
215
.
Girardi
,
G.
and
Ergün
,
A.T.
(
2013
), “
Systemic risk measurement: multivariate GARCH estimation of CoVaR
”,
Journal of Banking and Finance
, Vol. 
37
No. 
8
, pp. 
3169
-
3180
, doi: .
Im
,
K.S.
,
Hashem Pesaran
,
M.
and
Shin
,
Y.
(
2003
), “
Testing for unit roots in heterogeneous panels
”,
Journal of Econometrics
, Vol. 
115
No. 
1
, pp.
53
-
74
, doi: .
IMF
(
2025
),
Uzbekistan: Financial Sector Assessment Program (FSAP)
,
International Monetary Fund
,
Washington, DC
.
Jiang
,
L.
,
Yang
,
Y.
and
Li
,
S.
(
2020
), “
Measuring systemic risk in China’s banking sector: a CoVaR approach
”,
Finance Research Letters
, Vol. 
33
, 101192.
Khan
,
S.
,
Ahmad
,
N.
and
Abbas
,
S.
(
2021
), “
Measuring systemic risk in Indian banks: a quantile-based approach
”,
International Journal of Finance and Economics
, Vol. 
26
No. 
3
, pp. 
4201
-
4218
.
Levin
,
A.
,
Lin
,
C.-F.
and
Chu
,
C.-S.J.
(
2002
), “
Unit root tests in panel data: asymptotic and finite-sample properties
”,
Journal of Econometrics
, Vol. 
108
No. 
1
, pp.
1
-
24
, doi: .
Mendonça
,
H.F.
,
Galvão
,
D.J.C.
and
Loures
,
R.F.V.
(
2019
), “
Determinants of systemic risk in Brazil: a CoVaR approach
”,
Economic Modelling
, Vol. 
77
, pp. 
254
-
265
.
Pham
,
T.T.T.
,
Nguyen
,
D.T.
and
Vo
,
X.V.
(
2021
), “
Systemic risk and bank characteristics in ASEAN banking systems: evidence from CoVaR
”,
Borsa Istanbul Review
, Vol. 
21
No. 
4
, pp. 
329
-
343
.
Rizwan
,
M.S.
and
Ahmad
,
M.
(
2019
), “
Measuring systemic risk of Pakistani banks: a CoVaR approach
”,
Journal of Financial Economic Policy
, Vol. 
11
No. 
2
, pp. 
202
-
221
.
Verma
,
R.
,
Reddy
,
Y.V.
and
Ghosh
,
S.
(
2019
), “
Systemic risk in the Indian banking system: a CoVaR analysis
”,
Journal of Asian Economics
, Vol. 
63
, 101134.
Wen
,
F.
,
Xu
,
L.
,
Ouyang
,
G.
and
Kou
,
G.
(
2020
), “
Retail investor attention and systemic risk: evidence from China’s banking sector
”,
Journal of International Financial Markets, Institutions and Money
, Vol. 
65
, 101176.
Xu
,
H.
,
Li
,
Y.
and
Wang
,
S.
(
2021
), “
Time-varying systemic risk and bank heterogeneity: evidence from China using dynamic CoVaR
”,
Emerging Markets Review
, Vol. 
46
, 100746.
Published in International Trade, Politics and Development. Published by Emerald Publishing Limited. This article is published under the Creative Commons Attribution (CC BY4.0) license. Anyone may reproduce, distribute, translate and create derivative works of this article (for both commercial and non-commercial purposes), subject to full attribution to the original publication and authors. The full terms of this license may be seen at Link to the terms of the CC BY 4.0 licence.

or Create an Account

Close Modal
Close Modal