Research Article  Open Access
Xinyuan Chen, Inhi Kim, "Modelling RailBased Park and Ride with Environmental Constraints in a Multimodal Transport Network", Journal of Advanced Transportation, vol. 2018, Article ID 2310905, 15 pages, 2018. https://doi.org/10.1155/2018/2310905
Modelling RailBased Park and Ride with Environmental Constraints in a Multimodal Transport Network
Abstract
Environmental sustainability is a significant aspect in the sustainable development of modern urban cities, especially in the road transport system. As traffic demands increase, public transport requires more promotion to accommodate the increasing travel demands while maintaining the environmental quality. Public transport, however, is less attractive in vast suburb areas mainly due to its longer travel distance and waiting time. Therefore, this paper proposes a railbased ParkandRide (RPR) scheme to promote public transport in the multimodal transport network. To remedy the heterogeneous distribution of vehicle pollutants in the network, regulations in environmental sensitive districts are required and studied in this paper. To quantitatively evaluate and analyse this joint RPR and environmental regulation strategy in multimodal transport systems, this paper develops an environmental constrained combined modal split and traffic assignment (ECCMSTA) model. The proposed formulation adopts the concept of fixpoint to reformulate the nonlinear complementarity conditions associated with the combined modal split and user equilibrium conditions, which is subsequently incorporated into a VI formulated nonlinear complementarity conditions associated with environmental constraints. The proposed VI formulation can handle a general constraint structure, which enhances the modelling adaptability and flexibility. The strictly monotone and Lipschitz continuity properties of this model are rigorously proved, giving rise to efficient algorithms for the model. A customized projection based selfadaptive gradient projection (SAGP) algorithm is then developed. Numerical studies demonstrate that the ECMSTA model could enhance the behavioural modelling of network users’ travel decisions and assist in quantitatively evaluating the effectiveness of RPR schemes and environmental regulations.
1. Introduction
Environmental protection is an important component in the sustainable development of modern urban cities and closely related to the transport system. According to Hong Kong Air Pollutant Emission Inventory (2016), the amount of CO generated in the year of 2016 in Hong Kong is 58600 tones and among them, road transport contributes to 31500 tones and shares 54% of the CO emissions. Also, the US’s environmental protection agency reported motorized vehicles in the United States contribute to 86% of CO emissions, 62% of emissions, 12% of PM2.5 emissions, and 7% of SO_{2} emissions [1]. Kitthamkesorn et al. [2] pointed out that only about 20% of the world’s town residents enjoy good enough air quality as measured by the levels of emissions and over 600 million people in urban areas worldwide were exposed to dangerous levels of trafficgenerated emissions from congestion. In addition, these emissions have not only adverse effects on human health but also a source of greenhouse gas that contributes to global warming [3].
With these adverse effects, the promotion of alternative green transport modes is thus pressing for the governments and transportation authorities to accommodate the growing travel demand while maintaining the environmental quality. Public transport aiming to satisfy the travel demand of urban commuters safely, quickly, and affordably is therefore a significant alternative to private transport. However, in vast suburb areas which are far away from an urban centre, public transport becomes less attractive, mainly due to the unreliable time table, tortuous bus itinerary, much lower travel speed, and lack of door to door services [4]. Largescale deployments of public transport system in the remote areas are hardly seen as an economical available strategy. Park and Ride (P&R), therefore, becomes an important strategy to promote public transport since it enables travellers to drive when the road is not congested and to change a mode to transit accessing to CBD without traffic congestion. Since its first implementation in 1930s, P&R has been widely recognized as a significant traffic management strategy and becomes an indispensable part of majority modern urban cities.
P&R can be categorized into busbased P&R (BPR) and railbased P&R (RPR) based on its ride component [5]. In a metropolis with advanced rail system, railbased P&R is superior to busbased P&R, due to the faster speed, higher level of safety, comfortability, and its advantages in environmental protection (zeroemission).
Despite the fact that environmental deterioration could somehow be alleviated by promoting public transport through RPR schemes, the environmental requirements are usually not strictly met in certain districts mainly due to the heterogeneous distribution of pollution in the network [3]. To remedy this limitation, specific environmental regulations, therefore, are required to certain district, e.g., the city business district (CBD) with many tall buildings, to protect the environmental quality in this area.
Therefore, we aim to use RPR services to assist mode shift to reduce the overall vehicle travel demands in the road network and then apply areabased regulation scheme to deal with the environmental deterioration caused by heterogeneity of traffic flows. The environmental protection approaches are integrated as a combined strategy in this paper. Previous studies revealed that a systematic study is necessary to properly understand and evaluate the effectiveness of the environmental regulations, as well as the RPR scheme. Thus, this paper focuses on the modelling and solution algorithm development of the proposed environmental protection strategy.
2. Literature Review
Road network equilibrium approaches are commonly preferred for estimating vehicle emissions and noises [3], because the distribution of vehicle pollution heavily relies on the states of traffic flows [6]. Other than the overall pollution level, the spatial distribution level should be controlled since it is heterogeneous.
One line of the research developed models by setting environmental protection as objectives [3], modelling environmental considerations as an objective function. Rilett and Benedek [7] proposed the concept of emission optimum (EO) and emissionoriented user equilibrium (EOUE), in which the former one assumes travellers make route choices with the objective of minimizing the total emissions or noises, and the latter includes the toll charges to represent the impact on the environment. Later, Benedek and Rilett [8] made extension by modelling the assignment of vehicles based on equity which considers the objectives of neighbouring residents, rather than the system operator’s or drivers’ objectives. Sugawara and Niemeier [9] used the emissionsoptimized trip assignment model to estimate the maximum carbon monoxide emission reduction under varying congestion levels.
Although the overall environmental quality could somehow be improved by taking the environmental considerations as the objective, the pollution level is not always strictly met on each link and zone in the road network. Therefore, another line of the research focuses on modelling the emission as environmental constraints. Chen, Zhou and Ryu [10] developed a traffic assignment model by considering various physical and environment restrictions as side constraints, which is formulated as a VI problem and solved by a predictorcorrector (PC) algorithm. Li et al. [11] developed a toll design model by simultaneously considering congestion and environmental externalities. Recently, Xu, Chen, and Cheng [1] reformulated the environmental constrained traffic assignment model through a smooth gap function to enhance the modelling adaptability and flexibility.
As a rational extension of traffic assignment problem, Nagurney, Ramanujam, and Dhanda [12] extended the environmental constrained traffic assignment problem to the case of multimodal traffic network. However, the model presented in this paper simply follow the deterministic user equilibrium (DUE) principle for both modal split and traffic assignment. In order to evaluate the environmental regulation strategy in the multimodal transport network, a better behavioural model is needed. Therefore, the combined modal split and traffic assignment model (CMSTA) is adopted in this paper. In the literature, the CMSTA problems are widely adopted to evaluate the P&R schemes. Fernandez et al. [13] conducted an initial work to model P&R in deterministic case. Later, this work was extended to the case of asymmetric link travel time functions by García and Marín [14] and to the logitbased stochastic model by Li et al. [15] and Lam et al. [16]. Recently, Liu and Meng [5] provided a formulation for the busbased P&R in a network with the congestion charge. Pineda et al. [17] and Song, He, and Zhang [18] developed an integrated P&R facilities and transit service optimization model. Kitthamkesorn et al. [2] and Kitthamkesorn and Chen [19] considered the influence of “go green” transport modes on environmental protection by taking use of recent advances in logit model. Liu et al. [20] proposed a new concept of P&R, remote park and ride, and developed an exact solution method for the optimal design of parking facilities, following the nonlinear valid inequality technique [21, 22]. Some other studies focus on the P&R schemes in a linear corridor, e.g., Liu et al. [23], Wang and Du [24], and Du and Wang [25]. However, these studies mainly focus on congestion alleviation. As abovementioned, the environmental requirements in certain areas are not considered. Therefore, a more general modelling framework that incorporates environmental requirements is needed. To this end, this paper aims to build an environmental constrained combined modal split and traffic assignment (ECCMSTA) model to better evaluate and analyse environmental sustainable transportation development strategies.
This paper distinguishes itself from previous studies in threefold aspects. First, we propose ECCMSTA, CMSTA model with side constraints, which could handle the general environmental constraint structures, including linear or nonlinear, link based, or areabased constraints. Also, the general link/route cost structure is considered here. Second, by incorporating general side constraints into CMSTA, solving is much more complex compared to the traditional side constrained unimodal traffic equilibrium problems. To explicitly consider the general environmental constraints and make use of existing efficient algorithms, we propose a VI formulation by adapting the concept of fixpoint. This VI formulation is strictly monotone and Lipschitz continuous, permitting a number of existing efficient algorithms for its solution [10, 26]. Third, even though the solution methods for the traffic assignment problems with side constraints have been well studied with numerous achievements [1, 10, 12, 27], solutions for the multimodal traffic equilibrium with side constraints are relatively limited. Recently, Meng, Liu and Wang [28] proposed a CA incorporated selfadaptive PC algorithms to solve asymmetric stochastic user equilibrium (SUE) problem with elastic demand and link capacity constraints. Inspired by Meng, Liu, and Wang [28], we extend this approach to multimodal traffic equilibrium problem and develop a customized projection based algorithm with a selfadaptive step size scheme which permits the convergence under a much milder assumption of the underlying mapping.
The remainder of this paper is organized as follows: we present assumptions and an ECCMSTA model for the environmental regulations in the multimodal transport system with RPR scheme in the next section. A general VI formulation and a projection based solution algorithm are provided in Sections 3 and 4. Section 5 provides two numerical studies to validate the proposed methodology. Finally, Section 6 concludes this study.
3. Problem Description
Considering a specific environmental regulation strategy, to analyse its impact in the multimodal transport system, we build a model for the equilibrium flows following environmental requirements.
Given a multimodal transportation network represented by a strongly connected topology graph where is the set of nodes and is the set of links. The origins, destinations, and RPR facilities are scattered in the network and represented by certain nodes, . The set of links can be further categorized into three subsets: set of autolinks ; set of rail links ; and set of P&R links ; .
3.1. Modelling Environmental Requirements
As abovementioned, the environmental requirements are not strictly met in the road network. Thus, it is necessary to explicit considering environmental protection thresholds to prevent environmental deterioration in certain zones. In the literature, the environmental constraints can be expressed into a general form:where is a component of and represents a set of environmental sensitivity elements. In the case of link based environmental protection strategy, represents a link which is assumed in this paper. is the vector of link flows; is a function to measure the amount of vehicle pollution and environmental thresholds. Equations (1) ensure the amount of pollution no larger than predetermined thresholds. Note that (1) are general formulation that could measure different types of vehicle pollutants [1]. Herein, we make assumption as follows.
Assumption 1. The environmental threshold measure function , is a continuously differentiable and strictly monotone increasing function of the link flows .
Suppose the total number of environmentally sensitive areas in the network is . The inbound flows of these areas should be restricted according to environmental protection strategies. All the links in the area are grouped into set , where . Then, the specific link based environmental constraints in these areas could be expressed aswhere is the amount of vehicle pollutants per vehicle on link and is the environmental protection threshold on link imposed by the operators.
3.2. Network Assessment
Due to the fact that the travel decisions of a network user consist of modal choice and route choice, a combined modal split and traffic assignment model is needed. The travel impedance on an auto or rail link is clearly influenced by its flows. For the RPR links, each of them always connects an autolink and a rail link, so that the RPR flows could transfer. The total travel impedance on a RPR link includes the time spent on parking, waiting/boarding the dedicated bus as well as the bus invehicle travel time, which is clearly flowdependent. We then make the following assumption.
Assumption 2. The link travel time is a continuously differentiable, monotone increasing of the link flow .
Let denote the set of paths under ; denotes the flow on path for mode ; denotes the demand of mode ; and denotes the total travel demand between OD pair . We have the following flow conservation conditions:where (3) and (5) define demand/flow conservation conditions. Equations (4) and (6) are the nonnegativities. Equations (7) define the linkpath flow relationship. Constraints (1), (3)(7) define the feasible set of of where is the vector of link flows; is the vector of path flows; and is the vector of modal demands.
Let denote the set of OD pairs; denote the total travel demand for ; and denote the set of modes between OD pair . The multinomial logit (MNL) model is adopted to analyse the modal split pattern. The probability of choosing any mode between an OD pair equalswhere is the dispersion parameter and is the deterministic component of the utility of mode between OD pair , which can be taken as the equilibrium path travel time of mode . The satisfaction of MNL, , equals , which is concave function [28].
The user equilibrium path flow solution is then formulated by the following conditions:where is the Lagrange multiplier associated with the constraints (5), which represents the travel impedance of using mode . At equilibrium, the travel time of a used path equals . All commuters under mode have identical path travel time, and no traveller can unilaterally reduce his/her travel time by changing routes.
4. Variational Inequality Formulation
In this section, we propose two equivalent VI formulations for the proposed ECCMSTA model. In order to strengthen the flexibility and adaptability of the proposed model, we further make the following two assumptions.
Assumption 3. The link travel time function is asymmetric, i.e., is not only affected by its own flow , but also flows on some other links, which is expressed as .
Assumption 4. The path travel time is nonadditive to its links, i.e., , where if path uses link and 0 otherwise.
Then, the VI mode for ECCMSTA is formulated as follows:where is the travel time on route , where . represents the entropy term in travel cost function of mode , ; and are the optimal solution of this problem; and are the vector of optimal solution.
Let be the vector of and be the vector of , i.e., and . The VI formulation can be simplified as
Note that the inclusion of nonlinear environmental constraints makes the feasible set more complex in general. Existing solution algorithms for CMSTA cannot be directly adopted to solve it. Most of existing studies focus on solving the simplified linear link capacity constrained traffic equilibrium problems by a penalty based algorithm. However, due to the complexity of ECCMSTA (simultaneous achieving equilibrium for both route flow and modal demand), reformulating it into an integrated formulation or directly solving it by a penalty/gradient based algorithm is less effective. Therefore, this paper develops an alternative approach to reformulate the ECCMSTA which permits existing solution algorithms for its solution.
4.1. Alternative VI Formulation
Daganzo [29] built a fixpoint model in terms of link flows for the asymmetric SUE problem with fixed demand and showed availability of the MSA, which was extended by Cantarella [30] to the case of asymmetric SUE problem with elastic demand. Recently, Meng, Liu and Wang [28] further made extension by incorporating link capacity constraints. This problem was formulated as a VI problem by incorporating the fixpoint formulated asymmetric SUE problem with elastic demand. Following the same technique, we propose a fixed point reformulated CMSTA model which is expressed as follows:where is the equilibrium link flow solution, is the probability of mode choice defined in (8), and is the link usage probability defined byand is the probability of path perceived as the shortest one among the path set . In the deterministic case, equals one if path is the shortest path among alternatives and zero otherwise.
A set of generalized SUE conditions with elastic demand and capacity constraints was proposed by Meng, Liu, and Wang [28] for the optimal solution of SUE problem with elastic demand, asymmetric link travel time function and link capacity constraints. These conditions can be further extended to the multimodal transport system which is the problem addressed in this paper. Let denote the equilibrium link flow pattern of the addressed problem if and only if there exists a vector of Lagrangian multipliers corresponding to environmental constraints, denoted by . where is defined byand if path goes through and , otherwise. Based on the optimal Lagrangian multiplier , the ECCMSTA is identical to the CMSTA in terms of generalized link travel time functions. Equations (16)(18) imply that if there is no environmental constraint on any link , .
Note that the complementary slackness conditions (see (16)(18)) could be regarded as a nonlinear complementary problem (NCP) of the vector of Lagrangian multipliers . Due to the feasible set of is defined on the whole nonnegative orthant, this NCP problem can be equivalently formulated as a VI model. Hence, the alternative formulation of VI (11) is proposed as follows:where denotes the feasible set. is the dimensional vector functions, which measures vehicle pollution: is thus the optimal solution of the VI problem (21). Herein, is a link flow solution for the CMSTA with the generalized link travel time functions (20) (without environmental constraints). In other words, is a solution of the fixpoint model with the Lagrangian multiplier :where denotes the vector of generalized link travel time functions.
Note that the feasible set of is unbounded, which makes it difficult to shown the existence of a solution to VI model . Fortunately, according to Meng, Liu and Wang [28], one can always find an upper bound for . Therefore, we define a subset of which is bounded by . The restricted VI model is subsequently formulated asThe feasible set is a nonempty, closed, and convex set. To show the existence of solution to as well as the convergence of its solution algorithm, the following properties of vector function are proposed. Rigorous proofs of Propositions 5 and 6 are given in the Appendices A and B, respectively.
Proposition 5. Vector function is strictly monotone on ; namely,
Proposition 6. Vector function is Lipschitz continuous on ; namely,
Though we have Propositions 5 and 6, the existence of a solution to VI model can be ensured. We proceed to prove the equivalence between the VI model and the generalized CMSTA conditions and summarized in Proposition 7. Rigorous proofs are given in Appendix C.
Proposition 7. is a solution of if and only if and fulfil the generalized CMSTA conditions (15)(18).
5. Solution Algorithm
As the vector function of the VI model is an implicit function, it is burdensome to directly evaluate its gradient. On the other hand, with any given Lagrangian multiplier vector , the value of vector function can be evaluated by solving the CMSTA with generalized travel time with no constraints, which can be efficiently solved by a number of existing solution algorithms.
Projection type algorithms have been thoroughly studied and proven to be efficient in solving traffic equilibrium problems [31–33]. Among them, the gradient projection (GP) is commonly adopted due to its concise form and good performance, for which a recent survey can be found in Ryu, Chen, and Choi [34]. The selfadaptive gradient projection (SAGP) algorithm was originally proposed by He et al. [35], where the gradient projection operation is adopted to establish the descent direction, coupled with a selfadaptive step size scheme to permit the convergence under a much milder assumption of the underlying mapping, and thus can handle general cost mapping structures and generalized side constraints. Later, Han and Sun [36] adopted this method to solve the asymmetric traffic assignment problems. Chen, Zhou, and Xu [32] further used this method to solve the asymmetric traffic assignment problem with nonadditive route cost functions. Inspired by previous studies, this paper presents the SAGP algorithm specifically designed for solving the CMSTA with generalized travel time functions.
As for the Lagrangian multipliers , the projection operation on of the VI model is straightforward, as well. Propositions 5 and 6 imply the monotony and continuity of the VI model . Therefore, the projection based algorithm is also suitable to evaluate the Lagrangian multipliers solutions. To this end, this section proposes a nested SAGP algorithm to solve the ECCMSTA model which is presented as follows.
Step 0 (initiation). Set parameters of SAGP within the given range: , , , , , . Set the iteration counter: , and the initial solution which should be a feasible solution: .
Step 1 (selfadaptive scaling procedure). This step aims to find the largest step size that satisfies the given criteria in (27). To this end, let , and set nonnegative integer . We perform the following procedures, if the given criterion is satisfied, terminate and go to Step 2; otherwise, let , and perform the following procedures again.
Step 1.1. Given vector , evaluate the CMSTA with generalized travel time functions (23), using SAGP algorithm, and obtain the equilibrium flows . The solutions are subsequently used to evaluate the environmental measurement function .
Step 1.2. Perform the projection operation for to obtain a tentative vector , i.e., .
Step 1.3. Let . If the following inequality is satisfiedgo Step 2. Otherwise, let , update step size , and go back to Step 1.1.
Step 2 (selection of ). If the following condition is satisfied, let .
Step 3 (convergence test). If the convergence criterion is satisfied, terminate; otherwise, let iteration counter: , go back to Step 1.
Note that the CMSTA with generalized travel time functions in Step 1.1 is also solved by the SAGP algorithm, which is not repeated here for conciseness.
6. Numerical Examples
In this section, two experiments based on a linear network and the SiouxFalls network are provided to demonstrate the features of the proposed ECCMSTA and validity of the proposed solution algorithm.
It is well known that carbon monoxide (CO), carbon dioxide (CO_{2}), sulphur dioxide (SO_{2}), nitrogen oxides (), respirable suspended particulates (RSP or PM_{10}), fine suspended particulates (FSP or PM_{2.5}), and volatile organic compounds (VOC) are the main types of vehicular emission pollutants that are critical for air quality protection. Among others, CO is taken as an important indicator of the level of vehicle atmosphere pollution, mainly because vehicles contribute to the majority of CO emissions in the atmosphere; CO is one of the most critical pollutants among the various types of vehicular emissions; the emission functions of other pollutants are similar to that of CO [11]. In addition, this consideration has also been used in a number of recent studies, e.g., Yin and Lawphongpanich [37]; Nagurney, Qiang and Nagurney [38]; Chen, Zhou and Ryu [10]; Li et al. [11]; Xu, Chen and Cheng [1]; Yang, Ban and Ma [39]; Kitthamkesorn et al. [2]. Therefore, we use CO as an illustration to model vehicle emission. The nonlinear macroscopic formula of CO emission rate proposed by Wallace et al. [6] is used for analysis as shown below:where is the CO emission rate on link measured in grams per hour. is the length of link measured in kilometers, and is the link travel time measured in minutes. This nonlinear function can be shown to be convex and continuously differentiable. The Bureau of Public Road (BPR) type function is adopted for the link travel time, i.e., for road links, and for rail and RPR links. The dispersion parameter of MNL mode is set as .
6.1. A Linear Corridor
As shown in Figure 1, consider a linear corridor , with one rail line (dotted line) from suburb to the urban city area and another highway (solid line) connecting . The OD demand is assumed to be 500 per hour. One RPR site is constructed at , adjacent to the train station . Travellers transfer from the RPR site to the train station through the RPR link (link 5). Link 2 passes into an environmentally sensitive urban area, which exerts restriction on the total CO emissions on link 2. Using the CO emission rate function defined by (29), the nonlinear environmental constraint on link can be expressed aswhere is the CO emission threshold on link . Attributes of the network are given in Table 1.

