# Travel Time Reliability for Urban Networks: Modelling and Empirics.

1. Introduction

Urban travel times are well recognized to be variable due to time-varying demand, stochastic capacity, different control measures (e.g., signal control at intersections), and external factors (e.g., weather conditions). The most widely used travel time uncertainty modelling philosophy is based on the concept of probability. With the probability distributions, the reliability analysis of link/route/network travel times can be conducted. In the literature, travel time distribution estimation methods in the urban context can be categorized into three groups: analytical approach, statistical (or numerical) approximation, and simulation.

For the analytical approach, Zheng et al. [1, 2] proposed an analytical link and trip travel time distribution model for urban signalized arterials. The model explicitly accounts for the influence of stochastic traffic processes, variation of traffic demand, queue dynamics, and different traffic signal control schemes (including cycle time, green time, and offset) on travel time variability. The "spillback" phenomenon (queued traffic from one intersection affecting operations at adjacent intersections) is also taken into account. The advantage of the analytical model lies in its transferability and applicability in different traffic conditions. However, the main drawback is the complexity of the modelling process for network-wide travel time distribution estimation.

For the statistical (or numerical) methods, normal, lognormal, truncated normal/lognormal, Gamma, and Weibull distributions [3-5] are widely applied to describe travel time reliability (variability). Due to the complexity of urban traffic conditions, a single distribution model could not well represent the travel time distribution of an urban road. A mixture of statistical distribution models [6, 7] has been proposed to capture multistate traffic conditions and the influence of signal control at intersections for urban roads. Clark and Watling  employed the method of moments to calculate various travel time reliability measures (expected travel time, variance, skewness, and kurtosis). The values of these four measures were then fitted to a family of probability densities called "Johnson" curves. In order to capture both the vehicle-to-vehicle and day-to-day variability in travel time data, Kim and Mahmassani  proposed a Gamma-Gamma mixture distribution model. Based on their proposed model, the heterogeneity within these two variability types under different weather conditions can be described as well.

For the simulation method, Kim et al.  employed traffic simulation models to derive travel time distributions under different scenarios considering various demand and supply uncertainty factors, such as weather, traffic incidents, work zones, and traffic control. Chen and Zhou  proposed a simulation approach to derive network travel times by solving a traffic equilibrium problem considering demand uncertainty. The network travel time samples were further used to estimate the network travel time distribution and its statistics. The main advantage of the simulation approach lies in its flexibility. However, intensive computation is required to conduct the simulation.

Although the concept of travel time reliability is relatively new in the field of transportation compared to some other engineering fields, various travel time reliability measures have been proposed by researchers, for example, standard deviation, statistical range , buffer index , 90th or 95th percentile , skewness-width , travel time budget or total travel time budget , and mean-excess total travel time . The choice of the most suitable reliability model depends on the context and purpose for which reliability is evaluated and on the characteristic of the travel time distribution . Skewness is sometimes a relevant characteristic because the more extreme long travel times may deter travelers to choose a certain route and the economic value of extreme long travel times is important . In cases for which skewness in the travel time probability distribution is not an issue, the variance of the travel time is a logical reliability measure. Total travel time budget (TTTB) and mean-excess total travel time (METTT) are two risk measures proposed for assessing transportation network performance under uncertainty [11, 17]. TTTB is defined as the minimum travel time threshold that satisfies a certain reliability requirement given by decision-makers at certain confidence level. Xu et al.  define the METT as the conditional expectation of travel times exceeding the corresponding TTTB at a given confidence level. The numerical example shows that applying METTT in the network design problem under uncertainty lowers the construction cost and reduces the unacceptable risk in the distribution tail substantially while incurring a minor TTTB increase.

Up to now, most research focuses on modelling or measuring link/segment level or route level travel time reliability. Research on network wide travel time reliability is rather limited and mainly based on simulation. Mahmassani et al.  used both simulated data and field GPS trajectory data to model travel time reliability. Their investigation with field data is rather limited. Moreover, GPS data is only a sample of the total traffic on the network and no flow data is available in their study area. The question is to what extent the derived relationship between the travel time variability and the network flow or density deduced from GPS sample data can represent the real situation.

The network travel time reliability has been considered as an important aspect in network management, control, and design . In [1, 2], the authors have developed an analytical travel time distribution model for link and route levels. When it comes to the network level, the previously developed "pure analytical model" is unable to describe the travel time variability. Therefore, in this paper, a network travel time distribution model based on flexible probability densities known as Johnson curves is introduced in Section 2. Different from previous research on travel time reliability modelling with Johnson curves estimated based on the method of moments, we apply a percentile-point method to estimate Johnson curves. The capabilities of Johnson curves are demonstrated with field travel time data. In Section 3, we discuss the calculation of travel time reliability measures for urban networks and their connection with network traffic characteristics (e.g., network traffic density and flow). An illustration of travel time reliability based on empirical data for a real urban network of a city in China using ANPR data is provided as well. Finally, some discussion and conclusions are provided in the last section.

2. Network Travel Time Distribution

