# Improvement on Load-Induced Cascading Failure in Asymmetrical Interdependent Networks: Modeling and Analysis

Authors: Haiyan Han; Rennong Yang
Journal: Mathematical Problems in Engineering
 (2015)

---

## Abstract

Many real-world systems can be depicted as interdependent networks and they usually show an obvious property of asymmetry. Furthermore, node or edge failure can trigger load redistribution which leads to a cascade of failure in the whole network. In order to deeply investigate the load-induced cascading failure, firstly, an asymmetrical model of interdependent network consisting of a hierarchical weighted network and a WS small-world network is constructed. Secondly, an improved “load-capacity” model is applied for node failure and edge failure, respectively, followed by a series of simulations of cascading failure over networks in both interdependent and isolated statuses. The simulation results prove that the robustness in isolated network changes more promptly than that in the interdependent one. Network robustness is positively related to “capacity,” but negatively related to “load.” The hierarchical weight structure in the subnetwork leads to a “plateau” phenomenon in the progress of cascading failure.


---

## Body

## 1. Introduction

Complex network theory has been a successful tool in modeling and analysis of modern systems [1]. A variety of network models have been proposed to approximate realistic systems such as power grids [2], transportation systems [3], communication networks [4], and other systems. Most of them focus on single and isolated systems. Recently,interdependent network [5] provides a new insight in understanding the structure [6], percolation [7], spreading processes [8], evolution games [9], and robustness [10, 11] of complex systems. This kind of structure usually consists of two or more subnetworks, with the subnetworks working dependently on each other [12, 13]. For example, power grids and computer systems must depend on each other because computers require power grids to supply electricity. Power grids, in turn, rely on computer systems to control power transmission process.In 2010, Buldyrev et al. [14] found that interdependent networks have become significantly more vulnerable than their noninteracting counterparts under random attack. Gao et al. [5] reviewed the connectivity properties of “networks of networks” formed by interdependent random networks. Hu et al. [15] constructed a partially coupled network with both interdependent and interconnecting links. The interconnections satisfy “one to one” condition. They found that the change of interconnecting links leads to the change of the phase transition from second order to first order through hybrid phase transition. Huang et al. [16] studied the robustness of interdependent networks under targeted attack on high or low degree nodes. It provided a routine method to study the degree-based targeted attack problems in both single networks with dependency links [17, 18] and other general randomly connected and uncorrelated interdependent networks. Moreover, Parshani et al. [19] described the dynamic process of cascading failures on two partially interdependent networks.Currently, researchers focus on the problem of modelling an interdependent network model for analyzing the progress of cascading failure. Based on classical network models such as ER [20], WS [21], and BA [22], researchers studied well some symmetrical and asymmetrical networks [19, 23, 24], such as ER-ER, WS-WS, BA-BA, ER-WS, BA-ER, and BA-WS; but these models differ greatly from real-world systems [25, 26]. Taking interdependent networks like ground transportation network and airline network as an example, ground transportation network displays a hierarchical property. That is, stations at provincial level have higher capacity and more importance than municipal stations. However, the airline network can be seen as a single level network with small world property [27]. So, some other researchers [28–30] made their trials to construct framework of asymmetrical interdependent networks which are more authentic to reality.Moreover, load such as cargoes transported in the transportation network and electric stream in power grids can trigger cascading failures [31]. The load carried by the failed nodes or edges will not disappear but will flow to the remaining part of the network, which, in a possible way, will cause further failures. As far as we are concerned, not enough attention has been paid to the cascading failure induced by load redistribution in interdependent networks.In order to investigate the load-induced cascading failure in interdependent networks, we propose asymmetrical interdependent networks model in Section2 and the load-induced cascading failure model in Section 3. Section 4 simulates the cascading failure in the proposed model and compares it with a WS-WS symmetrical network when nodes and edges suffer from intentional attacks, respectively. Section 5 summarizes the contribution of this paper and identifies future research needs.


## 2. Network Model

