Comparing Dynamic User Equilibrium and Noniterative Stochastic Route Choice in a Simulation-Based Dynamic Traffic Assignment Model: Practical Considerations for Large-Scale Networks

Abstract

Simulation-based dynamic traffic assignment (DTA) models play a vital role in transportation planning and operations. While the widely studied equilibrium-seeking DTA including dynamic user equilibrium (DUE) often provides robust and consistent outcomes, their expensive computational cost for large-scale network applications has been a burden in practice. The noniterative stochastic route choice (SRC) model, as a nonequilibrium seeking DTA model, provides an alternative for specific transportation operations applications that may not require equilibrium results after all (e.g., evacuation and major network disruptions) and thus tend to be computationally less expensive, yet may suffer from inconsistent outcomes. While DUE is a widely accepted approach for many strategic planning applications, SRC has been increasingly used in practice for traffic operations purposes. This paper aims to provide a comparative and quantitative analysis of the two modeling approaches. Specifically, a comparison has been made at two levels: link-level flows and network-level congestion patterns. Results suggest that adaptive driving improves the quality of the SRC solution, but its difference from DUE still remains significant at the link level. Results have practical implications for the application of large-scale simulation-based DTA models for planning and operations purposes.

This work is licensed under http://creativecommons.org/licenses/by/4.0/ (the “License”). Notwithstanding the ProQuest Terms and Conditions, you may use this content in accordance with the terms of the License.

1. Introduction

Simulation-based dynamic traffic assignment (DTA) models are powerful tools for planning and operations of transportation networks. The two main components of a DTA model are equilibrium-seeking route choice and network loading [1]. The literature also clearly distinguishes between two other commonly used terms known as route search and route choice. Route search refers to the filtering of the most appropriate paths from thousands or millions of possibilities in a large-scale network [2]. Route choice is, on the other hand, the selection of certain routes from an origin to a destination in the presence of alternative routes. Route selection can depend on multiple factors; however, the most widely assumed behavior is that users seek to minimize their respective travel time. Route choice is considered as the core of the traffic assignment problem commonly following Wardrop’s principle based on the assumptions that all drivers are homogeneous and completely rational, and they have perfect information on all available paths and their costs [3–7]. While equilibrium-seeking DTA models often provide robust and consistent outcomes, their expensive computational cost in large-scale networks has made their applications in practice limited and an ongoing challenge.

A noniterative stochastic route choice (SRC) model, as a nonequilibrium seeking DTA model, provides an alternative approach for specific large-scale transportation operations applications that may not require equilibrium after all such simulation network-level evacuation or major disruptions. Thus, it tends to be computationally less expensive. While DUE is a widely accepted approach for many strategic planning applications, SRC has been increasingly used in practice for traffic operations optimization because of its lower computational cost without proper acknowledgement of its different nature and use case compared to equilibrium-seeking approaches. SRC models allow drivers to choose a route from a set of feasible paths given the utility function of the path and select a path without strictly considering the flows and travel times on other paths. Consequently, the path overlapping probability increases due to lack of iterative process and can increase travel time on certain paths. This problem is often resolved by giving the option to simulated drivers to alter their selected paths in order to decrease their travel times enroute. This behavior, known as adaptive driving, can increase the average network flow and reduces the formation of gridlock in the simulation [8].

Despite the existence of a vast literature on both simulation-based DUE and SRC, very few studies have provided a quantitative comparative analysis of the two modeling approaches and how selection of each could result in significantly different outcomes in practice. A recent study by [9] highlights that care is needed when using deterministic as opposed to stochastic day-to-day models for transport planning and control purposes. They argued that the two models may provide qualitatively similar outcomes only if the mean of the stochastic model behaves similarly to the deterministic model. In this paper, we aim to provide a comparative quantitative analysis of the solution characteristics of the DUE and SRC models applied to a large-scale network of Melbourne, Australia [10]. We also demonstrate the differences between the two modeling outcomes in a simple network traffic signal optimization problem. The paper also explores the impact of adaptive driving on the quality of the SRC model at both the link and network levels. It is very well known in the academic literature that the two models are fundamentally different and serve different purposes in transportation planning and operations. However, given the increasing use of SRC in practice for different applications, in this paper we aim to demonstrate with numerical evidence how use of each of the two modeling approaches affects the estimated traffic volumes and network wide congestion patterns and provides significantly different outcomes if used carelessly for traffic operations optimization.

2. Background

2.1. Simulation-Based Dynamic Traffic Assignment Models

Traffic assignment models including simulation-based dynamic traffic assignment models are extensively used for a variety of applications such as transportation planning and traffic operations and management [4, 5, 11–13]. In the presence of a reliable origin-destination (OD) demand matrix, a reasonable estimation of network and link traffic states can often be achieved in these models. The two main components of the dynamic traffic assignment models, namely, route choice and network loading, aim to capture vehicles’ movements dynamics in a network. Network loading mainly replicates vehicles’ movements in the network on different paths specified by the route choice model.

The route choice problem is the most critical aspect of the traffic assignment problem as it involves multiple behavioral factors especially in a large-scale network. Trying to capture the reasons behind the selection of paths has always been a challenging task due to the heterogeneity of drivers’ preferences. Even under the assumption of similar behavior of drivers, the model faces many challenges such as route search from a huge set of possible routes in a large-scale network and real-time access and provision of traffic conditions. Vehicles’ path-changing enroute due to nonrecurrent events such as incidents, vehicle breakdowns, and special events requires special modeling considerations as most existing models in practice apply pre-defined routes before vehicles’ departure [14]. Different equilibrium-seeking models, such as DUE, dynamic stochastic user equilibrium (DSUE), boundedly rational user equilibrium (BRUE) [15–17], multiclass fuzzy user equilibrium (MFUE) [18], prospect-based user equilibrium (PBUE) [19], and behavioral user equilibrium (BUE) [20], adopt unique behavioral assumptions of drivers [1]. On the other hand, drivers’ route choice using a probabilistic approach considering utilities of different alternative routes results in the so-called stochastic route choice (SRC) models that are not necessarily seeking an equilibrium solution [21].

2.1.1. Dynamic User Equilibrium (DUE)

