Assessment of Spatiotemporal Peak Shift of Intra-Urban Transportation Taking a Case in Bangkok, Thailand

: Reducing congestion has been one of the critical targets of transportation policies, par-ticularly in cities in developing countries suffering severe and chronic trafﬁc congestions. Several traditional measures have been in place but seem not very successful. This paper applies the agent-based transportation model MATSim for a transportation analysis in Bangkok to assess the impact of spatiotemporal transportation demand management measures. We collect required data for the simulation from various data sources and apply maximum likelihood estimation with the limited data available. We investigate two demand management scenarios, peak time shift, and decentralization. As a result, we found that these spatiotemporal peak shift measures are effective for road transport to alleviate congestion and reduce travel time. However, the effect of those measures on public transport is not uniform but depends on the users’ circumstances. On average, the simulated results indicate that those measures increase the average travel time and distance. These results suggest that demand management policies require considerations of more detailed conditions to improve usability. The study also conﬁrms that microsimulation can be a tool for transport demand management assessment in developing countries.


Introduction
The global passenger transportation demand was 51 trillion passenger-km in 2014, which is expected to more than double and reach 110 trillion passenger-km in 2060 [1]. Most of the increase is expected to emerge in developing countries where the population and economy are growing significantly. Many cities, both in industrialized and developing countries, have faced severe traffic congestion problems due to increased transportation demand. Reducing congestion has been one of the key targets of transportation policies for various purposes, including alleviating economic loss, reducing pollutants and energy consumption, and reducing travel time. In particular, cities in developing countries suffer severe and chronic traffic congestions caused by the surging travel demand and insufficient transportation infrastructure [2]. Many studies have demonstrated that the road infrastructure cannot be a single solution to relieve congestion [3][4][5][6]. Other studies suggest that a combination of various transportation demand management measures, including public transport provision [7], pricing [8], information service [9], coordination with land use [10,11], and active transport promotion [12,13], are needed to tackle congestion. These measures usually require public policies. Recently, the use of smart technologies in the transport sector has also been investigated intensively, such as bike-sharing [12,13], application of big data [14][15][16], implementation of mobility as a service platform and information services [17][18][19][20][21], and new urban transit system [22]. Bangkok, Thailand, is the target city in this study. It is known as one of the most severely traffic-congested megacities in the world [23]. Various countermeasures to alleviate the congestion have been implemented, including expansion of road capacity, improvement in public transport [24], but its congestion remains at a very high level. These measures were mainly related to infrastructure development and improvement, while the demand management measures were not successfully implemented. There are various studies that analyze important aspects of this situation, including automobile dependencies [25], a spatial demand-supply gap of public transportation [24], assessment of accessibility to the metro system [26], investigation of impact factors on transportation mode choice [27], obstacles towards sustainable transportation [23], and impact of telecommuting on traffic [28].
Most of these studies employed statistical approaches or static transportation models rooted in the four-step transportation model design [29]. After 2000, agent-based modeling approaches [30] have been applied in transportation research to answer various research questions considering the individuality of travelers and temporal transportation behavior. This highly disaggregated approach provides some insight into transportation research, spatial development, quality of life, and social equity, becoming more prominent in urban policies [31]. Nevertheless, this approach usually requires massively disaggregated activity data [32]. Therefore, applying agent-based models to developing countries is a challenge in terms of data limitation.
The purpose of this study is to apply the agent-based transportation model MAT-Sim [33] for transportation analyses in Bangkok to assess the impact of spatiotemporal transportation demand management measures, including working place dispersion and departure time shift. We conducted this research as a part of a Japan-Thailand bilateral international research project, and transportation and planning authorities in Thailand provided various transportation and spatial data. We also utilized unconventional data, including data compiled by private companies. We took advantage of existing statistics and data for the analysis and did not collect new data in the field. However, available data were not sufficient to generate the input for MATSim. To complement the data, we proposed a method to estimate a model to generate complementary data. With these generated data, we analyzed the impact of transportation policy measures on travel demand. Hereafter, Section 2 shows the features of data used in this study, and the needs for complementary data are discussed. Section 3 describes a method to generate the complementary data for a microsimulation. Section 4 shows the spatiotemporal demand shift policy scenario. Section 5 describes the analytical results. Section 6 is the discussion and conclusion.