Setting the environmental constraints at 500, we first examine the correctness of the results. As expected, the results in Table 2 satisfy the modal demand conversation condition and the MNL based modal choice condition. The Lagrangian multiplier on link 2 equals 1.73, permitting the emission on link 2 subject to the requirements.

We proceed to explore the influence of the environmental constraint on the flow shift of modal demands, link flows, and total system travel time by varying the emission constraint between 0 and 800. The corresponding modal demand shift pattern under environmental constraint is given in Figure 2. We can see that as increases (relax the CO emission constraints), the auto demand first increases and then becomes stable; meanwhile, the demands of rail and RPR show the crosscurrent. It makes sense because relaxing is similar to increase the road capacity. When the emission level is beyond 631.47, constraint (30) is not blinding anymore which means the PRP facilities are able to achieve the environmental protection requirements without additional regulations or restrictions. The dotted line indicates the evolution of Lagrangian multiplier on link 2, . As the environmental regulation can be strengthened, keeps increasing and trends to be infinity, when .
Detailed results of link flows, link travel time, and Lagrangian multiplier are shown in Table 3. As marked by bold and italic text in Table 3 the minimal system travel cost is achieved when is increased to around 631.47. In other words, setting is equivalent to have no CO emission constraint on link 2, which reduces the problem into CMSTA. Meanwhile, the Lagrangian multiplier reduces to zero as well indicating that there is no restrictions for vehicles which come into this area.