Figure1 shows an asymmetrical interdependent network model. N
e
t
w
o
r
k
A has a hierarchical and weighted structure, where nodes are assigned with weights and distributed into different levels, where N
e
t
w
o
r
k
B is a single level network. The numbers of nodes in N
e
t
w
o
r
k
A and N
e
t
w
o
r
k
B are equal. The coupling proportion is set to p
1. Nodes in different subnetworks randomly construct “one to one” connections; for example, node i
A in N
e
t
w
o
r
k
A merely couples with j
B in N
e
t
w
o
r
k
B.Figure 1
Model of asymmetrical interdependent networks, in which two subnetworks are partly coupled together with proportionp
1. The solid lines in N
e
t
w
o
r
k
A and N
e
t
w
o
r
k
B represent the internal edges, while the dashed lines connecting the two subnetworks represent the coupling edges.The model can be described mathematically byG
N
,
E
,
W, where N
=
N
A
+
N
B is the set of nodes from N
e
t
w
o
r
k
A and N
e
t
w
o
r
k
B. E
=
E
A
+
E
B
+
E
A
B denotes the set of edges, where E
A (E
B) denotes the internal edges of N
e
t
w
o
r
k
A (N
e
t
w
o
r
k
B) and E
A
B represents the coupling edges between the two subnetworks. Matrix W
=
W
A
W
A
B
W
A
B
T
W
B represents the adjacency matrix, where W
A and W
B represent the connections inside N
e
t
w
o
r
k
A and N
e
t
w
o
r
k
B, respectively, while W
A
B is the coupling adjacency matrix. The element ω
i
j of W represents the weight of connection between two neighboring nodes i and j. Node intensity s
i is defined to represent the total weights of all edges connected with node i:(1)
s
i
=
∑
j
∈
Γ
i
ω
i
j
,where Γ
i is the neighborhood of node i, including neighboring nodes in both N
e
t
w
o
r
k
A and N
e
t
w
o
r
k
B.As mentioned before, inN
e
t
w
o
r
k
A, nodes are divided into several levels. There are H levels in N
e
t
w
o
r
k
A and each node from a certain level has M subordinate nodes from the next level. There are no connections between subordinate nodes in the same level except for the top level, where all nodes are connected with each other. Internal edges connecting different levels have diverse weights representing the difference in importance of interlevel connections.The hierarchical structure is established as follows. Firstly, we establish a fully connected top level, the level marker is set toh
=
1, the number of nodes is n
1, and the edge weight is ω
1. Then, each node in top level establishes relations with M subordinate nodes from the next level, where the level marker is set to h
=
2. Nodes in this level do not connect with each other. We repeat the former step until h
=
H; then, N
e
t
w
o
r
k
A with hierarchical weight is completed. The weight of edges connecting level h and level h
+
1 is set to(2)
ω
h
,
h
+
1
=
ω
1
-
C
·
h
-
1
,where constant C controls the weightiness of internal edges which connect neighboring levels.N
e
t
w
o
r
k
B is evolved by the rule of WS small-world network because small-world network yields the shortest and most effective paths [32]. The weight of each internal edge is equivalent. It turns from regular network to random network by adjusting the reconnecting probability p
2
∈
0,1, where p
2
=
0 corresponds to the case of nearest-neighbor coupled network; the nodes are placed on a ring lattice with periodic boundary conditions and each node is initially connected to K nearest neighbors. When p
2
∈
(
0,1
), the WS small-world network model is created by rewiring a small fraction of the links with probability p
2 to nodes chosen at random. The average degree of node is still K. When p
2
=
1, it corresponds to a completely random network.The coupling proportion is set top
1
=
0.7 between N
e
t
w
o
r
k
A and N
e
t
w
o
r
k
B. The sizes of the subnetworks are set to N
A
=
N
B
=
340. In N
e
t
w
o
r
k
A, M
=
4, H
=
4, and n
1
=
4, so the numbers of subordinate nodes in the next three levels are 16, 64, and 256. Because C
=
2 and ω
1
=
10, the weights of internal edges that connect different levels are ω
1,2
=
8, ω
2,3
=
6, and ω
3,4
=
4 according to (2). N
e
t
w
o
r
k
B is a WS small-world network with p
2
=
0.1 and K
=
4. Due to the randomly “one to one” relationship, we set the weight of coupling edges to be equal to the maximum weight of internal edges connected to the coupled node in N
e
t
w
o
r
k
A.


