Abstract

Since the “Belt and Road” Initiative, air transport has become the second largest mode of transport after ocean shipping and has contributed more and more to international trade between the countries. To explore the hierarchical characteristics of the aviation network of the “Belt and Road,” analyze the relationship between levels, and identify the core layer of the network, k-core analysis based on the “degree” value was designed and performed; the data of airports and air routes were collected to construct the “Belt and Road” aviation network model; then, k-core decomposition was conducted to reduce the size of nodes layer by layer from the outside to the inside, and a network structure model with 19 levels was obtained; the relationship between the coreness and centrality of nodes in the network was investigated, and the changes in the attribute of networks at different levels and the connection between networks were explored; according to the structural characteristics of the network, the “Belt and Road” aviation network was divided into three categories: the core layer, the middle layer, and the detail layer. The results showed that the highest-level network obtained by k-core decomposition was used as the core layer of the “Belt and Road” aviation network, including 53 airports in 37 countries, all of which have high degree centrality and eigenvector centrality; the core layer network exhibits strong “small world” characteristics; there is little difference between the degree values of nodes, and the network has high stability. The status of China’s airports in the core layer network is not prominent; the four regional sectors such as East Asia, Southeast Asia, Central Asia, and the Middle East are closely connected regarding air transport. Europe is relatively weakly connected with three sectors: East Asia, Southeast Asia, and Central Asia; geographical factors are still the dominant factors determining the status of hub airports along the “Belt and Road.”

1. Introduction

Since 2013, air transportation has become the second largest mode of transport after ocean shipping. Compared with ocean shipping, air transportation is less limited by geographical environment, more timely and more accessible. In addition, the influence of geopolitical and security situation is also significantly less than that of ocean shipping. With its unique technical and economic advantages, air transportation has maintained a rapid growth momentum in recent years. Taking China as an example, air transport is growing by more than 20% per year between China and the other countries of the “Belt and Road” before the global spread of COVID-19. In 2019, air transport accounted for 26 percent of China’s total imports and exports. Among 167 international routes opened in 2019, 105 or 62.87 percent was accessible to countries along the “Belt and Road.” The aviation network between countries along the “Belt and Road” is becoming more and more perfect, contributing more and more to international trade. However, the aviation network of the “Belt and Road” is a complex network consisting of more than 1,000 general airports and tens of thousands of air routes; a few key airports play a dominant role in the service function and service ability of the whole aviation network. Therefore, identifying the core layer that plays a leading role through the hierarchical analysis of the aviation network is the key to further optimize the aviation network structure and provide a stronger support for the “Belt and Road” Initiative.

At present, abundant results have been achieved in the research of the characteristics of aviation network structure. Empirical research on the network topology revealed that the international air transport network is “scale-free” with “small world” characteristics [1]; study on the degree centrality, closeness centrality, and eigenvector topological centrality of airports in the aviation network showed that North America and East Asia focus on the development of regional hub airports while Western Europe, Southeast Asia, and the Middle East focus on building their main airports into global hubs [2]; although China’s aviation network also has the characteristics of “small world” and “scale-free” [3], it is different from the global aviation network in which the difference in nodes is more prominent [4]. Besides, the aviation network of the “Belt and Road” has also become a research focus. Plenty of scholars have explored the spatial pattern and evolution of the aviation network of the “Belt and Road” through empirical research and pointed out the unbalanced overall spatial distribution of the aviation network of the “Belt and Road” [5, 6]. The heterogeneity of nodes is a remarkable feature of aviation network [7]. A large number of scholars adopted centrality analysis to study the differences between nodes in the network. They calculated the “degree centrality, betweenness centrality, and closeness centrality” of nodes to study the status of different nodes in the network and investigated their driving factors [810]. The abovementioned methods can find the airports playing a central role in the network but fail to identify the core layer of the aviation network, making the research results difficult to be used to further optimize the aviation network structure. Therefore, many scholars have tried to study the core groups in the network through the cohesive subgroups analysis [11, 12], analyzed the cliques in the network and the relationship within and between cliques by such methods as cliques, n-clique, k-cluster, and k-core [13], and explored the “substructures” of the network [14]; a series of algorithms were used to reveal how the structure of the network is composed of substructures [15] and how some important substructures exert direct impact on the overall function of the network [16]. Therefore, using k-core analysis based on the “degree” value to study the hierarchical structure of the aviation network and identifying the core layer of the network by analyzing the relationship between the layers can provide theoretical guidance for the further improvement of the “Belt and Road” aviation network.

