Network models of financial systemic risk: A review
Fabio Caccioli
Paolo Barucca
Teruyoshi Kobayashi

November 2017
Discussion Paper No.1719

GRADUATE SCHOOL OF ECONOMICS
KOBE UNIVERSITYROKKO, KOBE, JAPANNetwork models of financial systemic risk: A review
Fabio Caccioli,1,2,3 Paolo Barucca,4,5 and Teruyoshi Kobayashi6∗
1

2

Department of Computer Science, University College London, London, UK
Systemic Risk Centre, London School of Economics and Political Science, London, UK
3
London Mathematical Laboratory, London, UK
4
Department of Banking and Finance, University of Zurich, Zurich, Switzerland
5
London Institute for Mathematical Sciences, London, UK
6
Graduate School of Economics, Kobe University, Kobe, Japan
October 31, 2017

Abstract
The global financial system can be represented as a large complex network in which banks, hedge funds
and other financial institutions are interconnected to each other through visible and invisible financial
linkages. Recently, a lot of attention has been paid to the understanding of the mechanisms that canlead to a breakdown of this network. This can happen when the existing financial links turn from beinga means of risk diversification to channels for the propagation of risk across financial institutions. In thisreview article, we summarize recent developments in the modeling of financial systemic risk. We focus
in particular on network approaches, such as models of default cascades due to bilateral exposures or
to overlapping portfolios, and we also report on recent findings on the empirical structure of interbank
networks. The current review provides a landscape of the newly arising interdisciplinary field lying at theintersection of several disciplines, such as network science, physics, engineering, economics, and ecology.

1Introduction

Since the global financial crisis of 2008–2009, many studies on financial systemic risk have been accumulated
to date. One of the most distinctive features of this newly arising field is its interdisciplinary nature, withresearchers having backgrounds in economics and finance, statistical physics, ecology, engineering, appliedmathematics, etc [1–3]. The study of financial systemic risk has attracted such a diversity of disciplines
because financial markets are complex systems, and the study of complex systems has traditionally been aninterdisciplinary field.
In the financial market there is a wide variety of market participants, such as commercial banks, insurancecompanies, hedge funds, individual investors, and central banks. These participants are interacting with each
∗ Corresponding author: kobayashi@econ.kobe-u.ac.jp1other by selling and buying financial assets, creating complex webs of financial liabilities, cross-asset holdings,
and correlations in asset returns. Financial systemic risk is, loosely speaking, the risk associated with theoccurrence of a breakdown of the financial system. Looking at systemic risk from the point of view of complex
systems means thinking of it as emerging from the interactions between different players that operate in thefinancial market. Moreover, individual participants react to the aggregate dynamics of the market that theycollectively create, so that to understand systemic risk the feedback loop between individual and collective
dynamics has to be accounted for. A key to understanding systemic risk is thus to uncover the mechanismsthat lie behind the micro-macro feedback. Because of the fact that many interactions that take place infinancial markets can be represented as a network of financial linkages between institutions, a significant
fraction of research in systemic risk has been devoted to the study of financial networks [4, 5], which is thefocus of this short review.
In this article, we provide a review of recent studies on financial systemic risk. Here, we focus on networkmodels of financial markets that has been developed outside traditional economics, although there are also
many studies of systemic risk in the literature of economics whose approaches are typically based on gametheory, finance and macroeconomic modeling. An advantage of modeling the financial system as a complex
network is that we can directly analyze complex feedback between micro- and macroscopic phenomena
without oversimplifying the structure of financial linkages. It is the structure of networks that plays anessential role in leading micro events to collective phenomena. Moreover, the empirical structure of financial
networks rarely takes a stylized form such as the Erdős-Rényi random graph and a star graph, but it takes
more complex structures such as multiplex, bipartite, core-periphery, and time-varying networks, dependingon the property of the financial linkage concerned. Such complex yet realistic network structures cannot betreated in the traditional economic models.
In section 2, we first explain basic clearing algorithms that are needed to compute the allocation of
debtor’s assets among creditors. In normal times, a bilateral credit contract is settled when a debtor repays
in full the amount of borrowed funds. However, if a debtor fails, then it is no longer straightforward to know
how to consistently allocate the debtor’s remaining assets across its creditors. It may seem reasonable toallocate the remaining assets proportionally to the amount of funds lent, but the losses that the creditors
incur may also cause their defaults. If such contagious defaults happen, the remaining assets of the failed
creditors should also be allocated proportionally to their creditors, which might in turn lead to a collapse of
the creditors’ creditors, and so on. Therefore, in the presence of contagious default cascades, calculating thefinal allocation of funds is essentially equivalent to computing a fixed point of an iterative map. We provide
a brief introduction to a widely known clearing algorithm, the Eisenberg-Noe algorithm [6].
The allocation of funds achieved by the Eisenberg-Noe’ clearing algorithm is economically sensible in thesense that the amounts of funds paid back to the creditors are endogenously determined in proportion to
the amounts of funds lent. In the real world, however, such an ideal allocation may no be feasible, especiallyin financial turmoil, because evaluating the assets of failed banks and negotiating among creditors takes along time. Therefore, many studies consider a fixed exogenous recovery rate (i.e. creditors receive a fixed
percentage of their payment), which is often set to zero to analyze worst case scenarios. An advantage ofimposing the zero-recovery assumption is that it allows us to analyze financial contagion by using a modelof social contagion that has been developed in network science. Pioneers of this approach are Gai and
Kapadia [7], who exploited the social cascade model of Watts [8]. We review some of the models of this type2

of default cascades in section 3.
Propagation of distress between borrowers and creditors can also occur before the borrower’s default,
because of credit quality deterioration. In section 4, we review a model, the so-called DebtRank [9], which
has been introduced to account for this situations, together with some of its extensions.
Stress does not only propagate between borrowers and lenders. In fact, losses can also propagate between
investors having common assets. For example, the devaluation of an asset that is commonly held by manybanks would simultaneously undermine the balance sheets of the banks holding the asset. If the loss of abank is so severe that the bank is unable to meet the minimum requirement for its capital ratio, then thebank will have to sell some of its assets. This liquidation has a negative impact on the prices of the assetsthat are being sold, which causes losses to other banks holding these assets. In this way, a cascade of defaults
could be triggered by the initial decline in the price of an asset, and fueled by the presence of overlapping
portfolios among banks. Modeling the hidden interbank linkages formed through cross asset holdings is thenimportant to understand systemic risk. Recent studies in this domain are introduced in sections 5.
In section 6, we summarize recent work on the structure of empirical interbank networks and its dynamics.
Whether the initial default of a bank can trigger a large-scale default contagion depends largely on the way
financial institutions are connected to each other. We introduce studies that examined the empirical structureof networks formed by interbank bilateral trades and explain the well-studied core-periphery structure and its
possible limitations. The necessity to analyze daily network dynamics is also discussed. Section 7 concludes.

2

Clearing algorithms

In the context of financial networks and systemic risk, clearing plays a central role as the settlementsof different transactions are entangled in a network of mutual commitments. Hence, it constitutes the
structural basis leading to observable liquidity problems, failed payments, losses, and insolvencies. In fullgenerality, financial transactions can be over-the-counter (OTC) contracts, where each pair of counterparties
has to settle its own contract, or can be collected and managed by Central-Clearing-Counterparties (CCPs)- absorbing all or some of the financial risks - or be mediated by a market via an order book dynamics.
Here, we introduce a well known clearing model, the one proposed by Eisenberg and Noe [6], which
led to a series of works on the problem of valuating systemic risk in financial networks. In a schematic
representation, a debt contract between two counterparties is represented by an amount Lab to be paid attime T by the financial institution A to financial institution B. Though this contract is simple in nature,
it is easy to realize the complexity that may arise by thinking that, in the case of conflicting financial
obligations, institution A might be unable, or even unwilling, to repay B in full at the prescribed time T .
It is therefore necessary to consider the seniority of the debt, i.e., the priority of its repayment with respect
to other obligations. Besides deterministic quantities, the counterparty B is interested in knowing ex-ante- before the maturity of the contract - the probability of default (PD) of A, i.e., the probability of A not
honoring the contract, and the loss given default (LGD), that is the credit that B will be able to recover
given the default event.
In this section, we focus on a deterministic clearing procedure [6] and in particular we will investigate
a specific case of an arbitrary number of financial institutions entangled in a mono-layer network of simple
bilateral debt contracts. All institutions in the model need to clear their payments at the time of maturity,

3

same for all bilateral contracts. Eisenberg and Noe provide a solid set of assumptions to study the basic
properties of the solution of this kind of clearing procedure. This methodology allows to compute cascades ofdefaults, reallocation of funds, and systemic effects when dealing with network of contracts. It is important
to stress that these models constitute only one mathematical aspect of the complex and multidisciplinary
problem of ensuring financial stability, involving political, legal, economic, financial, and infrastructural
aspects. Opacity and information asymmetry in markets cause a general discrepancy between credit risk
estimations by different financial institutions and the real risks. It is arguable, given the complexity ofsystem, that even a fully-fledged structural model would not be able to correctly compute credit risks. As
a consequence, opacity and financial complexity lead to a systematic inefficiency in risk taking, both in theconstruction of portfolios and in the issuing of credit [10].

2.1

The Eisenberg-Noe model

Following closely the work by Eisenberg and Noe [6], let us consider an economy composed of N financial
institutions, hereafter called banks for the sake of simplicity. Each bank has nominal liabilities to otherbanks that need to be settled at the same time. Such structure of liabilities can be represented with a N × N
matrix of non-negative real numbers L, where each entry Lij stands for the nominal liability of node i tonode j. Nominal liabilities are all non-negative because a debt contract Lij between bank i and bank j witha negative value would constitute an effective credit for bank i and therefore a debt for bank j, and would
simply appear as a positive value of the entry Lji . A further realistic assumption is the absence of nominal
claims of a bank against itself, which results in having null elements on the diagonal of the liabilities matrix.
Finally, each bank has a non-negative operating cash flow, ei , representing the net cash received by eachbank from the outside of the financial system under consideration. External liabilities can be introduced
either by setting a negative value to cash flows or by adding an extra bank in the liabilities matrix. Suchfictitious bank has zero cash flow, i.e. e0 = 0, and is supposed to receive an amount Li0 from each financial
institution i. Nevertheless, this choice is not entirely equivalent to having a negative ei , as the seniority the priority that a contract takes - of the liabilities in the matrix L, including of the ones towards node 0,
is lower than the one that would result from simply subtracting Li0 from the non-negative cash flow ei . In
the following we will take the latter choice and consider the possibility of an extra bank.
In this simplified picture, a financial system F is a pair of a non-negative liabilities matrix and a operating
cash flow vector, F = (L, e). This framework excludes the existence of many realistic characteristics offinancial systems, such as: multiple contracts between a given pair of banks, multiple levels of seniority, or
involving more than two banks, or also, different times to maturity. Further, it does not model in detail thestochastic features of cash flows, nor their correlation structure or the existence of common asset holdings
among different institutions.
Despite its specificity, this framework successfully deals with the problem of identifying a clearing vectorof payments, i.e. a vector that associates to each bank the total amount that it is able to repay to its debtors
given the financial system, (L, e), and, in doing so, it is able to introduce some crucial quantities and featuresof systemic events in financial systems, such as the dynamics of the contagion process and the conditions forthe existence of unique solutions.
Let us now enter in more detail into the definition of the clearing procedure as defined in [6]. It is useful4