2.1. Vehicle-to-Vehicle Travel Time Variability from Empirical Observations. Travel times that vehicles experience in the urban network can be influenced by many factors, such as fluctuations in traffic demand and supply, signal control at intersections, turning vehicles from cross streets, bus maneuvers at bus stops, parking vehicles along the roadside, crossing pedestrians and cyclists, and, of course, intermediate stops of the vehicles themselves, for example, for certain activities of the driver along the trip. As a result, urban travel times are rather variable given known traffic condition (traffic flow). Figure 1(a) illustrates the empirical travel time-flow relationship derived from local 5 min aggregated measurements of time-mean flow and median travel time for each 5 min on an urban arterial road "Renmin Road" in April 2015 in Changsha City, China. The region "A" illustrates the low demand condition and region "B" shows the high demand condition. It can be clearly seen in Figure 1(a) that there is no one-one correspondence relationship between travel time and flow over the whole range of traffic flow. Even in the condition of low traffic volume (region "A"), for a certain traffic flow, there is a big range of travel time corresponding to it. In the high demand condition, a large variation of travel times can be observed. Figure 1(b) illustrates travel times experienced by vehicles departing within a short time period (e.g., 1 min). The difference between the maximum travel time and minimum travel time can be as large as 221 seconds even within a small departure time period of 1 min.

2.2. Modelling Network Travel Time Distribution Using Johnson Curves. Due to a lot of uncertainties in urban transportation systems, travel times on both the link level and the route level are rather variable, which aggregately results in the variability of travel times for the network level. As shown in some empirical studies, the travel time distribution is unlikely to have a normal distribution but is heavily skewed with a long fat tail [15, 21, 22]. In this case, flexible probability densities known as Johnson curves  can be applied to model travel time data. The Johnson curves have been widely applied in the field of agriculture and environmental engineering and proved to be a powerful tool to model nonnormal and multi-modality data [24, 25]. One important property of Johnson curves is that parameters in the model can be estimated analytically or explicitly, which facilitate the generalization to different situations and data sets. In the following, the Johnson curves approach is discussed.

Let X denote the variable of intersection-normalized travel time which is defined as the travel time trip divided by the number of intersections on the trip. Jones et al.  and Mahmassani et al.  suggested that the use of distance-normalized travel time can help exclude the source of variability coming from trip distance and only captures variability due to speed variation. This is the case with trips without intersections, when the variance of travel time depends only on the variation in speed and congestion at bottlenecks. However, on urban arterials, a major part of travel time variance comes from signalized intersections. Therefore, we propose to normalize travel time on the number of intersections for a certain trip and focus on investigating the intersection-normalized travel time distribution for urban networks. This normalized travel time will be denoted by "travel time rate."

In the Johnson system, three types of probability density functions are used:

(1) The unbounded system (or [S.sub.U]): z = [gamma] + [eta] arcsinh((x - [epsilon])/[lambda])

(2) The bounded system (or [S.sub.B]): z = [gamma] + [eta] ln((x - [epsilon])/([lambda] + [epsilon] - x))

(3) The lognormal system (or [S.sub.L]): z = [gamma] + [eta], [epsilon] ln ((x - [epsilon])/[lambda])

where z is a standard normal variable and [gamma], [eta], [epsilon], and [lambda] are parameters of the Johnson curves. In the [S.sub.U] curves, the domain is the entire real line. By contrast, the domain of [S.sub.B] curves is defined over the open interval between [epsilon] and [epsilon] + [lambda]. [S.sub.L] curves can be further simplified by eliminating the parameter [lambda] as z = [gamma] + [eta]ln(x - [epsilon]). In the literature, Johnson curves were usually estimated by the method of moments [8, 27, 28], which requires calculating four moments (mean, standard deviation, skewness, and kurtosis) of a distribution and sometimes it is not straightforward to derive these moments. Slifker and Shapiro  proposed a method to estimate parameters of Johnson curves based on percentiles. Their method has shown to be an easier and more reliable procedure compared with methods based on moments . The procedure to estimate parameters is as follows.

Step 1. Estimate the percentiles in the travel time data corresponding to the cumulative probabilities [P.sub.-3z], [P.sub.-z], [P.sub.z], and [P.sub.3z] associated with four numbers -3z, -z, z, and 3z. The value of z should be chosen such that the percentiles can be reliably estimated. Slifker and Shapiro  suggested a choice of z = 0.524, which provides reliable estimates of four percentiles and works well for a variety of data sets.

Step 2. Estimate the percentile travel times at [x.sub.-3z], [x.sub.-z], [x.sub.z], and [x.sub.3z] based on the probabilities [P.sub.-3z], [P.sub.-z], [P.sub.z], and [P.sub.3z] and calculate the parameters m, n, and p as

[mathematical expression not reproducible]. (1)

