Abstract

A better understanding of financial contagion and systemically important financial institutions (SIFIs) is essential for the prevention and control of systemic financial risk. Considering the ripple effect of financial contagion, we integrate the relevant spatiotemporal information that affects financial contagion and propose to use the ripple-spreading network to simulate the dynamic process of risk contagion in China’s financial system. In addition, we introduce the smooth-transition vector autoregression (STVAR) model to identify “high” and “low” systemic risk regimes and set the relevant parameters of the ripple-spreading network on this basis. The results show that risk ripples spread much faster in high than in low systemic risk regimes. However, systemic shocks can also trigger large-scale risk contagion in the financial system even in a low systemic risk regime as the risk ripple continues. In addition, whether the financial system is in a high or low systemic risk regime, the risk ripples from a contagion source (i.e., a real estate company) spread first to the real estate sector and the banking sector. The network centrality results of the heterogeneous ripple-spreading network indicate that most securities and banks and some real estate companies have the highest systemic importance, followed by the insurance, and finally the diversified financial institutions. Our study provides a new perspective on the regulatory practice of systemic financial risk and reminds regulators to focus not only on large institutions but also on institutions with strong ripple capacity.

1. Introduction

Financial security is an important part of national security; preventing and defusing financial risk, especially holding the bottom line of no systemic financial risk, is of great significance to ensure the smooth operation of China’s economy and finance. In recent years, many financial institutions in China have been engaged actively in cross-financial business and conduct mixed financial services [1], which have strengthened the linkages between financial institutions and increased systemic risk [2]. Besides, the events of 2008 global financial crisis have reminded us that regulating financial institutions or markets in isolation and ignoring the connectedness between institutions or markets can no longer effectively prevent and control systemic financial risk [3]. Therefore, it is essential and urgent to explore the connectedness among financial institutions in China and examine the spillover and contagion effects of systemic risk by considering the financial system as a whole.

The global financial crisis in 2008 has aroused attention to the connectedness and risk contagion among financial institutions or markets and increased the emergence of related research methods. A commonly used measure of systemic financial risk is the conditional value-at-risk (CoVaR) and the delta conditional value-at-risk (ΔCoVaR) [4], which are based on a “bottom-up” perspective and reflect the risk contribution of financial institutions to the system [57]. Another commonly used measure is the marginal expected shortfall (MES) [3], which is based on an “up-bottom” perspective and focuses on the systemic risk contribution of individual financial institution when the system is in a crisis [8]. In addition, the conditional expected shortfall (CoES) [9, 10] and the SRISK risk index [11, 12], etc., have also been developed to measure systemic financial risk. Although methods such as CoVaR, MES, CoES, and SRISK are widely used, they ignore the connectedness in the financial system and take less account of the risk contagion effect. While the contagion effect of risks is at the heart of the outbreak of the systemic risk or financial crisis, when individual or local risk emerges, different institutions may be affected by risk sharing, risk spreading, and risk amplification mechanisms, resulting in chain reactions and even the collapse of the entire financial system.

Network theory, another important branch of systemic risk research, provides an effective and intuitive analytical tool for systemic financial risk by constructing network nodes and edges, which can not only display the risk contagion paths but also identify systemically important nodes [13], and thus helps the regulators to establish a stable financial system and an effective regulatory system. In general, there are two types of financial network analysis methods. The first method uses detailed information on interbank asset-liability exposures to construct financial networks and then determines the specific risk transmission mechanisms or simulates the impact of shocks on the network. Based on this method, scholars have established the interbank exposure network [14], interbank payment network [15, 16], and asset-liability network [17, 18]. Although this method allows for the identification of specific directions and paths of risk contagion, it relies on nonpublic data, which makes research limited. The second method uses financial market data, i.e., stock prices (e.g., return, volatility, and tail risk) to build financial networks and then assesses risk contagion through network topology indicators. Since financial market data are publicly available and forward-looking, they are widely used in the study of financial networks and systemic risk. Related studies mainly include correlation-based networks [1921], Granger causality networks [2224], volatility spillover networks [2527], and tail risk spillover networks [28, 29]. In addition, in order to take more risk information into account, some scholars have developed composite networks [30] and multilayer information spillover networks [3134]. In these financial networks, nodes denote financial institutions or markets, and edges denote relationships between nodes. By examining whether there are correlation, causality, or spillover effects between nodes, we can determine the risk contagion paths.