to introduce a few auxiliary variables, i.e. the total nominal obligations p̄i defined as:
p̄i =

N
X

Lij ,

(1)

j=0

and the relative liabilities matrix, quantifying the fraction of liabilities from bank i that a bank j is entitled
to receive in the case of full repayment, Lij , if p̄ > 0
i
p̄i
(2)
Πij =
0 otherwise.
Finally, we define the payment vector p, the vector of unknown quantities that we want to identify, that is
the amount that each bank is actually able to repay. By definition, all the elements of the payment vector
are less than or equal to the elements of the obligation vector p̄ and greater or equal to zero, or in formulas,p ≤ p̄, p ≥ 0. The clearing procedure may yield multiple solutions for p, but, by making some intuitive andnatural financial assumptions, it can be shown that the solution is unique.
The financial requirements for the Eisenberg-Noe clearing procedure are the following: (i) all the elementsof the payment vector are less than or equal to the available cash flow of the bank (Limited Liabilities), (ii)banks repay as much as they can, i.e. they are not allowed to keep cash in their balance as long as they
have not fully repaid all their liabilities; this can also be expressed that the equity has a lower seniority
with respect to interbank liabilities; (Absolute Priority) (iii) the individual payment of a given liability, i.e.
the effective value repaid to a bank, has to be proportional to the fraction of the total obligation that theliability represents, as given by the relative liabilities matrix Πij (Proportionality). More explicitly, for each
institution the ratio between the liability repaid to a given counterparty and the total amount repaid to all
counterparties has to equate the ratio between the nominal liability and the total amount of liabilities thatthe institution has.
From these simple assumptions, we can compute the net position of a bank assuming a given payment
P
vector p. In fact, the total assets of bank i will amount to ei + j Πji pj , whilst having a total obligationin the interbank market p̄i . Given assumption (ii), as long as the obligation is less than the total amount ofP
assets then the bank will repay in full, i.e. pi = p̄i , and will remain with a net position ei + j Πji pj − p̄i .
P
Otherwise, when ei + j Πji pj < p̄i , bank i will have to use all its assets to repay its counterparties, i.e.
P
pi = ei + j Πji pj , leaving the bank with a net position equal to zero and a total amount of unrepaid debt
Pequal to p̄i − ei − j Πji pj . Hence, if we require all banks to simultaneously satisfy the same relations, wederive the following system of equations:




X
pi = min ei +
Πji pj , p̄i∀ i = 1, . . . , N,
(3)


j

which defines a fixed-point problem for the payment vector, whose components also have to satisfy theconditions p̄i ≥ pi ≥ 0 for each bank i.

2.2

Convergence to the payment vector

The system, (3) yielding the solutions to the clearing procedure, is a set of non-linear equations, for which
the dependence on each component of the payment vector is piece-wise linear, monotone, bounded, and5

continuous. All these properties allow to demonstrate the existence and uniqueness of the solution, bymeans of the general Knaster-Tarski theorem [11].
To identify the solution, Eisenberg and Noe propose a simple iterative algorithm, where each iteration ismade of two-steps: first, an update that identifies the set of defaulted banks
D(p) = {i ∈ {1, ...N } |pi < p̄i }.

