Research Article | Open Access

# Formation and Propagation of Local Traffic Jam

**Academic Editor:**Cengiz Çinar

#### Abstract

Large scale traffic congestion often stems from local traffic jam in single road or intersection. In this paper, macroscopic method was used to explore the formation and propagation of local traffic jam. It is found that (1) the propagation of traffic jam can be seen as the propagation of traffic signal parameters, that is, virtual split and virtual green time; (2) for a road with endogenous flow, entrance location influences the jam propagation. With the same demand (upstream links flow and entrance flow), the upstream got more influence; (3) when a one-lane road is thoroughly congested, virtual signal parameters everywhere are the same as that at stop line; for a basic road, the virtual signals work in a cooperative manner; (4) phase sequence is one important parameter that influences traffic performances during peak hour where spill back of channelization takes place. The same phase plan for left-turn flow and through flow would be preferred; (5) signal coordination plays an important role in traffic jam propagation and hence effective network signal parameters should be designed to prevent jam from propagation to the whole network. These findings would serve as a basis for future network traffic congestion control.

#### 1. Introduction

Urban road network has been considered as the “life line” of urban daily life. It serves as the basis of urban economy. The assurance of its normal operation is very important. However, as the development of urbanization, the number of vehicles relatively exceeds the capacity of road network which often results in large scale traffic jam during peak hours. This leads to many problems such as pollution, and noise. Large scale traffic jam often stems from local jam which takes place in a road or single node. According to Wright and Roberg [1], there are three sources of traffic jam: a temporary obstruction, stochastic fluctuation in demand, and permanent capacity bottleneck. These types of congestion cannot be prevented efficiently because they are related tightly with road geometry, land usage, or other factors which cannot be changed easily. However, if we can control the propagation of traffic jam, then, network traffic jam will be prevented. The propagation process looks like “domino phenomena” which is characterized by connected set of events. Many literatures are focusing on this subject.

Simulation methods including microscopic and macroscopic types are among the most used methods. Microscopic tools include cellular automata [2] and car following models [3]. They are based on dynamic motion of single vehicle and hence can grasp the traffic performance of the whole network [1, 4, 5]; macroscopic tools often adopted a discrete form [6, 7] of LWR model which stemmed from Lighthill and Whitham [8] and Richards [9]. This method divided links into a series of consecutive cells and time into equal time steps. By checking and tracking basic traffic parameters (flow, density, and velocity) of every cell during each step, the dynamic nature of traffic jam can be derived. This kind of method greatly relies on underlying models and the simulation networks are often regular [10, 11] which lacks generality. Furthermore, no analytical conclusion can be obtained through this approach.

Some researchers tried to establish analytical results of traffic jam. Wright and Roberg [1] proposed a simple analytical model for incident-based jam growth and discussed the effect of the length of the channelized part of roads and stop line width assignment on jam formation; Michalopoulos et al. [12] formulated a rigorous but complex mathematical modeling for signalized links; Skabardonis and Geroliminis [13, 14] developed an analytical model based on kinetic wave theory to construct monitor methods for arterials. Based on the derivation of waves produced by traffic signal, vehicle trajectory was identified and hence travel time function was formulated; Daganzo [15] divided the network into connected neighborhoods and constructed the relationship between defined traffic variables and formulated the so-called “gridlock” condition. Through these researches we can get a deep understanding of traffic jam. If we examine the general road network, we can see that it is comprised of basic road (some may be simpler) as shown in Figure 1. The road is multilane structure controlled by traffic signal. The traffic demand may come from upstream links or entrance within roads. Compared with this situation, current researches are not enough to describe traffic flow operation on basic road segment, not to mention the network traffic jam propagation.

Some scholars used data collected from control system (mainly loop detectors) to monitor long queue and spillover [16, 17]. However, under congestion hour traffic jam involves various traffic flows that include multilane flow and endogenous flow from within the entrance, and these are difficult to distinguish from loop detector data. These factors may undermine the applicability of these methods.