Data
In general, official statistics for transportation in developing countries are insufficient for microsimulation transportation models. In transportation analysis, trip distribution data are essential information to capture the spatial travel pattern; however, it usually requires many samples in the survey. In addition, road network data are indispensable for traffic assignment. Recently, open data, such as OpenStreetMap, are available for transportation analysis in many cities. For the public transportation analysis, common format data, such as General Transit Feed Specification (GTFS), are used in many cities, and it is being utilized as input data for transportation research.
The target area in this study was Bangkok and its five surrounding provinces, namely, Nonthaburi, Samutprakarn, Pathum Thani, Samut Sakhon, and Nakhon Pathom. Figure 1 shows the target area. The total population was 15 million in 2017, and the total area is 7735 km 2 . Table 1 summarizes the data used in this study.  First, the essential data to capture the transportation behavior was the Household Travel Survey (HTS) provided by the Office of Transport and Traffic Policy and Planning (OTP). These data record daily personal travel behaviors, including the origin and destination zones, departure and arrival time, transportation mode with attributes of respondents, including age, gender, and household information. These data contain essential information for the agent-based simulation model. However, the number of samples was 38 thousand, which is only 0.2% of the population in the target area. These data were not sufficient to capture the metropolitan-wide travel pattern. Past study indicates that the trip generation and trip attraction volumes correlate with the population at the origin and economic activities at the destination, respectively [34]. In this study, we estimated the zonal trip attraction volume based on the number of employees at the destination zone and estimated the destination choice model combined with HTS data. Here we utilized the trip data for all purposes, not limited to commuting trips. We assumed that the number of activities other than work was also spatially correlated with the number of employees. Therefore, the trip concentration for all purposes except go-home was also estimated based on the number of employees.
Here the population by age classes and gender in Table 1 are given only for subdistrict level. For the Traffic analysis zone (TAZ), which are smaller zones than the subdistricts, the population was provided for total number only. This study estimated the population by attribute for TAZ as proportional to the belonging subdistrict population by attributes. We utilized MapFan DB Data Model to construct a road network model. A private company provided these data for the base data of car navigation systems. It included the road type and the number of lanes for each link. We obtained the GTFS data from OTP for the public transportation network and its schedule. The public transport modes included railway, metro, tram, bus, and ferry. Figure 2 shows the population and employment density of TAZ and the road and public transportation network of the target area. We utilized MapFan DB Data Model to construct a road network model. A private company provided these data for the base data of car navigation systems. It included the road type and the number of lanes for each link. We obtained the GTFS data from OTP for the public transportation network and its schedule. The public transport modes included railway, metro, tram, bus, and ferry. Figure 2 shows the population and employment density of TAZ and the road and public transportation network of the target area.  Figure 3 shows the analytical flow of this study. First, we estimated the travel time matrix (skim matrix) using MATSim. For the estimation, tentative origin-destination (OD) travel demand was generated based on the zonal population and employment. The homebased trip generation volume was assumed to equal to the population, and the volume of trip attraction was assumed to be proportional to the employment at the destination zone. We estimated a gravity model parameter so that the average Euclidean distance equaled the average trip length of the HTS data. Using this gravity model, we estimated the initial OD pattern and modified it to make it consistent with the assumed trip generation and attraction volumes using the Fratar method. Sampling from this OD pattern and departure time distribution in HTS data, we generated the initial travel plans of agents for MATSim. As a result, we calculated the skim matrices for car and public transport based on the simulated travel time, which represented the daily average of travel time between  Figure 3 shows the analytical flow of this study. First, we estimated the travel time matrix (skim matrix) using MATSim. For the estimation, tentative origin-destination (OD) travel demand was generated based on the zonal population and employment. The homebased trip generation volume was assumed to equal to the population, and the volume of trip attraction was assumed to be proportional to the employment at the destination zone. We estimated a gravity model parameter so that the average Euclidean distance equaled the average trip length of the HTS data. Using this gravity model, we estimated the initial OD pattern and modified it to make it consistent with the assumed trip generation and attraction volumes using the Fratar method. Sampling from this OD pattern and departure time distribution in HTS data, we generated the initial travel plans of agents for MATSim. As a result, we calculated the skim matrices for car and public transport based on the simulated travel time, which represented the daily average of travel time between zones. The skim matrix for non-motorized mode was obtained by dividing Euclidean distance between zones by the assumed walking speed of 5 km/h. Second, we estimated the logit model parameters for mode choice by age class and gender using HTS data and the skim matrices. The explanatory variable of this model was travel time by mode only. We calculated the expected minimum time between OD using estimated parameters of the mode choice model, and this representative time between zones was used as the explanatory variable for the destination choice model.

