# Analysing the Vulnerability of Public Transport Networks

Authors: Mary Luz Mouronte-López
Journal: Journal of Advanced Transportation
 (2021)

---

## Abstract

This paper analyses the robustness of specific public transport networks. Common attributes and which of them have more influence on the networks’ vulnerability are established. Initially, the structural properties of the networks in two graphical representations (L-Space and P-Space) are checked. Afterwards, the spread of problems (traffic jams, etc.) are simulated, employing a model based on a propagation and recovery mechanism, similar to those used in the epidemiological processes. Next, the size of the largest connected subset of stops of the network (giant component) is measured. What is shown is that the faults randomly happened at stops or links, also displaying that those that occurred in the highest weighted links spread slower than others. These others appear at stops with the largest level of betweenness, degree, or eigenvector centralities and PageRank. The modification of the giant component, when several stops and links are removed, proves that the removal of stops with the highest interactive betweenness, PageRank, and degree centralities has the most significant influence on the network’s integrity. Some equivalences in the degree, betweenness, PageRank, and eigenvector centrality parameters have been found. All networks show high modularity with values of indexQ close to 1. The networks with the highest assortativity and lowest average number of stops are the ones which a passenger can use to travel directly to their destination, without any change. The Molloy–Reed parameter is higher than 2 in all networks, demonstrating that high integrity exists in them. All stops were characterized by low k-cores ≤3.


---

## Body

## 1. Introduction

The question of vulnerability of a public transport network (PTN) against problems in nodes or links is a matter of high interest. All nodes and links may accidentally fail with a similar probability, but the severity of the damage caused is higher when it occurs in some nodes or links rather than in others. The PTN is a key issue in the main cities of the world, having to operate efficiently on a constant basis. They must display a high tolerance to problems, keep the possible delays to a minimum, and manage the duration of a journey in accordance with the trip planning. With the purpose of persuading passengers to make use of the PTN, some strategies should be implemented to ensure the continued running of networks, the prevention of problems before, and the mitigation of losses during and after the failure [1, 2].This paper analyses the vulnerability of some PTNs. The study firstly focuses on the analysis of the structural properties of the PTN (betweenness centrality, degree, eigenvector centrality, PageRank, average path length,k-core, modularity, and assortativity). Next, the fault propagation is simulated (e.g., traffic jams caused by accidents or roadworks), employing different selection criteria of nodes and links, in order to detect the stops in whose immediate vicinity congestion is more prone to happening. A procedure based on epidemic spread algorithms is used. The effect of the structure of the network on the failure propagation is also examined. Finally, the size of the largest connected subset of stops of the PTN (giant component, GC) is estimated, which is calculated when several stops and links are disconnected according to certain criteria. Depending on the value of this magnitude, a serious issue could occur due to the existence of inaccessible routes (construction work or redesign activities in networks), which cause unavoidable stops or important delays during a trip. The analysis allows the identification of the most critical stops and links of a PTN and establishes a subnetwork of highly connected stops.Ten networks are analysed. They are Auckland in New Zealand, which consists of 5223 stops and 318,896 links; EMT in Madrid, Spain, which has 4636 stops and 2,070,508 links; BC Transit in Vancouver, Canada, which consists of 3981 stops and 702,197 links; Kolumbus in Rogaland, Norway, which has 3828 stops and 487,432 links; AVL, CFL, RGTR, and TICE in Luxembourg which contains 1372 nodes and 340,684 links; STAR in Rennes, France, which consists of 1415 stops and 9,477,213 links; Thunder Bay Transit, in Ontario, Canada, which contains 825 nodes and 78,247 links; TransAntofagasta, in Chile, which has 650 nodes and 724,362 links; Linja-Karjala Oy, in Kuopio, Finland, which contains 551 nodes and 63,339 links; and finally, CIT Chambly-Richelieu-Carignan, in Quebec, Canada, which consists of 346 nodes and 9366 links. The reason for the selection was to have networks of various sizes (small, medium, and large) and topologies.Specifically, the goals of this research are as follows:(1)
Analysis of the robustness of some PTNs: detection of commonalities and key aspects.(2)
Discovery of properties on a networkʼs vulnerability that can be generalized.Several investigations exist about PTN, some of them describing a comprehensive survey of statistical properties of PTN based on the data of several cities, and at the same time, presenting models that reproduce the best part of their properties [3]. Other research detects the critical nodes in PTN in order to enable contingency mechanics in them.Berche et al. [4] studied the PTNʼs resilience applying different attack scenarios, which range from random failure to targeted destruction. Some nodes are removed according to certain operating characteristics. The candidate nodes to be eliminated were selected according to their highest degree, closeness, graph, stress, and betweenness centralities, as well as the largest clustering coefficient and next nearest neighbouring number. In [5], an explanation exists of a procedure for the detection of critical airports in the global air transport network (ATN), which is based on simulating an attack on specific airports using several adaptive selection criteria. Wandelt et al. [6] evaluated several strategies for dismantling of networks and identified large heterogeneities in their performance. The authors show that the use of the interactive betweenness creates far stronger attacks. Because of the high complexity of this algorithm, the authors recommend its use on small networks. After, in [7], the authors propose approximations to interactive betweenness in order to improve the computational cost of the algorithm. Other research analyses the effect that the addition of new links has on the robustness of the networks or design methodologies in order to detect the most vulnerable links in PTN [8].The robustness of a network can be estimated in terms of passenger welfare [9]. Additionally, the accessibility of the PTN is evaluated by several types/examples of research. Albacete et al. [10] compared and applied two location-based methods to analyse PTN accessibility. The results had considerable implications for transport policy making. In [11], a thorough assessment is carried out, focusing on the possible applications of an optimization-based scheme for the building of several timetables and proactive railway traffic management. This is carried out over a large network, employing stochastic disturbances. Wandelt et al. [12] proposed a framework to evaluate and improve the robustness of transportation systems by exploiting the existence of communities. The method is applied to several real-world transportation systems. Sun et al. [13] studied the resilience of cities when disruptions in their airports happen.Relationships between robustness and accessibility are also studied [14, 15]. Other authors suggested methods and spread models which analyse the performance of networks. He et al. [16] defined a mathematical model to evaluate the effect of the propagation of the risk of failure in multimodal transport networks. Baspinar and Koyuncu [17] proposed a new model for the air transportation network under stress and defined parameters in order to describe not only air sectors, but also airports, as well as flights. They used epidemic spreading processes by assuming that the characteristics of disease spreading and the delay propagation are similar. Akdere et al. [18] studied the reliable data dissemination in the context of wireless sensor network environments. It shows the applicability of epidemic spread algorithms in those environments and carries out a comparative performance analysis of several mechanisms in terms of message delivery rate, average message latency, and messaging overhead on the network. As other research suggests, Nekovee [19] presented a new epidemic algorithm for information dissemination in highly dynamic and intermittently connected vehicular ad hoc networks (VANET). It shows through realistic simulations of highway traffic that the proposed algorithm is suitable to achieve a reliable and efficient level of information transmission in a context with frequent network fragmentation and large density variations.This paper is organized as follows: firstly, resources and methods used in this work are explained in detail. It focuses on the used data and the developed software programs, which correspond with the design of the network and the employed parameters to characterize the networks’ topology. The propagation model and the study of evolution of the GC when certain elements are removed are also studied. Afterwards, the results, conclusions, and future projects are described and discussed.


## 2. Materials and Methods

### 2.1. Overview of Used Resources

Information of stops and routes of Auckland, EMT, BC Transit, Kolumbus, STAR, AVL, CFL, RGTR, TICE, Thunder Bay Transit, TransAntofagasta, Linja-Karjala Oy, and CIT networks was retrieved from public sources (see the Data Availability section for details). Several programs in R [20] and Python [21] were implemented, using the R.3.6.0 and 3.8.3 version, respectively. We use the Python networkx and pathpy packages and the igraph R package. The programs followed the typical development life cycle with stages of specification, design, and testing.


### 2.2. Overview of Used Methods

The structural properties of the networks are analysed. A new propagation algorithm, which was designed and undertaken in the networks in order to simulate the spread of a failure, is detailed. This is followed by a description of how the evolution of the GC is calculated, with the aim of studying the specific impact a problem might have on each network.


### 2.3. Study of Structural Properties

A PTN can be represented in two topological spaces L-Space [22, 23] and P-Space [2, 24]. In both spaces, the network is mapped as a graph G = (N; L), where N is the set of nodes symbolizing the stops and L is the set of links established between them. In L-Space, one node symbolizes one stop, and one link means a union between two consecutive stops, which tells us there is a link between two stops, if one stop is the successor of the other on a route. This space aims to show the geographical proximity between stops. In the P-Space [2, 24], one node represents one stop, and one link joins a pair of stops, if at least one route provides direct service between them. A link means that passengers can take at least one route for a direct trip between two stops. If travellers have to exchange routes, then the pair of stops is joined by more than one link. This space aims to show the transfer relationship between routes.In L-Space and P-Space, the average values of betweenness (<BC>) and eigenvector (<EC>) centralities, the PageRank <PR>, and degree (<k>) of nodes as well as the minimum distance between nodes distributions were estimated. Pagerank was measured considering a damping factor equal to 0.85. In order to determine if the degree distributions, similar in a way to what happens in other networks [3, 25], followed a power law function, this characteristic was also analysed. Next, all these magnitudes are defined:(i)
The degree of a nodei, ki, for an undirected graph, G, such as a PTN, is [26](1)ki=∑j=1NAij,whereAij is the element ij of the adjacency matrix, for example, Aij=1 if the node iis linked to node j and 0 otherwise.(ii)
The minimum distance between two nodesi,j in G, l, is the length of the shortest path between them.(iii)
The betweenness centrality of a nodeiin G, BCi is [27](2)BCi=∑u≠i≠wσu,wiσu,w,whereσu,w is the total number of shortest paths from node u to node w and σu,wi is the number of those paths that pass through i.(iv)
Regarding the eigenvector centrality of a nodei in G, ECi [28, 29]:λ1,λ2,λ3,,…,λN are the eigenvalues of the adjacency matrix A=Aij of G. Then, the largest eigenvalue of matrix A is λmaxwith an eigenvector e=e1,e2,…,eNT such that λmax∗ei=∑j=1NAij∗ej. The eigenvector centrality for node i denoted as ECi can be defined as(3)ECi=1λmax∑j=1NAij∗ej.(v)
Pagerank, PR, of a nodei in G is [28, 30](4)PRi=qN+1−q∑j:j⟶iPRjkoutj,i=1,2,3,…,N,where N is the number of nodes in G, PRj is the PageRank of a node j, and koutj is the outdegree of node j. We add PRj/koutj for all nodes with a link ending in i. In the case of the PTN, it is considered that G is an undirected graph; therefore, koutj=kj. q is the damping parameter, ∈ [0, 1].Regarding the modularity of the network, in L-Space, this was calculated employing the generate overlapping cluster generator (OCG) method [31]. The clusters were originally hierarchically joined together, optimizing the modularity of the partition, resulting in overlapping clusters. Similarly, to what happens in the protein networks, the uniqueness of node categorization prevents from disclosing the implication of some stops in various trip processes. In L-Space, the assortativity parameter [32] which determines the tendency of nodes being connected to other similar ones was also estimated. Additionally, the Molloy–Reed parameter (MRP) [33] was calculated to know the structural integrity of the networks.


