Next Article in Journal
Possibility and Time in Quantum Mechanics
Previous Article in Journal
Double-Matrix Decomposition Image Steganography Scheme Based on Wavelet Transform with Multi-Region Coverage
Previous Article in Special Issue
Enhancement of Cooperation and Reentrant Phase of Prisoner’s Dilemma Game on Signed Networks
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Analysis of the Structure and Dynamics of European Flight Networks

by
Matteo Milazzo
1,
Federico Musciotto
2,
Salvatore Miccichè
2,* and
Rosario N. Mantegna
2,3
1
Dipartimento di Fisica e Astronomia Ettore Majorana, Università degli Studi di Catania, Via S. Sofia 64, I-95123 Catania, Italy
2
Dipartimento di Fisica e Chimica Emilio Segrè, Università degli Studi di Palermo, Viale delle Scienze, Ed. 18, I-90128 Palermo, Italy
3
Complexity Science Hub Vienna, Josefstaedter Strasse 39, 1080 Vienna, Austria
*
Author to whom correspondence should be addressed.
Entropy 2022, 24(2), 248; https://doi.org/10.3390/e24020248
Submission received: 25 December 2021 / Revised: 30 January 2022 / Accepted: 5 February 2022 / Published: 8 February 2022
(This article belongs to the Special Issue Structure and Dynamics of Complex Socioeconomic Networks)

Abstract

:
We analyze structure and dynamics of flight networks of 50 airlines active in the European airspace in 2017. Our analysis shows that the concentration of the degree of nodes of different flight networks of airlines is markedly heterogeneous among airlines reflecting heterogeneity of the airline business models. We obtain an unsupervised classification of airlines by performing a hierarchical clustering that uses a correlation coefficient computed between the average occurrence profiles of 4-motifs of airline networks as similarity measure. The hierarchical tree is highly informative with respect to properties of the different airlines (for example, the number of main hubs, airline participation to intercontinental flights, regional coverage, nature of commercial, cargo, leisure or rental airline). The 4-motif patterns are therefore distinctive of each airline and reflect information about the main determinants of different airlines. This information is different from what can be found looking at the overlap of directed links.

1. Introduction

The air transportation system (ATS) is a socio-technical system analyzed as a complex network for many years [1,2]. The ATS is analyzed at different geographical scales (see, for example, studies covering the ATSs of China [3], Europe [4] and the U.S. [5]) and at different resolutions starting from the airport–flight network down to the network of the reference points used in the definition of flight routes (called navigation points) [6].
In the majority of studies, the ATS is investigated by setting up a flight network where nodes are airports and flights connecting airports are links. The flight networks have been investigated by considering them undirected and/or directed networks (in this last case, the direction of the links originates from the departing airport and ends up in the arrival airport), unweighted and/or weighted [7]. Several studies have considered the problem of the resilience of the ATS to failures and attacks [5,8,9,10]. Other studies have selected a subset of links (labeled as the “backbone” of the ATS) presenting statistical properties that are not consistent with a specific null hypothesis [11,12], making the ATS one of the first systems where statistically validated networks [13] have been investigated.
The ATS is a complex system composed of well-defined subunits. In fact, flights are operated by different airlines that compete and collaborate among them. Since 2010, the ATS has been analyzed by distinguishing the role of its subunits (i.e., by analyzing properties of flight networks of single airlines [5]). Moreover, the presence of different flight networks observed for different airlines made this system a natural candidate for the study of so-called multiplex, which are networks where nodes can have multiple kinds of relations called layers. In fact, the ATS was one of the first socio-technical systems described as a multiplex, where layers represent flights operated by different airlines [14,15].
Flight networks have been investigated from different perspectives and at different scales [16,17,18], for example, by considering basic network metrics, topology of the degree distribution, resilience to attack or failures, community detection of large clusters and computation and analysis of network motifs. Motifs are isomorphic subnetworks of a specified number of nodes and shape. Motifs were first investigated in studies of social networks [19]. In these earlier studies, motifs were primarily investigated as triads (i.e., as subnetworks of three nodes) and put in relation with the properties of the degree sequence. At the beginning of this century, such structures were also investigated in biological systems under the name of motifs [20]. By considering isomorphic motifs (i.e., subnetworks where the identity of the node is not taken into account when considering the shape of the subnetwork), there are 13 isomorphic for subnetworks with 3 nodes or 3-motifs. This number soon explodes when the subnetwork includes more nodes. For subnetworks with 4 nodes (or 4-motifs), one counts 199 isomorphic motifs [20].
Network motifs have been investigated in flight networks both in studies comparing the informativeness of network detection in several types of complex networks [21] and in studies fully focused on the static and dynamics characteristics of the flight networks [22,23,24,25]. In this study, we investigate the temporal evolution of 3-motifs and 4-motifs for the 50 European airlines with the highest number of flights in the European Civil Aviation Conference (ECAC) airspace in the year 2017. By investigating the number and temporal evolution of the 3- and 4-motifs, we are able to perform an unsupervised classification of the 50 airlines indicating that main differences among different airlines are due to their regional specialization (including the ability to perform intercontinental flights) and to their business model. We observe that the business model of each airline ranges between the two stylized models of hub-and-spoke and point-to-point business models [26,27]. In a hub-and-spoke model, one or more airports act as “hubs”, i.e., as special airports directly connecting all remaining airports. In a hub-and-spoke structure with a single hub, the network therefore has a star topology with the hub at the center of the star and all the other airports acting as leaves of the network. In the point-to-point structure, all the airports are equivalent and the network degree is characterized by pair interconnections between airports.
The main goal of our investigation is a reliable and effective classification of airlines. The classification is obtained by an unsupervised methodology that only takes into account the information about the airline flights. We hypothesize that the business models of each airline induce specific constraints on its flight network. These constraints are reflected in the motif occurrence of each airline. Our network analysis shows that European airlines present a heterogeneous profile distributed between the two boundaries of hub-and-spoke and point-to-point business models. The heterogeneity is clearly shown by using a measure of concentration of degree in the degree sequence. Specifically, as a measure of concentration, we use an adapted version of the Herfindal–Hirshman index [28,29]. For the sake of simplicity, in the remaining text, we will call this index by the more traditional, although imprecise, name of Herfindal index.The time evolution of motifs shows that the basic temporal unit of the flight schedule is the week. Differences in the degree concentration observed during winter and summer schedules are detected, but their amount is negligible for most airlines. Average values of the motif occurrences may therefore be a useful proxy of the average behavior of the airlines over a calendar year. By using average values of the 4-motifs occurrence, we are able to obtain an unsupervised classification of airlines. The obtained hierarchical clustering is showing that the presence of a given number of hubs together with the presence or absence of intercontinental flights characterizes groups of airlines. On the other hand, a hierarchical clustering based on a similarity measure estimated starting from the co-presence of the two airlines in the origin–destination flight is providing a poorly informative hierarchical clustering.
The paper is organized as follows. In Section 2, we discuss the data used in our analysis and the metrics and methods used to characterize flight networks. In Section 3, we present our results about the heterogeneity of the degree concentration and our results about the structure and time evolution of 3- and 4-motifs for the different airlines. Average 4-motif occurrences are used to perform an unsupervised clustering of the 50 airlines providing an informative hierarchical cluster. In Section 4, we discuss our results.