## 3. Cascading Failure Model

The coupling property makes the interdependent networks fragile when suffering from intentional attack [16, 18]. Here we focus on the cascading failure of interdependent network induced by load. We define the initial loads L
i
(
0
) of node i as a function of the node intensity s
i:(3)
L
i
0
=
λ
s
i
α
,where load parameters λ
>
0 and α
>
0 are adjustable to control the distribution of initial load. We can see that each node bears more or less some initial loads according to the node intensity and the loads on each node are a nonlinear function unless the parameter α is equal to 1. ML model [33] conjectures that the capacity of a node is proportional to its initial loads. There is a linear relationship between capacity C
i of node i and its initial load L
i
(
0
). Consider(4)
C
i
=
1
+
β
·
L
i
0
,where β
>
0 is the tolerance parameter. It is well studied that big β leads to excellent robustness. However, considering other perspectives such as cost, the capacity is impossible to be infinite, so it tends to a finite β. In addition, the linear relationship does not fit for many real-world networks, so a nonlinear “load-capacity” model is proposed as follows:(5)
C
i
=
L
i
0
+
β
·
L
i
0
θ
,where two tunable parameters β
>
0 and θ
>
0 are introduced. If θ
=
1, it decays to the linear “load-capacity” model. According to a local nearest redistribution strategy, the proportions of load distribution Π
j and the new added load Δ
L
i
→
j where the failed node i passes to its neighbor j are as follows[34]:(6)
Π
j
=
C
j
∑
n
∈
Γ
i
C
n
,
Δ
L
i
→
j
=
Π
j
L
i
=
C
j
∑
n
∈
Γ
i
C
n
L
i
,where n is the neighbor of the failed node i and Γ
i is the set of neighbors which node i connects. If the initial load of node j plus the load that node i transfers to exceeds its capacity (C
j
<
L
j
+
Δ
L
i
→
j), node j fails and leads to a new round of load redistribution. The process repeats until there is no overloaded node or the entire network is paralyzed. The evolving procedure is illustrated in Figure 2.Figure 2
The cascading failure under node attack on asymmetrical interdependent network. The initial topology of the asymmetrical network, shown in (a), has a set of nodes inN
e
t
w
o
r
k
A and N
e
t
w
o
r
k
B, labeled as A
1
,
A
2
,
…
,
A
8 and B
1
,
B
2
,
…
,
B
7, respectively. The internal edges are represented as solid black lines and the coupling edges are represented as dashed blue lines. Red circles represent the failed nodes and green circles represent the nodes that suffered from the load redistribution but are still active.
(a)
Stage 1
(b)
Stage 2
(c)
Stage 3
(d)
Stage 4At stage 1 in Figure2(a), node A
1 is attacked and fails; the load distributes to its neighbors A
2
,
A
3
,
A
4
,
A
5. Node A
2 fails because of overloading and it triggers the failure of B
1 because the load of B
1 exceeds its capacity after receiving some amount of load from A
2. Meanwhile, nodes A
3
,
A
4, and A
5 are still active because they are not overloaded. At stage 2 in Figure 2(b), due to the load redistribution of node B
1, node B
6 is overloaded, but nodes B
2 and B
5 are still active. B
6 distributes the load to its neighboring nodes which leads to the failure of nodes B
7 and A
8. In this way, the cascading failure propagates in both N
e
t
w
o
r
k
A and N
e
t
w
o
r
k
B. At stage 3 in Figure 2(c), the failures of nodes B
7 and A
8 trigger the failures of B
4 and A
5. There is no new failure in the interdependent networks in Figure 2(d) and the cascading progress ends at stage 4.We useR
s to express the relative scale damage caused by node failures:(7)
R
s
=
S
′
S
,where S (S
′) is the sum of initial node intensity before (after) node failures.Similar to the case of node failure, the initial load of edgeL
i
j
(
0
) is related to the node intensities of i and j, as defined in (8). The “load-capacity” model of edges can be expressed in (9). Figure 3 shows the load redistribution rule for edge failure. Consider(8)
L
i
j
0
=
s
i
·
s
j
γ
,where γ
>
0 is a tunable parameter to control initial load of edges. Consider(9)
C
i
j
=
L
i
j
0
+
β
·
L
i
j
0
θ
.Figure 3
The load redistribution triggered by an edge failure. The internal edges are represented as solid black lines, while the coupling edges are represented as dashed blue lines. The black arrowed lines represent the directions of load redistribution from the failed edge and the red solid line represents the edge which is overloaded.If a single internal edgee
i
A
h
A in N
e
t
w
o
r
k
A fails, the load on the broken edge will be redistributed bidirectionally to its neighboring internal and coupling edges. If the load on neighboring edges plus the extra load exceeds their capacity, this may trigger a recursive process of cascading failures, such as internal edge e
i
A
k
A and internal edge e
i
A
j
B shown in Figure 3. After coupling edge e
i
A
j
B fails, its load proportionally is reassigned to neighboring edges connected to node i
A in N
e
t
w
o
r
k
A and node j
B in N
e
t
w
o
r
k
B. Consider(10)
Δ
L
i
A
→
j
A
=
C
i
A
j
A
∑
a
∈
Γ
i
A
C
i
A
a
+
∑
b
∈
Γ
j
B
C
j
B
b
-
2
C
i
A
j
B
L
i
A
j
B
,
Δ
L
j
B
→
k
B
=
C
j
B
k
B
∑
a
∈
Γ
i
A
C
i
A
a
+
∑
b
∈
Γ
j
B
C
j
B
b
-
2
C
i
A
j
B
L
i
A
j
B
,where Γ
i
A (Γ
j
B) is the set of neighbors of node i
A (j
B) and L
i
A
j
B is the initial load of edge e
i
A
j
B. Δ
L
i
A
→
j
A is the additional load distributed to edge e
i
A
j
A. Δ
L
j
B
→
k
B is the additional load distributed to edge e
j
B
k
B. Only if C
m
n
>
L
m
n
+
Δ
L
m
→
n for any edge will there be no new edge failure in the system.We finally adopt the relative size of removed edge asR to describe the network robustness under edge attack:(11)
R
=
∑
i
j
E
i
j
E
,where E
i
j is the number of failed edges due to cascading failure after e
i
j is cut off, while E is the total number of edges in the whole system before cascading failures.