In summary, the research on systemic risk under financial network has yielded rich results. It is worth noting that, however, most of the relevant network analyses are static. Although relevant studies have used the rolling time window approach [27, 30, 34] to capture the dynamics of risk contagion, they are unable to model the dynamic ripple effects triggered by a contagion source, i.e., which financial institution is affected first and which is affected later when a financial institution suffers a loss or goes bankrupt. The ripple effect describes the gradual spreading of the effects caused by local events. In fact, the formation and evolution of many real-world complex systems depend to a large extent on the spreading of a few local events [35]. Suppose a failed (lost or bankrupt) financial institution triggers financial contagion, the contagion will spread around like ripples in a calm pond. When it reaches the nearest (highest associated) financial institution, the institution may be infected and trigger a new ripple-spreading process, which is analogous to a ripple in a pond reaching a stake and creating a new responding ripple due to the reflection effect. As the contagion continues to spread, more and more institutions may be infected. It is well known that many successful artificial intelligence (AI) techniques are actually inspired by certain natural system or phenomena [36, 37]. For example, genetic algorithms (GA) are inspired by the process of natural selection and evolution, particle swarm optimization (PSO) is inspired by the learning behaviors within populations, and artificial neural networks (ANNs) are inspired by animal brains, etc. These algorithms or their derived algorithms are also widely used in economic and financial fields [38, 39]. Following the common practice of learning from nature in the field of artificial intelligence, Hu et al. [35] developed a ripple-spreading network model (RSNM) which is inspired by the natural ripple-spreading phenomenon on the calm water surface and emphasized its application potential and flexibility in their paper. On this basis, Hu et al. [40] attempted to apply the genetic algorithm (GA) to tune the ripple-spreading related parameters and made it a great flexibility to study many real-world complex network systems.

Existing studies have shown that financial contagion has ripple effects [41]. Looking at the latest research, several scholars have applied the ripple-spreading network to model the contagion path of financial risk and identify systemically important financial institutions or markets. For example, Su et al. [42] proposed a ripple network-based collective spillover effect approach and identified the systemic importance of financial markets. Xu et al. [43] pointed out that the ripple effect is one of the most features of financial contagion and modeled the paths of China’s systemic risk contagion under different contagion sources. However, the results are far from sufficient, and there are still some issues that need to be further explored. For example, few literature studies have analyzed the ripple effect of financial contagion in the framework of regime switching. Financial markets may experience regime shifts and nonlinear risk contagion as a result of sudden structural changes due to selling behavior of common asset holders [44], investor panic [45], and asymmetric dependence on financial asset returns in both upward and downward phases of the market [46]. Feng et al. [47] pointed out that macroeconomic and financial variables often have sudden structural changes due to the exposure to external shocks such as policy changes and ignored the possible effects of regime shifts in the process of risk contagion which may lead to a significant bias. Therefore, it is necessary and urgent to explore the network connectedness and examine the dynamic ripple-spreading process of financial risk under different systemic risk regimes. In addition, relevant studies have examined financial contagion in the framework of a deterministic ripple-spreading network. However, the deterministic and uncertain factors coexist in the process of financial contagion, and how to balance this relationship deserves further exploration. In this paper, we mainly carry out the following innovative work.

First, we improve the basic algorithm of the ripple-spreading network to make it applicable to financial contagion and apply it to the analysis of financial contagion in China. The semideterministic ripple-spreading network model (SD-RSNM) can balance the deterministic and uncertain factors in the process of financial contagion. By observing the instantaneous state of the ripple-spreading network, we can identify which financial institutions are affected first and which are affected later.

Second, we introduce the smooth-transition vector autoregression (STVAR) model to identify the “high” and “low” state regimes of systemic risk and parameterize the ripple-spreading network on this basis. In terms of parameter specification, we take into account some important spatiotemporal factors affecting financial contagion, such as the magnitude of financial shock, risk amplification factor, risk resistance ability, and spreading speed of each financial institution.

