AbstractDisruptions due to either natural or anthropogenic hazards can significantly impact the operation of critical infrastructure networks (e.g., transportation systems) as they may instigate network-level (cascade) systemic risks, thus impacting the overall city resilience. Recent relevant studies demonstrated the need to quantify the resilience of city infrastructure networks following failures of one/some of their main components, considering both topological and operational network measures. Subsequently, focusing on robustness (a key resilience attribute) and on transit (a major critical infrastructure network), the current study develops a related quantification tool employing a hybrid approach that integrates complex network theoretic measures with data analytics, and specifically clustering and genetic algorithms. To demonstrate the practical utility of the developed tool, the robustness of the City of Minneapolis bus transit network is quantified under possible cascade failures represented by node (i.e., bus stop), link (i.e., route segment), and route failure scenarios. The robustness quantification of this transit network is facilitated by analyzing 43 topological and operational measures using a coupled map lattice model integrated with a direction-based passenger flow redistribution model. Absorptive capacity thresholds are subsequently identified under different passenger flow-to-route capacity ratios. Finally, the routes are categorized based on their influences on the network robustness using genetic algorithms coupled with K-means clustering. The developed approach aims at providing a better understanding of transit systems pre- and postdisruptions by identifying key components that control the network robustness and subsequently devising reliable systemic risk management strategies and recovery plans. Such strategies and plans are expected to facilitate city resilience planning through management of cascade failure risks attributed to natural and anthropogenic hazard events.IntroductionThe city economy, security, public health, and safety are enabled by critical infrastructure networks (e.g., power, transportation, and water) that typically comprise multiple complex, interdependent, and dynamic (i.e., evolve with time) subsystems (Salama et al. 2020). Therefore, a disruption of a single subsystem (e.g., due to cyberattacks, natural hazards, operational incidents) may instigate network-level cascade risks, known as systemic risks (Ezzeldin and El-Dakhakhni 2019) that may hinder different human activities and cause significant economic losses. Among the different infrastructure networks, transit is the backbone of the city transportation system as it provides the basic mobility and accessibility services (Litman 2014), supports economy and urban expansion, induces transit-oriented developments, and plays a vital role in emergency responses and evacuation schemes during major incidents (Humphrey 2008).Among the different public transportation modes (e.g., light rail, subway, commuter rail, rapid transit), transit networks are key in serving dense regions of the city (e.g., the downtown core) and for extending land use (Martin et al. 1998; Metropolitan Council 2009). For instance, in the United States, within 2019, the Minneapolis Metro transit agency served an average of 249,300 trips during the weekdays, with the urban local bus network contributing approximately 70% of this total number of served trips (APTA 2019). As such, disruptions to bus networks can cause a significant reduction in the overall transit system performance, thus impacting the reliability and quality of the service provided. Examples of recent disruptions that significantly influenced different urban transit networks in North America include (1) the 2020 demonstrations across the United States (e.g., in Minneapolis, Chicago, Miami, and Los Angeles) due to the death of George Floyd, which resulted in rerouting several bus lines, canceling hundreds of trips, reducing the bus frequency, interrupting and suspending the transit service for several days, and stranding thousands of commuters (Fox9 2020); and (2) the shutdown of the Toronto subway in 2020 due to a train derailment, which triggered massive disruptions at several bus stops causing a chaotic situation for the transit network in Toronto for nearly 4 h during the morning peak time. Such disruptions, as well as others not listed herein, have significantly impacted key components of urban transit networks (e.g., stations and routes) either directly or indirectly, and thus transit service. This situation highlights the critical need to quantify and enhance the resilience of urban local bus transit networks to natural and anthropogenic hazard-induced disruptions.There is a lack of consensus on the definition of resilience within the transportation engineering field (Hosseini et al. 2016; Sun and Guan 2016; Haggag et al. 2020). However, according to the US Department of Transportation, resilience of transportation networks is defined as “the ability to anticipate, prepare for, adapt and absorb perturbation and withstand, respond to, and recover rapidly from disruptions” (FHWA 2012). Therefore, in transportation networks, resilience can be related to different operational measures, including hours of congestion, travel time index, traffic flow, and the optimal spare capacity (Hollnagel et al. 2008; Weilant et al. 2019). The spare capacity in turn can be further categorized into three different subcapacities (Hosseini and Barker 2016; Weilant et al. 2019): (1) absorptive capacity defined as the ability of a system to absorb and withstand the impacts and consequences of disruptive events; (2) adaptive capacity defined as the ability of a system to adapt in response to disruptions in order to maintain normal functions; and (3) restorative capacity defined as the ability of a system to quickly recover from and restore normal operations following disruptions (Norris et al. 2008; Turnquist and Vugrin 2013). Resilience can also be represented by its 4R attributes: robustness, rapidity, redundancy, and resourcefulness (Bruneau et al. 2003). Robustness measures the system’s ability to withstand the stress induced by sudden disruptions with minimal performance loss; redundancy is the ability of a system to meet the required functions during disruptions through alternatives to the out-of-service components; resourcefulness represents the capacity to mobilize resources to address priorities during disruptions; and, finally, rapidity reflects the timely rebound from the disrupted state to the (near) fully operational one. Note that the current study focuses only on the robustness measure of resilience in the context of transportation networks.For a transportation network, when a single component (e.g., bus stop, route segment) becomes nonoperational, network passengers typically resort to other operational components, which can instigate a cascade failure scenario if such components reach their capacities. Robustness quantification can thus facilitate identifying the most critical components that can cause network-level collapse (i.e., posing systemic risks). Exploring such critical components, prior to a disruption, can facilitate developing effective risk mitigation strategies (e.g., rescheduling transit services and/or providing alternative travel routes for passengers) that can improve the network robustness, and thus its resilience to cascade failures. Therefore, several studies have attempted to relate the robustness to different operational measures in transportation networks. For example, Zhang and Miller-Hooks (2015) used the average travel time to identify vulnerable links and evaluate the reserve capacity in rapid transit networks. Cats and Jenelius (2015), Hosseini and Barker (2016), and Weilant et al. (2019) employed different capacity metrics and investigated their impacts on the robustness of transit networks. Numerous indices have been also suggested in other related studies, including the giant component size (Lordan et al. 2014; Jia et al. 2021), the average shortest path length (Yang et al. 2015; Hossain and Alam 2017), the number of transfer stations (Wang et al. 2017), the delays and travel time (De Los Santos et al. 2012), and the percentage of cancelled trips (Lusby et al. 2018). However, other studies (e.g., Derrible and Kennedy 2010; Abdelaty et al. 2020; Yassien et al. 2020; Alzoor et al. 2021) highlighted the importance of using indices that consider both operational and topological characteristics of the network in order to achieve more realistic and practical robustness representation. This is because the interplay between such network characteristics has been proven to control cascade failures and accelerate recovery following natural and anthropogenic hazards (Strogatz 2001; Boccaletti et al. 2006; Ezzeldin and El-Dakhakhni 2019). In this respect, the current study employs the fraction of active nodes (i.e., operating bus stops) as a robustness measure because it considers both the topological (i.e., node degree and strength) and operational (i.e., redistribution of the passenger flow on links) characteristics of transit networks during disruptions (i.e., over time), which is key for resilience planning.Unlike simple graph theory (Guze 2019; Żochowska and Soczówka 2020) and GIS-based network analysis techniques (Inanloo et al. 2016; Yang et al. 2018), complex network theory (CNT) represents a computationally efficient alternative that provides a unique set of measures to simulate the aggregated temporal behavior of complex systems during disruptions (Derrible and Kennedy 2011; Barabási 2013). CNT is thus ideal for integration with network-specific simulation techniques for systemic risk and robustness quantification. Within CNT, networks are simulated as a set of nodes (i.e., representing the network components) connected by links (i.e., representing the interdependence between these components). CNT has been applied in earlier studies to investigate the robustness of a wide spectrum of transportation networks (including transit networks). Such studies evaluated the robustness of transportation networks either statically (e.g., Wang et al. 2011; Zou et al. 2013) or dynamically (Sullivan et al. 2010; Cats and Jenelius 2015; Huang et al. 2015; Nian et al. 2019).Static robust assessment of transportation networks relies on investigating the topological characteristics of such networks with the assumption that the disruption of a single component (e.g., a stop or a route) does not trigger a redistribution of node demand throughout the network. In general, most previous studies (e.g., Zou et al. 2013; Ma et al. 2020) highlighted the limitation of such an assumption and recommended applying other sophisticated models in order to achieve realistic and practical results. Dynamic robustness assessment of transportation networks, in contrast, considers not only the network topology but also the redistribution of the passengers following the failure of some network components. Therefore, dynamic robustness assessment is preferred as it can be employed to evaluate the impact of triggering cascade failures to other components beyond those that failed initially. Various models have been developed to quantify the dynamic robustness of transportation networks. The coupled map lattice (CML) models, first introduced by Kaneko (1992) to model the spatiotemporal chaos and pattern formation in fluids, were adapted to simulate transportation networks considering cascade failures (Chazottes and Fernandez 2005; Huang et al. 2015; Sun et al. 2018). The basic assumption of the CML model is that the spatiotemporal chaotic failure is attributed to both the topology of the network and the redistribution of the link weights. Specifically, in the CML models, cascade failure simulation considers (1) the state, degree, and strength of nonoperational stops and their neighboring stops; and (2) the weights and capacities of links connected to such stops.Research Gaps and Study Goal and ObjectivesDespite adopting CML models within the transportation engineering field, most related studies simulated their underlying transit networks using undirected networks. Such a simulation approach can result in misleading conclusions as the ridership (e.g., passenger flow) may vary over directions, and centrality measures are significantly different between directed and undirected networks (White and Borgatti 1994; Barabási 2013). Only a few studies have applied the CML model to quantify the robustness of directed transit networks. In addition, only a limited number of studies investigated the robustness of transit networks under route failure scenarios. For example, Shen et al. (2019) developed a modified directed version of the CML model to assess the robustness of the Nanjing metro transit network. However, the cascade failure conducted in this study has been limited to node failure scenarios only. Zhang et al. (2019a) also investigated the impact of route failures in weighted bus transit networks considering the dynamic load redistribution and link prediction method. Although these previous studies facilitated analyzing dynamic cascade failures, link weights and capacities as well as directed centrality measures were not considered for critical routes, necessitating the need for further investigations to assess the robustness of directed transit networks based on different failure scenarios (i.e., nodes, links, and routes) and relating the robustness thresholds to operational transit measures through a practical mapping approach.To address this research gap, the main objectives of the current study are to develop operationalizable tools to quantify the robustness of transit networks under systemic risks and subsequently identify the corresponding critical network routes through a novel hybrid approach that utilizes CNT, genetic algorithm (GA), and cluster analysis. The current study also investigates the influence of the absorptive capacity on the transit network robustness considering both its topological and operational characteristics under different disruption scenarios. The overarching goal of the study is to empower transit managers with the tools necessary for network resilience planning, considering the complexity of simulating interdependence-induced (cascade) disruptions.To demonstrate the utility of the developed tools, the Minneapolis local bus transit network is investigated, where CNT is employed to extract the topological measures of the network, and a CML model is developed considering these measures. To facilitate analyzing the network robustness considering the route directions, a direction-based passenger flow redistribution model is integrated with the developed CML model. The robustness of the network is then quantified for different passenger flow-to-route capacity ratios under node, link, and route failure scenarios, and absorptive capacity thresholds are accordingly obtained. The robustness of the network under route failures is evaluated considering 43 different topological (e.g., closeness centrality of bus stops) and operational (e.g., number of bus stops) measures, and GA is then coupled with K-means cluster analysis to identify and categorize critical routes based on key network measures. Specifically, the application of GA and K-means clustering facilitates identifying the specific few metrics to use in order to accurately characterize critical routes under disruptions. Empowered by the results of the aforementioned analyses, transportation service providers and city officials can devise proactive systemic risk mitigation strategies for transit networks (e.g., developing adaptive transit schedules and managing the transit fleet by increasing the spare capacity on specific routes) when subjected to disruptive events; this is key to enhancing the overall resilience of city infrastructure systems.Methods and MeasuresThis section describes the hybrid approach used in the developed tool to assess the robustness of bus transit networks. The analysis within this tool is divided into three stages, as shown in Fig. 1. In Stage 1, the underlying network is simulated using the CNT, where bus stops and route segments are represented by nodes and links, respectively, and the passenger flow and route capacities are then assigned to each link based on automated passenger count (APC) data. In Stage 2, the network classification and topological characteristics are evaluated using several CNT-based measures (i.e., node in-degree, out-degree, and strength). In Stage 3, the network robustness is assessed under systemic risks using a directed CML model that considers the redistribution of passenger flow due to the failure of the different network components (i.e., nodes and links). The inputs of the directed CML model in Stage 3 are (1) the transit network topology (from Stage 1); and (2) the node in-degree, out-degree, and strength, in addition to the passenger flow (from Stage 2). The network robustness is also investigated under different passenger flow to route capacity ratios and the corresponding absorptive capacity thresholds are identified. Finally, binary GA and K-means clustering are applied to categorize the network routes based on both their failure impacts on the network and their operational and topological characteristics.Complex Network Theoretic MeasuresDisruptions to bus stops/routes can cause significant changes in the network topological characteristics such as the degree centrality (ki), degree distribution (Pk), and node strength (Si). The degree centrality, ki, is defined as the total number of links connected directly to node i. In transit networks, the degree centrality represents the number of route segments originated from or directed towards a specific stop along a route, and can thus reflect the stop size (e.g., terminals or stations) within the network. In directed networks, this measure is calculated as the summation of the outward, kiout, and inward, kiin, links connected directly to node i, where kiout and kiin are calculated as follows (Opsahl et al. 2010; Barabási 2016): (1) (2) where aij and aji = elements of the adjacency matrix representing the network given that aij≠aji and aii=ajj=0; i is the start stop (source node); j is the end stop (all other nodes connected to the source node); and N = total number of nodes in the network. Note that aij is equal to 1 when nodes i and j are directly connected and is equal to zero elsewhere.The degree distribution, Pk, represents the percentage of nodes with a degree k, and is calculated in directed networks as the fraction of nodes with kout (Pkout) and the fraction of nodes with kin (Pkin) (Albert and Barabási 2002). In transit networks, the degree distribution represents the fraction of stops connected to a specific number of route segments and can then indicate the frequency of stop types (e.g., terminals or stations) within the network. Pkout and Pkin are evaluated based on the number of nodes with kout (Nkout) and the number of nodes with kin (Nkin), respectively (Albert and Barabási 2002) as (3) (4) The node strength, Si, represents the total passenger flow over the route segments that are connected to stop i (Sun et al. 2018), and can thus be used to reflect the passenger circulation (i.e., number of passengers transferred to/from a stop) at each stop along a route. Si is evaluated as (Sun et al. 2018) (5) Si=∑j≠iN wij+∑j≠iN wjinwhere n = number of nodes connected to node i; and wij and wji = sectional passenger flows (i.e., the number of passengers transferred by bus between stops i and j) on the links connected to node i. Such links are assigned directions based on the sequence of the stops along the routes. It should be highlighted when a transit network is represented by a directed network, Si can be partitioned into out-strength (Siout=∑j≠iN wij/n) and in-strength (Siin=∑j≠iN wji/n) to reflect the passenger circulation from and to stop i, respectively.Coupled Map Lattice ModelThe modified CML model developed by Huang et al. (2015) is applied in the current study to investigate the cascading failure within a weighted, directed urban local bus transit network. In this model, the passenger flow redistribution is modified to match the characteristics of a directed network, where the state of node i (xi) is updated over time (t) according to (6) xi(t+1)=|(1−ξ1−ξ2)f(xi(t))+ξ1∑j=1j≠iNaijf(xj(t))k(i)+ξ2∑j=1j≠iNwijf(xj(t))s(i)|+Rwhere f(xi(t)) = local dynamic chaotic behavior of the node i at time t; ξ1 = topological coupling coefficient; ξ2 = flow coupling coefficient; wij = hourly sectional passenger flow (Pf) between nodes i and j; and R = external perturbation severity that reflects the effect and magnitude of the failure that may occur. Larger R values correspond to failures with higher impacts, which seldom occur (Xu and Wang 2005; Huang et al. 2015; Sun et al. 2018). The ξ1 and ξ2 also typically ranges between 0 and 1 and their summation should be less than 1 (Huang et al. 2015; Sun et al. 2018).In the current study, the application of the modified CML model is started by assigning each node a random xi value between 0 and 1, where higher values reflect worse transportation conditions (Shen et al. 2019). A logistic chaotic map of [f(xi)=4xi(1−xi)] is then used to evaluate the corresponding dynamic behavior pattern of nodes (Sun et al. 2018; Shen et al. 2019). An R value of 1 is subsequently assigned to failed nodes to enable investigating the network robustness under failures with minimal impacts (Sun et al. 2018), and ξ1 and ξ2 values are assumed to be constant and equal to 0.25 (Huang et al. 2015). Note that if both xi and f(xi) are within the open interval] 0,1 [and R=0, the network keeps a permanent normal state. In contrast, if a node/link fails at t=l, additional passenger flow is produced and redistributed over downstream connections. A cascade failure scenario may consequently occur if the redistributed passenger flow results in violating the route capacity (RC).Failure ScenariosIn the current study, the network robustness is quantified when its main components become nonoperational due to different failure scenarios. Specifically, the study initially investigates the network robustness under the random failure of a single component (i.e., a node or link), which occurs frequently due to minor road accidents or severe weather conditions. The network robustness is also evaluated under random route failure scenarios, which may occur due to major incidents such as floods or snowstorms.Node and Link FailuresIn the current study, a node (i.e., stop) fails when the passengers cannot be transferred by bus from (or to) this stop due to any disruption. In such cases, buses served by failed nodes (i.e., nonoperational stops) would deviate to adjacent streets to continue serving passengers along the intended routes. However, the details of these rerouting procedures are not considered in this study due to data restrictions. Instead, when node i fails at time t, this node and its associated links are removed from the network during the subsequent time steps. As shown in Fig. 2(a), the weight of the removed links is then redistributed over the downstream links (i.e., located within a walkable distance) using Eq. (7) (7) The in-degree, out-degree, and strength of the network nodes are subsequently updated according to Eqs. (3)–(5), respectively. Using such node parameters together with the updated Pf values, the state of all the nodes is re-evaluated using the CML model presented in Eq. (6). Additional links may thus fail due to the weight redistribution process when their updated Pf/RC ratios are larger than 1.0, and they are in turn removed from the network. Subsequent rounds of redistribution are applied that can trigger cascade failures. Note that both the network topology and link weights are updated during the same time step and prior to the following ones. The pseudocode for the node failure scenario is presented in Algorithm 1.A link failure occurs when passengers cannot be transferred by bus between stops (nodes) i and j. Similar to the node failure scenario, buses served by failed links are assumed to move through adjacent streets to continue serving passengers along the intended routes. However, as mentioned earlier, the details of the rerouting procedures are not considered in the current study due to data restrictions. When a link fails at time t, its weight (Pf) is transferred to the destination node and subsequently distributed over the connected active (operating) links using Eq. (7), as shown in Fig. 2(b). Similar to the node failure scenario, the degree and strength of the network nodes are updated according to Eqs. (3)–(5), and are subsequently used together with the new Pf values (i.e., updated link weights) to calculate the temporal node state (i.e., xi) based on Eq. (6). The weight redistribution process continues (cascade failure occurs) when the value of xi of any node exceeds 1.0 or when the Pf/RC ratio of any link is larger than 1.0. The pseudocode for the link failure scenario is presented in Algorithm 2.Under random node and link failure scenarios, a node/link is selected randomly to initiate the failure, and the fraction of active nodes (RN) (i.e., the fraction of nodes with xi<1) is used to represent the network robustness. A single value of RN is then estimated for each failure scenario (i.e., considering the cascade effect). As the network robustness typically decreases with the increasing number of inactive nodes/links, critical thresholds, fc-N and fc-L, are defined to indicate the fraction of inactive nodes and links, respectively, that minimize RN. Note that the values of fc-N and fc-L depend on the average Pf/RC across the network. A sensitivity analysis is therefore essential to investigate the impact of Pf/RC variability (which represents the disruptions to bus frequency or demand over the network) on the values of fc-N and fc-L as well as the network robustness. A network absorptive capacity can be subsequently quantified as the average Pf/RC ratio (as suggested by Weilant et al. 2019) at which the network can withstand minor changes in fc-N and fc-L values under node and link failure scenarios, respectively.Route FailureRoute failures are simulated in the current study based on two distinct scenarios: (1) all stops along the route are assumed to be completely nonfunctional, or (2) all segments along the route become nonoperational. In the first scenario, stops on the failed routes are completely removed from the network. Accordingly, links (i.e., bus route segments) directed inward to or outward from such stops are also removed from the network. Subsequently, the network robustness is evaluated based on both the percentage of failed nodes (Nf-N%) and the route failure impact ratio (RFIR), where the latter represents the ratio between the number of failed nodes (NFR) along the route and the total number of nodes (NTR) along the same route. In the second scenario, links along the failed routes are completely removed from the network, albeit without removing stops connecting such links. The network robustness is subsequently quantified using the percentage of failed links (Nf-L%). Note that during both scenarios a bus route cannot serve a specific bus line and is thus fully disconnected from other routes.Genetic Algorithm and Cluster AnalysesCluster analysis is typically used to group observations, in an unsupervised manner, based on their shared similarities. The K-means clustering is widely used for such a purpose (Likas et al. 2003; Zhang et al. 2019b); however, it requires a prior definition of the number of clusters (K). A range of K values is therefore assumed when K-means clustering is applied, and the performances of the resulting models are assessed. The percentage of variance explained is the typical measure of the model performance and can be implicitly expressed through several alternatives (e.g., within cluster sum of squares, silhouette score) (Pandit and Gupta 2011).In the current study, K-means clustering was applied to categorize transit routes based on their criticality (i.e., their impacts on the network robustness under disruptions) and their main characterizing attributes are identified accordingly. Due to the large number of measures specific to the network topology, operation, and robustness, a feature selection approach should be applied. The current study thus integrates K-means clustering with GA to identify the most significant important measures for categorizing network routes. The GA is used to identify the optimum set of measures for each predefined cluster model and the process is repeated iteratively (i.e., by employing different K values) to identify the optimal variable set that maximizes the percentage of variance explained (Sexp2). The K values were thus varied between 2 and 30, with the upper limit selected to ensure the homogeneity of the resulting clusters (i.e., the clusters include a similar number of observations). The GA-based clustering can be formulated as (8) maxziSexp2subject to: zi is binary∀  iwhere Sexp2=1−SSc/SST; SSc = within cluster sum of squared distance; and SST = sum of squared distance between observation pairs (Jones and Harris 1999). As the optimum variable set may differ based on the K value, the frequency of including a specific variable in the optimum set is used to reflect its importance for clustering the network routes. Subsequently, the optimum value of K and the optimum variable set were chosen based on assessing the performance of the optimal solution over the different K values. It should be highlighted that, when GA is used to select the optimal features for clustering, decision variables of the optimization problem are represented by a binary vector where a value of 1 is used to indicate that a specific feature is included for clustering, while a value of 0 is used to indicate that the feature is not included. For example, a vector of the decision variables in the form of [1 0 1 0 0 0 0 … 0 1] indicates that the first, third, and last features will be included for clustering, while all other features will not be included.Demonstration ApplicationIn the current study, the Minneapolis local bus transit network was utilized as a demonstration for the application of the developed robustness assessment tool. Minneapolis is the largest and most populous city in the state of Minnesota (US Census Bureau 2019). The transportation system in the city is characterized by a highly connected urban multimodal transit network that serves numerous attractions, destinations, and other critical facilities such as universities, hospitals, public schools, and recreational and shopping centers (Minneapolis City Council 2019). The transit network serves an average of 249,300 passengers per day, with the local bus transit networks contain 27 routes [Fig. 3(a)] and serve around 68% of the daily demands (Metropolitan Council 2019).Data Description and ProcessingThe automated passenger count (APC) dataset of the network is provided by the Metro Transit Agency under the Metropolitan Council Agreement (Minnesota Geospatial Information Office 2019). The APC dataset is filtered by excluding the days of extreme weather, unusual events, and abnormal travel patterns. Note that, within the APC dataset, the total boarding and alighting are unequal for different bus stops, routes, and provider levels. Therefore, the boarding and alighting volumes were balanced in the current study following the method described by Lu (2008) and Boyle (2009) through two stages: (1) the variation between the total boarding (ons) and the total alighting (offs) along the routes has been calculated (i.e., between 0.8 and 1.1) and compared to the allowable thresholds (i.e., between 0.7 and 1.3) according to the TCRP synthesis 77 guidelines (Boyle 2009); and (2) ons and offs are balanced based on factoring the average method such that negative loading errors of the cumulative passenger flow are avoided (i.e., correcting negative loads). Note that other methods can be used to balance boarding-alighting errors such as factor into the higher count (Lan 2015) and pseudostop (Furth et al. 2006). However, factoring the average method was utilized in the current study as Lu (2008) highlighted that the three methods typically yield similar balanced trips. The transit time-of-day analysis is then performed to convert the total average daily passenger flow into a peak hourly flow, which was subsequently used as link weights. Because the peak hour flow determines the required size of the transit fleet, the transit time-of-day analysis was performed based on a conversion factor of 0.14 [according to the TRB-National Highway research program (Martin et al. 1998)] in order to convert the total average daily passenger flow (i.e., available within the dataset) to peak hour flow for all purpose transit trips.Network DevelopmentThe topological configuration of transit networks can be defined through a space of changes (Kurant and Thiran 2006), a space of routes (Xu et al. 2007; Zhang et al. 2019b), or a space of stops (Kurant and Thiran 2006; Xu et al. 2007). The current study adopts the latter configuration as it captures the passenger flow variability over route segments as well as the operating frequency and capacity of bus transit routes (Zhang et al. 2018). In this configuration, bus stops are represented by nodes, while links simulate either the physical or operational connections between these stops. Accordingly, the Minneapolis local bus transit network is conceptualized as a weighted, directed network of 4,339 nodes and 4,664 links, as shown in Fig. 3(b).Result AnalysesThis section demonstrates the topological characteristics and cascade failure analysis results of the network, with the topological characteristics of the network are initially examined based on various centrality measures. The CML model is subsequently used to obtain the dynamic robustness of the network based on random node, link, and route failure scenarios that represent different disruptions to bus transit components, as discussed earlier. Such disruptions may induce additional demands, thus impacting transit stops, roadway sections, or transit corridors (Bureau of Transportation Statistics 2018). The network robustness towards route failures is quantified considering various centrality and operational measures. A sensitivity analysis is also conducted to evaluate the impact of Pf/RC ratio variations on the failure thresholds, and the corresponding absorptive capacities are obtained.Topological CharacteristicsThe static analysis of the network topology demonstrates that approximately 81% of the nodes have kout=1 and kin=1, as shown in Fig. 4. As can be seen also in the figure, the node degree distribution of the network was found to follow a power-law distribution (Pkout=0.725kout−3.219 and Pkin=0.72kin−3.22), which supports classifying it as a scale-free network (Barabási 2013). However, the power-law exponent values (i.e., 3.219 and 3.22) highlight that the network is of significant heterogeneity, may act as a random network under random failures (Barabási 2013), and may therefore fragment at a finite critical threshold (i.e., fc-N or fc-L) independent of the network size, as will be discussed next.Node and Link FailuresFig. 5 shows the robustness analysis results for the network under both node (i.e., nodes and their connected links are completely removed from the network) and link (i.e., links are completely removed from the network) failure scenarios. The results showed that the network experienced different fragmentation schemes under both failure scenarios. For example, at a fraction of removed nodes or links, f, of 0.3, the corresponding network robustness values are 0.38 and 0.61, for node and link failures, respectively, as shown in Fig. 5. This observation indicates that the influence of node failures on the network robustness is approximately double that of link failures. This significant difference in the network behavior is attributed to the fact that when a node fails (i.e., a bus stop becomes nonoperational), this node is removed from the network along with all of its connected links, where the number of such links can be represented by the network average (typically larger than 1.0). Therefore, when an f fraction of nodes is removed, an f fraction of links is removed from the network, which results in more damage to the network than the removal of an f fraction of links. However, Fig. 5 shows that both scenarios yield completely disconnected networks at similar critical threshold values of fc-N(act)=0.58 and fc-L(act)=0.63 when nodes and links are triggered, respectively. By comparing such values to the theoretical critical threshold value, fc(th), developed by Molloy et al. (1995) for scale-free networks and presented in Eq. (9), there are variations of 20% and 15% for node and link failures, respectively. These variations are mainly because fc(th) is based only on the network topology [i.e., the minimum degree of nodes (kmin) and the degree exponent of the network power-law distribution (γ)] (Molloy and Reed 1995; Barabási 2013), whereas fc-N(act) and fc-L(act) values presented in the current study depend on the developed CML model that considers both the topology and the passenger flow when systemic risks within bus transit networks are evaluated (9) fc(th)=1−1(|γ|−2)(|γ|−3)kmin−1,for  |γ|>3Capacity DisruptionsThe passenger flow to route capacity (Pf/RC) ratios are utilized in the current study to represent any disruptions that may occur due to either variation in the frequency of buses or excess demands with respect to the route capacities. In this respect, disruptions to the network capacity have been investigated and the corresponding network robustness is quantified under various Pf/RC ratios. The absorptive capacity thresholds of the network are also obtained for both node and link failure scenarios. Figs. 6(a and b) show the sensitivity of fc-N and fc-L to different Pf/RC ratios under node and link failure scenarios, respectively. As can be seen in both figures, the critical threshold values remain almost the same at initial ratios of Pf/RC (i.e., less than 0.35 and 0.45 for nodes and links, respectively); however, the threshold values are reduced significantly when the Pf/RC exceeds these initial ratios. Specifically, the network exhibits an initial Pf/RC ratio of 0.18 and corresponding critical threshold value of fc-N(act)=0.58 and fc-L(act)=0.63 for node and link failures, respectively. The absorptive capacity is defined in the current study as the Pf/RC ratio that the network can withstand with minor changes in the corresponding critical thresholds. Therefore, the absorptive capacity of the network is identified as 0.35 (i.e., 200% increase) and 0.45 (i.e., 250% increase) under node and link failures, with corresponding minor variations in fc-N(act) and fc-L(act), respectively. These results highlight the significant impact of node failures on transit network capacities compared to that of link failures. As can be also seen in Figs. 6(a and b), for all Pf/RC ratios, the fc-N(act) and fc-L(act) are lower than fc(th). For instance, at Pf/RC=0.09 (i.e., 50% less than the original value), the fc(act) values are 0.59 and 0.64 for node and link failures, respectively, which are smaller than fc(th). This finding emphasizes the importance of considering the Pf/RC ratio as well as the topological characteristics of the network in obtaining the critical absorptive capacity thresholds for transit networks.Route FailureAs every route encompasses a set of nodes and links, the robustness of transit networks under route failures highly depends on the characteristics of such components. Therefore, under each route failure scenario, 43 different variables are calculated to represent a wide range of operational and topological measures, and these are listed in “Notation” section. A preliminary analysis is conducted to investigate the correlation between such measures and the route failure indices (Nf-N%, Nf-L, and RFIR). This correlation analysis was investigated based on the value of the Pearson’s correlation coefficient as seen in Fig. 7. As can be seen in Fig. 7, no significant correlation is apparent between the different robustness indices and most of the topological and operational measures. For example, the RFIR has a strong direct correlation with only two topological measures (DRavg and PRRavg, with correlation coefficient value>0.7), and a strong inverse correlation with only two operational measures (NL and NN, with correlation coefficient value<−0.7). Conversely, the Nf-L%, and Nf-N% have a strong direct correlation with two operational measures only (NL and NN, with correlation coefficient value>0.7). These results highlight the notion of indices/measure interdependence, rather than their simple correlation, which supports that the robustness of the network is a function of the interplay between all topological and operational measures. The integrated GA-clustering analysis was subsequently applied to categorize the network routes considering the interplay between their topology and operational measures and the robustness indices. The GA was applied using a population size of 200, a crossover ratio of 0.2, maximum generations of 100, and a tournament selection with a size equals to 2. The integrated GA-clustering analysis was applied 29 times, each of which includes: (1) predefining a K value (between 2 and 30); and (2) applying the GA to identify the optimum variable set that can represent the whole dataset through that number of clusters (i.e., K). An optimum variable set is thus identified for each K value and is referred to as the GA solution. Fig. 8 shows the frequency of including each of the 43 measures in the GA solution. A high frequency indicates that the corresponding measure was frequently used for clustering, which reflects its importance for representing the whole dataset. For example, the top five measures of the highest importance are the maximum weighted degree, SRmax, RFIR, average degree, DRavg, average Eigen-centrality, ECRavg, and maximum weighted outdegree, SRout-avg. These measures appeared 21, 18, 17, 14, and 14 times in the GA solutions, respectively. Measures of the highest importance include topology and centrality measures, which support the limitation of characterizing bus transit routes based on operational measures only and highlight the need for including topology measures to assess the overall bus transit network robustness under route failures.The percentage of the variance explained by each GA solution increases with the number of clusters (i.e., K), as shown in Fig. 9, where the best cluster model is that with six variables and 27 clusters (describing approximately 100% of the variance). However, this is considered an infeasible model due to the large number of clusters required. As such, a model with a smaller number of clusters can be more practical. The use of six clusters with eight variables can represent approximately 96% of the variance, which can be increased by only 4% using a model with more clusters. Accordingly, a model with six clusters and eight variables is considered herein as the near-optimal cluster model. The definitions, physical representations, and mathematical formulations of the selected eight variables are presented in Table 1, while Fig. 10 shows the statistical behavior of each of these variables over the six clusters.Table 1. Selected variables measures and definitionsTable 1. Selected variables measures and definitionsName and variableDefinition and representationFormulaRoute failure impact ratio (RFIR)The ratio between the total number of failed nodes (NFR) when a specific route failed and the number of nodes along that route (NTR) and is used to identify the critical routes that can facilitate the cascading failure in the network.RFIR=NFR/NTRAverage betweenness centrality (BCR)The average fraction of shortest paths passing through each node k along the route R and is used to identify the routes containing the relatively important stops in terms of location.BCR=1NTR∑k∑l∑mρ(l,k,m)ρ(l,m)where k≠l≠mAverage closeness centrality (CCR)The distance between each node k on the route R and all other nodes in the network (Opsahl et al. 2010; Barabási 2013) the degree of accessibility to the stops belonging to the route R compared to other stops in the network.CCR=1NTR∑kN−1∑jdkjwhere dkj = length of the shortest path between nodes k and j.Average Eigen centrality (ECR)The average influence of each node k along the route R based on the importance of the neighboring nodes.ECR=1NTR∑k1λ∑jAkjDECjwhere DEC = right leading eigenvector; Akj = adjacency matrix; and λ = eigenvalue of Akj for which DEC exists.Maximum in-degree (DRin)The maximum number of links directed toward any node k over the route R and corresponds to the stops with the largest size (major stop like terminals, or stations) that exist along the route R based on the number of routes directed toward the stop.DRin=maxk∑jajkwhere ajk = adjacency matrix entry corresponding to the connection between nodes j and k.Maximum modularity class (QR)The maximum modularity class (Q) amongst the nodes located along the route R, where Q is used to measure the likelihood of dividing the network into modules, where higher modularity denotes that a group of edges are most likely located within the same community.Q=1L∑i·j[aij−DiinDjoutL]δ(ci,cj)where aij = adjacency matrix; L = number of links in the network; Diin = indegree of the node i; Djout = outdegree of the node j;ci (cj) = community to which the node i (j) belongs; and δij = Kronecker delta function.Maximum weighted out-degree (SRout)The maximum outward node strength along the route R and represents the maximum passenger flow (Pf) originated from that route.SRout=maxk∑jwkjwhere wkj = entry of the weighted adjacency matrix corresponding to the connection between nodes k and jMaximum weighted degree (SR)The maximum node strength along the route R and corresponds to the stop with the highest passenger circulation along the route.SR=maxk(∑jwkj+∑jwjk)Although the routes can be classified based on their levels of operability, this classification may be misleading as the failure of a route can cause the failure of several nodes/links based on the connectedness of this route to other network components. As such, the network routes have been classified based on their RFIR (as an operability/robustness indicator) together with the BCRavg, CCRavg, ECRavg, DRin-max, and QRmax (as connectedness indicators) into low (i.e., low robustness/low connectedness), moderate (i.e., moderate robustness/intermediate connectedness), and high (i.e., high robustness/high connectedness). It should be emphasized that such classification procedures can cause the intervals of a single measure to overlap with those of other classes. Specifically, Clusters 1 to 4 include routes with low RFIR values that span between 0.18 and 0.48. Cluster 5 represents routes with high RFIR values that are between 0.62 and 0.99, while Cluster 6 contains moderate impact routes with RFIR values that vary between 0.34 and 0.63. The results show that BCRavg and ECRavg have a significant impact on the severity of the failure compared to other connectedness indicators, thus highlighting the importance of the location, accessibility levels, and characteristics of the neighboring nodes in controlling the severity of failure propagation through the network. Moreover, routes with relatively high connectedness indicators and moderate to high passenger flow indicators can significantly impact the overall network vulnerability. In contrast, routes with relatively low connectedness indicators and high passenger flow indicators are characterized by low or moderate route failure impact ratio (i.e., their failures have low to moderate effects on the overall network vulnerability). Fig. 11 shows the spatial distribution of the six clusters over the network. Most of the routs have a low impact on the overall network robustness; however, routes of relatively higher failure impacts (i.e., those with moderate to high RFIR values) are highly connected and located in the core of the network. Most of the bus stops along such routes exhibit high connectedness (expressed through higher BCRavg, CCRavg and ECRavg) compared to other routes. Conversely, routes of lower failure impacts (i.e., those with low RFIR values) are located near the network boundaries. However, a small fraction of the bus stops along these routes is highly connected to other routes in the network.ConclusionsTransportation infrastructure networks (e.g., urban transit) provide fundamental accessibility and mobility services for users, and thus any disruptions to such complex and dynamic networks are associated with severe socioeconomic impacts. Subsequently, research efforts focusing on quantifying and enhancing their robustness under cascade failure scenarios are critically needed in order to ensure reliable transit service and transportation network resilience under systemic risks posed by cascade disruptions. In the current study, a novel analysis approach is developed to quantify the robustness of urban bus transit networks under systemic risks, where complex network theory (CNT), genetic algorithm (GA), and K-means clustering are integrated in a hybrid simulation approach. The approach aims at utilizing several topological and operational measures to evaluate the robustness of transit networks following the failure of their main components (i.e., stops, links, and routes). In this respect, a coupled map lattice (CML) model is developed considering several CNT-based topological measures and employing a direction-based passenger flow redistribution approach to facilitate evaluating the robustness of the transit networks under different passenger flow-to-route capacity (Pf/RC) scenarios. GA is subsequently coupled with K-means clustering to identify the primary measures influencing the network robustness as well as to categorize the network routes based on their topological and operational measures. The City of Minneapolis bus transit network is considered to demonstrate the utility of the developed approach.The results of the network robustness analysis showed that node (i.e., stop) failures can trigger larger fragmentation of the transit network compared to those due to link (i.e., bus route segments) failures. The results also demonstrated that currently developed theoretical robustness thresholds are not representative of bus transit networks as they are determined based on static topological characteristics only. The capacity disruption analysis highlighted that absorptive capacity thresholds of 0.35 and 0.45 can be used for node and link failures, respectively, suggesting that the considered bus transit network’s spare capacity enables absorbing an increase in passenger flow of up to 200% and 250% under the two failure scenarios, respectively. In addition, the results demonstrated the strong inverse correlation between the absorptive capacity and the network robustness, where the network can exhibit a remarkable robustness reduction when exceeding its absorptive capacity thresholds.The application of GA and K-means cluster analyses identified the set of measures that can be used to characterize routes in transit networks. These measures include a robustness indicator (i.e., route failure impact ratio), four connectedness indicators (i.e., the average betweenness centrality, closeness centrality, Eigen centrality, and maximum modularity class), and three passenger flow indicators (i.e., the maximum weighted in-degree, out-degree, and total degree). This finding supports the importance of integrating both topological and operational measures when the robustness of transit networks is evaluated under systemic risks. The results also showed that both the average betweenness centrality and the Eigen centrality have a significant impact on the severity of cascade failures in transit networks compared to other connectedness indicators. This highlights the importance of location, accessibility levels, and characteristics of the neighboring nodes in controlling the failure propagation throughout the network. In other words, the characteristics of the nonoperational stops along with their neighboring stops are key factors that can either amplify or hinder cascade failures within transit networks.The analysis results highlighted the importance of modeling bus transit networks as directed ones. This is because the different directions of the same route can have different failure impacts on the network and subsequently, they may not fall within the same cluster as demonstrated in the current study. In addition, the current study showed that routes with relatively high connectedness indicators and moderate-to-high passenger flow indicators can significantly impact the overall network vulnerability. In contrast, routes with relatively low connectedness indicators and high passenger flow indicators are characterized by low or moderate route failure impact ratios, which supports the importance of including both topological and operational measures to identify critical routes in transit networks.In general, the hybrid simulation approach developed in the current study facilitates understanding the behavior of transit networks under different component (i.e., stops, links, and bus routes) failure scenarios. Transit service providers and city officials can employ the developed approach to (1) evaluate the robustness of different transit network structures under what-if scenarios and subsequently select the optimal configuration; (2) enhance the transit network capacity through calculating the absorptive capacity under various operational conditions; (3) devise capacity-based recovery plans (i.e., rescheduling the transit service, increasing bus frequency, reducing headways) to manage the risk of cascade failures due to passenger flow redistribution; and (4) classify the network routes according to their contribution to the network robustness and subsequently prepare effective risk mitigation strategies (e.g., providing alternative routes for travelers and/or managing the transit fleet). The presented model empowers transit service providers and city officials with the necessary simulation tools for effective city systems resilience planning and management under different forms of interdependence-induced network-level disruptions.Algorithm 1. Robustness of the Minneapolis Local Bus Transit Network under Random Node FailureAlgorithm 1. Robustness of the Minneapolis Local Bus Transit Network under Random Node FailureInput: initial random xi(t)for all υi, ξ1, ξ2, wij, aij, k(i), s(i), f(xi(t))Output:RN1: Procedure Random node failure in the Minneapolis local bus transit network2:  for all υi∈Vido3:      ifxi>1 or k(i)=04:       xi=0, ki=05:     else6:       xi=xi and k(i)=k(i)7:     end if8:  υk← Randomly select Node k9:  xk(t)=xk(t)+R10:  f(xk)=4xk(1−xk);11:   for all υk∈Viwhereaik=112:    do Eq. (6)13:   end for14:  ifxi(t+1)<115:   ki=kiout+kiin16:   do Eq. (5)17:   do Eq. (7)18:   f(xi(t+1))=4xi(1−xi)19:  else20:   xi(t+1)=xi(t+1)+R21:    for all υi∈Viwhereail=122:     do Eq. (6)23:    end for24:  end if25:    V←RN=N−Nf/N;26: end for27: end procedureAlgorithm 2. Robustness of the Minneapolis Local Bus Transit Network under Random Link FailureAlgorithm 2. Robustness of the Minneapolis Local Bus Transit Network under Random Link FailureInput: initial random x(t)for all υi, ξ1, ξ2, wij, aij, k(i), s(i), f(xi(t))Output:RN1: Procedure Random link failure in Minneapolis local bus transit network2:  for all υi∈Vido3:     ifxi>1 or k(i)=04:      xi=0, ki=05:     else6:      xi=xi and k(i)=k(i)7:     end if8:   for all lki∈Li9:     all υi∈Vi10:   lki← Randomly select Link11:    do Eq. (6)12:   end for13:  ifxi(t+1)<114:   ki=kiout+kiin15:   do Eq. (5)16:   do Eq. (7)17:   f(xi(t+1))=4xi(1−xi)18:  else19:   xi(t+1)=xi(t+1)+R20:    for all υi∈Viwhereail=121:     do Eq. (6)22:    end for23:   end if24:     V←RN=N−Nf/N;25: end for26: end procedureData Availability StatementAll data, models, and code generated or used during the study appear in the published article.AcknowledgmentsThis work was supported by the Natural Science and Engineering Research Council of Canada (NSERC) through the Canadian Energy Infrastructure Resilience under Systemic Risk (CaNRisk)—Collaborative Research and Training Experience (CREATE) under Grant No. CREATE/482707-2016. The additional support through the INTERFACE Institute and the INViSiONLab is also acknowledged.NotationThe following symbols are used in this paper:Topological CharacteristicsAR(v)avg=average authority score;AR(v)max=maximum authority score;BCRavg=average betweenness centrality;BCRmax=maximum betweenness centrality;CRavg=average clustering;CRmax=maximum clustering;CCRavg=average closeness centrality;CCRmax=maximum closeness centrality;CCHRavg=average harmonic closeness centrality;CCHRmax=maximum harmonic closeness centrality;DRavg=average degree;DRmax=maximum degree;DRin-avg=average indegree centrality;DRin-max=maximum indegree;DRout-avg=average outdegree centrality;DRout-max=maximum outdegree;ECRavg=average Eigen centrality;ECRmax=maximum Eigen centrality;ER(v)avg=average eccentricity;ER(v)max=maximum eccentricity;HR(v)avg=average hub;HR(v)max=maximum hub;N(Din=1)%=percentage of nodes with in-degree = 1;N(Dout=1)%=percentage of nodes with out-degree =1;PRRavg=average page ranks;PRmax=maximum page ranks;QRavg=average modularity class;QRmax=maximum modularity class;SCRavg=average strong component number; andSCRmax=maximum strong component number.Operational CharacteristicsNL=number of links along the route;NN=number of stops along the route;(Pf/RC)Ravg=average Pf/RC;(Pf/RC)Rmax=maximum Pf/RC;SRavg=average weighted degree;SRmax=maximum weighted degree;SRout-avg=average weighted out-degree;SRout-max=maximum weighted out-degree;SRin-avg=average weighted in-degree; andSRin-max=maximum weighted in-degree.Robustness IndicatorsNf-L%=percentage of failed nodes based on link failure;Nf-N%=percentage of failed nodes based on node failure; andRFIR=route failure impact ratio.References Abdelaty, H., M. Mohamed, M. Ezzeldin, and W. El-Dakhakhni. 2020. “Quantifying and classifying the robustness of bus transit networks.” Transportmetrica A: Transp. Sci. 16 (3): 1176–1216. APTA (American Public Transportation Association). 2019. Public transportation ridership report: Fourth quarter 2019. Washington, DC: APTA. Barabási, A. L. 2016. Network science. Cambridge, UK: Cambridge University Press. Boyle, D. K. 2009. Passenger counting systems. Washington, DC: Transportation Research Board. Bruneau, M., S. E. Chang, R. T. Eguchi, G. C. Lee, T. D. O’Rourke, A. M. Reinhorn, M. Shinozuka, K. Tierney, W. A. Wallace, and D. Von Winterfeldt. 2003. “A framework to quantitatively assess and enhance the seismic resilience of communities.” Earthquake Spectra 19 (4): 733–752. Bureau of Transportation Statistics. 2018. Transportation statistics annual. Washington, DC: Bureau of Transportation Statistics. Cats, O., and E. Jenelius. 2015. “Planning for the unexpected: The value of reserve capacity for public transport network robustness.” Transp. Res. Part A: Policy Pract. 81 (Feb): 47–61. Chazottes, J. R., and B. Fernandez. 2005. Dynamics of coupled map lattices and of related spatially extended systems. Berlin: Springer. De Los Santos, A., G. Laporte, J. A. Mesa, and F. Perea. 2012. “Evaluating passenger robustness in a rail transit network.” Transp. Res. Part C: Emerging Technol. 20 (1): 34–46. Furth, P. G., B. Hemily, T. H. Muller, and J. G. Strathman. 2006. Using archived AVL-APC data to improve transit performance and management. Washington, DC: Transportation Research Board. Haggag, M., M. Ezzeldin, W. El-Dakhakhni, and E. Hassini. 2020. “Resilient cities critical infrastructure interdependence: A meta-research.” Sustainable Resilient Infrastruct. 12 (Nov): 1–22. Hosseini, S., and K. Barker. 2016. “Modeling infrastructure resilience using Bayesian networks: A case study of inland waterway ports.” Comput. Indus. Eng. 93 (Jan): 252–266. Hosseini, S., K. Barker, and J. E. Ramirez-Marquez. 2016. “A review of definitions and measures of system resilience.” Reliab. Eng. Syst. Saf. 145 (Jan): 47–61. Huang, A., H. M. Zhang, W. Guan, Y. Yang, and G. Zong. 2015. “Cascading failures in weighted complex networks of transit systems based on coupled map lattices.” Math. Problems Eng. 2015 (Jul): 1–15. Inanloo, B., B. Tansel, K. Shams, X. Jin, and A. Gan. 2016. “A decision aid GIS-based risk assessment and vulnerability analysis approach for transportation and pipeline networks.” Saf. Sci. 84 (Apr): 57–66. Kurant, M., and P. Thiran. 2006. “Extraction and analysis of traffic and topologies of transportation networks.” Phys. Rev. E: Stat. Nonlinear Soft Matter Phys. 74 (3): 036114. Litman, T. 2014. Evaluating public transit benefits and costs—Best practices guidebook. Victoria, BC: Victoria Transport Policy Institute. Lordan, O., J. M. Sallan, P. Simo, and D. Gonzalez-Prieto. 2014. “Robustness of the air transport network.” Transp. Res. Part E: Logist Transp. Rev. 68 (5): 155–163. Lu, D. 2008. Route level bus transit passenger origin-destination flow estimation using APC data: Numerical and empirical investigations. Columbus, OH: Ohio State Univ. Ma, F., W. Shi, K. F. Yuen, Q. Sun, X. Xu, Y. Wang, and Z. Wang. 2020. “Exploring the robustness of public transportation for sustainable cities: A double-layered network perspective.” J. Cleaner Prod. 265 (20): 121747. Martin, W. A., and N. A. McGuckin. 1998. Travel estimation techniques for urban planning. Washington, DC: Transportation Research Board. Metropolitan Council. 2009. Twin cities transit system performance evaluation. New York: Metropolitan Council. Metropolitan Council. 2019. Transit system performance evaluation-twin cities metropolitan region. New York: Metropolitan Council. Minneapolis City Council. 2019. Minneapolis 2040—The city’s comprehensive plan. Minneapolis: Minneapolis City Council. Minnesota Geospatial Information Office. 2019. “Minnesota Geospatial commons website.” Accessed October 1, 2021. Nian, G., F. Chen, Z. Li, Y. Zhu, and D. Sun. 2019. “Evaluating the alignment of new metro line considering network vulnerability with passenger ridership.” Transp. A: Transp. Sci. 15 (2): 1402–1418. Norris, F. H., S. P. Stevens, B. Pfefferbaum, K. F. Wyche, and R. L. Pfefferbaum. 2008. “Community resilience as a metaphor, theory, set of capacities, and strategy for disaster readiness.” Am. J. Commun. Psychol. 41 (1–2): 127–150. Salama, M., M. Ezzeldin, W. El-Dakhakhni, and M. Tait. 2020. “Temporal networks: A review and opportunities for infrastructure simulation.” Sustainable Resilient Infrastruct. 10 (Sep): 1–16. Shen, Y., G. Ren, and B. Ran. 2019. “Cascading failure analysis and robustness optimization of metro networks based on coupled map lattices: A case study of Nanjing, China.” Transportation 48 (2): 537–553. Sullivan, J. L., D. C. Novak, L. Aultman-Hall, and D. M. Scott. 2010. “Identifying critical road segments and measuring system-wide robustness in transportation networks with isolating links: A link-based capacity-reduction approach.” Transp. Res. Part A: Policy Pract. 44 (5): 323–336. Sun, D. J., and S. Guan. 2016. “Measuring vulnerability of urban metro network from line operation perspective.” Transp. Res. Part A: Policy Pract. 94 (3): 348–359. Sun, L., Y. Huang, Y. Chen, and L. Yao. 2018. “Vulnerability assessment of urban rail transit based on multi-static weighted method in Beijing, China.” Transp. Res. Part A: Policy Pract. 108 (5): 12–24. Wang, J., H. Mo, F. Wang, and F. Jin. 2011. “Exploring the network structure and nodal centrality of China’s air transport network: A complex network approach.” J. Transp. Geogr. 19 (4): 712–721. Wang, X., Y. Koç, S. Derrible, S. N. Ahmad, W. J. A. Pino, and R. E. Kooij. 2017. “Multi-criteria robustness analysis of metro networks.” Physica A 474 (Jan): 19–31. Weilant, S., A. Strong, and B. Miller. 2019. Incorporating resilience into transportation planning and assessment. Santa Monica, CA: RAND Corporation. Yang, Y., Y. Liu, M. Zhou, F. Li, and C. Sun. 2015. “Robustness assessment of urban rail transit based on complex network theory: A case study of the Beijing Subway.” Saf. Sci. 79 (Nov): 149–162. Yang, Y., S. T. Ng, F. J. Xu, and M. Skitmore. 2018. “Towards sustainable and resilient high density cities through better integration of infrastructure networks.” Sustainable Cities Soc. 42 (Jan): 407–422. Yassien, Y., M. Ezzeldin, M. Mohamed, and W. El-Dakhakhni. 2020. “Air transportation infrastructure robustness assessment for proactive systemic risk management.” J. Manage. Eng. 36 (4): 04020037. Zhang, J., F. Chen, and Q. Shen. 2019a. “Cluster-based LSTM network for short-term passenger flow forecasting in urban rail transit.” IEEE Access 7 (Dec): 147653–147671. Zhang, L., J. Lu, B. Fu, and S. Li. 2019b. “A cascading failures model of weighted bus transit route network under route failure perspective considering link prediction effect.” Physica A 523 (4): 1315–1330. Zhang, L., J. Lu, B. B. Fu, and S. Li. 2018. “A review and prospect for the complexity and resilience of urban public transit network based on complex network theory.” Complexity. Żochowska, R., and P. Soczówka. 2018. “Analysis of selected transportation network structures based on graph measures.” Zeszyty Naukowe. Gliwice, Poland: Politechnika Śląska.

Source link

Leave a Reply

Your email address will not be published.