Method
Here, for the walking travel, the Manhattan distance can be a better approximation of the travel distance than the Euclidean distance in practice. Obviously, the Manhattan distance is larger than the Euclidean distance when the distance between origin and destination is long, and therefore it could overestimate the non-motorized travel share for short trips. Even though the model parameters were estimated under the assumed distance to represent the observation, so the bias caused by the assumed distance was expected to be corrected to some extent. Table 2 shows the estimated parameters of the mode choice model where the parameters for all the segments met the sign condition except for females less than ten years old and males more than or equal to 70. The p-values were almost zero except for the age class less than 10, more than or equal to 70, and males of 50-59. So, most of the models were statistically significant, although McFadden's R2 of all models was quite low, i.e., the representability of the mode choice model was poor. This is because the transportation mode choice is actually affected not only by travel time but also by various factors, such as availability of car, possibility to accompanied drive with a family member, fare, reliability on timetable, safety for public transport, etc. The mode choice decision is, therefore, affected by each individual's conditions. In this study, most of the mode choice model parameters satisfied the sign condition and were statistically significant. Thus, we assumed that the average behavior could be reflected by the model parameters, which were applied for the transportation mode choice estimation. The other studies applied alternative approaches to estimating the modal share under the data limitation, for example, modeling based on the expert survey [35] and statistical model using travel time ratio [36]. They can be candidates for the modal split analysis under the data limitation. However, we used the discrete choice model that was highly compatible with agent-based modeling. Third, we estimated the destination choice model. Here the trip pattern of HTS survey and the estimation based on the gravity model were largely different due to the small sample of HTS. On the other hand, the HTS survey contained sample attributes information. To reflect both the information of HTS and crude estimate of OD pattern, we estimated the destination choice model parameters θ k , and ρ k that maximized the following log-likelihood function LL.
where q d is the number of trips of distance class d, q dk is the estimated number of trips of distance class d of traveler's attribute k by model, T ijk is expected minimum travel time between zones i and j of attribute k, N Ej is the number of the employee at zone j, N Rik is the number of residents of attribute k at zone i, Ω d is a set of OD pair belonging to distance class d, n z is the number of zones, Q ij is travel demand between i and j given by gravity model and Fratar method, σ 2 is the variance of OD demand model. By definition, LL1 is the log-likelihood of travel distance distribution, and LL2 is that of OD travel demand of the model. Table 3 indicates the estimates of the destination choice model parameters. The parameters indicate that the choice probability of the destination zone was higher for shorter the access time to the zone and larger numbers of the employee at the zone. Most of the parameters were statistically significant, while some parameters did not have an adequate estimation of variance; therefore, we could not determine the significance of those parameters. As formulated in the likelihood function, the parameters for all segments were estimated simultaneously. There might be correlations between the explanatory variables among segments that would cause the failure to estimate the variance of some parameters. Merging the age class may improve the estimation stability. The destination choice model was statistically significant according to the likelihood ratio.