The dynamic nature of network traffic jam (critical condition for its formation, propagation speed, propagation route, and so on) is complex in that they are influenced by many factors such as road physical characteristics, network topology, and control parameters. Understanding of these relationships is important and premise for further control and management. However, they were not fully accounted for in traditional methods. In this paper, with the consideration that network jam always stems from traffic jam formed within single road, attention was put on the traffic jam formation on signal controlled network comprised of basic road as in Figure 1 that is generally existed. After description of critical condition of jam formation for one-way link (Section 2.2), decomposition method is used to divide the road into two multilane segment (Section 2.3). Factors include insufficient split or excess demand (Section 3.1), endogenous flow (Section 3.2) and channelization spillback (Section 3.3) are discussed. The influences of road spillback on upstream links are then analyzed (Section 4). We end with some highlights for further researches.

#### 2. Preliminary Results

##### 2.1. Notations

: optimal flow rate, : optimal density, : jam density, , , , : green time, red time, cycle and split respectively, : virtual red time, , , : stopping wave, starting wave, and wave that emerges after and meet, : maximal queue length, , , : output flow rate of left-turn flow, through flow and right-turn flow. “” means output, , , : green time for left-turn flow, through flow and right-turn flow.

##### 2.2. Wave Dynamics under Simple Condition

Consider a one-lane road controlled by signal without endogenous flow as shown in Figure 2 which is the simplest condition in the road network (such as a road regulated by one-way traffic organization), assuming that fundamental diagram of traffic flow is parabolic function as shown in Figure 3. Upstream flow is (point A in Figure 2) and from parabolic function can be obtained (density of upstream flow is always smaller than optimum density under maximum flow; that is, point A should be on left side of the curve, otherwise stopping wave speed would be bigger than starting wave which denotes that queue would not disperse forever, which can be seen from Figure 3):

**(a)**

**(b)**

Formation and dispersion of a queue behind stop line at signal controlled intersection are shown in Figure 4(a). At first (the beginning of red time), stopping wave (red lines in Figure 4 which represent queue back) propagates upstream with velocity and a queue forms. When green time begins a starting wave (blue line in Figure 4) emerges and also propagates upstream with bigger speed (denoted as ). After starting wave catches up with stopping wave and queue has dispersed and a new wave (denoted by its speed ) forms. It takes for wave to pass through stop line. If green time, , is bigger than , then it can be concluded that traffic state would reproduce cycle by cycle. However, if is smaller than , the wave propagation profile is different (as shown in Figure 4(b)) cycle by cycle. At first, stopping wave will spread with speed which makes queue back of this cycle further from stop line than former cycle. As time elapses, queue back becomes further and further, so it is an unstable condition which represents oversaturated condition.

**(a)**

**(b)**

Based on the analysis above, some formulas can be got:

When holds (we can deduce that this equation determines a specific split ) a stable state forms (if decreases, it will result in oversaturation), that is,

The critical condition is defined such that the furthest queue back reaches road tail because larger or smaller split will result spillover of queue:

Equations (4a), (4b), and (4c) establishe the relationship between link length , cycle , split , and flow . Given split can be obtained as shown in embedded figure in Figure 5. When and are given, and (also ) can be deduced as shown in Figure 5. It can be seen that links with different lengths require different control parameters and upstream inflow to maintain their normal operation. Much more green time is needed during higher demand under the same road length. This is why unchanged green time usually causes longer queue especially during peak hours. Spillover caused by smaller green time will be discussed later.

##### 2.3. Decomposition of Complex Link Topology

The problems discussed above are mainly focusing on the ideal scenario of one-way links while the case in reality shows more topological complexity. We define the basic road shown in Figure 6(a) that always appears in urban road network. It is a road controlled by signal with a channelized section of length . Overall upstream traffic demand, is assumed to be uniformly distributed between two lanes, that is, . Proportion of left-turn flow, right-turn flow, and through flow is , , and (and ), respectively. The same symbol expression method is for signal parameters, that is, , , , and , so we can get

**(a)**

**(b)**

Decompose the basic road as shown in Figure 6(b). The basic road has been divided into two sections: upstream section and channelized section. Both are comprised of some single one-way links. Channelized sections are controlled by traffic signal (except right-turn section). The analysis hence is simplified that the method in Section 2 can be used in further analysis.

#### 3. Complex Circumstances

During daily operation, circumstances often hold more complex nature. This section deals with this problem by analyzing dynamic and insufficient split, endogenous flow of basic road, and the spillover of channelized section.

##### 3.1. Dynamic and Insufficient Split