2. Data and Methods

We investigate the flight networks of the 50 biggest commercial airlines flying over the European flight zone. Specifically, we consider all flights that occurred during the period from 1 January 2017 to 31 December 2017.
A flight network is a network where nodes are airports and links are flights that occurred in a given time interval. By considering that the flight occurs from a departing airport to an arrival airport, flight networks can be described as directed weighted networks (where the weight of a link is the number of flights that occurred from airport i to airport j in the chosen time interval). In this study, we considered flight networks as directed networks while we disregard the weights of the links. Networks are computed using daily and weekly time intervals.
Flight networks and their metrics of each airline are analyzed both in their time evolution and in their subunits. Specifically, we investigate the daily degree sequence of each airline for each day. In our analysis, we primarily focus on the concentration of the highest degree values on a limited set of airports usually described as “hubs”. This is performed by adapting the Herfindal index, i.e., a well-known measure of concentration, to the degree sequence. The subunits analysis is carried out by considering all isomorphic small networks with 3 or 4 nodes. These subnetworks are called motifs in the biological literature or triads or subnetworks in the social science literature.
We compare similarity between pairs of airlines both by considering the links, i.e., flights, they are performing on a specific day or week and by considering the motifs they present on a specific day or on average over the full year. Similarity between the airlines is therefore estimated and interpreted by extracting hierarchical trees from the selected similarity matrix.

2.1. Flight Data