Scenarios for Transport Policy Measures
Using the models estimated in the previous section, we estimated the OD travel demand by transport modes and generated agents' initial daily travel behavior. We generated the coordinates of origin and destination randomly within the geographical extent of the origin and destination zones. The agents were assumed to make two trips a day. The departure time of the two trips was sampled from the HTS. Here, we input these initial travel behaviors to MATSim by mode separately, and the mode choice module of MATSim was not applied afterward. In other words, the modal share was fixed to the initial estimation by the mode choice model and did not feedback the MATSim simulation results to the mode choice and destination choice behavior. As a result, each agent's travel time and distance were simulated for each demand management scenario.
In this study, the decentralization of employment places and the time shift of peak hours were analyzed as the demand management scenarios. Compared with the business as usual (BAU) case representing the current employment and departure time, we estimated the impact of these scenarios. The time shift impacts combined with mobility as a service (MaaS) system in Bangkok have already been analyzed in the other article [37], but we analyzed the impact of combination with the decentralization. To focus on the impact analysis, we fixed the other factors, such as population, road infrastructure, and public transport level of service, to the current level.
As a decentralization scenario, 30% of employees of the zones ranked as top 10% in terms of density (i.e., these were the CBD where the employment density is very high) moved into the zones ranked as top 10-20% (i.e., these are potential urban subcenters in the less-dense area that are being connected or soon to be connected by urban railway lines). Such employment decentralization will change the destination choice probability and relieve the traffic congestion in the city center. That will affect the trip length, travel time, and road speed. Figure 4 shows the location of the zone with changed comparing with the BAU case: decreased zones displayed in blue and increased zones in reds.
As a time shift scenario, we arbitrarily shifted the departure time of 30% of agents at peak time. The morning peak hour was 7:00-9:00 and the evening peak hour was 16:00-18:00. In this scenario, the sampled 30% of agents departing at the morning peak hour were forced to depart 2 h later, and the sampled 30% departing at the evening peak hour was made to depart 2 h earlier. This scenario will affect the road congestion and may change the travel time and shift the route. It was also expected to affect the public transport travel time and route because these factors were different time by time depending on the time schedule of public transport. Figure 5 shows the departure time distributions of BAU and time-shift scenarios.
terms of density (i.e., these were the CBD where the employment density is very high) moved into the zones ranked as top 10--20% (i.e., these are potential urban subcenters in the less-dense area that are being connected or soon to be connected by urban railway lines). Such employment decentralization will change the destination choice probability and relieve the traffic congestion in the city center. That will affect the trip length, travel time, and road speed. Figure 4 shows the location of the zone with changed comparing with the BAU case: decreased zones displayed in blue and increased zones in reds. As a time shift scenario, we arbitrarily shifted the departure time of 30% of agents at peak time. The morning peak hour was 7:00-9:00 and the evening peak hour was 16:00-18:00. In this scenario, the sampled 30% of agents departing at the morning peak hour were forced to depart 2 h later, and the sampled 30% departing at the evening peak hour was made to depart 2 h earlier. This scenario will affect the road congestion and may change the travel time and shift the route. It was also expected to affect the public

Results
We applied these scenarios with MATSim and compared their simulated travel times and distances. The sample rate was chosen to be 5% for road transport and 1% for public transport in the simulation, which is a reasonable approach in agent-based modeling to economize on runtime [38]. The total volumes indicated in the results were scaled back to 100% to be consistent with observed demand volumes. In this section, first, we validate the model, and then we present the simulated results of the scenarios.

Validation
The simulation results were compared with the observed data. Table 4 shows the key travel demand indices of the HTS 2017 and the simulation results produced by MATSim. The HTS estimation of car travel demand 2qw 14.9 million trips/day, and 9.1 million trips/day for motorcycles. Our estimation did not distinguish cars from motorcycles, and the total demand was calculated to be 17.5 million trips/day. Assuming the passenger car equivalent value (PCE) of a motorcycle was 0.33, the travel demand of the HTS was equivalent to 17.9 million PCE-trips/day. The difference of this volume to our estimate was 2%. The number of trips by public transport in our simulation was smaller than that in the HTS. In MATSim, the public transport module simulated the agent's behavior to use public transport and walk. If the travel time of public transport was too long, the agents did

Results
We applied these scenarios with MATSim and compared their simulated travel times and distances. The sample rate was chosen to be 5% for road transport and 1% for public transport in the simulation, which is a reasonable approach in agent-based modeling to economize on runtime [38]. The total volumes indicated in the results were scaled back to 100% to be consistent with observed demand volumes. In this section, first, we validate the model, and then we present the simulated results of the scenarios.