Usually oversaturation is developed gradually and queue back extends cycle by cycle, that is, domino phenomenon is the result of oversaturation extending in several cycles as shown in Figure 7(a). Each green signal tends to pull back queue back while each red signal prolongs queue back. Furthermore, the length pulled back or prolonged is proportional to time duration of traffic signal because in the figure slope of wave trajectory; that is, wave speed including starting wave and stopping wave is fixed.

Given green time of cycle (or red time ) queue length pulled back ( prolonged) can be computed from Figures 7(b1) and 7(b2):

After th red time (th green time ) location of queue back is where denotes the length of initial oversaturated queue. It can also be calculated easily by assuming that the traffic is undersaturated at first, as shown in Figure 8. Symbols and derivations are as follows:

Replace in (7) queue back can be calculated at any time. Given changeable and and upstream traffic demand , and hence each also can be determined. When queue back exceeds road length, that is, then spillover emerges ( in Figure 9 is the time when queue back exceeds road length).

However, spillover itself does not result domino phenomena because there may be no upstream traffic demand, this problem will be discussed in Section 4.

In order to describe jam propagation upstream, the method of SVS (speed of virtual signal) is used. The principle is that there is a virtual signal located everywhere. When the jam has not reached a specific location, the virtual signal is always green. When the location is occupied by a queue (i.e., jam propagates to this position), the virtual signal shows red. When the queue is cleared, the virtual signal turns to green. From this description, we can see that the propagation of jam can be transformed to the propagation of virtual signal parameter. We use split to explore the propagation of jam. If we obtain locations of , , shown in Figure 10, then we can derive the split dynamics along time and space (see Figure 11). If demand exceeds supply, at the end of each cycle, the residual vehicle number is (the queue that remains after past cycles is not included in this formula), so three locations for each cycle can be derived:

After obtaining , , for every cycle, together with Figure 11, split everywhere for any cycle can be determined.

Figures 12 and 13 present a simple numerical example for propagation of split. From Figure 12, for a specified location, the gradual change of virtual signal cycle by cycle can be easily seen. Figure 13 gives the split on space-time (expressed as cycle number) diagram. The slope can be seen as the speed of jam propagation.

In order to develop a measure of this speed, an index which incorporates stopping wave and starting wave is proposed later on. For a cycle of jam propagation, the increment of queue back can be got:

If the control parameters are constant, then we define speed of virtual signal (SVS) to express the severity of jam, that is, the propagation of congestion. During a cycle, queue back propagates , and the average speed is then

From the expression we can see that SVS is in linear relationship with split.

Figures 14 and 15 present the relationship between SVS and split, arriving flow rate, respectively. From the figure we can see that SVS is very sensitive to traffic state under saturation condition. The results above are based on parabolic fundamental, which is symmetrical about optimal flow. Figures 16 and 17 plot the SVS versus split and arriving flow under triangular fundamental diagram like Figure 3(b). The parameters are given beneath the figure.

Similar trend is also found for triangular fundamental diagram. However, it is evident from both figures that change of SVS is somewhat moderate compared with parabolic diagram. This phenomena stem from the different nature embedded in these two fundamental diagrams: wave speed for triangular diagram is constant when the density is below , but for parabolic fundamental diagram, wave speed (tangent line of specific point) is changing simultaneously with traffic demand.

Equation (12) assumes that the arriving flow is constant; however, it is not the case. But we can generalize it to dynamic condition where within one cycle the arriving is stable but different from cycle to cycle, that is, where is stop wave speed at cycle .

##### 3.2. Influence of Endogenous Flow

Traditionally it is assumed that flow is generated at nodes (i.e., intersections within the network or originations of a road network) as we did in network traffic analysis [18]. However, intersections inside road network actually do not generate flow (of course, from the viewpoint of downstream intersection, flow is “produced” at upstream intersection). In fact, most traffic stems from entrance within road (e.g., shopping malls, schools, and official departments about the road). It is necessary to take endogenous traffic flow into account. Due to complexity in analyzing multilane road, we only deal with one-way road with entrance inside which the flow is denoted as (shown in Figure 18). The basic assumption is major road priority; that is, vehicles in major road have priority to cross.

When a formed queue back has not reached entrance location, the upstream flow is ; when queue length exceeds (distance between stop line and endogenous flow entrance location) upstream flow becomes , so the stopping wave and starting wave propagate differently from that without endogenous flow (Figure 4(a)). Given and stopping wave profile and starting wave profile can be obtained using the same method. Denote stopping wave as (wave when queue length has exceeded ), (wave when queue length has not exceeded ) and , as the wave produced after and meet as shown in Figure 20.