(4)Secondly, the payment vector is updated by looking for the fixed-point of the following map:
p~0 = Λ(p~0 )(Π(Λ(p~0 )~p + (1 − Λ(p~0 ))p~¯ + ~e) + (1 − Λ(p~0 ))p~¯,(5)where Λ(p~0 ) is a diagonal matrix, such that Λ(p~0 )ii is one if i ∈ D(p~0 ), and zero otherwise. Equation (5)admits a unique fixed-point if the financial network is regular; the regularity condition is defined in full detailin [6], but for the sake of brevity we just notice that such condition is easily respected when each bank inthe system has a strictly positive equity.
The fixed-point p~∗ of the vector equation (5) then is used to update the set of defaulted banks D
accordingly. The two-steps are repeated until convergence to a set of defaulted banks D and payment vectorp satisfying (3), and under the regularity condition such payment vector is the only solution to (3).

2.3

Related literature and generalizations

The seminal paper by Eisenberg and Noe (EN, hereafter) ignited research and many generalizations havebeen proposed since its publication. Here we discuss only the ones that introduce the most important effects
and outline their consequences on the properties of the systemic losses.
Rogers and Veraart [12] introduce default costs in the system. In fact, when insolvency occurs it isunrealistic to assume - as EN do in their clearing model - the absence of additional costs: an insolvent
financial institution may (i) need to rapidly liquidate valuable external assets at a lower price than its
present market valuation, e.g. due to the price impact of a fire-sale, and (ii) need to withdraw their interbank
assets conceding their counterparty a discount related to the early repayment before maturity. With such
motivation the EN system is modified, introducing two factors α and β - both between zero and one -, that
discount respectively the external, and the interbank assets of an insolvent institution. The resulting systemof equations reads:p̄ , if p̄ < e + P Π p ,i
ii
j ji jpi =(6)Pαei + βΠij pj otherwise.
j

It is easy to recognize that when α = β = 1 the system of equations (6) is equivalent to (3). Otherwise,
these discount factors increase losses, introducing costs for bailouts that exceed the initial losses of thefinancial system. In fact, while EN simply redistributes losses across the financial system, as discussed in[13], Rogers and Veraart clearing procedure recognizes the existence of extra costs, which are exactly theones that financial regulatory institutions want to minimize. The motivation is that, while external losses
coming from economic shocks are driven by scarcely controllable complex economic dynamics, the possible
endogenous loss amplification due to the financial system interconnectedness could be avoided by regulators
by imposing specific policies.

6Both the EN and the Rogers-Veraart models are completely deterministic. They constitute a mechanism,
respectively, for the redistribution and amplification of losses in a financial system. Nevertheless, suchmechanism is triggered only by an actual insolvency, i.e. the liabilities have to exceed assets for a bank topropagate its losses. They are default contagion mechanisms. This is a limitation of the clearing frameworkthat could or may be overcome in two ways: first, directly account for a propagation that activates when acounterparty is in financial distress, as in DebtRank [9], and define a distress contagion mechanism; secondly,
account - before maturity - for the uncertainty on the value that the external assets will take at maturity,
i.e. introduce a probability distribution over external assets. In the latter, the propagation remains basedon a default contagion mechanism.
Nevertheless, the distribution on the external assets will include scenarios that on average will cause
expected losses, that would have been absent, if the external assets were taken at their present value before
maturity. This is the approach taken by Elsinger et al. [14]. A recent effort was made to consider the twokinds of mechanisms in a unified framework of network valuation [15] accounting for uncertainty, and default
costs in a compact form.

3

Cascades of bank defaults due to bilateral interbank exposures

3.1

The Gai-Kapadia model

Since the work of Gai and Kapadia [7], many researchers have developed network models of default cascades
in financial networks, especially interbank networks in which banks lend to and borrow from each other. The
basic structure of the Gai-Kapadia cascade model is heavily based on the Watts model of global cascades [8]that would occur on networks formed by social interactions between humans. In the Watts model, themechanism of how a node affects its neighbors is quite simple; a node gets “activated” (or “infected”) if
and only if at least a certain fraction R ∈ [0, 1] of its neighbors are activated. The Watts model of cascades
is therefore categorized as a linear threshold model or just a threshold model, which belongs to the class ofcomplex contagion.1 The main implication of the Watts model is that, on randomly connected networks,even a vanishingly small fraction of initial active nodes may lead a significant fraction of infinitely many
nodes to get activated as long as the network is not too sparse or too dense. This critical phenomenon iscalled a global cascade, and the analytic condition under which a global cascade may occur, called the cascadecondition, can be computed by exploiting a mean-field approximation.
In this section, we fist explain the basic properties of the threshold model developed by Watts [8]. Since thethreshold of activation for humans can be reinterpreted as the threshold of defaults for banks, understanding
the threshold model used in the social network literature is important to understand many existing models
of default contagion in interbank networks. Here, we explain two different approaches to calculating thesize of cascades, namely the tree-based approximation [16] and the generating function approach [7, 8].
While the Watts model assumes that edges are undirected, extending its framework to directed networks isstraightforward [7].
1 In contrast, simple contagion refers to a contagion process in which the probability of a node being affected by its neighboris exogenously given.

7

1
0.8

⇢ 0.6
0.4
0.2
0

cascade region
0

2

4

6

8

Figure 1: Threshold cascade model with a Poissonian degree distribution. Line denotes the value of ρ
calculated from Eq. (7), and circle represents the simulated average cascade size, averaged over 1000 runs
with N = 105 and R = 0.18. z is the mean degree, and seed fraction is ρ0 = 10−4 .
3.1.1

Tree-based approximation

Here, we briefly explain a tree-like approximation method for solving the model of social contagion [16].
The basic idea of a tree-like approximation is to calculate the average final fraction ρ of activated nodes byassuming that the network is locally tree-like. The probability of a randomly chosen node being active, orthe average size of global cascades, is calculated by the following equation:
!
∞
k
m
X
X
k
q m (1 − q)k−m Fρ = ρ0 + (1 − ρ0 )
pk
,
(7)
k
m
m=0
k=1

where q is the probability that a randomly chosen neighbor is active, ρ0 is the chance that a node is initially
active (i.e., a seed node), pk is the degree distribution, and k and m stand for the degree and the numberof active neighbors, respectively. At this point, the neighbors’ activation probabilities, which are considered
to be identical to the mean value, are regarded as independent since there is no cycle of influence at least
locally thanks to the assumption of a locally tree-like structure. In the simple Watts model, the response
function F (·) takes 1 if m/k > R, and 0 otherwise.
The probability q is given as
!
k−1
∞
m
X
k X k−1 m
q = ρ0 + (1 − ρ0 )pk
q (1 − q)k−1−m F
,
(8)
z m=0
k
m
k=1

where z denotes the mean degree, which is the connectivity parameter of the network. Note that one should
use the excess degree distribution kpk /z, as opposed to the degree distribution pk , to compute the averagefraction of active neighbors of a neighbor. This can be understood as a situation in which a “child node”
is influenced by its “parent nodes” and then the child node affects the “grandchildren”, and so forth. The
solution for q is then obtained as a fixed point of recursion equation (8), which in turn gives the solution for8the mean cascade size ρ from Eq. (7).2 Regardless of its simplicity, the tree-based method can predict the
final size of a global cascade very accurately (Fig. 1).
The parameter space within which a global cascade may occur is called the cascade region, and theanalytical condition for the parameters to be satisfied in the cascade region is called the cascade condition.
To see the derivation of the cascade condition, let us define S(q) as
!
∞
k−1
m
X
k X k−1 m
S(q) ≡.
(9)pk
q (1 − q)k−1−m F
z m=0
k
m
k=1

Gleeson and Cahalane [16] argue that if the derivative of the RHS of Eq. (8) (i.e., ρ0 + (1 − ρ0 )S(q)) nearq = 0 takes a value larger than one, then a vanishingly small initial seed ρ0 results in a large value of ρ. The
first-order cascade condition is therefore given by
  

∞
X
1
k(k − 1)(1 − ρ0 )pk F− F (0) > 1.
(10)
z
k
k=1

In fact, a comparison with simulation results reveals that this is not a very accurate condition for theparameter space (R, z). Therefore, they also propose the second-order cascade condition given by(C1 − 1)2 − 4C0 C2 + 2ρ0 (C1 − C12 − 2C2 + 4C0 C2 ) < 0,P∞
where it is assumed that ρ20 ≈ 0, and C` is defined such that S(q) = `=0 C` q ` and
!
!
∞ X
`
n
X
k−1
`
k
(−1)`−n pk F
C` ≡
.
z
k
`
n
n=0

(11)

(12)

k=`+1

Condition (11) states that the second-order approximation of Eq. (8), q = ρ0 + (1 − ρ0 )(C0 + C1 q + C2 q 2 ),has no solution near q = 0, because the existence of a positive root near q = 0 implies that a global cascade
is impossible.3 Gleeson and Cahalane [16] showed that the second-order cascade condition (11) well matches
the cascade region predicted by the numerical simulation.
3.1.2

Generating function approach

Now we explain the generating function approach to calculating the expected cascade size. In doing so, we
compute the probability that a randomly chosen node is vulnerable; we say a node is vulnerable if its degree,k, satisfies R ≤ 1/k. That is, if a node is vulnerable, then the node will get activated if at least one neighboris active. Let µk = P [R ≤ 1/k] denote the probability of a node having k edges being vulnerable. Since
the probability that a randomly chosen node has degree k is pk , the generating function of vulnerable nodedegree is given as
G0 (x) =∞Xµk pk xk .

(13)

k=0
2 When there are multiple fixed points, the smallest solution is selected as the valid solution. In fact, depending on the model
parameters, recursion equation (8) may exhibit a saddle-node bifurcation, which could cause a phase transition [16]. Such aphase transition can be observed in a variety of threshold models [17, 18].
3 To obtain the second-order condition, we can simply approximate S(q) near q = 0 up to the second order: S(q) =S(0) + S 0 (0)q + 21 S 00 (0)q 2 , where C0 = S(0), C1 = S 0 (0), and C2 = 21 S 00 (0).

9Generating function G0 (x) has information on all of the moments of the degree distribution only of vulnerable
nodes.4 The generating function for the excess degree distribution for the vulnerable nodes leads toP∞kµ p xk−1
G0 (x)
P k k
= 0
,
(14)
G1 (x) = k=1z
k=1 kpk
where G0 represents derivative. Note that G1 (1) is equal to the probability that a randomly chosen neighboris vulnerable.
We now introduce the generating function for the vulnerable cluster size:
H0 (x) =H1 (x) =

∞
X
n=0
∞
X

θn xn ,

(15)

θ̃n xn ,

(16)

n=0

where θn denotes the probability that a randomly chosen node belongs to a vulnerable cluster of size n, andθ̃n is the corresponding probability for a neighbor of a randomly chosen node. The generating function forthe probability that a randomly selected neighbor belongs to a vulnerable cluster should satisfy the following
self-consistency equation [8, 20]:
H1 (x) = 1 − G1 (1) + xG1 (H1 (x)).

(17)

The first term represents the probability that a neighbor is not vulnerable, and the second term corresponds
to the size distribution of vulnerable clusters to which a chosen neighbor belongs. Note that if a node
belongs to a vulnerable cluster of size n, then its neighbors must also belong to a vulnerable cluster ofsize n. Therefore, the generating function for the probability that a randomly chosen neighbor belongs to avulnerable cluster of size n depends on the second- and higher-order neighbor’s generating function, resulting
in a self-consistent determination (17) [8, 20]. Here, the assumption of a locally tree-like structure is needed
to obtain the generating function, in which case different neighbors belong to (locally) independent subsets
of a vulnerable cluster. Once H1 (x) is obtained, H0 (x) is computed as
H0 (x) = 1 − G0 (1) + xG0 (H1 (x)),(18)

We note that, roughly speaking, the procedure for calculating H1 (H0 ) corresponds to the derivation of fixed
point q (ρ) in the recursion equation (8) (Eq. (7)) in the Gleeson-Cahalane’s [16] tree-based method.5
The average vulnerable cluster size is given byhni = H00 (1)= G0 (1) +(G00 (1))2.
z − G000 (1)It follows that the cascade condition is expressed asXz < G000 (1) =k(k − 1)µk pk .

(19)

(20)

k=1
4 See Ch. 13 of Newman [19] for the basic explanation of generating functions.
5 This correspondence between the generating function approach and the tree-based method is not rigorous in the sense thatthe size of a vulnerable cluster is not identical to the average cascade size [8].

10Interbank asset

Interbankliability

AIB

LIB

External asset

DepositD

AE

CapitalK
Figure 2: Stylized balance sheet of a bank.
Note that since µk = F (1/k) for k > 0, Eq. (20) is equivalent to the previous cascade condition derived fromthe Gleeson-Cahalane’s method (Eq. (10)) as long as the threshold value is strictly positive (i.e., R > 0 and
F (0) = 0) and the seed fraction is sufficiently small (i.e., ρ0 → 0).
3.1.3

Financial contagion

In models of financial contagion, nodes and directed edges represent banks and lending-borrowing relationships, respectively, and an activation of a node is interpreted as a bank default. In fact, the Gai-Kapadia
model is isomorphic to the Watts model, the difference being that the former treats a directed random graphwhile the latter focuses on an undirected random graph. To see this, let us consider a stylized balance sheetof a bank (Fig. 2). Suppose that each bank may have two types of assets: interbank assets, AIB , and external
assets, AE (such as stocks, bonds, etc). On the liability side, there can be interbank liabilities, LIB , anddeposits from customers, D. Then, the solvency condition for bank i is given byE
IB
AIB
i + Ai − Li − Di > 0,(21)which is equivalent to saying that the net worth (or the capital) of a bank should be positive.
Now, consider the situation in which the amounts of loans extended from a bank to other banks areevenly distributed, so that each interbank exposure is simply written as AIB
i /|Ni |, where Ni denotes the setof borrowers to which bank i lends. We assume that the ratio of total interbank assets, AIBi , to net worth,
∀i
for
R
>
0.
Under
these
assumptions,theKi , is common across banks and is given by AIB
/K
=
1/R
i
i
default condition for bank i leads to
φi >

Ki
= R,
AIB
i

(22)

where φi is the fraction of bank i’s counterparties that have defaulted. The loss given default is assumed tobe 100% for simplicity; the lender would loose the full amount of funds lent to a defaulted bank. Eq. (22)states that bank i will default if the actual fraction of defaulted counterparties exceeds a constant threshold11

level R, which is essentially the same mechanism as that of social contagion in the Watts model (i.e., justreplacing R with R). Note that if bank i holds a sufficient amount of capital buffer such that R > 1, thenthere is no chance for bank i to (contagious) defalt simply because the capital buffer can fully absorb the
maximum possible losses. Since edges exposed to the risk of bank default are out-going edges (i.e., lending
to other banks), the mechanism of social contagion applies here in a straightforward manner as long as there
are no bidirectional edges.6 The trick here is that the volumes of interbank exposures, or edge weights, areevenly distributed between borrowers to which a bank lends. This greatly simplifies the analysis because
otherwise we have to replace φi in the default condition (22) with the total fraction of losses that bank i
incurs. A more general case of heterogeneous edge weights will be discussed in the following section.
Given the isomorphic property, it is straightforward to analyze cascades of bank defaults in the same
way with the methods developed by Gleeson and Cahalane [16] and Watts [8]. Due to its simplicity andanalytical tractability, over the past years the Gai-Kapadia model has been frequently used as a baseline
framework of more sophisticated models of financial contagion. We review some of these extensions in thenext section.

3.2

Extensions of the threshold cascade model

The Gai-Kapadia model spurred a flurry of research on financial contagion due to interbank exposures.
Before reviewing these works, let us summarize some of the important assumptions made in the simplest
version of the Gai-Kapadia model; i) loans are evenly distributed, ii) interbank lending forms a Erdős-Rényi
random graph, and iii) the risk of external assets is not considered. Recently, various models are proposed tomake the cascade model more realistic by relaxing these assumptions. In considering the possible extensionsof the threshold financial cascades, we can take advantage of the isomorphic property that characterizes theWatts and Gai-Kapadia models. That is, any models proposed as extensions of the Watts model could beapplied to the model of financial contagion as well.
3.2.1

Heterogeneous edge weights

When a heterogeneity of loan weights is introduced in the Gai-Kapadia model, the default condition is nolonger captured by Eq. (22). The condition (22) is valid only if the amounts of interbank loans that a bankhas lent to other banks are identical. Otherwise, the default condition cannot be expressed just by thefraction of defaulted borrowers, but rather expressed by the ratio of losses to total interbank assets:P
IBj∈Nidef Aij
> R,
(23)
AIB
i
P IB
IB
def
where AIB
is theij denotes the amount of funds lent from bank i to bank j (i.e., Ai =j Aij ), and NiIB
set of bank i’s borrowers that have defaulted. Note that assuming identical edge weights, Aij = AIB
/|N
i |,
i
def
recovers condition (22) since the LHS would reduce to |Ni |/|Ni | ≡ φi .
Under a generic default condition (23), the standard mean-field approximation will not be appropriate
since different borrowers have different weights, meaning that the number of defaulted banks itself is not
6 The presence of bidirectional edges, although they rarely appear in empirical financial networks [21], could affect the sizeof contagion [22, 23].