Validation
The simulation results were compared with the observed data. Table 4 shows the key travel demand indices of the HTS 2017 and the simulation results produced by MATSim. The HTS estimation of car travel demand 2qw 14.9 million trips/day, and 9.1 million trips/day for motorcycles. Our estimation did not distinguish cars from motorcycles, and the total demand was calculated to be 17.5 million trips/day. Assuming the passenger car equivalent value (PCE) of a motorcycle was 0.33, the travel demand of the HTS was equivalent to 17.9 million PCE-trips/day. The difference of this volume to our estimate was 2%. The number of trips by public transport in our simulation was smaller than that in the HTS. In MATSim, the public transport module simulated the agent's behavior to use public transport and walk. If the travel time of public transport was too long, the agents did not use public transport but traveled on foot. In the simulation, a substantial number of agents were estimated to travel only by walking, reflecting a bias of level of service estimation due to the small sampling. The total and average trip length of the car in our estimation was much longer than those of HTS. The simulated trip length was almost double the average Euclidean distance between ODs. This trip length reflects the network topology and the detouring behavior to avoid traffic congestion in the simulation. The simulated total trip length of public transport was much shorter than that of HTS. The simulation underestimated both the number of trips and the average trip length of public transportation as well. The car travel time of our estimate was much longer than the HTS, reflecting the estimation of the longer trip length. The average travel speed was almost the same level as HTS because the model was calibrated to fit with the average speed. The estimated travel time of public transport was relatively close to the HTS data, but the estimated travel speed was slower by about 20% because of the shorter trip length in our estimation. While these differences deserve further investigation, it should be noted that HTS was a comparatively small survey with a certain level of uncertainty. The resulting traffic flows appear to be sufficient for scenario analysis. Figure 6 plots the observed and estimated hourly average link speed. From the observed data panel (a), the speed declined during the morning peak hours. During the daytime, the traffic speed was fairly low and stable at the same level as the morning peak hour. In the evening peak-hour, the average speed became lower and reached the bottom between 19:00 and 20:00. After that time, the speed recovered to free flow toward midnight. While our estimation shown in panel (b) indicates that the average speed was too low in the morning and evening peak hours, the midday speed recovered nearly to the free-flow speed. Such a difference between the observed situation and the model estimation is probably because we did not have business trips or logistics in the simulation, resulting in the total traffic volume being underestimated. However, when considering the acceptable morning and evening peak pattern, we decided to calibrate the MATSim to fit with the daily average speed.
Even though there were biases in traffic flow and speed estimation, we applied this model to the sensitivity analysis of the scenarios. The discussion section will address the issues related to estimation bias. Figure 6. Hourly average link speed: (a) observed data by iTIC (Intelligent Traffic Information Center Foundation); (b) estimation in this study; iTIC collects the data from probes on taxis, buses, logistic vehicles, and mobile phones. We compiled the data using the raw data from 15 to 24 November 2020 in the TAZ region and calculate hourly average link speed on weekdays. Panel (b) shows the average speed trajectory of simulated road travel aggregated by departure time. Figure 7 shows the rate of change in total travel time and distance for the scenarios of time shift, decentralization, and their combination compared with the BAU. The travel time of car was reduced for all scenarios, while that of public transport increased. Car travel time was reduced by 10.5% by the time shift scenario, reflecting the alleviation of congestion by avoiding the peak hour. The decentralization scenario reduced travel time by 2.4% and alleviated the congestion by shifting the traffic concentration from the center to the near suburbs. The combined scenario reduced travel time by a car by 13.4%, which was more than the sum of the two scenarios. That implies synergies in the spatiotemporal transportation demand management measures. On the other hand, the travel time of public transport slightly increased for all scenarios. All scenarios had the same timetable for public transport. In other words, the supplied level of service was unchanged. The results indicate that a simple departure time shift would cause an increase in waiting time or a detour depending on the public transport service level. In this case, the average travel time increased by 1.3% by the time shift scenario. The decentralization scenario increased the travel time by 1.1%. This scenario did not necessarily move the employment to places that were accessible by public transport. The combined scenario increased travel time by 1.4%. Spatiotemporal demand management measures did not necessarily have a uniform effect on public transportation travel.

Effect of Policy Scenarios
Travel distances declined in the time shift scenario by 1.0% for cars and inclined by 0.7% for public transport. In the simulation, car users could choose shorter routes when road congestion was reduced. In contrast, the travel distance of public transport increased probably due to the timetable and change in routes. The decentralization scenario decreased the travel distance of cars by 0.3% and increased that of public transport by 1.2%. Dispersion of employment to near suburbs may reduce the distance between home and workplace or avoid the detour by congestion for car users. Meanwhile, decentralization may affect the route choice of public transport or increase the egress distance. The combined scenario decreased the travel distance of a car by 1.4% and increased that of public transport by 0.8%. Figures 8 and 9 show the spatial distributions of the rate of change in travel time by departure zones for car and public transport, respectively. Panel (a) in Figure 8 indicates that the time shift scenario decreased the car travel time for most zones but increased for a few zones. The decentralization scenario shown in panel (b) had many zones where the travel time increased. This result implies that the spatial dispersion of the urban core function possibly causes longer travel for some departure zones. From this panel, we can find that zones in the central area, and inner suburbs tended to decrease the travel time.
On the other hand, outer suburbs had relatively more zones that increased the travel time. Thus, possibly decentralization will increase the congestion in the outer area even though it reduces the congestion in the central area. The combined scenario shown in panel (c) decreased the travel time for more zones, and the decreasing rates were larger than the time shift solely.  As shown in Figure 8, there was no clear pattern in the spatial distribution of travel time change for public transport. The zones where the travel time increases or decreases seem to be mixed geographically, and many zones had a larger rate of change than car travel. Smaller differences were likely to be caused by the random variation of different model runs. Public transport users may be affected more by a change in departure time or destination place than car users because the travel time and optimal route are sensitive to the timetable and transit route connections. This result indicates that the spatiotemporal measures of transportation demand management have to be coordinated carefully for public transport operations.