Finally, we manage to identify systemically important financial institutions (SIFIs) based on the heterogeneous networks generated by the dynamic ripple-spreading processes under high and low systemic risk regimes. We find that most securities and banks and some real estate companies have the highest systemic importance in China’s financial system. In particular, the securities sector has the strongest ripple-spreading capacity and plays an intermediary role in the system network. The results remind financial regulators and government departments that systemic risk regulation should focus not only on large institutions but also on institutions with strong ripple effects.

The remainder of the article is structured as follows: Sections 2 and 3 describe the methodology and data, respectively. In Section 4, we apply our method to Chinese financial institutions and present the empirical analysis results. In Section 4, we provide a brief discussion of the findings.

2. Methodology

2.1. Basic Idea of the Ripple-Spreading Network

Considering the simultaneous existence of deterministic and uncertain factors in financial contagion, we derive the semideterministic ripple-spreading network model (SD-RSNM) based on Hu et al. [35] and Xu et al. [43] to simulate financial contagion in China’s financial system. In this model, there are two classes of parameters: contagion source related and network node related. The parameters of the contagion source include , , and . denotes the energy of the initial ripple of the contagion source, and it measures the magnitude of the systemic shock; denotes the ripple-spreading speed of the contagion source; denotes the distance between the contagion source and node . The parameters of the network nodes include , , , and (). denotes the risk amplifying factor of the node ; denotes the connection threshold, representing the risk resistance capability of the node ; is the distance between the node and . Once a node is infected, it will trigger a new ripple with energy and speed . In addition, the energy parameters, i.e., and , decay following the same function , where is a constant and is the ripple radius of the node at time . It is worth noting that, in the framework of SD-RSNM, the node follows the following activation principles: if and , then the node is activated by the node , while if and , then the node is activated with probability ; is the probability decay coefficient. Using these parameters, the specific simulation steps of SD-RSNM are described in Appendix A.

To visualize the basic principle of the risk ripple-spreading process, we give a simple example, as shown in Figure 1, where node 0 denotes the contagion source and nodes 1–3 denote normal network nodes, i.e., financial institutions. The energy of a ripple decays as it spreads; i.e., the strength of energy is reflected as the thickness of the ripple.

At , the contagion source, i.e., node 0, first triggers an initial ripple with energy and spreads out, but it has not yet reached any node.

At , the initial ripple triggered by node 0 reaches node 3. Since , node 3 is activated and a directed link from node 0 to node 3 is established. Then, node 3 generates a response ripple with energy that spreads out again.

At , the initial ripple triggered by node 0 reaches node 1. Since , node 1 is activated and generates the corresponding response ripple.

At , risk ripples triggered by contagion source disappear; i.e., the energy value decays to 0, while risk ripples triggered by nodes 1 and 3 continue to spread out.

At , the ripple triggered by node 1 reaches node 2; since , node 2 is activated by node 1 and generates the corresponding response ripple.

At , the ripple triggered by node 2 reaches node 1, since , node 1 is not activated; the ripple triggered by node 1 reaches node 0, since , node 0 is activated by node 1 and generates the corresponding response ripple. At the same time, the ripple triggered by node 1 reaches node 3, although , a dotted link is created from node 1 to 3, which is determined by chance. Then, node 3 generates a response ripple with energy .

Thus, a directed financial network consisting of four nodes and five edges is formed.

2.2. Parameter Specifications

Parameter specifications are at the heart of the ripple-spreading network, and in this paper, we set them in the context of systemic risk regime switching. The smooth-transition vector autoregression (STVAR) model proposed by Weise [48] is capable of examining asymmetric mechanisms in both high- and low-state regimes. The general form of the STVAR model can be expressed as follows:where is the m-dimensional endogenous variable in the period , is the q-dimensional exogenous variable in the period , is the lag order of the STVAR model, is the state variable, and is the transition function. Given and , the sample can be partitioned into two states, a low-state variable regime () and a high-state variable regime (); and are the k-order lag term coefficient matrices of the endogenous variables in the high- and low-state regimes, respectively. and are the exogenous variable coefficient matrices; and are the m-dimensional intercept vectors. is the transition function which takes the form of a logistic function. portrays the probability of a sample being partitioned into different state regimes; denotes the diagonal matrix; determines the degree of smoothing of the regime transition, the larger the value of , the faster the rate of regime transition; is the threshold parameter for regime partitioning.