Slifker and Shapiro  showed that the expression mn/[p.sup.2] can be used to determine the appropriate Johnson curve. The Johnson [S.sub.U] corresponds to mn/[p.sup.2] > 1, the Johnson [S.sub.B] corresponds to mn/[p.sup.2] < 1, and Johnson [S.sub.L] corresponds to mn/[p.sup.2] = 1. It should be noted that when using the value of mn/[p.sup.2] to select among Johnson curve types, [S.sub.U] and [S.sub.B] predominate since only in rare cases mn/[p.sup.2] will be exactly equal to 1.

Step 3. Estimate the Johnson curve parameters [gamma], [eta], [epsilon], and [lambda]. For the detailed calculation of these parameters, readers can refer to literature [24, 25].

Based on the above estimated parameters, travel time rate distributions using the Johnson curves can be derived straightforwardly. It is worthwhile noting that the circumstances in which mn/[p.sup.2] exactly equals 1 are rather rare. Therefore, [S.sub.U] and [S.sub.B] are more frequently chosen when the value of mn/[p.sup.2] is used to determine the type of Johnson curves. However, mn/[p.sup.2] could be less than or larger than 1 when samples are extracted from an [S.sub.L] distribution in some data sets . In this case, examination of the third (skewness) and fourth (kurtosis) moments of data in relation to the moment plane could help determine the proper type of Johnson curves when mn/[p.sup.2] is close to 1. For instance, as suggested by Mage , if the (third and fourth) moments lie within the [S.sub.B] region and mn/[p.sup.2] is larger than 1, the [S.sub.L] distribution is an appropriate choice.

2.3. Illustration with Field ANPR Data

2.3.1. Study Area and Data Description. In order to investigate how the Johnson curve system performs in terms of describing network travel time distribution (variability), we selected a study network located in the CBD of Changsha City, China. Figure 2 shows the selected network with 73 signalized intersections and 21 intersections with ANPR cameras. The data were collected for three days from 20 April 2015 to 22 April 2015. The time period of analysis considered in this study is between 7:00 a.m. and 7:00 p.m. Each day was divided into 30 min time intervals. OD travel times between intersections with ANPR can be estimated from camera records. Outliers such as double counting and mismatching were eliminated from the data set. Apart from these outliers, by comparing travel times between vehicles consecutively passing the intersection, abnormal link travel times between two adjacent cameras (e.g., parking along the road for a coffee) were identified and eliminated as well. The process can be described as follows.

Vehicles i -1, i, i + 1, ..., n sequentially pass intersections j, j + 1, ..., m. Let us denote [[??].sup.j.sub.i] as the time instant when vehicle i passes the upstream intersection j, [[??].sup.j.sub.i - 1] as the time instant when the immediate preceding vehicle i - 1 passes the upstream intersection j, [[??].sup.j.sub.i + 1] as the time instant when the immediate following vehicle i + 1 passes the upstream intersection j, [[??].sup.j + 1.sub.i - 1] as the time instant when vehicle i - 1 passes the adjacent downstream intersection j + 1, [[??].sup.j + 1.sub.i] as the time instant when vehicle i passes the adjacent downstream intersection j + 1 ,and [[??].sup.j + 1.sub.i + 1] as the time instant when vehicle i + 1 passes the adjacent downstream intersection j + 1; then the link travel time between intersections j and j+1 for vehicles i - 1, i, and i + 1 can be calculated as

[mathematical expression not reproducible]. (2)

The link travel time differences between vehicle i - 1 and i and vehicle i and i + 1 are calculated as

[mathematical expression not reproducible]. (3)

If [DELTA][[tau].sup.-] > [delta] and [DELTA][[tau].sup.+] > [delta], vehicle i is considered as "abnormal" vehicle whose travel time is much longer compared with "normal" preceding vehicle i - 1 and following vehicle i + 1. Then the link travel time of vehicle i is eliminated, where [delta] is the predefined threshold (the average cycle time of the intersections is about 3 min; therefore, "3 min" is chosen in our case).

The vehicle travel time rates were then calculated with individual travel times divided by their travel distance.

2.3.2. Results. We calculated the value of parameters [gamma], [eta], [epsilon], and [lambda] in the Johnson curve system according to the steps described in the previous subsection. The value of mn/[p.sup.2] is smaller than 1 for all cases, which indicates that the Johnson [S.sub.B] curves are appropriate to describe travel time rate distributions for our selected network.

Figures 3(a)-3(c) illustrate the comparison between the network travel time rate distributions estimated from ANPR data and the Johnson fitted and lognormal fitted curves during peak time period of 7:00-7:30 a.m. for three days. The Johnson curves show a better fit to the ANPR travel time data compared with the lognormal fit. Similar results can be identified in Figures 3(d)-3(f), where ANPR travel time rate distributions are presented along with Johnson and lognormal fitted probability curves during off-peak period of 10:00-10:30 a.m. We further applied the one-sample Kolmogorov-Smirnov (KS) test with four time periods each day (in total, 12 time periods). Table 1 presents the K-S test results. From Table 1, it can be concluded that the null hypothesis that the travel time rates follow the Johnson curve (SB) distribution cannot be rejected for most of the time periods (except time periods 5:00 p.m.-5:30 p.m. on day 2015-0421) at the significance level of 5% since the P value is larger than 0.05. The better fitting performance of the Johnson curve model is due to the fact that the Johnson system allows describing a wide spectrum of traffic conditions as usually existing in the urban traffic system. The Johnson model can flexibly adapt its shape to the observed travel time data. We calculated the Akaike information criterion (AIC) for both the Johnson curve model and the lognormal model as well. AIC is a measure of the relative quality of statistical models for a given set of data. It is defined as