2. The k-Core Analysis of the Network

2.1. The k-Core of the Network

Seidman proposed the k-core analysis method in 1983, believing that the minimum degree can be used to determine the groups with different cohesiveness in the network and realize the research on the hierarchical structure of the network [17]. k-core is a concept of subgroups based on degree value; if the network is abstracted as a graph G = (V, E) composed of a point set V and an edge set E, then k-core refers to a subgraph that satisfies the following condition; that is, each node in the subgraph is at least connected with k nodes in this subgraph [18]. The k-core of the network can be obtained by recursively removing all nodes with a degree value less than k and their edges until all nodes in the residual network have a degree value of at least k. The coreness of a node represents the deepest k-core that contains the node, and the maximum value of the coreness of nodes is defined as the coreness of the network. The set of nodes that exists in the k-core and not in the (k + 1)-core is defined as the “k-residual set” [19].

2.2. k-Core Decomposition

The k-core decomposition can reveal the hierarchical characteristics of the network [20]. K-core decomposition is realized by recursively removing all nodes with a degree value less than or equal to K in the network; since nodes with small degree values are generally located in the outer layer of the network and nodes with large degree values are often located in the inner layer of the network, the k-core decomposition is conducted layer by layer from the outside to the inside [21].

The process of k-core decomposition is described in Figure 1; in a specific network, the first degree value of each node in the network is calculated; then, according to the order of the degree values from small to large, the corresponding nodes and edges in the network are removed in turn. As shown in Figure 1, three nodes such as , and are isolated points with a degree value of 0. Therefore, these three nodes are removed firstly. The remaining nodes and their edges form a new subgraph, namely, 1-core; then, the nodes with a degree value of 1 in the network and their corresponding edges are removed () to obtain 2-core; the nodes with a degree value of 2 () and their corresponding edges in the network are removed to obtain 3-core; after removing the nodes with a degree value of 3 and their corresponding edges, all nodes and edges do not exist, and the decomposition process is terminated; the subgraph corresponding to the 3-core can be regarded as the core layer of the network.

In the process of k-core decomposition, the nodes are removed together with their corresponding edges. Therefore, although some nodes have different initial degree values during decomposition, they may eventually be kept at the same level. For example, the initial degree value of the node is 3, but when the node and its corresponding edge are removed, the degree value of is 2; thus, it is kept in the 2-core.

3. k-Core Analysis of the “Belt and Road” Aviation Network

3.1. Data Collection and Network Construction