From Figure 19 formulas below can be got:

For case I, critical condition can be formulated just as in (3a) and (3b):

For cases II and III, critical conditions as follows:

The virtual red time in Figure 20 and (16a), (16b), and (16c) is obtained from geometric relationship. Figure 21 gives the relationship between critical split and different entrance locations from (15a), (15b), (16a), (16b), and (16c). We can see that upstream location gets more influence in that once queue propagates exceed entrance endogenous flow cannot enter. Moreover, propagation of jam (i.e., SVS under this condition) also holds such trend because of the same tendency of propagation speed of stopping and starting wave.

##### 3.3. Channelized Section Spillover

Most researches do not take spillover of channelized section into account. But at peak hours due to length constraint the channelized section often cannot accommodate excessive vehicles which make spillover of channelized section inevitable. Following Wright and Roberg [1], we assume that traffic flow running on the upstream section was mixed uniformly (i.e., velocity, density, and flow are uniformly distributed in the section) and lane changing are executed instantly at interface between channelized section and upstream section. This is acceptable from the viewpoint of system description although some details we do not concern with are lost. Road layout is as that in Figure 6. From fundamental traffic theory we can easily get flow-density relationship as shown in Figure 22. Parabolic curve beneath is flow-density relationship of left-turn channelized section (right-turn or through channelized section is the same) and the upper curve is that of upstream section.

When there exists channelized section spillover the influence can be considered as “virtual red time,” principle of which is described in Figure 23: once there is spillover a queue will propagate upstream and block upstream section result of which is the same as signal control. A virtual signal is set at the interface between two sections as shown in Figure 23. Virtual red time equals to time difference between stopping wave and starting wave spreading over the interface. Now the analysis of spillover with consideration of channelized section can be taken by two steps:(I) analyze supply-demand of channelized section to obtain block time (note that queues of two channelized sections, that is, through section and left-turn section may be overflowed at the same time or may be not, so block time of these spillovers may be separated, overlapped, or overlapped partly);(II)after virtual time (interface block time) is obtained, analyze supply-demand of upstream section. If upstream section queue back exceeds road length, then upstream intersection will be blocked.

Due to the fact that both channelized sections may block the interface we first deal with only one section spillover and then two.

###### 3.3.1. One Section Spillover

Since results of spillover of the two channelized section are the same, here we simply assume that only left-turn section spillover takes place. When , that is, left-turn traffic flow equals through flow . It can be deduced from Figure 22 that starting wave and stopping wave keep its speed after spread over interface. It is condition A in Figure 24. This condition degrades to one-way link spillover which has been discussed in Section 2.

When (i.e., ), upstream section holds more flow averagely than channelized section. The wave propagation profile should be the pattern as in Figure 24. Its stopping wave speed is bigger than and is smaller than . When green time extends over , the queue will be cleared completely. From the figure critical conditions for domino formation are deduced:

When (i.e., ), upstream section holds less flow averagely than channelized section, so the wave propagation profile should be the pattern as in Figure 24. Its stopping wave speed is smaller than and is greater than . When green time extends over then the queue will be cleared completely. Critical conditions for domino formation are deduced:

From geometry and Figure 24, we can deduce virtual red time for interface when critical condition (i.e., (17) and (18)) is not reached: where is virtual red for interface. From it split can also be determined: . Since , from (19) we can see that virtual red time increases against channelization section length, and opposite relationship holds for split and channelization section length.

Now critical condition of (4a), (4b), and (4c) can be used to determine supply-demand of upstream section. Note that there is no channelization spill back for through traffic flow. Since virtual split of interface is bigger than (from (19) relationship that is evident and also ), for condition A and C in Figure 24, if critical condition of (18) has not reached, demand of upstream section is always satisfied because

For condition B, if the critical relationship of (4a), (4b), and (4c) is not hold, traffic jam will form and propagate upstream.

If channelization spillover of only one direction happens, capacity of the other two directions may be wasted, which imply the importance of well-designed phase sequences. From (4a), (4b), and (4c), it implies that the equation below should hold to prevent wasted capacity:

Equations (21a) and (21b) assure that the demand can be satisfied and (21c) guarantees queue back would not reach interface such that block cannot happen.