6.2. SiouxFalls Network
Another numerical example built from the SiouxFalls network (Figure 3) is used to further access the performance of the proposed model and solution algorithm for solving the environmental constrained CMSTA problem. This study network consists of 30 nodes and 76 links. The area highlighted by the red ellipse is the CBD district. It has four rail lines between the suburb areas and the city loop. Detail attributes of this network are tabulated in Appendix D. We focus on the travel demands in commuting hours, where the majority of travel demands origin from the suburbs and go into business districts. There are 24 OD pairs and the OD demand data are given in Appendix E. Four sites (nodes 3, 6, 19, and 14) which locate close to rail lines are selected as RPR car parks where RPR users park their car and transfer to rail mode.
The column generation scheme is adopted in this case and incorporated into our proposed solution algorithm to generate paths in evaluating the equilibrium flows. In the multimodal transport network, we first examine the effect of the proposed RPR scheme with the CMSTA model and access the emissions in the road subnetwork. 5 links in the CBD district are subsequently identified as the critical links, and be further grouped into the link set . As shown in Figure 4, the provision of RPR facilities significantly influence travellers’ mode choice as well as route choice behaviour. 14.03% of previous autodrivers and 5.51% of previous train riders shift to RPR mode, due to the utility provided by RPR services. The amount of emissions on critical links also significantly reduced with the provision of RPR schemes. This phenomenon indicates RPR schemes could alleviate the environmental deterioration by redistributing the flow and modal demand in the network and reflects the rationality of modelling users modal choice and route choice decisions in a combined model, i.e., CMSTA.
We proceed to exert environmental constraints on these critical links. Here, we consider five scenarios of the thresholds and summarize in Table 4. Parameters in the SAGP algorithm are set as follows: , , and . In order to reflect the convergence trend of the SAGP algorithm apparently, we take the following logarithmic value of the maximum absolute error as a performance index, namely,The SAGP algorithm is then used to solve the environmental constrained CMSTA in terms of the five scenarios provided in Table 4. Figure 5 depicts the logarithmic value of the maximum absolute error versus the number of iterations used for solving scenarios 1, 2, 3, and 4. It clearly shows that the algorithm can converge to a solution of VI model . In scenario 5, the solution converges much faster than other scenarios because only link 22 is constrained.