AIC = 2 * k - 2 ln (L), (4)

where k is the number of parameters in the model; L is the maximized value of the likelihood function for the model. From (4), it can be seen that AIC rewards goodness of fit (as assessed by the likelihood function), but it also includes a penalty that is an increasing function of the number of estimated parameters.

Given a set of candidate models for the data, the preferred model is the one with the minimum AIC value. For more detailed information about AIC, readers can refer to Akaike's work . Table 2 shows the AIC values for both models in different time periods. The lower AIC values of the Johnson curve model further confirm that the Johnson curve model outperforms the lognormal distribution model.

Besides, as can be seen from Figure 3, network travel time distributions are wide and skewed with a long fat right tail during different time periods, which suggest that the networktravel time rate is rather variable (unreliable). In the following section, we will look at network travel time reliability measures and how these measures vary under different traffic conditions.

3. Network Travel Time Reliability Metrics and Their Connection with Network Traffic Characteristics

3.1. Network Travel Time Reliability Measures. Based on travel time distributions as proposed in the previous section, any travel time reliability metrics (e.g., standard deviation, percentiles, skewness, and buffer index) can be derived. In the past decades, travel time reliability metrics are commonly used to quantify link or route travel time reliability, whereas very little research has been dedicated to investigate the network travel time reliability. In this section, we investigate the network travel time reliability (in terms of intersection-weighted standard deviation and intersection-weighted skewness of travel time rate) and its connection with network traffic characteristics under different traffic conditions.

A method that is frequently used in literature to characterize traffic states in a network is the macroscopic fundamental diagram (MFD) [31-33]. In the MFD, the traffic state is characterized by production (sum of distance-weighted traffic flow) and accumulation (sum of distance-weighted traffic density). Let us denote the length of a link j by [l.sub.j], the outflow of link j by [q.sub.j], the traffic density (vehicles per length of the road) of link j by [k.sub.j], and the total number of links in the network by n. Then, the weighted average network density during time interval [DELTA]t is defined as

[mathematical expression not reproducible]. (5)

Since we use ANPR data, the traffic density cannot be derived directly. We calculate the traffic density from the total vehicle travel time on link j and the link length [l.sub.j]. We use the methodology proposed by Edie  applied on single links, based on trajectory analysis . The average traffic density of link j [k.sup.j] during time interval [DELTA]t is estimated as the total time spent on the link during the time interval.

There is a methodological problem in the choice of the frame in space and time for the definition of the traffic density. The time the vehicle i spent on the link j [t.sub.ij] for the first few vehicles passing the end of the link after [t.sub.1] is not entirely spent during the time period [DELTA]t as shown in Figure 4. The situation is represented in a rather extreme way to show the problem. The first four vehicles registered by the camera have a travel time which has not been spent completely in the time frame [DELTA]t. The next four vehicles are not registered in time frame [DELTA]t although they spend a part of the time in the space-time frame [DELTA]x x [DELTA]t. Although the detailed trajectories are not known, we can apply the Three-Dimensional Trajectory method (Saberi et al. ). To be accurate we have to correct the travel time of the first vehicles passing in [DELTA]t, that is, the vehicles for which the travel time is longer than the time in [DELTA]t at which they are registered. Similarly, we have to add at the end of the timeframe the times spend by the vehicles that entered [DELTA]x but did not reach the end of the road section within the timeframe.

The travel time of vehicle i on link j is given by the difference in passing time at the stop line of the link j [[??].sup.t.sub.ij] and the link [mathematical expression not reproducible]. The part of the travel time that should be allocated to the time interval [DELTA]t is calculated as

[mathematical expression not reproducible]. (6)

In an analogue way, the vehicles that are not counted at the end of the time frame can be detected in the next time frame and partly allocated to [DELTA]t. The travel time that has to be added to the time spent in the period [DELTA]t is calculated as

[mathematical expression not reproducible]. (7)

In practice, if the time intervals [DELTA]t are chosen longer than the travel time on the link the boundary problem gives only a marginal influence of the calculation of the time spent. Furthermore, At is chosen long enough to have little influence of the signal cycle time (which was not a common cycle time for all intersections). This avoids the time-space window problems as discussed by Courbon and Leclercq 