DUE aims to achieve the state of user equilibrium in the network for each time interval in the presence of time-dependent demand. The process is typically carried out in an iterative manner for each origin-destination-time (ODT) pair in the network. At each iteration, the set of shortest paths is updated, and path flows are calculated using, for instance, the method of successive averages (MSA) as one of the most widely used techniques [1, 22]. To measure the quality of the DUE solution, the relative gap metric is commonly used for each ODT that quantifies between the total cost of paths used by vehicles and the total cost of the shortest path used by all the vehicles [13]. A relative gap can be defined and expressed as follows to measure the proximity to an equilibrium solution [23]:$$\begin{array}{}\text{(1)}& \text{Rgap}n=\frac{{{\displaystyle \sum }}_t{{\displaystyle \sum }}_{r,s\in \delta }{{\displaystyle \sum }}_{p\in P_{rs}t}{f}_{rsp}^{n}t{\tau }_{rsp}^n-{\theta }_{rs}^{n}t}{{{\displaystyle \sum }}_t{{\displaystyle \sum }}_{r,s\in \delta }{d}_{rs}t{\theta }_{rs}^{n}t},\end{array}$$where $f_{rsp}^{n}t$ represents flow on path $p$ from OD pair r from s at time t and iteration n whereas ${\tau }_{rsp}^n-{\theta }_{rs}^{n}t$ quantifies the difference in cost experienced by vehicle following path with cost ${\tau }_{rsp}^n$ in comparison to the shortest path cost ${\theta }_{rs}^{n}t$, $d_{rs}t$ represents the demand at ODT pair of r, s and t, $P_{rs}t$ is the set of all available paths for the ODT, and $\delta$ is the set of all paths for the OD in network independent of time t. When the relative gap reduces below an acceptable threshold, the network in question is considered to reach the equilibrium state.

MSA is a widely used iterative method that tends to slowly converge to the equilibrium state especially in a large-scale network. Therefore, the whole simulation process can be quite computationally expensive even for mesoscopic simulations. This is a major concern for large-scale DUE applications. Some other methods in the literature are known to converge at an improved rate with a smaller number of iterations like WMSA (weighted-MSA), gradient-based methods, simple travel time responsive method, and alternating direction (AD) methods [24–26]. However, the level of convergence in user equilibrium widely depends on the complexity, size of the network, and the criterion for convergence [25] making it difficult to specify which method is better in general. The MSA algorithm used in our study is the modified version proposed by [27] which is considered to be computationally efficient [26]. The WMSA is a variant of MSA, with the only difference being in the step size calculation, proposed by [24]. Where the MSA algorithm iterates by 1/n of the demand, the implemented WMSA tends to move by 2/(n + 1) where n is the number of iterations. This modification in the step size sequence allows the solution in some cases to converge quicker. Analytical DUE methods such as [28] tend to perform faster on small networks but are known to be impractical for large-scale networks [29]. Another concern is the degree of validity of the assumed equilibrium-seeking behavior of drivers in a real-world network, as empirical evidence has suggested the existence of the noncooperative or selfish equilibrium at the macroscopic level only [30].

2.1.2. Noniterative Stochastic Route Choice (SRC)

The assumption underlying DUE that every driver has perfect knowledge of the path costs to help choose the best alternative does not necessarily hold in practice as there is a possibility of drivers making imperfect decisions due to inadequate knowledge [3]. Also, unknown factors in relation to preference of some routes that cannot be applied to the general behavior makes the intrinsic nature of traffic extremely variable. This perception error can be incorporated into the SRC model by decomposing the path utility into a deterministic and a random term. The probit model assuming a normally distributed random term and the logit model assuming a Gumbel distribution are two widely used methods to solve the SRC problem [7]. The probit model is comparatively computationally expensive due to the absence of a closed-form formula [31]. On the other hand, the logit model has a closed-form structure for easy calculation while having two main drawbacks, one being the issue with highly overlapped path set due to its inability of finding path correlation and the other being the scaling problem of not identifying the heterogeneity in drivers’ perception errors. To overcome the path overlapping problem, modified versions of the multinomial logit model have been introduced, such as the C-logit (which will be used in this paper), path size logit (PSL), and implicit availability/perception (IAP), by introducing a new utility function that tends to choose the spatially separated routes. Scaling problem can be addressed by introducing a scaling factor $\theta$ [32]. C-logit presents probability $P_k$ for selecting the paths using$$\begin{array}{}\text{(2)}& P_k=\frac{e^{\theta V_k-{\text{CF}}_k}}{{{\displaystyle \sum }}_{l\in K_i}{e}^{\theta V_l-{\text{CF}}_l}},\end{array}$$where $V_k$ and $V_l$ are the perceived utility for selected path k and alternative paths l, respectively and ${\text{CF}}_k$ denotes the “commonality factor” of the path k that describes the degree of overlapping for current path to all the alternative paths chosen. This ${\text{CF}}_k$ is calculated using the cost of links $L_{lk}$ common to paths k and l, $L_l$, $L_k$ the link costs of path l and k, respectively (expressed in hours), as follows$$\begin{array}{}\text{(3)}& {\text{CF}}_k=\beta \cdot \ln {\displaystyle \sum }_{l\in K_i}{\frac{L_{lk}}{L_l^{1/2}{L}_k^{1/2}}}^{\gamma },\end{array}$$Here, $\beta$ and $\gamma$ are the two opposite parameters where $\beta$ displays the importance of utility $V_k$ and $\gamma$ with lesser influence plays an opposite role to $\beta$ and is usually in the range of [0, 2] [23]. SRC model characteristics are applicable to both iterative (e.g., DSUE, DSUE-SRDTC [33–35]) and noniterative techniques. Noniterative SRC is used for all the analysis here as the focus of the paper is on the comparative analysis of iterative user-equilibrium versus noniterative route choice models. Going forward, noniterative SRC would be referred to as SRC.

2.1.3. Adaptive Driving

Compared with the DUE, SRC is much less computationally expensive because it does not require multiple iterations to reach the equilibrium state. However, without seeking equilibrium and with the introduction of perception errors, the quality of the SRC solution cannot be guaranteed and the network might evolve into unexpected gridlock. It has been shown that the network flow can be improved by considering adaptive driving for a proportion of drivers [8, 36]. Here, adaptive driving refers to the ability to alter the assigned path at each time step by choosing the new shortest path according to the prevailing traffic conditions. Our goal in this paper is to study how adaptive driving affects the estimated network flows and thus the quality of the SRC solution in comparison to the DUE solution as it may have significant practical implications if used for traffic operations optimization.

3. Model Description

In this paper, we perform a mesoscopic simulation-based DTA on a subnetwork of Melbourne metropolitan area (Figure 1) covering the city center using an existing calibrated and validated model in AIMSUN [10]. The network consists of 4375 links and 1977 nodes. Links have multiple attributes including the number of lanes, capacity, and the free-flow speed. The geometric configuration of the network is adapted from the Victoria Integrated Transport Model (VITM) including 492 traffic zones. The static OD matrix included 2,173,306 vehicles for a four-hour morning peak period (6 AM–10 AM). A time-dependent OD matrix is then estimated based on the initial static OD matrix and observed counts from Sydney Coordinated Adaptive Traffic System (SCATS) loop detector data from over 1500 actuated signals. For more details on calibration and validation of the DTA model, please see [10, 37–39]. The supply and demand models were calibrated and validated. For calibration and validation details, please see [10]. AIMSUN uses a simplified car following model as the supply model [26, 40–42]. The simplification results in a link fundamental diagram (FD) that estimates traffic flow characteristics at a mesoscopic level (see for [23] further details). The resulting triangular flow-density FD is as follows:$$\begin{array}{c}\text{(4)}& q_{t}k=\min V\cdot k_t,\frac{1-L{k}_t}{R}=\min V\cdot k_t,\frac{1}{R}1-\frac{k}{k_{\text{jam}}},R=R_{\text{experiment}}\cdot R_{\text{factor}},\end{array}$$where V and $k_t$ denote the speed limit and the density level of a given link. The first term in the right-hand side of the equation represents the undersaturated regime of the traffic while the second term represents the oversaturated regime where L and R are the effective length (sum of the vehicle length and the clearance distance) and the reaction time in a link, respectively. The link flow is zero either for k = 0 or for $k_{\text{jam}}$ = 1/L and is maximum at the intersection of two regimes. Reaction time, here, is the production of the global reaction time and the reaction time factor of the link. The local reaction time factor brings more flexibility in the calibration process.

[figure omitted; refer to PDF]

Real observed traffic count data obtained from loop detectors and SCATS detectors at the signalized intersections is used for comparative analysis and validation of simulated solutions. We use a relative gap of less than 5% to ensure the DUE solution is properly converged based on previous studies in the literature [27, 43–45]. The literature suggests that the relative gap has a high dependency on the size of the network. For a small size network, a very low relative gap close to zero is possible. However, for medium to large-scale networks, approaching a zero relative gap is very difficult. Hence, most practitioners working with large-scale networks tend to choose appropriate relative gap of around 5–10%. AIMSUN provides three different methods to calculate DUE, namely, MSA (method of successive averages), WMSA (weighted-MSA), and gradient-based. Each method has a different approach to converge to the appropriate relative gap (5%); hence their computation performance varies accordingly (see Figure 2). The MSA and WMSA processing times for each iteration are very similar as the underlying basic mechanism is the same with only the difference being in step size calculation, whereas the gradient-based method is more costly per iteration as compared to the other two methods. The second and important factor in the processing time is the number of iterations that is required to converge to the desired solution, in which WMSA performs the best with only 36 iterations as opposed to its competitors, MSA and gradient based with 120 and 121 iterations. Despite improved convergence of WMSA, DUE still took 3214.45 seconds which suggests that the noniterative SRC still outperforms DUE. All the computation metrics are recorded while using the same platform with specifications as follows: processor: Intel (R) Core (TM) i9-9900 CPU @ 3.10 GHz 3.10 GHz, RAM: 64.0 GB.

[figures omitted; refer to PDF]

For SRC, we use the C-logit model and consider different proportions of drivers with adaptive driving. C-logit is used with the default values of scale factor as 1 and commonality factor parameters, $\beta$ and $\gamma$, as 0.15 and 1, respectively. These parameters are fixed for all the OD pairs in the matrix throughout the simulation. The model takes into consideration 5 paths initially and, from those, chooses up to three most suitable paths by calculating the shortest path(s) from each OD pair using initial path travel times. At each simulation time interval, paths are recalculated using experienced average link travel times. Adaptive drivers also recalculate their shortest path at each time interval to check for the possible alternative routes for dynamic rerouting. The distribution of the number of paths used by vehicles in the DUE (see Figure 3) suggests considering three alternative paths in the path assignment of the SRC should be a reasonable assumption. To further strengthen our analysis, we have conducted analysis of running noniterative SRC with three, five, and seven paths with same the number of adaptive drivers and compared the results against real data (see Figure 4). The results using mean absolute error (MAE) suggest that the noniterative SRC achieves more accurate outcomes when three paths are used. When root mean square error (RMSE) and root mean square normalized (RMSN) errors are used, no significant difference is observed in the model outcomes when 3, 5, or 7 paths are used.

[figure omitted; refer to PDF]

The SRC with 40% to 80% adaptive drivers produces the closest result to the DUE, in terms of individual link flows, as indicated by MAE, RMSE, and RMSN. Interestingly, a number of links in the network experience flows larger than 1000 veh/hr in the DUE, which, however, exhibit a rather small flows in the SRC with very few adaptive drivers. This is because when the adaptive driving proportion is low, most drivers simply follow their initially assigned paths even if they encounter severe congestion enroute. With a higher proportion of adaptive drivers, flows are improved as also observed and reported in [8]. A comparison between the observed link volumes and the simulated volumes of both DUE and SRC is performed to further validate the results (see Figure 7). Results suggest that an increase in adaptive driving brings the solution closer to the observed data and produces relatively closer model outcomes to real data than DUE. SRC with 100 % adaptive driving estimates the closest link volumes to real-world observations. However, this does not necessarily conclude that SRC with high adaptive driving is qualitatively better than DUE as both models are characteristically different.

[figures omitted; refer to PDF]

4.2. Network-Level Comparison

Traffic congestion can be measured using a variety of methods. In this section, we use link speeds to measure the congestion level by comparing the observed speed with the maximum speed (i.e., the speed limit) [46, 47]. For each link in the network, we obtain the simulated speed $v_{i}t$ for each time interval $t$ and the maximum speed $v_i^{\text{max}}$. To quantify congestion as a function of speed, we define$$\begin{array}{}\text{(6)}& r_{i}t=\frac{v_{i}t}{v_i^{\text{max}}},\end{array}$$where $r_{i}t$ is the ratio of the link speed $v_{i}t$ over the maximum speed limit $v_i^{\text{max}}$. We then classify links into congested links with $c_i=1$ and uncongested links with $c_i=0$ using a threshold $\rho$ [47]:$$\begin{array}{}\text{(7)}& c_{i}t=\begin{array}{cc}1,& r_{i}t<\rho ,\\ 0,& r_{i}t\ge \rho ,\end{array}\end{array}$$where $\rho$ represents a user-defined congestion level. Here, we consider $\rho =0.2$ and $\rho =0.5$ (i.e., 20% and 50% of the maximum speed) as two congestion levels. With $\rho =0.2$, we are capturing links that are severely congested. On the other hand, with $\rho =0.5$ (Figure 8), any link with a speed less than half of the maximum speed is considered as congested. Consequently, $Ct$ is defined as the fraction of congested links in the network:$$\begin{array}{}\text{(8)}& Ct=\frac{{\displaystyle {\sum }_1^N{c}_{i}t}}{N},\end{array}$$where $N$ is the total number of links in the network and $Ct$ curve represents the overall congestion pattern in the network at any time t which evolves throughout the simulation.

[figures omitted; refer to PDF]

The main difference between the DUE and SRC with low percentage of adaptive drivers appears at the onset of congestion as well as in the recovery period. With 60% and 80% of adaptive drivers, the network manages to recover more rapidly similar to DUE as compared to the other adaptiveness levels. However, even with 80% of drivers being adaptive, the SRC model still produces slightly different congestion patterns in the network as compared to the DUE. Traffic congestion obtained from DUE solution tends to peak earlier as compared to SRC and recovers more rapidly as well.

Results also suggest that propagation and recovery rates of congested links, measured as the gradient of the curves shown in Figure 9, in the DUE model are smaller as compared to the SRC model with low percentage of adaptive drivers. Adaptive driving does not necessarily eliminate the propagation of congestion in the network, but it helps the network to recover from the gridlock state more rapidly than without adaptive drivers. Figure 10 shows and compares the spatial distribution of the congested links in the network under both the DUE and SRC scenarios. As expected, a higher dissimilarity is observed in the presence of low adaptiveness in SRC, and vice versa.

[figures omitted; refer to PDF]

[figures omitted; refer to PDF]

5. Application in Network Traffic Signal Optimization

Network traffic signal optimization is a well-studied problem in the literature [48–51]. Among the more commonly used methods to solve the problem are heuristic and metaheuristic approaches including genetic algorithms, neural networks, feedback control, and simulation-based optimization [50, 52–54]. We are specifically interested in formulating the problem as a bi-level program similar to what has been proposed in [49] where the upper-level problem is signal timing optimization and the lower-level problem is the DTA (equilibrium seeking DUE or nonequilibrium seeking SRC). In fact, the literature on network traffic signal optimization can be grouped into two, in which the first group of studies applies an equilibrium route choice model to solve the signal control problem [55–58] while the second group applies a nonequilibrium route choice model [59–65]. Table 1 provides a summary of the relevant studies in the literature classified based on the network size and utilized traffic assignment model. The aim here is not to discuss or advocate whether, for network traffic signal optimization, the use of an equilibrium route choice is a more reasonable and realistic approach or a nonequilibrium approach. Rather, we only aim to provide a comparative analysis to shed light on how use of each of the models affects the solution.

Table 1

A summary of literature on simulation-based traffic signal optimization studies.

Traffic assignment	Network size	Summary
Deterministic user equilibrium (DUE)	Local (<5 signals or <5 available routes)	Very few studies are found carrying out signal optimization using DUE in relatively small networks with few routing options including [66] that used DUE to calculate the reserve capacity of the network. Meanwhile, [67, 68] implemented DUE on a local network with isolated signal or few signals as well as with a real network of multiple signalized intersections.
	Global (>5 signals or >5 available routes)	DUE is the popular choice of traffic assignment for signal optimization in relatively large networks [67–70]. These studies yield efficient techniques for signal optimization based on DUE traffic assignment with some providing a comparison between DUE with SO, concluding significant advantage of DUE for adaptive signals.
Stochastic user equilibrium (SUE)	Local (<5 signals or <5 available routes)	SUE is more prominent in localized signal optimization [57, 68, 71, 72]. The study reported in [58] has provided a comparison of the performance of deterministic and stochastic user equilibrium techniques over different networks and observed improved results for SUE in local networks.
	Global (>5 signals or >5 available routes)	Very few studies are found in the literature using SUE in relatively large networks including [72] that has recommended the use of DUE on global traffic signal optimization.
System optimal (SO)	Local/global	References [67, 73, 74] used SO in both large and small networks. Reference [67] argued that DUE is advantageous over SO in designing adaptive signals in a network.
Stochastic route choice (SRC)	Local (<5 signals or <5 available routes)	Use of SRC mostly prevails in the studies on a single freeway or highway with few alternate routes available [65, 75–79].
	Global (>5 signals or >5 available routes)	A few studies are found in the literature that used relatively large networks and SRC [80–82]. In [82], the selection of route choice model was arbitrary as the main focus of the study on the optimization part of the problem rather than the simulation.

Here, we formulate a simple network traffic signal optimization problem in which we consider the signal cycle length as the only decision variable which is assumed to be homogeneous across the selected traffic signals in the study network. To reduce the computational burden of the optimization in this case study, we select 19 traffic signals in the CBD area only to be included in the optimization problem (see Figure 11) while the remainder of the traffic signals in the network are actuated with a pre-determined and fixed cycle length. Note that the simulations are run on the complete network of Melbourne as shown earlier in Figure 1 while the optimization and traffic measurements are only conducted for the CBD area. The goal of this exercise is to analyze the impact of DUE or SRC selection on the developed congestion patterns and the estimated optimal solution.$$\begin{array}{c}\text{(9)}& \underset{\omega }{\min }\frac{{{\displaystyle \sum }}_{t=1}^{\tau }ct}{N},\text{subject\hspace{0.17em}to}\\ ct=\text{DTA}\omega ,\end{array}$$where $\omega$ represents signal cycle length, $\tau$ represents the study time period, and $N$ is number of time intervals between 1 and $\tau$. The upper-level problem (equation (6)) minimizes the average number of congested links in the network over the study time period, while the lower-level problem (equation (7)) represents the simulation-based DTA (SRC or DUE model) with variable cycle length $\omega$ for the selected traffic signals.

[figure omitted; refer to PDF]

To solve the signal optimization problem, involving a simulation-based DTA model, a variety of simulation-based optimization (SBO) methods could be applied. Without loss of generality, SBO refers to the optimization of an objective function subject to certain constraints, both of which are evaluated through stochastic computer simulations. In transportation network design and analysis, SBO is frequently applied in the presence of “black-box” simulation [83–86]. For an overview of different SBO methods, refer to [87]. In this paper, we adopted a naïve response surface method where a parametric model mimicking the simulation input-output mapping was constructed based on results of selected sample points. Specifically, we randomly selected 7 sample points representing different cycle lengths between 60 and 150 sec, for each of which simulations were performed and the objective function was evaluated. Note that SRC was run for five replications with 30% adaptive drivers. Admittedly such a naïve approach lacks a mathematical rigor without accounting for any infill sample points as one would see in the literature [88], but it suffices to highlight the difference between the two modeling approaches as the key emphasis in this paper.

Figure 12 summarizes the analysis results for both DUE and SRC models. The formulated objective function value, namely, the average fraction of congested links in the study subnetwork, varies for different cycle length values. However, the variations for both SRC and DUE models show a similar pattern. When SRC is used, the optimal cycle length is estimated to be 126 seconds while for the DUE the optimal cycle length is 93 seconds. The observed and quantified difference demonstrates how selection of the route choice model affects the optimization solution. Therefore, considerable care is needed when selecting the route choice model and comparing the results of DUE and SRC in traffic planning and operations applications.

[figures omitted; refer to PDF]

6. Conclusions

Finding dynamic equilibrium is computationally expensive in large-scale networks, but it often provides robust and consistent model outcomes given the equilibrium-seeking behavior. The SRC model introduces a random perception error in the path assignment without seeking the equilibrium state of the network, thereby resulting in a more probabilistic DTA model outcomes. Given the increasing use of SRC model in practice and its lower computational costs, in this paper we aimed to provide a quantitative comparison between the two modeling approaches when applied to a large-scale network and used in traffic operations optimization. Due to the absence of the equilibrium-seeking behavior, the SRC with low adaptive driver’s population tends to produce more congestion in the network and, hence, lower link flows. However, by introducing adaptive driving, the quality of the SRC solution in terms of overall congestion rate is shown to be improved as measured by the increased link flows as well as the faster recovery from congestion. With increasing the proportion of adaptive drivers in the network, the SRC model performs more closely to both DUE and real-world observations in terms of link flows and congestion patterns. However, if the whole population of drivers behave adaptively, the estimated link flows drift away from DUE. These results have practical implications for the application of large-scale DTA models for planning and operations purposes. A simplified case study on network traffic signal optimization was performed and revealed the practical implications of use of each of the route choice models. While DUE is widely used in academic literature in many applications including network design problems, the increasing use of SRC in practice and ongoing discussions among practitioners on the necessity of equilibrium (or lack thereof) in large-scale modeling studies have motivated us to provide evidence on the significant practical differences that each modeling approach provides. We hope the provided discussions and quantitative comparisons in this paper help practitioners with selection of the most appropriate route choice model for their specific purposes.

References

[1] L. Zhang, "Behavioral foundation of route choice and traffic assignment," Transportation Research Record: Journal of the Transportation Research Board, vol. 2254 no. 1,DOI: 10.3141/2254-01, 2012.

[2] P. H. L. Bovy, E. Stern, "Route choice: wayfinding in transport networks," Studies in Operational Regional Science, vol. 9, 1990.

[3] J. N. Prashker, S. Bekhor, "Route choice models used in the stochastic user equilibrium problem: a review," Transport Reviews, vol. 24 no. 4, pp. 437-463, DOI: 10.1080/0144164042000181707, 2004.

[4] J. G. Wardrop, "Road paper. some theoretical aspects of road traffic research," Proceedings of the Institution of Civil Engineers, vol. 1 no. 3, pp. 325-362, 1952.

[5] C. F. Daganzo, Y. Sheffi, "On stochastic models of traffic assignment," Transportation Science, vol. 11 no. 3, pp. 253-274, DOI: 10.1287/trsc.11.3.253, 2008.

[6] C. Gawron, "An iterative algorithm to determine the dynamic user equilibrium in a traffic simulation model," International Journal of Modern Physics C, vol. 9 no. 3, pp. 393-407, 1998.

[7] X. Lai, J. Li, "Modelling stochastic route choice behaviours with a closed-form mixed logit model," Mathematical Problems in Engineering, vol. 2015 no. 2,DOI: 10.1155/2015/729089, 2015.

[8] M. Saberi, H. S. Mahmassani, A. Zockaie, "Network Capacity, Traffic instability, and adaptive driving: findings from simulated urban network experiments," EURO Journal on Transportation and Logistics, vol. 3 no. 3-4, pp. 289-308, DOI: 10.1007/s13676-013-0040-2, 2014.

[9] A. Mahmoodjanlou, M. L. Hazelton, K. Parry, "Apples versus oranges? comparing deterministic and stochastic day-to-day traffic assignment models," Transportmetrica B, vol. 7 no. 1, pp. 1426-1443, DOI: 10.1080/21680566.2019.1627953, 2019.

[10] S. Shafiei, Z. Gu, M. Saberi, "Calibration and validation of a simulation-based dynamic traffic assignment model for a large-scale congested network," Simulation Modelling Practice and Theory, vol. 86 no. 2017, pp. 169-186, DOI: 10.1016/j.simpat.2018.04.006, 2018.

[11] S. C. Dafermos, F. T. Sparrow, "The traffic assignment problem for a general network," Journal of Research of the National Bureau of Standards B, vol. 73 no. 2, pp. 91-118, 1969.

[12] M. Ben-Akiva, A. De Palma, K. Isam, "Dynamic network models and driver information systems," Transportation Research Part A: General, vol. 25 no. 5, pp. 251-266, 1991.

[13] Y.-C. Chiu, J. Bottom, M. Mahut, "Dynamic traffic assignment: a primer," 2011. Transportation Research Circular, No. E-C153

[14] E. Frejinger, Route Choice Analysis: Data, Models, Algorithms and Applications, 2008.

[15] H. S. Mahmassani, G.-L. Chang, R. Herman, "Individual decisions and collective effects in a simulated traffic system," Transportation Science, vol. 20 no. 4, pp. 258-271, 1986.

[16] H. S. Mahmassani, G.-L. Chang, "On boundedly rational user equilibrium in transportation systems," Transportation Science, vol. 21 no. 2, pp. 89-99, 1987.

[17] Y. Lou, Y. Yin, S. Lawphongpanich, "Robust congestion pricing under boundedly rational user equilibrium," Transportation Research Part B: Methodological, vol. 44 no. 1, pp. 15-28, 2010.

[18] M. Miralinaghi, Y. Lou, Y. Hsu, R. Shabanpour, Y. Shafahi, "Multiclass fuzzy user equilibrium with endogenous membership functions and risk‐taking behaviors," Journal of Advanced Transportation, vol. 50 no. 8, pp. 1716-1734, 2016.

[19] H. Xu, Y. Lou, Y. Yin, J. Zhou, "A prospect-based user equilibrium model with endogenous reference points and its application in congestion pricing," Transportation Research Part B: Methodological, vol. 45 no. 2, pp. 311-328, 2011.

[20] L. Zhang, "Agent-based behavioral model of spatial learning and route choice," Proceedings of the Transportation Research Board 85th Annual Meeting, .

[21] J. B. Lesort, A. Nuzzolo, F. Russo, A. Vitetta, "A modified logit route choice model overcoming path overlapping problems. specification and some calibration results for interurban networks," Transportation and Traffic Theory, pp. 697-711, .

[22] J. Almond, "Traffic assignment with flow-dependent journey times," .

[23] Aimsun, Aimsun Next 20 User’s Manual, Aimsun Next Version 20.0.3, . Barcelona, Spain, 2021, http://aimsun.com.aimsun.20.0/doc/UsersManual/Intro.html

[24] H. X. Liu, X. He, B. He, "Method of successive weighted averages (MSWA) and self-regulated averaging schemes for solving stochastic user equilibrium problem," Networks and Spatial Economics, vol. 9 no. 4, pp. 485-503, DOI: 10.1007/s11067-007-9023-x, 2007.

[25] M. Carey, Y. E. Ge, "Comparison of methods for path flow reassignment for dynamic user equilibrium," Networks and Spatial Economics, vol. 12 no. 3, pp. 337-376, DOI: 10.1007/s11067-011-9159-6, 2012.

[26] J. Barceló, Fundamentals of Traffic Simulation, 2010.

[27] M. Florian, M. Mahut, N. Tremblay, "Application of a simulation-based dynamic traffic assignment model," European Journal of Operational Research, vol. 189 no. 3, pp. 1381-1392, DOI: 10.1016/j.ejor.2006.07.054, 2008.

[28] X. Ban, H. X. Liu, "A link-node complementarity model and solution algorithm for dynamic user equilibria with exact flow propagations," Transportation Research Part B: Methodological, vol. 42 no. 9, pp. 823-842, DOI: 10.1016/j.trb.2008.01.006, 2008.

[29] B. N. Janson, "Dynamic traffic assignment for urban road networks," Transportation Research Part B, vol. 25 no. 2–3, pp. 143-161, DOI: 10.1016/0191-2615(91)90020-J, 1991.

[30] M. Yildirimoglu, O. Kahraman, "Searching for empirical evidence on traffic equilibrium," PLoS One, vol. 13 no. 5,DOI: 10.1371/journal.pone.0196997, 2018.

[31] T. Yai, S. Iwakura, S. Morichi, "Multinomial probit with structured covariance for route choice behavior," Transportation Research Part B: Methodological, vol. 31 no. 3, pp. 195-207, 1997.

[32] S. Pravinvongvuth, A. Chen, "Adaptation of the paired combinatorial logit model to the route choice problem," Transportmetrica, vol. 1 no. 3, pp. 223-240, 2005.

[33] P. C. A. Vythoulkas, "Dynamic stochastic assignment model for the analysis of general networks," Transportation Research Part B: Methodological, vol. 24 no. 6, pp. 453-469, 1990.

[34] M. Du, H. Tan, A. Chen, "A faster path-based algorithm with Barzilai-Borwein step size for solving stochastic traffic equilibrium models," European Journal of Operational Research, vol. 290 no. 3, pp. 982-999, 2021.

[35] J. Long, W. Y. Szeto, Q. Shi, Z. Gao, H.-J. Huang, "A nonlinear equation system approach to the dynamic stochastic user equilibrium simultaneous route and departure time choice problem," Transportmetrica A: Transport Science, vol. 11 no. 5, pp. 388-419, 2015.

[36] C. F. Daganzo, V. V. Gayah, E. J. Gonzales, "Macroscopic relations of urban traffic variables: bifurcations, multivaluedness and instability," Transportation Research Part B: Methodological, vol. 45 no. 1, pp. 278-288, DOI: 10.1016/j.trb.2010.06.006, 2011.

[37] S. Shafiei, M. Saberi, A. Zockaie, M. Sarvi, "Sensitivity-based linear approximation method to estimate time-dependent origin–destination demand in congested networks," Transportation Research Record: Journal of the Transportation Research Board, vol. 2669 no. 1, pp. 72-79, DOI: 10.3141/2669-08, 2017.

[38] S. Shafiei, M. Saberi, H. L. Vu, "Integration of departure time choice modeling and dynamic origin–destination demand estimation in a large-scale network," Transportation Research Record: Journal of the Transportation Research Board, vol. 2674 no. 9, pp. 972-981, DOI: 10.1177/0361198120933267, 2020.

[39] S. Shafiei, M. Saberi, H. L. Vu, "Nonlinearity in time-dependent origin-destination demand estimation in congested networks," Proceedings of the 2019 IEEE Intelligent Transportation Systems Conference (ITSC), .

[40] M. Mahut, "Speed-maximizing car-following models based on safe stopping rules," Proceedings of the 78th Annual Meeting of the Transportation Research Board, .

[41] M. Mahut, "Behavioral car following models," 2000. Ph. D. thesis

[42] M. A. Mahut, "Discrete flow model for dynamic network loading," 2001. Ph. D. dissertation

[43] J. Auld, O. Verbas, M. Stinson, "Agent-based dynamic traffic assignment with information mixing," Procedia Computer Science, vol. 151 no. 2018, pp. 864-869, DOI: 10.1016/j.procs.2019.04.119, 2019.

[44] G. Gentile, "Solving a dynamic user equilibrium model based on splitting rates with gradient projection algorithms," Transportation Research Part B: Methodological, vol. 92, pp. 120-147, DOI: 10.1016/j.trb.2016.02.005, 2016.

[45] H. Sbayti, C. C. Lu, H. S. Mahmassani, "Efficient implementation of method of successive averages in simulation-based dynamic traffic assignment models for large-scale network applications," Transportation Research Record: Journal of the Transportation Research Board, vol. 2029 no. 1, pp. 22-30, DOI: 10.3141/2029-03, 2007.

[46] M. Saberi, H. Hamedmoghadam, M. Ashfaq, "A simple contagion process describes spreading of traffic jams in urban networks," Nature Communications, vol. 11 no. 1,DOI: 10.1038/s41467-020-15353-2, 2020.

[47] D. Li, B. Fu, Y. Wang, "Percolation transition in dynamical traffic network with evolving critical bottlenecks," Proceedings of the National Academy of Sciences, vol. 112 no. 3, pp. 669-672, 2015.

[48] D. I. Robertson, "Transyt a traffic network study tool," RRL Report LR 253, Road Resarch Laboratory Crowthorne, 1969.

[49] D. Sun, R. F. Benekohal, S. T. Waller, "Bi-level programming formulation and heuristic solution approach for dynamic traffic signal optimization," Computer-Aided Civil and Infrastructure Engineering, vol. 21 no. 5, pp. 321-333, DOI: 10.1111/j.1467-8667.2006.00439.x, 2006.

[50] R. Mohebifard, A. Hajbabaie, "Optimal network-level traffic signal control: a benders decomposition-based solution algorithm," Transportation Research Part B: Methodological, vol. 121, pp. 252-274, DOI: 10.1016/j.trb.2019.01.012, 2019.

[51] K. Jang, H. Kim, I. G. Jang, "Traffic signal optimization for oversaturated urban networks: queue growth equalization," IEEE Transactions on Intelligent Transportation Systems, vol. 16 no. 4, pp. 2121-2128, DOI: 10.1109/TITS.2015.2398896, 2015.

[52] M. Zlatkovic, X. Zhou, "Integration of signal timing estimation model and dynamic traffic assignment in feedback loops: system design and case study," Journal of Advanced Transportation, vol. 49 no. 6, pp. 683-699, DOI: 10.1002/atr.1295, 2015.

[53] B. Park, C. J. Messer, T. Urbanik, "Traffic signal optimization program for oversaturated conditions: genetic algorithm approach," Transportation Research Record: Journal of the Transportation Research Board, vol. 1683 no. 1, pp. 133-142, DOI: 10.3141/1683-17, 1999.

[54] D. Srinivasan, M. C. Choy, R. L. Cheu, "Neural networks for real-time traffic signal control," IEEE Transactions on Intelligent Transportation Systems, vol. 7 no. 3, pp. 261-272, DOI: 10.1109/TITS.2006.874716, 2006.

[55] N. Gartner, "Area traffic control and network equilibrium," Traffic Equilibrium Methods, 1976.

[56] N. H. Gartner, "OPAC: a demand-responsive strategy for traffic signal control," IFAC Proceedings Volumes, vol. 23 no. 2, pp. 241-244, DOI: 10.1016/S1474-6670(17)52677-4, 1990.

[57] Y. Sheffi, W. B. Powell, "Optimal signal settings over transportation networks," Journal of Transportation Engineering, vol. 109 no. 6, pp. 824-839, DOI: 10.1061/(ASCE)0733-947X(1983)109:6(824), 1983.

[58] E. Cascetta, M. Gallo, B. Montella, "Models and algorithms for the optimization of signal settings on urban networks with stochastic assignment models," Annals of Operations Research, vol. 144 no. 1, pp. 301-328, DOI: 10.1007/s10479-006-0008-9, 2006.

[59] T.-Y. Hu, H. S. Mahmassani, "Day-to-Day evolution of network flows under real-time information and reactive signal control," Transportation Research Part C: Emerging Technologies, vol. 5 no. 1, pp. 51-69, 1997.

[60] B. G. Heydecker, C. Cai, C. K. Wong, "Adaptive dynamic control for road traffic signals," Proceedings of the IEEE International Conference on Networking, Sensing and Control, ICNSC 2007, .

[61] B. G. Heydecker, "Objectives, stimulus and feedback in signal control of road traffic," Journal of Intelligent Transportation Systems, vol. 8 no. 2, pp. 63-76, 2004.

[62] M. Smith, R. Mounce, "A splitting rate model of traffic Re-routeing and traffic control," Procedia-Social and Behavioral Sciences, vol. 17, pp. 316-340, 2011.

[63] T. Toledo, T. Balasha, M. Keblawi, "Optimization of actuated traffic signal plans using a mesoscopic traffic simulation," Journal of Transportation Engineering Part A: Systems, vol. 146 no. 6,DOI: 10.1061/JTEPBS.0000363, 2020.

[64] J. Stevanovic, A. Stevanovic, P. T. Martin, T. Bauer, "Stochastic optimization of traffic control and transit priority settings in VISSIM," Transportation Research Part C: Emerging Technologies, vol. 16 no. 3, pp. 332-349, DOI: 10.1016/j.trc.2008.01.002, 2008.

[65] A. Stevanovic, P. T. Martin, J. Stevanovic, "VisSim-based genetic algorithm optimization of signal timings," Transportation Research Record: Journal of the Transportation Research Board, vol. 2035 no. 1, pp. 59-68, DOI: 10.3141/2035-07, 2007.

[66] G. Ziyou, S. Yifan, "A reserve capacity model of optimal signal control with user-equilibrium route choice," Transportation Research Part B: Methodological, vol. 36 no. 4, pp. 313-323, DOI: 10.1016/S0191-2615(01)00005-4, 2002.