Table 5 gives the resultant link emissions and corresponding Lagrangian multiplier for each scenario. We can observe that the emissions on each link are bounded by the required thresholds (some link emissions slightly excess the requirements caused by the computation error). The corresponding multiplier takes positive value only on the links with emission bounded by the threshold. These two phenomena tally with (16)(18), which numerically verify the effectiveness of the VI model and the SAGP algorithm. In addition, from scenario 1 to scenario 5, it can be seen that, by tightening the thresholds, the corresponding Lagrangian multipliers become larger, which indicates that it would be more costly to restrict the emission to a much lower level.
 
LM denotes Lagrangian multiplier. 
Table 6 further investigates the modal shift trend from scenario 1 to scenario 5. As the values of threshold getting smaller, travellers who use automode become less. The environmental regulations in road network could significantly promote the modal shift from auto to “green” travel modes (i.e., Train and RPR). Taking scenario 1, for example, the automode share is depressed about 11.11%, from 53.21% to 42.10%, where 5.38% of travel demand shifts to rail mode and 5.72% rest shifts to RPR mode. This phenomenon indicates that environmental regulations significantly influence not only travellers’ route choice behaviour but also mode choice behaviour.

To further see the impacts of demand level on the Lagrangian multipliers, a sensitivity test is carried out on the dispersion parameter on the modal choice function (8). Without loss of generality, the first scenario of threshold setting is taken for this test, and then the value of is changed from 0.5 to 0.9 with an interval of 0.1. The Lagrangian multipliers obtained by the SAGP algorithm in these five cases are presented in Table 7. We noted that a larger value of results in a larger level of autodemand level and the road network congestion level is higher accordingly. Table 7 indicates the generalized CMSTA conditions can still be fulfilled. As the demand level ascends, the corresponding Lagrangian multipliers also increases, meaning that it is more difficult to maintain the threshold constraints at a higher demand level.
 
