How Many and Where to Locate Parking Lots? A Space–time Accessibility-Maximization Modeling Framework for Special Event Trafﬁc Management

Special event trafﬁc planning and management needs to accommodate high trafﬁc demand volume and special distribution patterns with dramatic structural deviations from the normal conditions. To provide sufﬁcient transportation service supply that matches non-typical demand needs, this paper explains how to systematically optimize the locations of park-and-ride stations, the number of additional parking lots, and the bus rapid transit schedules. The goal is to maximize the number of travelers who can complete their activity tours within a reasonable travel time budget. Based on a space–time network construct, this paper formulates a network design problem to maximize the system-wide transportation accessibility from different origins to activity locations at special event sites. A linear integer programing model is proposed to formulate the joint optimization of the location and capacity of parking lots associated with mega-event sites. Illustrative and real-world examples are used to examine the effectiveness and practical usefulness of the proposed modeling framework.


Motivation and Problem Illustration
For regular days, transportation demand can be represented as a repeated stochastic process with similar within-day dynamic patterns. A majority of transportation planning efforts have been devoted to accurately predicting and managing regular OD demand patterns, especially in the context of long-range transportation planning. On the other hand, it is extremely important to fully recognize various structural deviations of complex OD traffic desires from the normal condition, which could be caused by special events, severe weather conditions, dramatic responses of travelers under traffic incidents, as well as other demand management strategies.
Specific event traffic management needs to accommodate unusually high traffic demand volume and special spatial distribution patterns due to sports games, concerts, holidays, or other major reactional activities. In recent years, many regional planning organizations are also actively involved in the medium-term transportation planning of mega events, for example, major trade shows such as World Expo, large sport events such as Olympics Games, to boost the regional economic attractiveness and competitiveness.
A mega event can generate an extraordinarily high concentration of traffic superimposed on urban traffic networks for a few weeks. As an example, an Olympic game could not only attract daily travelers of 500,000 but also requires logistics personnel of 200,000 in size, which leads to almost 1.5 million trips on each day [1]. Accordingly, the transportation organization/authorities need to carefully plan and manage a coordinated multi-model infrastructure and service network. This will mitigate the potential traffic congestion while fully maintaining the accessibility to event sites. In this paper, we focus on how to optimize the location and capacity of additional parking lots. This will be particularly associated with park-and-ride stations, to allow a large number of visitors to successfully complete a trip chain within a reasonable travel time budget.
To improve the accessibility to the mega-event site, we will determine the locations and capability of park-and-ride (P&R) stations, which is critically important because they can facilitate the interchange between the private/lower occupancy mode to the public/higher occupancy mode. This will further help to complete trip chains through a sustainable multi-modal service network provision [2]. In recent Olympic Games, ranging from Atlanta in 1996 to London in 2012, P&R mode has demonstrated its capability in integrating an accessible public transport system to a well-planned and comprehensive transportation system for managing complex traffic demand at those mega events [3].
Planning and managing intermodal traffic demand/supply systems for special events is theoretically and practically challenging in its own right. To clearly illustrate the essential modeling elements and complexity of this problem, we use Fig. 1 to describe the problem of joint optimizing parking lot location and capacity for traffic demands in a special event. As shown in Fig. 1, visitors have trip desires to travel from multiple origins (O1, O2, O3, O4) to the mega-event site (D). Visitors can reach destinations through driving links (L1, L2, L3, L4), transit links (T1, T2), or P&R mode that combines the use of driving and transit links. Park-and-ride lots P1, P2, P3 are used to connect the road network and the transit network.
The final parking lot on link L4 allows driving-only travelers to park their cars and then walk to the mega-event site. In a complex user-equilibrium context, a traveler needs to make a departure time/mode/route decision to minimize his/her own traveling disutility that involves travel time/delay, transit and parking fares. From a transportation system management perspective, when the detailed transit fare has not been finalized yet, it is better to consider system-wide objectives such as minimizing the total travel time or maximizing the space-time accessibility to the event site, that is, the number of travelers who can complete their tours within a reasonable travel time budget. In this case, there are still a wide range of traveling options available and many of them are dependent on the locations and capacity of parking lots, traffic condition on the (driving) road network, as well as the capacity and schedule of transit services. For example, a visitor from origin O1 can drive along the road network O1 ? L1 ? L2 ? L3 ? L4 ? D. If link L3 is congested or parking lot on link L4 is saturated, he/she might consider an intermodal option through route O1 ? L1 ? L2 ? P2 ? T2 ? D by parking the car at P&R lot P2. If the capacity at parking lot P2 is still not sufficient, the visitor will need to consider driving a short distance to P&R location P1, through route O1 ? L1 ? P1 ? T1 ? T2 ? D, or taking a transitonly route through route O1 ? P1 ? T1 ? T2 ? D.
For a more complex situation, in a real-world large-scale network, large numbers of visitors travel from multiple origins to the special event site. They dynamically choose their routes to the destination, and strongly interact with each other. Under some emergent or unusual situations (e.g., accident, bad weather condition), congestion will form unexpectedly and travelers may dynamically switch their routes or be forced to change route informed by traffic managers. If there is sufficient information provision (such as VMS, radio, mobile phone Apps), travelers may dynamically change their routes, departure time, or even travel mode when they meet some heavy congestion, e.g., Gao et al. [4] discussed the crowd's dynamic route choice behavior during evacuations with sufficient information. Besides, changing mode can of course improve the accessibility to the event. In this paper, we consider more realistic but potentially complicated dynamic route/mode choice behaviors.