## 4. Results and Discussion

### 4.1. Analysis of Network Robustness under Node Attack

A similar phenomenon of the isolated status and interdependent status is thatR
S has positive relations with the capacity parameters β and θ. This means that high capacity benefits the robustness of N
e
t
w
o
r
k
A. Additionally, it is notable that the curve of R
S abruptly emerges when β approaches a critical value in Figure 4(a), and then it rapidly grows to 1; but, in Figure 4(b), after its abrupt emergence, the robustness indicator then shows a slow increase to 1. In order to find the detailed difference between the two statuses, we fix θ at 0.2, 0.4, and 0.6 separately to compare the curves of R
S in isolated status and interdependent status in Figure 4(c). It can be seen that the solid lines are above the corresponding dashed ones because the interdependent relationship leads to a weaker robustness of N
e
t
w
o
r
k
A. Furthermore, the robustness indicator R
S in isolated status undergoes a short “plateau” state before it grows to 1. In comparison, in the interdependent status, the robustness undergoes a hybrid phase transition.Figure 4
The relationship between capacity and robustness ofN
e
t
w
o
r
k
A under node intentional attack. (a) N
e
t
w
o
r
k
A in isolated status. (b) N
e
t
w
o
r
k
A in interdependent status. (c) A comparison of N
e
t
w
o
r
k
A in both isolated status and interdependent status. Capacity parameters β
∈
[
0,2
) and θ
∈
[
0.2,1.6
]. Load parameter α
=
0.8. In (c), the critical value β
C is labeled with green dot. The solid lines denote the isolated status, while the dashed lines denote the interdependent status. Simulation results are averaged over 1000 independent trials.
(a)
(b)
(c)The robustness ofN
e
t
w
o
r
k
A has a positive relation with the capacity parameter β in Figure 5. However, it has a negative relationship with load parameter α as shown in Figures 5(a) and 5(b). We fix α at 1.2, 1.4, and 1.6 separately to compare the changes of robustness in both statuses. We also find that the robustness of N
e
t
w
o
r
k
A in isolated status is better than that in interdependent status. In addition, the trend of each curve in Figure 5(c) is similar to that in Figure 4(c). Because of the hierarchical weight structure in N
e
t
w
o
r
k
A, nodes from the same level have equal initial load and capacity according to (1) and (3). If a node fails, the proportion of load redistribution is divided by the level to which the neighboring nodes belong. The failure probabilities of neighboring nodes from the same level are equivalent, while those of the nodes from different levels are not. That finally leads to the sharp change of R
S when most neighbors fail or the short plateau without large scale of failures in the neighborhood. The hierarchical weight structure in N
e
t
w
o
r
k
A also leads to a cascade of failures in the isolated status unlike a regular first-order phase transition. Meanwhile, the interdependent relationship leads to the failures from a first-order phase transition to a second-order phase transition.Figure 5
The relationship between load and robustness ofN
e
t
w
o
r
k
A under node intentional attack. (a) N
e
t
w
o
r
k
A in isolated status. (b) N
e
t
w
o
r
k
A in interdependent status. (c) A comparison of N
e
t
w
o
r
k
A in both isolated status and interdependent status. Capacity parameters β
∈
[
0,2
) and θ
=
0.8. The critical value β
C is labeled with green dot. The solid lines denote the robustness of N
e
t
w
o
r
k
A in isolated status, while the dashed lines denote the robustness of N
e
t
w
o
r
k
A in interdependent status. Simulation results are averaged over 1000 independent trials.
(a)
(b)
(c)