Furthermore, (1) can be equated to the following form:where , , , and .

Following Caggiano et al. [49], the model can be simplified as follows:where , is the normalized state variable, i.e., , is the original observation, is the mean, and is the standard deviation. It should be noted that, in this paper, we use systemic risk measured by CoVaR [50] as the state variable and the risk indicator, i.e., historical volatility [26] of each financial institution as endogenous variables in the model.

On the basis of the coefficient estimation and regime partitioning of the STVAR model, following Diebold et al. [25], we can obtain the variance contribution under high and low systemic risk regimes, i.e., how much of the future uncertainty of variable is due to shocks in variable :where denotes the proportion of changes in caused by the shocks of the endogenous variable and describes the risk spillover intensity from to , is the number of endogenous variables, is the forecast period, is the standard deviation of the error term for the th equation, is the covariance matrix for the error vector , is the H-step moving average coefficient matrix, and is the selection vector, with one being the th element and zeros otherwise. Since the shocks between variables are not orthogonal, the entries of each row in the variance decomposition matrix do not add up to 1. Hence, we normalize it according to the row summation approach and obtain

Furthermore, we can set the network parameters based on the results of the regime transformation.

Connection threshold : We assume that the more external shocks an institution is exposed to, the more vulnerable it is to financial contagion. It is similar to in the connectedness method proposed by Diebold et al. [25]. Thus, based on (7), the connection threshold for institution under high and low systemic risk regimes, i.e., and , can be expressed as follows:

Amplifying factor : We use the market capitalization of each financial institution to specify under high and low systemic risk regimes. For example, the average market capitalization of ICBC under a high-risk regime is billion yuan, so we set ; the average market capitalization of ICBC under a low-risk regime is billion yuan, so we set .

Spreading speed : Following Xu et al. [43], we use the average turnover rate to specify under different systemic risk regimes. For example, the average turnover rate of ICBC under a high systemic risk regime is 0.081, so we set ; the average turnover rate of ICBC under a low systemic risk regime is 0.066, so we set .

Market distance : is determined by the inverse of the volatility correlation coefficient between institutions and ; i.e., the higher the correlation between two financial institutions, the shorter their distance and the easier it is for risk ripples to reach. Thus, the market distance between institutions and under the high- and low-risk regimes, i.e., and , can be expressed aswhere is the correlation coefficient calculated based on history volatility data between institutions and . In addition, if , then .

2.3. Heterogeneous Network

Given a lager enough , the ripple-spreading process will eventually form a stable fully connected network with nodes and edges , which is not conducive to analyzing the systemic importance of financial institutions. Following Xu et al. [43], we use the variance of node degrees to identify the heterogeneous network of the ripple-spreading process under different systemic risk regimes. For every time instant , let be the adjacency matrix of the instantaneous network in the ripple-spreading process.

The variance of node degrees at can be expressed as follows:where denotes the node degree of the node , i.e., , and denotes the average of node degrees, i.e., ; when takes the maximum value, the network is of best heterogeneity and appropriate for analyzing systemically important financial institutions (SIFIs).

3. Data

Limited to data availability, we select 55 financial institutions listed before 2011 as the main research object, specifically including 11 real estate companies, 16 banks, 14 securities, 4 insurance, and 10 diversified financial institutions, as shown in Table 1. In terms of sample selection, we include some real estate companies (collectively, financial institutions) in our research sample due to their financial-like attributes. The research sample covers the periods from January 4, 2011, to February 10, 2023. The daily historical volatility [26] of each institution can be expressed aswhere , and are the logs of daily high, low, opening, and closing prices, respectively, and the data are obtained from the Wind database.

