1. Introduction
Global climate change has induced rainfall extremes in terms of frequency, duration, pattern, and intensity in many cities of developing countries [
1]. As a result, urban flooding is one of the major challenges due to intensified rainfall events and inadequate urban engineering infrastructures [
2,
3]. For example, in the Bangkok metropolitan region (BMR), average total annual rainfall measured by the Thai Meteorological Department exceeded 2000 mm and caused road floods. In the BMR, maximal road flood depths are recorded by detectors installed to monitor flood conditions on roads and tunnels, and collect cumulative average rainfall measurements. When average rainfall accumulated at less than 60 mm/h, 14 flooded areas occurred, while rainfall at over 60 mm/h produced 56 flooded areas, which negatively impacted everyday activities [
4]. Drivers modify their travel patterns by taking alternative routes, leaving for destinations at different times or canceling the trip altogether. To better understand and minimize the traffic impact of different flood conditions, flood-related data must be integrated with road usage statistics.
The development and implementation of strategies that minimize the effects of adverse flood conditions require comprehensive knowledge of how flood events impact traffic operations, and how best to assess flood-related effects for a given scenario. Previous research integrated flood simulation with traffic analysis to model different viewpoints of the interplay between these two spatiotemporal phenomena [
5]. This integration was used to build a road-network model to assign trip paths, consider the effects of road closures, and evaluate travel delays and vehicle volume redistribution in a given flash-flood scenario by considering the influence of floods on mobility, mainly focusing on severe flood conditions. The relationship between flood events and traffic operations was assessed by a traffic simulation model that could handle complex roadway geometrics, traffic control devices, and vehicle configurations.
However, if key parameters within a traffic model could be changed under various flood conditions, this would greatly assist in developing improved flood-responsive traffic management strategies. The MFD or the relationship between traffic density and traffic variables [
6] was determined to assess the service level of the road network and serve as a vehicular access control measure in urban areas with complex road networks [
7,
8]
MFD determinants as the dynamic nature of traffic flow in urban networks were copiously studied [
9,
10,
11,
12,
13]. Suwanno et al. [
14] used precipitation intensity and the number of nonrainy days with Bluetooth-enabled traffic data installed at junctions to analyze the main MFD parameters decreasing in each rain condition. These parameters were applied to the delay function of the static traffic model to measure traffic efficiency. Results showed that increased traffic concentration had significant impact on the flow characteristics of the static traffic model of an urban road network. However, the intensity of the rain affected the shape of the MFD as a weak correlation factor because rainfall alone does not consider fast flow paths and areas where deep puddles completely block the roads. Moreover, using Bluetooth data collection in previous work showed a lack of transportation mode in upstream and downstream traffic, which is able to cut the outlier from only abnormal speed. Probe data are utilized in this study to eliminate this imprecision of the data of mode of transportation and driving speed.
Using the MFD concept, the volume–density relationship of a large network was used to generate a performance diagram and evaluate traffic-flow characteristics. As mentioned by Knoop and Hoogendoorn [
15], the Greenshields relationship was used to estimate MFD parameters for road surfaces under varying flood depths. A road-network model was then built to determine optimal travel routes by considering how the effects of road degradation impacted travel delays and vehicle distribution in a given road flood.
An accurate assessment of the impact of different flood conditions on MFD shape has not been previously undertaken for large road networks. The development of empirical relationships between depth of floodwater and vehicle speed reduction cannot analyze impact at the network macroscale that requires the relevant parameters of speed, flow, and density. To resolve this problem, changes in MFD parameters with flood conditions were assessed, and results were applied to generate a traffic model to analyze impact on the whole network. However, this study only considered MFD shape comparisons for a single reservoir. Future investigations should examine flow exchanges in multireservoirs to assess the principle of MFD-based modeling over multiple regions. The main idea is to protect a given area by limiting inflow at its periphery to maintain maximal throughput. In practice, inflow limitation can be achieved by tuning traffic-light timings or road pricing.
2. Data Collection and Extraction
The case study was for the Sukhumvit district, located in the eastern suburban area of Bangkok. This area has high flood risk because it is flat, and the gravity discharge of stormwater takes time. Pumps are used to move water to the main drainage system [
16,
17].To construct the MFD, taxi data that were collected during 2019 were used as open data (historical raw vehicles and mobile probe data in Thailand) provided by [
18]. Data from 5 a.m. to 12 a.m. on weekdays were included in the study with no freeway transport. Only taxis probes occupied by passengers were considered because they reflected the behavior of regular vehicles. In total, 78,078 taxis were recorded as registered in Bangkok in March 2019 [
19], and a taxi penetration rate per day of 22.45% was used to estimate the MFDs [
20]. The road flood monitoring system monitors flood conditions on roads and underpasses in the Bangkok area [
19]. The system collects data at 5-minute intervals. The MS SQL database stores the data for a long time period, making it a very useful tool for road-flood forecasting. A total of 94 sensors were installed on important roads throughout urban Bangkok, with 8 sensors in tunnels to measure and display real-time flood graphs.
This study estimated the MFD using vehicle trajectories from taxi probes that covered the Bangkok area. A taxi characteristically drives randomly with wide distribution, and the destination is determined by the passengers. This conforms to the characteristics of resident travel and reports their locations as timestamps at every uplink interval. Taxi trajectories were reconstructed on a time–space diagram as shown in
Figure 1.
For vehicle trajectories, formulae from Edie’s definitions [
21] based on time–space diagram presented in
Figure 1 were used to extract average flow (
) (Veh/h), average density
(Veh/km), and average speed
(km/h) in the network to obtain the MFD for a rectangular region in time and space with dimensions of length of time intervals (
T) and total length of all roads in the network (
L). Flow and density were defined on the basis of total distance travelled (
) and total time spent (
), as shown below for section
i.
where,
: distance traveled by vehicle k in section I, and
: time spent by vehicle k in section i.
Then, flow and density are defined as
where,
: flow by a vehicle in section i,
: density by a vehicle in section i,
: number of lanes in section i, and
: length of section i.
The flow and density of full traffic are estimated as
where
Lastly, the area average flow (
Q) and density (
K) are calculated by averaging section variables across an area according to the following equations.
4. Development of the Traffic Model
A dynamic traffic assignment (DTA) model implements a dynamic user equilibrium in the time-dependent nature of demand and network characteristics to produce more realistic traffic conditions on a large scale. Mesoscopic models, on the other hand, capture traffic dynamics in less detail, but are faster and easier to apply and calibrate than microscopic models are. Therefore, they are most suitable for modelling large networks, while microscopic models are usually applied to smaller areas [
22].
This study considered a mesoscopic traffic simulation as a dynamic equilibrium-based model to determine a balance between large-scale low-level details in macromodels, and small-scale high-fidelity in micromodels [
23]. This model can capture traffic-queue build-up and dissipation across lanes through upstream intersections caused by traffic controls (signals or signs), flow capacity, and lane-changing behavior [
23,
24,
25,
26] to generate and track individual vehicles or packets of vehicles. The model assesses the movement of these vehicles according to macroscopic relationships as speed, density, and flow, and considers queuing and spillback due to subject link capacity and downstream link queuing capacity.
Vehicles are simulated at the microscopic level to determine path costs. In these models, the planner specifies the level of detail for vehicles, time, and network inputs. The mesoscopic model consists of the three main elements of traffic demand, network definition, and traffic controls. Traffic demand is represented by a time-dependent origin–destination (OD) matrix for each class of modeled vehicles. The representation of the road network is at the level of individual lanes and includes the definition of turn pockets. The permitted lanes for each intersection movement and lane prohibitions by vehicle class are specified where necessary. The network definition does not require lane width or detailed intersection geometry, such as a turning radius. Traffic controls contain all relevant information regarding the different signaling plans on the intersections.
However, the existing Bangkok traffic model, developed in 2017 and called the extended Bangkok urban model (eBUM), is a macroscopic traffic model, developed by the Office of Transport, Traffic Policy, and Planning [
27]. Before applying the eBUM model to analyze and forecast transport and traffic situations in the mesoscopic traffic model, we updated the traffic model in 2019 with the most relevant current traffic conditions by using a four-step method. Then, we extracted the road network for the case study of the Sukhumvit area from the updated eBUM model and imported it into a baseline scenario as a new database. The matrix resulting from the transversal assignment was also imported into the new database and used as the demand assigned in the subarea scenario to perform the requisite adjustments. The subarea model consisted of 88 zones, 207 nodes, 766 links, and 16 traffic signals, as shown in
Figure 3.
This subarea network was the base scenario for the dynamic model debugging and calibration process. The models simulated the peak hour at full demand, with half an hour for the network to warm up and another half an hour to cool down after peak, with 80% of the peak-hour demand. Vehicles in the model were of the private car class. The model was run with a maximum of 200 iterations and a 1% relative gap. A 15-min assignment interval, equivalent to 8 intervals for each assigned OD matrix, was used. Up to 20 paths were searched in one assignment interval. Simulation results were collected every 15 min after the warm-up period for one hour.
In a DTA model, a relative gap of less than 2% is considered to be good, but is not easily achieved in highly congested networks (Chiu et al., 2010). However, the relative gap only gives an idea of the stability of the current solution and does not quantify the solution with respect to the underlying equilibrium (user-optimal) objective; this indicates that the algorithm is no longer improving the results, but does not indicate how well the final solution satisfies the desired objective of equilibrating travel times on alternative paths. Calibration is required to ensure that the model can replicate traffic volumes, travel times, and other travel patterns in the network. This process involves the calibration of both the demand and supply parts of the model. This was achieved by using model volumes and observed counts at the important road segments in the network.
Figure 4 shows a scatterplot graph that compares model flow and observed flow for the morning peak hours between 7:30 and 8:30. Results indicated that the model was well-calibrated, with
of 0.975.
Network simulation results are based on link results aggregated to totals or averages for the entire network, as shown in
Table 3, ‘In count’, we measured the number of vehicles that entered the physical space of the network and passed the entrance connectors during the time interval. ‘Out count’ measured the number of vehicles that exited the physical space of the network during the time interval and exited from the exit connectors, and VHT measured the total vehicle hours of travel during the time interval inside the network. This measure did not include waiting time on virtual links, while VKT measured total vehicle kilometers traveled during the interval. As with VHT, this does not include vehicles on virtual links.
Figure 5 shows the relevant outputs as the number of vehicles waiting to enter the network and the number of vehicles remaining in the network at the end of the simulation. The blue line changed suddenly at the third interval (7:45), reflecting the start of the peak-hour interval until after the last peak interval demand. Over 25,000 vehicles were waiting to enter the road network. These vehicles were stuck on virtual links and could not enter the road network to complete their trips.
Application of MFD Parameters to the Traffic Model
THe systematic consideration of the influence of patterns in a supply and demand structure on a global scale is a requisite for the successful and effective transfer of policy interventions between networks with different structures. Thus, the influence of supply and demand on traffic conditions was investigated. Demand was assessed as the magnitude and distribution of movements of travelers between different locations when aggregated across the geographical area served by a network. The main study interest concerned the performance of the transport system when placed under extreme pressure as an average weekday morning peak period for private cars and assumed that modal travel behavior did not change much during road floods. Supply was assessed as physical infrastructure, composed of roads and intersections, and the capacity to efficiently facilitate the movements of travelers by systematically changing input MFD parameters. Flood depths are displayed in
Table 4. Maximal flood depth was considered to be 30 cm.
A mesoscopic model was used with three components of traffic flow, as shown below [
25]. (1) Car-following models describe the steady-state properties of traffic on a roadway. Steady-state properties can be observed when traffic is moving at a constant speed, and thus also at a constant flow and density. (2) Gap acceptance is the amount of time available for a lower-priority vehicle to execute its movement and clear the point of conflict with a higher-priority vehicle. The relative wait is the difference between the delay that two vehicles incur at an intersection (not including delays due to red traffic signals) as the next vehicles to exit by their respective lanes. This delay is the amount of time spent waiting for acceptable gaps in higher-priority movements. (3) Lane changing is the movement of vehicles on individual lanes of a roadway. How drivers utilize the lanes on the road can significantly impact delays and how these delays propagate through the network. Vehicle trajectories within links are implicitly modeled. Each driver must choose the lane by which to enter and exit a link just before actually arriving at the link. Once onto the link, the choice cannot be reconsidered. The rules used to model drivers’ lane-choice behaviors are complex, combining a look-ahead procedure with local lane-choice rules. The look-ahead feature captures the behavior of drivers familiar with the roadway and recurrent congestion patterns along their usual paths.
Following the above, the application of the parameters in
Table 4 that define the supply part of the car-following model, which can vary, are response time and relative length, as shown in Equation (
5). The effective length (
L) of a vehicle as a predefined value was set to be 6.25 m, while the driver’s response time (
R) was set at 1.25 s [
25], both the microscopic traffic parameters of this model. Any modifications to these two parameters results in a change in the jam density and maximal flow rate of each link. Free-flow speed is the other important component that defines the capacity of the links, and the set values were modified according to flood levels.
where,
: position of the following vehicle at time t;
: position of the leading vehicle at time t;
: free speed of the roadway,
L: effective length of the following vehicle; and
R: response time of the following vehicle.
For the mesoscopic traffic model, a triangular diagram was chosen. This can be described by only three parameters, thereby minimizing the required input data. The triangular diagram is presented in
Figure 6 below.
Values of the three macroscopic traffic-flow parameters for a specific vehicle type and a specific roadway can be determined from the free-flow speed of the link, the effective length (
L), and response time (
R) of the vehicle type, as shown in Equation (
6) below.
where,
: maximal possible flow rate (veh/h/lane);
: jam density (veh/hm/lane); and
: speed (km/h) at which shock waves move through a platoon of traffic against the direction of flow, for a specific vehicle type.
VKT and VHT during the peak hour were selected as the two major performance measures. VKT is regarded to be a benchmark of network utilization, i.e., how much traffic can be served in the system. When demand increased, more vehicles traveled in the network and VKT increased; however, in the peak hour, this did not hold because VKT and VHT were determined by vehicles that had completed their trips. Due to network congestion, the number of exited vehicles (vehicles that had been served) stabilized when the system throughput reached capacity, while demand continued to increase.
Figure 7 shows the number of vehicles entering and exiting from the network. The number of vehicle entries and exits in flooding scenarios shows a decrease from the baseline scenario to the high-depth scenario. This occurs because the number of trips by entering vehicles reduces, while for exiting vehicles, congestion in the network was high, and vehicle movement decreased. Most of the vehicles could not reach their destination on time.
Figure 8 compares the base scenario (dry situation) and flood scenarios. The network indicator evaluation clearly explains that, when floods occur, VHT (
Figure 8) increases, and VKT (
Figure 8b) in the network decreases because the number of vehicles entering and exiting the network decreases.
5. Traffic Management during Flooding
The base dry scenario, and the highest flood of 20–30 cm scenario, both with and without traffic management measures, were simulated in the model. Results were compared to investigate the improvement of road-network performance after the implementation of traffic management measures for (1) an optimized signal system management, and (2) when traffic demand on arterial roads was induced onto the expressways. Road link performance in the study area was also evaluated.
5.1. Optimization of Signal System Management
The optimization of signal system management by synchronizing multiple signalized intersections generates a “green wave”, allowing for platooning vehicles to travel along the arterial road with minimal or no stops [
28]. Synchronizing provides incentives for vehicles on arterial roads to speed and catch up with the platoon, so that they can pass through the intersection without any stops [
29].
This study considered the mesoscopic model (Dynameq) that can model different intersection controls. Dynameq allows for signal offsets as inputs and provides improved computational efficiency compared to a microscopic model in a regional network. This is particularly important for simulation-based DTA. However, Dynameq has no program packet to cater for multiple coordinated intersections. Therefore, the microsimulation model (Vissim) was used to design a working layout for junction and traffic-signal design, and coordinate signal timings for Dynameq. The Dynameq model was used to evaluate how traffic was assigned to the local road network, considering delays as a result of flooding. Results of traffic flows and signal groups at each signalized node were input into the Vissim model to attain coordinated signal timings, which were then fed back into the Dynameq model to produce a final set of traffic flows for the network.
Figure 9 Vissim output using genetic algorithms for signal timing optimization for delay minimization [
30,
31]. Average delay of all vehicles passing through the signal heads at the intersection can be automatically determined from the node evaluation in Vissim. To optimize the signal stages, the highest delay times in all the stages were determined. Stages with minimal and maximal times were adopted as the best and worst stages. The worst stages’ green time was adjusted by deducting the time from the best stage. If the time could not be deducted from the best stage, the next best stage was selected for deduction. Similarly, the remaining stages were optimized for signal timing.
5.2. Use of Elevated Metropolitan Expressway with Free Tolls
The concept of an expressway with free tolls during flooding was suggested by [
32] as an application of the MFD congestion pricing strategy (generalized cost) to improve mobility and relieve congestion in cities by reducing delays and the length of the rush hour. The DTA solution refers to assignments being in approximate equilibrium, where each driver attempts to minimize travel time or generalized cost, including intersection delays (total cycle time) and link cost as an operating cost including fuel cost and driving cost (value of time) multiplied by the time in the network to arrive at an operating cost. We added a perception penalty to the generalized cost for each road type, while the toll penalty was only used in scenarios that included congestion pricing. The toll penalty was calculated by dividing the toll cost assigned to a link by the assumed value of time to achieve a time cost in seconds of paying a toll. To be consistent with the dynamic model, the current model assumed that the value of time (fuel+distance cost, FC) was 86 s per km [
27]. Trips were assigned to the network using the generalized cost (GC) Equation (
7).
where,
: generalized cost for movements;
: left-turn penalty (10 s);
: right-turn penalty (30 s);
: facility penalty (fuel and distance cost);
: toll penalty.
Generalized cost in the expressway for this scenario is the same as another road type that provides a simplified description of travelers’ choices the effect of flood in their decisions to improve the city’s mobility. If the free toll is implemented during floods, road users are likely to take an elevated road instead of a road with the possibility of flooding. Therefore, the implementation of the free toll improves network performance and ensures the unobstructed flow of vehicles during floods.
5.3. Results of Traffic Management
The simulation period to evaluate how vehicles entered and exited the network during floods covered only peak hours. The numbers of vehicles entering and exiting, as shown by the model, were consistent for all scenarios.
Dry-scenario results in
Table 5 were assessed with traffic management measures. The elevated metropolitan expressway with free tolls (1) showed a slight change in the number of vehicles entering (in count) and exiting (out count) the road network by 5.34 and 6.94, respectively, but almost no change in optimizing the signals (2). However, the total number of vehicles in the network at the end of the interval (traveling) showed a high increase in all three measures. VHT increased in optimizing the signals (2) by 9.92%, as the in count was allocated to travel at intersections to minimize delays in the overall road network.
Table 6 shows the results of traffic management measures during flooding with the improved effectiveness of the flood depth of 20–30 cm scenario. Using the elevated metropolitan expressway with free tolls (1) showed a slight increase in the number of vehicles entering and exiting the road network, thus slightly improving VHT as vehicles chose alternative routes. By optimizing the signals (2), the number of vehicles exiting the network increased by 22% and improved VHT by a 15% reduction.
Overall, network efficiency in the dry scenario improved when traffic measures were taken, and there was only an increased VHT for optimization of the signal measures. In the flood scenario, the optimization of the signal condition showed better results than those of the use of the elevated metropolitan expressway with free tolls. Signal optimization impacted the overall network, whereas the use of the elevated metropolitan expressway with free tolls applied only to limited roads (expressway). When the elevated metropolitan expressway with free tolls and the optimization of the signal network were applied together (1 + 2), the number of vehicles exiting the network or completing trips increased by 38%, VHT decreased by 27%, and VKT increased by 30%. Overall, a substantial improvement in the network was achieved.
Results showed that VKT after applying traffic management measures by all vehicles showed improved performance of the whole network. VHT was reduced with a higher vehicle out count. The decrease in VHT indicated a better traffic situation. Similar to VHT, VKT is a meaningless criterion to reflect the change in traffic demand. Nevertheless, when combined with VHT, VKT indicated the extent to which drivers changed their route choices, and this affected network performance. Sometimes, other criteria such as average speed and delay are considered. However, evaluation results are mostly represented by travel time and travel distance.