12

informative. One obvious way to analyze such a more general environment is to rely on numerical simulations.
Hurd and Gleeson [24], Hurd [25] and Unicomb et al. [26], however, proposed alternative approximation
methods to compute the solution of cascade dynamics. Hurd and Gleeson [24, 25] consider a situationin which edge weights w are random variables, whose CDFs Gkk0 (w) depend on the degrees of nodes inboth sides of an edge k and k 0 . They show that by imposing additional assumptions, one can obtain theanalytical solutions for the mean cascade size. Unicomb et al. [26] extend Gleeson’s [27, 28] approximate
master equation, showing that an increased weight heterogeneity will reduce the size of cascades.
3.2.2

Non-Erdős-Rényi networks

The Erdős-Rényi random graph [29] is probably the most widely used network structure in the analyticalmodels of global cascades, yet no empirical financial networks exhibit the Erdős-Rényi structure. For example, it has been shown that the degree distribution of interbank networks follows fat-tail distributions
such as a power-law distribution and a log-normal distribution [30–32]. Moreover, although the most analytical methods assume a locally tree-like structure, empirical networks have local clusters, and edges havea degree-degree correlation called assortativity [19].
It is natural to think that the tree-based method could be extended to the configuration graphs with
arbitrary degree distribution pk as long as the network has a locally tree-like structure. A possible problem isthat the tree-like assumption might no longer hold true once a more realistic network structure is considered.
Even in the configuration model, for example, the clustering coefficient can be large when a fat-tail degree
distribution is assumed [19].7 Fortunately, recent studies develop several ways in which the presence of
local cycles would not affect the accuracy of the analytical solutions. Melnik et al. [33] provide someconditions under which the tree-like approximation works “unreasonably” well even in networks with highlevels of clustering. Radicci and Castellano [34] develop an alternative technique based on a message-passing
algorithm, which gives an accurate approximation in networks with local clusters. Ikeda et al. [35] also showthat the presence of local clusters will enhance the chance of global cascades.
Another possible departure from the Erdős-Rényi graph is that there are negative degree-degree correlations (or disassortativity) in real-world financial networks [36, 37]. That is, banks with high degrees arelikely to trade with low-degree banks. Dodds and Payne [38], Payne et al. [23, 39] and Hurd et al. [40] study
the effect of (dis)assortativity on the level of systemic risk, allowing for an arbitrary degree distribution.
They show that assortativity of financial linkages strongly affects the expected cascade size.
We can extend the tree-based method to even a multiplex structure of interbank networks. A multiplex
network is a networked system consisting of multiple layers, on each of which nodes are connected by edges.8
Interbank networks may exhibit a multiplex structure when banks trade different types of assets. For
example, if there is a seniority in interbank assets (i.e., difference in the risk of loans), a “monoplex” model
would no longer suffice. Brummitt and Kobayashi [42] generalized the Gai-Kapadia model in a way thatallows for different seniority levels, where different risk assets are traded in different layers. They consider
a general case in which there are M seniority levels, showing that the cascade condition is generally given
7 In the configuration model, the clustering coefficient will approach 0 as the network size goes to infinity if the second

moment of the degree distribution takes a finite value. See Eq. (13.47) of [19] for details.
8 The essential difference between multiplex networks and multilayer networks is that in the latter, different layers can havedifferent nodes while each layer has the same set of nodes in the former [41].

13

by the trace of the Jacobian of M many recursion equations. Of course, variety of interbank assets is not
limited to seniority levels. There can be other “layers” in which banks trade long- and short-term assets,
foreign exchange exposures, derivatives, and etc. [43, 44]. The multiplex structure of financial networks,
however, is a relatively pre-matured research area in the sense that analytical models of financial contagion
are still scarce.
3.2.3

Risk of external assets

The only contagion channel considered in the Gai-Kapadia model is a cascade of repayment failures in theinterbank market. In reality, however, this channel is just a part of the source of systemic risk. The riskof devaluation in external assets is one of the major concerns not only for the portfolio management of
individual banks, but also for the systemic risk of the entire financial system. In recent years lots of works
have been done on the contagion channel through overlapping portfolios, in which a fall in the price ofan asset will simultaneously affect many banks holding the same (or a correlated) asset [45–48]. Such asimultaneous shock to multiple banks has a potential to accelerate the traditional contagion process throughinterbank exposures. We will explain these studies in detail in sections 4 and 5.
In the Gai-Kapadia model, it is considered that the risk of external assets plays a role in initiating acontagion. Suppose that returns of external assets held by different banks are independent, and the priceof an asset held by bank i falls. If the devaluation of the asset is so large that the default condition (22) issatisfied (due to a reduction in Ki ), then it may cause bank i’s creditors to default, initiating a contagion
process. However, the role of external assets in reality is not that simple because the actual external assets
are correlated, and the volatility of assets makes the health of balance sheets differ from bank to bank.
To take into account these more realistic situations, many studies conduct simulations to understand theimpact that a correlation in external assets has on systemic risk [49]. Kobayashi [50] provides a simple way togeneralize response function F to include the possibility that the value of external assets follow a probability
distribution, which in fact corresponds to the Watts model in which the threshold for contagion is a randomvariable [16].

4Distress propagation due to credit quality deterioration

The importance of counterparty default contagion for practical purposes has been challenged both theoretically and empirically. From the theoretical point of view Glasserman et al. [51] proved for instance that,within the Eisenberg-Noe framework, the contribution of contagion to the default probability of a bank isalways small, and Battiston et al. [52] showed that this is the case because of a “conservation of losses” thatis implicitly embedded in the Eisenberg-Noe algorithm, which prevents it from amplifying exogenous shocks.
From the empirical point of view, contagion analysis of real interbank systems have shown domino effects
triggered by the failure of a small number of banks are unlikely to occur in practice [53]. On the other hand
it was shown that networks of interbank exposures can significantly amplify distress propagation in presenceof other contagion channels, such as for instance fire sales and overlapping portfolios [48]Beyond its interaction with other contagion mechanisms, another reason why networks of interbank
exposures can be important is the following: Models of contagion due to counterparty default risk assume14

that losses propagate from borrowers to lenders only after the default of a borrower. However, in practice,losses could occur even in absence of default, because of credit quality deterioration [54]. Consider thesituation in which bank i is exposed to bank j, which suffers a large loss. After the loss, the probability thatj defaults has increased, and therefore the expected cash flow associated with the exposures between i and jis reduced. If interbank assets were to be marked to market, this would mean that the value of the interbank
asset of i that is associated with its exposure to j is reduced. The idea of accounting for the propagation of
distress before defaults led to the introduction of DebtRank [9].

4.1

DebtRank

Let us consider a system of N banks, and let us denote by Wij the interbank exposure of i towards bankj, by Aext
the external (non interbank) assets of bank i, and by Li its total liabilities. DebtRank is aidiscrete-time map that describes the evolution of the equity of all banks after a shock hits the system. In
the DebtRank dynamic banks can be in two states: active or inactive. An active bank is a bank that willpass distress to its creditors if subject to a loss, and a bank becomes inactive after it has passed distress to
its creditors once. This does not necessarily mean that the bank has defaulted, nor that the bank cannot
suffer additional losses, but simply that further losses will not be transmitted to its creditors. If we denotei (t)by hi (t) = Ei (0)−E
the relative loss of equity of bank i at time t, and by A(t) the set of active banks at
Ei (0)
time t, the DebtRank dynamic reads




X Wij
hi (t + 1) = min 1, hi (t) +hj (t)(24)
Ei (0)j∈A(t)

A(t + 1)={i | hi (t) > 0, hi (t − 1) = 0}(25)

The meaning of the above dynamic is the following: the loss experience by bank i between time 0 and timet + 1 is its loss up to time t plus the new losses that are transmitted by its active counterparties. The
contribution to the loss of i due to counterparty j is proportional to the level of distress of j (the factor
hj (t)) and to the exposures of i towards j relative to its equity (the factor Wij /Ei (t)). The matrix withelements Wij /Ei (t) has been named matrix of interbank leverage ([55]) because it represents the percentage
loss of equity of i that corresponds to a 1% devaluation of its exposure to j.
In the original DebtRank paper, Battiston et al. [9] present a study of US commercial banks, and theyshow that their algorithm can effectively be used to rank banks in terms of their systemic importance.
By showing that relatively small banks can be among the most systemically important, and because of theanalogy between DebtRank and centrality measures in networks, they introduced into the debate on systemicrisk the idea that some banks might be “too central to fail”.

4.2

Extensions

According to the above formulation, because nodes become inactive after they propagate distress once, losses
can flow through a cycle in the network only once. To account for further rounds of propagation, Bardoscia

15

et al. [56] derived, from the iteration of the balance sheet identity, the following modified dynamic


N


X
Wij
hi (t + 1) = min 1, hi (1) +hj (t) ,
Ei (0)

(26)

j=1

(27)
where hi (1) is the initial exogenous shock that affects bank i, and it is assumed that hi (0) = 0 for alli = 1, . . . N . This formulation makes it easier to understand the stability of a system with respect to a small
perturbation. In particular, if the largest eigenvalue of the matrix of interbank leverages is larger than one,shocks will be amplified by the network and lead to the default of some banks in the system.
The underlying assumption of DebtRank is that losses propagate from borrowers to lenders linearly: an
x% devaluation of the equity of the borrower leads to an x% devaluation of the interbank asset of the lender.
This assumption can be relaxed by considering dynamics of the form


N


X
Wij
hi (t + 1) = min 1, hi (1) +f (hj (t)) ,(28)
Ei (0)
j=1

with f (x) a function that maps the interval [0, 1] into the positive real semiaxis. For instance, Bardoscia etal. [57] considered the following function
f (x) = xe−α(x−1) ,(29)