In the STVAR model with the systemic risk index [50] as the state variable, the results of the grid point search show that the smoothing coefficient of the transition function is estimated to be 10 and the parameter takes the value 1.126. Figure 2(a) shows the results of the transition function. It can be found that the transition function shows a smooth and asymptotic trend, suggesting a nonlinear relationship of asymptotic evolution of risk contagion as the state of systemic risk changes. Furthermore, the distribution of the high and low states of systemic risk over the sample period can be obtained, as shown in Figure 2(b). The shaded areas in Figure 2(b) mark the high systemic risk regime. It can be found that the high systemic risk regime covers the period of “the market liquidity crisis caused by the money shortage in China’s banking sector in 2013,” “the stock market crash in China in 2015-2016,” and the outbreak of the public health event, i.e., “COVID-19” in early 2020. In particular, the aggregation characteristics of systemic risk are the most obvious during the period of “the stock market crash in China in 2015-2016.” Since the end of 2014, China’s stock market has seen explosive growth and high market sentiment, with investors leveraging into the market through brokerage financing and over-the-counter matching, resulting in an over-leveraged stock market and a serious market valuation bubble. However, by the second half of 2015, the market trend took a sharp turn for the worse, and the departure of leveraged funds accelerated the decline of the stock market, which eventually led to the outbreak of stock disasters and the high fluctuation of systemic risk. Compared with the actual situation, the model constructed in this paper can well identify the evolution of China’s systemic risk in various periods, indicating that the logistic smooth-switching model adopted in this paper can effectively identify various types of crisis events, and the estimation results of the model are reasonable and reliable. Finally, it can be estimated that, in the sample range from 2011 to 2023, China’s financial system has about 20% of the time in a state of high systemic risk.

On the basis of the above systemic risk regime identification, we give the behavioral parameters of the ripple-spreading network under high and low systemic risk regimes, respectively, as shown in Tables 2 and 3. The market distance, i.e., , , is measured by the inverse of the correlation coefficient of historical volatility between financial institutions. In addition, in China’s financial system, the real estate sector occupies an important position, so we set a real estate company as the contagion source to explore the risk ripple-spreading processes. To ensure sufficient network links, the energy of initial ripple of contagion source is set as .

4. Results and Discussion

4.1. Ripple-Spreading Process under High Systemic Risk Regimes

Given the model-related parameters mentioned above, we perform the risk ripple-spreading simulation processes under high systemic risk regimes according to the algorithm in Appendix A.

Figure 3 shows the dynamic ripple-spreading processes for 55 financial institutions, by giving some heat maps at t = 20, 24, 28, 32, 36, 40, 44, 48, 52, 56, 60, and 64. The rows of the heat maps represent the ripple-spreading effect from a given financial institution to other financial institutions. For example, the first row of Figure 3(a) represents the ripple-spreading effect from the real estate company, i.e., RE01, to other institutions. In addition, in order to better distinguish the ripple-spreading effects within and across sectors, we use dotted lines to partition the heat map by sectors, i.e., the real estate, banking, securities, insurance, and diversified financial sectors. The five diagonal areas of each subfigure in Figure 3 indicate the risk contagion within sectors, and nondiagonal areas indicate the risk cross-contagion between sectors. For example, the top-left corner of Figure 3(a) represents the risk contagion within the real estate sector; the bottom-left corner of Figure 3(a) represents risk cross-contagion from the diversified financial sector to the real estate sector.

From these heat maps, it can be found that the number of contagion links increases gradually over time. By observing the instantaneous state of the ripple-spreading network, we can identify which financial institutions are affected first and which are affected later. On the whole, the contagion triggered by the real estate company, i.e., RE01, first spreads to the real estate and banking sectors and then to the insurance and securities sectors and thus triggers the widespread contagion within the financial system. Specifically, Figures 3(a) and 3(b) show that the financial contagion originating from the source node, i.e., RE01, first reaches the real estate sector (RE02) and the banking sector (ABC, BOC, CCB, BCM, etc.) and then spreads to the insurance sector (CPIC, CLIC, etc.) and the securities sector (GYS). By t = 24, as shown in Figure 3(b), all banks are directly affected by the risk ripple of the contagion source. These findings suggest that real estate companies have a significant impact on the financial system, especially on the banking sector. Real estate companies are linked with banks through massive amounts of debt, so when real estate risks occur, the banking sector is affected more deeply and broadly. As risk ripples continue to spread, the cross-contagion occurs among institutions beyond the source node, and the network density increases gradually, as shown in Figures 3(c)3(l). By t = 28, as shown in Figure 3(c), the risk ripple from CCB reaches most banks such as ICBC, ABC, BOC, and BCM. Then, the risk ripple from CCB reaches the securities sector (CMS, CJS, CITIC, etc.) and the insurance sector (CLIC, CPIC, and PAIC), as shown in Figure 3(d). Meanwhile, the risk ripples from securities institutions (CITIC and GYS) begin to spread to the financial system. It is worth noting that, by t = 32, all diversified financial institutions are not yet affected by risk ripples and do not send any links to the system, suggesting that diversified financial institutions are more distant from other financial institutions.