### 4.2. Analysis of Network Robustness under Edge Attack

When the edge with maximum initial loadm
a
x
(
L
i
j
(
0
)
) in N
e
t
w
o
r
k
A suffers from intentional attack, the relationship among load, capacity, and robustness is analyzed based on the results in Figures 6 and 7. Under edge attack, the robustness indicator R also shows positive correlations with the capacity parameters β and θ but a negative correlation with load parameter γ. The results shown in Figures 6 and 7 also prove that the robustness can be enhanced by high capacity as well as low initial load. As shown in Figures 6(a) and 6(b) and Figures 7(a) and 7(b), the enhancement of robustness is more prompt in isolated status than in interdependent status. In Figures 6(c) and 7(c), the solid line denoting each value of θ and γ is above the corresponding dashed ones which also indicates a weaker robustness of N
e
t
w
o
r
k
A in interdependent status. The hierarchical weight structure in N
e
t
w
o
r
k
A also leads to a plateau on the curves of R in Figure 7(c). Because the initial load and capacity of internal edge L
i
j which connect the neighboring levels are equal, according to (1), (2), and (8). Edges connecting the neighboring levels have equal initial load and capacity. If one internal edge fails, the failure probabilities of edges connecting the same neighboring levels are equivalent, while those of the edges connecting different levels are not. So it leads to the sharp change of R when most edges fail or the plateau without large scale of edge failures.Figure 6
The relationship between capacity and robustness ofN
e
t
w
o
r
k
A under edge intentional attack. (a) N
e
t
w
o
r
k
A in isolated status. (b) N
e
t
w
o
r
k
A in interdependent status. (c) A comparison of N
e
t
w
o
r
k
A in both isolated status and interdependent status. Capacity parameters β
∈
[
0,2
) and θ
∈
[
0.2,1.6
]. Load parameter γ
=
0.8. The solid lines denote the isolated status, while the dashed lines denote the interdependent status. Simulation results are averaged over 1000 independent trials.
(a)
(b)
(c)Figure 7
The relationship between load and robustness under edge intentional attack in isolated status (a) and in interdependent status (b). (c) A comparison ofN
e
t
w
o
r
k
A in both isolated status and interdependent status. Capacity parameters β
∈
[
0,2
) and γ
∈
[
0.2,1.6
]. Capacity parameter θ
=
0.8. The solid lines denote the robustness of N
e
t
w
o
r
k
A in isolated status, while the dashed lines denote the robustness of N
e
t
w
o
r
k
A in interdependent status. Simulation results are averaged over 1000 independent trials.
(a)
(b)
(c)