Our dataset comprises all the flights that, even partly, cross the ECAC airspace for the entire 2017 year. Data were obtained by EUROCONTROL (http://www.eurocontrol.int, accessed on 4 February 2022), the European public institution that coordinates and plans air traffic control for all of Europe.
Specifically, we obtained access to the Demand Data Repository (DDR) from which one can obtain all flights followed by any aircraft in the ECAC airspace. Data about flights contain several types of information. In the present study, we just focus on the origin–destination of each flight crossing the ECAC airspace at a given time.
By considering that our focus is on the specific characteristics of airlines, in the present study, we investigate flights of the major 50 airlines performing flights in the ECAC airspace in 2017. In our set, we do not consider Air Berlin because this airline ceased operations on 27 October 2017. Since 2016, Germanwings has been a lease operator for its sister company Eurowings. In our set, we are not considering Germanwings flights. The selected airlines have performed 65.7 % of the total number of flights of 2017, which corresponds to approximately 3000 flights per company per month on average. The list of the 50 airlines is provided in Appendix A. The large majority of airlines are commercial airlines. There are 24 flag carrier airlines, 14 low cost carrier (LCC) airlines, 6 regional airlines, 2 leisure airlines, 2 scheduled airlines, 1 cargo airline and 1 rental airline.

2.2. Herfindal Index

The Herfindal index [28] has been introduced in the economic literature in order to measure the amount of competition among industrial firms. As such, it has also been used as an indicator of concentration, as large firms usually contribute more to the Herfindal index than smaller ones. In the context of complex networks, the Herfindal index can be defined as
H   =   i   =   1 N d i 2 m 2
where d i is the degree of node i and 2 m is twice the number of directed links.

2.3. Motifs Detection

The investigation of subnetworks of fixed size (also called motifs) has a long history. Originally investigated as triads and put in relation with the properties of the degree sequence in the investigation of social networks [19], they were then also introduced in biology where the term “motif” was used for the first time [20].
In network analysis, a motif of size k is a structure of k nodes not necessarily all linked between each other, as, for example, in Figure 1. Motifs are different from cliques. A clique is defined in undirected networks, and it is a subgraph such that every two distinct vertices are adjacent.
For size k = 3, there are 13 isomorphic 3-motifs. In Figure 1, we are showing all of them together with the classification scheme used in [20]:
Isomorphic 3-motifs present unidirectional links (as in the case of motifs labeled as 6, 12, 36, 38 and 98), bidirectional links (as in the case of motifs 78 and 238) and both types of links (as in the case of motifs 14, 46, 74, 102, 108 and 110).
The number of isomorphic 4-motifs is 199 and therefore much larger than 13. As for the 3-motifs, we use the classification of [20]. For the shape of each 4-motif, one can consult the motifs dictionary that can be downloaded from the website of Uri Alon laboratory.
Network motif analysis can be performed by computational or analytical approaches. In our investigation, we considered a computational approach as it allows for the exact count of network motifs. Computational approaches usually follow a three-step procedure that can be summarized as follows:
  • Search and enumerate occurrences of a topology with fixed size in the observed network;
  • Classify topologies by their isomorphic classes;
  • Calculate statistical significance for each isomorphic classes comparing occurrences with those in random ensemble.
In particular, we considered the mfinder [30] software developed by Uri Alon laboratory.

2.4. Average Linkage Clustering Analysis

We assess the similarity between each pair of the n airline by estimating the correlation between the average occurrence of each 4-motif of each airline. The average is computed over the 365 days of the year. To take into account the large interval of values observed for the different 4-motifs, we use the Spearman correlation coefficient. Therefore, by starting from the matrix of records obtained by averaging the occurrence of each 4-motif, we estimate a correlation matrix and we use the correlation ρ i j as a measure of similarity between airlines i and j.
From the correlation values, we compute a distance according to the relation d i j   =   2 ( 1     ρ i j ) . This distance is therefore used to extract a hierarchical tree with the method of the average linkage.
The average linkage cluster analysis is a hierarchical clustering procedure [31,32]. The procedure gives as an output a rooted tree or dendrogram. In this procedure, at each step, when two elements or one element and a cluster or two clusters p and q merge in a wider single cluster t, the distance d t r between the new cluster t and any cluster r is recursively determined as the average distance between any element of t and any other element of cluster r.

3. Results

3.1. Herfindal Index

Our first analysis determines the daily flight network of each investigated airline. The day is defined as the calendar day at European Central Time. For illustrative purposes, we show the networks of the nine biggest airlines on day 1 September 2017 in Figure 2. This day has been retrospectively selected as an example of a day with routinely operational activities.
For each flight network, we extract the degree sequence by considering the network as a directed network. The average values over the year of the number of nodes N (i.e., number of airports where airlines flight), the number of direct links E (i.e., the number of distinct origin destination flights), minimum degree, median degree, mean degree, maximum degree, standard deviation of the degree and Herfindal index are shown in Table 1. The metrics shown in Table 1 are quite basic and standard with the exception of the adaptation of the Herfindal index as an indicator of concentration in the degree sequence observed in one or more of the nodes.
Given the definition of Equation (1), a pure hub-and-spoke setting of flights would imply a Herfindal index of 0.25 for large values of N. This is what we observe (see Table 1) as average yearly value for Brussels Airlines (BEL), Aeroflot (AFL), KLM, Iberia Airlines (IBE) and Finnair (FIN). Networks of these airlines are very close to a pure hub-and-spoke setting. Other airlines show lower values of the average Herfindal index. The values observed range from 0.213 for Austrian Airlines to 0.016 for Ryanair, showing a high variability of the underlying flight network structure. For the sake of compactness, in Table 1, we show only the yearly average values of the selected indicators. To assess the degree of variability of the Herfindal index, we show in Figure 3 the daily profile of this index for the top ten airlines in number of flights. They are Ryanair (RYR), Lufthansa (DLH), Turkish Airlines (THY), EasyJet (EZY) Air France (AFR), Scandinavian Airlines (SAS), British Airways (BAW), KLM (KLM), Vueling Airlines (VLG) and Alitalia (AZA). Time dynamics of the Herfindal index is detectable for several airlines but fluctuations are quite limited and primarily reflect a weekly or intra-weekly periodicity. Some airlines also show detectable winter–summer dynamics. Examples are THY, AZA and VLG. Horizontal dashed line is the expected values of the Herfindal index for networks with only bidirectional links and with K pure hubs and all the remaining (large) number of leaves only flying to a single hub for K ranging from one (top dashed line) to five (bottom dashed line). In particular, KLM networks are compatible with a network structure having a single hub (i.e., Schipol airport), Lufthansa (DLH) networks are compatible with a two hub network (prominent Lufthansa hubs are Frankfurt and Munich airports). Vueling (VLG) and Scandinavian Airlines (SAS) have a pattern compatible with three or more hubs and/or with a prominent section of the flight network based on point-to-point flight circulation, whereas the Herfindal index of Ryanair and EasyJet have pretty low values, manifesting the poor relevance of the hub-and-spoke structure in their flight networks.
Our analysis can therefore confirm that flight network characteristics are deeply related to the business organization of each airline with a prominent role played by the choice of a hub-and-spoke versus a point-to-point structure and with a role played by the number of hubs characterizing the flight network.
In the next section, we investigate 3-motifs to better characterize similarity and differences among the flight networks of airlines.

3.2. 3-Motifs

We have computed the number of 3-motifs present on daily flight networks for all 50 airlines. In Figure 4, we show a color code map of the occurrence of the 13 isomorphic 3-motifs for the 9 largest airlines.
The occurrence of each 3-motif presents large variability among the different types of motifs and is correlated with properties of the flight networks such as number of nodes, number of links, number of bidirectional links and topology structure of the network. The most common 3-motif is motif 78. This type of motif is clearly manifesting that a hub-like structure and bidirectional links are essential ingredients of all flight networks. The other 3-motif with all bidirectional link, i.e., 3-motif 238, is significantly present in airlines presenting flight networks with a pronounced point-to-point structure, such as Ryanair, EasyJet and Vueling, or airlines having more than a single hub, such as Lufthansa, Turkish and Scandinavian Airlines. The 3-motifs with only unidirectional links are poorly observed (see average occurrence values of 3-motifs 6, 12, 36, 38 and 98). Some of the 3-motifs with mixed types of links are significantly present (for example, 3-motifs 14 and 74), while others are rather poorly expressed (as in the case of 3-motifs 102 and 108).
The profile of occurrence of the 3-motifs in different airlines is certainly informative. However, the number of 3-motifs is somewhat limited and therefore it is useful to consider motifs of larger size. In the next section, we investigate the occurrence of 4-motifs.

3.3. 4-Motifs

3.3.1. Daily Occurrence of 4-Motifs

We compute the occurrence of all 4-motifs for the daily flight networks of the 50 biggest airlines. In Figure 5, we show a color code map of the occurrence of the 199 isomorphic 4-motifs for the 9 largest airlines.
The profile of 4-motifs is richer than the one of the 3-motifs. Occurrences of the 4-motifs span about 5 orders of magnitude. For this reason, in Figure 5, we show the decimal logarithm to provide a comprehensive overview of the results. Airlines characterized by the presence of a single hub such as KLM present only a very limited number of 4-motifs with occurrence different from zero. Airlines with a business model closer to a point-to-point structure such as Ryanair and EasyJet present a higher number of observed 4-motifs. The other airlines characterized by a different number of hubs present an intermediate behavior between the two extremes. In addition to the presence or absence of a given motif at a given day, Figure 5 also shows a time variation of the occurrence of a given motif. To investigate the main frequencies associated with this time variation, we compute the periodogram of the occurrence of a set of 4-motifs. Specifically, we consider the twelve 4-motifs with the highest occurrence averaged over all considered days. In Figure 6, we show the power spectrum of the time evolution of the occurrence of the top twelve 4-motifs of Ryanair. For all 4-motifs, frequency peaks are detected for f = 0.14 day 1 and for its second and third harmonics. The main frequency f = 0.14 day 1 corresponds to a weekly cycle and the second and third harmonics correspond to a bi-weekly or three-weekly cycle. Therefore, the main underlying periodicity is the week periodicity as already observed in the estimation of the Herfindal index (see periodicity observed in Figure 3).
Tracking in details of the occurrence for the 199 different 4-motifs is impractical. For this reason, we first consider the ten 4-motifs with highest occurrence in the nine biggest airlines. Specifically, we rank these 4-motifs, and the rank is obtained by considering the daily occurrence of each motif averaged over all days of the year. Labels of these motifs are listed on Table 2 according to their rank for the 9 biggest airlines.
The link configuration of these 4-motifs is shown in Figure 7. In the figure, we show on the left 4-motifs composed by unidirectional links, whereas on the right we have motifs with only bidirectional links. The 4-motifs with both unidirectional and bidirectional links are shown in the middle of the figure. It is worth noting that 4-motifs with only unidirectional links (i.e., 4-motifs 14, 28, 280 and 2184) are only observed for KLM in the top 10. KLM is one of the airlines with an almost pure hub-and-spoke structure and the flight concurring to this type of motif is, in the majority of the cases, an intercontinental flight. All the other airlines have the top 10 4-motifs presenting a high number of bidirectional links. The 4-motif with the highest occurrence for all the top 9 airlines (and indeed the top 4-motif for 49 of 50 airlines) is 4-motif 4382. These motifs present three bidirectional links originating from the same node. As for 3-motif 78, the largest occurrence of this motif reflects the fact that at least one important airport is used as a hub by the airline generating the network. The other 4-motifs composed by only bidirectional links (i.e., 4-motifs 4698, 4958, 13260 and 13278) are compatible with a point-to-point structure or with a hub-and-spoke structure in the presence of at least two hubs. In fact, these 4-motifs are not observed for KLM and are observed at the highest rank for more oriented point-to-point airlines such as Ryanair, EasyJet and Vueling. They are also present when more than one hub is present as, for example, in the case of Lufthansa or Scandinavian Airlines. The ranking of the 4-motifs can therefore be used to evaluate the similarity of flight airline networks and we investigate this possibility in the next section.

3.3.2. Similarity of 4-Motif Profile

We use the information about the 4-motifs occurrence to obtain a categorization of airlines by using the methodology of Section 2.4. It is worth recalling here that, given this specific purpose, it is not necessary for us to maintain the information about the specific airports that is present in a motif. In fact, since we are interested in extracting a clusterization of airlines by using the structural information about the 4-motifs, only the isomorphic motifs will be relevant for us.
The result of our analysis is shown in Figure 8.
The hierarchical tree of Figure 8 is highly informative with respect to the clustering of groups of airlines. One airline markedly distinct from all others is NetJets Transportes Aéreos, S.A. (NJE). This airline is the only airline of the set providing rental of jets and therefore observing it distinct from all the others indicates that the observed flights of this rental company have 4-motifs that are quite distinct from the ones of all other airlines. An inspection of the hierarchical tree indicates the presence of clusters of airlines presenting a certain similarity among them and a degree of dissimilarity from the other airlines. Here, we wish to comment about some of them. One cluster is the cluster of KLM, Aeroflot (AFL) and Brussels Airlines (BEL). These three airlines are airlines with a single large hub as testified by a Herfindal index very close to 0.25 (see Table 1). Another cluster comprises Delta Air Lines (DAL), American Airlines (AAL) and United Airlines (UAL). These three airlines are American airlines primarily performing intercontinental flights. A large cluster is composed by Vueling Airlines (VLG), Volotea (VOE), Norwegian Air International (IBK), EasyJet (EZY), Ryanair (RYR), Eurowings (EWG), Wizz Air (WZZ), Norwegian Air Shuttle (NAX) and Scandinavian Airlines (SAS). These are all airlines with several hubs and/or with a point-to-point business model. Another cluster comprises Royal Air Maroc (RAM), Pegasus (PGT) and Turkish Airlines (THY). These airlines are primarily serving Middle East destinations and airlines are headquartered in Middle East countries. Another distinct cluster comprises Qatar Airways Company Q.C.S.C. (QTR), Austrian Airlines (AUA), European Air Transport Leipzig (BCS) and Lufthansa (DLH). With the exception of Qatar Airways, the airlines of this cluster are all based in central Europe. In fact, Lufthansa and European Air Transport Leipzig are German airlines (Lufthansa is the second largest commercial airline in Europe and European Air Transport Leipzig is the largest cargo company in Europe by number of flights) and Austrian Airlines is a subsidiary of the Lufthansa Group.

3.4. Airline Networks Overlap

It is worth estimating whether similarity between 4-motif occurrences could just be due to overlap between the links of the flight network of airlines. We rule out, to a large extent, this possibility by investigating the degree of overlap between all pairs of airline networks. Our investigation is conducted by estimating the Jaccard measure J ( n e t 1 , n e t 2 ) between each pair of airline networks n e t 1 and n e t 2 . The Jaccard similarity is defined as
J ( n e t 1 , n e t 2 )   =   | E 1 E 2 | | E 1 E 2 | ,
where | E 1 E 2 | is the number of directed links appearing in both flight networks, and | E 1 E 2 | is the number of links that appear in at least one of the two networks.
To take into account weekly variability of flight schedules, we have performed this analysis by considering the weekly schedule of each airline. The results obtained at the daily level are showing a degree of similarity of the same order or less. In Figure 9, we show the average linkage hierarchical tree obtained by using the Jaccard measure as a similarity measure. The hierarchical tree is poorly informative and only a very limited number of small clusters can be highlighted. This is in marked contrast with what we have obtained in the previous section when the similarity measure between airlines was obtained from the analysis of 4-motifs. The hierarchical tree shown in Figure 9 is representative of hierarchical trees obtained for all weeks of 2017.
By summarizing, we are the first to use the Herfindal index to characterize each airline operating in a given period. Moreover, by using the Herfindal index together with 3- and 4-motif analysis, we are able to achieve an unsupervised classification of airlines, clarifying the main characteristics of each airline.

4. Discussion and Conclusions

In the present study, we have analyzed the structure and dynamics of flight networks of 50 airlines performing most of the flights that occurred in the European airspace in 2017. Our analysis of directed flight networks shows that the degree concentration of the different networks is quite heterogeneous among the different airlines. We have been able to quantify this heterogeneity by using an adapted version of a classic measure of concentration, i.e., the Herfindal index. The Herfindal index provides a simple and reliable estimation of the closeness of the airline network to the reference models classified as hub-and-spoke and point-to-point. It can also be informative about the number of main hubs that are present in a network with a hub-and-spoke structure and multiple hubs. It is worth noting that the European ATS presents a very heterogeneous set of airline companies. In other words, business optimization performed at the level of a single airline generates different business models that eventually coexist in the global system.
The time evolution of the different networks presents a basic time cycle that is a weekly cycle. This basic timescale is evident both from the analysis of the time evolution of the Herfindal index and from the analysis of the time evolution of the occurrence of 3-motifs and 4-motifs. The summer–winter cycle primarily detected in the number of flights occurring daily or weekly does not significantly affect the long-term time evolution of the Herfindal index and the occurrence of 3-motifs and 4-motifs. These indicators are therefore more related to the type of business model followed by the airline than to the specific origin–destination links or number of flights operated in a given time interval.
In summary, an unsupervised classification based on hierarchical clustering and obtained by using a correlation coefficient between the occurrence profile of 4-motifs of airline networks as similarity measure is highly informative with respect to the properties of the different airlines (for example, the number of main hubs, their participation to intercontinental flights, their regional coverage, their nature of commercial, cargo, leisure or rental airline). The 4-motifs are therefore distinctive of the airlines and reflect information about the main determinants of the different airlines. Information is distinct from that originating from the overlap of the same directed links.
Such results indicate that a reliable and effective classification of airlines can be obtained by an unsupervised methodology that only takes into account the information about the airline flights. This is an important result given that, currently, the characterization of airlines and their business model has become a fundamental part of modern air transportation systems. An appropriate airline categorization is important not only for the practitioners but because it also influences the passengers’ perception. An indubitable advantage of our approach is that it is flexible as it may directly reflect any positioning of an airline within the general landscape of airlines, due to any change in its business model as reflected within its flight plans.

Author Contributions

Conceptualization, M.M., S.M. and R.N.M.; methodology, M.M., F.M., S.M. and R.N.M.; formal analysis, M.M. and S.M. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Italian Ministry of Research, grant number 2017WZFTZP.

Data Availability Statement

Restrictions apply to the availability of these data. Data were obtained from EUROCONTROL.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

The table below contains the list of airlines considered in this work and some information about their business model. Such data were obtained from the following links (accessed on 4 February 2022):
Table A1. Fifty airlines performing the largest number of flights in the ECAC airspace in 2017. Airlines are ordered by the total number of flights performed in 2017. The last column gives the number of flights performed by the airline in 2017.
Table A1. Fifty airlines performing the largest number of flights in the ECAC airspace in 2017. Airlines are ordered by the total number of flights performed in 2017. The last column gives the number of flights performed by the airline in 2017.
AirlineICAO CodeType# of Flights
RyanairRYRLCC815,878
LufthansaDLHFlag561,437
Turkish AirlinesTHYFlag555,720
EasyJetEZYLCC538,533
Air FranceAFRFlag369,991
Scandinavian AirlinesSASFlag330,515
British AirwaysBAWFlag304,112
KLMKLMFlag274,822
Vueling AirlinesVLGLCC219,839
AlitaliaAZAFlag215,033
Wizz AirWZZLCC186,713
PegasusPGTLCC179,942
FlybeBEERegional169,299
SwissSWRFlag156,173
TAP PortugalTAPFlag148,363
Austrian AirlinesAUAFlag141,305
Norwegian Air ShuttleNAXLCC138,010
EurowingsEWGLCC133,738
FinnairFINFlag128,253
WideroeWIFRegional127,047
AeroflotAFLFlag122,480
LOT Polish AirlinesLOTFlag117,962
Iberia AirlinesIBEFlag105,055
Norwegian Air InternationalIBKLCC104,835
Qatar Airways Company Q.C.S.C.QTRFlag99,298
Air EuropaAEAScheduled92,521
Brussels AirlinesBELFlag92,279
EmiratesUAEFlag90,894
Royal Air MarocRAMFlag90,626
Aer LingusEINFlag87,330
Air NostrumANERegional82,254
Hop!HOPRegional79,385
European Air Transport LeipzigBCSCargo76,697
TUI AirwaysTOMLeisure74,834
Jet2.comEXSLCC69,338
Transavia HollandTRALCC68,758
United AirlinesUALFlag67,796
Ukraine InternationalAUIFlag64,963
Delta Air LinesDALLCC64,552
Aegean AirlinesAEEFlag59,762
Olympic Air S.A.OALRegional58,135
SunExpress TKSXSLCC55,761
EasyJet Switzerland SAEZSLCC54,120
Air BalticBTIFlag52,722
American AirlinesAALFlag52,118
VoloteaVOELCC49,722
NetJets Transportes Aereos, S.A.NJERental49,338
Condor FlugdienstCFGLeisure46,257
LogonairLOGRegional45,688
Thomas Cook Airlines LimitedTCXScheduled45,230

References

  1. Barrat, A.; Barthelemy, M.; Pastor-Satorras, R.; Vespignani, A. The architecture of complex weighted networks. Proc. Natl. Acad. Sci. USA 2004, 101, 3747–3752. [Google Scholar] [CrossRef] [Green Version]
  2. Guimera, R.; Mossa, S.; Turtschi, A.; Amaral, L.A.N. The worldwide air transportation network: Anomalous centrality, community structure, and cities’ global roles. Proc. Natl. Acad. Sci. USA 2005, 102, 7794–7799. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Li, W.; Cai, X. Statistical analysis of airport network of China. Phys. Rev. E 2004, 69, 046106. [Google Scholar] [CrossRef] [Green Version]
  4. Burghouwt, G.; de Wit, J. Temporal configurations of European airline networks. J. Air Transp. Manag. 2005, 11, 185–198. [Google Scholar] [CrossRef]
  5. Wuellner, D.R.; Roy, S.; D’Souza, R.M. Resilience and rewiring of the passenger airline networks in the United States. Phys. Rev. E 2010, 82, 056101. [Google Scholar] [CrossRef] [Green Version]
  6. Gurtner, G.; Vitali, S.; Cipolla, M.; Lillo, F.; Mantegna, R.N.; Miccichè, S.; Pozzi, S. Multi-scale analysis of the European airspace using network community detection. PLoS ONE 2014, 9, e94414. [Google Scholar] [CrossRef] [Green Version]
  7. Zanin, M.; Lillo, F. Modelling the air transport with complex networks: A short review. Eur. Phys. J. Spec. Top. 2013, 215, 5–21. [Google Scholar] [CrossRef] [Green Version]
  8. Fleurquin, P.; Ramasco, J.J.; Eguiluz, V.M. Systemic delay propagation in the US airport network. Sci. Rep. 2013, 3, 1159. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  9. Sun, X.; Wandelt, S. Complementary strengths of airlines under network disruptions. Saf. Sci. 2018, 103, 76–87. [Google Scholar] [CrossRef]
  10. Oliveira, I.M.; Carpi, L.C.; Atman, A.P.F. The Multiplex Efficiency Index: Unveiling the Brazilian air transportation multiplex network—BATMN. Sci. Rep. 2020, 10, 13339. [Google Scholar] [CrossRef]
  11. Serrano, M.Á.; Boguná, M.; Vespignani, A. Extracting the multiscale backbone of complex weighted networks. Proc. Natl. Acad. Sci. USA 2009, 106, 6483–6488. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Marcaccioli, R.; Livan, G. A Pólya urn approach to information filtering in complex networks. Nat. Commun. 2019, 10, 745. [Google Scholar] [CrossRef] [PubMed]
  13. Tumminello, M.; Miccichè, S.; Lillo, F.; Piilo, J.; Mantegna, R.N. Statistically validated networks in bipartite complex systems. PLoS ONE 2011, 6, e17994. [Google Scholar] [CrossRef] [Green Version]
  14. Cardillo, A.; Zanin, M.; Gómez-Gardenes, J.; Romance, M.; del Amo, A.J.G.; Boccaletti, S. Modeling the multi-layer nature of the European Air Transport Network: Resilience and passengers re-scheduling under random failures. Eur. Phys. J. Spec. Top. 2013, 215, 23–33. [Google Scholar] [CrossRef] [Green Version]
  15. Cardillo, A.; Gómez-Gardenes, J.; Zanin, M.; Romance, M.; Papo, D.; Del Pozo, F.; Boccaletti, S. Emergence of network features from multiplexity. Sci. Rep. 2013, 3, 1344. [Google Scholar] [CrossRef] [Green Version]
  16. Zhang, S.; Derudder, B.; Witlox, F. Dynamics in the European Air Transport Network, 2003–9: An Explanatory Framework Drawing on Stochastic Actor-Based Modeling. Networks Spat. Econ. 2016, 19, 643–663. [Google Scholar] [CrossRef] [Green Version]
  17. Zhang, S.; Hu, Y.; Tang, X.; Fuellhart, K.; Dai, L.; Witlox, F. Modeling the Evolutionary Mechanism of China’s Domestic Air Transport Network. Sustainability 2020, 12, 6295. [Google Scholar] [CrossRef]
  18. Zhang, S.; Derudder, B.; Witlox, F. The determinants of full-service carriers airfares in European hub-to-hub markets. Eur. J. Transp. Infrastruct. Res. 2014, 14, 4. [Google Scholar]
  19. Wasserman, S.; Faust, K. Social Network Analysis: Methods and Applications; Cambridge University Press: Cambridge, UK, 1999; pp. 1–819. [Google Scholar]
  20. Milo, R.; Shen-Orr, S.; Itzkovitz, S.; Kashtan, N.; Chklovskii, D.; Alon, U. Network motifs: Simple building blocks of complex networks. Science 2002, 298, 824–827. [Google Scholar] [CrossRef] [Green Version]
  21. Tran, N.T.L.; DeLuccia, L.; McDonald, A.F.; Huang, C.H. Cross-disciplinary detection and analysis of network motifs. Bioinform. Biol. Insights 2015, 9, BBI-S23619. [Google Scholar] [CrossRef] [PubMed]
  22. Jin, Y.; Wei, Y.; Xiu, C.; Song, W.; Yang, K. Study on structural characteristics of China’s passenger airline network based on network motifs analysis. Sustainability 2019, 11, 2484. [Google Scholar] [CrossRef] [Green Version]
  23. Agasse-Duval, M.; Lawford, S. Subgraphs and motifs in a dynamic airline network. arXiv 2018, arXiv:1807.02585. [Google Scholar]
  24. Wei, Y.; Jin, Y.; Ma, D.; Xiu, C. Impact of colored motif characteristics on the survivability of passenger airline networks in China. Phys. Stat. Mech. Appl. 2021, 566, 125658. [Google Scholar] [CrossRef]
  25. Yang, H.; Le, M.; Wang, D. Airline Network Structure: Motifs and Airports’ Role in Cliques. Sustainability 2021, 13, 9573. [Google Scholar] [CrossRef]
  26. Fotheringham, A.S.; O’Kelly, M.E. Spatial Interaction Models: Formulation and Applications; Kluwer: Normwel, MA, USA, 1989. [Google Scholar]
  27. Derudder, B.; Devriendt, L.; Witlox, F. Flying where you don’t want to go: An empirical analysis of hubs in the global airline network. Tijdschr. Voor Econ. Soc. Geogr. 2007, 98, 307–324. [Google Scholar] [CrossRef]
  28. Hirschman, A.O. The Paternity of an Index. Am. Econ. Rev. 1964, 54, 761–762. [Google Scholar]
  29. Hall, M.; Tideman, N. Measures of concentration. J. Am. Stat. Assoc. 1967, 62, 162–168. [Google Scholar] [CrossRef]
  30. Mfinder Tool Guide. Available online: https://www.weizmann.ac.il/mcb/UriAlon/sites/mcb.UriAlon/files/uploads/NetworkMotifsSW/mfinder/mfindermanual.pdf (accessed on 4 February 2022).
  31. Anderberg, M.R. Cluster Analysis for Applications; Academic Press: New York, NY, USA, 1973. [Google Scholar]
  32. Coronnello, C.; Tumminello, M.; Lillo, F.; Miccichè, S.; Mantegna, R.N. Sector identification in a set of stock return time series traded at the London Stock Exchange. Acta Phys. Pol. B 2005, 36, 2653–2679. [Google Scholar]
Figure 1. Isomorphic 3-motifs and related codes according to the classification given in [20]. The arrow indicates a flight from origin (tail) to destination (head). Bidirectional flights are indicated with a double head arrow.
Figure 1. Isomorphic 3-motifs and related codes according to the classification given in [20]. The arrow indicates a flight from origin (tail) to destination (head). Bidirectional flights are indicated with a double head arrow.
Entropy 24 00248 g001
Figure 2. Airline networks on day 1 September 2017. Each panel refers to one of the nine biggest airlines. The airline name is indicated on top of the panel.
Figure 2. Airline networks on day 1 September 2017. Each panel refers to one of the nine biggest airlines. The airline name is indicated on top of the panel.
Entropy 24 00248 g002
Figure 3. Herfindal index of flight networks (defined as in Equation (1)) as a function of the calendar day for the top 10 airlines Ryanair (RYR), Lufthansa (DLH), Turkish Airlines (THY), EasyJet (EZY) Air France (AFR), Scandinavian Airlines (SAS), British Airways (BAW), KLM (KLM), Vueling Airlines (VLG) and Alitalia (AZA).
Figure 3. Herfindal index of flight networks (defined as in Equation (1)) as a function of the calendar day for the top 10 airlines Ryanair (RYR), Lufthansa (DLH), Turkish Airlines (THY), EasyJet (EZY) Air France (AFR), Scandinavian Airlines (SAS), British Airways (BAW), KLM (KLM), Vueling Airlines (VLG) and Alitalia (AZA).
Entropy 24 00248 g003
Figure 4. Color code representation of the base 10 logarithm of the daily occurrence of each isomorphic 3-motif (x-axis) as a function of the time (y-axis). The 3-motifs are indicated with the labels defined in [20]. A zero occurrence is indicated with a −1 value. Each panel refers to one of the nine biggest airlines. The airline name is indicated on top of the panel.
Figure 4. Color code representation of the base 10 logarithm of the daily occurrence of each isomorphic 3-motif (x-axis) as a function of the time (y-axis). The 3-motifs are indicated with the labels defined in [20]. A zero occurrence is indicated with a −1 value. Each panel refers to one of the nine biggest airlines. The airline name is indicated on top of the panel.
Entropy 24 00248 g004
Figure 5. Color code representation of the base 10 logarithm of the daily occurrence of each isomorphic 4-motif (x-axis) as a function of the time (y-axis). The 4-motifs are ranked according to the labels defined in [20]. Some of the most populated 4-motifs are indicated in the x-axis. A zero occurrence is indicated with a −1 value. Each panel refers to one of the nine biggest airlines. The airline name is indicated on top of the panel.
Figure 5. Color code representation of the base 10 logarithm of the daily occurrence of each isomorphic 4-motif (x-axis) as a function of the time (y-axis). The 4-motifs are ranked according to the labels defined in [20]. Some of the most populated 4-motifs are indicated in the x-axis. A zero occurrence is indicated with a −1 value. Each panel refers to one of the nine biggest airlines. The airline name is indicated on top of the panel.
Entropy 24 00248 g005
Figure 6. Power spectrum of the time evolution of the occurrence of the top 12 4-motifs of Ryanair. The label of the 4-motif is on top of each panel. Frequency peaks are detected for f = 0.14 and for its harmonics f = 0.28 and f = 0.43 day 1 . The main frequency f = 0.14 day 1 corresponds to a weekly cycle and the second and third harmonics correspond to a bi-weekly or three-weekly cycle.
Figure 6. Power spectrum of the time evolution of the occurrence of the top 12 4-motifs of Ryanair. The label of the 4-motif is on top of each panel. Frequency peaks are detected for f = 0.14 and for its harmonics f = 0.28 and f = 0.43 day 1 . The main frequency f = 0.14 day 1 corresponds to a weekly cycle and the second and third harmonics correspond to a bi-weekly or three-weekly cycle.
Entropy 24 00248 g006
Figure 7. Shape of the 10 top 4-motifs observed for the 9 largest airlines (see Table 2). The 4-motifs are labeled according to [20]. Left column groups 4-motifs with only unidirectional links. Right column groups 4-motifs with only bidirectional links. Center column groups 4-motifs with both unidirectional and bidirectional links.
Figure 7. Shape of the 10 top 4-motifs observed for the 9 largest airlines (see Table 2). The 4-motifs are labeled according to [20]. Left column groups 4-motifs with only unidirectional links. Right column groups 4-motifs with only bidirectional links. Center column groups 4-motifs with both unidirectional and bidirectional links.
Entropy 24 00248 g007
Figure 8. Average linkage hierarchical tree obtained by using the relation d i j = 2 ( 1 ρ i j ) as distance, where ρ i j is the Spearman correlation coefficient between the vectors of the average occurrence of each 4-motif of airlines i and j. The average is computed over the 365 days of the year.
Figure 8. Average linkage hierarchical tree obtained by using the relation d i j = 2 ( 1 ρ i j ) as distance, where ρ i j is the Spearman correlation coefficient between the vectors of the average occurrence of each 4-motif of airlines i and j. The average is computed over the 365 days of the year.
Entropy 24 00248 g008
Figure 9. Average linkage hierarchical tree obtained by using as distance the relation d i j = 2 ( 1 J ( n e t i , n e t j ) ) , where J ( n e t i , n e t j ) is the Jaccard measure of the overlap of directed links observed in weekly flight network i and j. The hierarchical tree is obtained for the 20th week of 2017 (from Monday 15 May to Sunday 21 May).
Figure 9. Average linkage hierarchical tree obtained by using as distance the relation d i j = 2 ( 1 J ( n e t i , n e t j ) ) , where J ( n e t i , n e t j ) is the Jaccard measure of the overlap of directed links observed in weekly flight network i and j. The hierarchical tree is obtained for the 20th week of 2017 (from Monday 15 May to Sunday 21 May).
Entropy 24 00248 g009
Table 1. Basic metrics of the flight networks of 50 largest airlines by number of flights performed. The degree of each airport is obtained as the average value of the degree observed for each day of the year. N is the number of nodes, m is the number of direct links, Min, Median, Mean, Max and Std are the minimum, median, mean, maximum and standard deviation of the average degree distribution. The Herfindal index is obtained from the average degree distribution as explained in the main text. Airlines are ordered according to the value of the Herfindal index.
Table 1. Basic metrics of the flight networks of 50 largest airlines by number of flights performed. The degree of each airport is obtained as the average value of the degree observed for each day of the year. N is the number of nodes, m is the number of direct links, Min, Median, Mean, Max and Std are the minimum, median, mean, maximum and standard deviation of the average degree distribution. The Herfindal index is obtained from the average degree distribution as explained in the main text. Airlines are ordered according to the value of the Herfindal index.
AirlineN2mMinMedianMeanMaxStdHerfindal
BEL65.08239.651.002.003.68119.5214.570.253
AFL63.87248.191.172.003.89123.3015.180.251
KLM143.69538.271.002.003.75268.5022.240.251
IBE60.39227.551.002.003.77110.8414.010.243
FIN78.79308.371.002.003.92147.9516.430.234
AUA82.96342.751.012.004.13155.9916.900.213
TAP76.26327.581.002.004.30144.4416.400.202
QTR67.19266.831.002.003.97117.9914.150.201
BTI36.78160.371.062.004.3767.6210.840.192
SWR78.68351.811.002.004.47148.6616.760.190
AUI54.36233.901.012.004.2998.6413.130.189
UAE63.07247.101.002.003.92103.4312.800.183
LOT66.62288.131.032.004.32119.2914.470.182
AFR156.12657.341.002.004.21270.1021.760.176
OAL26.62125.881.892.004.7248.199.040.174
RAM79.28327.751.002.004.13132.7814.750.172
EIN52.90238.281.072.004.4892.7212.460.164
BAW162.38660.511.002.004.07244.3720.010.155
AEE50.90236.931.012.004.5979.5811.080.141
PGT89.83459.421.012.005.10160.0017.080.138
THY245.361274.181.002.005.18449.5129.370.135
EZS42.87221.231.082.205.1561.8811.280.133
AZA77.52410.591.012.005.29138.8116.150.132
AEA50.64237.181.012.004.6977.3710.930.126
DLH175.871015.781.002.275.77289.7426.450.124
SXS33.31199.471.043.585.9146.179.710.113
ANE38.43198.791.272.175.1354.178.820.104
TRA54.48286.731.012.075.2475.5810.870.098
UAL31.36201.631.163.476.3947.358.790.091
LOG19.34109.381.623.335.6219.334.900.089
BCS59.85288.981.231.864.6260.628.460.080
NAX67.06401.901.022.215.9987.3812.480.079
HOP35.22214.971.062.186.0840.507.840.076
DAL35.56210.471.092.735.8736.187.610.075
SAS94.26652.471.092.626.89110.6316.920.074
AAL23.42157.981.554.866.6122.765.690.074
CFG41.51234.101.002.205.5344.997.830.074
TCX33.97222.881.003.566.0933.876.820.070
VLG86.65689.101.073.607.90153.7017.310.067
EWG71.22465.411.002.486.3873.0611.450.062
EXS34.71332.781.165.878.8235.448.220.059
VOE35.79230.881.144.215.8421.074.780.054
IBK51.23381.921.003.327.4643.999.310.050
TOM46.92370.301.004.557.5940.198.450.050
BEE46.38429.761.315.359.2443.379.740.045
WIF43.98325.591.695.577.2925.355.850.039
EZY110.221714.641.087.3215.44157.9721.530.027
WZZ100.23895.021.005.338.8760.0811.020.026
NJE90.32262.931.001.982.8816.032.770.022
RYR172.473183.891.009.6418.33206.0924.720.016
Table 2. Rank of the 10 top 4-motifs observed for the 9 biggest airlines. The rank is obtained by considering the daily occurrence of each motif averaged over all days of the year. The shape of these 4-motifs are shown in Figure 7.
Table 2. Rank of the 10 top 4-motifs observed for the 9 biggest airlines. The rank is obtained by considering the daily occurrence of each motif averaged over all days of the year. The shape of these 4-motifs are shown in Figure 7.
RankRYRDLHTHYEZYAFRSASBAWKLMVLG
1438243824382438243824382438243824382
2469828628646984374437428643744374
34958437443744958286495843742864958
443744958282437428246984682282286
5442813,2783013,278437028628243704698
613,278282437044283013,2784698304446
728630495813,260469843703028013,278
84064698469828649584428437028282
946824370444640646824446495421844370
104446444628495028049504702144950
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Milazzo, M.; Musciotto, F.; Miccichè, S.; Mantegna, R.N. Analysis of the Structure and Dynamics of European Flight Networks. Entropy 2022, 24, 248. https://doi.org/10.3390/e24020248

AMA Style

Milazzo M, Musciotto F, Miccichè S, Mantegna RN. Analysis of the Structure and Dynamics of European Flight Networks. Entropy. 2022; 24(2):248. https://doi.org/10.3390/e24020248

Chicago/Turabian Style

Milazzo, Matteo, Federico Musciotto, Salvatore Miccichè, and Rosario N. Mantegna. 2022. "Analysis of the Structure and Dynamics of European Flight Networks" Entropy 24, no. 2: 248. https://doi.org/10.3390/e24020248

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop