Quantification of the Hysteresis of Macroscopic Fundamental Diagrams and Its Relationship with the Congestion Heterogeneity and Performance of a Multimodal Network

Recent studies have observed hysteresis loops in the macroscopic fundamental diagram (MFD). In particular, for the same network density, higher network flows occur during congestion onset than during congestion offset. To evaluate management strategies using the MFD, investigating the relationship between the size of these loops and network performance is needed. The existing literature has mainly discussed correlating loop width (difference in density) and height (capacity drop) with congestion heterogeneity, but has failed to prove a relationship between the capacity drop and traffic conditions. Moreover, quantification of the MFD loop in complex multimodal networks has not been investigated. The objective of this paper covers these aspects. We simulated the Sioux Falls network with different mode-share ratios (car and bus users) based on a multi-agent simulation, MATSim. We investigated the relationships between MFD loop size and congestion heterogeneity (standard deviation of density) and network performance (average passenger travel time), and found that both were directly correlated with loop width, while weakly correlated with loop height. Moreover, we divided the MFD loop into two parts according to congestion onset and offset periods and found that the heights of the two parts had opposite effects. Accordingly, we show why the relationship between capacity drop and congestion heterogeneity is not found in the literature. We also found that network performance inversely affected the height of part of the loop while the height of its other part increased with an increase in congestion heterogeneity. These results help to evaluate network performance in the presence of MFD hysteresis, leading to elaborated management decisions. How to cite this paper: Hemdan, S., Wahaballa, A.M. and Kurauchi, F. (2018) Quantification of the Hysteresis of Macroscopic Fundamental Diagrams and Its Relationship with the Congestion Heterogeneity and Performance of a Multimodal Network. Journal of Transportation Technologies, 8, 44-64. https://doi.org/10.4236/jtts.2018.81003 Received: November 28, 2017 Accepted: December 25, 2017 Published: December 28, 2017 Copyright © 2018 by authors and Scientific Research Publishing Inc. This work is licensed under the Creative Commons Attribution International License (CC BY 4.0). http://creativecommons.org/licenses/by/4.0/

between the size of these loops and network performance is needed.The existing literature has mainly discussed correlating loop width (difference in density) and height (capacity drop) with congestion heterogeneity, but has failed to prove a relationship between the capacity drop and traffic conditions.Moreover, quantification of the MFD loop in complex multimodal networks has not been investigated.The objective of this paper covers these aspects.We simulated the Sioux Falls network with different mode-share ratios (car and bus users) based on a multi-agent simulation, MATSim.We investigated the relationships between MFD loop size and congestion heterogeneity (standard deviation of density) and network performance (average passenger travel time), and found that both were directly correlated with loop width, while weakly correlated with loop height.Moreover, we divided the MFD loop into two parts according to congestion onset and offset periods and found that the heights of the two parts had opposite effects.Accordingly, we show why the relationship between capacity drop and congestion heterogeneity is not found in the literature.We also found that network performance inversely affected the height of part of the loop while the height of its other part increased with an increase in congestion heterogeneity.These results help to evaluate network performance in the presence of MFD hysteresis, leading to elaborated management decisions.

Introduction
Traffic congestion is a severe matter worldwide as it causes delays and air pollution.Thus, several traffic-management strategies have been developed to control traffic congestion.However, evaluating the efficiency of any traffic management strategy prior to its implementation is a challenge due to the wide variety of control variables that make the joint optimization of all traffic-flow measures difficult.This is, for instance, due to the large data sets that need to be processed and the large solution space needed.Alternatively, the macroscopic fundamental diagram (MFD), proposed by [1] and [2], may have the potential to provide a simpler way to assess traffic control impacts on transportation networks.MFD represents a network's traffic operations by a scatter plot relating the average network density and flow over a certain period.Instead of a detailed description of the traffic state, one variable, the number of vehicles in an area, may help to describe the traffic state.Geroliminis and Daganzo [2] showed that the concept worked in practice for road traffic in the center of a major city.It was claimed, though, that the concept of the MFD worked only for homogeneous traffic conditions.
Buisson and Ladier [3] used empirical observations of different freeways and showed that the MFD may have a significant amount of scatter; that is, multiple flows may be observed for a given value of density, and the flows during the onset of congestion differ significantly from those during the dissipation of congestion (later referred to as hysteresis).Mazloumian et al. [4] found similar MFD hysteresis loops in a simulation of an urban grid network.They found, for different estimations of spatial distribution of congestion (represented by the standard deviation of the densities), that the MFD had different shapes.Subsequently, many researchers (e.g., [5] [6] [7] [8]) have confirmed the presence of MFD hysteresis loops using empirical data and simulation results.Gayah and Daganzo [8] tested the hysteresis on MFD in a simplified two-bin network and concluded that the existence of a MFD hysteresis loop can be expected even on a symmetric network with uniform demand.Thus, in real-world conditions of city-scale networks, a hysteretic type of MFD may exist.Qian [9] simulated installing ramp meters to restrict vehicle entry to a highway when its flow reached a certain threshold.Qian concluded that hysteresis did not affect recognition of the patterns of a MFD and that the critical density and the maximum flow could be recognized readily.Nevertheless, few research efforts have been made concerning the relationship between MFD hysteresis loops and traffic conditions, as detailed in Section 2.3.Hemdan et al. [10] used MFD to qualify the effect of route, Journal of Transportation Technologies mode, and departure-time choices on multimodal network performance based on simulation results.However, a quantitative evaluation of MFD loops was not provided.Quantifying MFD hysteresis loops and correlating them to different traffic characteristics may help in comparing the effects of different management strategies on network performance, leading to more appropriate management decisions.
Recently, Saberi and Mahmassani [11], using loop detector data on freeways, found that the heterogeneous spatial distribution of network occupancy was independent of both the width of a hysteresis loop (difference in density) and its height (capacity drop); however, it was correlated with the area of the hysteresis loop.Shi and Lin [12] used loop detector data from urban expressways and found that the mean speed correlated with the width of the hysteresis loop.
However, consistent with the results of [11], it was not correlated with capacity drop.Orfanou et al. [13] correlated capacity drop with hysteresis loop shape based on vehicle trajectory data.They also found that the capacity drop could not fully describe the form or duration of the hysteresis loop.These studies and others tried to illustrate the main characteristics of MFD hysteresis loops; however, they failed in interpreting the loose relationship between capacity drops and the spatial heterogeneity of density.In this study, we continue the previous work on quantifying and characterizing MFD hysteresis loops.We attempted to find a relationship between loop height/width and the heterogeneity of congestion over the network.Additionally, the hysteresis loop was correlated with a network performance indicator, represented by the average passenger travel time consumed.We illustrate the reasons why capacity drop appearing in the MFD loop and spatial heterogeneity of density are insignificantly correlated in the literature and how they can be correlated.
The next section reviews the relevant literature with an overview of using MFDs as an indicator of network performance under different traffic control strategies.Section 3 illustrates the simulation tool used in this paper and the development of MFDs for different network characteristics and travel choices.Section 4 presents the results of the proposed quantification of the MFD hysteresis loop and its relationships.Section 5, finally, presents a summary and conclusions.

Literature Review
The MFD, developed earlier by [14] and reintroduced by [1] [2] [15] [16] is a representation of a network's traffic operations by a scatter plot for the average flow on a network and its average density.Because our focus is on comparing network performance under different control strategies in the presence of a MFD hysteresis loop, the next subsections summarize the relevant research.

Control Strategy Evaluation Based on MFDs
MFDs have been used to represent traffic evolution and to evaluate network Journal of Transportation Technologies performance.Geroliminis and Daganzo [2], using data collected from loop detectors and GPS-equipped taxis in Yokohama, proved the existence of a MFD for a city network.Later, many researchers investigated MFDs for urban networks and freeways, based on simulations (e.g., [17] [18] [19]), and empirical data (e.g., [2] [3] [6] [20]).
MFDs have been used to evaluate the performance of transport networks subject to a zone-based routing system and cordon pricing (e.g., [21] and [22]).
Zheng et al. [23] studied a time-dependent area-based pricing and have updated pricing value based on the critical density appear in MFD.Their focus was considering users heterogeneity and the traveler's behavioral adaptation to the pricing values.Zhang et al. [24] studied the impacts of various adaptive traffic signaling methods on signalized intersections with two-way arterial streets, while [19] performed a simulation on an idealized grid network.They found that adaptive traffic signals showed consistent MFDs when the network was moderately congested.Ji et al. [17] evaluated the effect of providing four ramp meters in one direction on the "A10 West" highway, Amsterdam.They found that ramp metering had a direct effect on MFD shape.The reduced inflow resulted in higher flows and lower accumulation with ramp metering, whereas without ramp metering, accumulation increased significantly while flow decreased (resulting in heavy congestion).Qian [9] simulated installing ramp meters on urban roads feeding the A10 highway to restrict vehicle feeding when the flow on the A10 reached a certain threshold.It was found that, for the A10 highway, the maximum density during the congestion period decreased, while the maximum flow value remained the same.For urban roads, congestion branches did not appear on the MFD before the implementation of ramp metering.However, the curve with ramp metering showed flow reduction at high densities, and the maximum density increased.Qian concluded that hysteresis did not restrict recognizing MFD patterns and that important values, such as the critical density and the maximum flow, could be recognized readily from MFDs.He suggested that the size of the drop in MFD loops when congestion resolved (the difference between the values of flow accompanying a certain value of density) could indicate the severity of congestion.
Other researchers (e.g., [25] [26] [27]) have used MFDs to analyze multimodal network performance.Zheng and Geroliminis [25] used MFDs to evaluate the effects of road space distribution on the performance of a congested multimodal network by optimization and simulation.The objective of their optimization was to minimize the total passenger hours traveled (PHT) by reallocating road space among modes.A simulation was performed on a bi-modal city case study to investigate the performance under a dynamic space allocation strategy.They concluded that PHT could be minimized by a dynamic space allocation strategy that used the bus lane space better during the off-peak period and served more passengers during the peak period.Geroliminis et al. [27] developed a three-dimensional MFD (3D-MFD) for a bi-modal network based on simulation Journal of Transportation Technologies data.They related the number of cars and buses to the total flow in the network and found that the network's vehicle flow decreased monotonically with the number of buses serving in the network, whereas passenger throughput was maximized at a non-zero accumulation of buses.
In summary, analyses using MFDs for evaluating the effects of different traffic management strategies suggest that networks can be evaluated based on this approach.According to changes in critical density and maximum flow values, the performance of networks can be qualified under different traffic management strategies.Additionally, the size of the drop in MFD hysteresis loops when congestion resolves may reflect congestion severity.

Hysteresis Phenomena
Early investigation of hysteresis phenomena belongs to [28] and [29].They illustrated the hysteresis of the speed-density relationship according to acceleration and deceleration.Detailed illustration of traffic hysteresis characteristics can be found in [30].Recently, Ji et al. [17] performed a simulation of a combination of an urban arterial and an urban freeway and found that during congestion onset and offset, the MFD was scattered.Buisson and Ladier [3] used empirical data on Toulouse and showed that the MFD might have a hysteresis.They summarized the following possible reasons for such hysteresis: • Differences between highways, surface streets, city center, and suburban roads.
• When demand is not distributed uniformly, congestion may occur in some parts of the network.
• The appearance and disappearance of congestion can also be a reason for successive values of mean flow.
• Differences in data-measurement locations: some loops are located near traffic signals, and some are much farther away.
For urban networks, Mazloumian et al. [4] showed by simulation that the variance in density over different locations was an important factor when determining total network performance.They concluded that, "An inhomogeneity in the spatial distribution of car density increases the probability of spillover, which substantially decreases the network flow."Later, Geroliminis and Sun [6] confirmed the findings of [4] using empirical data from the Yokohama metropolitan area.Daganzo et al. [7] showed that turning at intersections was among the key reasons for the drop in the performance with unevenly spread congestion.Gayah and Daganzo [8] tested the effects of driver adaptivity on MFD hysteresis for a homogeneous square-grid network simplified as a two-bin network.They concluded that, when congestion resolved, it did not do so instantaneously at all locations; however, it did resolve completely from one side of the queue, causing clockwise hysteresis loops in the MFD.They also concluded that these clockwise hysteresis loops should be expected, even under conditions of perfect network symmetry with uniform demand, and that they will diminish when drivers avoid congestion by an adaptive re-routing process.They suggested that the presence of highly adaptive drivers might explain why no hysteresis loop existed in Yokohama, and they expected the Yokohama network to exhibit a clockwise hysteresis loop if it was subjected to disturbances.Knoop et al. [21] concluded that congestion resolution would increase the spatial variance of the accumulation and thus (relatively) decrease the performance.To consider the spatial inhomogeneity, [31] proposed a generalized MFD (GMFD) that related the average flow to both the average density and the inhomogeneity of density (the standard deviation of the densities in different links).They showed that weighted flow was a continuous function of the weighted density and the inhomogeneity of density and they described hysteresis patterns using this function.Their results were consistent with many other studies (e.g., [4]) and MFDs may have a scatter in which multiple flows may be observed for a given value of density.
In response to the appearance of MFD hysteresis loops, some researchers had developed clustering algorithms to divide the network into homogeneous sub-networks to get low scattered MFD (e.g., [32] [33]).Overview of relevant research can be found in [34].Qian [9] suggested that the size of the drop in MFD loops when congestion resolves (the difference between the values of flow accompanying a certain value of density) could reflect congestion severity.Additionally, [21] concluded that the hysteresis loops in a MFD were an effect by themselves.Nevertheless, there has been limited research concerning the relationship between MFD hysteresis loops and traffic conditions, as detailed below.

Characterization of MFD Hysteresis Loops
Saberi and Mahmassani [20] investigated the effect of the spatial and temporal distribution of congestion on the form and characteristics of MFD hysteresis loops based on loop detector data for Portland and found a scattered pattern in the MFD when the spatial distribution of link densities was inhomogeneous.
Also, [8] characterized hysteresis loops by their direction of formation: "clockwise", "counter-clockwise", and "figure-eight".Later, [11] investigated hysteresis and capacity drop phenomena in three different freeway networks in the US using loop detector data.They identified two types of capacity drop: capacity drop while the network was loading, due to the network's inability to sustain its throughput at its peak value for a relatively long time; and capacity drop associated with network instability during a reloading process such as the afternoon loading period after an incomplete recovery from the morning peak period.
They observed two loop shapes on both freeways.A loop was associated with unstable recovery of a network in which the average network flow constantly decreased when the average network occupancy decreased.Another loop shape was observed when the network was recovering, and the average network flow remained roughly unchanged when the average network occupancy decreased.
They tried to correlate the size of the hysteresis loop and the standard deviation of network occupancy during the recovery phase.They found a very weak cor-Journal of Transportation Technologies relation between the standard deviation of network occupancy and both the width and height of the hysteresis loop.However, they found that the hysteresis loop area correlated with a congestion homogeneity index (the difference between the number of congested and uncongested links divided by the total number of links in the network) with a correlation coefficient of 0.77.Shi and Lin [12] used loop detector data to investigate the characteristics of MFDs for an urban expressway network in Shanghai, China.They explored factors that affected the characteristics of hysteresis loops.Using two days' data, resulting in four different hysteresis loops, they found a direct correlation between the width of the hysteresis loop and the mean speed, meaning that when congestion of the network was more severe, wider hysteresis loops formed.However, they also found no correlation between the capacity drop of a hysteresis loop and the mean speed.
Although these studies ( [11] and [12]) used different variables (standard deviation of network occupancy and mean speed, respectively), there was a direct inverse relationship between these variables when looking at a specific regime (loading or recovery period).In the loading process, the flow and density increased rapidly, the standard deviation of density was relatively low, and the average speed was relatively high (congestion had not happened).When the congestion dissipated, the densities of many links decreased.However, some high-density links, still in congestion, led to higher standard deviation of density among links.Moreover, flows on high-density links were low with congestion and the flows of low-density links were definitely low.Meanwhile, the average network speed in the recovery period was lower because of congestion, which also led to lower network flow.This means that both studies reflect the same result, that capacity drop was independent of the traffic situation.
Orfanou et al. [13] correlated capacity drop, acceleration, and spacing to hysteresis shape based on vehicle trajectory data observed by video cameras on the I-80 freeway in San Francisco.They identified determinants that affected the shape of hysteresis and related them to two driver-behavior classes: aggressive and timid.They classified the loops into two categories and concluded that aggressive behavior by the following vehicle resulted in a counter-clockwise hysteresis loop, whereas timid behavior resulted in a clockwise hysteresis loop.The spacing and acceleration of the following vehicle at the end of the hysteresis were found to be dominant, while the capacity drop between the onset and offset of hysteresis was the least influential factor and could not fully describe the hysteresis loop in terms of its form or duration.Recently, [35] used MFDs to compare traffic performance on various idealized hierarchical urban networks consisting of local and arterial streets using microsimulation.They differentiated arterials from local streets by additional travel lanes and additional green time at intersections.Using the size and shape of hysteresis loops, they found that the presence of arterials had a significant impact on the spatial distribution of congestion in the network.They subjectively compared networks based on the width and height of hysteresis loops (a superior network has a smaller width and height).However, they did not investigate the quantitative relationship between the size of hysteresis loops and congestion inhomogeneity.
Like MFD hysteresis loop, [36] proved the existence of a hysteresis loop in the relationship between average travel time and its variance (for journeys with the same average travel time, the variance for departures as congestion dissipates are higher than the variance for departures during the onset of congestion).Few research efforts studied the relationship between MFD loops and travel time variability loops.Gayah et al. [37] and Yildirimoglu et al. [38] studied the average travel time and its day-to-day variability relationship and analyzed the observed hysteresis loops.However, they did not investigate the relationship between the sizes of the two loops.Some promising correlations have found by [39].They found strong correlations between loops sizes and their components, specifically the differences in flow, density, and the average and standard deviation of the travel time.
In summary, this review shows that characterizing MFD hysteresis loops and their relationships to traffic characteristics is still an ongoing research topic.Few methodologies have been employed and different assumptions have been made regarding the formation of hysteresis loops and the factors that influence them.
Little attention has been paid, however, to quantifying MFD hysteresis loops and their relation with congestion inhomogeneity, network performance, and other traffic characteristics.Those who have attempted to do so have used limited sample sizes and studied specific loop shapes and have found, in the main, relationships between loop size and heterogeneity, not network performance.Additionally, the literature has provided no relation between the capacity drop on hysteresis loops and traffic conditions.Moreover, quantifying hysteresis loops in complex multimodal networks has not been investigated.Closing these gaps in the literature is the objective of this paper.

Developing MFDs for Different Multimodal Network
Characteristics and Mode Choices

Simulation of a Multimodal Network
Vehicle and passenger travel behavior can be estimated based on the MATSim framework, which uses an iterative approach for agent-based dynamic traffic assignment.Integrated simulation of private and public transport based on a queuing model allows time-dependent calculation of travel times, accounting for spillover effects and the direct interaction of private and public transport.Agents alter their behavior from iteration to iteration, based on a co-evolutionary algorithm to try to find optimal routes, modes, and departure times and, thus, maximize the total utility of their daily activity schedule.Following each iteration of the queue-based network assignment, the choices of each agent are evaluated and scored, allowing agents to select more successful options for execution of their schedule in the next iteration.The selection of travel alternatives from the choice set of each agent is performed based on a random utility model, which, after a number of iterations, leads to the convergence of individual and total utilities and, thus, to an agent-based SUE [40].The extended Sioux Falls network proposed by [41] was used.The MATSim simulation was performed using the integrated public transport simulation, in which buses share road space with cars and are affected by the congestion.Individual mode choices for agents can be explored and used for MFD drawing.
The extended Sioux Falls network [41]

Characteristics of Multimodal Network and Mode-Choice Scenarios
We tested different supply characteristics and different mode-choice combinations to consider a variety of MFD loops.The supply characteristics of scenarios 1, 2, and 3, shown in Table 1, represent homogeneous networks where all links have the same number of lanes, with three different capacity levels.Different mode-choice combinations are developed in Scenarios 4, 5, and 6.The simulation framework presented here enables mode choices for some individuals: some agents can change between using a car, bus, or walking.In this work, we considered only two modes: cars and buses.We simulated a whole day period from  referring to [41] and the marginal disutility for the car was varied between 0 and −1.The different mode-share ratios are shown in scenarios 4, 5, and 6 (Table 1).
Rieser et al. [42] concluded that less flow fluctuation and a more stable and faster convergence process could be achieved by allowing only some of the population to generate new plans.He stated that a typical set-up in MATSim is to allow re-planning for 10% of the agents.Mode choice re-planning probabilities were thus permitted for 10% of agents in this paper, and agent selections were based on agent-based SUE, as illustrated in Section 3.1.Bus fleet size and headway remained constant for all cases.The MFD represents the relationship between the average network flow and the average network density.The MFDs for different scenarios are shown in Figure 1.We used the MATSim outputs to calculate the values of the average network flow, the average network density, the standard deviation of density among links and the average passenger travel time using the following equations.showed almost the same correlation trend and accuracy as with the actual one.
Thus, for simplicity, we determined the next set of relationships using the approximated area size (by multiplying the width and height of the loop).Different relationships were found between the average standard deviation of density (for the recovery phase only and considering both the loading and recovery phases) and loop metrics.Table 2 shows the R 2 values for those relationships.Table 4 shows R 2 values for the relationships between the average passenger travel time and loop metrics.A direct correlation was found between the average passenger travel time and loop size and width.Also, a weak relationship was found between the average passenger travel time and overall loop height, which is the same result as the relationship between the standard deviation of density and whole loop height (Table 2).However, by dividing the loop into two parts, a clear inverse relationship was found between the average passenger travel time and the height of Part 1 of the loop (R 2 = 0.701).However, a weak correlation was found between the average passenger travel time and the height of Part 2 of the loop (R 2 = 0.510) (Figure 6).Comparing Figure 4 and Figure 6 shows that   flow correlated with an increase in ∆q 1 .While ∆q 2 occurs during congestion dissipation, as some link densities decrease and some links are still in congestion, this leads to a higher standard deviation of density among links that is correlated with the ∆q 2 increase.From Figure 5, increasing the number of car users increases the density and severity of congestion, as indicated by the increase in the width of the MFD loop.The scenarios of "78% Cars" and "65% Cars" (where 78% and 65% of travelers use cars, respectively) lead to more severe congestion and larger ∆k 1 .Consequently, a higher average passenger travel time occurs.However, ∆q 1 decreases according to its inverse relationship to average passenger travel time; this is accompanied by a decrease in traffic flow and a large increase in density (large ∆k 1 ).This increase in density can be eliminated by reducing the number of car users.The scenario of "58% Cars" results in light congestion, indicated by a decrease in traffic flow and a small increase in density, accompanied by an increase in ∆q 1 and, consequently, a lower average passenger travel time.An additional decrease in the percentage of car users, while keeping the bus fleet the same, would allow the network to operate with nearly free flow: see the results for the "52% Cars" and "55% Cars" modes.

Summary and Conclusions
MFDs that relate the average network density and flow have been used in the literature for network performance evaluation.However, it is difficult to attain low scattered MFDs in real-world city-scale networks: instead, recent studies have suggested more complex MFDs that exhibit a hysteresis loop.In particular, for the same network density, a higher network flow occurs during congestion onset than during congestion offset.Investigating the relationship between these loop sizes and network performance is of great importance in comparing the effect of different management strategies.This study investigated quantifying MFD hysteresis loops and the relationship between the size of this loop and both congestion heterogeneity and network performance in multimodal networks.
We estimated MFDs for traffic conditions associated with different mode-share ratios (car and bus users) on the Sioux Falls network, based on a multi-agent transport simulation.Mode choice was based on a random utility model with an agent-based stochastic user equilibrium.The goal was to compare traffic conditions associated with different mode-share ratios, based on MFDs, and, consequently, evaluate the effect of different management strategies on network performance in the presence of a hysteresis loop.First, we evaluated a quantification method for MFD loops.Previous research [11] and [12] has approximated hysteresis loop area as the product of its width and height.To ensure that this approximation does not affect result accuracy, we calculated the actual loops' areas for a sample of the developed MFDs and compared the results with the approximation.Second, we correlated loop dimensions with the standard deviation of density (representing congestion heterogeneity) and the average passenger travel time (as an indicator of multimodal network performance).Moreover, we decomposed the problem into its components by dividing loop area into two parts, according to congestion onset and offset times (associated with the maximum flow point and loop start points during network recovery).We tested the effect of partitioning loop area, based on the important results found by [43], who concluded that a higher spatial spread of density occurred during the starting point of the hysteresis loop.
By correlating with the standard deviation of density, the approximated area showed almost the same correlation trend and accuracy as did the actual area.We conclude that a hysteresis loop's width (difference in density), its height (capacity drop), and its area (by multiplying the loop's width and height) are representative metrics for MFD loops.Different relationships were found between the average standard deviation of density and loop metrics.It is concluded that the average standard deviation of density should be calculated over both loading and recovery periods, and not for the recovery period alone, especially for loops covering a long period.Dividing loop area into two parts according to congestion onset and offset times (upper and lower parts, respectively) provides some important insights.We consider, to some extent, that the partitioning criteria for a loop's area suggested in this paper are helpful for understanding a number of relationships that have not been found in the literature.It was found that the correlation between the standard deviation of density and the loop's lower part is stronger than its correlation with the overall loop.That is, the lower part of the loop alone can represent the heterogeneity effect.We also found that the correlation between the standard deviation of density and the heights of both loops' parts was stronger than the correlation with overall loop height.This is because the two heights have opposite effects on the standard deviation of den-

Figure 3 .
Figure 3. Correlation between (a) actual versus approximated loop area, (b) the standard deviation of density versus actual loop area, and (c) the standard deviation of density versus approximated loop area.
the correlation between the standard deviation of density and the height of the Part 2 loop is stronger than its correlation with the height of the Part 1 loop, whereas the correlation between the average passenger travel time and the height of the Part 1 loop is stronger than its correlation with the height of the Part 2 loop.This means that network performance inversely affects the height of the Part 1 loop whereas the height of the Part 2 loop increases with the increase of the spatial spread of density.Network performance inversely affects the capacity drop while the network is loading, while the inhomogeneous spatial distribution of congestion increases the flow reduction during network unloading with congestion dissipation.As discussed in detail in Section 4.2, the capacity drop in the Part 1 loop occurs due to network loading, although it is unable to sustain its flow at its maximum value for long.More loading leads to a decrease in network

Figure 6 .
Figure 6.Relationship between average passenger travel time and (a) height of part 1 loop and (b) height of part 2 loop.
consists of two road types; urban and highway roads, with a total of 334 links and 282 nodes.The bus network consists of five bus lines with a 5-min headway operation.In total, 56,904 workers commute, and 27,206 persons perform secondary activities, with a total of 168,220 daily trips.Every agent performs only one round trip per day (home -workhome or home -secondary -home).The time window for an activity and its typical duration affect how peaked or spread out the demand is in the peak-hours, and this will reflect the degree of congestion in the network.The work activity opens at 8 a.m. and closes at 6 p.m., with a typical duration of 9 h (latest work start time is 9 a.m.), while secondary trips can be taken between 8 a.m. and 8 p.m. and completed within 1 h during this period.MATSim allows different travel choices for individuals.The simulation presented here enables mode choices for some individuals according to the different traffic management strategies described in the next subsection.

Table 1 .
Scenarios of different network characteristics and different mode-choice combinations (Sioux Falls network).
a.m. to 11:59 p.m. Mode-choice rate was controlled by changing the disutility for traveling by car.The marginal disutility for the bus mode was set by

Table 4 .
Summary of R 2 values for relationships between network performance and loop metrics.