###### 3.3.2. Two Section Spillovers

Sometimes during peak hour left-turn demand and through demand may both exceed respective capacity which results in spillover of both sections. In such case the problem may be complicated but can be explained conveniently. Spillover of two sections is just as we set two signals in the interface at the same time. When both virtual signals show green then overall virtual signal is green otherwise the overall signal is red. One case is shown in Figure 25.

According to the signal assignment of left-turn flow and through flow, there will be three cases as shown in Figure 26. Apparently for the same split, different “offset” will produce different spillover virtual signal which implies the importance of phase sequence under spillover condition. We can see that during congestion hours (where channelized section spill back take place); it would be better to put through flow and left-turn flow in the same phase.

#### 4. Impact of Spillover on Upstream Intersection

Spillover is the precondition of congestion propagation. However, spillover itself is not enough for domino phenomena formation because if there is not enough upstream demand congestion would not spread. Thus the sufficient and necessary condition for domino phenomena is first, spillover of downstream road and second, there is continuous demand upstream.

When a queue forms and spill back and blocks upstream intersection, there are basically two events to be analyzed: first, how long it will block upstream intersection and second, its influence on upstream flows. When queue back reaches upstream intersection at time, say, (as shown in Figure 27) then it blocks the intersection till the starting wave catches up with it.

We denote that the time duration of the intersection is blocked (i.e., blocked time) as , then from the viewpoint of upstream signal phases there are three types (using a typical four-phase intersection shown in Figure 28 as an example) of block. Take phase A to describe as following.(I)Block fully happens at green time. Then the phase will lose effective green time which will decrease capacity respectively.(II)Block fully happens at red time. Since no matter whether the intersection is blocked, this movement always has not got RoW (Right of Way), so this type of block holds no impact.(III)Block happens partially at green time and partially at red time (let us neglect amber time because it will not harm the conclusion). This block will encroach some green time which will reduce some capacity.

Due to reduced capacity caused by downstream block, specific phase upstream will be under oversaturation temporarily and long queue forms which may be block upstream intersection of this crossroad.

The block influences upstream link flow differently. In Figure 28 the most possible flow blocked by downstream jam is highlighted in red line for each phase. Green time for these flows then may be wasted hence new jam may form at upstream links. This can be seen as domino phenomena of traffic jam in the whole network.

From the description above we can see that badly designed signal may contribute to network jam propagation, which suggests the importance of network signal coordination during peak hours.

#### 5. Conclusion

Due to the complexity of urban traffic jam, many researches were taken under simplified assumption or used simulated methods. In this paper, urban traffic jam formation is studied. From the research we can see that many factors (such as road geometry, flow structure, and signal settings) contribute to the traffic jam formation and propagation which should be took into account in the control of jam during peak hours.

However, limited by the macroscopic method, we cannot look into deeply mutual influence of different flow, (left-turn flow, right-turn flow, and through flow) which may contribute to the formation and propagation of traffic jam. This may be solved by using microscopic analysis methods. Furthermore, formation of congestion is only prelude of large scale urban traffic jam. Propagation of urban traffic jam (denominated as “domino phenomena” here) itself is influenced by many factors. Some further researches could be highlighted.(1)The calibration of fundamental diagram: the relationship between flow, density, and speed is a basic problem of traffic engineering field and exerts great influence on jam propagation (severity and propagation speed). A thorough examination of the influence is needed.(2)The asssumption of uniformly mixed flow at upstream section shoule be revised. In reality, the distribution may be not averagely. Whether this feature prompts the propagation f traffic jam can be included in our model with a coefficient that needs to be explored closely. (3)Isolated intersection signal design at peak hours. From this research we can see that traditional control principle of minimizing delay may be ineffective in dealing with peak hour flow due to the spillover of long queue. Till now the control strategy that takes this condition into account is absent.

#### Acknowledgments

This work was supported by the State Key Development Program for Basic Research of China (Grant no. 2012CB725402); the National Science Foundation for Postdoctoral Scientists of China (Grant no. 2012M521175); Excellent Postdoctoral Science Foundation of Zhejiang Province (Bsh1202056).

#### References