### 4.3. Comparison of Asymmetrical and Symmetrical Models

We compare between robustness of the asymmetrical network model we propose and that of a symmetrical network model (WS-WS). The configuration of WS subnets in the symmetrical network model is the same as that of the WS subnet in the asymmetrical network model. The robustness of network under different initial conditions is shown in Figures8 and 9.Figure 8
Comparison of node attack in symmetrical and asymmetrical network models under different initial conditions. (a)α
=
0.8 and θ
=
0.6. (b) α
=
0.6 and θ
=
0.8. The red lines without triangles denote the robustness in WS-WS network model, while solid lines with triangles denote the robustness in N
e
t
w
o
r
k
B. Simulation results are averaged over 1000 independent trials.
(a)
(b)Figure 9
Comparison of edge attack in symmetrical and asymmetrical network models under different initial conditions. (a)γ
=
0.8 and θ
=
0.6. (b) γ
=
0.6 and θ
=
0.8. The red lines without triangles denote the robustness in WS-WS network model, while solid lines with triangles denote the robustness in N
e
t
w
o
r
k
B. Simulation results are averaged over 1000 independent trials.
(a)
(b)We test intentional node attacks on the studied models under two initial conditions in Figure8. It can be seen that the robustness of the symmetrical model is stronger than that of the asymmetric model. Additionally, no matter in the symmetrical or the asymmetrical case, the robust levels of N
e
t
w
o
r
k
A and N
e
t
w
o
r
k
B under the same conditions are very close to one another.Additionally, we test intentional edge attacks on the two models under different initial conditions in Figure9. Similar to the cases on node attack, the robust levels in the symmetrical model are better than those in the asymmetric model. However, what is notable is that R
-
β curves of each subnetwork in the asymmetrical model deviate from each other, while those of the symmetrical model remain resemblant. N
e
t
w
o
r
k
A is less vulnerable than N
e
t
w
o
r
k
B.


## 4.1. Analysis of Network Robustness under Node Attack

A similar phenomenon of the isolated status and interdependent status is thatR
S has positive relations with the capacity parameters β and θ. This means that high capacity benefits the robustness of N
e
t
w
o
r
k
A. Additionally, it is notable that the curve of R
S abruptly emerges when β approaches a critical value in Figure 4(a), and then it rapidly grows to 1; but, in Figure 4(b), after its abrupt emergence, the robustness indicator then shows a slow increase to 1. In order to find the detailed difference between the two statuses, we fix θ at 0.2, 0.4, and 0.6 separately to compare the curves of R
S in isolated status and interdependent status in Figure 4(c). It can be seen that the solid lines are above the corresponding dashed ones because the interdependent relationship leads to a weaker robustness of N
e
t
w
o
r
k
A. Furthermore, the robustness indicator R
S in isolated status undergoes a short “plateau” state before it grows to 1. In comparison, in the interdependent status, the robustness undergoes a hybrid phase transition.Figure 4
The relationship between capacity and robustness ofN
e
t
w
o
r
k
A under node intentional attack. (a) N
e
t
w
o
r
k
A in isolated status. (b) N
e
t
w
o
r
k
A in interdependent status. (c) A comparison of N
e
t
w
o
r
k
A in both isolated status and interdependent status. Capacity parameters β
∈
[
0,2
) and θ
∈
[
0.2,1.6
]. Load parameter α
=
0.8. In (c), the critical value β
C is labeled with green dot. The solid lines denote the isolated status, while the dashed lines denote the interdependent status. Simulation results are averaged over 1000 independent trials.
(a)
(b)
(c)The robustness ofN
e
t
w
o
r
k
A has a positive relation with the capacity parameter β in Figure 5. However, it has a negative relationship with load parameter α as shown in Figures 5(a) and 5(b). We fix α at 1.2, 1.4, and 1.6 separately to compare the changes of robustness in both statuses. We also find that the robustness of N
e
t
w
o
r
k
A in isolated status is better than that in interdependent status. In addition, the trend of each curve in Figure 5(c) is similar to that in Figure 4(c). Because of the hierarchical weight structure in N
e
t
w
o
r
k
A, nodes from the same level have equal initial load and capacity according to (1) and (3). If a node fails, the proportion of load redistribution is divided by the level to which the neighboring nodes belong. The failure probabilities of neighboring nodes from the same level are equivalent, while those of the nodes from different levels are not. That finally leads to the sharp change of R
S when most neighbors fail or the short plateau without large scale of failures in the neighborhood. The hierarchical weight structure in N
e
t
w
o
r
k
A also leads to a cascade of failures in the isolated status unlike a regular first-order phase transition. Meanwhile, the interdependent relationship leads to the failures from a first-order phase transition to a second-order phase transition.Figure 5
The relationship between load and robustness ofN
e
t
w
o
r
k
A under node intentional attack. (a) N
e
t
w
o
r
k
A in isolated status. (b) N
e
t
w
o
r
k
A in interdependent status. (c) A comparison of N
e
t
w
o
r
k
A in both isolated status and interdependent status. Capacity parameters β
∈
[
0,2
) and θ
=
0.8. The critical value β
C is labeled with green dot. The solid lines denote the robustness of N
e
t
w
o
r
k
A in isolated status, while the dashed lines denote the robustness of N
e
t
w
o
r
k
A in interdependent status. Simulation results are averaged over 1000 independent trials.
(a)
(b)
(c)


## 4.2. Analysis of Network Robustness under Edge Attack

When the edge with maximum initial loadm
a
x
(
L
i
j
(
0
)
) in N
e
t
w
o
r
k
A suffers from intentional attack, the relationship among load, capacity, and robustness is analyzed based on the results in Figures 6 and 7. Under edge attack, the robustness indicator R also shows positive correlations with the capacity parameters β and θ but a negative correlation with load parameter γ. The results shown in Figures 6 and 7 also prove that the robustness can be enhanced by high capacity as well as low initial load. As shown in Figures 6(a) and 6(b) and Figures 7(a) and 7(b), the enhancement of robustness is more prompt in isolated status than in interdependent status. In Figures 6(c) and 7(c), the solid line denoting each value of θ and γ is above the corresponding dashed ones which also indicates a weaker robustness of N
e
t
w
o
r
k
A in interdependent status. The hierarchical weight structure in N
e
t
w
o
r
k
A also leads to a plateau on the curves of R in Figure 7(c). Because the initial load and capacity of internal edge L
i
j which connect the neighboring levels are equal, according to (1), (2), and (8). Edges connecting the neighboring levels have equal initial load and capacity. If one internal edge fails, the failure probabilities of edges connecting the same neighboring levels are equivalent, while those of the edges connecting different levels are not. So it leads to the sharp change of R when most edges fail or the plateau without large scale of edge failures.Figure 6
The relationship between capacity and robustness ofN
e
t
w
o
r
k
A under edge intentional attack. (a) N
e
t
w
o
r
k
A in isolated status. (b) N
e
t
w
o
r
k
A in interdependent status. (c) A comparison of N
e
t
w
o
r
k
A in both isolated status and interdependent status. Capacity parameters β
∈
[
0,2
) and θ
∈
[
0.2,1.6
]. Load parameter γ
=
0.8. The solid lines denote the isolated status, while the dashed lines denote the interdependent status. Simulation results are averaged over 1000 independent trials.
(a)
(b)
(c)Figure 7
The relationship between load and robustness under edge intentional attack in isolated status (a) and in interdependent status (b). (c) A comparison ofN
e
t
w
o
r
k
A in both isolated status and interdependent status. Capacity parameters β
∈
[
0,2
) and γ
∈
[
0.2,1.6
]. Capacity parameter θ
=
0.8. The solid lines denote the robustness of N
e
t
w
o
r
k
A in isolated status, while the dashed lines denote the robustness of N
e
t
w
o
r
k
A in interdependent status. Simulation results are averaged over 1000 independent trials.
(a)
(b)
(c)