Literature Review
In general, there are a number of practical guidelines used in programing an event-based transportation management plan [5], which could have three major components; namely traffic management plan, transit plan, and travel demand management initiatives. These plans are critically  Fig. 1 Parking lot location and capacity optimization problem in an intermodal network needed to encourage event attendees to choose public transit as the primary mode especially in a mega event such as Olympic Games [6]. P&R is a key bridge to connecting the spatially distributed potential attendees to the public transit and then to the event site. There are a number of related studies in two categories, namely (i) survey-based empirical investigation and (ii) model-based optimization analysis. In the first category, a number of studies examine the practical effectiveness of P&R facilities, for example; Meek et al. [7] suggested that P&R might increase the average travel distance due to low load factors on dedicated buses, public transport abstraction, and trip generation. The second category of studies focus on modeling travelers' parking searching behavior and finding an optimal solution for designing and operating parking facilities. There are mainly five categories of research lines, e.g., parking searching and pricing, traditional network design problems, park-and-ride facility network design problems, intermodal simulation and assignment, and general space-time network. More technical details are provided in the references listed in Table 1.
The research on travelers' parking searching behavior aims to find the user-equilibrium solutions under different parking facilities' constraints [8][9][10][11][12][13][14][15][16]. Based on the route choice behavior under parking facilities, researchers have extended the traditional network design problems [17][18][19][20][21][22] to park-and-ride facilities' network design and location problems [23][24][25][26]. Considering the optimal location and pricing of a P&R facility simultaneously in a linear monocentric city, Wang et al. [23] aimed to find a deterministic mode choice equilibrium solution with objective functions as profit maximizing and social cost minimizing. Liu et al. [24] proposed an improved model based on deterministic continuum equilibrium that can be formulated through a super-network approach.
Another category of studies aim to simulate and implement the dynamic network loading process in both single-modal road network and multi-modal transit network [27][28][29][30][31][32][33][34][35][36][37][38][39]. These approaches can be applied to the intermodal infrastructure and service network design for special event management which has its own unique characteristics. Compared to the common system-optimal objective function that minimizes total travel time for a given OD demand, the special event organizers typically want to attract more visitors from different origins to attend the event in a reasonable time interval. Therefore, the space-time concept has been included to study the activity characteristics [40][41][42][43][44][45][46][47][48], e.g., finding the optimal trip chain under activity and time budget constraints, calculating the travel accessibility.
In this case, the space-time accessibility to the special event site is more relevant or important to achieving the overall management goal, compared to the simple mobility measure. While there are a wide range of studies (e.g., [40][41][42][43] examining accessibility-oriented strategies, a few researchers recently started systematically incorporating accessibility/connectivity measures in a network design modeling framework. To name a few, Santos et al. [44] introduced a transportation network design problem based on equity and accessibility. The activity-based network design problem studied by Kang et al. [45] aims to minimize both the network design costs and activity-related disutility using a bi-level model. In this paper, we extend a recent network design model by Tong et al. [48] that maximizes space-time accessibility, based on a spacetime prism analysis framework established by Miller [42].
Modeling the parking searching behavior under different pricing policies is a complicated and difficult task, this paper will only consider the simple but most important parking lot capacity constraint to describe the parking searching behavior without loss of generality. The remainder of this paper will be organized as follows. First, we formulate the related network design problem within a space-time network to improve individual visitors' accessibility of reaching desired activities at the special event from different origins. This is followed by a linear integer programing model for joint optimization of the location and capacity of parking lots. After presenting a Lagrangian decomposition-based solution approach [48,49], we use illustrative and real-world examples to examine the effectiveness and practical usefulness of the proposed modeling framework.

Model Formulation
In this study, we consider an intermodal urban transportation network, where travelers reach their destinations via different travel modes, such as transit (bus, Bus rapid transit (BRT), subway), driving (private car, HOV), and combined mode (park and ride). This multi-modal network can be modeled as a multi-layer network [31] for dynamic traffic assignment with integrated management strategies, where travelers using different modes have mode-specific penalty or costs and need to consider service or road capacity constraints at different layers.
The space-time accessibility measure for a special event needs to consider a trip chain that starts from a particular location, reaching the activity location, and then return back home within a travel time budget. In this case, if the road network is too congested or parking lots have limited capacity, then it is difficult to enable a large number of travelers to complete their activity chains. Thus, it is important to formulate an integrated network design model to maximize individuals' accessibility through parking lot location and capacity optimization while subject to the construction budget constraint and various flow and parking capacity constraints.

Space-Time Accessibility for Travelers
Attending Special Events The first modeling task is to build an intermodal spacetime network. The analysis system horizon is discretized into T time intervals, and each time interval represents r minutes or seconds. To better characterize the accessibility measure for a trip chain, we consider a traveling time budget TTB. If a person needs to spend T w minutes at the event location, then there is only a remaining time of TTB À T w for traveling activities between the origins and the event location. Similar to the space-time prism concept [48], we define the activity location at the special event site is accessible if it can be reached within time budget TTB, under actual traffic conditions (e.g., free flow or congestion) and parking lot availability restriction in the intermodal network.
To construct an equivalent network flow model with essential flow balance conditions at each vertex, as shown in Fig. 2, we add virtual traveling arcs to connect the origin node and the final time t 0 þ TTB. Within a network flow formulation framework, the virtual arc connects from the origin at the start time to the same origin node at the end time to allow all solutions to be feasible. For example, if the driver arrives at 18:30, which is not accessible within the total time budget (an unfeasible solution), and he/she is forced to reach destination directly through the virtual arcs.
By introducing waiting arcs at the origin and activityperforming arcs at the activity location, a feasible tour is generated within the space-time prism (with both total time budget constraint and activity-performing constraint). More specially, the examples in Fig. 2 consider the total activity time T w ¼ 3 h, and the total time budget is T ¼ 10 h. In Fig. 2a), a visitor makes the trip by the driving-only mode, he/she will depart at 8:00 and return origin at 15:00 when the traffic condition is reasonably good.
Under congested road traffic conditions, using transit lines and P&R facilities can improve the space-time accessibility. In Fig. 2b), the green line illustrates an example of an improved accessibility with better parking lot locations and sufficient capacity. The visitor can return back to the origin at 17:00 by P&R trip mode. However, the parking lot capacities and locations also affect the individual's accessibility goal. If the visitor has to drive a much longer distance to reach the parking lot, e.g., the reddashed line driving to parking lot PL #2 as shown in Fig. 2b), or the capacity of the parking lot is not sufficient, e.g., the blue-dashed line driving to parking lot PL#1, as shown in Fig. 2b), he/she still has difficulties to reach the special event site within the travel time budget (e.g., only arrive at the node N at 18:00 within the total time budget). Table 2 lists the notations for the key parameters and variables. The variables' definitions and formulations are extended from the framework by Tong et al. [48].