[67] H. Yu, R. Ma, H. M. Zhang, "Optimal traffic signal control under dynamic user equilibrium and link constraints in a general network," Transportation Research Part B: Methodological, vol. 110, pp. 302-325, DOI: 10.1016/j.trb.2018.02.009, 2018.

[68] A. Karoonsoontawong, S. T. Waller, "Application of reactive Tabu search for combined dynamic user equilibrium and traffic signal optimization problem," Transportation Research Record: Journal of the Transportation Research Board, vol. 2090 no. 1, pp. 29-41, DOI: 10.3141/2090-04, 2009.

[69] F. Teklu, A. Sumalee, D. Watling, "A genetic algorithm approach for optimizing traffic control signals considering routing," Computer-Aided Civil and Infrastructure Engineering, vol. 22 no. 1, pp. 31-43, DOI: 10.1111/j.1467-8667.2006.00468.x, 2007.

[70] J. Guo, Y. Kong, Z. Li, W. Huang, J. Cao, Y. Wei, "A model and genetic algorithm for area-wide intersection signal optimization under user equilibrium traffic," Mathematics and Computers in Simulation, vol. 155, pp. 92-104, DOI: 10.1016/j.matcom.2017.12.003, 2019.

[71] N. Nihan, M. Hamed, G. Davis, "Interactions between driver information, route choice, and optimal signal timing on a simple network," Journal of Advanced Transportation, vol. 29 no. 2, pp. 163-182, DOI: 10.1002/atr.5670290203, 1995.

[72] Z. Yang, S. Lu, X. Liu, "Combined traffic signal control and route guidance: multiple user class traffic assignment model versus discrete choice model," Proceedings of the Multiconference on “Computational Engineering in Systems Applications”, .

[73] H. R. Varia, S. L. Dhingra, "Dynamic optimal traffic assignment and signal time optimization using genetic algorithms," Computer-Aided Civil and Infrastructure Engineering, vol. 19 no. 4, pp. 260-273, DOI: 10.1111/j.1467-8667.2004.00353.x, 2004.

[74] A. S. Abdelfatah, H. S. Mahmassani, "A simulation-based signal optimization algorithm within a dynamic traffic assignment framework," j, vol. 1 8, pp. 428-433, DOI: 10.1109/itsc.2001.948695, 2001.

[75] A. Z. Stevanovic, P. T. Martin, "Assessment of the suitability of microsimulation as a tool for the evaluation of macroscopically optimized traffic signal timings," Journal of Transportation Engineering, vol. 134 no. 2, pp. 59-67, DOI: 10.1061/(ASCE)0733-947X(2008)134:2(59), 2008.

[76] A. Stevanovic, J. Stevanovic, P. T. Martin, "Optimizing signal timings from the field: VISGAOST and VISSIM?ASC/3 software-in-the-loop simulation," Transportation Research Record: Journal of the Transportation Research Board, vol. 2128 no. 1, pp. 114-120, DOI: 10.3141/2128-12, 2009.

[77] J. Mulandi, A. Stevanovic, P. T. Martin, "Cross-evaluation of signal timing optimized by various traffic simulation and signal optimization tools," Transportation Research Record: Journal of the Transportation Research Board, vol. 2192 no. 1, pp. 147-155, DOI: 10.3141/2192-14, 2010.

[78] A. Stevanovic, C. Kergaye, J. Stevanovic, "Evaluating robustness of signal timings for varying traffic flows," Transportation Research Record: Journal of the Transportation Research Board, vol. 2259 no. 1, pp. 141-150, DOI: 10.3141/2259-13, 2011.

[79] B. B. Park, J. D. Schneeberger, "Evaluation of traffic signal timing optimization methods using a stochastic and microscopic simulation program," Research Report UVACTS-5-0-4, University of Virginia, 2003. Charlottesville, VA, USA

[80] L. Zheng, X. Xue, C. Xu, B. Ran, "A stochastic simulation-based optimization method for equitable and efficient network-wide signal timing under uncertainties," Transportation Research Part B: Methodological, vol. 122, pp. 287-308, DOI: 10.1016/j.trb.2019.03.001, 2019.

[81] L. Zheng, C. Xu, P. J. Jin, B. Ran, "Network-wide signal timing stochastic simulation optimization with environmental concerns," Applied Soft Computing Journal, vol. 77, pp. 678-687, DOI: 10.1016/j.asoc.2019.01.046, 2019.

[82] N. R. Sabar, L. M. Kieu, E. Chung, T. Tsubota, P. E. Maciel de Almeida, "A memetic algorithm for real world multi-intersection traffic signal optimisation problems," Engineering Applications of Artificial Intelligence, vol. 63, pp. 45-53, DOI: 10.1016/j.engappai.2017.04.021, 2017.

[83] L. Zheng, X. Xue, C. Xu, B. Ran, "A stochastic simulation-based optimization method for equitable and efficient network-wide signal timing under uncertainties," Transportation Research Part B: Methodological, vol. 122, pp. 287-308, 2019.

[84] C. Osorio, V. Punzo, "Efficient calibration of microscopic car-following models for large-scale stochastic network simulators," Transportation Research Part B: Methodological, vol. 119, pp. 156-173, 2019.

[85] X. M. Chen, C. Xiong, X. He, Z. Zhu, L. Zhang, "Time-of-day vehicle mileage fees for congestion mitigation and revenue generation: a simulation-based optimization method and its real-world application," Transportation Research Part C: Emerging Technologies, vol. 63, pp. 71-95, 2016.

[86] Z. Gu, S. T. Waller, M. Saberi, "Surrogate‐based toll optimization in a large‐scale heterogeneously congested network," Computer‐Aided Civil and Infrastructure Engineering, vol. 34 no. 8, pp. 638-653, 2019.

[87] S. Amaran, N. V Sahinidis, B. Sharda, S. J. Bury, "Simulation optimization: a review of algorithms and applications," Annals of Operations Research, vol. 240 no. 1, pp. 351-380, 2016.

[88] X. Chen, L. Zhang, X. He, C. Xiong, Z. Li, "Surrogate‐based optimization of expensive‐to‐evaluate objective for optimal highway toll charges in transportation network," Computer‐Aided Civil and Infrastructure Engineering, vol. 29 no. 5, pp. 359-381, 2014.