From t = 44 to t = 64, as shown in Figures 3(g)3(l), we can find that the cross-contagion between financial institutions gradually permeates the entire financial system and the securities sector sends the most risk ripples to the system network. In addition, there are as many as 1,735 links at t = 52, which is significantly higher than that at t = 40, i.e., 603, indicating that risk ripples among financial institutions spread very fast. After the contagion channels are fully established in the early stage, the rapid risk contagion begins in the later stage. Therefore, it is wise for regulators to take timely measures to block the paths of risk contagion before the contagion channels are fully established.

4.2. Ripple-Spreading Process under Low Systemic Risk Regimes

Given the model-related parameters mentioned above, we perform the risk ripple-spreading simulation processes under low systemic risk regimes according to the algorithm in Appendix A.

For a better comparative analysis with the ripple-spreading network under high systemic risk regimes, Figure 4 similarly shows heat maps of the ripple-spreading process under the low systemic risk regime at t = 20, 24, 28, 32, 36, 40, 44, 48, 52, 56, 60, and 64. From these heat maps, we find that the contagion triggered by the contagion source, i.e., RE01, first spreads to the real estate sector and then to the banking, securities, and insurance sectors and thus triggers the cross-contagion beyond the source node. Specifically, Figures 4(a)4(c) show that the financial contagion originating from the source node, i.e., RE01, first reaches the real estate sector (RE02, RE03, RE05, etc.) and then spreads to the banking sector (ICBC, ABC, CCB, etc.), the securities sector (CJS, CITIC, GFS, etc.), and the insurance sector (CPIC, PAIC). By t = 28, as shown in Figure 4(c), most banks are directly affected by the source node, which is similar to the ripple-spreading network under high systemic risk regimes; i.e., the real estate sector has the most significant impact on the banking sector. Moreover, by t = 28, the number of links in the system network is 32, which is smaller than the links under the high systemic risk regime, suggesting that financial institutions are relatively more distant from each other under low systemic risk regimes. As the risk ripple-spreading process goes on, we find that the cross-contagion occurs beyond the source node, as shown in Figures 4(d)4(f). By t = 32, the cross-contagion occurs within the securities sector, and some securities institutions such as CITIC begin to send risk ripples or links to the banking sector. By t = 36, the risk cross-contagion occurs within the banking sector. In addition, by t = 40, we can find that most of the links constructed by the diversified financial institutions are mostly from securities institutions, while the diversified financial institutions do not issue any links to the system network. This implies that the diversified financial institutions are the main receiver of risk ripples in the Chinese financial system. From t = 44 to t = 64, as shown in Figures 4(g)4(l), we find that the risk cross-contagion relationships become more complex and the network density increases gradually. Taken together, although risk ripples spread slower in low than in high systemic risk regimes, systemic shocks can also trigger large-scale risk contagion within the financial system even in low systemic risk regime as risk ripples spread.

4.3. Network Centrality Analysis of Heterogeneous Networks

In the above analysis, we have studied the dynamic ripple-spreading processes in China’s financial system under high and low systemic risk regimes. Furthermore, we discuss the network centrality of heterogeneous ripple-spreading networks and identify the systemic importance of financial institutions (SIFIs). According to (10), each process can generate a stable heterogeneous network. We set the upper time limit to t = 500. Figure 5 shows the dynamic change trend of network heterogeneity. It can be seen that network heterogeneity increases sharply from 0 to a maximum value and then decreases gradually as risk ripples spread. The network heterogeneity in high and low systemic risk regimes reaches its maximum value at t = 49 and t = 60, respectively.