### 2.4. Simulation of Faults Propagation

In the L-Space, the propagation of failures in each network is simulated using an infection and recovery mechanism, which consists of the following steps:Calculate the total number of nodes in the network (N),Initialize the infection rate (IR = 0),Infect the first node,Calculate the total number of infected nodes in the network (µ′),Initialize time (t = 1),while (IR < 0.60) and (t < T):{Obtain all uninfected nodes in the network,uinii=1i=µ, µ being the total number of uninfected nodes.Fori = 1 to µ{Obtain the neighbour nodes ofuininnuini−jj=1j=γi, γi being the total number of neighbours of uiniForj = 1 to γi:{If RandomBinomial (µ′, pI) > 0 nnuini−j is still infectedElseTo disinfectnnuini−j}}Calculate how many nodes are infected (δ),Fori = 1 to δ:{If RandomBinomial then (1,pR) > 0 thenTo disinfectini, pR being the recuperation probability,ElseTo infectini}Calculate the total number of infected nodes in the network (µ′)Calculate IR as:IR=µ′/N}where RandomBinomial (n, p) is a generator of random numbers, which are calculated from the binomial distribution specified by the number of trials n and the probability of success for each trial p. The np.random.binomial function in python is used.Several simulations were carried out until a timeT and without reaching 60% percent of the total of infected nodes in each transport network. pI and pR were modified in the range 0 to 0.1 in steps of 0.0025 (values were considered within a range of 0-0.1, because with pI > 0.1, all networks very quickly reached 60% infection (they were saturated in a short time for all groups). It was checked if some differences existed when the propagation process started in nodes with specific characteristics, such as nodes with the highest betweenness (B) and eigenvector centralities (E), nodes with the largest PageRank (P) and degree (D), and nodes randomly selected (A). Pagerank was estimated considering a damping factor equal to 0.85. Once the simulation was concluded, the existing differences between groups (B, E, P, D, A) were analysed. The following procedure was carried out:(1)
The normality of distributions was tested using the D’Agostino test [34] with a significance level α = 0.05. This test provides very effective results when it is applied to large size samples [34, 35]. The considered hypothesis was as follows:(i)
H0: “the samples came from a normal distribution.”(ii)
Ha: “the samples did not come from a normal distribution.”Ifp value ≤ α, H0 is rejected and Ha is taken, else H0 is accepted.(2)
If normality existed, the homoscedasticity of variances would be studied in each distribution using the Breusch–Pagan test [36] with a significance level of α = 0.05. The following hypothesis should be taken into account:(i)
H0: “the variance was constant.”(ii)
Ha: “the variance was not constant.”Again, ifp value ≤ α, H0 should be rejected and Ha should be accepted; otherwise, Ha should be taken.(3)
If there was no normality or homoscedasticity in the distributions, the Kruskal–Wallis test [37] should be executed in order to detect whether the population distributions were identical, or at least one differed from the rest. The used significance level was α = 0.05. The hypothesis employed was as follows:(i)
Null hypothesis: “the groups were from identical populations.”(ii)
Alternative hypothesis: “at least one of the groups comes from a different population than the others.”(4)
In the next step, the Wilcoxon rank sum test [37] was executed with a significance level α = 0.05 in order to realize a pairwise comparison between groups. The hypothesis was as follows:(i)
Null hypothesis: “both groups had the same distribution.”(ii)
Alternative hypothesis: “both groups have a different distribution.”(5)
If there are normality and homoscedasticity, the analysis of variance (ANOVA) method was realized.


### 2.5. Analysis of Evolution of the GC

The size of the GC describes the largest fraction of overall nodes such that any pair of them is linked through a path. It evaluates the largest extension that a route (or a union of them) provides in terms of the available stops a traveller can reach from an origin inside the GC.In the L-Space, for each network, the size of GC was calculated as nodes and links were removed according to the following criteria: links (criterion 1) or nodes (criterion 2) randomly eliminated; nodes removed in descending order of betweenness centrality (criterion 3), degree (criterion 4), eigenvector centrality (criterion 5), and PageRank (considering a damping factor equal to 0.85) (criterion 6); and the highest weighted links removed (criterion 7) following a decreasing sequence. In this criterion, a weighted graph was built where the weight of each link,wij, represented the number of links between the nodes i, j and the highest number of links between two nodes in the entire network. Finally, the interactive betweenness calculation was carried out (criterion 8). This criterion, instead of removing the nodes in decreasing order of the static betweenness, recomputes the betweenness after the removal of a node.It is relevant to analyse which of the used criteria allow us to achieve an 80% reduction in the size of GC in less time. It is true that when the network is fragmented, the stops belonging to the GC are linking a fraction of the PTN, which can still be functional in certain cases. However, it is also relevant to observe the evolution of size over time in order to estimate the attack tolerance. Several commonalities and differences between PTN can be detected. The relevance of studying the evolution of PTN is also corroborated by the interdependencies between PTN and the accessibility to other infrastructures such as hospitals, business centers, schools, and airports. Therefore, the fluctuations in the availability of stops and routes might greatly impact urban life.


## 2.1. Overview of Used Resources

Information of stops and routes of Auckland, EMT, BC Transit, Kolumbus, STAR, AVL, CFL, RGTR, TICE, Thunder Bay Transit, TransAntofagasta, Linja-Karjala Oy, and CIT networks was retrieved from public sources (see the Data Availability section for details). Several programs in R [20] and Python [21] were implemented, using the R.3.6.0 and 3.8.3 version, respectively. We use the Python networkx and pathpy packages and the igraph R package. The programs followed the typical development life cycle with stages of specification, design, and testing.


## 2.2. Overview of Used Methods

The structural properties of the networks are analysed. A new propagation algorithm, which was designed and undertaken in the networks in order to simulate the spread of a failure, is detailed. This is followed by a description of how the evolution of the GC is calculated, with the aim of studying the specific impact a problem might have on each network.


## 2.3. Study of Structural Properties

A PTN can be represented in two topological spaces L-Space [22, 23] and P-Space [2, 24]. In both spaces, the network is mapped as a graph G = (N; L), where N is the set of nodes symbolizing the stops and L is the set of links established between them. In L-Space, one node symbolizes one stop, and one link means a union between two consecutive stops, which tells us there is a link between two stops, if one stop is the successor of the other on a route. This space aims to show the geographical proximity between stops. In the P-Space [2, 24], one node represents one stop, and one link joins a pair of stops, if at least one route provides direct service between them. A link means that passengers can take at least one route for a direct trip between two stops. If travellers have to exchange routes, then the pair of stops is joined by more than one link. This space aims to show the transfer relationship between routes.In L-Space and P-Space, the average values of betweenness (<BC>) and eigenvector (<EC>) centralities, the PageRank <PR>, and degree (<k>) of nodes as well as the minimum distance between nodes distributions were estimated. Pagerank was measured considering a damping factor equal to 0.85. In order to determine if the degree distributions, similar in a way to what happens in other networks [3, 25], followed a power law function, this characteristic was also analysed. Next, all these magnitudes are defined:(i)
The degree of a nodei, ki, for an undirected graph, G, such as a PTN, is [26](1)ki=∑j=1NAij,whereAij is the element ij of the adjacency matrix, for example, Aij=1 if the node iis linked to node j and 0 otherwise.(ii)
The minimum distance between two nodesi,j in G, l, is the length of the shortest path between them.(iii)
The betweenness centrality of a nodeiin G, BCi is [27](2)BCi=∑u≠i≠wσu,wiσu,w,whereσu,w is the total number of shortest paths from node u to node w and σu,wi is the number of those paths that pass through i.(iv)
Regarding the eigenvector centrality of a nodei in G, ECi [28, 29]:λ1,λ2,λ3,,…,λN are the eigenvalues of the adjacency matrix A=Aij of G. Then, the largest eigenvalue of matrix A is λmaxwith an eigenvector e=e1,e2,…,eNT such that λmax∗ei=∑j=1NAij∗ej. The eigenvector centrality for node i denoted as ECi can be defined as(3)ECi=1λmax∑j=1NAij∗ej.(v)
Pagerank, PR, of a nodei in G is [28, 30](4)PRi=qN+1−q∑j:j⟶iPRjkoutj,i=1,2,3,…,N,where N is the number of nodes in G, PRj is the PageRank of a node j, and koutj is the outdegree of node j. We add PRj/koutj for all nodes with a link ending in i. In the case of the PTN, it is considered that G is an undirected graph; therefore, koutj=kj. q is the damping parameter, ∈ [0, 1].Regarding the modularity of the network, in L-Space, this was calculated employing the generate overlapping cluster generator (OCG) method [31]. The clusters were originally hierarchically joined together, optimizing the modularity of the partition, resulting in overlapping clusters. Similarly, to what happens in the protein networks, the uniqueness of node categorization prevents from disclosing the implication of some stops in various trip processes. In L-Space, the assortativity parameter [32] which determines the tendency of nodes being connected to other similar ones was also estimated. Additionally, the Molloy–Reed parameter (MRP) [33] was calculated to know the structural integrity of the networks.


## 2.4. Simulation of Faults Propagation

In the L-Space, the propagation of failures in each network is simulated using an infection and recovery mechanism, which consists of the following steps:Calculate the total number of nodes in the network (N),Initialize the infection rate (IR = 0),Infect the first node,Calculate the total number of infected nodes in the network (µ′),Initialize time (t = 1),while (IR < 0.60) and (t < T):{Obtain all uninfected nodes in the network,uinii=1i=µ, µ being the total number of uninfected nodes.Fori = 1 to µ{Obtain the neighbour nodes ofuininnuini−jj=1j=γi, γi being the total number of neighbours of uiniForj = 1 to γi:{If RandomBinomial (µ′, pI) > 0 nnuini−j is still infectedElseTo disinfectnnuini−j}}Calculate how many nodes are infected (δ),Fori = 1 to δ:{If RandomBinomial then (1,pR) > 0 thenTo disinfectini, pR being the recuperation probability,ElseTo infectini}Calculate the total number of infected nodes in the network (µ′)Calculate IR as:IR=µ′/N}where RandomBinomial (n, p) is a generator of random numbers, which are calculated from the binomial distribution specified by the number of trials n and the probability of success for each trial p. The np.random.binomial function in python is used.Several simulations were carried out until a timeT and without reaching 60% percent of the total of infected nodes in each transport network. pI and pR were modified in the range 0 to 0.1 in steps of 0.0025 (values were considered within a range of 0-0.1, because with pI > 0.1, all networks very quickly reached 60% infection (they were saturated in a short time for all groups). It was checked if some differences existed when the propagation process started in nodes with specific characteristics, such as nodes with the highest betweenness (B) and eigenvector centralities (E), nodes with the largest PageRank (P) and degree (D), and nodes randomly selected (A). Pagerank was estimated considering a damping factor equal to 0.85. Once the simulation was concluded, the existing differences between groups (B, E, P, D, A) were analysed. The following procedure was carried out:(1)
The normality of distributions was tested using the D’Agostino test [34] with a significance level α = 0.05. This test provides very effective results when it is applied to large size samples [34, 35]. The considered hypothesis was as follows:(i)
H0: “the samples came from a normal distribution.”(ii)
Ha: “the samples did not come from a normal distribution.”Ifp value ≤ α, H0 is rejected and Ha is taken, else H0 is accepted.(2)
If normality existed, the homoscedasticity of variances would be studied in each distribution using the Breusch–Pagan test [36] with a significance level of α = 0.05. The following hypothesis should be taken into account:(i)
H0: “the variance was constant.”(ii)
Ha: “the variance was not constant.”Again, ifp value ≤ α, H0 should be rejected and Ha should be accepted; otherwise, Ha should be taken.(3)
If there was no normality or homoscedasticity in the distributions, the Kruskal–Wallis test [37] should be executed in order to detect whether the population distributions were identical, or at least one differed from the rest. The used significance level was α = 0.05. The hypothesis employed was as follows:(i)
Null hypothesis: “the groups were from identical populations.”(ii)
Alternative hypothesis: “at least one of the groups comes from a different population than the others.”(4)
In the next step, the Wilcoxon rank sum test [37] was executed with a significance level α = 0.05 in order to realize a pairwise comparison between groups. The hypothesis was as follows:(i)
Null hypothesis: “both groups had the same distribution.”(ii)
Alternative hypothesis: “both groups have a different distribution.”(5)
If there are normality and homoscedasticity, the analysis of variance (ANOVA) method was realized.


## 2.5. Analysis of Evolution of the GC

The size of the GC describes the largest fraction of overall nodes such that any pair of them is linked through a path. It evaluates the largest extension that a route (or a union of them) provides in terms of the available stops a traveller can reach from an origin inside the GC.In the L-Space, for each network, the size of GC was calculated as nodes and links were removed according to the following criteria: links (criterion 1) or nodes (criterion 2) randomly eliminated; nodes removed in descending order of betweenness centrality (criterion 3), degree (criterion 4), eigenvector centrality (criterion 5), and PageRank (considering a damping factor equal to 0.85) (criterion 6); and the highest weighted links removed (criterion 7) following a decreasing sequence. In this criterion, a weighted graph was built where the weight of each link,wij, represented the number of links between the nodes i, j and the highest number of links between two nodes in the entire network. Finally, the interactive betweenness calculation was carried out (criterion 8). This criterion, instead of removing the nodes in decreasing order of the static betweenness, recomputes the betweenness after the removal of a node.It is relevant to analyse which of the used criteria allow us to achieve an 80% reduction in the size of GC in less time. It is true that when the network is fragmented, the stops belonging to the GC are linking a fraction of the PTN, which can still be functional in certain cases. However, it is also relevant to observe the evolution of size over time in order to estimate the attack tolerance. Several commonalities and differences between PTN can be detected. The relevance of studying the evolution of PTN is also corroborated by the interdependencies between PTN and the accessibility to other infrastructures such as hospitals, business centers, schools, and airports. Therefore, the fluctuations in the availability of stops and routes might greatly impact urban life.


## 3. Results and Discussion

### 3.1. Study of Structural Properties

The PageRank (with a damping parameter equal to 0.85), betweenness and eigenvector centralities, degree, and minimum distance characteristics were analysed in the L-Space and P-Space. Figures1 and 2 show the L-Space for all studied networks. The average values of degree, PageRank, eigenvector, and betweenness centralities for this space can be observed in Table 1. Those correspond to degrees, which are similar in all networks. The average eigenvector and PageRank centralities are lower than 0.002, having very small values in all networks. Therefore, the failures randomly happening in a node will have a low probability of arriving at any other. The same happens with respect to betweenness centrality, presenting small magnitudes (a mean value minor than 0.025), which denotes a slight average in how a node plays a bridging role in a network. Table 1 also shows the value of MRP in all networks, proving the existence in them of a GC (MRP >2).Figure 1
In L-Space, for EMT (a), BC Transit (b), Kolumbus (c), Auckland (d), and TransAntofagasta (e) networks.
(a)(b)(c)(d)(e)Figure 2
In L-Space, for Linja-Karjala Oy (a), AVL, CFL, RGTR, and TICE (b), STAR (c), CIT (d), and Thunder Bay Transit (e) networks.
(a)(b)(c)(d)(e)Table 1
Structural properties in L-Space. Average Betweenness Centrality (<BC>), average degree (<k >), average Eigenvector Centrality (<EC>), average PageRank (with a damping factor = 0.85) (<PR>), average path Length (<l>), diameter (d), density, and MRP.
NetworkNumber of nodesNumber of links<BC><k><EC><PR><l>dDensityMRPAuckland5,223318,8960.006552.360710.002170.000192.695231510.000452.66878EMT4,6362,070,5080.003832.711400.001810.000221.00000610.000583.21782BC Transit3,981702,1970.002332.288940.001670.000259.11427830.000582.54653Kolumbus3,829487,4320.002582.595460.002160.000262.89389980.000683.24314STAR1,415661,1500.006642.816960.005590.0007110.0000370.001993.54340AVL, CFL, RGTR, and TICE1,372340,6840.007062.810500.005640.000721.00000340.002053.79201Thunder Bay Transit82578,2470.034992.176970.005600.0012129.79920.002642.35412TransAntofagasta650724,3620.024362.963080.009020.0015416.78525510.004563.46729Linja-Karjala Oy55163,3990.016112.558980.014970.001814.50476290.004653.18865CIT3469,3660.027142.248550.017000.002892.33333290.006512.75804In L-Space, an equivalence between the degree distributions of the PTN of some cities was pointed out by Von Feber et al. [3], where the authors showed that this distribution followed a power law. This characteristic was checked for the analysed networks, fitting their degree distribution to the function:(5)Pk∼k−α.The obtained values for the fitting are shown in Table2. It is observed, according to the R square value, that those obtained for the fitting were not good (R squared is not close to 1). Figure 3 depicts, in L-Space, the degree distribution for the analysed networks. It is clear that all of them have a low value of k (kmax ≤ 17) compared with the total links. Therefore, the failures occurring at random on a node will have a high probability of having a small impact. Even so, these networks might be vulnerable to disruptions since a very high amount of nodes is joined to very few other nodes (low degree).Table 2
Fitting in L-Space of the network degree distribution to a power law functionPk∼k−α.
L-SpaceNetworkαR squareLinja-Karjala Oy0.1580.600Kolumbus0.4590.330EMT0.6220.239Auckland0.6030.237BC Transit0.4120.161AVL, CFL, RGTR, and TICE0.2360.156CIT0.3020.156Thunder Bay Transit0.4060.068TransAntofagasta0.1820.044STAR0.1580.041Figure 3
Degree distribution in linear (x-axis) and linear scale (y-axis) for the following networks: EMT (a), BC Transit (b), Kolumbus (c), Auckland (d), TransAntofagasta (e), Linja-Karjala Oy (f), AVL, CFL, RGTR, TICE (g), STAR (h), CIT (i), and Thunder Bay Transit (j).
(a)(b)(c)(d)(e)(f)(g)(h)(i)(j)It is clear from the above that the values of betweenness, eigenvector, and PageRank centralities in the analysed networks ensure that they are highly robust against random failures.The lowest average minimum distance between nodes shows a relevant variability in all networks. TransAntofagasta and EMT as well as AVL, CFL, RGTR, and TICE networks have the lowest value. The highest value corresponds to Thunder Bay Transit network.The modularity of the network is also evaluated in L-Space, employing the generate overlapping cluster generator (OCG) method [31]. It can be noted in Table 3 that all networks show a high value of Q with magnitudes near to 1, which demonstrates a strong community structure in these networks. Romano et al. [38] showed that network efficiency is modularity dependent, with the highest values of transmission occurring at intermediate levels and low values of modularity. High values tend to negatively influence transmission in networks [38]. Therefore, the high modularity in PTN could help to slow down the propagation of a failure in the network.Table 3
Q Modularity computed using the overlapping cluster generator (OCG) method in L-Space.
NetworkQBC Transit0.99750Auckland0.99610EMT0.98550Thunder Bay Transit0.99100Kolumbus0.98730STAR0.97020AVL, CFL, RGTR, and TICE0.96880CIT0.95880TransAntofagasta0.93910Linja-Karjala Oy0.93960Table4 depicts the values of the assortativity coefficient [32]. All networks show assortative higher than 0.4, demonstrating that a significant chance exists in which fractions of stops will randomly join with stations of the same degree. The nodes appear to join with a specific preference. A failure in a node with the highest degree will not have much negative impact, since other nodes with the highest degree are still connected to ones with a similar degree.Table 4
Assortativity coefficients in L-Space.
NetworkAssortativityKolumbus0.90599STAR0.89463BC Transit0.86067Auckland0.84123Linja-Karjala Oy0.81336AVL, CFL, RGTR, and TICE0.72469EMT0.67211TransAntofagasta0.58252Thunder Bay Transit0.45575CIT0.42243The highestk-core in all networks was 3 or 2, demonstrating that the largest subgraphs exist where every node is connected to at least 3 or 2 other nodes within the network. Most nodes show a small and similar k-core (≤2).With respect to P-Space, Table5 shows the average values of eigenvector and betweenness centralities, PageRank, degree, minimum distance between nodes, diameter, and density.Table 5
Structural properties in P-Space. Average Betweenness Centrality (<BC>), average degree (<k>), average Eigenvector Centrality (<EC>), average PageRank (with a damping factor = 0.85) (<PR>), average path Length (<l>), diameter (d), and density.
NetworkNumber of nodesNumber of links<BC><k><EC><PR><l>dDensityAuckland5,2238,079,0950.00039101.285080.023860.000192.6914370.01940EMT4,01025,263,1080.0006547.669830.046850.000253.5923870.01189BC Transit3,98115,737,7290.0001770.767140.034150.000252.6548160.01778Kolumbus3,82911,317,7140.0002878.983550.025840.000262.9688970.02063STAR1,4159,477,2130.0013743.112360.018620.000711.4827660.03049AVL, CFL, RGTR, and TICE1,3724,254,6560.0012048.425700.017910.000731.0000060.03532Thunder Bay Transit8252,014,7790.0013375.510300.021970.001212.0940.09164TransAntofagasta65031,171,3070.00108199.747690.036660.001541.7012330.30778Linja-Karjala Oy5511,167,6820.0018064.268600.035000.001811.4761940.11685CIT346135,7400.0027433.934780.034450.002891.0000040.103578The degree distributions of the analysed networks were also fitted to a power law function according to equation (5). The results are shown in Table 6. It can be noted that, as in the case of L-Space, in accordance with the obtained R square value, the fitting is not good. P-Space allows us to analyse interconnections between routes and transfers in the network. In this space, k reveals how many stops a person can travel through by using only one route. The network with the highest value average k, <k>, is TransAntofagasta. The obtained values for eigenvector and PageRank centralities are very low and similar in all networks, reinforcing the idea exposed in the analysis of L-space that failures randomly occurring in a node have a low probability of quickly arriving at any other one. Many possible transfer nodes exist, making the values for betweenness centrality very small and achieving high robustness against random failures on a node. The PTNs with lower <l> are the best connected in this space, having a minor level of vulnerability. Figure 4 shows, in P-Space, the degree distribution for the analysed networks.Table 6
Fitting in P-Space of the network degree distribution to a power law function.
NetworkαR squareEMT0.5130.180Kolumbus0.1690.166Auckland0.2040.077BC Transit0.1770.068Thunder Bay Transit0.2750.047TransAntofagasta0.1630.044STAR0.1180.037CIT0.1090.036AVL, CFL, RGTR, and TICE0.0610.022Linja-Karjala Oy0.0150.001Figure 4
Degree distributions in P-Space in linear (x-axis) and linear scale (y-axis) for the following networks: EMT (a), BC Transit (b), Kolumbus (c), Auckland (d), TransAntofagasta (e), Linja-Karjala Oy (f), AVL, CFL, RGTR, and TICE (g), STAR (h), CIT (i), and Thunder Bay Transit (j)
(a)(b)(c)(d)(e)(f)(g)(h)(i)(j)


### 3.2. Simulation of Faults Propagation

Failure propagation was simulated in L-Space using an infection and recovery algorithm.(1)
Firstly, a total of 100 experiments for each selection criterion and duple (pI, pR) were undertaken in each network. In each experiment, the simulation was implemented until t = T or until the percentage of infected nodes was higher or equal to 60%. The values of t were modified in steps of 10 units.(2)
The second step was to check the normality of distributions of the infection rate (IR), in all groups (A, B, D, E, P) using the D’Agostino test for anyT, which took different values until reaching a level propagation ≥60% of the network. The results show that p value was not higher than 0.05 in all cases. Since not all distributions presented normality, it was not necessary to study the homoscedasticity of variances applying the Breusch–Pagan Test. The Kruskal–Wallis test could be used for the detection of differences between groups.(3)
The third step was using the Kruskal–Wallis test to verify for eachT, whether all groups (A, B, D, E, P) were identical or by contrast, one of the groups tended to give observations that were different from those of other groups. Again, the p value obtained was lower than 0.05. The null hypothesis had to be rejected because at least one group had a different distribution.(4)
The fourth step was carried out utilising the Wilcoxon rank sum test in order to compare paired groups; we also applied the Bonferroni correction, both with a significance levelα = 0.05.In all analysed duos (pI, pR), for EMT, there was a significant difference (p value ≤ 0.05) between all groups except between groups D-E. With respect to BC Transit network, there were dissimilarities (p value ≤ 0.05) between some groups except between groups D-E, D-P, and E-P. Regarding the Kolumbus network, a relevant similarity is found between groups B-E and P-D. In the Auckland network, similarities happened between groups P-D, P-B, and D-B. For pI = 0.05 and pR = 0.0025, in T = 100, Figure 5 shows the results of the propagation algorithm in EMT, BC Transit, Kolumbus, and Auckland networks. For the same values of pI and pR and T = 65, Figure 6 shows the outcomes for STAR and AVL, CFL, RGTR, and TICE networks. Figure 7 depicts the results for TransAntofagasta, Linja-Karjala Oy, CIT, and Thunder Bay Transit networks. For those values of pI, pR, and T, Tables 7–16 show the obtained p value in the Wilcoxon rank sum test in each network.Figure 5
In L-Space for (a) EMT, (b) BC Transit, (c) Kolumbus, and (d) Auckland networks.pI = 0.05, pR = 0.0025, and T = 150 (boxplot diagram for A, B, D, E, and P groups).
(a)(b)(c)(d)Figure 6
In L-Space for (a) STAR and (b) AVL, CFL, RGTR, and TICE networks.pI = 0.05, pR = 0.0025, and T = 65 (boxplot diagram for A, B, D, E, and P groups).
(a)(b)Figure 7
In L-Space for (a) TransAntofagasta, (b) Linja-Karjala Oy, (c) CIT, and (d) Thunder Bay Transit.pI = 0.05, pR = 0.0025, and T = 65 (boxplot diagram for A, B, D, E, and P groups).
(a)(b)(c)(d)Table 7
In L-Space, for BC Transit network,pI = 0.05, pR = 0.0025, and T = 150, p-value of the Wilcoxon rank sum test (resultant p-value per Bonferroni correction is indicated in brackets).
ABDEB<5.70e − 11 (8.60e − 10)D1.10e − 14 (1.60e − 13)1.60e − 09 (2.30e − 08)E2.40e − 14 (3.50e − 13)3.70e − 09 (5.50e − 08)0.90000 (1.00000)P2.20e − 14 (3.30e − 13)1.30e − 08 (2.00e − 07)0.27000 (1.00000)0.34000 (1.00000)Table 8
In L-Space, for Kolumbus network,pI = 0.05, pR = 0.0025, and T = 150, p-value of the Wilcoxon rank sum test (resultant p-value per Bonferroni correction is indicated in brackets).
ABDEB<1.40e − 13 (2.00e − 12)D0.00473 (0.071)<2.00e − 16 (1.10e − 15)E3.00e − 14 (4.50e − 13)0.44174 (1.000)<2.00e − 16 (1.00e − 15)P0.00554 (0.083)<2.00e − 16 (1.10e − 15)0.99123 (1.000)<2.00e − 16 (1.00e − 15)Table 9
In L-Space, for EMT network,pI = 0.05, pR = 0.0025, and T = 150, p-value of in the Wilcoxon rank sum test (resultant p-value per Bonferroni correction is indicated in brackets).
ABDEB<2.00e − 16 (<2.00e − 16)D<2.00e − 16 (<2.00e − 16)0.00027 (0.0041)E<2.00e − 16 (<2.00e − 16)1.40e − 06 (2.00e − 05)0.30986 (1.00000)P<2.00e − 16 (<2.00e − 16)<2.00e − 16 (2.00e − 16)<2.00e − 16 (<2.00e − 16)<2.00e − 16 (<2.00e − 16)Table 10
In L-Space, for Auckland network,pI = 0.05, pR = 0.0025, and T = 150, p-value of the Wilcoxon rank sum test (resultant p-value per Bonferroni correction is indicated in brackets).
ABDEB<2.00e − 16 (<2.00e − 16)D<2.00e − 16 (<2.00e − 16)0.49000 (1.00000)E<2.00e − 16 (<2.00e − 16)7.40e − 14 (1.10e − 12)3.60e − 12 (5.30e − 11)P<2.00e − 16 (<2.00e − 16)0.54000 (1.00000)0.96000 (1.00000)1.00e − 12 (1.50e − 11)Table 11
In L-Space, for STAR network,pI = 0.05, pR = 0.0025, and T = 65, p-value of the Wilcoxon rank sum test (resultant p-value per Bonferroni correction is indicated in brackets).
ABDEB1.20e − 10 (1.20e − 09)D5.70e − 12 (5.70e − 11)0.056 (0.56)E2.90e − 14 (2.90e − 13)0.752 (1.00000)0.007 (0.07000)P1.20e − 12 (1.20e − 11)0.193 (1.00000)0.689 (1.00000)0.01900 (0.19000)Table 12
In L-Space, for AVL, CFL, RGTR, and TICE networks,pI = 0.05, pR = 0.0025, and T = 65, p-value of the Wilcoxon rank sum test (resultant p-value per Bonferroni correction is indicated in brackets).
ABDEB2.30e − 14 (2.3e − 13)D2.00e − 16 (4.8e − 16)1.70e − 07 (1.7e − 06)E< 2.00e − 16 (3.1e − 16)9.00e − 06 (9.0e − 05)0.19 (1.00000)P< 2.00e − 16 (<2e − 16)0.19 (1.00000)2.70e − 07 (2.7e − 06)4.30e − 05 (0.00043)Table 13
In L-Space, for TransAntofagasta network,pI = 0.05, pR = 0.0025, and T = 65, p-value of the Wilcoxon rank sum test (resultant p-value per Bonferroni correction is indicated in brackets).
ABDEB0.0013 (0.008)D1.80e − 11 (1.10e − 10)5.50e − 07 (3.30e − 06)E1.70e − 14 (1.70e − 13)2.90e − 14 (2.90e − 13)<2.00e − 16 (<2.00e − 16)P6.50e − 12 (3.90e − 11)1.50e − 07 (8.90e − 07)0.5348 (1.00000)<2.00e − 16 (<2.00e − 16)Table 14
In L-Space, for Linja-Karjala Oy network,pI = 0.05, pR = 0.0025, and T = 65, p-value of the Wilcoxon rank sum test (resultant p-value per Bonferroni correction is indicated in brackets).
ABDEB1.10e − 13 (1.10e − 12)D2.40e − 16 (2.40e − 15)0.11600 (1.00000)E1.50e − 15 (1.50e − 14)0.11000 (1.00000)0.59500 (1.00000)P<2.00e − 16 (3.20e − 16)0.65500 (1.00000)0.16400 (1.00000)0.09800 (0.98000)Table 15
In L-Space, for CIT network,pI = 0.05, pR = 0.0025, and T = 65, p-value of the Wilcoxon rank sum test (resultant p-value per Bonferroni correction is indicated in brackets).
ABDEB1.30e − 13 (1.30e − 12)D<2.00e − 16 (7.50e − 16)1.00000 (1.00000)E<2.00e − 16 (1.40e − 15)0.29000 (1.00000)0.20000 (1.00000)P<2.00e − 16 (7.30e − 16)0.63 (1.00000)0.56 (1.00000)0.41 (1.00000)Table 16
In L-Space, for Thunder Bay Transit network,pI = 0.05, pR = 0.0025, and T = 65, p-value of the Wilcoxon rank sum test (resultant p-value per Bonferroni correction is indicated in brackets).
ABDEB3.00e − 13 (3.00e − 12)D<2.00e − 16 (1.50e − 15)0.47900 (1.00000)E<2.00e − 16 (5.60e − 16)0.02700 (0.27000)0.08300 (0.83000)P<2.00e − 16 (9.30e − 16)0.16900 (1.00000)0.44800 (1.00000)0.27900 (1.00000)For STAR networks, in all studied duos (pI, pR) there was a similarity (p value > 0.05) between groups B-D, B-E, B-P, and D-P. In AVL, CFL, RGTR, and TICE network there is an equivalence between groups E-D and B-P.With respect to the TransAntofagasta network, a similarity exists between groups D-P. In Linja-Karjala Oy and CIF networks, there is a difference between groups A-B, A-D, A-E, and E-P. Additionally, dissimilarities exist between B-E groups in the Thunder Bay Transit network.In all analysed networks, according to the analysis carried out, the A category presents differences with respect to other groups, showing that those failures randomly happening on a node have a distinct behaviour to those originated by target attacks.In all networks, forpI = 0.05 and pR = 0.0025, Tables 17 and 18 show the time in which IR reaches 60%.Table 17
ForpI = 0.05 and pR = 0.0025, average propagation time in which IR reached 60% in each group (A, B, D, E, P) and globally.
ABDEPAllAuckland542.13386.85389.85408.83389.37415.24EMT261.31171.75162.83162.06195.33202.99BC Transit>700>700>700>700>700>700Kolumbus552.99458.23>700489.01>700>700STAR116.7593.5199.7495.0997.8899.67AVL, CFL, RGTR, and TICE121.3784.4290.4389.5883.7391.51TransAntofagasta>150145>150101.16>150105.57Linja-Karjala Oy127.2298.64102.54102.1297.50103.78CIT140.00115.65117.98119.56120.18120.01Thunder Bay Transit>150>150>150>150>150>150Table 18
ForpI = 0.05 and pR = 0.0025, median propagation time in which IR reached 60% in each group (A, B, D, E, P) and globally.
ABDEPAllAuckland545385390409388399EMT249.50169164161193181BC Transit>700>700>700>700>700>700Kolumbus547455>700486>700486.50STAR112.509498959999.67AVL, CFL, RGTR, and TICE116.0084899081.5088TransAntofagasta>150149>150100>150101Linja-Karjala Oy1289710310397102CIT138113117119120119Thunder Bay Transit>150>150>150>150>150>150It can be noted that in all analysed networks, random failures take the longest propagation time, demonstrating, as was previously covered in the section labelled study of structural properties, that these networks are less vulnerable to such failures. Depending on the network, some categories (D, E, P, or B) show similar propagation times. This seems to point out that in some cases, there is an existence of a certain amount of correlation between betweenness, degree, PageRank, and eigenvector centralities. Research shows that topology and density significantly affect the correlation between centrality measures and suggested that they could also affect the robustness [39].The networks that showed reduced failure propagation times were AVL, CFL, RGTR, and TICE; STAR; Linja-Karjala Oy; and TransAntofagasta. However, those that presented the highest times were in the order of BC Transit, Kolumbus, Auckland, EMT, and finally Thunder Bay Transit. In the magnitude of these times, there is an influence of the diameter presented by the network in L-space.


### 3.3. Analysis of Evolution of the GC

It is necessary to know the sensitivity of the network to the elimination of stops and links according to certain criteria in order to investigate which of them has a higher impact. For criteria 1–7, the nodes are removed in blocks of five elements and the percentage of the reduction in the GC is calculated. The applied criteria are the removal of random links (criterion 1) and the removal of random nodes (criterion 2). Other strategies such as the removal of nodes in the descending order of betweenness centrality (criterion 3), degree (criterion 4), eigenvector centrality (criterion 5), PageRank (with a damping factor equal to 0.85) (criterion 6), and highest weight links (criterion 7) are also applied. Using this criterion, a weighted graph was built, in which the weight of each link,wij, represented the portion of links taken from the total number of links between i and j. The interactive betweenness calculation (criterion 8) was also used. Figures 8 and 9 show, in all networks, the obtained results for the eight criteria.Figure 8
Evolution of the GC in EMT (a), BC Transit (b), Kolumbus (c), and Auckland (d) networks.
(a)(b)(c)(d)Figure 9
Evolution of the GC in AVL, CFL, RGTR, and TICE (a), STAR (b), TransAntofagasta (c), Linja-Karjala Oy (d), CIT Chambly-Richelieu-Carignan network (e), and Thunder Bay Transit (f).
(a)(b)(c)(d)(e)(f)In Table19 and Figures 8 and 9, it can be observed that in all networks, the four criteria that most quickly achieve 80% of the reduction of the GC after removing a number of nodes in the range [0%, 20%] have the highest interactive betweenness, PageRank, degree, and static betweenness centrality. And among them, in concordance with [6], in most networks, the elimination of nodes with the highest betweenness interactively calculated obtained an 80% reduction in the GC more rapidly, with a minor % of nodes eliminated. The criteria in which 80% was reached more slowly were those with a removal of random links, random nodes, the highest weighted links, and eigenvector centrality.Table 19
Percentage of removed nodes with which the 80% of reduction in the GC is achieved.
NetworkMethod% of removed elementsNetworkMethod% of removed elementsAucklandPagerank centrality4.98EMTInteractive betweenness centrality8.35Degree7.47Pagerank centrality12.62Interactive betweenness centrality8.67Degree13.70Betweenness centrality16.28Betweenness centrality24.38Random links28.79Random links47.02Weighted links56.85Eigenvector centrality51.02Random nodes61.99Random nodes52.42Eigenvector centrality64.66Weighted links60.22BC TransitPagerank centrality4.52KolumbusPagerank centrality8.89Degree5.79Degree12.29Interactive betweenness centrality9.55Interactive betweenness centrality15.04Betweenness centrality15.58Betweenness centrality19.74Random links30.19Random links43Eigenvector centrality59.55Weighted links49.75Weighted links68.94Random nodes64.17Random nodes86.50Eigenvector centrality73.86STARInteractive betweenness centrality4.53AVL, CFL, RGTR, TICEInteractive betweenness centrality3.57Pagerank14.18pagerank6.93Degree15.60Degree9.85Betweenness centrality17.73Betweenness centrality10.22Eigenvector centrality45.74Random links53.51Random links55.53Eigenvector centrality39.05Weighted links64.07Random nodes64.25Random nodes78.49Weighted links66.49TransAntofagastaInteractive betweenness centrality6.63Linja-Karjala OyInteractive betweenness centrality5.27Pagerank17.83Pagerank13.64Degree22.48Betweenness centrality17.27Betweenness centrality30.23Degree19.09Random links45.31Eigenvector centrality37.27Weighted links45.83Random links52.48Random nodes65.12Random nodes63.51Eigenvector centrality79.84Weighted links73.76CITInteractive betweenness centrality2.90Thunder Bay TransitInteractive betweenness centrality2.31Pagerank5.80Pagerank4.88Degree8.70Degree7.93Betweenness centrality13.04Betweenness centrality18.29Random links38.96Random links18.99Random nodes58.70Random nodes61.11Eigenvector centrality62.32Eigenvector centrality64.63Weighted links74.03Weighted links96.97The information above highlights the importance of taking special protective measures at some stops. The networks are more robust to failures in random links, random nodes, the highest weighted links, and eigenvector centrality than in the nodes with the highest interactive betweenness, PageRank, static betweenness centrality, and degree. However, when the faults occur randomly, those happening in links produce a higher impact. There is consistency with [40], which proves that links are more robust than nodes when the problems happen in highest degree nodes or in links, for instance, with a weight based on passengers’ flow between stops.Regardingk-core feature, it does not seem to influence the level of impact that a failure originating from a node will have on the network. This is in line with some research points, in reference to social networks, where in rumour spreading models, only the k-core of a node does not determine its propagation capacity [41].Additionally, if some of the topological parameters in L-Space, which are depicted in Table1, are seen in conjunction with the results shown in Tables 17 and 18, it can be noted that the networks with the highest value of density and lowest value of diameter are between those that obtain 60% in the IR more quickly.


## 3.1. Study of Structural Properties

The PageRank (with a damping parameter equal to 0.85), betweenness and eigenvector centralities, degree, and minimum distance characteristics were analysed in the L-Space and P-Space. Figures1 and 2 show the L-Space for all studied networks. The average values of degree, PageRank, eigenvector, and betweenness centralities for this space can be observed in Table 1. Those correspond to degrees, which are similar in all networks. The average eigenvector and PageRank centralities are lower than 0.002, having very small values in all networks. Therefore, the failures randomly happening in a node will have a low probability of arriving at any other. The same happens with respect to betweenness centrality, presenting small magnitudes (a mean value minor than 0.025), which denotes a slight average in how a node plays a bridging role in a network. Table 1 also shows the value of MRP in all networks, proving the existence in them of a GC (MRP >2).Figure 1
In L-Space, for EMT (a), BC Transit (b), Kolumbus (c), Auckland (d), and TransAntofagasta (e) networks.
(a)(b)(c)(d)(e)Figure 2
In L-Space, for Linja-Karjala Oy (a), AVL, CFL, RGTR, and TICE (b), STAR (c), CIT (d), and Thunder Bay Transit (e) networks.
(a)(b)(c)(d)(e)Table 1
Structural properties in L-Space. Average Betweenness Centrality (<BC>), average degree (<k >), average Eigenvector Centrality (<EC>), average PageRank (with a damping factor = 0.85) (<PR>), average path Length (<l>), diameter (d), density, and MRP.
NetworkNumber of nodesNumber of links<BC><k><EC><PR><l>dDensityMRPAuckland5,223318,8960.006552.360710.002170.000192.695231510.000452.66878EMT4,6362,070,5080.003832.711400.001810.000221.00000610.000583.21782BC Transit3,981702,1970.002332.288940.001670.000259.11427830.000582.54653Kolumbus3,829487,4320.002582.595460.002160.000262.89389980.000683.24314STAR1,415661,1500.006642.816960.005590.0007110.0000370.001993.54340AVL, CFL, RGTR, and TICE1,372340,6840.007062.810500.005640.000721.00000340.002053.79201Thunder Bay Transit82578,2470.034992.176970.005600.0012129.79920.002642.35412TransAntofagasta650724,3620.024362.963080.009020.0015416.78525510.004563.46729Linja-Karjala Oy55163,3990.016112.558980.014970.001814.50476290.004653.18865CIT3469,3660.027142.248550.017000.002892.33333290.006512.75804In L-Space, an equivalence between the degree distributions of the PTN of some cities was pointed out by Von Feber et al. [3], where the authors showed that this distribution followed a power law. This characteristic was checked for the analysed networks, fitting their degree distribution to the function:(5)Pk∼k−α.The obtained values for the fitting are shown in Table2. It is observed, according to the R square value, that those obtained for the fitting were not good (R squared is not close to 1). Figure 3 depicts, in L-Space, the degree distribution for the analysed networks. It is clear that all of them have a low value of k (kmax ≤ 17) compared with the total links. Therefore, the failures occurring at random on a node will have a high probability of having a small impact. Even so, these networks might be vulnerable to disruptions since a very high amount of nodes is joined to very few other nodes (low degree).Table 2
Fitting in L-Space of the network degree distribution to a power law functionPk∼k−α.
L-SpaceNetworkαR squareLinja-Karjala Oy0.1580.600Kolumbus0.4590.330EMT0.6220.239Auckland0.6030.237BC Transit0.4120.161AVL, CFL, RGTR, and TICE0.2360.156CIT0.3020.156Thunder Bay Transit0.4060.068TransAntofagasta0.1820.044STAR0.1580.041Figure 3
Degree distribution in linear (x-axis) and linear scale (y-axis) for the following networks: EMT (a), BC Transit (b), Kolumbus (c), Auckland (d), TransAntofagasta (e), Linja-Karjala Oy (f), AVL, CFL, RGTR, TICE (g), STAR (h), CIT (i), and Thunder Bay Transit (j).
(a)(b)(c)(d)(e)(f)(g)(h)(i)(j)It is clear from the above that the values of betweenness, eigenvector, and PageRank centralities in the analysed networks ensure that they are highly robust against random failures.The lowest average minimum distance between nodes shows a relevant variability in all networks. TransAntofagasta and EMT as well as AVL, CFL, RGTR, and TICE networks have the lowest value. The highest value corresponds to Thunder Bay Transit network.The modularity of the network is also evaluated in L-Space, employing the generate overlapping cluster generator (OCG) method [31]. It can be noted in Table 3 that all networks show a high value of Q with magnitudes near to 1, which demonstrates a strong community structure in these networks. Romano et al. [38] showed that network efficiency is modularity dependent, with the highest values of transmission occurring at intermediate levels and low values of modularity. High values tend to negatively influence transmission in networks [38]. Therefore, the high modularity in PTN could help to slow down the propagation of a failure in the network.Table 3
Q Modularity computed using the overlapping cluster generator (OCG) method in L-Space.
NetworkQBC Transit0.99750Auckland0.99610EMT0.98550Thunder Bay Transit0.99100Kolumbus0.98730STAR0.97020AVL, CFL, RGTR, and TICE0.96880CIT0.95880TransAntofagasta0.93910Linja-Karjala Oy0.93960Table4 depicts the values of the assortativity coefficient [32]. All networks show assortative higher than 0.4, demonstrating that a significant chance exists in which fractions of stops will randomly join with stations of the same degree. The nodes appear to join with a specific preference. A failure in a node with the highest degree will not have much negative impact, since other nodes with the highest degree are still connected to ones with a similar degree.Table 4
Assortativity coefficients in L-Space.
NetworkAssortativityKolumbus0.90599STAR0.89463BC Transit0.86067Auckland0.84123Linja-Karjala Oy0.81336AVL, CFL, RGTR, and TICE0.72469EMT0.67211TransAntofagasta0.58252Thunder Bay Transit0.45575CIT0.42243The highestk-core in all networks was 3 or 2, demonstrating that the largest subgraphs exist where every node is connected to at least 3 or 2 other nodes within the network. Most nodes show a small and similar k-core (≤2).With respect to P-Space, Table5 shows the average values of eigenvector and betweenness centralities, PageRank, degree, minimum distance between nodes, diameter, and density.Table 5
Structural properties in P-Space. Average Betweenness Centrality (<BC>), average degree (<k>), average Eigenvector Centrality (<EC>), average PageRank (with a damping factor = 0.85) (<PR>), average path Length (<l>), diameter (d), and density.
NetworkNumber of nodesNumber of links<BC><k><EC><PR><l>dDensityAuckland5,2238,079,0950.00039101.285080.023860.000192.6914370.01940EMT4,01025,263,1080.0006547.669830.046850.000253.5923870.01189BC Transit3,98115,737,7290.0001770.767140.034150.000252.6548160.01778Kolumbus3,82911,317,7140.0002878.983550.025840.000262.9688970.02063STAR1,4159,477,2130.0013743.112360.018620.000711.4827660.03049AVL, CFL, RGTR, and TICE1,3724,254,6560.0012048.425700.017910.000731.0000060.03532Thunder Bay Transit8252,014,7790.0013375.510300.021970.001212.0940.09164TransAntofagasta65031,171,3070.00108199.747690.036660.001541.7012330.30778Linja-Karjala Oy5511,167,6820.0018064.268600.035000.001811.4761940.11685CIT346135,7400.0027433.934780.034450.002891.0000040.103578The degree distributions of the analysed networks were also fitted to a power law function according to equation (5). The results are shown in Table 6. It can be noted that, as in the case of L-Space, in accordance with the obtained R square value, the fitting is not good. P-Space allows us to analyse interconnections between routes and transfers in the network. In this space, k reveals how many stops a person can travel through by using only one route. The network with the highest value average k, <k>, is TransAntofagasta. The obtained values for eigenvector and PageRank centralities are very low and similar in all networks, reinforcing the idea exposed in the analysis of L-space that failures randomly occurring in a node have a low probability of quickly arriving at any other one. Many possible transfer nodes exist, making the values for betweenness centrality very small and achieving high robustness against random failures on a node. The PTNs with lower <l> are the best connected in this space, having a minor level of vulnerability. Figure 4 shows, in P-Space, the degree distribution for the analysed networks.Table 6
Fitting in P-Space of the network degree distribution to a power law function.
NetworkαR squareEMT0.5130.180Kolumbus0.1690.166Auckland0.2040.077BC Transit0.1770.068Thunder Bay Transit0.2750.047TransAntofagasta0.1630.044STAR0.1180.037CIT0.1090.036AVL, CFL, RGTR, and TICE0.0610.022Linja-Karjala Oy0.0150.001Figure 4
Degree distributions in P-Space in linear (x-axis) and linear scale (y-axis) for the following networks: EMT (a), BC Transit (b), Kolumbus (c), Auckland (d), TransAntofagasta (e), Linja-Karjala Oy (f), AVL, CFL, RGTR, and TICE (g), STAR (h), CIT (i), and Thunder Bay Transit (j)
(a)(b)(c)(d)(e)(f)(g)(h)(i)(j)


## 3.2. Simulation of Faults Propagation

Failure propagation was simulated in L-Space using an infection and recovery algorithm.(1)
Firstly, a total of 100 experiments for each selection criterion and duple (pI, pR) were undertaken in each network. In each experiment, the simulation was implemented until t = T or until the percentage of infected nodes was higher or equal to 60%. The values of t were modified in steps of 10 units.(2)
The second step was to check the normality of distributions of the infection rate (IR), in all groups (A, B, D, E, P) using the D’Agostino test for anyT, which took different values until reaching a level propagation ≥60% of the network. The results show that p value was not higher than 0.05 in all cases. Since not all distributions presented normality, it was not necessary to study the homoscedasticity of variances applying the Breusch–Pagan Test. The Kruskal–Wallis test could be used for the detection of differences between groups.(3)
The third step was using the Kruskal–Wallis test to verify for eachT, whether all groups (A, B, D, E, P) were identical or by contrast, one of the groups tended to give observations that were different from those of other groups. Again, the p value obtained was lower than 0.05. The null hypothesis had to be rejected because at least one group had a different distribution.(4)
The fourth step was carried out utilising the Wilcoxon rank sum test in order to compare paired groups; we also applied the Bonferroni correction, both with a significance levelα = 0.05.In all analysed duos (pI, pR), for EMT, there was a significant difference (p value ≤ 0.05) between all groups except between groups D-E. With respect to BC Transit network, there were dissimilarities (p value ≤ 0.05) between some groups except between groups D-E, D-P, and E-P. Regarding the Kolumbus network, a relevant similarity is found between groups B-E and P-D. In the Auckland network, similarities happened between groups P-D, P-B, and D-B. For pI = 0.05 and pR = 0.0025, in T = 100, Figure 5 shows the results of the propagation algorithm in EMT, BC Transit, Kolumbus, and Auckland networks. For the same values of pI and pR and T = 65, Figure 6 shows the outcomes for STAR and AVL, CFL, RGTR, and TICE networks. Figure 7 depicts the results for TransAntofagasta, Linja-Karjala Oy, CIT, and Thunder Bay Transit networks. For those values of pI, pR, and T, Tables 7–16 show the obtained p value in the Wilcoxon rank sum test in each network.Figure 5
In L-Space for (a) EMT, (b) BC Transit, (c) Kolumbus, and (d) Auckland networks.pI = 0.05, pR = 0.0025, and T = 150 (boxplot diagram for A, B, D, E, and P groups).
(a)(b)(c)(d)Figure 6
In L-Space for (a) STAR and (b) AVL, CFL, RGTR, and TICE networks.pI = 0.05, pR = 0.0025, and T = 65 (boxplot diagram for A, B, D, E, and P groups).
(a)(b)Figure 7
In L-Space for (a) TransAntofagasta, (b) Linja-Karjala Oy, (c) CIT, and (d) Thunder Bay Transit.pI = 0.05, pR = 0.0025, and T = 65 (boxplot diagram for A, B, D, E, and P groups).
(a)(b)(c)(d)Table 7
In L-Space, for BC Transit network,pI = 0.05, pR = 0.0025, and T = 150, p-value of the Wilcoxon rank sum test (resultant p-value per Bonferroni correction is indicated in brackets).
ABDEB<5.70e − 11 (8.60e − 10)D1.10e − 14 (1.60e − 13)1.60e − 09 (2.30e − 08)E2.40e − 14 (3.50e − 13)3.70e − 09 (5.50e − 08)0.90000 (1.00000)P2.20e − 14 (3.30e − 13)1.30e − 08 (2.00e − 07)0.27000 (1.00000)0.34000 (1.00000)Table 8
In L-Space, for Kolumbus network,pI = 0.05, pR = 0.0025, and T = 150, p-value of the Wilcoxon rank sum test (resultant p-value per Bonferroni correction is indicated in brackets).
ABDEB<1.40e − 13 (2.00e − 12)D0.00473 (0.071)<2.00e − 16 (1.10e − 15)E3.00e − 14 (4.50e − 13)0.44174 (1.000)<2.00e − 16 (1.00e − 15)P0.00554 (0.083)<2.00e − 16 (1.10e − 15)0.99123 (1.000)<2.00e − 16 (1.00e − 15)Table 9
In L-Space, for EMT network,pI = 0.05, pR = 0.0025, and T = 150, p-value of in the Wilcoxon rank sum test (resultant p-value per Bonferroni correction is indicated in brackets).
ABDEB<2.00e − 16 (<2.00e − 16)D<2.00e − 16 (<2.00e − 16)0.00027 (0.0041)E<2.00e − 16 (<2.00e − 16)1.40e − 06 (2.00e − 05)0.30986 (1.00000)P<2.00e − 16 (<2.00e − 16)<2.00e − 16 (2.00e − 16)<2.00e − 16 (<2.00e − 16)<2.00e − 16 (<2.00e − 16)Table 10
In L-Space, for Auckland network,pI = 0.05, pR = 0.0025, and T = 150, p-value of the Wilcoxon rank sum test (resultant p-value per Bonferroni correction is indicated in brackets).
ABDEB<2.00e − 16 (<2.00e − 16)D<2.00e − 16 (<2.00e − 16)0.49000 (1.00000)E<2.00e − 16 (<2.00e − 16)7.40e − 14 (1.10e − 12)3.60e − 12 (5.30e − 11)P<2.00e − 16 (<2.00e − 16)0.54000 (1.00000)0.96000 (1.00000)1.00e − 12 (1.50e − 11)Table 11
In L-Space, for STAR network,pI = 0.05, pR = 0.0025, and T = 65, p-value of the Wilcoxon rank sum test (resultant p-value per Bonferroni correction is indicated in brackets).
ABDEB1.20e − 10 (1.20e − 09)D5.70e − 12 (5.70e − 11)0.056 (0.56)E2.90e − 14 (2.90e − 13)0.752 (1.00000)0.007 (0.07000)P1.20e − 12 (1.20e − 11)0.193 (1.00000)0.689 (1.00000)0.01900 (0.19000)Table 12
In L-Space, for AVL, CFL, RGTR, and TICE networks,pI = 0.05, pR = 0.0025, and T = 65, p-value of the Wilcoxon rank sum test (resultant p-value per Bonferroni correction is indicated in brackets).
ABDEB2.30e − 14 (2.3e − 13)D2.00e − 16 (4.8e − 16)1.70e − 07 (1.7e − 06)E< 2.00e − 16 (3.1e − 16)9.00e − 06 (9.0e − 05)0.19 (1.00000)P< 2.00e − 16 (<2e − 16)0.19 (1.00000)2.70e − 07 (2.7e − 06)4.30e − 05 (0.00043)Table 13
In L-Space, for TransAntofagasta network,pI = 0.05, pR = 0.0025, and T = 65, p-value of the Wilcoxon rank sum test (resultant p-value per Bonferroni correction is indicated in brackets).
ABDEB0.0013 (0.008)D1.80e − 11 (1.10e − 10)5.50e − 07 (3.30e − 06)E1.70e − 14 (1.70e − 13)2.90e − 14 (2.90e − 13)<2.00e − 16 (<2.00e − 16)P6.50e − 12 (3.90e − 11)1.50e − 07 (8.90e − 07)0.5348 (1.00000)<2.00e − 16 (<2.00e − 16)Table 14
In L-Space, for Linja-Karjala Oy network,pI = 0.05, pR = 0.0025, and T = 65, p-value of the Wilcoxon rank sum test (resultant p-value per Bonferroni correction is indicated in brackets).
ABDEB1.10e − 13 (1.10e − 12)D2.40e − 16 (2.40e − 15)0.11600 (1.00000)E1.50e − 15 (1.50e − 14)0.11000 (1.00000)0.59500 (1.00000)P<2.00e − 16 (3.20e − 16)0.65500 (1.00000)0.16400 (1.00000)0.09800 (0.98000)Table 15
In L-Space, for CIT network,pI = 0.05, pR = 0.0025, and T = 65, p-value of the Wilcoxon rank sum test (resultant p-value per Bonferroni correction is indicated in brackets).
ABDEB1.30e − 13 (1.30e − 12)D<2.00e − 16 (7.50e − 16)1.00000 (1.00000)E<2.00e − 16 (1.40e − 15)0.29000 (1.00000)0.20000 (1.00000)P<2.00e − 16 (7.30e − 16)0.63 (1.00000)0.56 (1.00000)0.41 (1.00000)Table 16
In L-Space, for Thunder Bay Transit network,pI = 0.05, pR = 0.0025, and T = 65, p-value of the Wilcoxon rank sum test (resultant p-value per Bonferroni correction is indicated in brackets).
ABDEB3.00e − 13 (3.00e − 12)D<2.00e − 16 (1.50e − 15)0.47900 (1.00000)E<2.00e − 16 (5.60e − 16)0.02700 (0.27000)0.08300 (0.83000)P<2.00e − 16 (9.30e − 16)0.16900 (1.00000)0.44800 (1.00000)0.27900 (1.00000)For STAR networks, in all studied duos (pI, pR) there was a similarity (p value > 0.05) between groups B-D, B-E, B-P, and D-P. In AVL, CFL, RGTR, and TICE network there is an equivalence between groups E-D and B-P.With respect to the TransAntofagasta network, a similarity exists between groups D-P. In Linja-Karjala Oy and CIF networks, there is a difference between groups A-B, A-D, A-E, and E-P. Additionally, dissimilarities exist between B-E groups in the Thunder Bay Transit network.In all analysed networks, according to the analysis carried out, the A category presents differences with respect to other groups, showing that those failures randomly happening on a node have a distinct behaviour to those originated by target attacks.In all networks, forpI = 0.05 and pR = 0.0025, Tables 17 and 18 show the time in which IR reaches 60%.Table 17
ForpI = 0.05 and pR = 0.0025, average propagation time in which IR reached 60% in each group (A, B, D, E, P) and globally.
ABDEPAllAuckland542.13386.85389.85408.83389.37415.24EMT261.31171.75162.83162.06195.33202.99BC Transit>700>700>700>700>700>700Kolumbus552.99458.23>700489.01>700>700STAR116.7593.5199.7495.0997.8899.67AVL, CFL, RGTR, and TICE121.3784.4290.4389.5883.7391.51TransAntofagasta>150145>150101.16>150105.57Linja-Karjala Oy127.2298.64102.54102.1297.50103.78CIT140.00115.65117.98119.56120.18120.01Thunder Bay Transit>150>150>150>150>150>150Table 18
ForpI = 0.05 and pR = 0.0025, median propagation time in which IR reached 60% in each group (A, B, D, E, P) and globally.
ABDEPAllAuckland545385390409388399EMT249.50169164161193181BC Transit>700>700>700>700>700>700Kolumbus547455>700486>700486.50STAR112.509498959999.67AVL, CFL, RGTR, and TICE116.0084899081.5088TransAntofagasta>150149>150100>150101Linja-Karjala Oy1289710310397102CIT138113117119120119Thunder Bay Transit>150>150>150>150>150>150It can be noted that in all analysed networks, random failures take the longest propagation time, demonstrating, as was previously covered in the section labelled study of structural properties, that these networks are less vulnerable to such failures. Depending on the network, some categories (D, E, P, or B) show similar propagation times. This seems to point out that in some cases, there is an existence of a certain amount of correlation between betweenness, degree, PageRank, and eigenvector centralities. Research shows that topology and density significantly affect the correlation between centrality measures and suggested that they could also affect the robustness [39].The networks that showed reduced failure propagation times were AVL, CFL, RGTR, and TICE; STAR; Linja-Karjala Oy; and TransAntofagasta. However, those that presented the highest times were in the order of BC Transit, Kolumbus, Auckland, EMT, and finally Thunder Bay Transit. In the magnitude of these times, there is an influence of the diameter presented by the network in L-space.


## 3.3. Analysis of Evolution of the GC

It is necessary to know the sensitivity of the network to the elimination of stops and links according to certain criteria in order to investigate which of them has a higher impact. For criteria 1–7, the nodes are removed in blocks of five elements and the percentage of the reduction in the GC is calculated. The applied criteria are the removal of random links (criterion 1) and the removal of random nodes (criterion 2). Other strategies such as the removal of nodes in the descending order of betweenness centrality (criterion 3), degree (criterion 4), eigenvector centrality (criterion 5), PageRank (with a damping factor equal to 0.85) (criterion 6), and highest weight links (criterion 7) are also applied. Using this criterion, a weighted graph was built, in which the weight of each link,wij, represented the portion of links taken from the total number of links between i and j. The interactive betweenness calculation (criterion 8) was also used. Figures 8 and 9 show, in all networks, the obtained results for the eight criteria.Figure 8
Evolution of the GC in EMT (a), BC Transit (b), Kolumbus (c), and Auckland (d) networks.
(a)(b)(c)(d)Figure 9
Evolution of the GC in AVL, CFL, RGTR, and TICE (a), STAR (b), TransAntofagasta (c), Linja-Karjala Oy (d), CIT Chambly-Richelieu-Carignan network (e), and Thunder Bay Transit (f).
(a)(b)(c)(d)(e)(f)In Table19 and Figures 8 and 9, it can be observed that in all networks, the four criteria that most quickly achieve 80% of the reduction of the GC after removing a number of nodes in the range [0%, 20%] have the highest interactive betweenness, PageRank, degree, and static betweenness centrality. And among them, in concordance with [6], in most networks, the elimination of nodes with the highest betweenness interactively calculated obtained an 80% reduction in the GC more rapidly, with a minor % of nodes eliminated. The criteria in which 80% was reached more slowly were those with a removal of random links, random nodes, the highest weighted links, and eigenvector centrality.Table 19
Percentage of removed nodes with which the 80% of reduction in the GC is achieved.
NetworkMethod% of removed elementsNetworkMethod% of removed elementsAucklandPagerank centrality4.98EMTInteractive betweenness centrality8.35Degree7.47Pagerank centrality12.62Interactive betweenness centrality8.67Degree13.70Betweenness centrality16.28Betweenness centrality24.38Random links28.79Random links47.02Weighted links56.85Eigenvector centrality51.02Random nodes61.99Random nodes52.42Eigenvector centrality64.66Weighted links60.22BC TransitPagerank centrality4.52KolumbusPagerank centrality8.89Degree5.79Degree12.29Interactive betweenness centrality9.55Interactive betweenness centrality15.04Betweenness centrality15.58Betweenness centrality19.74Random links30.19Random links43Eigenvector centrality59.55Weighted links49.75Weighted links68.94Random nodes64.17Random nodes86.50Eigenvector centrality73.86STARInteractive betweenness centrality4.53AVL, CFL, RGTR, TICEInteractive betweenness centrality3.57Pagerank14.18pagerank6.93Degree15.60Degree9.85Betweenness centrality17.73Betweenness centrality10.22Eigenvector centrality45.74Random links53.51Random links55.53Eigenvector centrality39.05Weighted links64.07Random nodes64.25Random nodes78.49Weighted links66.49TransAntofagastaInteractive betweenness centrality6.63Linja-Karjala OyInteractive betweenness centrality5.27Pagerank17.83Pagerank13.64Degree22.48Betweenness centrality17.27Betweenness centrality30.23Degree19.09Random links45.31Eigenvector centrality37.27Weighted links45.83Random links52.48Random nodes65.12Random nodes63.51Eigenvector centrality79.84Weighted links73.76CITInteractive betweenness centrality2.90Thunder Bay TransitInteractive betweenness centrality2.31Pagerank5.80Pagerank4.88Degree8.70Degree7.93Betweenness centrality13.04Betweenness centrality18.29Random links38.96Random links18.99Random nodes58.70Random nodes61.11Eigenvector centrality62.32Eigenvector centrality64.63Weighted links74.03Weighted links96.97The information above highlights the importance of taking special protective measures at some stops. The networks are more robust to failures in random links, random nodes, the highest weighted links, and eigenvector centrality than in the nodes with the highest interactive betweenness, PageRank, static betweenness centrality, and degree. However, when the faults occur randomly, those happening in links produce a higher impact. There is consistency with [40], which proves that links are more robust than nodes when the problems happen in highest degree nodes or in links, for instance, with a weight based on passengers’ flow between stops.Regardingk-core feature, it does not seem to influence the level of impact that a failure originating from a node will have on the network. This is in line with some research points, in reference to social networks, where in rumour spreading models, only the k-core of a node does not determine its propagation capacity [41].Additionally, if some of the topological parameters in L-Space, which are depicted in Table1, are seen in conjunction with the results shown in Tables 17 and 18, it can be noted that the networks with the highest value of density and lowest value of diameter are between those that obtain 60% in the IR more quickly.


## 4. Conclusions

An efficient PTN must be robust to provide a good quality service, even when it faces faults or disaster. This study analyses several PTNs, applying a methodology based on the study of structural properties, simulations of faults propagation, and the examination of the evolution of the GC. The analysis aims to detect the most vulnerable nodes and links when some faults occurred in them. The influence of the topological properties of the networks over their robustness has also been checked. This could allow preventative action in certain nodes and links against traffic congestion and unattainable stops, occurring as a result of construction work and catastrophes.Problems that happened at random nodes and over the strongest links spread slower than other types of failure. Those nodes whose failures spread more quickly depended on the network, but they were those with the highest betweenness centrality, degree, or PageRank. All networks show high modularity >0.8 and a Molloy–Reed parameter >2. Additionally, all nodes were characterized by lowk-core.The AVL, CFL, RGTR, and TICE networks were that presented the highest-level propagation in a timeT, for all analysed pairs (pI, pR), followed in descending order by STAR, Linja-Karjala Oy, and TransAntofagasta. In the majority of the networks, 80% reduction in the size of GC with the minor percentage of removed nodes is reached for those nodes with the highest betweenness calculated interactively. In a small number of networks, this was achieved with the removal of nodes with the highest PageRank. The networks that showed the highest amount of sensitivity were in sequential order of Thunder Bay Transit; CIT; and AVL, CFL, RGTR, and TICE, with 2.31%, 2.90%, and 3.57% of deleted nodes, respectively. In all networks, random failures took the longest time to spread across all networks, which proves, in conjunction with some structural properties (k, PageRank, and eigenvector in L-Space), that the networks are less vulnerable to such failures.This research can be continued by analysing passenger flow in the PTN as a dynamically changing characteristic. Some models that take into account the individual passengers’ behaviour, with respect to the choice of routes and connections during a trip, can be implemented in order to dynamically establish the possible consequences of the congestion. It is also possible to study how the passengers modified their travel behaviour when considering their previous experience. Additionally, the PTN vulnerability can be examined in other graph representations, as well as considering different link weights, such as trip time or financial costs. In addition to the above, in a catastrophic situation, the postdisaster circumstances can be examined, completing the analysis of the evolution of GC with a study of the access to critical facilities (i.e., emergency services).



---

*Source: 5513311-2021-03-11.xml*