Impact Estimation of Unplanned Urban Rail Disruptions on Public Transport Passengers: A Multi-Agent Based Simulation Approach

Once unplanned urban rail disruptions occur, it is essential to evaluate the impacts on public transport passengers since impact estimation results enable transit agencies to verify whether alternative transit services have adequate capacity to evacuate the affected rail passengers and to adopt effective emergency measures in response to the disruptions. This paper focuses on estimating the impacts of unplanned rail line segment disruptions on rail passengers as well as original bus passengers, as the latter are overlooked in existing studies. A method of identifying affected rail passengers based on passenger tap-in time is proposed, which is helpful for evaluating the scale and origin-destination distribution of the affected passengers. Passengers’ response behaviors are analyzed and modeled in a multi-agent simulation system. The system realizes the simulation of the multimodal evacuation process, in which a rule-based logit model is employed to describe passengers’ travel selection behavior and the Monte Carlo method is utilized to address the issue of uncertainty in passengers’ travel selection. In particular, the original bus passengers are integrated into the simulation and interact with rail passengers. Finally, some indicators assessing the impacts on rail passengers and bus passengers are presented, and a case study based on the Ningbo urban rail transit network is conducted.


Introduction
Urban rail transit (URT) shows a trend of large-scale and rapid development for cities worldwide. Due to continuous expansion, high-intensity operation and aging equipment, unplanned service disruptions occur frequently, causing negative impacts on a large number of passengers as well as the transportation system, not limited to the URT system. What exactly is the influence of rail disruption on passengers and the transportation system? It is essential for transit agencies to identify the impacts and adopt appropriate measures to minimize the negative effects of rail disruptions.
The robustness (resilience, redundancy, vulnerability, etc.) of the URT service partly depends on the topological structure of the URT network, so some studies focused on the effects on the URT system under random disruptions based on graph theory or complex network theory [1][2][3][4][5][6]. The identical nodes or links were identified by analyzing the topological characteristics, such as node degree, betweenness, centrality, etc. For example, Yang et al. [1] assessed the robustness of the Beijing URT network in face of unplanned disruptions. A new node importance evaluation index based on the weighted sum of betweenness and node degree was presented to measure node significance. Jian and Guan [3] 2 of 25 proposed a methodology to identify the most influential lines and measure the URT network vulnerability from a line operation perspective with consideration of the disruption probability. Chen et al. [6] developed a demand-impedance indicator and an effective path betweenness index to evaluate the importance of rail stations.
The previous studies belong to the static impact analysis of the URT service disruptions. The dynamic passenger flow was not taken into consideration. In practice, once a disruption occurs, affected passengers would make travel adjustments, including mode shift, resetting original and terminal points, and travel route selection. The travel adjustment choice of the affected passengers at the micro level is finally reflected in the redistribution of passenger flow on the transportation network at the macro level. Thus, some studies were dedicated to understanding the mechanism of passengers' behaviors and constructing passenger reassignment models under unplanned rail disruptions [7][8][9][10][11][12][13][14][15]. Generally, these works were conducted by means of a passenger flow survey. For instance, Teng and Liu [7] made a passenger behavior survey and a stated preference survey in Shanghai. They found that passengers' primary choice in response to the interruption is to make a detour in the URT network and that the temporary shuttle bus routes and existing bus lines are also welcomed by passengers as alternative modes. Then a multinomial logit model was built to assign the passenger flow under section interruption in the URT system. An online survey of rail users was conducted by Currie and Muir [8] to understand passenger behaviors, perceptions and priorities during unplanned rail disruptions. An interesting observation was obtained that over two-thirds of passengers choose the replacement buses as their substitute transit mode even though it takes a long time to wait for the bus service to operate. According to Zhu et al. [9], income played a vital role in determining travel reactions to disruptions. Nguyen-Phuoc et al. [10] conducted semi-structured interviews with 30 commuters from Melbourne, Australia, and found the most important factors influencing the mode shift when public transport ceases in the short term, e.g., car access, travel time, travel cost, trip importance, which were also observed by Adelé et al. [11] and Li et al. [12].
On the basis of survey results, logit models were commonly employed to make travel choice behavior analyses. Teng and Liu [7] built a multinomial logit model to assign the passenger flow under section interruption in the URT system. Considering uncertain disruption duration, Li et al. [14] developed a nested logit model to explore metro passengers' travel plan choice behavior under unplanned service disruptions. Wang et al. [15] established a nested logit model following random regret minimization principles to estimate passenger travel choice behaviors under metro emergency context. Dai et al. [16] also proposed a nested logit model for the decision-making of affected passengers in a metro emergency evacuation.
Apart from survey-based approaches, some scholars used data-driven methods to explore the influence of unplanned rail disruptions on the URT passenger flows. The historical passenger flow data collected from the automated fare collection (AFC) system under normal states and disrupted states were compared and analyzed, hence capturing the passenger volume changes and estimating the impacts of rail disruptions [17][18][19][20][21][22]. For instance, Silva et al. [17] used network-wide data obtained from smart cards in the London transport system to predict future traffic volumes and to estimate the effects of disruptions due to unplanned closures of stations or lines. Sun et al. [18] introduced a Bayesian model to identify disruptions based on automated fare collection (AFC) data, and the effects of disruptions were modeled concerning the affected passengers and the length of delay in the URT network. Through analyzing mobile phone records before and after a collision accident of a metro line in Shanghai, Duan et al. [19] studied the passenger evacuation process and the impacts of the accident on commuters. Based on smart card data, Eltved et al. [20] proposed a method to analyze the effects of long-term disruptions on passenger travel behavior by distinguishing travel behavior changes prior to the disruption and after the disruption.
Deterministic methods or aggregated methods cannot reflect the heterogeneity of passengers' behaviors under unplanned rail disruptions, while simulation is a better approach to model the heterogeneous passengers. On one hand, agent-based simulation was widely used [23][24][25][26][27][28][29]. Leng and Corman [24] applied agent-based micro-simulation (MATSim) to imitate large-scale passengers' behaviors during public transport disruptions. A within-day replanning method was integrated into MATSim to study passengers' route choices. The research revealed the importance of the issue time of information, which greatly impacts passenger satisfaction with public transport disruptions. A similar study was conducted by Muller et al. [26], where MATSim was used to simulate agents' choice to find an alternative route after a shutdown of an underground line. Based on simulation results with different policies, the effect of the disruption was analyzed, including travel time, number of line changes, and line usage. Li et al. [27] developed a discrete-event simulation model with parallel computing to estimate the effects of metro emergencies. Zargayounaa et al. [28,29] presented a multi agent simulation model to evaluate the impact of real-time information provision on passenger travel times. Results showed that the impact was positive until a threshold was reached.
On the other hand, Monte Carlo (MC) is another simulation method to analyze the uncertainty of passengers' behaviors. Seger and Kisgyörgy [30] utilized MC simulation to quantify uncertainty in traffic assignment on a transport network. Similarly, the MC technique was applied by Bhat [31] to evaluate the variations in travel mode choice in a multinomial logit model. To examine the day-to-day evolution of network traffic flow, Wang and Li [32] constructed a mixed logit model and designed an MC simulation algorithm to solve the model. Huang et al. [33] proposed an MC-based method to evaluate travel time reliability under the impact of traveler information.
However, there are some limitations in the previous studies. First, the existing works failed to provide a method to accurately identify and estimate the affected passenger flow. Second, concerning passengers' modal choice behavior, the existing studies only considered passengers' individual choices based on the utility of each alternative mode, such as cost, travel time, and transfer times. Actually, due to the limited capacity of transit modes, the choices of passengers who access the URT network in advance would affect those of passengers who enter the network later. However, the complex interactions among different passengers and the transit system were neglected in previous research. Third, most previous studies, model-based, or simulation-based, only performed the impact analysis of unplanned rail disruptions on affected rail passengers and rail networks, e.g., passenger volume, delay, travel time, and rail usage. As a matter of fact, the impacts are not limited to the interior of the rail system but extend to the whole transportation system. For instance, buses are a relatively popular choice for affected rail passengers as an alternative transportation mode to complete their travel [8,34]. When a large number of rail passengers shift to the bus system, the original bus passengers are certain to be affected. To the best of our knowledge, no other study focused on the impacts of unplanned rail disruptions on bus passengers.
To fill the research gap, this paper proposes a passenger behavior simulation framework under urban rail disruptions based on the Monte Carlo method and logit model. Affected rail passengers along with conventional bus passengers are modeled as agents that interact with each other and with the transport system. The heterogeneity of passengers' travel choices can be simulated in a multi-modal public transit network, including the urban rail, bus, taxi, shared bikes and walking. On the basis of simulation, the impacts of the disruptions on rail passengers as well as bus passengers can be analyzed. The results can enable transit agencies to evaluate the passenger evacuation capacity of the present public transit system under unplanned rail disruptions. They can further assist transit management departments in taking appropriate measures to minimize the negative effects on commuters.
The rest of the paper is organized as follows. Section 2 proposes a method to identify the affected passengers under unplanned rail disruptions. Section 3 builds a multi-agent simulation framework to model the travel adjustment behavior of the affected passengers based on the logit model and Monte Carlo method. Then the impact evaluation indices of disruptions are introduced in Section 4. Section 5 discusses a case study. Finally, the conclusions and future research are summarized in Section 6.

Identification of Affected Rail Passengers
This study focuses on unplanned URT disruptions that result in the closure of rail links between two adjacent turnover stations in both directions [35]. Note that the interrupted section may contain transfer stations or not. The former case, i.e., containing transfer stations, is obviously a more complicated situation, which is studied in-depth in this paper. As shown in Figure 1, Line 1 and Line 2 are both bidirectional rail lines. An unplanned disruption causes service breakdown on a segment of Line 1 including the transfer station A, while Line 2 still functions normally.
The rest of the paper is organized as follows. Section 2 proposes a method to identify the affected passengers under unplanned rail disruptions. Section 3 builds a multi-agent simulation framework to model the travel adjustment behavior of the affected passengers based on the logit model and Monte Carlo method. Then the impact evaluation indices of disruptions are introduced in Section 4. Section 5 discusses a case study. Finally, the conclusions and future research are summarized in Section 6.

Identification of Affected Rail Passengers
This study focuses on unplanned URT disruptions that result in the closure of rail links between two adjacent turnover stations in both directions [35]. Note that the interrupted section may contain transfer stations or not. The former case, i.e., containing transfer stations, is obviously a more complicated situation, which is studied in-depth in this paper. As shown in Figure 1, Line 1 and Line 2 are both bidirectional rail lines. An unplanned disruption causes service breakdown on a segment of Line 1 including the transfer station A, while Line 2 still functions normally. The passengers affected by the unplanned disruption refer to those whose space-time travel path in the URT network coincides with the occurrence place and duration of the disruption simultaneously. Under the scenario mentioned before, the affected passengers can be divided into two types in terms of whether they need to transfer in the interrupted line section. In this section, we analyze the time domain of the two types of affected passengers entering the URT system, respectively. On this basis, the affected passenger volume can be estimated using historical automatic fare collection (AFC) data. To identify the affected passengers by their travel route, it is assumed that passengers follow the shortest path to complete their trip.    The passengers affected by the unplanned disruption refer to those whose space-time travel path in the URT network coincides with the occurrence place and duration of the disruption simultaneously. Under the scenario mentioned before, the affected passengers can be divided into two types in terms of whether they need to transfer in the interrupted line section. In this section, we analyze the time domain of the two types of affected passengers entering the URT system, respectively. On this basis, the affected passenger volume can be estimated using historical automatic fare collection (AFC) data. To identify the affected passengers by their travel route, it is assumed that passengers follow the shortest path to complete their trip.

Notations
S: the set of rail stations, S = {S 1 , S 2 , . . . , S n }. S o : the origin station of a passenger. S d : the destination station of a passenger. S h : the transfer station in the interrupted section. S w : the two turnover stations adjacent to the interrupted section. S z : the first interrupted station in the forward direction of a passenger's path. S l : the last interrupted station in the forward direction of a passenger's path. S k : the transfer station outside the interrupted section and nearest to the interrupted section in a passenger's detour paths.
T interrupt : the time of the occurrence of the rail disruption. T: the passenger tap-in time. t interrupt : the time duration of the rail disruption. T recover : the time of the recovery of the rail disruption, i.e., T recover = T interrupt + t interrupt . Affected passengers who need not transfer in the interrupted section can be classified into four categories in terms of the spatial relationship between passengers' origins and destinations and the interrupted section (see Table 1). Based on the passenger tap-in time, we can identify whether a certain passenger is influenced by the disruption or not. With this regard, the affected tap-in time domain is determined for each category as follows.

Passenger category
The first category : t up = 0 and t down = 0 The sec ond category : t up = 0 and t down > 0   Affected passengers who need not transfer in the interrupted section can be classified into four categories in terms of the spatial relationship between passengers' origins and destinations and the interrupted section (see Table 1). Based on the passenger tap-in time, we can identify whether a certain passenger is influenced by the disruption or not. With this regard, the affected tap-in time domain is determined for each category as follows.  T : the passenger tap-in time. Affected passengers who need not transfer in the interrupted section can be classified into four categories in terms of the spatial relationship between passengers' origins and destinations and the interrupted section (see Table 1). Based on the passenger tap-in time, we can identify whether a certain passenger is influenced by the disruption or not. With this regard, the affected tap-in time domain is determined for each category as follows.

Affected Non-Transfer Passengers
Affected passengers who need not transfer in the interrupted section can be classified into four categories in terms of the spatial relationship between passengers' origins and destinations and the interrupted section (see Table 1). Based on the passenger tap-in time, we can identify whether a certain passenger is influenced by the disruption or not. With this regard, the affected tap-in time domain is determined for each category as follows.

Affected Non-Transfer Passengers
Affected passengers who need not transfer in the interrupted section can be classified into four categories in terms of the spatial relationship between passengers' origins and destinations and the interrupted section (see Table 1). Based on the passenger tap-in time, we can identify whether a certain passenger is influenced by the disruption or not. With this regard, the affected tap-in time domain is determined for each category as follows.   Affected passengers of the first category are those whose origins and destinations are both in the interrupted section (i.e., t up = 0 and t down = 0). For passengers who enter the URT network before T interrupt , the affected time domain is T interrupt − t metro o,d < T < T interrupt ; for those who enter the network later than T interrupt and depart from the network earlier than T recover , the affected time domain is T interrupt ≤ T ≤ T recover − t metro o,d ; and for those who enter the network after T interrupt and leave the network after T recover , the affected time Similarly, the affected time domain of the other three categories can be analyzed. The results are shown in Table 1.

Affected Transfer Passengers
As analyzed to non-transfer passengers, affected passengers who need to transfer in the interrupted section can also be segmented into four categories. Note that the first category is inexistent for affected transfer passengers. The other three categories are explained in Table 2.

Passenger category
The second category : t up = 0 and t down > 0 The third category : t up > 0 and t down = 0

Affected Transfer Passengers
As analyzed to non-transfer passengers, affected passengers who need to transfer in the interrupted section can also be segmented into four categories. Note that the first category is inexistent for affected transfer passengers. The other three categories are explained in Table 2.

Data Preparation
(1) The URT network diagram We define the URT network as a directed graph is the rail node set, corresponding to rail stations, and  Table 1.

Affected Transfer Passengers
As analyzed to non-transfer passengers, affected passengers who need to transfer in the interrupted section can also be segmented into four categories. Note that the first category is inexistent for affected transfer passengers. The other three categories are explained in Table 2.

Data Preparation
(1) The URT network diagram We define the URT network as a directed graph is the rail node set, corresponding to rail stations, and (2) The URT disruption data The data on the rail disruption including the interrupted section and some time-related parameters, i.e., int errupt T , int errupt t , cov re er T , are assumed to be known.

Passenger category
The fourth category : t up > 0 and t down = 0 (transfer from Line 1 to Line 2) The fourth category : t up > 0 and t down > 0(transfer from Line 2 to Line 1) As analyzed to non-transfer passengers, affected passengers who need to transfer in the interrupted section can also be segmented into four categories. Note that the first category is inexistent for affected transfer passengers. The other three categories are explained in Table 2. We define the URT network as a directed graph is the rail node set, corresponding to rail stations, and categories can be analyzed. The results are shown in Table 1.

Affected Transfer Passengers
As analyzed to non-transfer passengers, affected passengers who need to transfer in the interrupted section can also be segmented into four categories. Note that the first category is inexistent for affected transfer passengers. The other three categories are explained in Table 2. We define the URT network as a directed graph is the rail node set, corresponding to rail stations, and is the rail node set, corresponding to rail stations, and E = {e 1 , e 2 , . . . e n } is the rail arc set, representing links between adjacent stations.
(2) The URT disruption data The data on the rail disruption including the interrupted section and some time-related parameters, i.e., T interrupt , t interrupt , T recover , are assumed to be known.

(3) The shortest paths of all the OD in the URT network
The Dijkstra algorithm is applied to acquire the shortest paths of all the OD in the URT network, which constitute the database of PathData. Correspondingly, the travel time of each shortest path is stored into the database. Based on the shortest path of each OD, we can determine its affected category and time domain that are also laid in PathData.
(4) The AFC data After data cleaning, the AFC data of a certain day are stored in the database of OriginalURTData. Each piece of data represents a passenger trip, including the passenger ID, the tap-in station ID, the tap-in time, the tap-out station ID, and the tap-out time.

Identification Process
Based on the AFC data and the affected tap-in time domains, affected passengers can be identified. The basic idea is that for a piece of data, we clarify its affected category and compare its tap-in time with the corresponding affected time domain. However, it is tedious and time-consuming to make the comparison for every piece of AFC data, so we simplify the identification process in such a way that the data are sorted in the light of the tap-in time and then some data are excluded using the upper and lower limits of the time domains. The detailed process is listed as follows: Step 1: Compare the affected tap-in time domains of all OD pairs in PathData. Find out the smallest lower limit T in min and the largest upper limit T in max .
Step 2: Sort the data in OriginalURTData according to the tap-in time. Store the data with the tap-in time within T in min and T in max into URTData.
Step 3: Divide the data in URTData into several groups based on the OD. As a result, data with the same OD are included in the same group. Then sort the data in each group according to the tap-in time.
Step 4: For data in each group, find the affected time domain in PathData. Then identify the data within the domain and store them into AffectedURTData.

Trip Plan Selection
Under unplanned rail service disruptions, affected passengers are forced to make adjustments to their trip plans. Figure 2 summarizes the possible passenger behaviors in response to the disruptions. There are five options: (1) changing to an alternative route in the remaining URT network, (2) changing to other modes, (3) using rail and other modes with the latter functions as a "bridge" to restore connectivity between disrupted stations, (4) waiting until the URT system is restored, and (5) delaying or canceling the trip. Other alternative modes include walking, shared bike, bus, and taxi. Note that the response of altering the origin or destination is beyond our consideration. Here Option (4) is valid for passengers who have already entered the URT system before T interrupt , while Option (5) is valid to passengers whose T is later than T interrupt .
(4) The AFC data After data cleaning, the AFC data of a certain day are stored in the database of Orig-inalURTData. Each piece of data represents a passenger trip, including the passenger ID, the tap-in station ID, the tap-in time, the tap-out station ID, and the tap-out time.

Identification Process
Based on the AFC data and the affected tap-in time domains, affected passengers can be identified. The basic idea is that for a piece of data, we clarify its affected category and compare its tap-in time with the corresponding affected time domain. However, it is tedious and time-consuming to make the comparison for every piece of AFC data, so we simplify the identification process in such a way that the data are sorted in the light of the tap-in time and then some data are excluded using the upper and lower limits of the time domains. The detailed process is listed as follows: Step 1: Compare the affected tap-in time domains of all OD pairs in PathData. Find out the smallest lower limit min in T and the largest upper limit max in T .
Step 2: Sort the data in OriginalURTData according to the tap-in time. Store the data with the tap-in time within min in T and max in T into URTData.
Step 3: Divide the data in URTData into several groups based on the OD. As a result, data with the same OD are included in the same group. Then sort the data in each group according to the tap-in time.
Step 4: For data in each group, find the affected time domain in PathData. Then identify the data within the domain and store them into AffectedURTData.

Trip Plan Selection
Under unplanned rail service disruptions, affected passengers are forced to make adjustments to their trip plans. Figure 2 summarizes the possible passenger behaviors in response to the disruptions. There are five options: (1) changing to an alternative route in the remaining URT network, (2) changing to other modes, (3) using rail and other modes with the latter functions as a "bridge" to restore connectivity between disrupted stations, (4) waiting until the URT system is restored, and (5) delaying or canceling the trip. Other alternative modes include walking, shared bike, bus, and taxi. Note that the response of altering the origin or destination is beyond our consideration. Here Option (4) is valid for passengers who have already entered the URT system before int errupt T , while Option (5) is valid to passengers whose T is later than int errupt T .  The options for affected passengers of different categories are somewhat different from each other. As shown in Table 3, Options (2), (4) and (5) are available to every passenger category; Options (3) is not optional to only the first category, while Option (1) is available to only the fourth category.  (1) is no longer available to them. Considering that it takes a while for the disruption information dissemination if passengers obtain the information before they reach S k , they can resume their travel by detouring in the degraded URT network; otherwise, they miss the detour opportunities due to lack of information. The tap-in time domains in which Option (1) is available to passengers of the fourth category are displayed in Table 4.  (1) is available to passengers of the fourth category.

Passengers of the Fourth Category The Tap-In Time Domain
Non-transfer passengers Transfer passengers From Line 1 to Line 2 In order to determine passengers' decisions on trip plan selection, the following rules are established: (1) Affected passengers would select Option (1) if feasible detour routes exist in the degraded URT network. This assumption has been proven to be true in reality by Duan et al. [19]. Feasible detour paths shall meet the following two conditions: a. the travel time of the path should not exceed a predetermined limit t limit detour .
b. the total transfer times of the path should be no more than a predetermined limit F trans f er . f trans f er ≤ F trans f er .
(2) Transit agencies can take passenger flow guidance measures to direct passengers with short travel distances to choose slow traffic. With this in mind, we assume that if there is just one station between a passenger's current located station and his destination station, he will use a shared bike when there are available shared bikes. Additionally, if a passenger's current located station and his destination station are adjacent stations, passengers will choose walking. (3) The proportion of travelers canceling or delaying their trips is supposed as a given constant. (4) The duration of service disruption is predictable, and it is announced to passengers.
Passengers prefer to wait until the system is restored, i.e., select Option (4), if the announced recovery time is less than a predetermined value (denoted as t wait−recover .)

Candidate Trip Chains
Based on the above rules, affected passengers choosing Options (1), (4) and (5) are definite on some specific conditions. However, Options (2) and (3) are uncertain options that heterogeneous passengers may have distinct choices. Comparatively, Option (3) is more complicated in that for different categories of affected passengers, candidate trip chains of Option (3) exhibit differences. Herein, the trip chain means that travel comprises several trip legs, and the transit modes for adjacent legs are different. Especially, Option (2) can be regarded as a trip chain with only one leg.
We further define the starting and ending stations of the trip leg in Options (2) (2) and (3) can be derived (see Figure 3). Trip chain 3 belongs to Option (2), while the other five trip chains belong to Option (3). In Figure 3, "pre" and "after" represent the trip leg before and after disruption, respectively. Furthermore, Table 5

The Logit Model and Monte Carlo Method
Note that "other" in the trip chains means the transit modes other than rail include walking, shared bicycle, bus and taxi, so with distinct transit modes, a trip chain can evolve into various trip chains. Thus, passengers of each category in either Group (i) or (ii) have more than one choice. Here a classical logit model is constructed to describe the trip chain selection for affected passengers [36] (see Equation (3)).
Given an OD pair, r od U represents the utility of the rth trip chain. As seen in Equation (4), it consists of the generalized cost of the rth trip chain and a random item. The

The Logit Model and Monte Carlo Method
Note that "other" in the trip chains means the transit modes other than rail include walking, shared bicycle, bus and taxi, so with distinct transit modes, a trip chain can evolve into various trip chains. Thus, passengers of each category in either Group (i) or (ii) have more than one choice. Here a classical logit model is constructed to describe the trip chain selection for affected passengers [36] (see Equation (3)).
Given an OD pair, U r od represents the utility of the rth trip chain. As seen in Equation (4), it consists of the generalized cost of the rth trip chain and a random item. The generalized cost is composed of the money cost and the cost of travel time, which can be divided into in-vehicle time, transfer time and waiting time (see Equation (5)).
where α denotes the transfer penalty coefficient; β denotes the waiting penalty coefficient; and δ denotes the value of time.
Deterministic methods ignore the randomness in passengers' travel adjustment behavior, which may result in estimation bias [22]. In this study, the Monte Carlo method is utilized to address the issue of uncertainty in passengers' trip chain selection. Specifically, for each passenger who has to select a trip chain, we perform a sampling based on the logit model, i.e., generate a random number representing the index of a trip chain according to the proportions of trip chains [37][38][39]. In this study, the Monte Carlo method can be taken as a subroutine of the multi-agent simulation introduced in the next section.

Multi-Agent Simulation Model
Passengers' trip chain selection depends on the availability of the transport modes. Specifically, passengers can only choose transport modes with residual capacity. Note that the residual capacity at a certain time is influenced by passengers' selection behavior before that time. In other words, passengers' decision-making on trip chain selection interacts with each other. To better reflect this interaction, a multi-agent simulation model is developed to simulate the process of passengers' travel adjustment behavior in a multimodal transport system with dynamic changes.

Multi-Agent System
Considering the complexity and dynamics of the multi-modal transport system, it is preferable to split the system into several independent and interrelated parts [40]. Each part is regarded as a type of agent, including the passenger agent, vehicle agent and station agent. Information is able to be exchanged among these agents. The logical framework of the multi-agent simulation system is shown in Figure 4. The basic idea is that first, the database is constructed, then the simulation experiment is carried out, and finally, the impacts of an unplanned rail disruption are estimated and output.
The basic data, including the infrastructure data, passenger data, and vehicle data, are predetermined as input information for the simulation system. The infrastructure data mainly comprises information on the disrupted rail stations and related bus stations. The passenger data include attributes of affected rail passengers and original bus passengers, e.g., the origin station, the destination station and the tap-in time. The vehicle in the simulation system consists of trains, buses, taxis and shared bikes. The train and bus data mainly involve the capacity, the load factor, the remaining capacity and the timetable, while the taxi and shared bike data mainly contain the number of available vehicles, the turnover rate and the located station.
Considering the complexity and dynamics of the multi-modal transport system, it is preferable to split the system into several independent and interrelated parts [40]. Each part is regarded as a type of agent, including the passenger agent, vehicle agent and station agent. Information is able to be exchanged among these agents. The logical framework of the multi-agent simulation system is shown in Figure 4. The basic idea is that first, the database is constructed, then the simulation experiment is carried out, and finally, the impacts of an unplanned rail disruption are estimated and output. The basic data, including the infrastructure data, passenger data, and vehicle data, are predetermined as input information for the simulation system. The infrastructure data mainly comprises information on the disrupted rail stations and related bus stations. The passenger data include attributes of affected rail passengers and original bus passengers, e.g., the origin station, the destination station and the tap-in time. The vehicle in the simulation system consists of trains, buses, taxis and shared bikes. The train and bus data mainly involve the capacity, the load factor, the remaining capacity and the timetable, while the taxi and shared bike data mainly contain the number of available vehicles, the turnover rate and the located station.

Passenger Agent
As the fundamental simulation unit, the passenger agent connects and interacts with other agents in the multi-agent simulation model. Passenger agents are classified into two types in this paper: rail passenger agents and bus passenger agents, representing the affected rail passengers and the original bus passengers in the real world, respectively.
For the rail passenger agents, the simulation model mainly simulates the process of train matching, the process of trip chain selection, and the process of using alternative modes. The logical flow chart of the rail passenger agent simulation is displayed in Figure 5. As analyzed before, passengers selecting Options (1), (4) and (5) can be determined by rules, so they are first picked up from the affected passengers. In order to reduce the simulation complexity, this study mainly focuses on the simulation of passengers in the interrupted section. Thus, passengers not entering the interrupted section, e.g., those choosing trip chains 3 and 4 in Group (ii) are also excluded beforehand. Then other passengers are put into the main simulation program. Train matching determines S o−new and the time passengers arrive at S o−new based on the train timetable and passenger tap-in time. Trip chain selection is realized by the logit model and Monte Carlo method stated in Section 3.2.2. Additionally, the passenger agents choosing the bus as their substitute mode will be integrated with bus agents and proceed with the waiting, boarding and alighting process.
Bus passenger agents have the function of waiting for bus, boarding the bus and alighting the bus (see Figure 6). The bus passenger agent is matched with a bus agent according to the time passengers reach the station, the bus timetable, the bus residual capacity and its sequence in the waiting queue. The waiting time of each bus passenger agent can be determined and recorded for the subsequent impact analysis.

Vehicle Agent
As the carrier for the passenger agent, the vehicle agent is to realize the function of simulating the boarding and alighting process by interacting with the passenger agent and the infrastructure agent. Vehicle agents can be divided into train agents, bus agents, taxi agents and share bike agents. Trip chain selection is realized by the logit model and Monte Carlo method stated in Section 3.2.2. Additionally, the passenger agents choosing the bus as their substitute mode will be integrated with bus agents and proceed with the waiting, boarding and alighting process.

Int. J. Environ. Res. Public Health 2022, 19, x FOR PEER REVIEW
Bus passenger agents have the function of waiting for bus, boarding the alighting the bus (see Figure 6). The bus passenger agent is matched with a bu according to the time passengers reach the station, the bus timetable, the bus resi pacity and its sequence in the waiting queue. The waiting time of each bus pa agent can be determined and recorded for the subsequent impact analysis.

Vehicle Agent
As the carrier for the passenger agent, the vehicle agent is to realize the fun simulating the boarding and alighting process by interacting with the passeng and the infrastructure agent. Vehicle agents can be divided into train agents, bus The train agent mainly realizes the departure function and arrival function based on the train timetable. When a train arrives at a station, it will dwell for a while. Subsequently, rail passenger alighting and boarding events will be triggered. In addition, when the scheduled departure time is up, the departure function is called. As depicted in Figure 7, the bus agent also performs arrival and departure functions. It is worth noting that both bus passenger agents and rail passenger agents who choose the bus as the substitute mode to pass through the interrupted section are involved in the bus operation simulation. Unlike train agents and bus agents, one taxi agent or shared bike agent represents several taxis and bikes at a certain station in the real world, so both train agents and bus agents have the attribute of residual capacity indicating the number of available carriers.

Impact Estimation of Unplanned Rail Disruptions
In evaluating the impacts of unplanned rail disruptions, our main concern is the impacts on passengers, including rail passengers and non-rail passengers. Note that the existing parallel bus system is one of the first alternatives for affected rail passengers to resume their trips. The original bus travelers are certain to be influenced when a large number of rail passengers shift to the bus system. Thus, for the impacts on non-rail passengers, we focus on the delay time of bus travelers caused by rail disruptions. As indicated before, our simulation focuses on passengers in the interrupted section, so here only these passengers are taken into the following impact estimation.

Estimation of the Impacts on Rail Passengers
(1) Affected rail passenger flow As noted before, affected rail passengers refer to their travel paths intersecting with the rail disrupted section in space and the rail disruption occurrence period in time. The affected rail passenger flow, denoted as affacted metro N , is the total number of affected rail passengers. This index can be estimated using the approach proposed in Section 2.
(2) The stranded time Herein, the stranded time of the ith affected rail passenger, denoted as i stranded t , is defined as the time duration the passenger is detained outside a rail station until he is taken away by a certain transit mode. In particular, for passengers who choose walking or

Impact Estimation of Unplanned Rail Disruptions
In evaluating the impacts of unplanned rail disruptions, our main concern is the impacts on passengers, including rail passengers and non-rail passengers. Note that the existing parallel bus system is one of the first alternatives for affected rail passengers to resume their trips. The original bus travelers are certain to be influenced when a large number of rail passengers shift to the bus system. Thus, for the impacts on non-rail passengers, we focus on the delay time of bus travelers caused by rail disruptions. As indicated before, our simulation focuses on passengers in the interrupted section, so here only these passengers are taken into the following impact estimation.