LM denotes Lagrangian multiplier. 
We further compare the networklevel performance before and after implementing the set of environmental restraints. The first scenario of threshold setting is still taken for this test. Without loss of generality, the total travel time (TTT) and total vehicle distances (TVD) travelled are used in this study to measure the networkwide performance. From Table 8 we can observe that implementation of environmental constraints decreases both TTT and TVD (13.89% decrease of TTT and 13.19% decrease of TVM), which implies that some travellers must be shifting their modes in order to benefit the air quality of the CBD area.

7. Conclusions
Environmental protection is an important consideration in planning a transportation system. Environmental constraint is a useful means to explicitly reflect environmental protection requirements. In this paper, we studied a combined environmental protection strategy in the multimodal transport network: a railbased parkandride scheme to assist mode shift to reduce the overall vehicle travel demands in the road network and an areabased regulation scheme to deal with the environmental deterioration caused by heterogeneity of traffic flows. To quantitatively evaluate the effectiveness of the proposed scheme, we built an environmental constrained combined modal split and traffic assignment (ECCMSTA) model. Specifically, we developed an equivalent VI reformulation for the ECCMSTA. Different from the explicit formulation as a VI or NCP, the proposed reformulation adapts the concept of fixpoint to reformulate the CMSTA conditions, which is, later, incorporated into the VI formulated nonlinear complementarity conditions associated with the nonlinear environmental constraints. The proposed VI reformulation can handle a general nonlinear environmental constraint structure and a general link and route cost structure, enhancing the modelling adaptability and flexibility, and in addition permits a number of efficient algorithms for its solution. A customized projection algorithm with selfadaptive step size scheme is adopted in this paper. Numerical examples demonstrated the ECCMSTA model has the potential to enhance the behavioural modelling of network users’ travel decisions and to assist in quantitatively evaluating the effectiveness of RPR schemes and environmental regulations.
Several directions are worthy of further investigation. First, in this paper, we use CO as an indicator to model vehicle emission. Further study will take other types of emission into consideration, e.g., the CO_{2}. Second, a number of efficient solution algorithms proposed by recent studies, e.g., Ryu, Chen, and Choi [34]; Ryu, Chen, and Choi [40], should be explored to further improve the solution efficiency. Third, we plan to test our model on larger networks with more real world applications. For very largescale problems, existing algorithms may be not efficient and specific algorithms should be developed.
Appendix
A. Proof of Proposition 5
Vector function is strictly monotone on ; namely,
Proof. For any two distinct nonnegative Lagrangian multiplier vectors and , let and denote the corresponding CMSTA link flow solutions, and similarly, let , , , and denote the corresponding CMSTA path flow and modal demand solutions. We havewhere defines the linkpath incidence, , , , , , , , , and .
According to the definition of CMSTA conditions, the modal travel time is taken as the equilibrium path travel time of mode , it follows thattogether with the flow conservation condition (5), we haveSubtracting (A.9) and (A.10) yields thatnamely,The monotonicity of function (8) implies that According to (A.13) and (A.14), we thus obtainAs and , (A.15) can be equivalently written asRearranging the left side of (A.16), it follows that According to the definition of generalized link travel time function, (20), we haveTogether with Assumption 1, is a strictly monotone function with respect to link flow solutions ; it follows thatwhich implies is strictly monotone with respect to on .
B. Proof of Proposition 6
Vector function is Lipschitz continuous on ; namely,
Proof. Following the approach presented by Meng, Liu, and Wang [28], it is easily to conclude that is continuously differentiable in a neighbourhood of according to the implicit function theorem [41].
Since is continuously differentiable on , its Jacobian matrix is thus continuous on . The twonorm of is therefore bounded from the nonempty closed and convex set ; namely, there is a positive constant such that According to the Theorem of Ortega and Rheinboldt [41] (meanvalue theorem), it can be seen thatSimilarly, according to Assumption 1, is a continuous differentiable function w.r.t. . The Jacobian matrix is thus continuous on which is also nonempty closed and convex. Thus, there is a positive constant such that and we haveAccording to (B.3) and (B.5), we thus haveThis completes the proof.
C. Proof of Proposition 7
is a solution of if and only if and fulfil the generalized CMSTA conditions (15)(18).
Proof.
Necessary Condition. Suppose is a solution of , and we now show that fulfil the generalized CMSTA conditions.
We first demonstrate that . Assume that there is at least on element with . Then, we define a feasible vector . Substituting in the VI model with vector , it follows thatnamely, Now, we proceed to demonstrate by using an apagogical approach. Note that, according to the monotony increasing of , there exists a upper bound of link flow that, for any flow , we haveThen, we proceed to show . Suppose at least one flow loaded path on the road network uses link . According to the DUE conditions, for any flow loaded paths, they have identical path travel time. Denote one of these paths as and its follows thatTherefore, its generalized path travel time and the corresponding modal travel time fulfil the condition:According to the monotonicity of function (8), we haveIn other words,However, (C.10) contradicts (C.3), then the hypothesis is incorrect. In other words,To further prove the necessary conditions, we construct the following two specific vectors for any link :where is a given parameter. Substituting vector of VI model with these two particular vectors, respectively, yields thatSince , (C.14) implies thatIn accordance with the fact that , we thus haveHence, (C.17) in conjunction with (C.15) implies thatIn other words, and fulfil the generalized CMSTA conditions.
Sufficient Condition. Suppose fulfil the generalized CMSTA conditions (15)(18); we have the following conditions:In addition, we have . Therefore, according to (C.19) and (C.20), it follows thatThis completes the proof.
D. Attributes of Study Network
See Table 9.
 
T is tail node; H is head node; Fft denotes free flow time; Capa denotes capacity. 
E. OD Travel Demand of Study Network
See Table 10.

Data Availability
No data were used to support this study.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
References
 X. Xu, A. Chen, and L. Cheng, “Reformulating environmentally constrained traffic equilibrium via a smooth gap function,” International Journal of Sustainable Transportation, vol. 9, no. 6, pp. 419–430, 2015. View at: Publisher Site  Google Scholar
 S. Kitthamkesorn, A. Chen, X. Xu, and S. Ryu, “Modeling mode and route similarities in network equilibrium problem with gogreen modes,” Networks and Spatial Economics, vol. 16, no. 1, pp. 33–60, 2016. View at: Publisher Site  Google Scholar
 W. Y. Szeto, X. Jaber, and S. C. Wong, “Road network equilibrium approaches to environmental sustainability,” Transport Reviews, vol. 32, no. 4, pp. 491–518, 2012. View at: Publisher Site  Google Scholar
 X. Chen, Z. Liu, I. Kim, and J. Leng, “Optimal Sitting and Sizing of the Remote ParkandRide Scheme in the Multimodal Transport Network,” in Proceedings of the Transportation Research Board 97th Annual Meeting, Washington DC, USA, 2018. View at: Google Scholar
 Z. Liu and Q. Meng, “Busbased parkandride system: a stochastic model on multimodal network with congestion pricing schemes,” International Journal of Systems Science, vol. 45, no. 5, pp. 994–1006, 2014. View at: Publisher Site  Google Scholar  MathSciNet
 C. E. Wallace, K. Courage, D. Reaves, Schoene. G., and G. Euler, TRANSYT7F user's manual.In.
 L. R. Rilett and C. M. Benedek, “Traffic assignment under environmental and equity objectives,” Transportation Research Record, no. 1443, pp. 92–99, 1994. View at: Google Scholar
 C. M. Benedek and L. R. Rilett, “Equitable traffic assignment with environmental cost functions,” Journal of Transportation Engineering, vol. 124, no. 1, pp. 16–22, 1998. View at: Publisher Site  Google Scholar
 S. Sugawara and D. Niemeier, “How much can vehicle emissions be reduced? Exploratory analysis of an upper boundary using an emissionsoptimized trip assignment,” Transportation Research Record Journal of the Transportation Research Board, no. 1815, pp. 29–37, 2002. View at: Google Scholar
 A. Chen, Z. Zhou, and S. Ryu, “Modeling physical and environmental side constraints in traffic equilibrium problem,” International Journal of Sustainable Transportation, vol. 5, no. 3, pp. 172–197, 2011. View at: Publisher Site  Google Scholar
 Z.C. Li, W. H. K. Lam, S. C. Wong, and A. Sumalee, “Environmentally sustainable toll design for congested road networks with uncertain demand,” International Journal of Sustainable Transportation, vol. 6, no. 3, pp. 127–155, 2012. View at: Publisher Site  Google Scholar
 A. Nagurney, P. Ramanujam, and K. K. Dhanda, “A multimodal traffic network equilibrium model with emission pollution permits: complicance vs noncompliance,” Transportation Research Part D: Transport and Environment, vol. 3, no. 5, pp. 349–374, 1998. View at: Publisher Site  Google Scholar
 E. Fernandez, J. de Cea, M. Florian, and E. Cabrera, “Network equilibrium models with combined modes,” Transportation Science, vol. 28, no. 3, pp. 182–192, 1994. View at: Publisher Site  Google Scholar
 R. García and A. Marín, “Network equilibrium with combined modes: models and solution algorithms,” Transportation Research Part B: Methodological, vol. 39, no. 3, pp. 223–254, 2005. View at: Publisher Site  Google Scholar
 Z.C. Li, W. H. K. Lam, S. C. Wong, D.L. Zhu, and H.J. Huang, “Modeling parkandride services in a multimodal transport network with elastic demand,” Transportation Research Record: Journal of the Transportation Research Board, no. 1994, pp. 101–109, 2007. View at: Google Scholar
 W. H. Lam, Z. Li, S. C. Wong, and D. Zhu, “Modeling an elasticdemand bimodal transport network with parkandride trips,” Tsinghua Science and Technology, vol. 12, no. 2, pp. 158–166, 2007. View at: Publisher Site  Google Scholar  MathSciNet
 C. Pineda, C. E. Cortés, P. JaraMoroni, and E. Moreno, “Integrated traffictransit stochastic equilibrium model with parkandride facilities,” Transportation Research Part C: Emerging Technologies, vol. 71, pp. 86–107, 2016. View at: Publisher Site  Google Scholar
 Z. Song, Y. He, and L. Zhang, “Integrated planning of parkandride facilities and transit service,” Transportation Research Part C: Emerging Technologies, vol. 74, pp. 182–195, 2017. View at: Publisher Site  Google Scholar
 S. Kitthamkesorn and A. Chen, “Alternate weibitbased model for assessing green transport systems with combined mode and route travel choices,” Transportation Research Part B: Methodological, vol. 103, pp. 291–310, 2017. View at: Publisher Site  Google Scholar
 Z. Liu, X. Chen, Q. Meng, and I. Kim, “Remote parkandride network equilibrium model and its applications,” Transportation Research Part B: Methodological, vol. 117, pp. 37–62, 2018. View at: Publisher Site  Google Scholar
 S. Wang, Q. Meng, and H. Yang, “Global optimization methods for the discrete network design problem,” Transportation Research Part B: Methodological, vol. 50, pp. 42–60, 2013. View at: Publisher Site  Google Scholar
 Z. Liu, S. Wang, and Q. Meng, “Optimal joint distance and time toll for cordonbased congestion pricing,” Transportation Research Part B: Methodological, vol. 69, pp. 81–97, 2014. View at: Publisher Site  Google Scholar
 T. Liu, H. Huang, H. Yang, and X. Zhang, “Continuum modeling of parkandride services in a linear monocentric city with deterministic mode choice,” Transportation Research Part B: Methodological, vol. 43, no. 6, pp. 692–707, 2009. View at: Publisher Site  Google Scholar
 D. Z. W. Wang and B. Du, “Reliabilitybased modeling of parkandride service on linear travel corridor,” Transportation Research Record: Journal of the Transportation Research Board, vol. 2333, pp. 16–26, 2013. View at: Publisher Site  Google Scholar
 B. Du and D. Z. W. Wang, “Continuum modeling of parkandride services considering travel time reliability and heterogeneous commuters – A linear complementarity system approach,” Transportation Research Part E: Logistics and Transportation Review, vol. 71, no. 0, pp. 58–81, 2014. View at: Publisher Site  Google Scholar
 Z. Liu, Q. Meng, and S. Wang, “Variational inequality model for cordonbased congestion pricing under side constrained stochastic user equilibrium conditions,” Transportmetrica A: Transport Science, vol. 10, no. 8, pp. 693–704, 2014. View at: Publisher Site  Google Scholar
 T. Larsson and M. Patriksson, “Side constrained traffic equilibrium models— analysis, computation and applications,” Transportation Research Part B: Methodological, vol. 33, no. 4, pp. 233–264, 1999. View at: Publisher Site  Google Scholar
 Q. Meng, Z. Liu, and S. Wang, “Asymmetric stochastic user equilibrium problem with elastic demand and link capacity constraints,” Transportmetrica A: Transport Science, vol. 10, no. 4, pp. 304–326, 2014. View at: Publisher Site  Google Scholar
 C. F. Daganzo, “Stochastic network equilibrium with multiple vehicle types and asymmetric, indefinite link cost Jacobians,” Transportation Science, vol. 17, no. 3, pp. 282–300, 1983. View at: Publisher Site  Google Scholar  MathSciNet
 G. E. Cantarella, “A general fixedpoint approach to multimode multiuser equilibrium assignment with elastic demand,” Transportation Science, vol. 31, no. 2, pp. 107–128, 1997. View at: Publisher Site  Google Scholar
 A. Chen, X. Xu, S. Ryu, and Z. Zhou, “A selfadaptive Armijo stepsize strategy with application to traffic assignment models and algorithms,” Transportmetrica A: Transport Science, vol. 9, no. 8, pp. 695–712, 2011. View at: Publisher Site  Google Scholar
 A. Chen, Z. Zhou, and X. Xu, “A selfadaptive gradient projection algorithm for the nonadditive traffic equilibrium problem,” Computers & Operations Research, vol. 39, no. 2, pp. 127–138, 2012. View at: Publisher Site  Google Scholar
 Z. Zhou, A. Chen, and S. C. Wong, “Alternative formulations of a combined trip generation, trip distribution, modal split, and trip assignment model,” European Journal of Operational Research, vol. 198, no. 1, pp. 129–138, 2009. View at: Publisher Site  Google Scholar
 S. Ryu, A. Chen, and K. Choi, “Solving the combined modal split and traffic assignment problem with two types of transit impedance function,” European Journal of Operational Research, vol. 257, no. 3, pp. 870–880, 2017. View at: Publisher Site  Google Scholar  MathSciNet
 B. S. He, H. Yang, Q. Meng, and D. R. Han, “Modified GoldsteinLevitinPolyak projection method for asymmetric strongly monotone variational inequalities,” Journal of Optimization Theory and Applications, vol. 112, no. 1, pp. 129–143, 2002. View at: Publisher Site  Google Scholar  MathSciNet
 D. Han and W. Sun, “A new modified GoldsteinLevitinPolyakprojection method for variational inequality problems,” Computers Mathematics with Applications, vol. 47, no. 12, pp. 1817–1825, 2004. View at: Google Scholar
 Y. Yin and S. Lawphongpanich, “Internalizing emission externality on road networks,” Transportation Research Part D: Transport and Environment, vol. 11, no. 4, pp. 292–301, 2006. View at: Publisher Site  Google Scholar
 A. Nagurney, Q. Qiang, and L. S. Nagurney, “Environmental impact assessment of transportation networks with degradable links in an era of climate change,” International Journal of Sustainable Transportation, vol. 4, no. 3, pp. 154–171, 2010. View at: Publisher Site  Google Scholar
 X. Yang, X. J. Ban, and R. Ma, “Mixed equilibria with common constraints on transportation networks,” Networks and Spatial Economics, vol. 17, no. 2, pp. 547–579, 2017. View at: Publisher Site  Google Scholar  MathSciNet
 S. Ryu, A. Chen, and K. Choi, “Solving the stochastic multiclass traffic assignment problem with asymmetric interactions, route overlapping, and vehicle restrictions,” Journal of Advanced Transportation, vol. 50, no. 2, pp. 255–270, 2016. View at: Publisher Site  Google Scholar
 J. M. Ortega and W. C. Rheinboldt, Iterative Solution of Nonlinear Equations in Several Variables, Academic Press, Thailand, 1970. View at: MathSciNet
Copyright
Copyright © 2018 Xinyuan Chen and Inhi Kim. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.