- C. Wright and P. Roberg, “The conceptual structure of traffic jams,”
*Transport Policy*, vol. 5, no. 1, pp. 23–35, 1998. View at: Publisher Site | Google Scholar - M. Schreckenberg, L. Neubert, and J. Wahle, “Simulation of traffic in large road networks,”
*Future Generation Computer Systems*, vol. 17, no. 5, pp. 649–657, 2001. View at: Publisher Site | Google Scholar - M. Treiber, A. Hennecke, and D. Helbing, “Congested traffic states in empirical observations and microscopic simulations,”
*Physical Review E - Statistical Physics, Plasmas, Fluids, and Related Interdisciplinary Topics*, vol. 62, no. 2 B, pp. 1805–1824, 2000. View at: Google Scholar - N. J. Linesch and R. M. D'Souza, “Periodic states, local effects and coexistence in the BML traffic jam model,”
*Physica A: Statistical Mechanics and its Applications*, vol. 387, no. 24, pp. 6170–6176, 2008. View at: Publisher Site | Google Scholar - T. Nagatani, “The physics of traffic jams,”
*Reports on Progress in Physics*, vol. 65, no. 9, pp. 1331–1386, 2002. View at: Publisher Site | Google Scholar - C. F. Daganzo, “The cell transmission model: A dynamic representation of highway traffic consistent with the hydrodynamic theory,”
*Transportation Research Part B*, vol. 28, no. 4, pp. 269–287, 1994. View at: Google Scholar - C. F. Daganzo, “The cell transmission model, part II: Network traffic,”
*Transportation Research Part B*, vol. 29, no. 2, pp. 79–93, 1995. View at: Google Scholar - M. J. Lighthill and G. B. Whitham, “On kinematic waves. II. A theory of traffic flow on long crowded roads,”
*Proceedings of the Royal Society. London. Series A. Mathematical, Physical and Engineering Sciences*, vol. 229, pp. 317–345, 1955. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet - P. I. Richards, “Shock waves on the highway,”
*Operations Research*, vol. 4, pp. 42–51, 1956. View at: Publisher Site | Google Scholar | MathSciNet - J. Long, Z. Gao, X. Zhao, A. Lian, and P. Orenstein, “Urban traffic jam simulation based on the cell transmission model,”
*Networks and Spatial Economics*, vol. 11, no. 1, pp. 43–64, 2011. View at: Publisher Site | Google Scholar | MathSciNet - P. Roberg and C. R. Abbess, “Diagnosis and treatment of congestion in central urban areas,”
*European Journal of Operational Research*, vol. 104, no. 1, pp. 218–230, 1998. View at: Google Scholar - P. G. Michalopoulos, D. E. Beskos, and J.-K. Lin, “Analysis of interrupted traffic flow by finite difference methods,”
*Transportation Research. Part B. Methodological. An International Journal*, vol. 18, no. 4-5, pp. 409–421, 1984. View at: Publisher Site | Google Scholar | MathSciNet - A. Skabardonis and N. Geroliminis, “Real-time estimation of travel times along signalized arterials,” in
*Proceedings of the 16th International Symposium on Transportation and Traffic Theory*, 2005. View at: Google Scholar - A. Skabardonis and N. Geroliminis, “Real-time monitoring and control on signalized arterials,”
*Journal of Intelligent Transportation Systems: Technology, Planning, and Operations*, vol. 12, no. 2, pp. 64–74, 2008. View at: Publisher Site | Google Scholar - C. F. Daganzo, “Urban gridlock: Macroscopic modeling and mitigation approaches,”
*Transportation Research Part B: Methodological*, vol. 41, no. 1, pp. 49–62, 2007. View at: Publisher Site | Google Scholar - X. Wu, H. X. Liu, and D. Gettman, “Identification of oversaturated intersections using high-resolution traffic signal data,”
*Transportation Research Part C: Emerging Technologies*, vol. 18, no. 4, pp. 626–638, 2010. View at: Publisher Site | Google Scholar - N. . Geroliminis and A. Skabardonis, “Queue spillovers in city street networks with signal-controlled intersections,” in
*Proceedings of the 89th Annual Meeting of the Transportation Research Board*, Washington, DC, USA, 2010. View at: Google Scholar - Y. Sheffi,
*Urban Transportation Networks: Equilibrium Analysis with Mathematical Programming Methods*, Prentice-Hall, Englewood Cliffs, NJ, USA, 1985.

#### Copyright

Copyright © 2013 Hong-sheng Qi et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.