In this paper, an aviation network model was established using the airports with international air routes in 65 countries along the “Belt and Road” as nodes and the air routes among these countries as edges. A total of 522 airports and 6,822 international air routes were collected. The data come from the website (http://www.variflight.com). With 522 airports as the rows and columns of the matrix, the 0-1 adjacency matrix A of the “Belt and Road” aviation network was constructed based on whether there are air routes between countries.

3.2. k-Core Decomposition of the Aviation Network of “Belt and Road”

The adjacency matrix A was input into the UCINET software for k-core analysis so that the coreness of airports in the aviation network can be calculated; the calculation shows that the maximum coreness is 19. As the coreness k increases, the extinction of nodes and edges in each layer of the network is shown in Table 1.

It can be seen from Table 1 that when k = 2 and k = 3, the change in node is the greatest, indicating that there are high-density subgroups in these two levels; when k = 4 and k = 5, the remaining nodes drop to 8%, and the subsequent sequence maintains a relatively stable state, suggesting that in the aviation network of “Belt and Road,” high-density subgroups of nodes are mainly distributed on the periphery of the network, which are manifested as the state-level and regional aviation network subsystems; when k > 6, the extinction of nodes is relatively stable, indicating that the high-density subgroups have disappeared. The subsequent network hierarchy is balanced, and the main structure of the network appears layer by layer. Generally, the extinction of edges is relatively stable, slightly increasing with the value of k, demonstrating that the peripheral nodes have a relatively small number of air routes.

3.3. Analysis of the Correlation between the Coreness and Centrality of the Nodes in “Belt and Road” Aviation Network
3.3.1. Centrality Indicators of Node

According to different standards, the centrality of nodes can be divided into four types: degree centrality, betweenness centrality, closeness centrality, and eigenvector centrality.(1)Degree Centrality. If a node is directly connected to many other nodes, then it is located at a relatively central position [22]. Relative degree centrality is shown in equation (1), where n is the number of nodes in the network and dx is the degree of node x. (2)Betweenness Centrality. If a node is on many lines, then it has the ability to control the connection of other nodes, thereby further affecting the function of the entire network [23]. Relative betweenness centrality is defined in the following equation:where represents the number of shortest paths between node i and j and is the number of shortest paths between node i and j passing through node d.(3)Closeness Centrality. The closer the node is to other nodes, the less it depends on other nodes, and the higher the closeness centrality the node has [24]. Relatively closeness centrality is expressed by equation (3), where dij is the number of upper edges of the shortest distance between nodes i and j.(4)Eigenvector Centrality. Eigenvector centrality needs to use “factor analysis” to determine the “dimensions” of the distance between nodes. The corresponding position of each node in each dimension is called an “eigenvalue,” and a series of such values are defined as eigenvectors [25]. For the network adjacency matrix, its element aij means the status contribution of node i to j. Let x represent the eigenvector centrality, and the expression is shown in equation (4). That is, the centrality of a node is a function of the centrality of other nodes connected to this node

3.3.2. The Correlation between the Coreness and Centrality of Nodes

Coreness and centrality are both important indicators to determine the core nodes of the aviation network. The adjacency matrix A of the aviation network of “Belt and Road” was input into the UCINET software to calculate the coreness and 4 centrality indexes of all 522 nodes, as shown in Figure 2.

With the change in the coreness of nodes, the closeness centrality and betweenness centrality indexes of the nodes are evenly distributed. The index of Pearson correlation between the coreness, closeness centrality, and betweenness centrality of nodes is 0.428 and 0.427, respectively, indicating that coreness has a weak correlation with closeness centrality and betweenness centrality; the Pearson correlation between the coreness, degree centrality, and eigenvector centrality is 0.783 and 0.861, respectively, indicating that coreness has a strong positive correlation with degree centrality and eigenvector centrality; as the coreness of nodes increases, both the degree centrality and the eigenvector centrality show an increasing trend; the larger the coreness of a node, the higher the probability of high degree centrality and eigenvector centrality of the node. A node with small coreness has a relatively low-degree centrality and eigenvector centrality.

4. Analysis of the k-Core Network Structure of the “Belt and Road” Aviation Network

4.1. k-Core Networking of the “Belt and Road” Aviation Network

The “Belt and Road” aviation network can be divided into 19 levels by k-core decomposition. The navigation nodes at each layer and the air routes between them constitute 19 networks of different sizes, denoted as ki networks, where i = 1, …, 19; the structure of some networks is shown in Figure 3.

4.2. Comparative Analysis of the Properties of ki Network

The average distance, clustering coefficient, and degree distribution of the network are three important indicators to reveal the properties of the network.

4.2.1. Average Distance

The distance of the network refers to the number of edges contained in the shortest path between any two nodes; the average distance is the mean of the distance between all pairs of nodes in the network, as shown in equation (5); the larger the value, the greater the network depth and the worse the convenience.where is the number of nodes and lij is the number of edges with the shortest distance between nodes.

4.2.2. Clustering Coefficient

If a node i has di edges connected to other nodes, the clustering coefficient Ci refers to the ratio of the number of edges that actually exist to the number of edges that may exist among di nodes, as shown in equation (6). The clustering coefficient C of the entire network is the mean of the clustering coefficient Ci of all nodes i. The clustering coefficient is an indicator to reflect the degree of aggregation among nodes in the network.

The 19 networks obtained from the k-core decomposition of the “Belt and Road” aviation network were input into the UCINET software to calculate the average distance and clustering coefficient of each network, and their values are displayed in Figure 4.

t can be seen from Figure 4 that the clustering coefficient of the network in each layer shows a gradual increase. The clustering coefficient of the k1 network is the smallest, at 0.417, and the clustering coefficient of the k19 network is the largest, reaching 0.640; in general, the clustering coefficient of the network in each layer is relatively small compared with the clustering coefficient of China’s aviation network and the global aviation network [26]; this reflects that the aggregation among the airports in the “Belt and Road” aviation network is relatively weak, and the status and role of hub airports along the “Belt and Road” are not significant. The average distance of the network in each layer exhibits an obvious downward trend. The average distance of the k1 network is the largest, at 2.796, which means two nodes in the k1 network can reach each other after an average of 2.796 transfers; the average distance of the k19 network is the smallest, at 1.505. The decrease in the average distance reflects that the degree of convenience of the network in each layer is continuously improved with the increasing degree of decomposition.

4.2.3. Degree Distribution

The degree distribution represents the probability that the degree value of a randomly selected node is exactly k. For the convenience of research, the cumulative degree distribution function is often used to reflect the degree distribution. It represents the probability distribution of nodes whose degree value is not less than k.

When the degree distribution conforms to the power-law distribution, the cumulative degree distribution function conforms to the power law with a power exponent of γ−1 [27], as shown in equation (7), where γ is also called the degree distribution index. The smaller the value of γ, the greater the difference between the degree value of the nodes in the network and the stronger the scale-free property of the network. In the log-log coordinate system, the power-law distribution corresponds to a straight line [28].

Seven representative networks after k-core analysis were selected to calculate the node degree value and degree distribution function and the cumulative degree distribution curves in the log-log coordinate system, as shown in Figure 5.

In the log-log coordinate system, the cumulative degree distribution curve of each network corresponds approximately to a straight line; thus, the above cumulative degree distribution curves conform to power-law distribution; fitting was performed on the above distribution curves, and the values of a and γ are shown in Table 2.

It can be seen from Table 2 that as the value of k increases, the degree distribution index γ shows an increasing trend, indicating that with the increasing degree of network decomposition, the difference between the degree values of nodes in the network is continuously reduced, and the scale-free attribute of the network is constantly weakened. The value of γ in k1 and k4 is less than 2, revealing that there is a considerable proportion of nodes functioning as hubs in the network, and the hub nodes have a great impact on the network function; the value of γ in the k7, k10, k13, and k16 networks is between 2 and 3, and the network structure is similar to Internet router network and Movie actors network [29], suggesting that the connection between nodes in the network is not strong, and there are a small number of hub nodes, but the impact of hub nodes on network functions is further decreased. The value of γ in the k19 network exceeds 3, which indicates that the network is similar to a random network, where the difference between the degree values of nodes in the network is not large. The realization of network functions does not depend on nodes with certain characteristics. Each node has a relatively limited influence on network functions, and the network is highly stable.

4.3. Analysis of the Connection of the Network in Each Layer

The analysis of the connection of the network in each layer can more in-depth reflect the law of connection between nodes and help to identify the core layer of the “Belt and Road” aviation network. Denote “k-residual set” as , i = 1, …, 19. The connection of the network in each layer can be expressed by the interconnection between “k-residual sets”; the percentage of the number of interconnected edges between “k-residual sets” in the total number of edges of the nodes in the “k-residual sets” is displayed in Table 3.

Table 3 shows that the proportion of connections between the “k-residual set” of each layer and is the maximum, reflecting that is at the core; since the maximum coreness of the network is 19, and k19 are the same network. Therefore, it can be seen that the k19 network has the highest percentage of connections with the networks in other layers.

5. Research on the Levels of the “Belt and Road” Aviation Network

5.1. Division of the “Belt and Road” Aviation Network

According to the structural parameters such as the number of nodes and edges in , the largest connected subgraph, clustering coefficient, and network density, the categories of the levels of the network can be divided [30]. In , since many nodes are isolated without connected to other nodes by edges, the nodes connected by edges in were regarded as “subgraphs,” and the subgraph with the largest number of nodes is called the largest connected subgraph. The calculated structural parameters of the network are shown in Table 4.

Based on the method for network level division proposed [31] by Xiao et al. combining the variation characteristics of the largest connected subgraph, clustering coefficient, and network density in different cores, the “Belt and Road” aviation network was divided into three categories: core layer network, middle layer network, and detail layer network.

5.1.1. Core Layer Network

Since and k19 are the same network, it can be regarded as the core network in the “Belt and Road” aviation network. Table 4 shows that this network has the largest clustering coefficient and network density, and 53 nodes in this network are connected to each other. Since coreness has a strong correlation with degree centrality and eigenvector centrality, the nodes of this layer have larger degree centrality and eigenvector centrality and play a greater role in connection. Since the network at this layer has a larger clustering coefficient and a smaller average distance, it exhibits stronger “small world” features and a higher degree of convenience; this network has the largest degree distribution index γ, indicating that there is little difference between node degree in the network, and the network has high stability. In addition, it can be seen from Table 3 that this network has the largest proportion of connections to networks at other layers, revealing that this network occupies a core position in the “Belt and Road” aviation network structure.

5.1.2. Middle Layer Network

The network composed of the nodes in , ,…, can be regarded as the middle layer network in the “Belt and Road” aviation network. Table 4 shows that in , , …, , the clustering coefficient and network density have increased significantly, and a larger connected subgraph appears in the network.

5.1.3. Detail Layer Network

It can be seen from Table 4 that in , ,…, , not only the clustering coefficient is 0 but the network density is also small, and there are few nodes in the largest connected subgraph, reflecting that most of the nodes are in an isolated state. Therefore, the nodes in , , …, constitute the detail layer network.

5.2. Analysis of the Core Layer Network in the “Belt and Road” Aviation Network
5.2.1. Composition of the Core Layer Network in the “Belt and Road” Aviation Network

The core network includes 53 airports in 37 countries; the country with the largest number of airports is the United Arab Emirates (5), followed by China and Saudi Arabia (4 each); the distribution is shown in Figure 6. The 37 countries covered by the core layer account for 56.92% of the total number of countries along the “Belt and Road,” 86.45% of the GDP, and 93.74% of the total population along the “Belt and Road.” They are highly representative.

5.2.2. Analysis of the Cohesive Subgroups of the Core Layer Aviation Network of “Belt and Road”

Cohesive subgroup analysis can reveal the actual or potential interrelationship between nodes in the network. Cliques are the most closely related cohesive subgroups (Ren et al. 2000); in the network structure consisting of 53 airports in the core layer, although some airports are not prominent from a global perspective, they may play a crucial role in the connection of the core layer; therefore, identifying cliques based on the strength of the connection between nodes is of great significance for effectively enhancing the overall connection of the aviation network.

UCINET software was used to analyze cliques in the core layer network. With a minimum of 3 nodes as the standard, a total of 857 cliques were identified, of which the largest clique contains 13 nodes and the smallest clique contains 4 nodes. Through the analysis of the identified cliques, the following characteristics of the nodes in the core layer network were discovered:Although there are many cliques identified, a large number of “overlapping” nodes exist; for example, Qatar’s Doha (QA-DOH) appears in 582 cliques and is the node overlapped most. In addition, airports that “overlap” frequently include Turkey’s Chistanbul (TR-IST), UAE’s Dubai (AE-DXB), and Thailand’s Bangkok (AE-DXB), which appear in 560, 374, and 371 cliques, respectively; this indicates that these airports play an important role in the connection of the core layer network, and they are also the most important hub airports in the “Belt and Road” aviation network.China’s 4 airports such as Beijing Capital (CN-PEK), Shanghai Pudong (CN-PVG), Guangzhou (CN-CAN), and Kunming (CN-KMG) appear in 183, 67, 41, and 28 cliques, respectively; except for Beijing, the number of cliques where the other three airports are engaged is much lower than the average value of 135, indicating that their status in the core layer is not prominent.The “Belt and Road” route can be divided into five major sectors: East Asia, Southeast Asia, Central Asia, the Middle East, and Europe. By analyzing the composition of airports in cliques, it can be seen that there exists a relatively close connection between the airports within each major regional sector, and a large number of clique members are from the same sector; the analysis of some large cliques showed that the four sectors such as East Asia, Southeast Asia, Central Asia, and the Middle East are closely related, and the airports in these four sectors often appear in the same clique. Europe only has a close connection with the Middle East sector, while its connection with the three sectors such as East Asia, Southeast Asia, and Central Asia is weak.Factors such as economy, population, and land do not have a significant impact on the status of airports in the network. Airports in China, India, Russia, and other countries do not have many cliques, while airports in Qatar, Turkey, UAE, and Thailand have a large number of cliques due to their unique geographical location, showing that geographical location is currently the main factor determining the status of the nodes in “Belt and Road” aviation network.

5.3. Analysis of the Middle Layer Network in the “Belt and Road” Aviation Network

The middle layer network includes 92 airports in 39 countries; the country with the largest number of airports is China (15), followed by Russia (9) and India (6 each); the distribution is shown in Figure 7.

The middle layer network contains 92 airports, among which the airports with high degree are Moscow Vnukovo International Airport (RU-VKO), Bangkok Don Mueang International Airport (TH-DMK), Indonesia Denpasar International Airport (ID-DPS), and Tan Son Nhat International Airport (VN-SGN) in Vietnam. Airports in the middle layer network have a high K-core value. These airports are mostly secondary aviation hubs in the country where they are located, and they have a certain degree of influence in the network. Compared with the airports in the core layer network, the airports in the middle layer network are connected to relatively single objects by air routes, and the clustering coefficient is not high, suggesting that these airports do not tend to form clusters.

5.4. Analysis of the Detail Layer Network in the“Belt and Road” Aviation Network

The detail layer network consists of 460 airports in 64 countries; the distribution is shown in Figure 8. Despite a large number of nodes, this network is the most sparse, with only 221 routes, accounting for 3.24% of the total number of routes of the “Belt and Road” aviation network. Among them, 258 airports have no connection with other airports at the same layer by air routes; this shows that they are on the outer edge of the “Belt and Road” aviation network, and there are a large number of isolated nodes or plenty of relatively independent subgroups formed by a large quantity of low-degree nodes.

The node airports in the detail layer network only exist in the “Belt and Road” aviation network by connecting with the airports in the core layer network or the middle layer network. Therefore, they are the radiation and expansion of the core layer network and the middle layer network to the periphery. Although the airports in the network at this layer play a small role, they effectively increase the transportation service coverage of the entire aviation network and are also important components of the “Belt and Road” aviation network.

6. Conclusion

In this paper, the data on airports and air routes were collected to construct the “Belt and Road” aviation network model; k-core decomposition was conducted, and a network structure model with 19 levels was obtained; the relationship between the coreness and centrality of nodes in the network was studied, the attribute variation of the network at different levels was analyzed, and the connection between these networks was explored; finally, based on the structural characteristics of the 19-layer network, the “Belt and Road” aviation network was divided into three categories: the core layer, the middle layer, and the detail layer networks, and their respective characteristics were investigated.

The following conclusions were obtained in this research:(1)The k-core analysis based on the “degree” value can distinguish different cohesive groups in the aviation network hierarchically, to realize a hierarchical network; through the k-core decomposition, the group size is reduced from the outside to the inside, and finally a relatively important core layer is obtained, which can effectively reduce the complexity of the aviation network.(2)The highest-level network obtained by k-core decomposition can be used as the core layer of the “Belt and Road” aviation network, which includes 53 general aviation airports in 37 countries and shows strong “small world” characteristics; the difference between the degree values of nodes is not large, suggesting that the network has high stability.(3)In the future, the connection between Europe and East Asia, Southeast Asia, and Central Asia should be strengthened, and the density of air routes should be increased; geographical factors are the dominant factors determining the status of the airport, indicating that air transport mainly plays a role in stimulating personnel exchange along the “Belt and Road,” and its role in economic and trade development needs to be further intensified.

Data Availability

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

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This study was supported by National Social Science Foundation of China (16BGL014).