
Identifying influential spreaders in complex networks is a crucial issue that can help control the propagation process in complex networks. An aviation network is a typical complex network, and accurately identifying the key city nodes in the aviation network can help us better prevent network attacks and control the spread of diseases. In this paper, a method for identifying key nodes in undirected weighted networks, called weighted Laplacian energy centrality, was proposed and applied to an aviation network constructed from real flight data. Based on the analysis of the topological structure of the network, the paper recognized critical cities in this network, then simulation experiments were conducted on key city nodes from the perspectives of network dynamics and robustness. The results indicated that, compared with other methods, weighted Laplacian energy centrality can identify the city nodes with the most spreading influence in the network. From the perspective of network robustness, the identified key nodes also have the characteristics of accurately and quickly destroying network robustness.
Citation: Shuying Zhao, Shaowei Sun. A study on centrality measures in weighted networks: A case of the aviation network[J]. AIMS Mathematics, 2024, 9(2): 3630-3645. doi: 10.3934/math.2024178
[1] | Yike Wang, Gaoxia Wang, Ximei Hou, Fan Yang . Motif adjacency matrix and spectral clustering of directed weighted networks. AIMS Mathematics, 2023, 8(6): 13797-13814. doi: 10.3934/math.2023706 |
[2] | Ling Zhou, Zuhan Liu . Blow-up in a p-Laplacian mutualistic model based on graphs. AIMS Mathematics, 2023, 8(12): 28210-28218. doi: 10.3934/math.20231444 |
[3] | Xiaodie Lv, Yi Liu, Yihua Zhong . A novel method to solve the optimization problem of uncertain network system based on uncertainty theory. AIMS Mathematics, 2023, 8(3): 5445-5461. doi: 10.3934/math.2023274 |
[4] | Bin Zhou, Xiujuan Ma, Fuxiang Ma, Shujie Gao . Robustness analysis of random hyper-networks based on the internal structure of hyper-edges. AIMS Mathematics, 2023, 8(2): 4814-4829. doi: 10.3934/math.2023239 |
[5] | Fan Yang, Xiaohui Ai . Exponential stability of stochastic complex networks with multi-weights driven by second-order process based on graph theory. AIMS Mathematics, 2024, 9(4): 9847-9866. doi: 10.3934/math.2024482 |
[6] | Yongyan Zhao, Jian Li . Integrating artificial intelligence with network evolution theory for community behavior prediction in dynamic complex systems. AIMS Mathematics, 2025, 10(2): 2042-2063. doi: 10.3934/math.2025096 |
[7] | Gonzalo Contreras-Aso, Cristian Pérez-Corral, Miguel Romance . Uplifting edges in higher-order networks: Spectral centralities for non-uniform hypergraphs. AIMS Mathematics, 2024, 9(11): 32045-32075. doi: 10.3934/math.20241539 |
[8] | Patarawadee Prasertsang, Thongchai Botmart . Improvement of finite-time stability for delayed neural networks via a new Lyapunov-Krasovskii functional. AIMS Mathematics, 2021, 6(1): 998-1023. doi: 10.3934/math.2021060 |
[9] | Wanchun Fan, Yan Jiang, Songyang Huang, Weiguo Liu . Research and prediction of opioid crisis based on BP neural network and Markov chain. AIMS Mathematics, 2019, 4(5): 1357-1368. doi: 10.3934/math.2019.5.1357 |
[10] | Muhammad Ahmad, Muhammad Faheem, Sanaa A. Bajri, Zohaib Zahid, Muhammad Javaid, Hamiden Abd El-Wahed Khalifa . Optimizing SNARK networks via double metric dimension. AIMS Mathematics, 2024, 9(8): 22091-22111. doi: 10.3934/math.20241074 |
Identifying influential spreaders in complex networks is a crucial issue that can help control the propagation process in complex networks. An aviation network is a typical complex network, and accurately identifying the key city nodes in the aviation network can help us better prevent network attacks and control the spread of diseases. In this paper, a method for identifying key nodes in undirected weighted networks, called weighted Laplacian energy centrality, was proposed and applied to an aviation network constructed from real flight data. Based on the analysis of the topological structure of the network, the paper recognized critical cities in this network, then simulation experiments were conducted on key city nodes from the perspectives of network dynamics and robustness. The results indicated that, compared with other methods, weighted Laplacian energy centrality can identify the city nodes with the most spreading influence in the network. From the perspective of network robustness, the identified key nodes also have the characteristics of accurately and quickly destroying network robustness.
Many real-world networks can be regarded as complex networks, such as social networks [1] and transportation networks [2]. The development of the complex systems theory has provided a strong theoretical foundation for the analysis of complex networks in real life. Studies in [3,4,5] have shown that the failure of a few nodes in a network often leads to complex cascading failure effects, even causing widespread network paralysis. Therefore, the identification of key nodes in networks has become a popular and important research topic. Currently, the study of network centrality is mostly based on the perspective of network topology [6,7,8]. The topology of a network contains a lot of functional information about the network. Starting from considering different degrees of topological structure information, the recognition method is extended from locality to globality. The simplest local method is degree centrality [9], which is simple and efficient, and suitable for large-scale networks. Degree centrality applied to weighted networks is known as strength centrality [10], which defines the strength of a node in terms of the total weight of its connection. The most classic global centrality methods are closeness centrality [11] and betweenness centrality [12]. Nodes with high closeness centrality have the best perspective to observe the information propagation in the network, while nodes with high betweenness centrality can directly influence the flow of information in the network. When applying these two methods to weighted networks, the calculation of the shortest paths was defined correspondingly as the weighted shortest paths in [13]. In recent years, Ma et al. [14] has applied graph energy to identify central nodes in networks. In [15], we proposed the application of Laplacian energy centrality on undirected unweighted networks, which showed promising results. However, this method has not been extended to weighted graphs.
The aviation network is a complex network formed by cities as nodes and airlines between them as edges. The aviation network serves as a platform for social communication and economic development. Events such as adverse weather conditions, public health emergencies, and terrorist attacks can result in airport closures or flight suspensions, impacting the transportation efficiency of the entire aviation network and causing significant economic losses. Therefore, identifying key nodes in the aviation network has become a crucial issue. The theoretical analysis methods of complex networks have effectively matched with the study of the topological structure of networks [16,17,18,19]. In [20,21,22,23], a global aviation transportation network model was established and the analysis of this model shows that the air transport network exhibits small world characteristics and scale-free properties. Key city nodes in transportation networks often have a significant impact on the operational performance of the network. Shen et al.[19] used principal component analysis (PCA) to identify key transportation cities in the transportation network and improve its performance. Lordan et al.[24] studied the robustness of the global aviation transportation network using complex network theory, proposed a method to determine key airports and intentionally attack the network based on different importance indicators, and compared the differences in the corresponding evaluation results. This method can accurately identify the central node and open some paths to future research in this area. Mo et al. [25] introduced classical centrality methods to analyze airport nodes based on the analysis of the aviation network structure, showing that traditional centrality methods can effectively identify the airport's structural system, but the paper did not consider the weight of the network. Li et al.[26] applied the concept of minimum connected dominating sets to complex networks. This method can simultaneously identify both key nodes and edges in the network, which is of practical significance for studying network resilience and backbone network construction. This method is highly integrated with the actual situation of air transportation, but the experimental data in this paper has limitations and the applicability of the theory is not well verified. Lou et al. [27] analyzed the changes in the robustness of network structures after attacks on different airline companies. However, the aforementioned literature only focused on network robustness and did not verify the identified key nodes from the perspective of network dynamics.
Based on the above research, we will extend Laplacian energy centrality to weighted networks and apply it to identify key nodes in the Chinese aviation network in this paper. We validate the identified key nodes from the perspectives of network propagation dynamics and network robustness. The rest of this paper is structured as follows: Section 2 provides a detailed exposition of the methods employed for network construction, the topology of the network, and several classical centrality measures. Following that, Section 3 focuses on the principles of the weighted Laplacian energy centrality method. In Section 4, we validate the weighted Laplacian energy centrality from the perspectives of influence and robustness. We give a conclusion in the last section.
In this section, we first introduce the way to construct the Chinese aviation network. Subsequently, the constructed aviation network and its topological structure are introduced, with the analysis of the network's topological structure being able to uncover various information about the network. Finally, other centrality measures that serve as comparisons to the methods proposed in this paper are listed.
The flight data between various cities of China from January 1–7, 2022 is obtained by an app called 'Ctrip'. The aviation network is abstracted as G=(V,E,W), where the node set V={v1,v2,v3,…,vn} represents all cities with air transport. If a city has two or more airports, merge them into one node. The edge set E represents the connections between cities; if there are flights between two cities vi and vj, then an edge vivj exists between these two cities. The weight of an edge vivj, denoted by wi,j, is the total number of round-trip flights between these two cities vi and vj from January 1–7, 2022, excluding the self-loop formed by routes between two different airports in the same city. The weight set W contains the weights of all edges. The aviation network constructed in this paper consists of 198 nodes and 2379 edges. The topological structure of the flight data on January 1, 2022, is shown in Figure 1.
The node degree is an important indicator that reflects the statistical characteristics of the interconnections between nodes in a network. The degree ki of a node vi is defined as the sum of the number of edges connected to the node vi. The degree of a node refers to the number of connections the node has in a network. The five largest cities in terms of node degrees in the aviation network shown in Figure 1 are Beijing, Shanghai, Chengdu, Guangzhou, and Shenzhen. It can be seen that these are the major hub cities in China.
Degree distribution[28] is the proportion of nodes with degree k in the network, denoted as p(k). Figure 2 shows the degree distribution of the aviation network depicted in Figure 1. From Figure 2, it can be observed that the degree distribution of this network approximates a power law distribution. This is consistent with the fact that most airports have very few routes, while a few large hub airports directly connect to hundreds of other airports.
In an unweighted network, the average path length is defined as the average of the shortest paths between any two nodes in the network, calculated as follows:
L=1N(N−1)∑i≠jdij, | (2.1) |
where N represents the number of nodes in the network and dij represents the length of the shortest path between nodes vi and vj. A lower L value indicates better connectivity in the network.
Calculating the shortest path between nodes in a weighted network requires consideration of weights, where the weight between two neighboring nodes indicates the length of the path between the two nodes. According to the definition of the average shortest path, a smaller average shortest path indicates better network connectivity, but in the aviation network, if we use the number of flight routes as the weight between two adjacent nodes, larger values indicate more flight routes, which means better connectivity. Thus in this paper, we normalize the weights of the network, and then use the method of taking the reciprocal of the weights to make the weights between strong connections smaller than those between weak connections. This normalization method was introduced in [29] as follows:
Wi,j=wi,j∑vivj∈E(G)wi,jm, | (2.2) |
where Wi,j represents the normalized weight value, and wi,j denotes the original weight value between node vi and node vj before normalization. m represents the total number of edges in network G.
The calculation method for the weighted shortest path length swij [30] between node vi and node vj is as follows:
swij=1Wi,c+…+1Wd,j, | (2.3) |
where c and d refer to the nodes traversed by the shortest path between node vi and node vj, then the average path length of a weighted network with N nodes in this paper is calculated as follows:
Lw=1N(N−1)∑i≠jswij. | (2.4) |
The average path length of the Chinese aviation network in Figure 1 is calculated to be 2.5097. This relatively small value indicates that we require only a few transfers to reach any city, which satisfies the transportation needs of our airline at the current stage.
The clustering coefficient of a network refers to the proportion of actual connections among the neighbors of a node in the network. It is commonly used to measure the density of connections between nodes and the strength of community structures within a network. Let Ei be the number of edges that exist among the neighbors of the node vi. The clustering coefficient Ci for a node vi is defined as Ei divided by the maximum number of possible edges, that is,
Ci=2Eiki(ki−1). | (2.5) |
A network with a relatively high average clustering coefficient means that the nodes in that network are more closely connected to each other, forming a stronger community structure. The clustering coefficient C of the network refers to the average value of the clustering coefficients of all nodes in the network. Note that C∈[0,1]. The clustering coefficient of the network shown in Figure 1 is calculated to be C=0.7365, which indicates that the network is more closely connected, and because the weighted average path length of this network is relatively small, we conclude that the aviation network in our country exhibits characteristics of small-world networks.
In weighted networks, the importance of a node vi is determined by the sum of the weights of the edges incident to node vi. The formula for strength centrality[10] is as follows:
DC(vi)=∑vj∈Niwi,j, | (2.6) |
where wi,j represents the weight between nodes vi and vj and Ni represents the neighbor sets of vi.
Betweenness centrality[30] takes this position by giving higher centrality values to the nodes that fall within the shortest path of many pairs of nodes. In simpler terms, nodes with higher betweenness centrality often have a direct influence on the flow of information in the network. In this study on the aviation network, a shortest path between two nodes is determined by the minimum weighted shortest path lengths between them. The betweenness centrality of a node vi is represented by the following equation:
BC(vi)=∑k≠j≠iδkj(i)δkj, | (2.7) |
where δkj represents the number of weighted shortest paths between vk and vj, and δkj(i) represents the number of weighted shortest paths between vk and vj that pass through the node vi.
The definition of closeness centrality[13] states that the closer a node is to the rest of the nodes in the network in terms of the weighted average distance, the faster information can spread throughout the network. The normalized value of closeness centrality essentially represents the inverse of the distance and can be expressed using the following equation:
CC(vi)=N−1∑i≠jswij, | (2.8) |
where N is the number of nodes in the network and swij is defined before.
The eigenvector centrality[31] suggests that the importance of a node depends on both the number of its neighboring nodes and the importance of those neighboring nodes. The more important the neighboring nodes that are connected to a particular node, the more important that node is considered. The calculation formula is as follows:
EC(vi)=xi=cN∑j=1aijxj, | (2.9) |
where c is a proportionality constant, x=(x1,x2,…,xN)T, and after multiple iterations to reach a stable state, it can be written in matrix form as follows:
x=cAx. |
Here, A=(aij) refers to the adjacency matrix of the network, and x is the eigenvector corresponding to the eigenvalue c−1 of matrix A.
In this section, we will extend a centrality method, called Laplacian energy centrality [15], to weighted networks. For this, we first give some definitions.
Let G=(V,E,W) be a weighted network with n nodes and m edges, without loops and multi-edges, where the node set is V(G)={v1,v2,v3,…,vn} and the edge set is E(G)={(e1,e2,…,em)}. Any edge e=(vivj) in E(G) has a weight value wi,j which is contained in set W(G). It is clear that wi,j=wj,i and wi,i=0. Then the adjacency matrix of G is defined below:
A(G)=(0w1,2…w1,nw2,10…w2,n⋅⋅…⋅wn,1wn,2…0). |
For each row i, we define its sum as si=∑ni=1wi,j=∑vj∈Niwi,j, where Ni is the set of neighboring nodes of node vi. si represents the weighted sum of node vi. The degree matrix of G is defined by D(G)=diag{s1,s2,…,sn}.
The Laplacian matrix of the weighted network G is defined as L(G)=D(G)−A(G). Some well-known properties of the Laplacian matrix L(G) are listed as follows:
● L(G) is symmetric, singular, and positive semi-definite;
● All eigenvalues are real and nonnegative;
● The smallest eigenvalue is always 0.
The eigenvalues of L(G) can be arranged as μ1(G)≥μ2(G)≥⋯≥μn(G)=0. Using the spectral features of L(G), the third Laplacian energy of G is defined as:
E3L(G)=n∑i=1μ3i(G). |
The third Laplacian energy centrality LCw(vi) of a node vi in the weighted network G (the unweighted version was defined in [15]) is then defined as
LCw(vi)=E3L(G)−E3L(G−vi). |
In this paper, we call the above centrality weighted Laplacian energy centrality. For the reason of choosing the third power of eigenvalues, readers can refer to the paper [15]. Next, we give an expression of LCw(vi) by the local information of the node vi.
Theorem 3.1. (weighted Laplacian energy centrality): Let G be a weighted network with n nodes. For a node vi∈V(G), we have
LCw(vi)=s3i+∑vj∈N(i)(3siw2i,j+3wi,js2j−2w3i,j+3wi,jn∑k=1w2j,k)+6Δwi, | (3.1) |
where si is the strength of node vi and Δwi is the sum of the weights of the triangles containing the node vi (the weight of a triangle means the product of the weights of its three edges).
Proof. Let Δw be the sum of weights of all triangles in G. From the definition of E3L(G), we have
E3L(G)=tr((D−A)3)=tr(D3−D2A−DAD+DA2−AD2+ADA+A2D−A3). |
Since tr(D2A)=tr(DAD)=tr(AD2)=0, tr(DA2)=tr(ADA)=tr(A2D)=∑nk=1sk∑nj=1w2k,j and tr(A3)=6Δw, from the above equation, we obtain
E3L(G)=n∑k=1(s3k+3skn∑j=1w2k,j)−6Δw. |
Therefore,
LCw(vi)=E3L(G)−E3L(G−vi)=∑vj∈Ni(s3j−(sj−wi,j)3+3sjn∑k=1w2j,k−3(sj−wi,j)(n∑k=1w2j,k−w2i,j))+s3i+3sin∑j=1w2i,j++6Δwi=∑vj∈Ni(wi,j(s2j+(sj−wi,j)2+sj(sj−wi,j))+3wi,jn∑k=1w2j,k+3(sj−wi,j)w2i,j)+s3i+3sin∑j=1w2i,j+6Δwi=s3i+∑vj∈Ni3siw2i,j+∑vj∈Ni(3wi,js2j+w3i,j−3sjw2i,j+3wi,jn∑k=1w2j,k+3sjw2i,j−3w3i,j)+6Δwi=s3i+∑vj∈Ni(3siw2i,j+3wi,js2j−2w3i,j+3wi,jn∑k=1w2j,k)+6Δwi. |
If we only consider the case that all weights of edges in the network are nonnegative, then from the above theorem with the fact that 3siw2i,j−2w3i,j≥0, we conclude that LCw(vi)≥0 for any node vi. Moreover, by the above theorem, one can easily obtain an expression of Laplacian energy centrality for unweighted networks, which is given in [15]. In the rest of the paper, LCw is simplified to LC for short.
The performance of the LC is evaluated by network dynamics and robustness for assessing central nodes in the network. We will present experimental results showcasing the performance of weighted Laplacian energy centrality and a series of other methods under these metrics.
In this paper, we assess the importance of identified key nodes by evaluating their significance from two perspectives: The influence of central nodes on information propagation and the degree of change in network robustness after being subjected to attacks. The SIR model (Susceptible-Infected-Recovered, which will be defined in next section) is used to simulate information propagation, while the primary indicators of network robustness include the relative size of the maximum connectivity subgraph and the network efficiency.
To analyze the failure process of the network under dynamic conditions, the most important thing is to study the propagation dynamics model, which can effectively analyze the dynamic connection between various factors in the network, and derive the law followed by the system under dynamic conditions. In recent years, some scholars have used the SIR model [32] to simulate the spread of disease in the global aviation network, and have given a strategy to curb the spread and provide a set of risk assessment systems. The SIR model can be used to simulate the spread of disease in the aviation network and the impact of flight delays.
In the SIR model, the nodes can be in one of three states: susceptible (S), infected (Ⅰ), or recovered (R). In the initial stage, the states of individual nodes in the network are established. At each time step, nodes in the infected state (Ⅰ) attempt to infect susceptible nodes (S) with an infection rate β, while recovering to the immune state (R) with a certain probability γ. The recovered nodes become immune and cannot be infected or infect others. The propagation process concludes when there are no nodes in the infected state (Ⅰ) present in the network, then we have
{dSdt=−βSIndIdt=βSIn−γIdRdt=γI | (4.1) |
with n=S+I+R, where S, I and R represent the number of susceptible nodes, infected nodes, and recovered nodes at time t, respectively.
The robustness of a complex network refers to the ability of the network to resist damage when it suffers different degrees of damage[33]. The robustness of the aviation network studied in this paper refers to the ability of the network to maintain its overall transportation function when natural disasters, public health emergencies, terrorist attacks, and other emergencies trigger route disruptions or airport closures, specifically whether it can reach the final destination at the time of the attack and ensure transportation efficiency as much as possible while ensuring that it reaches the destination. Consider that an attack can change the structure and transportation efficiency of the network. The following section introduces two of the most classical measures of maximum connected subgraph relative size and weighted network efficiency to measure the robustness of networks.
The maximum connectivity subgraph is the maximum connectivity component split after a network is attacked, which is an index reflecting the connectivity of the network. The relative size of the maximum connectivity subgraph S is defined as the ratio of the number of nodes in the maximum connectivity subgraph of a network after it is attacked to the total number of nodes in the original network, which is calculated as follows:
S=N′N, | (4.2) |
where N′ is the number of nodes in the largest connected component after the attack, and N is the number of nodes in the original network.
Network efficiency refers to the effect of structural changes on the shortest path distance between nodes after a network has been attacked. For unweighted networks, the shortest path refers to a path with the least number of edges between two nodes vi and vj is the number of edges of the shortest path length dij between two nodes. If there is no path between two nodes in the network, the shortest path length is infinite, and the efficiency between two nodes is expressed by the inverse of the shortest path length between two nodes, and its inverse is 0 which does not affect the result of the calculation. The efficiency of the entire network is defined as the average efficiency between all nodes denoted by E. The calculation is as follows:
E=1N(N−1)∑i≠j1dij. | (4.3) |
In the case of a weighted network, the calculation of the shortest path length between two nodes considers the weights assigned to the edges, denoted as swij. Therefore, the weighted network efficiency Ew is defined as:
Ew=1N(N−1)∑i≠j1swij. | (4.4) |
This paper simulates the spread process of the aviation network using the SIR model. The initial state of the network consists of the top 10 ranked nodes obtained through a certain method, as shown in Table 1, which serves as the initial infected node, and the remaining nodes are susceptible. In this aviation network, if we use an infection threshold βc[34] as the probability of infection, then the probability of infection between any two adjacent nodes is the same. However, since the edges of the network have weights, the edge with a large weight value means that the more round-trip flights between the two cities, then the probability of propagation of information between the two cities should be different according to the weight value of the edge, and the infection probability of the edge with a large weight value βij should also be larger. According to this principle, we set the infection rate between two adjacent nodes vi and vj as follows:
βij=wi,jwmax, | (4.5) |
where wmax is the maximum weight among all edges in the aviation network.
Ranking | DC | BC | CC | EC | LC |
1 | Beijing | Beijing | Beijing | Shanghai | Shanghai |
2 | Shanghai | Shanghai | Shanghai | Beijing | Beijing |
3 | Shenzhen | Chengdu | Shenzhen | Shenzhen | Shenzhen |
4 | Chengdu | Shenzhen | Guangzhou | Chengdu | Chengdu |
5 | Guangzhou | Guangzhou | Chengdu | Guangzhou | Guangzhou |
6 | Chongqing | Kunming | Chongqing | Chongqing | Chongqing |
7 | Kunming | Hangzhou | Hangzhou | Hangzhou | Hangzhou |
8 | Hangzhou | Hohhot | Kunming | Kunming | Kunming |
9 | Xiamen | Urumqi | Nanjing | Nanjing | Nanjing |
10 | Zhengzhou | Chongqing | Zhengzhou | Harbin | Harbin |
The recovery rate γ of the network is set to 1. The SIR propagation experiments based on each method are repeated 100 times to take the average value. The increase in the number of nodes that have been infected in the network as time increases is plotted in Figure 3. It can be observed that the number of infected nodes in the network increases rapidly with time during t<5, while it gradually stabilizes after t>5, in which the top 10 ranked nodes derived from the LC method as the initial infected nodes when propagation is stabilized can infect the network with the maximum number of nodes. This demonstrates that the top 10 ranked nodes derived from the LC method have the most influence and can maximize the influence of the nodes in the network. It can be seen that the top ten nodes of EC and LC are the same. Due to the issue of propagation probability in the propagation process, the SIR model has produced different results, but still ranks first and second. Protecting the most influential hub cities in the aviation network identified by the LC method can effectively stop the spread process of similar disease outbreaks.
In Figure 4, the X-axis represents the percentage of attacked nodes, that is, the ratio of the number of nodes that have been removed from the network to the initial number of nodes in the network, and the Y-axis represents the five evaluation methods for assessing the network stickiness metrics: The average degree <k>, the average clustering coefficients C, the ratio of maximally connected subgraphs S, the global efficiency Ew, and the average path length Lw. The networks all adopt two types of modes: random attacks and intentional attacks.
The trend of the average degree of the network <k> in the aviation network under random and intentional attack modes are shown in Figure 4a. The change of the average degree can reflect the robustness and anti-jamming ability of the network. As can be seen in Figure 4a: In the random attack mode, <k> shows a slow decreasing trend throughout the process. In the intentional attack mode, <k> shows an exponential decrease. It can be seen from the trend that the LC method of attack minimizes <k> the fastest compared to other methods.
The trends of the clustering coefficient C of the Chinese aviation network in random and intentional attack modes are shown in Figure 4b. It can be seen that in the random attack mode, the value of C changes more smoothly, indicating that the impact of random interference on the local transportation efficiency is not too great at the initial stage. Compared to random attack, in intentional attack mode, the value of C decreases sharply and collapses rapidly, and the air transportation industry enters a paralyzed state rapidly, indicating that intentional attack has a large impact on local transportation efficiency. It can be seen that the LC method is very close to DC and EC in attacking network effectiveness, but is generally superior to BC and CC.
The trends in the relative size of the maximum connected subgraph (S) in the Chinese aviation network under random and intentional attack modes are shown in Figure 4c. In Figure 4c, S exhibits a linear and continuously decreasing trend throughout the process in the random attack mode. It shows that the aviation is resistant to random attacks. In the intentional attack mode, S decreases dramatically. At %n=20, the decreasing trend slows down, at which the network also nearly collapses, and at %n=25, S is almost equal to 0. It can be seen from the changing trend that the LC method of attack minimizes S the fastest compared to other methods.
The trends of global efficiency Ew of in the Chinese aviation network under random and intentional attack modes are shown in Figure 4d. In the random attack mode, Ew shows a flat trend throughout, and the global efficiency does not decrease with the attack of the nodes, indicating that the random attack has little impact on the robustness of the network. While in the intentional attack mode, a variety of methods lead to a sharp drop in Ew. Between %n=8 and %n=19, the experimental difference between the various methods is not large, and Ew of the LC method drops to its lowest value at %n=19 then remains stable. The trend shows that the LC method can minimize the global efficiency of the network faster.
The trend of the average shortest path Lw is shown in Figure 4e. For the aviation network, the smaller Lw means fewer node cities need to transit in the air transportation process. In random attack mode, the change of the value of Lw is smooth with slight fluctuation. When %n=13.5, the value of Lw increases slightly, but is not very different from the initial value of the average shortest path of the network, and it has remained stable. These changes show that random attacks have little effect on the general convenience of the aviation network. In the case of an intentional attack, there is a sharp increase and then a sharp decrease, at %n=13.5. The fastest increase of Lw is caused by the attack of the LC method, which shows that the intentional attack has a significant impact on the overall convenience of the aviation network in the initial stage, because the node cities removed in the initial stage are the central cities of navigation. The intercity air transportation needs to make multiple transits, so the value of the node cities will increase dramatically as the nodes are removed. After two sharp changes, the remaining nodes of intercity transportation by the removed nodes of the impact of the nodes become smaller or even unaffected, that is, the number of transshipments will also be reduced to the average level, so after %n=23, the convenience of the entire network is gradually completely lost. The above changes indicate that intentional attacks have a large impact on the overall convenience of the aviation network. It can be seen that in the initial stage, the LC method breaks the network at the fastest speed, and in the later stage, the BC method can lead the other methods to make the average shortest path of the network maximize the words; the LC method is also superior to the other methods.
In this paper, we extended the Laplacian energy centrality to weighted networks and applied this method to a real dataset, a Chinese aviation network, which was transformed into a weighted undirected complex network. Based on the analysis of the topological characteristics of the aviation network, the central nodes of the network were identified using weighted Laplacian energy centrality. The results were validated from the perspectives of network dynamics and robustness.
The results of the topological analysis indicated that the weighted aviation network exhibits the characteristics of a small-world network with a low clustering coefficient. Experimental findings demonstrated that the weighted Laplacian energy centrality accurately identifies key hub nodes in the Chinese aviation network. In dynamic-based experiments, when the information propagation process simulated by SIR tends to stabilize, the results showed that weighted Laplacian energy centrality infects the most nodes at the end of propagation. It showed that weighted Laplacian energy centrality has the maximum impact on information propagation within the network compared to other methods. In robustness analysis experiments, we divided the attack modes into intentional attacks and random attacks. In attack scenarios, the results showed that the sequence of key nodes obtained by weighted Laplacian energy centrality can effectively and quickly destroy the network and weighted Laplacian energy centrality consistently outperforms other methods in terms of the average degree, clustering coefficient, global efficiency, maximum connected subgraph, and average shortest path of the aviation network. These results are valuable for identifying critical nodes in the aviation network, optimizing its layout, protecting crucial network nodes, and minimizing losses caused by unexpected events.
The authors declare they have not used Artificial Intelligence (AI) tools in the creation of this article.
The authors thank two anonymous referees for valuable comments which have considerably improved the presentation of this paper.
The authors declare no conflicts of interest.
[1] |
M. Ouyang, Z. Pan, L. Hong, L. Zhao, Correlation analysis of different vulnerability metrics on power grids, Physica A, 396 (2014), 204–211. https://doi.org/10.1016/j.physa.2013.10.041 doi: 10.1016/j.physa.2013.10.041
![]() |
[2] |
B. S. Kerner, Criticism of generally accepted fundamentals and methodologies of traffic and transportation theory: A brief review, Physica A, 391 (2013), 5261–5282. https://doi.org/10.1016/j.physa.2013.06.004 doi: 10.1016/j.physa.2013.06.004
![]() |
[3] |
D. Chen, H. Gao, L. Lü, T. Zhou, Identifying influential nodes in large-scale directed networks: The role of clustering, PLoS One, 8 (2013), e77455. https://doi.org/10.1371/journal.pone.0077455 doi: 10.1371/journal.pone.0077455
![]() |
[4] |
D. Chen, L. Lü, M. Shang, Y. Zhang, T. Zhou, Identifying influential nodes in complex networks, Physica A, 391 (2012), 1777–1887. https://doi.org/10.1016/j.physa.2011.09.017 doi: 10.1016/j.physa.2011.09.017
![]() |
[5] |
B. Michele, C. Davide, V. Simone, Efficiency of attack strategies on complex model and real-world networks, Physica A, 414 (2014), 174–180. https://doi.org/10.1016/j.physa.2014.06.079 doi: 10.1016/j.physa.2014.06.079
![]() |
[6] |
S. P. Borgatti, Identifying sets of key players in a social network, Comput. Math. Organ. Theory, 12 (2006), 21–34. https://doi.org/10.1007/s10588-006-7084-x doi: 10.1007/s10588-006-7084-x
![]() |
[7] |
T. Wen, D. Pelusi, Y. Deng, Vital spreaders identification in complex networks with multi-local dimension, Knowl. Based Syst., 195 (2020), 105717. https://doi.org/10.1016/j.knosys.2020.105717 doi: 10.1016/j.knosys.2020.105717
![]() |
[8] |
J. Zhao, Y. Song, F. Liu, Y. Deng, The identification of influential nodes based on structure similarity, Connect. Sci., 33 (2021), 201–218. https://doi.org/10.1080/09540091.2020.1806203 doi: 10.1080/09540091.2020.1806203
![]() |
[9] |
L. Freeman, Centrality in social networks conceptual clarification, Soc. Networks, 1 (1978), 215–239. https://doi.org/10.1016/0378-8733(78)90021-7 doi: 10.1016/0378-8733(78)90021-7
![]() |
[10] |
A. Barrat, M. Barthelemyt, R. Pastor-Satorrast, A. Vespignani, The architecture of complex weighted networks, PNAS, 101 (2004), 3747–3752. https://doi.org/10.1073/pnas.0400087101 doi: 10.1073/pnas.0400087101
![]() |
[11] | G. Sabidussi, The centrality index of a graph, Psychometrika, 31 (1966), 581–603. |
[12] |
L. Freeman, A set of measures of centrality based upon betweenness, Sociometry, 40 (1977), 35–41. https://doi.org/10.2307/3033543 doi: 10.2307/3033543
![]() |
[13] |
T. Opsahl, F. Agneessens, J. Skvoretzc, Node centrality in weighted networks: Generalizing degree and shortest paths, Soc. Networks, 32 (2010), 245–251. https://doi.org/10.1016/j.socnet.2010.03.006 doi: 10.1016/j.socnet.2010.03.006
![]() |
[14] |
Y. Ma, Z. Cao, X. Qi, Quasi-Laplacian centrality: A new vertex centrality measurement based on Quasi-Laplacian energy of networks, Physica A, 527 (2019), 121130. https://doi.org/10.1016/j.physa.2019.121130 doi: 10.1016/j.physa.2019.121130
![]() |
[15] |
S. Zhao, S. Sun, Identification of node centrality based on Laplacian energy of networks, Physica A, 609 (2023), 128353. https://doi.org/10.1016/j.physa.2022.128353 doi: 10.1016/j.physa.2022.128353
![]() |
[16] |
S. Bansal, J. Sen, Network assessment of Tier-Ⅱ Indian cities airports in terms of type, accessibility, and connectivity, Transp. Policy, 124 (2022), 221–232. https://doi.org/10.1016/j.tranpol.2021.05.009 doi: 10.1016/j.tranpol.2021.05.009
![]() |
[17] |
R. W. Daniel, R. Soumen, M. D. S. Raissa, Resilience and rewiring of the passenger airline networks in the United States, Phys. Rev. E, 82 (2010), 056101. https://doi.org/10.1103/PhysRevE.82.056101 doi: 10.1103/PhysRevE.82.056101
![]() |
[18] | A. Reggiani, S. Signoretti, P. Nijkamp, A. Cento, Network measures in civil air transport: A case study of lufthansa, 1 Eds., Berlin: Springer Press, 2009. https://doi.org/10.1007/978-3-540-68409-1-14 |
[19] |
J. Shen, H. Zong, Identification of critical transportation cities in the multimodal transportation network of China, Physica A, 628 (2023), 129174. https://doi.org/10.1016/j.physa.2023.129174 doi: 10.1016/j.physa.2023.129174
![]() |
[20] | F. Gao, Y. Dang, Analysis on distribution property of an international air transport network, Sci. Sci. Manage. S. T., 7 (2009), 75–79. |
[21] |
R. Guimera, L. Amaral, Modeling the world-wide airport network, Eur. Phys. J. B, 38 (2004), 381–385. https://doi.org/10.1140/epjb/e2004-00131-0 doi: 10.1140/epjb/e2004-00131-0
![]() |
[22] |
J. Hu, Y. Wang, X. He, Analysis and application of global aviation network based on complex network, Comput. Sci., 48 (2021), 321–325. https://doi.org/10.11896/jsjkx.200900112 doi: 10.11896/jsjkx.200900112
![]() |
[23] | W. Liu, M. Han, Z. Xie, Connectivity characteristics and community identification of worldcity network based on global airline, Econ. Geogr., 40 (2020), 34–40. |
[24] |
O. Lordan, J. Sallan, P Simo, Robustness of the air transport network, Transp. Res. Part E, 68 (2014), 155–163. https://doi.org/10.1016/j.tre.2014.05.011 doi: 10.1016/j.tre.2014.05.011
![]() |
[25] |
H. Mo, F. Jin, Y. Liu, J. Wang, Network analysis on centrality of airport system, Sci. Geol. Sin., 30 (2010), 204–212. https://doi.org/10.13249/j.cnki.sgs.2010.02.204 doi: 10.13249/j.cnki.sgs.2010.02.204
![]() |
[26] |
J. Li, X. Wen, M. Wu, F. Liu, S. Li, Identification of key nodes and vital edges in aviation network based on minimum connected dominating set, Physica A, 87 (2020), 123340. https://doi.org/10.1016/j.physa.2019.123340 doi: 10.1016/j.physa.2019.123340
![]() |
[27] | X. Luo, J. Wen, J. Zhong, Structural characteristics and robustness analysis of state-owned airline networks, Aeronaut. Comput. Tech., 51 (2021), 55–59. |
[28] |
D. J. Watts, S. H. Strogatz, Collective dynamics of 'small-world' networks, Nature, 393 (1998), 440–442. https://doi.org/10.1038/30918 doi: 10.1038/30918
![]() |
[29] | X. Feng, H. Jia, Aviation network robustness considering node failure and edge failure, J. Beijing Jiaotong Univ., 45 (2021), 84–92. |
[30] |
M. E. J. Newman, Y. Liu, Scientific collaboration networks. Ⅱ. Shortest paths, weighted networks, and centrality, Phys. Rev. E, 64 (2001), 016132. https://doi.org/10.1103/PhysRevE.64.016132 doi: 10.1103/PhysRevE.64.016132
![]() |
[31] |
P. Bonacich, Power and centrality: A family of measures, Am. J. Sociol., 92 (1987), 1170–1182. https://doi.org/10.1086/228631 doi: 10.1086/228631
![]() |
[32] |
K. Dietz, Infectious diseases of humans: Dynamics and control, Ann. Inter. Med., 117 (1992), 179. https://doi.org/10.1016/0169-4758(92)90018-W doi: 10.1016/0169-4758(92)90018-W
![]() |
[33] |
V. Latora, M. Marchiori, Efficient behavior of small-world networks, Phys. Rev. Lett., 87 (2001), 198701. https://doi.org/10.1103/PhysRevLett.87.198701 doi: 10.1103/PhysRevLett.87.198701
![]() |
[34] |
C. Castellano, R. Pastor-Satorras, Thresholds for epidemic spreading in networks, Phys. Rev. Lett., 105 (2010), 218701. https://doi.org/10.1103/PhysRevLett.105.218701 doi: 10.1103/PhysRevLett.105.218701
![]() |
1. | Haijiang Li, Xin Zhang, Peng Jia, Qianqi Ma, Research on the Pattern and Evolution Characteristics of Global Dry Bulk Shipping Network Driven by Big Data, 2025, 13, 2077-1312, 147, 10.3390/jmse13010147 |
Ranking | DC | BC | CC | EC | LC |
1 | Beijing | Beijing | Beijing | Shanghai | Shanghai |
2 | Shanghai | Shanghai | Shanghai | Beijing | Beijing |
3 | Shenzhen | Chengdu | Shenzhen | Shenzhen | Shenzhen |
4 | Chengdu | Shenzhen | Guangzhou | Chengdu | Chengdu |
5 | Guangzhou | Guangzhou | Chengdu | Guangzhou | Guangzhou |
6 | Chongqing | Kunming | Chongqing | Chongqing | Chongqing |
7 | Kunming | Hangzhou | Hangzhou | Hangzhou | Hangzhou |
8 | Hangzhou | Hohhot | Kunming | Kunming | Kunming |
9 | Xiamen | Urumqi | Nanjing | Nanjing | Nanjing |
10 | Zhengzhou | Chongqing | Zhengzhou | Harbin | Harbin |
Ranking | DC | BC | CC | EC | LC |
1 | Beijing | Beijing | Beijing | Shanghai | Shanghai |
2 | Shanghai | Shanghai | Shanghai | Beijing | Beijing |
3 | Shenzhen | Chengdu | Shenzhen | Shenzhen | Shenzhen |
4 | Chengdu | Shenzhen | Guangzhou | Chengdu | Chengdu |
5 | Guangzhou | Guangzhou | Chengdu | Guangzhou | Guangzhou |
6 | Chongqing | Kunming | Chongqing | Chongqing | Chongqing |
7 | Kunming | Hangzhou | Hangzhou | Hangzhou | Hangzhou |
8 | Hangzhou | Hohhot | Kunming | Kunming | Kunming |
9 | Xiamen | Urumqi | Nanjing | Nanjing | Nanjing |
10 | Zhengzhou | Chongqing | Zhengzhou | Harbin | Harbin |