Math Formulation
The proposed transportation network design problem considers a set of nodes and links. In the constructed space-time network, the space time arcs A are defined by a 5-index matrix x i;j;t;s ðpÞ, and the transportation space-time arcs have zero inaccessibility costs c i;j;t;s p ð Þ, while the c i;j;t;s p ð Þ = 1 for virtual arcs associated with the virtual arc for the inaccessible activity chain is illustrated in Fig. 2. The binary variable y i;j represents whether the parking lot (i, j) is selected to be constructed.
Given travel time budget TTBðpÞ, origin o p , departure time s p for each agent p, the problem under consideration aims to minimize the system-wide space-time inaccessibility measure, subject to a number of capacity and construction budget constraints. Then, the inaccessibility minimizing, i.e., accessibility maximizing, objective function is reformulated as an integrated optimization problem for joint locating parking lots and corresponding parking space capacity.
Objective function The objective function to be minimized is now formally stated as Eq. (1).
That is, if a traveler p can perform the activity at the special event site d within the given time budget, then its resulting activity chain-related inaccessibility cost is zero, otherwise, the cost is 1 by traveling through the virtual arcs. Equivalently, this objective function wants to minimize the total number of inaccessible trip chains across all agents who have desires to make the trips. It should be mentioned that, the cost matrix in the objective function [1] is for inaccessibility computation, which could be quite different with the travel time. We will use s i,j to describe the free flow or experienced time-dependent travel time for link (i,j). Travel time is embedded in the arcs (i,j,t,s = t?s i,j ) in the space-time network.
There are 6 types of constraints in our model, namely the space-time flow balance constraints, activity-performing constraints at special event site, traffic flow capacity constraints, park-and-ride facility constraints, transit service constraints, and total construction budget constraint.
(1) Space-time flow balance constraints X j;s ð Þ2Q: i;j;t;s ð Þ2A x i;j;t;s p ð Þ À X j;s ð Þ2Q: j;i;s;t ð Þ2A (2) Activity-performing constraints at special event site X The space-time activity-performing arcs at special event site and virtual traveling arcs for passenger p.
(3) Traffic inflow and spatial capacity constraints A spatial queue mesoscopic traffic flow model is applied to describe the traffic flow constraints. That is, the total inflow r constraint: The number of cumulative arrival and departure agents on the link ði; jÞ can be represented as: Here, for link (i,j), s i;j is the free flow or experienced travel time, L i;j is the length, and K jam i;j is the jam density. For more details on queue-based dynamic traffic simulation, please refer to Zhou et al. [48].
(4) Park-and-ride facility constraints (i) Capacity associated with cars arriving and departing at parking lots: the cumulative number of arrival agents minus the cumulative number of departure agents could not exceed the space capacity of the parking lot.
Transit service constraints: transit vehicles can only drive on specific space-time arcs defined by the given transit schedule. The number of passengers served by a transit schedule space-time arc should not exceed the passenger-carrying capacity of a transit v.
(6) Total construction budget constraints X The total construction cost for selected parking lots should not exceed the total construction budget.

Lagrangian Decomposition-Based Heuristic Solution Algorithm
In this section, a Lagrangian decomposition (LD) algorithm is introduced to reformulate and further decompose the integer linear programing problem [48,49]. This relaxation scheme allows parking lots to be overloaded, but penalizes it in the objective function. According to recursive calculation, the penalties will be updated iteratively, finally giving the marginal benefit of increasing capacity at those lots. First, we have the original problem P1.  (6), which links two sets of variables X and Y, we have a combination of two relatively easy-to-solve problems: P X as a set of constrained time-dependent routing problem for passenger agents embedded in a multi-modal dynamic traffic assignment program subject to constraints (2,3,4,5,7,8), and P Y as a knapsack problem subject to the total construction budget constraint (9). Interested readers are referred to Tong et al. [48] for detailed descriptions of LD solution approach with a similar modeling framework.
The dualized Problem P2 can be written as This problem is to minimize the total number of inaccessible trip chains across all agents and the total penalty cost under P&R capacity constraints.
Based on a dual optimization perspective, the Lagrange multipliers p i;j;t;s can also be interpreted as shadow price associated with P&R capacity constraints. The solution steps of the proposed algorithm can be listed as follows: Step 1: Initialization Set iteration number k ¼ 0; the set of available parking lot locations are given in terms of links in set E P ; and total construction budget TCB.
Choose positive values to initialize the set of Lagrangian multipliers p i;j;t;s : Step 2: Solve decomposed dual problems Step 2.1: Solve sub-problem P X using an enhanced multi-modal dynamic traffic assignment (DTA) simulator with a time-dependent least cost path algorithm and find a path solution XðpÞ for each agent p.
A spatial queue-based traffic flow simulator, such as DTALite [39], is used to ensure the traffic inflow and spatial capacity constraints (4) and cumulative flow count definitional constrain (5) and transit service capacity constraint (8) are all satisfied for all traveling agents through a multimodal dynamic network loading (DNL) program. Specifically, in a transportation network, a node is connected to different incoming links and outgoing links, and each link has two buffers in DNL, namely entrance buffer and exit buffer to facilitate traveling agents' transfers between links. These two buffers on each link are commonly implemented as first-in-and-first out (FIFO) queues. When the required link inflow and outflow capacities are available, an agent can move from the exit buffer of an upstream link to the entrance buffer of the downstream buffer.
To handle the remaining constraint sets (2, 3), namely space-time flow balance constraints and activity-performing constraints, with the dualized objective function in Eq. (10), a time-dependent routing problem with a set of constraints is solved for passenger agents. Specifically, each arc i; j; t; s ð Þin the available P&R locations (i.e., E P ) has an additional cost of p i;j;t;s for the dualized capacity constraints, which is equivalent to the estimated travel time penalties when agents use park-and-ride facilities. The space-time flow balance constraints are satisfied automatically in the routing algorithm, and the remaining activityperforming constraints at special event site can be handled through a simple decomposition to two trips. One from the origin to the special event site, and the other from the special event site back to the origin. Under road traffic congestion or large travel time penalty associated with p i;j;t;s at the parking lots, the travel time budget constraint TTBðpÞ might not be satisfied, and the routing algorithm will minimize the number of inaccessible virtual arcs as the total disutility. It should be remarked that, even there are optional P&R capacity available, some travelers (with the goal of accessibility maximization) could still select driving-only mode, if the related road traffic condition is less congested.
Step 2.2: Solve sub-problem Py, for example, using a dynamic programing algorithm, to find a value for Y. The dual cost p i;j;t;s associated with the dualized P&R capacity constraints will encourage the decision makers to select the most cost-effective P&R station location/capacity allocation option to maximize the total dual cost (i.e., profit) to be collected through a knapsack modeling framework, which is equivalent to maximize P i;j;t;s ð Þ2A:ði;jÞ2E P p i;j;t;s Â È cap P i;j Â y i;j h i g subject to constraint (9).
Calculate primal, dual, and gap values of P2.
Step 3: Update Lagrangian multipliers Update Lagrangian multipliers p i;j;t;s using subgradient the step length at iteration k.
Step 4: Termination condition test If k is less than a predetermined maximum iteration value, or the gap is smaller than a predefined toleration gap, terminate the algorithm; otherwise k ¼ k þ 1 and go back to Step 2.

3-Corridor Network for Testing Model Choice Sensitivity of Multi-Modal Network Loading
A 3-corridor network is constructed, as demonstrated in Fig. 3, to illustrate the proposed accessibility-maximization model. Node 1 represents the original location of travelers and node 12 represents the special event or megaevent location. The red lines represent a roadway system allowing passenger cars. The green lines represent dedicated BRT lines, which assign the right of the way to BRT vehicles only. The driving link 4-8-7 and 6-9-7 allow passenger cars to drive to the parking lots and connect to BRT stations. Three demand types are considered in this intermodal network, which includes driving alone (demand type 1), BRT (demand type 2), and Park & Ride (demand type 3). For simplicity, the total demand is assumed to be fixed, including 510 driving-only users and P&R users, 800 transit users, and 560 vehicles per hour background vehicles. During the trip, only en-route users can change their trip modes (P&R, driving only).
To investigate the effect of different levels of background traffic (driving only) on roadway network, a comparison of base and high driving-alone demand is conducted. Using DTALite [39], a queue-based mesoscopic traffic simulator, and the proposed LD algorithm, the problem is then solved and convergence is achieved by 20 iterations. The two demand scenarios are listed in Table 3. It can be seen that with the capacity of parking lot increasing, the number of pedestrians choosing P&R mode will increase from 50 to 396, and the average travel time can reduce to 147.9 min. The accessibilities of different time budgets are also listed in Table 3. The large parking lot capacity could also lead to a better accessibility, e.g.,   Figure 4 shows the relative gap and LB (lower bound) and UB (upper bound) with different scenarios of each iteration with our algorithm. Typically, the lower bound estimates are improved significantly during the first few iterations, and the duality gap between the upper bound and the lower bound estimates can be reduced dramatically to a relatively small difference after a certain number of iterations. In scenario 1, P&R is located on link 8 ? 7, and the inaccessible agents is about 242; while in scenario 2, P&R is located on link 9 ? 7, the inaccessible agents decrease to 208. It seems that the P&R station is beneficial to be located on link 9 ? 7. Here, the locations of P&R facilities are decision variables in this paper, and the variables can be regarded as space-discrete variables. For larger networks, we can also predetermine some alternative patterns to reduce the feasible solution space and make the problem easier to solve.

