1. Introduction
Globally there are 1.25 million deaths caused by traffic accidents every year [1]. In Colombia, there were more than 46,000 injuries for this cause in 2017. Specifically, public transportation (buses) is responsible for 3,100 of these cases [2], which is a worrying indicator. As a response, governments are turning to policies that encourage safety in public transit.
Bus lanes are a relatively recent alternative in the "push and pull" strategies [3] for promoting travelers to use public transport instead of a private car. Considering the future relevance of such types of lanes in the public transit planning, management, and operations, it is crucial to evaluate them in terms of various components such as operation, safety, environmental externalities among others. In terms of operation, there is a great amount of research developed [4, 5]. Some literature deal with the definition of different bus lane typologies and the way they operate [6, 7]. Other studies analyzed improving operation by evaluating some attributes like bus average speed and queue lengths before and after the implementation of these lanes [8]. However, there is limited literature on safety assessment in this kind of infrastructure and they are mainly based on accident-injury data. On the other side, studies of road safety were done at intersections [9], road links [10] or roundabouts [11], among others, but there is still a gap in this assessment type on bus lanes and public transportation.
The research described in this paper consists in the road safety analysis of a preferential bus lane, also known as preferential corridor, using microsimulation with modeling software VISSIM 9 and the Surrogate Safety Measures Assessment Model (SSAM 3.0). This study includes transit and behavioral driver parameters that give major accuracy and certainty to the calibration and validation of the model. The operating conditions and road safety are characterized in a section of the 7th Avenue in Bogota which has a preferential bus lane. The Swedish Traffic Conflicts Technique (STCT) is used to validate the conflict number of observation-based data with microsimulation results through SSAM as a methodology of quantitative risk analysis in road safety, in the infrastructure type described above. Three control countermeasures are proposed to mitigate the risk situations that are detected.
The rest of the paper is structured as follows: First, a literature review of bus lane operation and evaluation of road safety in different traffic interactions, the use of conflicts technique and surrogate safety measures. Then, a description of the study case is presented, including a deep explanation of the bus lane type studied and the theoretical-practical operation conditions. The above is followed by a description of the methodology, detailing the calibration process used. Finally, the results of microsimulation and calibration are presented, including results of three different control countermeasures in hypothetical conditions. The results are discussed in terms of operation and road safety, based on the reduction of the conflicts number and safer surrogate safety measures.
2. Literature review
2.1 Bus lanes
Different types of bus lanes were developed around the world with particular operational characteristics. Dedicated bus lanes with or without physical barriers are commonly used for BRT systems. The term “intermittent bus lane” was introduced in [6] and defined as a lane which acts as a dedicated bus lane only in presence of buses. In other cases, it is open to mixed traffic. Other authors [12] refer to this as “bus priority”, a more general term that involves multiple measures related to a strategy for improving bus operation and enhance environment for users [13]. In this case, operating conditions must be clearly explained at each research study.
In terms of operation of bus lanes, there is a great amount of literature. Different attributes of a corridor with an Intermittent Bus Lane in Lisbon, Portugal, were studied by [8]. They found increases in bus average speed and no significant variation in parameters like other vehicle speeds and queue lengths. In other study [14], a fitting between pre-evaluated data in a Paramics model and after-constructed data on an urban expressway bus lane in Beijing, China, was found. Additionally, they reported improved capacity in the same study area. Other research presents study cases of different cities around the world [15] including operation monitoring [16] and signaling analysis [17], among other topics.
Focusing on road safety dedicated bus lanes, limited literature was found. [18] studied the crash risk between two-wheeler drivers and buses in bus lanes, considering a bodily injury accidents-based risk measure. Furthermore, [19] research about maneuvers to analyze close interactions between bicyclists and buses in bus lanes. [20], conducted a study through empirical analysis of accident data, and used an approximation based on microsimulation. A first approach in this field [20], was developed through the implications of an exclusive bus lane. The research was conducted in Melbourne, modelling three scenarios (with AIMSUN) where two of them have an exclusive bus lane and the third one considers mixed traffic. The calibration of this model was based on the two-step calibration model proposed by [21] and used the deceleration rate required to avoid a crash (DRAC) and a crash potential index (CPI) as surrogate measures for assessing road safety.
2.2 Traffic conflict techniques
Studies in road safety include estimation of measures that make the assessment possible. Owing to the lack of crash data to apply studies of road safety assessment, researchers have focused on replicate the behavior of drivers with surrogate events. The approximation used in this research are traffic conflicts, which [22] are developed through the Swedish Traffic Conflicts Technique (STCT).
A traffic conflict is defined as an incident where, if two road users on collision trajectory maintain their speed and trajectory constant, the crash between them will be imminent. Based on this concept, there are different physical variables, used as indicators of how close a collision would have happened, as well as the severity of the traffic conflict. Such indicators are known as surrogate safety measures. To avoid the collision, it is necessary that, at least one of the two road users perform an evasive action. The user who first performs this action is called the "relevant user" [9].
As described in [23], many different surrogate safety measures or indicators are used in the road safety assessment. According to [23] , the most common measures used are: TTC (time to collision), defined as the time remaining to the collision if the invaded vehicle maintain the speed it currently has. TTC can be calculated in each moment of the risky situation using the speed and the remaining distance to the collision. In the same sense, TA (time to accident) is a critical value of the TTC function that occurs when the invaded vehicle initiates an evasive action; PET (post-encroachment time) is the real time measured between the moment of the trajectory invasion and the moment when the invaded vehicle reaches the eventually collision point, and DR (deceleration rate) is the initial deceleration at the beginning of the conflict. A combination of these measures were used to develop some traffic conflict techniques as "a framework for observed-based safety studies" [23].
2.3 The Swedish Traffic Conflict Technique - STCT
The STCT, first used by [22], estimates the severity of a conflict based on two surrogate measures:
Time to accident (TA): Time remaining to the collision, foreseen when the relevant user initiates the evasive action.
Conflicting speed (CS): Speed of the relevant user at the moment when the evasive action starts.
Figure 1 shows the diagram used to distinguish between serious and minor conflicts. Conflicts with severity level 26, or higher, are considered serious; however, there are some indications in the literature that suggests that the limit should be reduced one or two levels minor [24]. This information is collected in the field, through the recurrent observation of videos. The method allows the identification of events that can be considered as conflicts; these are analyzed accordingly using STCT framework, to estimate their severity and then decide if they should be included in the safety assessment.
2.4 STCT and microsimulation
The use of simulation models including traffic variables, behavioral parameters and transit operation to assess road safety, using the trajectory files of simulation models through the open tool Surrogate Safety Assessment Model (SSAM), is proposed by [25]. Therefore, SSAM aims to assess traffic safety through the number of simulated conflicts which are complemented with a characterization depending on the severity of some surrogate measures. These indicators are: time to collision (TTC), post encroachment time (PET), deceleration rate (DR), maximum speed (MaxS) and speed differential (DeltaS), as shown in Figure 2. These variables are reported by SSAM and their validity depends on the calibration of the simulation model, both by traffic parameters and by conflicts.
3. Case study
The study was conducted in a corridor of an important avenue in Bogota, Colombia, limited to approximately 700 m length of the road. Two important aspects make attractive this part of the city to be selected for our study: First, it has a preferential bus lane and, second, the sector presented the highest number of reported crash data, from 2014 to June 2018 [26]. The preferential bus lane was implemented in 2014 and consists in an external lane for public transport (buses) with visual segregation through horizontal and vertical signalization, as shown in Figure 3. Private transport can use the preferential lane just in case of right turns to access to intersections or properties. Through observation in situ, it was detected that the road users often use the bus lane in a different way than intended, this is that private vehicles often take up this lane causing a risky operation.
A possible source of traffic conflicts is not only the operation of the preferential bus lane, also the existence of bus stops for different routes and local public transport system’s components. Additionally, there are some important elements which made the study area’s operation more complex in terms of driving behavior. Such elements are entrances and exits to an educational institute, a hospital, and a gas station, among others. The general location of the study area and roads that currently have a preferential lane are shown in Figure 4. Figure 5 shows the study area with the bus lane the components related to the complex operation, and the location of traffic accidents during the last 5 years. The study period was determined as the hour with the highest number of observed conflicts. The operation was video recorded for 12 continuous hours, during a working day and a Saturday, finding the riskiest period to be from 07:00 to 08:00 on the working day.
4. Methodology
The methodology is divided into various steps. In the next section, the development of the base simulation model and the parameters used for its calibration are explained. The first step consists in the development of VISSIM simulation model. The model is carried out to represent the traffic and behavior conditions observed in field, and to evaluate the road safety of the corridor, according to [21]. After that, it is necessary to calibrate the model through some traffic variables as speed and volumes, and through safety variables as traffic conflicts. To validate the model from the point of view of road safety, the STCT is used to identify serious conflicts, in the field, looking for the simulation model to reproduce the operation observed in terms of the simulated conflicts. Subsequently, the calibrated and validated base model is characterized through traffic and safety variables. Finally, under the calibrated environment, control countermeasures are proposed and simulated, assessing the improvement of road safety.
4.1 Development of VISSIM simulation model
The base model was developed in VISSIM 9. The network was modeled based on a satellite image taken from official sources [27]. The network includes a scheme of 7th Avenue from 39th Street to 45th Street; links and connectors which represent real conditions, as shown in Figure 6. The model includes the preferential bus lane and stops for the local public transport system. Figure 7 shows priority rules and conflict areas that were designed and configured to avoid vehicle overlapping that could lead to the emergence of virtual crashes.
The initial model calibration consisted in the comparison between simulated and field-observed traffic variables. We evaluated volumes and speeds [21] to ensure that the simulation model correctly represents the number of vehicles circulating on the network and the speed at which they circulate. In the base model, volumes and free-flow speeds were evaluated by means of the statistical GEH [28].
At the end of the second calibration level, the VISSIM software trajectory files are extracted to be used in SSAM. This tool is used to filter simulated traffic conflicts and to obtain their surrogate measures necessary to calibrate the base model. The process has two stages [21]. The first corresponds to the comparison of the observed and simulated number of conflicts, validated by calculating the mean absolute percentage error (MAPE), as is shown in Equation (1):
where,
= field observed number of conflicts.
= simulated number of conflicts.
n = number of simulations runs. In this case n=10.
The initial parameters that are configured in SSAM are TTC, PET and conflict type by angle, as defined in [29]. The initial values by default of TTC and PET are 1.5s and 5s, respectively. After analyzing the trajectory files, filters are applied to the surrogate measures, according to field observations to calibrate and validate the number of conflicts. The filters are defined in this order: TTC, PET, DR and MaxD. If the expected results are not obtained, behavior parameters of the model are adjusted, until reaching validation criteria.
4.2 Control countermeasures proposed
Once the base model is validated, different scenarios are proposed and modeled trying to improve road safety. Three possible operating scenarios are proposed as engineering countermeasures for the better operation of the studied preferential bus lane. The first scenario imposes a strict application of the regulations [30]. Private cars are forced to avoid the bus lane, and only used them on the cases explained previously. The second scenario includes a barrier for the same bus lane, a physical segregation with an open space only at the bus stops to allow overtaking. In Madrid (Spain) there are segregated bus lanes, as shown in Figure 8, with an operation similar to the proposed scenario.
For the third scenario, the proposal is totally different. In this case, a physically segregated internal lane is recommended, which imply left ascending and descending of passengers. Figure 9 shows the scheme for this scenario.
5. Results and discussion
5.1 Calibration results
As mentioned above, the calibration and validation are reached based on three parameters: speeds, volumes and traffic conflicts. Free-flow speeds were measured in field for each type of vehicle: Taxis, light duty cars, buses and motorcycles. The statistical parameter used to validate the model was the GEH index, whose results are satisfactory when GEH < 5 at least in the 85% of the screen lines, for urban areas and a count comparison at individual turning movements and / or directional link counts [28]. Free flow speeds were calculated in the model using the data collection point option in two points of the road where free flow speeds were expected to occur. Comparison with field-measured speeds allowed the calculation of the GEH statistical. Results are shown in Table 1. On the other hand, GEH statistical was computed for volumes of fifteen different specific vehicle types. Results were obtained for eight control points in the corridor, each one for four time periods between 07:00am, 07:15, 07:30, 07:45 and 08:00, getting 480 GEH statistical values. Only 16 of those were values above 5, which was the accepted threshold found in the literature. This reveals a 96.7% of fitting values for the GEH statistical calibration.
5.2 Scenarios comparison
Comparing the base scenario with the three-proposed safety improved scenarios, two operational variables (latent demand and queues) and two road safety-related variables (amount of traffic conflicts and surrogate safety measures), are discussed. Evaluation of operational variables is important in this context, given that we aim to improve road safety without negatively affecting the operation.
Latent demand
The latent demand is calculated in control points where the main inputs and outputs affect the operation of the network. These points correspond to the approaches of the main intersections (39th and 45th Streets), before the traffic lights. The latent demand remains the same in the W-E and E-W directions, because those are not affected by the bus lane operation. When considering S-N direction, we report an increase in latent demand in scenario 3, at both intersections, which is maintained throughout the study hour. On the other hand, the N-S direction presents a slightly different behavior. When vehicles approach from 45th Street, capacity decreases along time in the proposed scenarios. However, the latter occurs much more in scenario 3. The number of vehicles going south from 39th Street has the same characteristics, but the reduction in capacity is lower. Table 2 and Table 3 show the latent demand in these approaches, for each scenario, in 15-minute intervals. In scenarios 0, 1 and 2, for S-N direction, we obtained a similar total latent demand. In the first two scenarios, there is a greater demand compared to the base scenario. In contrast, scenario 3 presents a considerable increase in latent demand when compared to the other scenarios. These results show that the lower the latent demand on the network, the greater capacity to support the flow.
For N-S direction, all the proposed scenarios have a higher latent demand value when compared to scenario 0. The demand of scenarios 1 and 2 is close to the baseline scenario. In scenario 3, as in the S-N direction, the demand compared to the base scenario, indicates a lower capacity of the network to support the flow.
Demand of public transport vehicles
Figure 10 shows the demand of public transport vehicles at the intersection approaches, for the three different scenarios. We report that the demand behaves totally different compared with the demand of total vehicles and we observe homogeneity between scenarios. Scenario 3 shows a behavior contrary to that of vehicles, in general, since the demand is greater in the S-N direction on 39th Street than in the N-S direction on 45th Street.
Queue length
Queue analysis allow us to evaluate if applied countermeasures are generating an increase in traffic congestion or an improvement in operating conditions. For this, queue lengths were measured at the entrance and exit of the study area. Maximum and mean queue lengths were analyzed to have a deeper understanding of the incidence of countermeasures.
Figure 11 shows that there is no relevant variation in the maximum queue length for the S-N direction, for any of the scenarios proposed. However, we report important increases in mean queue lengths if the countermeasures were applied, implying travel time increases, more congestion and worsening of level of service. N-S direction of the corridor shows similar behavior except for a great increase in mean and maximum queue length for the left-side segregated bus lane scenario, as shown in Figure 12.
Traffic conflicts
To validate the model, in terms of road safety, it is necessary to identify and describe the conflicts observed in the study area. The conflicts measurements were made from video data by a qualified team who identified and characterized the conflicts using STCT. This team of 8 undergraduate and graduate students in Civil Engineering received direct training from the author of the STCT, Christer Hydén, for the identification and characterization of traffic conflict, as well as in microsimulation of traffic conflicts in VISSIM and SSAM. Figure 13 shows the 18 lateral conflicts, that were identified in the field and classified as serious events. These conflicts will be compared with simulated conflicts.
Table 4 shows the number of conflicts extracted from the microsimulation models after ten different random runs of each scenario. In each scenario one simulation was selected as the most representative (highlighted in Table 4]. The selected simulation in the base scenario was the one with the minor MAPE, implying a minor difference between observed and simulated conflicts. The analysis and comparison between observed and simulated conflicts are made based on lateral conflicts since these were the ones that we observed in the field. It is worth noting that in all proposed scenarios there is a clear decrease in the number of lateral conflicts, with the first scenario being the safest. The mean-square error (MSE) for the mean TTC, between the 10 runs, was calculated to identify the most representative run of the 3 scenarios, also highlighted in Table 4.
Total: Total number of conflicts, L: lateral conflicts
Surrogate safety measures
The relationship between simulated surrogate measures, and the implications they have on the number and severity of traffic conflicts, is analyzed by means of linear regressions TTC vs PET, TTC vs MaxS, TTC vs DR and TTC vs Conflict angle, for the most representative run. This run, in scenarios 1, 2 and 3, was the one with the lowest mean-square-error, which were the third, ninth and ninth respectively as shown in Table 4. Also, the relationship between each of the variables and the probability of occurrence or severity of the conflict was taken into account [31].
Simulated conflicts obtained using SSAM, for the three countermeasures proposed, were obtained based on the initial conditions established in Table 5 from the base scenario. In the same way, the characteristics of the field-observed conflicts are used to improve the thresholds of each measure, also improving the simulation of conflicts.
Previous research, [32-37] have characterized conflicts through surrogate safety measures denoting a level of risk, severity or speed. Table 6 shows the risk thresholds for five surrogate measures analyzed.
Figure 14a shows linear regressions between TTC vs PET, the most commonly studied variables when working with SSAM. For the baseline scenario, as well as for alternative scenarios 2 and 3, a linear regression analysis can be done, given that data was obtained for more than two conflicts. However, for scenario 1, it is not possible to generate a trend by analyzing only data from two conflicts.
In [37], the relationship between “headway” and TTC in car following situation is exposed, that was later adapted as a relationship between PET and TTC in the same operational conditions in [31]. The latter generated the classification shown in Table 7. Although no classification was found for lateral conflicts, this is a good approximation to understand the possible events that can occur based on TTC and PET values and their relationship with the risk.
As seen in Figure 14a, there are no conflict data in the upper left part of the graph. This area corresponds to high PET values and low TTC values, which in practice are difficult to find. This is because conceptually [23] PET is the time that elapses from the encroachment of the trajectory until the invaded vehicle reaches the conflict point, and TTC is the theoretical value of the time to an eventual collision under given conditions of speed and distance. Then, if the time after the encroachment is high, it is difficult that along this there is a critical value under a theoretical TTC. It can be observed that, in general, the application of countermeasures results in the existence of conflicts that have higher values of TTC and PET, implying a lower probability of occurrence. In Figure 14a to Figure 14d, linear models were used to verify if the trend presented by the surrogate measures of the simulated conflicts corresponds to the expected trend according to the theory.
In the case of the evaluations carried out in the model for TTC vs MaxS, shown in Figure 14b, there is no high correlation between the variables with values of R2 that do not reach 0.5. However, in the base scenario and scenario 1, a slight negative trend is shown in the slope, while in scenario 3, a positive trend is observed. Associating MaxS with the severity of the conflict, we report that the three proposed operational configurations considerably reduce the number of conflicts with MaxS greater than 5 [m/s], which implies less possible severity on them. However, among the 3 scenarios proposed, no great difference is observed.
In the case of TTC vs DR, only for the case of scenario 2, a certain correlation is observed with a slight negative slope in the regression line, as shown in Figure 14c. The deceleration parameter increases in magnitude for the proposed countermeasures, which would imply a higher probability of occurrence of a conflict. This may be since, in a more orderly flow, due to segregation or the strict application of the norm, higher speeds may be presented at the moment of trajectory encroachment, implying greater decelerations. This analysis contrasts with what was found for the variable MaxS, in which a decrease in the maximum speed was reported in the conflicts for the proposed scenarios.
No correlation was found between TTC vs. Conflict angle [Figure 14d). Although, no literature has been found regarding this analysis, one can detect a large decrease in the angle in which conflicts occur. This is related to the decrease in lane changes of the vehicular interaction that is expected to occur in all the proposed scenarios.
Scenario analysis
According to Table 8, the first proposed scenario would not have improvement over traffic operating conditions. However, in terms of road safety one can expect great benefits, since the number of conflicts decreases substantially from 17 to 2. As for the severity of those conflicts, it is observed that it is similar, except in terms of deceleration, which increases, and in terms of speed, which decreases. For an improvement of the first scenario, stronger regulations are required, defining stop time of private vehicles and taxis on the preferential lane. This could also improve the operation of the preferential lane and avoid conflicts as a result of the transit of buses to the second lane. Additionally, the need to control on-street parking is evident, since, as observed in the field, this regulation is not presented.
Note: HR: High Risk, LR: Low Risk, MR: Moderate Risk, HS: High Severity, LS: Low Severity, HSC: High Speed Conflict, LSC: Low Speed Conflict (Variable characterization based on Table 6].
For the second scenario, which considers the physical segregation of the external preferential lane, except at bus stops to allow overtaking, the number of conflicts decreases from 17 to 4. Latent demand and public transport decreases in all senses with respect to the base scenario, except for the latent demand from north to south that increases. For this scenario, queues and deceleration increased in both directions compared to the base scenario, but speeds decrease.
The third proposed scenario has a lower capacity to support the flow because it presents a higher latent demand, which causes queue lengths to increase considerably. Of the three proposed scenarios, this is the one that reports the greatest number of conflicts when compared to the base scenario (6 conflicts). For this case, the deceleration increases, and the speed decreases when compared to scenario 0. Considering that physical segregation is proposed for the second and third scenarios, we expect higher interaction between vehicles, because of the decrease in capacity of mixed lanes.
6. Conclusions
We propose and validate a methodology to link field studies of traffic conflicts and the microsimulation of such events, showing that it is also possible to validate a simulation modeling based on traffic and safety variables. Under the calibrated environment, three scenarios of operational changes were simulated seeking to improve safety indicators. After this procedure, different parameters were established, and we propose them to be used as a starting point in future local road safety studies based on conflicts.
The implementation of the preferential lane of 7th Avenue has a positive influence on road safety conditions along this corridor, reducing conflicts and potential accidents. However, only when the preferential lane operates according to the existing regulations (scenario 1) improvements are observed. We did not observe such behavior in the field. There were 18 conflicts observed in the field and only 2 can be observed when simulating the operation established in the regulations, corresponding to scenario
The configuration of alternative scenarios of exclusive and/or preferential bus lanes also improve road safety, in terms of traffic conflicts and surrogate safety measures. The above is concluded from the number of conflicts obtained in scenario 2 (4) and scenario 3 (6). The mean values of the surrogate measures studied indicate a safer general operation in scenario 1, validating the regulation. This is caused by the reduction on the number of encroaching vehicles of the bus lane; therefore, fewer events between public transport and private vehicles are expected. Additionally, an adjustment is recommended in the normative measures, including regulating the stopping time and stopping zones of private vehicles and taxis letting for further research the study of acceptable stop times.
For the initial calibration process, there were some parameters avoided due to high uncertainty in data collection. Queue lengths are considered volatile parameters [28] and also not appropriate for this study because of the difficulty in counting the number of vehicles in queue or measuring its length when using video data. Those difficulties were mainly due to high distance between the data collection point and the intersections observed. However, the simulated queue lengths were used to evaluate the operational effect of the proposed countermeasures. In this sense, it was evidenced that the use of the bus lane under the regulations does not represent a deterioration of the operation. Travel times were also discarded owing to the difficulties to identify the vehicles circulating along the corridor as a consequence of the camera quality and the distance. The hypothesis of high uncertainty in these variables were confirmed with statistical analysis such as standard deviation which shows a high variability between the data collected. For this reason, these variables were not considered from the beginning of the project. The variables used represent the reliability and certainty of the model. The analyses were carried out with the most representative simulation run, which was the one where the greatest number of conflicts is presented from lane 1 to lane 2.
Some other limitations were found in the research process. Obtaining traffic conflicts information from video data was done using the STCT. In this case it was possible to measure values of Conflicting Speed (CS) and distance to obtain Time to Accident (TA). However, the measures obtained from SSAM were minimum TTC, PET, DR, among others. Software SSAM is not able to measure TA, because it cannot detect where an evasive manoeuvre has initiated. So, second step calibration was performed comparing field-measured TA values with TTCmin calculated by SSAM.