In this subsection, we consider four network centrality indicators, i.e., degree centrality (DC), closeness centrality (CC), eigenvector centrality (EC), and betweenness centrality (BC), based on heterogeneous networks selected above under different risk regimes, to identify and analyze SIFIs. These four network centrality indicators measure the systemic importance of financial institutions from different perspectives. In general, the higher the network centrality of a node, the higher its systemic importance in the system network. The implications of these four centralities are presented in Appendix A. Figures 69 show the measurement results of four network centrality indicators for each financial institution under high and low systemic risk regimes. The results are given in the form of heat maps, where darker colors indicate higher centrality values and higher systemic importance of financial institutions. Due to space constraints, only some of the English abbreviations of financial institutions are shown in heat maps, and their order from left to right is consistent with Table 1. As can be seen, whether the financial system is in a high- or low-risk regime, nodes with higher network centrality are mostly distributed in the securities, banking, and real estate sectors, and nodes with lower network centrality are mostly distributed in the diversified financial sector.

Specifically, Table 4 shows the top 15 SIFIs ranked based on four network centralities under high and low systemic risk regimes. Overall, the top 15 SIFIs are concentrated in the banking and securities sectors. However, there are some differences in the ranking of financial institutions under high and low risk regimes. In particular, the systemic importance characteristics of some securities institutions in the high systemic risk regime are particularly significant. For example, securities institutions such as GYS, HZS, IS, and GFS all appear in the top 15 lists for four centralities under high systemic risk regimes and only appear in the top 15 lists for one or two centralities under low systemic risk regimes. In addition, some real estate companies, such as RE02, RE09, and RE10, appear in the top 15 lists, while diversified financial institutions do not appear in the top 15 lists. Therefore, we conclude that most securities and banks and some real estate companies have the highest systemic importance, and the diversified financial institutions have the lowest importance. Furthermore, we average the network centrality indicators for each financial institution by sector to obtain a systemic importance ranking for each sector, and the results are shown in Table 5. It can be found that the systemic importance of the real estate sector is higher in high than in low systemic risk regimes, and supervisors should pay particular attention to the risk contagion capacity of the real estate sector when the financial system is under pressure. In addition, whether the financial system is in a high or low systemic risk regime, the securities sector has the highest “degree centrality” and “betweenness centrality,” which indicates that the securities sector not only has significant risk ripple linkage capacity but also has significant risk ripple intermediation effects. The spreading of risk ripples relies heavily on information transmission from the securities sector, thus providing a new entry point for blocking risk contagion.

5. Conclusion

In this paper, we propose using the ripple-spreading network model to reveal the spatiotemporal evolutionary characteristics of systemic risk contagion in China’s financial system. Compared to existing financial networks such as correlation networks and spillover networks, the ripple-spreading network provides a new tool for modeling the dynamic process of how financial contagion spreading from the contagion source to the whole financial system. On the one hand, the dynamic ripple-spreading process can reveal which nodes are affected first and which are affected later. On the other hand, we can identify SIFIs based on the generated heterogeneous networks.

As for the dynamic ripple-spreading processes, we find that risk ripples spread much faster in high than in low systemic risk regimes. However, systemic shocks can also trigger large-scale risk contagion within the financial system even in low systemic risk regimes as the risk ripples continue. Excessive network connectedness among institutions can amplify financial shocks through contagion effects. Overall, whether the financial system is in a high- or low-risk regime, the risk ripples from the contagion source (i.e., a real estate company) spread first to the real estate sector and the banking sector. On the one hand, institutions belonging to the same sector share similarities in terms of their business scope, business pattern, investment pattern, risk management, and financial regulation, etc., which may make it easier for the cross-contagion to occur within the sector. On the other hand, the real estate companies are linked with banks through massive amounts of debt, so when real estate risks occur, the banking sector is affected more deeply and broadly. In addition, the diversified financial institutions have fewer risk interactions in the early stage, and all the links established by the diversified institutions are mostly from securities, while it issues fewer links to the system network. This implies that the diversified institutions are the main receiver of risk ripples in the Chinese financial system.