where α ≥ 0. This function represents a one-parameter family of distress propagation rules that interpolates
between the threshold model used in [7], which is recovered in the limit α → ∞, and the linear rule of
DebtRank, which corresponds to α = 0. By performing contagion analyses on a system of European banks,
they explored the dependence of the model on the parameter α, showing the existence of different regimes
for what concerns the amplification of distress.
Bardoscia et al. [58] also considered the case of a non-linear propagation of distress, and considered f (x)to be increasing and convex. They analyzed the largest eigenvalue of the interbank leverage matrix, and they
showed the existence of trajectories in the space of networks that can turn a system from stable to unstablethrough processes that are normally believed to increase the stability of financial markets, namely market
integration and diversification. To show this they considered the hypothetical case in which the network ofinterbank contract is a directed acyclic graph, which is stable. They then considered a situation in whichlinks are randomly added between banks in the system, but in such a way that, every time a link is added,
the total amount of lending of each bank is preserved. This implies that banks are on average increasing
their diversification. They showed that, through this process of increasing diversification, it is possible for an
initially stable network to become unstable. They argue that the instability of a network under these typeof dynamics is due to the emergence of peculiar cyclical structures in the network of interbank exposures.
DebtRank is also at the basis of the stress testing framework proposed by Battiston et al. [55], who
propose a framework based on the following steps: 1) application of an exogenous shock to the systemand estimation of direct losses; 2) propagation of distress through the DebtRank dynamic Eq. (25) andestimation of second-round losses; 3) further (third-round) losses caused by banks liquidating a common
asset to target their initial leverage. Through the application of this stress testing framework, Battiston et16

al. [55] find that the second-round effects (due to DebtRank) and the third-round effects (due to leveragetargeting) dominate the first-round losses (direct losses due to the exogenous shock). This finding has
potential implications for regulators, as it implies that stress tests that do not account for network effects
can significantly underestimate systemic risk.
In relation to policy making, an interesting work has been presented by Thurner and Polenda [59] andPolenda and Thurner [60], who, using DebtRank as a tool to measure systemic risk, showed how taxation
policies that do take into account the impact of interbank contracts on systemic risk can effectively promote
systemic stability while not reducing the volume of interbank lending.

5Overlapping portfolios and price mediated contagion

Section 3 discussed contagion due to counterparty default risk. Here we discuss a different contagion mechanism, which is associated with the fact that stress propagates between investors that hold common assets.
The idea of this contagion mechanism is the following: Consider the simple situation in which two banks iand j invest in a common asset x. Suppose now that bank i is under stress, and that in order to reduce its
risk exposure it has to liquidate part of its position on asset x. Because of market impact, the tendency of
prices to react to trading activity, the liquidation procedure causes a devaluation of the asset, whose pricewill drop. When assets are marked to market, this devaluation causes a loss to bank j. We then see that,even in the absence of direct contracts between i and j (such as those associated with interbank loans considered previously), stress can propagate from i to j through the intermediation of the price of the common
asset x. Similarly to the case of counterparty default risk, the question to be asked is therefore “how does
the pattern of overlapping portfolios between banks, which can be modeled as a bipartite network, affectsystemic risk“? A pictorial representation of a simple network of overlapping portfolios is shown in Fig. 3.
The effect of losses due to common asset holdings and fire sales has first be studied by Cifuentes etal. [61] in the context of the Eisenberg-Noe model. In their paper, Cifuentes et al. [61] consider a systemof banks that are interacting through a network of interbank lending relationships and in which all banks
are investing in one common external asset. Banks are subject to a capital constraint, and thus need toliquidate part of their investment in the common asset if they face a loss. The authors present numericalsimulations on a system of 10 banks to show the response of the system to the initial default of a bank, andthey study the effect of changing the average connectivity of the interbank network. They find that there isa non-monotonic relationship between the number of connections in the network and the number of observed
defaults. A similar setting is also explored by Gai and Kapadia [7], where the interbank lending network is
however modeled as a directed Erdős-Rényi network and counterparty default contagion is modeled throughthe threshold dynamics discussed in section 3, as well as in Nier et al. [62]. May and Arinaminpathy [63]build on the model of Nier et al. [62], of which they present a mean-field solution, by considering that banks
interact through different asset classes and accounting for contagion between those asset classes.
The papers mentioned above considered the effect of fire sales of one or few asset classes, but their focus
was the study of the stability of a system as a function of the properties of the network of interbank loans.
More recently, the focus shifted towards the study of the network of overlapping portfolios itself and on how
its shape affects systemic stability. The network of overlapping portfolios is usually modeled as a bipartite
network, where two types of nodes exist (banks and assets) and a link can only connect a bank to an asset,

17

BanksAssets

1

A
2

B

3

4

C

5

Figure 3: Pictorial representation of a network of overlapping portfolios. A bank is connected to the assets in
its balance sheet. Stress can propagate between banks with common assets. For instance: if bank 1 is under
distress and liquidates its portfolio, asset A would be devalued. This would cause a loss to bank 2, whichmight in turn need to liquidate its investment portfolio. This liquidation would cause assets A to be further
devalued and asset B to be devalued as well, thus causing losses to banks 3, 4 and 5 and the consequent
devaluation of asset C.
meaning that the bank is investing in that asset. If we consider a system of N banks and M assets, we candescribe the structure of the system in terms of the matrix Q, where the element Qia is the number of sharesof asset a held by bank i, and we also denote by pa the price of asset a. Apart from the network, there are
two main ingredients that are needed to define a model. The first is the response of the bank to its losses,
the second is the response of the asset to its liquidation. If we consider a dynamic that occurs over discrete
time steps t = 1, 2, . . ., the response of bank i can be defined in terms of the map
Qia (t) = fia [Qia (t − 1), Ai (t − 1), Ei (t − 1)] ,(30)

where we have denoted by Ai (t) the value of the assets of bank i at time t and by Ei (t) its equity, while theresponse of asset a can be expressed in terms of the map
pa (t) = ga [{Qia (t)}] ,(31)

where we denote by {Qia (t)} the set {Q1a (t), Q2a (t), . . . , QN a (t)}.
In the literature on network models of overlapping portfolios, two choices are common for what concern
the response of banks: Either banks are passive until they default, at which point they liquidate their entire
portfolio, or they target a certain level of leverage, defined as the ratio between the mark to market value of
their assets and their equity: In the following we briefly discuss some of these models. Although with some18

differences, the models we discuss are all very similar in their ingredients, but the analysis quite different intheir focus.

5.1

Threshold dynamics

Huang et al. [46] consider the situation in which a bank is passive until its default, and it liquidates its entire
portfolio when it defaults, so that
Q (0), if E (t − 1) ≥ 0
iai
fia [Qia (t − 1), Ai (t − 1), Ei (t − 1)] =.
(32)
0, if Ei (t − 1) < 0
In Huang et al. [46] asset prices are assumed to respond to liquidation asPia (0) − Qia (t)]i [QPga [{Qia (t)}] = pa (0) 1 − α,i Qia (0)(33)

where α ≥ 0 is a parameter related to the market impact associated with asset a. The above expression
means that the value of the asset at time t depends linearly on the fraction of its shares (relative to the total
number of shares held in the system) that has been liquidated up to that time.
Huang et al. [46] perform an empirical analysis concerning the situation of US commercial banks in 2007.
They consider data for 7,846 commercial banks and 13 asset classes, and they perform stress tests by reducing
the value of one asset class from pa (0) to (1 − ξ)pa (0), with 0 ≤ ξ ≤ 1. They then compute the number
banks that survive the cascading process triggered by the initial shock. They find that abrupt transitions
occur in the number of surviving banks as a function of the parameters α and ξ and that devaluation ofcommercial real estate loans are responsible for the failure of commercial banks during the subprime crisis.
Quite interestingly, Huang et al. [46] also perform an empirical validation of their model by comparing thebanks that their model predicts should fail with those banks that actually failed between 2008 and 2011.
Their analysis of false positive and true positive rates shows that there is predictive power in the model.
A similar model of overlapping portfolios is the one of Caccioli et al. [47], who also consider a bipartite
network of banks and assets and the map (32) for the update of banks positions on the assets. The rule forthe devaluation of assets is however linear in the log returns9 , as the one of [7, 61]. This can be written asP[Qia (0)−Qia (t)]α i Pi Qia (0)ga [{Qia (t)}] = pa (0) 1 − e.
(34)
Caccioli et al. [47] study the stability of the system in the limit when the number of banks and assets islarge. In particular they identify the conditions under which a small initial perturbation such as the initialbankruptcy of a bank or the devaluation of an asset can lead to a global cascade of bankruptcies. They showthat the model can be described in terms of a branching process. In particular, using the approximationthat the network is a tree, they define a transfer matrix Π, whose element Πij represents the probability
that the only failure of bank j triggers the failure of bank i:
"M
#


X
P
Πij = prob
Qia pa (0) 1 − e−αQja / k Qka > Ei .
(35)
a=1
9 Caccioli et al. [47] also test their results with respect to other choices of the market impact function, such as the linear

function (33) and a square root market impact law.

19

25
20
15
10
5
5

10

15

Figure 4: Region of instability for the cascade model of overlapping portfolios for a bipartite Erdős-Rényi
network. Within the red region the system displays global cascades.
The stability of the system as a function of the model parameters can at this point be assessed by studying
the largest eigenvalue of Π. In the paper, Caccioli et al. [47] provide results for networks in the bipartite
Erdős-Rényi ensemble, and, similarly to the case of counterparty default risk [7], they find the existence of a
non-monotonic relation between average diversification and probability of observing a global cascades. They
also show the existence of a critical value of leverage below which, independently on network connectivity,
the system is always stable with respect to the initial shock. This is shown in Fig. 4, which presents anillustration of the unstable region as a function of the average diversification of banks (i.e. the average
degree of banks in the network of overlapping portfolios) and the leverage of banks, which is defined as themarket value of the bank’s investment portfolio divided by its equity. The figure refers to the same setting
considered in [47], with a bipartite Erdős-Rényi networks, and under the assumption that all banks have thesame leverage. The effect of heterogeneous degree distributions on this model is studied in Banwo et al. [64]by means of numerical simulations.

5.2

Leverage targeting

Caccioli et al. [47] also consider the effect of relaxing the assumption that banks are passive investors, andlook at what happens if banks decide to target their initial leverage over the dynamics. If a bank suffers aloss, its leverage will go up, therefore a rebalancing is needed to reduce risk. They consider the situation inwhich banks react to losses at time t by liquidating a fraction
λi Ei (t)∆Ai (t) = γAi (t) 1 −(36)
Ai (t)of their investment. In the above formula, λi is the target leverage of bank i, while the parameter γ ∈ [0, 1]
determines how quickly the bank tries to reach its target. ∆Ai (t) is the total value of the assets liquidated by20

the bank at time t. To know how many shares of a given asset a are sold, one needs to divide by the number
ki of different assets in the bank’s portfolio (it is assumed here that banks liquidate the same fraction ofeach asset), and by the current price of asset a. This leads to the following response functionQia (t − 1) − γ Ai (t−1) 1 − λi Ei (t−1) , if Ei (t − 1) ≥ 0
ki pa (t−1)
Ai (t−1)
fia [Qia (t − 1), Ai (t − 1), Ei (t − 1)] =. (37)
0, if E (t − 1) < 0iThey find that the attempt of banks to reduce their individual risk through preemptive liquidation ends up
significantly widening the region of parameter space where global cascades can be observed.
A stress test framework based on shock propagation due to leverage targeting is the one proposed by
Greenwood et al. [65]. In this case banks estimated their losses at time t and then reduce the size of theirinvestment in order to restore their initial leverage at time t + 1, which corresponds to the map given inequation (37) with γ = 1. Deleveraging causes prices to be devalued, and consequently mark-to-market losses
for banks, that will further need to deleverage, and so on. The market impact function considered by [65] isa linear function like the one in equation (32). They try to disentangle the effect of different factors on thelosses produced by fire sale contagion, showing that the contribution of a bank to aggregate deleveraging is
the higher if a bank is more connected, bigger, more leveraged and more exposed to the initial shock. They
also introduce the concept of indirect vulnerability of a bank with respect to a given asset, as measured by
the loss of the bank’s equity due to the deleveraging associated with the asset. This has to be contrasted
with the notion of direct vulnerability, which is the direct exposure of the bank towards the asset. They test
the model on the largest 90 banks in the European Union for the period 2009-2011. Through a regression
analysis on those banks that are publicly traded, they find that direct and indirect vulnerabilities have thesame explanatory power on banks’ returns.
The model of Greenwood et al. [65] is at the basis of the framework developed in Duarte and Eisenbach [66]to measure aggregate vulnerability and systemic importance of banks. Greenwood et al. [65] also show that
their measure of aggregate vulnerability due to spillover effects grows well before the crisis of 2008, and they
are able to disentangle the contribution to aggregate vulnerability due to the increase of leverage, systemsize and concentration of investments in illiquid assets.
Cont and Schaanning [67] consider a dynamic which is in between that of a passive investor and a leveragetargeting one. In their model, they account for the fact that there is usually a buffer between the leverageof a bank and the maximum leverage allowed by regulation. This is done so that they are not forced toliquidate their position because of a relatively small loss. They consider then the bank as a passive investor
until its loss makes the bank break its leverage constraint. When this happens, the bank develerages inorder to reach its target, which is however below the maximum allowed by the regulatory regime (so thata small buffer is restored). They also distinguish between marketable securities, that can be liquidated andare subject to market impact, and illiquid assets that are not marketable and therefore cannot be liquidated.
Deleveraging only involves marketable securities.
Cont and Schaanning [67] consider data collected by the European Banking Authority on 51 European
banks, and they compare the outcome of stress tests performing with their model vs. target leveraging. In
terms of price changes, they consider both a linear and square root market impact. They find that there are
significant differences in the losses estimated with the two models. Quite interestingly, they also introduce
a matrix of overlaps between the portfolios of different banks, where each asset is weighted by its liquidity,

21

and they show that, although many financial institutions have zero overlap between their portfolios, they are
all connected by second order overlaps. This means that stress tests that do not account for second round
losses can significantly underestimate systemic risk.
Although there may be randomness in the shocks that hit the system, or in the construction of the
network of overlapping portfolios, the dynamics described above are all deterministic. A stochastic dynamic
is for instance the one considered in Corsi et al. [68]. They consider a system of N banks investing in m
randomly chosen assets out of a universe of M possible assets, where m is however computed as the optimalvalue of diversification that corresponds to banks maximizing their profit conditional on a VaR constraint,
which is equivalent to leverage targeting [69]. The value of the assets is then the sum of a linear market impact term that depends on banks trading plus a stochastic component, which is in turn the sum of a common
factor and an idiosyncratic component. As for the other models we discuss here, the relative compositionof the portfolio does not change over time, but at each time banks change the volume of their investment
on the portfolio in order to maintain their target leverage. Corsi et al. [68] show that, upon increasing
diversification, the system goes from a stable regime where time series of asset returns are stationary to an
unstable regime where they are characterized by bubbles and bursts.

6

Empirical structure of interbank networks

There are many works that aim to characterize the structure of real-world financial networks. Measuring
and analyzing the structure of financial networks has a two-fold objective: On one hand, knowledge of thestructure of financial networks gives insights on how local risks would spread over the entire network throughfinancial linkages. This class of studies aims to measure the systemic risks of particular network topologies
that could emerge and remain unchanged for a certain period of time (i.e., static structure). On the other
hand, since the topology changes over time, knowledge about the dynamical transition patterns could allow
us to predict how systemic risk would evolve over time. In this section, we provide a brief review of these
two lines of research that study the static and dynamic structures of interbank networks.

6.1

Static structure

6.1.1

Interbank networks in different countries

Over the past decade, the topology of interbank networks has been examined in many countries. These studies
include Boss et al. [30] for Austria, Upper and Worms [70] for Germany, Degryse et al. [71] for Belgium,
Van Lelyveld and Liedorp [72] for Netherland, Iori et al. [31] and Bargigli et al. [43] for Italy, Wells [73] andLangfield et al. [74] for the UK, Furfine [75] for the US, Cont et al. [32] for Brazil, Martinez-Jaramillo et
al. [76] for Mexico, and Imakubo and Soejima [77] for Japan.
While in some countries bilateral transactions data are available,10 in many countries the aggregatebalance-sheet data (e.g., total amount of loans) are the only source of information for bilateral trades. In
such cases, one needs to estimate the interbank network structure using a suitable estimation method. A
10 One example is the data for the Italian interbank market, in which banks trade through an online platform called e-MID.

The e-MID data are commercially available from e-MID S.p.A based in Milan (http://www.e-mid.it/).

22

widely used one is the maximum entropy (ME) method. The ME method estimates the network structure by
maximizing the entropy of interbank linkages, which implies that the total interbank lending is distributed to
all the possible borrowers as evenly as possible. A disadvantage of the ME method is that it is likely that theestimated network is much denser than the actual one. Mistrulli [78] argues that the ME method may over- or
underestimate the risk of default contagion. To overcome the problem, recently more sophisticated methods
are also proposed [79]. Anand et al. [80] compare the accuracy of several existing estimation methods byapplying them to various empirical networks.
A more direct way to extract information of bilateral transactions is to use interbank payment data.
Since payment flows contain information on interbank settlements and transfers, one could filter out the
information of bilateral interbank loans. This approach is taken by, among others, Furfine [75], Demiralp etal. [81], and Imakubo and Soejima [77].
6.1.2

Core-periphery structure

It has been argued that the structure of interbank networks at certain point in time is best described as acore-periphery structure [82]. A core-periphery structure is formed by two groups: core and periphery. The
core and peripheral nodes are distinguished as follows; the core forms a subgraph of the entire network inwhich nodes are connected densely to each other. Peripheral nodes are connected to the core nodes but not
to other peripheral nodes. This is expressed by a block adjacency matrix as
"
# "
#
CC CP
1
CP
A=
≈
,
(38)
PC PP
PC
0
where CC denotes the submatrix representing the connectivity among core nodes, and PC represents theconnectivity between the core and peripheral nodes. PC is identical to CP since we consider undirected
graphs. In general, detecting core-periphery structure in a directed graph is a challenging problem, andmost of the existing methods for core-periphery detection are developed for undirected graphs. In the purecore-periphery structure, we should have CC = 1 (i.e., a complete graph) and PP = 0 (i.e., there is nolink between peripheral nodes). Classifying all the nodes into core and periphery is a nontrivial task, and a
popular way to do this is to find core nodes so as to minimize the difference between the empirical adjacency
matrix and the ideal core-periphery block matrix (38) [82–84].
6.1.3

Empirical networks: core-periphery vs. bipartite structure

We summarize the works estimating the core-periphery structure in interbank networks in Table 1. Theseempirical studies differ in terms of the data and their time scales. For example, Fricke and Lux [84] andBarucca and Lillo [85, 86] used the data for interbank transactions, but the former studied quarterly aggregate
networks and the latter analyzed daily networks. Imakubo and Soejima [77] studied the transactions datafiltered out of interbank payment data, and other studies are based on the regulatory data reported byfinancial institutions to the financial authorities.
While many empirical works on core-periphery structure are essentially based on the standard detection
method described above, a more flexible and widely used approach to detecting block structure, called thestochastic block model (SBM) [87], has also been used. The SBM is a probabilistic model of random graphs

23June 23, 2014
— July 4, 2014

All ties
April 23, 2001—May 7, 2001

June 12, 2007—June 25, 2007

Figure 5: Visualization of Italian interbank networks, e-MID. The networks are aggregated over 10 businessdays. Red and black circle denote Italian and other foreign banks, respectively. Visualization is done bypython-igraph with the Kamada-Kawai algorithm [91].
with a flexible block structure; nodes are assigned to different blocks and each pair of nodes is linked witha probability depending on the nodes’ blocks. This model can generate arbitrary block structures, such as
core-periphery, modular, and bipartite structures. In fact, Barucca and Lillo [85, 86] employ this approachand find that the two-block structure that best represents the e-MID overnight money market is bipartite(i.e. borrowers and lenders) at the daily resolution.
Even in the case of aggregate networks, the plausibility of the stylized core-periphery structure as acharacteristic of interbank networks can be controversial. As is evident from Eq. (38), the previous works
implicitly assumed that an empirical network consists of a single core block and a single peripheral block.
This suggests that even if there were no such a standard core-periphery structure in the empirical network,
the estimation method classifies each node as either a core node or a peripheral node. In fact, recent studies
argue that the seemingly core-periphery structure might just come from a heterogeneous degree distribution
or consist of two cores [88–90]. As an example, visualization of Italian interbank networks aggregated over 10
business days is presented in Fig. 5. It appears that there are densely connected core nodes at the center ofthe network while there are also (seemingly) peripheral nodes that are not linked to each other yet connectedto the core. However, in 2007, Italian banks and other foreign banks seem to form two core-like groups,
which would make it difficult to extract a stylized pure core-periphery structure. It might be reasonable toinfer that there are two core-periphery structures in the network [89].

ignificant ties

6.2

Interbank network dynamics at the daily scale

6.2.1

Why daily scale?

Although the majority of empirical works are based on static and aggregate networks, the granularity of thedata can be crucial for recognizing the heterogeneous behavior of financial institutions and for capturing the24Reference
Fricke and Lux [84]Craig and von Peter [83]Veld and van Lelyveld [92]Langfield et al. [74]
Imakubo and Soejima [77]Barucca and Lillo [85, 86]

Data typetransactions data
regulatory data
regulatory data
regulatory data
settlement datatransactions data

Country
Italy
Germany
Netherlands
UK
Japan
Italy

Resolution
quarterly
quarterly
quarterly
annualmonthly
daily

Year
2015
2014
2014
2014
2010
2016

Table 1: Empirical studies on the core-periphery structure in interbank networks.
functioning of a market at its inherent time-scales. In interbank markets, most of the bilateral transactions
are overnight, meaning that the relationship between two banks as a lender and a borrower lasts only for
one day or shorter, depending on the time when the loan contract is made. In the Italian interbank market,
for example, more than 86% of transactions are overnight lending in the period between 2000–2015 [21].
If the issue of interest is to understand the interconnectedness of financial risks, then an aggregate
network of overnight lending relationships contains irrelevant information because different edges formed
in different days in fact do not exist at the same time. Many researchers studied aggregate networks notnecessarily because it conveys information about the interconnected risk structure, but rather because it
would reveal meaningful information about the structure of long-term relationships among banks. Anotherpossible reason for why aggregate networks attract attention is that daily networks can be much sparser
and noisier than aggregate networks, and they appear to change their structure from day to day in a purely
random manner [93, 94].
6.2.2

Daily network dynamics

Barucca and Lillo [85, 86] show that at the daily resolution, the structure of interbank networks in the e-MID
interbank market is not characterized by a core-periphery structure, but rather characterized by a bipartite
structure or more general community structures. Lowering the time resolution (such as weeks or months)
tends to increase the likelihood that a core-periphery structure is detected, yet a non-negligible fraction ofsuch networks are still best characterized by a bipartite structure. Kobayashi and Takaguchi [21] reinforce
their result by showing that the bipartivity of daily interbank networks have been increasing over the past
decade.11
If we regard interbank markets as dynamic systems in which the structure of bilateral exposures changes
every day, an interesting question to ask is whether the daily dynamics are just random or there are robust
and time-invariant properties. Kobayashi and Takaguchi [21] find that the daily market activity represented
by combination (N, M ), where N and M are the numbers of active banks and edges, respectively, is strictly√
ruled by a superlinear relationship N ∝ M 1.5 , or hN i ∝ M , independently of the structure and the size ofdaily networks (Fig. 6a). They also find several daily dynamical patterns in the e-MID market, such as thepower-law distribution of transaction duration (Fig. 6b) and a tent-shaped distribution of weight growth.
11 Bipartivity is a measure of bipartite structure proposed by Estrada and Rodrı́guez-Velázquez [95], which takes 0.5 if the

network is complete and one if purely bipartite.

25

(a)10(b) 100

3

data
slope = 2.9

CCDF

10-1

2

10

-2

10

10-3
10-4
10-5
100

1

10

101

102

101

102

duration (days)

103

Figure 6: Scaling laws in the Italian overnight interbank networks. (a) Superlinear relationship between thenumbers of nodes (N ) and edges (M ) between September 4, 2000 and December 31, 2015 (3922 businessdays). Each dot corresponds to a day. (b) Complementary cumulative distribution function (CCDF) of
transaction duration (in terms of the number of business days) of bank pairs between 2010–2015.
Interestingly, these properties are ubiquitous in social networks formed via human interactions such as phone
calls and face-to-face interactions [96–98].

7

Discussion

In this article, we reviewed recent works studying financial systemic risk based on network approaches. While
we tried to cover as many research topics as possible, these are obviously not exhaustive. In particular, twoimportant research areas that have not been discussed in the current article are the prediction and the control
of systemic risk.
As in other fields of sciences, there are generally three steps for the study of systemic risk to mature asa scientific research field: The first step is to understand and model the mechanisms behind the real-world
phenomena. The second is to forecast the future state of the system. The third and last step is to control
the system to avoid the occurrence of undesired phenomena. Most of the studies discussed in the current
article are still at the first stage. Researchers from various fields have just started working together since thelate 2000’s to develop models, such as the models of interbank default cascades through bilateral exposuresand overlapping portfolios, that can be used to describe the real-world phenomena analytically.
In recent years, however, a growing number of researchers are tackling the controllability of systemic risk
by simulating possible policy tools that could be taken by the financial regulators. We still need further
studies at every step to deepen our understanding of the complexity of financial networks and reduce systemicrisk. We hope that this review article will encourage researchers from various areas of sciences to join in thischallenging research field.

26

Acknowledgements
F.C. acknowledges support of the Economic and Social Research Council (ESRC) in funding the Systemic Risk Centre (ES/K002309/1). P.B. acknowledges financial support from FET Project DOLFINS nr.
640772. T.K. acknowledges financial support from the Japan Society for the Promotion of Science Grants
no. 15H05729 and 16K03551.

References[1] May R. M., Levin S. A., Sugihara G. (2008) Complex systems:
451(7181), 893–895.

Ecology for bankers. Nature[2] Schweitzer F., Fagiolo G., Sornette D., Vega-Redondo F., Vespignani A., White D. R. (2009) Economic
networks: The new challenges. Science 325(5939), 422–425.
[3] Battiston S., Farmer J. D., Flache A., Garlaschelli D., Haldane A. G., Heesterbeek H., Hommes
C., Jaeger C., May R., Scheffer M. (2016) Complexity theory and financial regulation. Science351(6275), 818–819.
[4] Haldane A. G., May R. M. (2011) Systemic risk in banking ecosystems. Nature 469(7330), 351–355.
[5] Fouque J.-P., Langsam J. A. (2013) Handbook on Systemic Risk. Cambridge University Press.
[6] Eisenberg L., Noe T. H. (2001) Systemic risk in financial systems. Management Science 47(2), 236–249.
[7] Gai P., Kapadia S. (2010) Contagion in financial networks. Proceedings of the Royal Society A 466(2120),
2401–2423.
[8] Watts D. J. (2002) A simple model of global cascades on random networks. Proceedings of the National
Academy of Sciences USA 99(9), 5766–5771.
[9] Battiston S., Puliga M., Kaushik R., Tasca P., Caldarelli G. (2012) Debtrank: Too central to fail?
financial networks, the fed and systemic risk. Scientific Reports 2, 541.
[10] Battiston S., Caldarelli G., May R. M., Roukny T., Stiglitz J. E. (2016) The price of complexity infinancial networks. Proceedings of the National Academy of Sciences USA 113(36), 10031–10036.
[11] Tarski A. (1955) A lattice-theoretical fixpoint theorem and its applications. Pacific Journal of Mathematics 5(2), 285–309.
[12] Rogers L. C., Veraart L. A. (2013) Failure and rescue in an interbank network. Management Science
59(4), 882–898.
[13] Visentin G., Battiston S., D’Errico M. (2016) Rethinking financial contagion. arXiv:1608.07831.
[14] Elsinger H., Lehar A., Summer M. (2006) Risk assessment for banking systems. Management Science
52(9), 1301–1314.

27

[15] Barucca P., Bardoscia M., Caccioli F., D’Errico M., Visentin G., Battiston S., Caldarelli G. (2016)
Network valuation in financial systems. arXiv:1606.05164.
[16] Gleeson J. P., Cahalane D. J. (2007) Seed size strongly affects cascades on random networks. PhysicalReview E 75(5), 056103.
[17] Lee K.-M., Brummitt C. D., Goh K.-I. (2014) Threshold cascades with response heterogeneity in multiplex networks. Physical Review E 90(6) 062816.
[18] Kobayashi T. (2015) Trend-driven information cascades on random networks. Physical Review E
92(6), 062823.
[19] Newman M. E. J. (2010) Networks: An Introduction. Oxford University Press.
[20] Callaway D. S., Newman M. E., Strogatz S. H., Watts D. J. (2000) Network robustness and fragility:
Percolation on random graphs. Physical Review Letters 85(25), 5468.
[21] Kobayashi T., Takaguchi T. (2017) Social dynamics of financial networks. arXiv:1703.10832.
[22] Boguñá M., Serrano M. A. (2005) Generalized percolation in random directed networks. Physical Review
E 72(1), 016106.
[23] Payne J., Harris K., Dodds P. (2011) Exact solutions for social and biological contagion models on
mixed directed and undirected, degree-correlated random networks. Physical Review E 84(1), 016110.
[24] Hurd T. R., Gleeson J. P. (2013) On Watts’ cascade model with random link weights. Journal of Complex
Networks 1(1), 25–43.
[25] Hurd T. R. (2016) Contagion!: Systemic Risk in Financial Networks. Springer.
[26] Unicomb S., Iñiguez G., Karsai M. (2017) Threshold driven contagion on weighted networks.
arXiv:1707.02185.
[27] Gleeson J. P. (2011) High-accuracy approximation of binary-state dynamics on networks. Physical Review Letters 107(6), 068701.
[28] Gleeson J. P. (2013) Binary-state dynamics on complex networks: Pair approximation and beyond.
Physical Review X 3(2), 021004.
[29] Erdős P., Rényi A. (1959) On random graphs I. Publicationes Mathematicae 6, 290–297.
[30] Boss M., Elsinger H., Summer M., Thurner S. (2004) Network topology of the interbank market. Quantitative Finance 4(6), 677–684.
[31] Iori G., De Masi G., Precup O. V., Gabbi G., Caldarelli G. (2008) A network analysis of the Italianovernight money market. Journal of Economic Dynamics and Control 32(1), 259–278.
[32] Cont R., Moussa A., Santos E. B. (2013) Network structure and systemic risk in banking systems. In:
Fouque J.-P., Langsam J. A. (eds) Handbook on Systemic Risk, Cambridge University Press, New York.

28

[33] Melnik S., Hackett A., Porter M. A., Mucha P. J., Gleeson J. P. (2011) The unreasonable effectivenessof tree-based theory for networks with clustering. Physical Review E 83(3), 036112.
[34] Radicchi F., Castellano C. (2016) Beyond the locally treelike approximation for percolation on real
networks. Physical Review E 93(3), 030302.
[35] Ikeda Y., Hasegawa T., Nemoto K. (2010) Cascade dynamics on clustered network. Journal of Physics:
Conference Series 221(1), 012005.
[36] Soramäki K., Bech M. L., Arnold J., Glass R. J., Beyeler W. E. (2007) The topology of interbank
payment flows. Physica A 379(1), 317–333.
[37] Bech M. L., Atalay E. (2010) The topology of the federal funds market. Physica A 389(22), 5223–5246.
[38] Dodds P. S., Payne J. L. (2009) Analysis of a threshold model of social contagion on degree-correlated
networks. Physical Review E 79(6), 066115.
[39] Payne J., Dodds P., Eppstein M. (2009) Information cascades on degree-correlated random networks.
Physical Review E 80(2), 026125.
[40] Hurd T. R., Gleeson J. P., Melnik S. (2017) A framework for analyzing contagion in assortative banking
networks. PLOS ONE 12(2), 1–20.
[41] Kivelä M., Arenas A., Barthelemy M., Gleeson J. P., Moreno Y., Porter M. A. (2014) Multilayer
networks. Journal of Complex Networks 2(3), 203–271.
[42] Brummitt C. D., Kobayashi T. (2015) Cascades in multiplex financial networks with debts of different
seniority. Physical Review E 91(6), 062813.
[43] Bargigli L., Di Iasio G., Infante L., Lillo F., Pierobon F. (2015) The multiplex structure of interbank
networks. Quantitative Finance 15(4), 673–691.
[44] Poledna S., Molina-Borboa J. L., Martnez-Jaramillo S., van der Leij M., Thurner S. (2015) The multilayer network nature of systemic risk and its implications for the costs of financial crises. Journal of
Financial Stability 20, 70–81.
[45] Beale N., Rand D. G., Battey H., Croxson K., May R. M., Nowak M. A. (2011) Individual versus
systemic risk and the regulator’s dilemma. Proceedings of the National Academy of Sciences USA
108(31), 12647–12652.
[46] Huang X., Vodenska I., Havlin S., Stanley H. E. (2013) Cascading failures in bi-partite graphs: Modelfor systemic risk propagation. Scientific Reports 3, 1219.
[47] Caccioli F., Shrestha M., Moore C., Farmer J. D. (2014) Stability analysis of financial contagion due tooverlapping portfolios. Journal of Banking & Finance 46, 233–245.
[48] Caccioli F., Farmer J. D., Foti N., Rockmore D. (2015) Overlapping portfolios, contagion, and financial
stability. Journal of Economic Dynamics and Control 51, 50–63.

29[49] Kobayashi T. (2013) Network versus portfolio structure in financial systems. European Physical Journal
B 86(10), 434.
[50] Kobayashi T. (2014) A model of financial contagion with variable asset returns may be replaced with asimple threshold model of cascades. Economics Letters 124, 113–116.
[51] Glasserman P., Young H. P. (2015) How likely is contagion in financial networks? Journal of Banking& Finance 50, 383–399.
[52] Battiston S., D’Errico M., Visentin G. (2016) Rethinking financial contagion. arXiv:1608.07831.
[53] Upper C. (2011) Simulation methods to assess the danger of contagion in interbank markets. Journalof Financial Stability 7(3), 111–125.
[54] Glasserman P., Young H. P. (2015) Financial networks. Department of Economics Discussion Paper
753, University of Oxford.
[55] Battiston S., Caldarelli G., D’Errico M., Gurciullo S. (2016) Leveraging the network: A stress-test
framework based on DebtRank. Statistics & Risk Modeling 33(3-4), 117–138.
[56] Bardoscia M., Battiston S., Caccioli F., Caldarelli G. (2015) DebtRank: A microscopic foundation forshock propagation. PLOS ONE 10(6), e0130406.
[57] Bardoscia M., Caccioli F., Perotti J. I., Vivaldo G., Caldarelli G. (2016) Distress propagation in complex
networks: the case of non-linear DebtRank. PLOS ONE 11(10), e0163825.
[58] Bardoscia M., Battiston S., Caccioli F., Caldarelli G. (2017) Pathways towards instability in financial
networks. Nature Communications 8, 14416.
[59] Thurner S., Poledna S. (2013) Debtrank-transparency: Controlling systemic risk in financial networks.
Scientific Reports 3, 1888.
[60] Poledna S., Thurner S. (2016) Elimination of systemic risk in financial networks by means of a systemicrisk transaction tax. Quantitative Finance 16(10), 1599–1613.
[61] Cifuentes R., Ferrucci G., Shin H. S. (2005) Liquidity risk and contagion. Journal of the European
Economic Association 3(2-3), 556–566.
[62] Nier E., Yang J., Yorulmazer T., Alentorn A. (2007) Network models and financial stability. Journal of
Economic Dynamics and Control 31(6), 2033–2060.
[63] May R. M., Arinaminpathy N. (2010) Systemic risk: The dynamics of model banking systems. Journalof the Royal Society Interface 7(46), 823–838.
[64] Banwo O., Caccioli F., Harrald P., Medda F. (2016) The effect of heterogeneity on financial contagion
due to overlapping portfolios. Advances in Complex Systems 19(08), 1650016.
[65] Greenwood R., Landier A., Thesmar D. (2015) Vulnerable banks. Journal of Financial Economics
115(3), 471–485.

30

[66] Duarte F., Eisenbach T. M. (2015) Fire-sale spillovers and systemic risk. Staff Report no.645, Federal
Reserve Bank of New York.
[67] Cont R., Schaanning E. F. (2017) Fire sales, indirect contagion and systemic stress testing. Working
Paper 2/2017, Norges Bank.
[68] Corsi F., Marmi S., Lillo F. (2016) When micro prudence increases macro risk: The destabilizing effectsof financial innovation, leverage, and diversification. Operations Research 64(5), 1073–1088.
[69] Shin H. S. (2010) Risk and Liquidity, Oxford University Press, Oxford.
[70] Upper C., Worms A. (2004) Estimating bilateral exposures in the German interbank market: Is therea danger of contagion? European Economic Review 48(4), 827–849.
[71] Degryse H., Nguyen G., et al (2007) Interbank exposures: An empirical examination of contagion riskin the Belgian banking system. International Journal of Central Banking 3(2), 123–171.
[72] van Lelyveld I., Liedorp F. (2006) Interbank contagion in the Dutch banking sector: A sensitivity
analysis. International Journal of Central Banking 2(2), 99–133.
[73] Wells S. (2004) Financial interlinkages in the United Kingdom’s interbank market and the risk ofcontagion. Bank of England Quarterly Bulletin 44(3), 331.
[74] Langfield S., Liu Z., Ota T. (2014) Mapping the UK interbank system. Journal of Banking & Finance
45, 288–303.
[75] Furfine C. (2003) Interbank exposures: Quantifying the risk of contagion. Journal of Money, Credit,and Banking 35(1), 111–128.
[76] Martinez-Jaramillo S., Alexandrova-Kabadjova B., Bravo-Benitez B., Solórzano-Margain J. P. (2014)
An empirical study of the Mexican banking systems network and its implications for systemic risk.
Journal of Economic Dynamics and Control 40, 242–265.
[77] Imakubo K., Soejima Y., et al (2010) The transaction network in Japan’s interbank money markets.
Monetary and Economic Studies 28, 107–150.
[78] Mistrulli P. E. (2011) Assessing financial contagion in the interbank market: Maximum entropy versus
observed interbank lending patterns. Journal of Banking & Finance 35(5), 1114–1127.
[79] Mastrandrea R., Squartini T., Fagiolo G., Garlaschelli D. (2014) Enhanced reconstruction of weighted
networks from strengths and degrees. New Journal of Physics 16(4), 043022.
[80] Anand K., van Lelyveld I., dm Banai, Friedrich S., Garratt R., Haaj G., Fique J., Hansen I., Jaramillo
S. M., Lee H., Molina-Borboa J. L., Nobili S., Rajan S., Salakhova D., Silva T. C., Silvestri L., de Souza
S. R. S. (2017) The missing links: A global study on uncovering financial network structures from partial
data. Journal of Financial Stability, in press.
[81] Demiralp S., Preslopsky B., Whitesell W. (2006) Overnight interbank loan markets. Journal of Economics and Business 58(1), 67–83.

31

[82] Borgatti S. P., Everett M. G. (2000) Models of core/periphery structures. Social Networks 21(4), 375–
395.
[83] Craig B., von Peter G. (2014) Interbank tiering and money center banks. Journal of Financial Intermediation 23(3), 322–347.
[84] Fricke D., Lux T. (2015) Core–periphery structure in the overnight money market: evidence from thee-mid trading platform. Computational Economics 45(3), 359–395.
[85] Barucca P., Lillo F. (2016) Disentangling bipartite and core-periphery structure in financial networks.
Chaos, Solitons & Fractals 88, 244–253.
[86] Barucca P., Lillo F. (2017) The organization of the interbank network and how ECB unconventional
measures affected the e-MID overnight market. Computational Management Science, in press.
[87] Peixoto T. P. (2017) Bayesian stochastic blockmodeling. arXiv:1705.10225.
[88] Rombach M. P., Porter M. A., Fowler J. H., Mucha P. J. (2014) Core-periphery structure in networks.
SIAM Journal on Applied mathematics 74(1), 167–190.
[89] Kojaku S., Masuda N. (2017) Finding multiple core-periphery pairs in networks. arXiv:1702.06903.
[90] Kojaku S., Masuda N. (2017) Core-periphery structure requires something else in the network.
arXiv:1710.07076.
[91] Kamada T., Kawai S. (1989) An algorithm for drawing general undirected graphs. Information Processing Letters 31(1), 7–15.
[92] Veld D. in ’t, van Lelyveld I. (2014) Finding the core: Network structure in interbank markets. Journalof Banking & Finance 49, 27–40.
[93] Finger K., Fricke D., Lux T. (2013) Network analysis of the e-MID overnight money market: the informational value of different aggregation levels for intrinsic dynamic processes. Computational Management
Science 10(2-3), 187–211.
[94] Musmeci N., Battiston S., Caldarelli G., Puliga M., Gabrielli A. (2013) Bootstrapping topologicalproperties and systemic risk of complex networks using the fitness model. Journal of Statistical Physics
151(3-4), 720–734.
[95] Estrada E., Rodrı́guez-Velázquez J. A. (2005) Spectral measures of bipartivity in complex networks.
Physical Review E 72(4), 046105.
[96] Cattuto C., Van den Broeck W., Barrat A., Colizza V., Pinton J.-F., Vespignani A. (2010) Dynamics
of person-to-person interactions from distributed RFID sensor networks. PLOS ONE 5(7), 1–9.
[97] Starnini M., Baronchelli A., Pastor-Satorras R. (2013) Modeling human dynamics of face-to-face interaction networks. Physical Review Letters 110(16), 168701.

32

[98] Schläpfer M., Bettencourt L. M., Grauwin S., Raschke M., Claxton R., Smoreda Z., West G. B., Ratti
C. (2014) The scaling of human interactions with city size. Journal of the Royal Society Interface
11(98), 20130789.

33