## 4.3. Comparison of Asymmetrical and Symmetrical Models

We compare between robustness of the asymmetrical network model we propose and that of a symmetrical network model (WS-WS). The configuration of WS subnets in the symmetrical network model is the same as that of the WS subnet in the asymmetrical network model. The robustness of network under different initial conditions is shown in Figures8 and 9.Figure 8
Comparison of node attack in symmetrical and asymmetrical network models under different initial conditions. (a)α
=
0.8 and θ
=
0.6. (b) α
=
0.6 and θ
=
0.8. The red lines without triangles denote the robustness in WS-WS network model, while solid lines with triangles denote the robustness in N
e
t
w
o
r
k
B. Simulation results are averaged over 1000 independent trials.
(a)
(b)Figure 9
Comparison of edge attack in symmetrical and asymmetrical network models under different initial conditions. (a)γ
=
0.8 and θ
=
0.6. (b) γ
=
0.6 and θ
=
0.8. The red lines without triangles denote the robustness in WS-WS network model, while solid lines with triangles denote the robustness in N
e
t
w
o
r
k
B. Simulation results are averaged over 1000 independent trials.
(a)
(b)We test intentional node attacks on the studied models under two initial conditions in Figure8. It can be seen that the robustness of the symmetrical model is stronger than that of the asymmetric model. Additionally, no matter in the symmetrical or the asymmetrical case, the robust levels of N
e
t
w
o
r
k
A and N
e
t
w
o
r
k
B under the same conditions are very close to one another.Additionally, we test intentional edge attacks on the two models under different initial conditions in Figure9. Similar to the cases on node attack, the robust levels in the symmetrical model are better than those in the asymmetric model. However, what is notable is that R
-
β curves of each subnetwork in the asymmetrical model deviate from each other, while those of the symmetrical model remain resemblant. N
e
t
w
o
r
k
A is less vulnerable than N
e
t
w
o
r
k
B.


## 5. Conclusion

Interdependent network is a representative abstraction of realistic systems, where failures in one system may trigger a cascade of failures in the other system through interactions among systems. For instance, load redistribution caused by node or edge failures in one subnetwork may trigger an iterative progress of cascading failure in interdependent networks. In this paper, firstly, we propose an interdependent network model considering asymmetrical property of the two subnetworks. Secondly, a nonlinear “load-capacity” model is also proposed for nodes and edges, respectively, to model the cascading failures. Thirdly, we conduct a systematic series of experiments to analyze the robustness of interdependent network under attack and make comparisons with the situation of well-studied symmetrical WS-WS network model. (1) Interdependent networks are more fragile compared with isolated ones. (2) The hierarchical weight structure leads to a short “plateau” phenomenon in isolated status unlike the regular first-order phase transition in nonhierarchical isolated networks. Under node attack, only if the capacity reaches a critical valueβ
C can the robustness abruptly emerge and rapidly approach 1. Furthermore, the interconnections lead to the robustness curve from a first-order phase transition to a second-order phase transition in interdependent networks. (3) Under edge attack, the hierarchical weight structure also causes a plateau in an iterative process of cascading failures in isolated status. The change of robustness in isolated network is more prompt than that in interdependent network. (4) Compared with symmetrical WS-WS network model, the robustness of the asymmetrical model is weaker. Under node attack, no matter in the symmetrical or the asymmetrical case, the two subnetworks have similar robustness. However, under edge attack, subnetworks in the asymmetrical model perform differently, while those of the symmetrical model remain resemblant.



---

*Source: 194568-2015-10-25.xml*