Estimation of the Impacts on Rail Passengers
(1) Affected rail passenger flow As noted before, affected rail passengers refer to their travel paths intersecting with the rail disrupted section in space and the rail disruption occurrence period in time. The affected rail passenger flow, denoted as N a f f acted metro , is the total number of affected rail passengers. This index can be estimated using the approach proposed in Section 2.
(2) The stranded time Herein, the stranded time of the ith affected rail passenger, denoted as t i stranded , is defined as the time duration the passenger is detained outside a rail station until he is taken away by a certain transit mode. In particular, for passengers who choose walking or shared bike to resume their trips, their stranded time is set to be 0. The total number of affected rail passengers in the interrupted section is assumed to be N a f f acted−in metro . Consequently, the total stranded time can be formulated as follows: The average and maximum stranded time can also be formulated as Equations (7) and (8), respectively.
Compared to the delay time commonly defined as the increase in travel time under disruptions, the stranded time indicating the additional waiting time for alternative transit services is more concerned by affected rail passengers. On the other hand, transport authorities are also more interested in the stranded time which suggests the evacuating speed of affected passengers than the travel time increase [41].
(3) Stranded rail passenger flow Due to the limited capacity of available transport services, some affected rail passengers may be stranded outside a rail station until the disrupted rail service is recovered. Let N stranded−in metro be the number of these stranded passengers in the interrupted section and the corresponding proportion can be calculated as follows: The larger α stranded−in metro is, the greater the gap between the residual capacity of the available transit modes and the scale of affected rail passengers. It suggests that other substitute transport services are required to evacuate stranded passengers, such as the bus bridging service.  The total delay time can be estimated by the following equation:

Estimation of the Impacts on Original Bus Passengers
where N stop is the number of bus stops and V j is the number of bus passengers at the jth bus stop during the rail disruption. The average and maximum delay time can be obtained by Equations (12) and (13).
(2) Affected bus passenger flow Suppose bus headway to be h bus . If t i,j wait−bus − h bus > 0, it reveals that the ith bus passenger is affected by the rail disruption. On this basis, a binary variable η i,j is introduced to indicate whether the ith bus passenger at the jth bus stop in the interrupted section is affected (see Equation (14)).
Subsequently, the total number of bus passengers in the interrupted section who are affected by the disruption can be estimated by Equation (15).

N a f f acted−in bus
(

3) Stranded bus passenger flow
Since affected rail passengers take up some bus resources, several affected bus passengers would be unable to board buses and are stranded at bus stops until the disrupted rail service is recovered. Let ω i,j be a binary variable indicating whether the ith bus passenger at the jth bus stop in the interrupted section is finally stranded. If yes, ω i,j = 1; otherwise,

Case Study
In this section, part of the Ningbo URT network is employed to conduct a case study. As shown in Figure 8, the test network consists of 3 bi-directional rail lines and 75 stations, of which 3 are transfer stations.

Scenario Description
Two scenarios of unplanned rail disruption are designed here to make a comparison. As seen in Table 6, both of the interrupted sections in the two scenarios contain the transfer station S 10 . Compared with Scenario B, Scenario A is a comparatively smaller disruption with fewer interrupted stations. Additionally, the disruption in Scenario A happens in non-peak hours of a working day, while that in Scenario B arises in peak hours. The historical AFC data of a typical weekday are stored in the database of OriginalURTData and applied to make the impact estimation.

Case Study
In this section, part of the Ningbo URT network is employed to conduct a case study. As shown in Figure 8, the test network consists of 3 bi-directional rail lines and 75 stations, of which 3 are transfer stations.   As indicated before, this study focuses on the simulation of passengers in the interrupted section. Therefore, three bus lines parallel to the interrupted section are integrated into the simulation system, and the impact analysis on the original bus passengers also targets these three bus lines. Related data on these bus routes, including data on passengers as well as buses, are stored in the database of the simulation system. Referring to the existing research [8,22,42], some parameters in the simulation system are assigned as listed in Table 7.

Results Analysis
Since the simulation involves a stochastic trip chain selection process by the Monte Carlo method, the results of each simulation are distinct. We perform the simulation 50 times. It appears that there are small differences among those simulation results. For instance, regarding the average stranded time for rail passengers, the simulation results of Scenario B are between 14.97 and 15.53, with a standard deviation of 0.15. Thus, the average value of 50 simulation results is taken as the estimation result of each indicator. Some proportion-related data are given using a certain simulation result.

Impact Estimation Results for Rail Passengers
(1) Estimation of affected rail passenger flow Using the identification method described in Section 2.3, we obtain the estimated results of N a f f acted metro for scenarios A and B as 6959 and 22,734, respectively. Figure 9 depicts more specific information on the affected rail passenger flow composition, i.e., the number of passengers and OD pairs in the four categories. It can be concluded that passengers of the fourth category take the majority of the total affected rail passengers. Therefore, emergency buses should be dispatched not only to the interrupted rail stations but also to stations of other uninterrupted rail lines. This finding is in line with the observations of Duan et al. [17].  Figure 9 depicts more specific information on the affected rail passenger flow composition, i.e., the number of passengers and OD pairs in the four categories. It can be concluded that passengers of the fourth category take the majority of the total affected rail passengers. Therefore, emergency buses should be dispatched not only to the interrupted rail stations but also to stations of other uninterrupted rail lines. This finding is in line with the observations of Duan et al. [17]. The results of trip plan selection and trip chain selection for affected rail passengers are shown in Figure 10 and Table 8, respectively. The result of Scenario A is illustrated using the inner circle in Figure 10, while that of Scenario B is illustrated using the outer circle. It can be seen that the outcomes of the two scenarios have small differences. Options (1) and (3) take the majority in both scenarios, indicating that most affected rail passengers would continue to use the degraded URT network to resume their travel [17]. The results of trip plan selection and trip chain selection for affected rail passengers are shown in Figure 10 and Table 8, respectively. The result of Scenario A is illustrated using the inner circle in Figure 10, while that of Scenario B is illustrated using the outer circle. It can be seen that the outcomes of the two scenarios have small differences. Options (1) and (3) take the majority in both scenarios, indicating that most affected rail passengers would continue to use the degraded URT network to resume their travel [17].  To further analyze passengers in the interrupted section, the outbound passengers at each interrupted station needing to be evacuated by a certain transportation mode other than rail are described in Figure 11. As can be seen, most affected rail passengers depart from the URT system at the transfer station in the interrupted section as well as the turnover stations. This suggests that emergency buses should be dispatched to these stations to evacuate passengers and avoid secondary accidents caused by passenger gathering. Furthermore, Figure 11 shows the modal split results of the "other" mode in the trip chains. Additionally, a special type of passengers who are stranded until the disrupted rail service is recovered is also added into Figure 11. Based on the preset rules listed in Section 3.1, there are many passengers walking to their destination stations. It relieves the pressure of passenger evacuation on other transit modes, e.g., buses, taxis. This confirms that passenger flow guidance is extremely important in response to unplanned disruptions.   Trip chain 1  3  18  13  43  257  861  299  1028  Trip chain 2  38  103  1498  4025  Trip chain 3  26  93  137  564  65  291  82  394  Trip chain 4  370  1167  319  1368  Trip chain 5  161  784  263  874  Trip chain 6  698  2175 To further analyze passengers in the interrupted section, the outbound passengers at each interrupted station needing to be evacuated by a certain transportation mode other than rail are described in Figure 11. As can be seen, most affected rail passengers depart from the URT system at the transfer station in the interrupted section as well as the turnover stations. This suggests that emergency buses should be dispatched to these stations to evacuate passengers and avoid secondary accidents caused by passenger gathering. Furthermore, Figure 11 shows the modal split results of the "other" mode in the trip chains. Additionally, a special type of passengers who are stranded until the disrupted rail service is recovered is also added into Figure 11. Based on the preset rules listed in Section 3.1, there are many passengers walking to their destination stations. It relieves the pressure of passenger evacuation on other transit modes, e.g., buses, taxis. This confirms that passenger flow guidance is extremely important in response to unplanned disruptions. (2) Estimation of the stranded time and stranded rail passenger flow As illustrated in Table 9, in Scenario B, more than 35% of affected rail passengers do not be taken away by other transit modes until the disrupted rail service resumes normal operations. This results in an excessively long waiting time for passengers detained at interrupted stations. As can be seen, the maximum stranded time reaches 58.97 min. Furthermore, the distribution of stranded time is depicted in Figure 12. Passengers with stranded time exceeding 10 min account for approximately 3% in Scenario A, while the proportion in Scenario B is about 47%. Additionally, even more than 3000 passengers detained for over 30 min. Obviously, the alternative transit modes are not able to evacuate affected rail passengers efficiently and effectively in Scenario B, indicating that other emergency buses are in urgent demand. (2) Estimation of the stranded time and stranded rail passenger flow As illustrated in Table 9, in Scenario B, more than 35% of affected rail passengers do not be taken away by other transit modes until the disrupted rail service resumes normal operations. This results in an excessively long waiting time for passengers detained at interrupted stations. As can be seen, the maximum stranded time reaches 58.97 min. Furthermore, the distribution of stranded time is depicted in Figure 12. Passengers with stranded time exceeding 10 min account for approximately 3% in Scenario A, while the proportion in Scenario B is about 47%. Additionally, even more than 3000 passengers detained for over 30 min. Obviously, the alternative transit modes are not able to evacuate affected rail passengers efficiently and effectively in Scenario B, indicating that other emergency buses are in urgent demand.

Impact Estimation Results for Original Bus Passengers
As anticipated, with a large number of rail passengers shifting to the bus system, the original bus passengers are more or less influenced. As seen in Table 10, bus passengers

Impact Estimation Results for Original Bus Passengers
As anticipated, with a large number of rail passengers shifting to the bus system, the original bus passengers are more or less influenced. As seen in Table 10, bus passengers in Scenario B are affected to a greater extent than those in Scenario A in terms of the estimation results of delay time and the affected bus passenger flow. Specifically, the average time passengers wait to board buses increases by about 11 min in Scenario B. The affected bus passengers number over 1000, taking up approximately 57% of all boarding passengers in the interrupted section. Moreover, there are still 577 bus passengers who are unable to board buses and are stranded at bus stops until the disrupted rail service is recovered. In contrast, there are no bus passengers detained in Scenario A. Figure 13 presents the delay time distribution of the original bus passengers. The results further confirm the fact that bus passengers in Scenario B are greatly affected. The bus passengers with a delay time of more than 10 min account for about 34%. Additionally, there are over 100 bus passengers delayed over 30 min.  To further investigate the impacts on bus passengers at different bus stops, Tables 11 and 12 are plotted illustrating the passenger waiting time with passenger arrival time at stops for Scenarios A and B, respectively. In Scenario A, the waiting time of most passengers is within the bus headway except for passengers at the transfer station S10. Comparatively, for passengers at more stops, including the transfer station, turnover stations and intermediate stations, wait times for buses exceeds bus headway in Scenario B, especially for passengers at stops in the upstream direction. These most affected bus stops need to be paid more attention by transit agencies. Note that in these bus stops, the passenger waiting time follows the rule that it gradually increases with the time of arrival and then displays a linear decrease after it reaches the maximum value. It is because the continuous accumulation of arriving passengers results in incremental waiting time for passengers. Passengers whose arrival time is prior to the inflection point can board buses before cov re er T . Their wait time gradually increases. Passengers whose arrival time is after the inflection point are finally detained at stops until the disrupted rail service is recovered. Their waiting time gradually declines based on the assumption that the affected rail passengers would return to the rail system after cov re er T and the detained bus passengers are able to board the first arriving buses. It also warns us that with regard to the bus stop with a large number of stranded passengers, the passengers reaching the stop after cov re er T may still be affected. To further investigate the impacts on bus passengers at different bus stops, Tables 11 and 12 are plotted illustrating the passenger waiting time with passenger arrival time at stops for Scenarios A and B, respectively. In Scenario A, the waiting time of most passengers is within the bus headway except for passengers at the transfer station S 10 . Comparatively, for passengers at more stops, including the transfer station, turnover stations and intermediate stations, wait times for buses exceeds bus headway in Scenario B, especially for passengers at stops in the upstream direction. These most affected bus stops need to be paid more attention by transit agencies. Note that in these bus stops, the passenger waiting time follows the rule that it gradually increases with the time of arrival and then displays a linear decrease after it reaches the maximum value. It is because the continuous accumulation of arriving passengers results in incremental waiting time for passengers. Passengers whose arrival time is prior to the inflection point can board buses before T recover . Their wait time gradually increases. Passengers whose arrival time is after the inflection point are finally detained at stops until the disrupted rail service is recovered. Their waiting time gradually declines based on the assumption that the affected rail passengers would return to the rail system after T recover and the detained bus passengers are able to board the first arriving buses. It also warns us that with regard to the bus stop with a large number of stranded passengers, the passengers reaching the stop after T recover may still be affected. accumulation of arriving passengers results in incremental waiting time for passengers. Passengers whose arrival time is prior to the inflection point can board buses before cov re er T . Their wait time gradually increases. Passengers whose arrival time is after the inflection point are finally detained at stops until the disrupted rail service is recovered. Their waiting time gradually declines based on the assumption that the affected rail passengers would return to the rail system after cov re er T and the detained bus passengers are able to board the first arriving buses. It also warns us that with regard to the bus stop with a large number of stranded passengers, the passengers reaching the stop after cov re er T may still be affected. accumulation of arriving passengers results in incremental waiting time for passengers. Passengers whose arrival time is prior to the inflection point can board buses before cov re er T . Their wait time gradually increases. Passengers whose arrival time is after the inflection point are finally detained at stops until the disrupted rail service is recovered. Their waiting time gradually declines based on the assumption that the affected rail passengers would return to the rail system after cov re er T and the detained bus passengers are able to board the first arriving buses. It also warns us that with regard to the bus stop with a large number of stranded passengers, the passengers reaching the stop after cov re er T may still be affected.

Conclusions
Evaluating the impacts of unplanned urban rail disruptions on the transportation system enables transit agencies to determine the severity of the disruptions and adopt effective measures to mitigate the negative effects and ensure passengers' safety. Previous studies paid much attention to evaluating the impacts on the rail system itself, but little attention to the potential effects on the entire public transport system. To fill the gap, this study proposes a simulation-based approach to estimate the impacts on public transport travelers, including rail travelers as well as bus travelers. Firstly, we present a method of identifying affected rail passengers and estimating the scale and origin-destination distribution of affected passengers. Then a multi-agent simulation model is developed to simulate the multimodal transport evacuation process. Original bus passengers are integrated into the simulation system and interact with the rail passengers. On the basis, some indicators evaluating the impacts on rail passengers and bus passengers are presented respec-

Conclusions
Evaluating the impacts of unplanned urban rail disruptions on the transportation system enables transit agencies to determine the severity of the disruptions and adopt effective measures to mitigate the negative effects and ensure passengers' safety. Previous studies paid much attention to evaluating the impacts on the rail system itself, but little attention to the potential effects on the entire public transport system. To fill the gap, this study proposes a simulation-based approach to estimate the impacts on public transport travelers, including rail travelers as well as bus travelers. Firstly, we present a method of identifying affected rail passengers and estimating the scale and origin-destination distribution of affected passengers. Then a multi-agent simulation model is developed to simulate the multimodal transport evacuation process. Original bus passengers are integrated into the simulation system and interact with the rail passengers. On the basis, some indicators evaluating the impacts on rail passengers and bus passengers are presented respec-

Conclusions
Evaluating the impacts of unplanned urban rail disruptions on the transportation system enables transit agencies to determine the severity of the disruptions and adopt effective measures to mitigate the negative effects and ensure passengers' safety. Previous studies paid much attention to evaluating the impacts on the rail system itself, but little attention to the potential effects on the entire public transport system. To fill the gap, this study proposes a simulation-based approach to estimate the impacts on public transport travelers, including rail travelers as well as bus travelers. Firstly, we present a method of identifying affected rail passengers and estimating the scale and origin-destination distribution of affected passengers. Then a multi-agent simulation model is developed to simulate the multimodal transport evacuation process. Original bus passengers are integrated into the simulation system and interact with the rail passengers. On the basis, some indicators evaluating the impacts on rail passengers and bus passengers are presented respec-

Conclusions
Evaluating the impacts of unplanned urban rail disruptions on the transportation system enables transit agencies to determine the severity of the disruptions and adopt effective measures to mitigate the negative effects and ensure passengers' safety. Previous studies paid much attention to evaluating the impacts on the rail system itself, but little attention to the potential effects on the entire public transport system. To fill the gap, this study proposes a simulation-based approach to estimate the impacts on public transport travelers, including rail travelers as well as bus travelers. Firstly, we present a method of identifying affected rail passengers and estimating the scale and origin-destination distribution of affected passengers. Then a multi-agent simulation model is developed to simulate the multimodal transport evacuation process. Original bus passengers are integrated into the simulation system and interact with the rail passengers. On the basis, some indicators evaluating the impacts on rail passengers and bus passengers are presented respec-

Conclusions
Evaluating the impacts of unplanned urban rail disruptions on the transportation system enables transit agencies to determine the severity of the disruptions and adopt effective measures to mitigate the negative effects and ensure passengers' safety. Previous studies paid much attention to evaluating the impacts on the rail system itself, but little attention to the potential effects on the entire public transport system. To fill the gap, this study proposes a simulation-based approach to estimate the impacts on public transport travelers, including rail travelers as well as bus travelers. Firstly, we present a method of identifying affected rail passengers and estimating the scale and origin-destination distribution of affected passengers. Then a multi-agent simulation model is developed to simulate the multimodal transport evacuation process. Original bus passengers are integrated into the simulation system and interact with the rail passengers. On the basis, some indicators evaluating the impacts on rail passengers and bus passengers are presented respec-

Conclusions
Evaluating the impacts of unplanned urban rail disruptions on the transportation system enables transit agencies to determine the severity of the disruptions and adopt effective measures to mitigate the negative effects and ensure passengers' safety. Previous studies paid much attention to evaluating the impacts on the rail system itself, but little attention to the potential effects on the entire public transport system. To fill the gap, this study proposes a simulation-based approach to estimate the impacts on public transport travelers, including rail travelers as well as bus travelers. Firstly, we present a method of identifying affected rail passengers and estimating the scale and origin-destination distribution of affected passengers. Then a multi-agent simulation model is developed to simulate the multimodal transport evacuation process. Original bus passengers are integrated into the simulation system and interact with the rail passengers. On the basis, some indicators evaluating the impacts on rail passengers and bus passengers are presented respec-

Conclusions
Evaluating the impacts of unplanned urban rail disruptions on the transportation system enables transit agencies to determine the severity of the disruptions and adopt effective measures to mitigate the negative effects and ensure passengers' safety. Previous studies paid much attention to evaluating the impacts on the rail system itself, but little attention to the potential effects on the entire public transport system. To fill the gap, this study proposes a simulation-based approach to estimate the impacts on public transport travelers, including rail travelers as well as bus travelers. Firstly, we present a method of identifying affected rail passengers and estimating the scale and origin-destination distribution of affected passengers. Then a multi-agent simulation model is developed to simulate the multimodal transport evacuation process. Original bus passengers are integrated into the simulation system and interact with the rail passengers. On the basis, some indicators evaluating the impacts on rail passengers and bus passengers are presented respectively. Finally, a case study based on the Ningbo URT network is conducted and some practical results are acquired.
Whether simulation or impact assessment, this study only focuses on passengers in the interrupted section, since they are influenced by the disruptions most immediately and to the greatest extent. As analyzed in Section 2, the affected rail passengers can spread all over the URT network. Correspondingly, affected bus passengers are also not limited to the interrupted section. If the research scope is expanded to the whole public transport network, much more bus lines and bus passengers would be added to the simulation system. The issue becomes more complicated, and it is our ongoing work. Concurrently, the scale of rail disruptions in the study can be enlarged to an entire line or even a partial network.