As for the identification for SIFIs, we find that most securities and banks and some real estate companies are the most systemically important financial institutions in China’s financial system. Although we set the real estate company, i.e., RE01, as the contagion source, securities institutions exhibit the strongest risk ripple-spreading ability as the risk ripples continue, sending the most links to the system network. Especially, the systemic importance characteristics of securities institutions in high systemic risk regimes are particularly significant. To some extent, the results are not consistent with the previous studies, which noted that large financial institutions such as banks and insurances are the most SIFIs [51, 52]. For example, Wang et al. [51] noted that banks and insurance institutions in China contribute more to systemic risk than securities institutions. Chen et al. [52] pointed out that SIFIs are concentrated in the banking and insurance sectors. The reason for this is that the methods used in these studies focus mainly on assessing the systemic risk contribution of financial institutions. The higher the systemic risk contribution of a financial institution, the higher the level of systemic importance. However, the measurement of systemic risk contribution usually takes into account the impact of the size of financial institutions, so that large-scale financial institutions such as banks and insurance are identified as SIFI. Our paper mainly focuses on the network correlation of financial contagion, rather than the systemic risk contribution. Financial institutions with the highest network correlation tend to have the highest systemic importance. At the same time, this inconsistency reminds financial regulators and government departments that systemic risk regulation should focus not only on large institutions but also on institutions with strong ripple-spreading effects.

Finally, it is worth noting that the channel mechanism of financial contagion is very complex, this paper only analyzes the dynamic ripple-spreading processes of risk triggered by contagion source under high and low systemic risk regimes and does not involve the exploration of specific contagion channels or mechanisms. In the future, a more comprehensive integration of the factors affecting financial contagion will provide a better understanding of the contagion process of risks.

Appendix

A. Simulation Steps of the Ripple-Spreading Network

Given the ripple-spreading network behavior parameters, the new dynamic ripple-spreading network model can be mathematically described as follows:Step 1. Initialize the current time instant, i.e., ; initialize the current point energy of contagion source as ; initialize the current point energy of each network node as . Assume contagion source and each node have a ripple with a current radius of 0, i.e., , .Step 2. If the stopping criteria are not satisfied, do the following:Step 2.1. Let ,Step 2.2. Update the current radius and point energy of contagion source as , , where is the spreading speed of contagion source, i.e., the change in the radius of a ripple during one time instant, is a function defining how the point energy decays as the ripple spreads out. A typical decaying function can be defined as follows:where is a decaying coefficient and is the mathematical constant. Clearly, has an important influence on the decaying speed of ripples and will therefore affect the final network topology. In this paper, following Xu et al. [43], we set .Step 2.3. Check which new nodes are reached by the ripples of contagion source. Suppose is the distance between the contagion source and node . If and , then node is activated by contagion source, and thus, a link from the contagion source to node is established. Node generates a responding ripple with initial energy and . In this step, we consider the uncertainty characteristics of financial risk contagion; i.e., if and , then node generates responding ripple with the following probability:where is the probability decay coefficient. Obviously, the lower the ripple energy, the lower the probability of generating node behavior.Step 2.4. If , then update the current radius and energy of the ripple starting from node in a similar way to the ripple from the contagion source, i.e., ; Step 2.5. Check which new nodes are reached by the ripples of nodes. If and , then node is activated by node , and thus, a directed link from node to node is established. Node generates a responding ripple with initial energy and . Likewise, we consider the uncertainty characteristics of financial risk contagion; i.e., if and , then node generates responding ripple with the probability shown in equation (A.2).

Finally, we can stop the simulation in Step 2 by setting an upper time limit.

B. Networks Centrality Measures

The degree centrality of node can be expressed as follows:where is the number of nodes and is the maximum out- or in-degree. If there is a directed link from node to , ; otherwise, .

The closeness centrality for node is calculated as follows:where denotes the length of a shortest directed path from to .

Eigenvector centrality is a global measure of network centrality. It assigns a relative score to each node in the network, and in the contribution of a given node’s score, connections to nodes with high scores are larger than connections to nodes with low scores. The relative centrality score for node can be defined aswhere is a constant. It can be written as an eigenvector equation: . Typically, each eigenvector will correspond to a different eigenvalue . Only the solution corresponding to the largest eigenvalue is required by the centrality measure.

The betweenness centrality for node is calculated as follows:where , is the number of shortest directed paths linking and and is the number of shortest directed paths linking and that contain node .

Data Availability

The data used to support the findings of this study are available from the corresponding authors upon request.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This work was supported by the National Statistical Science Research Project (Grant no. 2023LY098).