Discussion and Conclusions
In this study, we applied MATSim to the Bangkok metropolitan region in Thailand. Even with access to the data provided by government authorities through the governmentsupported project and complemented by private companies' data, the input information was insufficient for an agent-based model. Therefore, we developed an estimation method of travel patterns and synthesized the input data for MATSim. On the other hand, we used GTFS data for the entire Bangkok metropolitan area for public transportation analyses, which is still not applied in many studies in transportation research in developing countries.
Applying this model, we analyzed the impact of employment decentralization and peak time shift on transportation. These spatiotemporal peak shift measures were effective for road transport to alleviate congestion and reduce travel time. However, the effect of those measures on public transport was not uniform but different by users' circumstances. On average, the simulated results indicate that those measures increased the average travel time.
The effects of peak time shift measures on public transport were different depending on the timetable and routes. There was no clear spatial pattern of the travel time change. These results suggest that demand management policy measures for public transport require consideration of more detailed conditions to improve usability. For instance, the employment dispersion should be incorporated with transit-oriented development plans to improve the access and egress to and from public transport. The time shift measures may need to be tailored to reflect the individual travel circumstances to reduce the travel time and congestion for the public transport system. Potentially, MaaS systems may be an effective tool, especially for the demand management of public transport.
This study conducted the analysis for the base year population in 2015. In practical transportation planning, population forecasts are required. Some studies projected the depopulation of Asian cities in the latter half of this century [39]. We need to reflect on the future population and its spatial distribution in the simulation. The spatial distribution of the urban population is also known to be affected by the transport system [3,40]. For future studies, we need to incorporate this interaction between transport and land use. In the decentralization scenario, we arbitrarily dispersed the employment from the center to near suburbs. However, more precise scenarios are required to coordinate with transit and urban development policies.
The model output can be used not only for estimation of travel time and distance but also for the other various analysis, such as environmental load, energy efficiency and consumption [41], effective land use and land value, value capture for transportation investments, and social disparities. For instance, based on a statistical model, urban structure and road congestion reduction were demonstrated to have only a limited impact on CO 2 emission reduction in the transport sector [42] using aggregated analysis. Microsimulation models can analyze more detailed spatiotemporally policy measures that cannot be analyzed by conventional models. Such microsimulation models may provide in-depth knowledge in policy practice. In particular, agent-based models are known to be superior in pricing studies and equity analyses.
Another study integrated MATSim with an agent-based land-use model [43] to analyze whether SDGs may be achieved by urban and transport policies [44]. Our study will incorporate the land-use model to reflect the population dynamics in Bangkok and will be able to analyze urban and transportation sustainability.
As mentioned above, this model has some considerable biases. Our analysis did not include business travel and goods transport. That a major reason for the bias in road speed estimation, congestion, subsequent detour trips, trip length, and travel time. That requires improving the accuracy and stability of mode choice and destination choice models. Furthermore, motorcycles were currently not represented in the simulation explicitly. Given the relevance of this mode in Bangkok, we may need to distinguish motorcycles from cars in the traffic simulation appropriately. It should also be noted that potentially induced travel was not represented here. It is conceivable that some trips that were suppressed due to high levels of congestion would be added if congestion was reduced. The model was not able to quantify this latent demand, as it was not observed in current traffic volumes.
There are various urban problems in Bangkok, including air pollution, rapid aging of society, securing mobilities for vulnerable people, suburban development on transportation, gentrification by inner-city development, the fiscal balance of the urban railway project, and induced urban development by transportation provision. Our model can be applied to a wide range of urban transport policy problems combined with the land-use model. Data limitation, however, requires various assumptions in modeling, which may explain at least in part the estimation bias. Our future studies will address these modeling issues.