Real-World Testing Case Study
In this paper, we use a real-world testing case study for further examining the practical usefulness of the proposed methodology. Specifically, we consider International Horticultural Exposition 2019, Beijing, China (hereafter referred to as Beijing Expo 2019) which will be held in the Yanqing district, Beijing, from 29th April to 7th October 2019. Beijing Expo 2019 is expected to have more than 100 official exhibitors, more than 100 other exhibitors and more than 16 million visitors as an initial traffic demand estimation, and the expected range of potential visitors will be 34-37 million according to additional surveys.
As shown in Fig. 5, the proposed analysis methodology is applied to the Yanqing District study area to examine the effectiveness of the parking lot location and capacity allocations on traffic flow conditions/mobility and spacetime accessibility. In this network, there are 66 OD zones, 1519 nodes, 3299 links with 22 different multi-modal link types. To meet the dramatically increased traffic desires for this mega event, transit lines and parking lot facilities need to be well designed. At present, there are 2127 transit lines and 8 planned P&R facilities.
According to the proposed LD-based solution algorithm, we generate several representative scenarios with suggested/optimized parking plot locations and capacity allocation schemes. That is, Scenario 1 considers only transit and driving network without P&R or BRT; Scenario 2 considers one new BRT line, but no P&R facilities to transfer to the BRT line; Scenarios 3-6, under different total construction budget TCB, 5-8 parking lots can be built to allow visitors to transfer to the BRT line. To investigate the effectiveness of BRT line and parking lot facilities on the travelers' space-time accessibility, 142,318 agents (visitors) are assumed to attend to the special event on a peak-hour period, and they can travel by three modes, i.e., driving only, transit only (traditional transit and BRT), and P&R.
The preliminary numerical results based on our initial OD demand estimates with limited survey data are listed in Table 4. It can be seen that if all agents travel by drivingonly mode, the average travel time is 74.9 min, and only 23 % of travelers can reach their accessibility goal. If BRT lines are built, the travel time can be reduced to 59.1 min, and the overall accessibility ratio will increase to 47 % significantly. When 5 parking lots are built, the transit and BRT lines can take more than 60 % passengers, and the average travel time can be dramatically reduced to 19.2 min, and the tour-level or agent-based accessibility ratio will then approximately increase to 83 %. As a result, optimized P&R location and capacity allocation and BRT facilities are extremely helpful in terms of increasing visitors' space-time accessibility goals. It should be remarked that, if more parking lots (6 or more) are built for additional locations, the overall accessibility could not be further improved significantly based on our simulation results.

Conclusions
In an intermodal traffic network, travelers' accessibility measure highly depends on the traffic flow condition, parkand-ride facility, and transit schedule. This paper describes a space-time representation of an intermodal traffic network. It also integrates traffic flow, transit schedule, and park-and-ride facility to establish a linear integer programing network design model to maximize the tourist's  accessibility to a mega event within a multi-modal network. A Lagrangian relaxation and decomposition solution approach is developed to efficiently solve this problem. The numerical results show that building parking lots at proper locations with an optimized parking space capacity could be useful in terms of allowing more intermodal transfers, reduce individual's travel time, and improve the special event visitors' accessibility. In practice, the megaevent park-and-ride lots planning problem should be integrated with long-term transit development planning decisions to fully utilize/balance the potentially unused capacity of the facilities. Future research will also focus on algorithm enhancements, OD demand calibration/validation, as well as valuating potential applications with other network design models.