[k.sub.j] ([DELTA]t) = [[SIGMA].sup.m.sub.i = 1] (t[t.sub.'ij]/[l.sub.j])/[DELTA]t ([DELTA]t = [t.sub.2] - [t.sub.1]), (8)

where t[t'.sub.ij] is the corrected travel time of vehicle i on link j passing the stop line at the end of link j between [t.sub.1] and [t.sub.2]; m is the total number of vehicles leaving link j during time interval [DELTA]t united with the vehicles that entered the link j in [DELTA]t; n is the total number of links in the network; [t.sub.1] and [t.sub.2] are the start and end of the time interval.

The weighted (average) network flow during time interval [DELTA]t is defined as

[mathematical expression not reproducible], (9)

where [q.sub.j]([DELTA]t) is determined as the number of vehicles passing the end of link j during [DELTA]t. Apart from the network flow and density, three reliability related quantities which are intersection-weighted mean travel time rate, intersection-weighted standard deviation of travel time rate, and intersection-weighted skewness of travel time rate distributions are used in this study.

As we proposed in the previous section, the travel time rate of trip i can be calculated as the trip travel time [t.sub.i] divided by the number of intersections ni during this trip:

[[tau].sub.i] = [t.sub.i]/[n.sub.i]. (10)

The intersection-weighted mean travel time rate during time interval [DELTA]t is defined as

[mathematical expression not reproducible], (11)

where [n.sub.i] is the number of intersections during trip i, [[tau].sub.i] is the travel time rate of trip i, and N is the total number of trips observed in the network during time period [DELTA]t.

The intersection-weighted standard deviation of travel time rate during time interval [DELTA]t is defined as

[mathematical expression not reproducible]. (12)

Skewness is a measure of the asymmetry of the probability distribution of a random variable about its mean. In terms of reliability, the more skewed the travel time distribution is, the less reliable travel times become. If the value of skewness is zero, the travel time is normally distributed. If the value of skewness is larger than zero, the travel time distribution is right-skewed with a longer tail to the right. Otherwise, the travel time distribution is left-skewed with a longer tail to the left. The weighted skewness of travel time rate for a network during time interval At is defined as

[mathematical expression not reproducible], (13)

where [mu]([DELTA]t) is the mean travel time rate during time interval [DELTA]t; [sigma]([DELTA]t) is the standard deviation of travel time rate during time interval [DELTA]t.

3.2. Case Study with Field ANPR Data. The selected study area consisting of 54 links is shown in Figure 2 in the previous section. Again we use the three-day (between 7:00 a.m. and 7:00 p.m. from 20 April until 22 April 2015) ANPR data from Changsha City, China. The traffic flow can be derived directly since each ANPR camera records the plates number of all vehicles on lanes for straight-on and left turning movements passing the stop line. The traffic density however cannot be obtained directly from ANPR data. We apply (6) to estimate traffic density for each link. Figure 5 illustrates the macroscopic fundamental diagram for the selected CBD network based on three days' observations with time intervals of 5 min (Figure 5(a)) and 15 min (Figure 5(b)), respectively. This macroscopic fundamental diagram reveals that the observed maximum weighted traffic volume in this network is around 400 veh/h/lane for the weighted traffic density around 30 veh/km/lane. No gridlock happened in this region from the observed data. That does not mean that the network operates below capacity, but it is an indication that the traffic police actively takes measures to reduce inflow and increase outflow on links at the moment that spillback threatens to occur. Of course, this form of active traffic management has an influence on travel times and travel time reliability.

Figure 6 illustrates the relation between the weighted average travel time rate and the weighted standard deviation of travel time rate for every 5 min. Figures 6(a)-6(c) are intersection-weighted travel times and Figures 6(d)-6(f) are distance-weighted travel times. A clear linear trend can be observed with both intersection-normalized and distance-normalized approaches, though the travel time data appear to be more scattered with a slightly lower correlation coefficient ([R.sup.2] = 0.7167 and 0.5938) on 21 April. Higher correlations can be observed with intersection-weighted approach compared with the distance-weighted method for all three days' travel times. The magnitude of the slope which reflects how much the weighted standard deviation will increase when the weighted average travel time rate increases is between 0.46 and 0.65 for intersection-weighted approach and between 0.35 and 0.62 for distance-weighted approach within these three days' observations. The larger this value is, the more vulnerable to congestion induced reliability deterioration the network becomes.

Figure 7 shows the relation between the weighted average travel time rate and the network density, where the upper three figures (Figures 7(a)-7(c)) are intersection-weighted travel time rates and the lower three ones (Figures 7(d)-7(f)) are distance-weighted travel time rates. A relatively strong positive correlation exists between the weighted average travel time rate and the network density for all three days with [R.sup.2] > 0.55 for intersection-weighted method and [R.sup.2] > 0.48 for distance-weighted approach. The comparison between the intersection-weighted method and the distance-normalized method illustrates that closer correlations can be observed with the former than the later. Lower correlation coefficients can be observed for the weighted standard deviation versus the network density as shown in Figure 8. Nevertheless, the positive linear trend still exists, which indicates that the travel time variability increases when the network density increases. This is consistent with the argument in literature that the network travel time variability has a linear relation with the network density with simulated data .

Apart from investigating network travel time variability in terms of the standard deviation, the intersection-weighted skewness of ANPR travel time data is calculated as well according to (15). Figure 9 shows the relation between the weighted skewness and the network density. There is no strong correlation between the network density and skewness, though a weak trend, indicating that the skewness decreases with the increase of network density, can be observed. Nevertheless, we can observe that most of the values of skewness between 7:00 a.m. and 7:00 p.m. are larger than zero, which indicates that the travel time distribution is right-skewed with a longer right tail. This is consistent with travel time distributions derived in Section 2 with Johnson curves.

The relations between travel time reliability and network flow are illustrated in Figure 10. It is observed that the relationship between travel time reliability (e.g., the intersection-weighted standard deviation of travel time rate) and network flow rate is not monotonic anymore. The intersection-weighted standard deviation of travel time rate increases with network flow rate in low traffic demand condition until the maximum flow rate, beyond which the network flow rate drops while the standard deviation continues to increase. The results shown in Figure 10 are quite consistent with those derived from simulated data in .

4. Discussion and Conclusions

This paper proposes a network-level travel time distribution model by applying a flexible system of Johnson curves. The proposed model is able to describe a variety of travel time distributions with properly estimating model parameters based on the observed percentiles in the travel time data. We validate the performance of the proposed model using ANPR data collected in the CBD area of Changsha City, China. The results indicate that the travel time data can be better modelled with the proposed model than the widely used lognormal model in the transportation field. The reason behind is that the Johnson curve model is able to better capture the skewness in the data set with parameters describing the tail of the distribution (e.g., with the third and fourth moments or four percentile points), while in the lognormal distribution, only the mean and standard deviation of the logarithm are used to characterize the distribution of the data set.

We further investigated the network travel time reliability measures and their connection with network traffic characteristics in terms of network density and network flow, which are known as the macroscopic fundamental diagram. Three network travel time reliability measures which are the weighted average travel time rate, the weighted standard deviation of travel time rate, and the weighted skewness of travel time distributions are discussed, and the relation between these reliability measures and the network density is illustrated using the field ANPR data. The weighting is done with respect to both the number of intersections on a trip and the trip distance. Results show that the weighted average travel time rate has a positive linear relation with the weighted standard deviation of travel time rate for different time periods under time-varying demand. A positive linear trend can be clearly observed between the weighted standard deviation (or the weighted average travel time rate) and the network density. The correlations between reliability measures and network states are higher with the intersection-weighted approach than those using the distance-weighted approach. The investigation of the relation between the weighted skewness and the network density shows that there is a weak negative correlation between the former and the latter based on the available ANPR data. The value of skewness is larger than 0 for most observation periods between 7:00 a.m. and 7:00 p.m., which indicates that travel time distributions are right-skewed with a longer tail to the right.

The findings in this paper provide some insights into the network travel time reliability. The empirical relation between the network travel time reliability and network traffic characteristics (e.g., network density) can be possibly applied to assess the traffic management and control measures for reliability purposes; for example, by applying the relations between network flow and network density, as well as network reliability (network standard deviation) and network density, network control strategies can be optimized to improve network reliability while maintaining required network throughput.

https://doi.org/10.1155/2017/9147356

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This work is supported by the National Science Foundation of China (NSFC) under Projects 51308475, 61673321, and 71671147, the Fundamental Research Funds for the Central Universities under Contract 2682016CY07, Chengdu Science and Technology Commission under Contract 2015-RK00-00190-ZF, and Changsha Science and Technology Commission under Projects K1001010-11 and K1106004-11.

References

 F. Zheng, H. van Zuylen, and Y. Pu, "Model for predicting distribution of link travel times for urban signalized roads," Transportation Research Record: Journal of the Transportation Research Board, vol. 2311, pp. 59-72, 2012.

 F. Zheng and H. Van Zuylen, "The development and calibration of a model for urban travel time distributions," Journal of Intelligent Transportation Systems: Technology, Planning, and Operations, vol. 18, no. 1, pp. 81-94, 2014.

 W. Pu, "Analytic relationships between travel time reliability measures," Transportation Research Record: Journal of the Transportation Research Board, vol. 2254, pp. 122-130, 2011.

 S. Arroyo and A. L. Kornhauser, "Modeling travel time distributions on a road network," in Proceedings of 84th TRB annual meeting, Washington DC, USA, 2005.

 H. Al-Deek and E. B. Emam, "New methodology for estimating reliability in transportation networks with degraded link capacities," Journal of Intelligent Transportation Systems: Technology, Planning, and Operations, vol. 10, no. 3, pp. 117-129, 2006.

 F. Guo, H. Rakha, and S. Park, "A multi-state travel time reliability model," in Proceedings of the TRB 89th annual meeting, Washington DC, USA, 2010.

 P. Chen, K. Yin, and J. Sun, "Application of finite mixture of regression model with varying mixing probabilities to estimation of urban arterial travel times," Transportation Research Record: Journal of Transportation Research Board, vol. 2442, pp. 96-105, 2014.

 S. Clark and D. Watling, "Modelling network travel time reliability under stochastic demand," Transportation Research Part B: Methodological, vol. 39, no. 2, pp. 119-140, 2005.

 J. Kim and H. S. Mahmassani, "Compound Gamma representation for modeling travel time variability in a traffic network," Transportation Research PartB: Methodological, vol. 80, pp. 40-63, 2015.

 J. Kim, H. Mahmassani, P. Vovsha, Y. Stogios, and J. Dong, "Scenario-based approach to analysis of travel time reliability with traffic simulation models," Transportation Research Record: Journal of the Transportation Research Board, vol. 2391, pp. 56-68, 2013.

 A. Chen and Z. Zhou, "The a-reliable mean-excess traffic equilibrium model with stochastic travel times," Transportation Research PartB: Methodologies, vol. 44, no. 4, pp. 493-513, 2010.

 J. Bates, J. Polak, P. Jones, and A. Cook, "The valuation of reliability for personal travel," Transportation Research Part E: Logistics and Transportation Review, vol. 37, no. 2-3, pp. 191-229, 2001.

 T. Lomax, D. Schrank, S. Turner, and R. Margiotta, Selecting Travel Time Reliability Measures, Texas Transportation Institute, Texas, USA, 2003.

 FHWA, Travel Time Reliability: Making It There On Time, All The Time, Federal Highway Administration.

 J. W. C. van Lint, H. J. van Zuylen, and H. Tu, "Travel time unreliability on freeways: why measures based on variance tell only half the story," Transportation Research Part A: Policy and Practice, vol. 42, no. 1, pp. 258-277, 2008.

 H. K. Lo, X. W. Luo, and B. W. Y. Siu, "Degradable transport network: travel time budget of travelers with heterogeneous risk aversion," Transportation Research Part B: Methodological, vol. 40, no. 9, pp. 792-806, 2006.

 X. Xu, A. Chen, L. Cheng, and H. K. Lo, "Modeling distribution tail in network performance assessment: a mean-excess total travel time risk measure and analytical estimation method," Transportation Research PartB: Methodological, vol. 66, pp. 32-49, 2014.

 E. A. I. Bogers, Traffic Information and Learning in Day-to-Day Route Choice [Dissertation, thesis], Transport & Planning, Delft University of Technology, 2009.

 H. Mahmassani, T. Hou, and J. Dong, "Characterizing travel time variability in vehicular traffic networks: deriving a robust relation for reliability analysis," Transportation Research Record: Journal of Transportation Research Board, no. 2315, pp. 141-152, 2012.

 H. Li, Reliability-based Dynamic network Design with Stochastic Networks. Doctoral dissertation [Doctoral, thesis], Delft University of Technology, 2009.

 M. Fosgerau and A. Karlstrom, "The value of reliability," Transportation Research Part B: Methodological, vol. 44, no. 1, pp. 38-49, 2010.

 S. Susilawati, M. A. P. Taylor, and S. V. C. Somenahalli, "Distributions of travel time variability on urban roads," Journal of Advanced Transportation, vol. 47, no. 8, pp. 720-736, 2013.

 N. L. Johnson, "Systems of frequency curves generated by methods of translation," Biometrika, vol. 36, pp. 149-176, 1949.

 J. F. Slifker and S. S. Shapiro, "The johnson system: Selection and parameter estimation," Technometrics, vol. 22, no. 2, pp. 239-246, 1980.

 N. R. Farnum, "Using Johnson curves to describe non-normal process data," Quality Engineering, vol. 9, no. 2, pp. 329-336, 1996.

 E. G. Jones, H. S. Mahmassani, R. Herman, and C. M. Walton, "Travel time variability in a commuting corridor: implications for electronic route guidance," in Proceedings of the the First International Conference on Applications of Advanced Technologies in Transportation Engineering, 1989.

 I. D. Hill, R. Hill, and R. L. Holder, "Algorithm AS 99: fitting Johnson curves by moments," Applied Statistics, vol. 25, no. 2, p. 180, 1976.

 I. D. Hill, "Algorithm AS 100: Normal-Johnson and Johnson_ Normal transformations," Applied Statistics, vol. 25, no. 2, p. 190, 1976.

 D. T. Mage, "An explicit solution for sbparameters using four percentile points," Technometrics, vol. 22, no. 2, pp. 247-251, 1980.

 H. Akaike, "A New Look at the Statistical Model Identification," IEEE Transactions on Automatic Control, vol. 19, no. 6, pp. 716-723, 1974.

 N. Geroliminis and C. F. Daganzo, "Existence of urban-scale macroscopic fundamental diagrams: some experimental findings," Transportation Research Part B: Methodological, vol. 42, no. 9, pp. 759-770, 2008.

 C. F. Daganzo and N. Geroliminis, "An analytical approximation for the macroscopic fundamental diagram of urban traffic," Transportation Research Part B: Methodological, vol. 42, no. 9, pp. 771-781, 2008.

 Y. Ji, N. Geroliminis, and J. Luo, "Empirical observations of congestion propagation and dynamic partition with probe data for large scale systems," in Proceedings of the 93rd TRB annual meeting, Washington DC, USA.

 L. C. Edie, "Discussion of traffic stream measurements and definitions," in Proceeding of the 2nd International Symposium on the Theory of Traffic Flow, J. Almond, Ed., pp. 139-154, Organisation for Economic Co-operation and Development, 1963.

 M. Saberi, H. S. Mahmassani, T. Hou, and A. Zockaie, "Estimating network fundamental diagram using three-dimensional vehicle trajectories: extending edie's definitions of traffic flow variables to networks," Transportation Research Record: Journal of the Transportation Research Board, vol. 2422, pp. 12-20, 2014.

 T. Courbon and L. Leclercq, "Cross-comparison of macroscopic fundamental diagram estimation methods," in Proceedings of the Procedia Social and Behavioral Sciences, vol. 20, pp. 417-426, September 2011.

 H. Mahmassani, T. Hou, and M. Saberi, "Connecting networkwide travel time reliability and the network fundamental diagram of traffic flow," Transportation Research Record: Journal of the Transportation Research Board, no. 2391, pp. 80-91, 2013.

Fangfang Zheng, (1, 2) Xiaobo Liu, (1, 2) Henk van Zuylen, (1, 2, 3, 4) Jie Li, (4) and Chao Lu (1, 2)

(1) School of Transportation and Logistics, Southwest Jiaotong University, No. 111 Erhuanlu Beiyiduan, Chengdu 610031, China

(2) The National United Engineering Laboratory of Integrated and Intelligent Transportation, Southwest Jiaotong University, No. 111 Erhuanlu Beiyiduan, Chengdu 610031, China

(3) Transport and Planning Department, Delft University of Technology, P. O. Box 5048, 2600 GA Delft, Netherlands

(4) Civil Engineering College, Hunan University, Lushan South Road, Changsha, Hunan Province 410082, China

Correspondence should be addressed to Xiaobo Liu; xiaobo.liu@swjtu.cn

Received 9 March 2017; Revised 17 May 2017; Accepted 4 June 2017; Published 9 July 2017

Caption: FIGURE 1: Variability of travel times from empirical data collected on an urban arterial road "Renmin Road" on 20 April 2015 (flow and travel time are both measured in 5 min aggregation; the outliers in travel time have been filtered out) in Changsha City, China.

Caption: FIGURE 2: The selected network in Changsha City, China.

Caption: FIGURE 3: Network travel time rate distribution from ANPR data and fitted density curves at morning peak-hour period ((a)-(c)) and morning off-peak period ((d)-(f)).

Caption: FIGURE 4: Trajectories at the boundary of the time frame.

Caption: FIGURE 5: Network flow-density diagram for the study area in Changsha City (7:00 a.m.-7:00 p.m.).

Caption: FIGURE 6: Weighted average travel time rate versus weighted standard deviation of travel time rate.

Caption: FIGURE 7: Weighted average travel time rate versus network density.

Caption: FIGURE 8: Weighted standard deviation of travel time rate versus network density.

Caption: FIGURE 9: Intersection-weighted skewness versus network density.

Caption: FIGURE 10: Intersection-weighted standard deviation of travel time rate versus network flow.
```TABLE 1: P value of Kolmogorov-Smirnov (KS) test with significance
level of 0.05.

P value

Fitted distributions        Time period        Day 1   Day 2   Day 3

Johnson curve           7:00 a.m.-7:30 a.m.    0.456   0.334   0.284
10:00 a.m.-10:30 a.m.   0.312   0.265   0.202
3:00 p.m.- 3:30 p.m.    0.258   0.168   0.126
5:00 PM-5:30 p.m.     0.359   0.032   0.243
Lognormal fit           7:00 a.m.-7:30 a.m.    0.392   0.036   0.122
10:00 a.m.-10:30 a.m.   0.087   0.032   0.071
3:00 p.m.- 3:30 p.m.    0.036   0.073   0.048
5:00 p.m.-5:30 p.m.    0.033   0.059   0.076

Day 1: 2015-04-20; day 2: 2015-04-21; day 3: 2015-04-22.

TABLE 2: AIC values for the Lognormal model and the Johnson-curve
model.

AIC value

Fitted distributions        Time period        Day 1   Day 2   Day 3

Johnson curve           7:00 a.m.-7:30 a.m.    15823   16988   17154
10:00 a.m.-10:30 a.m.   16343   16772   16952
Lognormal               7:00 a.m.-7:30 a.m.    16128   18214   17993
10:00 a.m.-10:30 a.m.   16766   17525   17629
```
Title Annotation: Printer friendly Cite/link Email Feedback Research Article Zheng, Fangfang; Liu, Xiaobo; van Zuylen, Henk; Li, Jie; Lu, Chao Journal of Advanced Transportation Report Jan 1, 2017 7331 Suspended Load Path Tracking Control Strategy Using a Tilt-Rotor UAV. Dynamic Path Planning of Emergency Vehicles Based on Travel Time Prediction. Mathematical models Transportation research Travel time (Traffic engineering) Urban transportation