Information Sharing, Emission, and Safety in Sustainable Freight Transportation
View this Special IssueResearch Article  Open Access
Multimodal Green Logistics Network Design of Urban Agglomeration with Stochastic Demand
Abstract
This study investigates a multimodal green logistics network design problem of urban agglomeration with stochastic demand, in which different logistics authorities among the different cities jointly optimize the logistics node configurations and uniform carbon taxes over logistics transport modes to maximize the total social welfare of urban agglomeration and consider logistics users’ choice behaviors. The users’ choice behaviors are captured by a logitbased stochastic equilibrium model. To describe the game behaviors of logistics authorities in urban agglomeration, the problem is formulated as two nonlinear bilevel programming models, namely, independent and centralized decision models. Next, a quantumbehaved particle swarm optimization (QPSO) embedded with a Method of Successive Averages (MSA) is presented to solve the proposed models. Simulation results show that to achieve the overall optimization layout of the green logistics network in urban agglomeration the logistics authorities should adopt centralized decisions, construct a multimode logistics network, and make a reasonable carbon tax.
1. Introduction
The logistics network design problem (LNDP) has been widely studied [1–3]. Generally, this problem comprises two subproblems: (i) a location problem, namely, how to decide the locations of logistics nodes (such as logistics parks, distribution centers, and logistics terminals), and (ii) an allocation problem, namely, how to route the flow of goods to origin–destination (OD). Since the construction of urban logistics nodes has the characteristics of a large investment, long period, and great risk, the rationality and feasibility of the plan need to be ensured with the objective guidance of theory and method. Meanwhile, with the development of society and the improvement of people’s living standards, traffic congestion and transportationrelated environmental issues have become the focus of attention for scholars and governments [4–6]. Against this unsettling backdrop, the traditional single mode of transport does not meet the needs of modern cities. Therefore, how to design a multimodal logistics network with high efficiency, safety, and environmental friendliness is an important issue to be solved urgently in the development of urban agglomeration logistics.
To locate logistics nodes, as addressed in the conventional facility location problem (FLP) [7], the multimodal logistics network design problem (MLNDP) of urban agglomeration needs to address the rational allocation of limited resources between different cities in urban agglomerations. Because of the limited resources in urban agglomerations, there is a great game among each local government in the resources (e.g., politics, economy, and population), and this game will derive many gamedisordered problems under the restriction of the government’s political and economic behavior. Local governments usually establish the logistics network based on the situation of their own precinct to fit the regional development with little coordination and communication with other surrounding regions, resulting in the lack of effective integration of the basic logistics resources [8]. Therefore, game theory is used to analyze the game relationship among city governments in urban agglomeration, which is of great significance to the rational allocation of resources.
To route the flow of goods, multimodal transportation is a hybrid of transportation that involves two classic transportation services: the single mode of transport and combined transport [9]. Compared with the single transportation mode, the multimodal transportation can make full use of existing logistics networks of infrastructure and the advantage of various modes of transport, integrating transport capacity resources, meeting the transport demand, and achieving the target (e.g., time, cost, and profit) [10, 11]. In the multimodal transportation, there is no effective connection between each mode of transportation, and the design of the multimodal transportation plan needs to be further improved and optimized [12]. However, urban traffic congestion and transportationrelated environmental issues (e.g., carbon dioxide (CO_{2}), nitrogen oxide, and sulfur oxide) are the bottlenecks that restrict the sustainable development of cities [13, 14]. Consequently, there is an inevitable trend toward the future development of urban agglomeration green logistics systems to develop multimodal and integrative transportation networks.
1.1. Literature Review
We classify the literature related to our research into three branches: (1) urban agglomeration logistics, (2) green logistics network design for facility location, and (3) logistics network design with choice behaviors.
1.1.1. Urban Agglomeration Logistics
The urban agglomeration logistics system is an important part of the urban agglomeration economy as it plays an important and supporting role in the sustainable development and evolution of the urban agglomeration. Mullen and Marsden [15] explored the adaptive adjustment mechanism of regional logistics and regional economy. Li et al. [16] emphasized the relationship between traffic investment and traffic efficiency in urban agglomeration. Jiang et al. [17] analyzed the bidirectional relationship between multimodal transportation investment and economic development. However, some scholars studied the spatial and temporal distribution characteristics of urban agglomeration logistics [18, 19]. Lindsey et al. [20] studied the relationships among freight transport, economic market drivers, and industrial space demand. Kumar et al. [21] investigated spatial patterns of transportation and logistics cluster in the US regions, applying spatial cluster and econometric analyses. The above studies mainly focus on the linkage mechanism between logistics and economy, demand and spatial and temporal distribution, and few focus on the green logistics system of urban agglomeration. Table 1 summarizes the main features of the proposed models in the reviewed studies.
 
S: single, M: multiple, De: deterministic, St: stochastic, UE: user equilibrium, SUE: stochastic user equilibrium, GA: Genetic Algorithm, and MOGA: Multiobjective Genetic Algorithm. SEAMO2: Simple Evolutionary Multiobjective Optimization 2, FLAA: Facility Linear Approximation Algorithm, GLS: Genetic Local Search, and B&B: BranchandBound. FW: Frank–Wolfe Algorithm, QPSO: Quantum behaved Particle Swarm Optimization, and MSA: Method of Successive Averages. 
1.1.2. Green Logistics Network Design for Facility Location
Green logistics aims to restrain the environmental damage (e.g., greenhouse gas emissions, noise, and accidents) caused by logistics activity and develop a sustainable balance between economic, social, and environmental objectives [22]. Facility location modeling is a strategic planning design approach that selects the optimal set of facilities from a set of potential facilities [23]. As a branch of network design, the green logistics network design problem originates from the classic facility location problem (FLP) [7, 24] and then extends it to the locationallocation problem (LAP) [25–27].
In the background of the vigorously developing green economy, Turken et al. [24] investigated the impact of carbon tax on plant capacity and location decisions of a firm. Yang et al. [26] studied the lowcarbon network design problem of the thirdparty logistics distribution system. To solve locationallocation problem with a carbon emissions constraint, Rezaee et al. [2] developed a twostage stochastic programming model, and Gao et al. [28] proposed five multidimensional mixedinteger nonlinear programming models. In practice, the urban agglomeration logistics network may comprise multiple urban logistics networks, which are independently and separately managed by local governments (or logistics authorities) with different targets. Nevertheless, the literature always supposed that the network is totally designed by a single management body, which may not suit the design of the logistics network in urban agglomeration.
1.1.3. Green Logistics Network Design under Route Choice Models
In general, route choice models should evaluate the utility of each route and route the flow of passengers or goods from origin to destination. Although they have been widely used to solve urban passenger transportation problems by researchers and engineers [29–33], using these models to solve urban cargo transportation problems has received little attention in the literature.
The logistics network design problem based on route choice behavior can be characterized as a bilevel programming model or an equivalent mathematical program with equilibrium constraints (MPEC) model. At the upper level, the logistics authority determines the number, location, and capacity of logistics facilities (nodes). At the lower level, shippers and carriers follow the user equilibrium (UE) principle or stochastic user equilibrium (SUE) principle [34–36]. Yamada et al. [37] presented a bilevel programming model to simulate the multimodal freight transport network design problem, where the lowerlevel model is a multimodal multiclass user equilibrium model. To address the intermodal hubandspoke network design problem for multiple stakeholders and multitype containers, Meng and Wang [38] established an MPEC model and used a variational inequality to describe the operators’ route choice behavior. Wang and Meng [39] extended the research work of Meng and Wang [38] by considering congestion effects and piecewise linear cost functions and solved it with a solution algorithm based on nonlinear optimization and branchandbound global optimization.
However, the route choice models of these studies did not consider the carbon emissions cost, perception error, and elastic demand. In addition, they mainly focus on static, determined network design issues. As strategic decision planning, the design of a logistics network in urban agglomeration involves longterm planning of a between 5 and 20year timeline, where the logistics demand is often uncertain with the development of urban evolution and industrial structure optimization. Hence, it is more realistic and reasonable to design the logistics network with demand uncertainty and stochastic route choice.
1.2. Objectives and Contributions
In view of the abovementioned realistic problems and specific characteristics of logistics development in urban agglomerations, the objective of this paper is to coordinate the distribution of various types of logistics nodes (e.g., logistics parks, distribution centers and logistics terminals) in the urban agglomeration logistics network from the overall perspective of urban agglomeration and fully consider the low carbon requirement together with demand uncertainties. The main contributions of this paper are as follows. First, to characterize different decisionmaking behaviors among multiple local authorities, the multimodal logistics network design problem is formulated as two nonlinear bilevel programming models. At the upper level, the logistics authority of each city attempts to maximize the social welfare including producer surplus and consumer surplus. At the lower level, the logistics users’ route choice decisions follow logitbased stochastic user equilibrium (SUE) with elastic demand under logistics demand scenario. Second, a heuristic solution algorithm that is a combination of quantum behaved particle swarm optimization (QPSO) and Method of Successive Averages (MSA) is developed to solve the proposed bilevel programming model. Third, the optimal number, scale, and location of logistics nodes, analysis of the impact on the distribution of city group logistics network of infrastructure investment budget, and carbon tax are illustrated with an example.
The remainder of this article is arranged as follows. The basic considerations of this paper, including general assumptions and network representation, are described in Section 2. A nonlinear bilevel programming formulation is proposed in Section 3, and the solution methods are presented in Section 4. A numerical example is presented to illustrate the availability of previous models in Section 5. Finally, the conclusion and future studies are discussed in Section 6.
2. Basic Considerations
2.1. General Assumptions
To facilitate the presentation of essential ideas without the loss of generality, the following basic assumptions are made.
A1. For the simplicity of expression, the urban agglomeration is assumed to be a city that is a set of all single cities, and the planning period is assumed to be one week.
A2. In the urban logistics system, the logistics nodes investment and subsidy are determined by the logistics authority.
A3. The disutility of each service route is measured by transport time, transport cost, and CO_{2} emission taxes (if any). Logistics users select their logistics service routes which are associated with their own perceptions of service disutility.
A4. The urban agglomeration comprises several cities, and each city has one logistics authority. In decentralized decisionmaking, the decisionmaking of logistics departments in urban agglomeration is completely independent.
2.2. Network Representation
To model logistics service of the urban agglomeration, we first represent the demand network, multimodal logistics physical network, and multimodal logistics service network.
Demand Network. Let be the set of logistics demand origin nodes and let be the set of logistics demand destination nodes, where is the set of logistics nodes (or transfer nodes) including existing and potential logistics parks, distribution centers, and logistics terminals. Denote by the set of logistics demand origindestination (O–D) pairs. For a given logistics OD pair, Figure 1(a) shows the different types of logistics demands such as industrial demand (K1), commercial demand (K2), and agricultural demand (K3). These demands are served by the multimodal logistics network, as shown in Figures 1(b)–1(e).
(a) Logistics user demand network
(b) Multimodal logistics physical network
(c) Representation of the virtual subnetwork of logistics transfer node 1
(d) Multimodal logistics service supernetwork
(e) Logistics service route
Multimodal Logistics Physical Network. Suppose that a multimodal logistics physical network needs to be designed, which is made up of a set of logistics nodes and a set of logistics links or arcs of different transport modes. Each logistics node or link is provided by a logistics operator.
We denote the multimodal logistics physical network as a directed network , where is the set of logistics nodes and is the set of logistics links. Let denote the set of logistics transfer nodes, and let be the set of modes. Each link is represented by a triplet , where are the starting point and ending point of link , respectively, and is the transportation mode on the link.
Analogously, for any logistics transfer node , let denote the set of logistics links pointing into transfer node , and let denote the set of logistics links stemming out of transfer node . For each logistics link , we make a copy of transfer node as its ending point. Let denote the set of heads (copies) of all logistics arcs in after copying. Similarly, For each logistics link , we make a copy of transfer node as its ending point. Let denote the set of tails (copies) of all logistics arcs in after copying. Thus, the set of virtual arcs at transfer node is denoted by , and the virtual subnetwork of logistics transfer node can be represented by a directed graph .
Figure 1(b) shows an example of a multimodal logistics physical network , where and We suppose and Then, the virtual subnetworks at each logistics transfer node can be generated. Figure 1(c) shows the example of the virtual subnetwork at logistics transfer node 1, where and .
Multimodal Logistics Service Network. A hypernetwork is used to construct a multimodal logistics service network, where represents the type of route transport modes, denotes the set of nodes after transfer node copies, denotes the set of logistics links, and denotes the set of transfer arcs. Based on the example of multimodal logistics physical network in Figure 1(b) and the virtual transfer subnetworks, the logistics service supernetwork is then generated and contains highway subnetworks and railway subnetworks (Figure 1(d)).
Feasible Sets of Route. Although in the multimodal logistics service supernetwork there are some routes from the logistics demand origin to the destination, only the routes that satisfy certain conditions, such as cost, time, or CO_{2} emissions, are called feasible routes. Thus, for any logistics OD pair , let denote the set of logistics feasible service routes connecting in the multimodal logistics service supernetwork, and let be one if service route uses transport mode and zero otherwise.
3. Model Formulation
3.1. The General Stochastic Bilevel Programming Model
Stochastic bilevel programming combines the characteristics of stochastic programming and bilevel programming and introduces random scenarios to describe the uncertainties involved in the model. The most commonly used stochastic bilevel programming model is the expected bilevel programming model, which is described in the following mathematical form [40]:where response function is implicitly defined by
Obviously, the expected bilevel programming model comprises two submodels, (U0), which is defined as an upper level and (L0), which is a lower level. , and denote the decision vector, objective function, and constraint set of the upperlevel decisionmakers or system managers, respectively. is the expectation operator with respect to random scenarios . , , and denote the decision vector, objective function, and constraint set of the lowerlevel decisionmakers or users, respectively.
Assume that the random scenarios have a finite number of states,. Let index be its possible realizations and let be their respective probabilities. Then, we can now express the extensive form of the stochastic bilevel program as follows:where response function is implicitly defined by
3.2. Decentralized Decision in Multiple Cities
Since the urban logistics system is exclusively managed by a logistics authority but serves for all logistics users (i.e., shippers and carriers) equally, the urban logistics network design is typically formulated as a Stackelberg game that can well characterize the interactions between the logistics authority and users.
As shown in Figure 2, a bilevel program is used to model the leaderfollower behaviors between logistics authority and users. At the upper level, the logistics authority attempts to maximize the social welfare by planning the investment capacity of logistics nodes and determining the carbon tax on transportation service links. At the lower level, the logistics users’ reactions and choice decisions to the urban logistics network design scheme will be assumed to follow logitbased stochastic user equilibrium (SUE) under the logistics demand scenario. Moreover, the main notations are provided in Appendix A.
3.2.1. UpperLevel Model of Green Logistics Network Design for a Single City
As the upperlevel decision maker, the logistics authority aims to maximize the total social welfare of its own city. It is well known that the total social welfare comprises the consumer surplus and producer surplus in scenario , which can be expressed aswhere
Eq. (6) formulates the consumer surplus, and is the inverse function of the logistics demand function. Consumer surplus is the extra benefit logistics users’ gains when the costs they actually pay are less than what they would be prepared to pay. Eqs. (7) and (8) formulate the producer surplus of all logistics nodes and logistics arcs, respectively. Producer surplus is a measure of producer welfare. Taking (7) as an example, the first term is transfer profits, the second term is subsidy revenues, and the last is park construction costs. Then, the upperlevel model of green logistics network design in single city is formulated assubject towhere can be obtained by solving the upperlevel model.
The objective function (9) formulates the expected social welfares. Constraint (10) represents the establishment or expansion capacity restraint of the logistics nodes; constraint (11) denotes the constraint of carbon tax; and constraint (12) represents the construction investment constraint of the logistics nodes.
3.2.2. LowerLevel Model of Green Logistics Network Design
Travel/Transfer Time. To capture the difference in attributes of different modes of transport, for each transport link , we consider the following service time function:where is the link freeflow transport service time, and is the average transport time interval. For HGVs or LGVs, the Bureau of Public Roadstype (BPR) function can be adopted to estimate transport service time. For railways or waterways, we consider service time function as a function of link freeflow transport service time and departure interval time [11, 35]. Similarly, for each virtual transfer arc, we use the following service time function:where is the arc freeflow transfer service time and and are impedance parameters.
Route Utility and Flow. According to A3, each logistics service route is associated with a given actual cost (disutility), which can be expressed aswhere is the value of time. , , and represent the transportation cost, logistics service time, and CO_{2} emission cost on service route between OD pair in city under logistics demand scenarios , respectively, which are expressed as
It is worth noting that the transportation cost (time) of each route includes the transport cost (time) on links and the transfer cost (time) at parks. Due to variations in perception, the route service disutility is perceived differently by each logistics user, and thus, the perceived disutility of each route is treated as a random variable. If the random variable can be considered to obey Gumbel distribution [41], then the path flow on route between OD pair in single city and scenario can be given bywhere represents the sensitivity of route selection disutility and denotes the logistics demand function. Eq. (19) is the most widely used flow assignment method in traffic planning and represents the flow assigned to each feasible route for each logistic demand [41]. To capture the logistics users’ responses to logistics service disutility, we assume that the elastic demand function between a generic OD pair is a monotonically decreasing function of the OD service cost between this OD pair. The generic elastic demand function [42] is expressed as follows:where represents the sensitivity to the expected service cost and is the expected minimum perceived service cost between OD pair in single city and scenario . In the case of stochastic user equilibrium (SUE) assignment with elastic demand [42], the expected minimum service cost between an OD pair under logistics demand scenario in city could be expressed as
Equivalence Model. According to the above analysis, for any given logistics authority decision , the lowerlevel model of green logistics network design in a single city can be formulated as the following equivalent minimization program:subject towhereConstraint (23) is the flow conservation constraint. Constraint (24) is the nonnegativity constraints of service route flows. Constraint (25) defines the relationship between the route flow and path flow. Constraint (26) defines the relationship between the route flow and logistics nodes flow, respectively.
Proposition 1. The minimization program (22)(26) is equivalent to the logitbased stochastic user equilibrium (SUE) assignment with elastic demand under logistics demand scenario in city . The proof of Proposition 1 is given in Appendix B.
3.3. Centralized Decision in Urban Agglomeration
In the previous subsection, we study the problem that several logistics authorities perform logistics network design decisions independently. However, such independence might be adverse to the overall benefit of urban agglomeration, since it does not explore the full potential of the investment to improve the need of all logistics users on the entire logistics network. We propose a centralized decision model in which logistics authorities cooperate with each other to maximize the total social welfare in urban agglomeration.
The upperlevel model of green logistics network design in urban agglomeration is formulated assubject to
The lower level program is still the traffic assignment problem described in Section 3.2.2.
4. Solution Methods
4.1. The Quantum Behaved Particle Swarm Optimization (QPSO)
The QPSO was introduced by Sun et al. [43], combining quantum theory based on the particle swarm optimization (PSO). Different from that in PSO, particles in QPSO have no velocity vectors; in addition, parameters that need to be adjusted are far fewer. It has been widely used in theory and engineering practice because of its global convergence together with comparative simplicity. In the QPSO, the particles are updated with the four following equations:where the denotes the particle’s position. denotes the mean best position of all the particles’ best positions. The and are the particle’s personal best position and the global best position, respectively. is the expansioncontraction coefficient, where is the current iteration number and is the maximum number of iterations; in general, , .
4.2. Hybrid QPSO Algorithm for Bilevel Decision Problem
To penalize the candidate solutions violating construction investment constraints (12), we first define the evaluation function as the sum of the objective function and penalty terms as follows:subject to (10),(11),(13)(26), and
where is the positive variable penalty coefficient.
The above model can then be solved by a solution algorithm based on the QPSO and Method of Successive Averages (MSA) hybrid algorithm (Figure 3). The detailed process of the above hybrid algorithm is described as follows.
Step 1. External Loop Initialization.(1)Set the QPSO parameters: population size () and related parameters().(2)Set the termination criterion for the external loop operation: maximum iteration (T).(3)Set the external loop iteration counter and randomly generate particle’s initial position , in which is sampled randomly in the feasible space. If , let the personal best (pbest) position of each particle
Step 2. Inner Loop Operation (Solve the LowerLevel Flow Assignment).
For a given logistics authority decision , use the method of successive averages (MSA) to solve the lowerlevel flow assignment problem for each scenario and obtain optimal solutions
Step 3. Calculate the mean of the personal best (pbest) positions of all particles by using (28).
Step 4. Update the expansioncontraction coefficient by using (33).
Step 5. For each particle, calculate the personal fitness value and pbest position .
If , update
Step 6. Update the global best (gbest) position
Step 7. For each particle, generate the new local attractor point by using (32).
Step 8. For each particle, update the particle’s position by using (34).
Step 9. Check the termination criterion for the external loop operation. If (), then stop the calculation and go to Step 10; otherwise, let , and go to Step 2.
Step 10. Output the final optimal solutions.
In the upper model, as shown in Figure 4, each particle is represented as one string comprising continuous numbers, where is the number of entities vector ; that is, , where is a set of consecutive positive integers.
The method of successive averages used in Step 2 is described in the following substeps.
Step 2.1. Inner Loop Initialization.(1)Set the termination criterion for the external loop operation ().(2)Set the initial link flow and calculate the route flow by using (19).(3)Determine effective service route set (4)Set the inner loop iteration counter ().
Step 2.2. Calculate the link flow and node flow by using (25) and (26), respectively.
Step 2.3. Calculate the link service time and node service time by using (13) and (14), respectively.
Step 2.4. Calculate the disutility by using (15).
Step 2.5. Calculate the auxiliary route flow by using (19).
Step 2.6. Update route flow by using method of successive average (MSA).
Let
Step 2.7. Check the termination criterion for the inner loop operation.
Letif , then stop the calculation and go to Step 2.8; otherwise, let and go to Step 2.2.
Step 2.8. Return to the inner loop optimal solutions
5. Numerical Examples
As shown in Figure 5, an example multimodal logistics network of urban agglomeration is used to illustrate the availability of the proposed bilevel model and QPSO algorithm. The solution procedure is coded by MATLAB R2012a and runs on a desktop Lenovo G3240 with an Intel Pentium 3.10 GHz and 4 GB RAM.
5.1. Main Data and Parameter Values
The network comprises 19 nodes, 42 arcs, and six logistics demand OD pairs. The logistics demand for each OD pair, the basic data of each logistics node, and each arc in the network are provided in Appendix C. To examine the effect of logistics demand uncertainty on logistics network design decisions, we consider three logistics demand scenarios. The total construction budget for urban agglomeration is 15,000 $/week. The sensitivity parameters and are 0.8 and 0.001, respectively. The value of time is 8 $/h. The unit investment subsidy for each logistics node is 0.5 $/week. The average emissions of expressway, railway, and waterway are 0.132, 0.022, and 0.016 kg/tonkm, respectively. The maximum carbon tax on arc is 1 $/kg. Unless otherwise specified, these input data are considered unchanged in the following analysis.
5.2. Numerical Results and Discussions
For comparison, the donothing model and the decentralized and centralized decision models are considered in this paper. In the donothing model and decentralized decision model, the logistics nodes budget for each single city is set as 0 and 5,000 $/week, respectively. However, in the centralized decision model, the budget is designed under a total budget constraint, 15,000 $/week. The optimal solutions and performance comparison for these three models are given in Tables 2 and 3, respectively.


Table 2 clearly shows that under the given investment budget five logistics parks located at (X1X4, X6) are constructed in the centralized decision model, and the total processing capacity is 3,658 tons/week. While the number of logistics parks increases to six (X1X6) in the decentralized decision model, the total processing capacity decreases to 3,552 tons/week. By comparing the centralized decision model with the decentralized decision model, it shows that in the decentralized decision model the investment budget for the city 1 is reasonable, while the investment budget for the city 2 is insufficient, and the investment budget for the city 3 is in surplus. Similarly, as shown in Table 3, among the three design models, the centralized decision model has the best performance, followed by the decentralized decision model. Compared with the donothing model, the centralized decision model and the decentralized decision model can reduce the average carbon emission rate by approximately 45% after the construction of a multimode network.
Given the above, to achieve overall optimality in the layout of urban agglomeration logistics networks, each logistics authority in the urban agglomeration should make joint decisions, allocate the investment budget, and plan the quantity and scale of the logistics nodes rationally. At the same time, the carbon emissions generated during the logistics and transportation can be greatly reduced by constructing a multimodal network and setting a reasonable carbon tax.
5.3. Impact Analysis of Budget Allocation on Network Performance
As seen from Table 2, the investment budget of city one is reasonable. To analyze the impact of budget distribution on network performance in the decentralized decision model, we set $/week and carried out numerical experiments of nine budget allocation scenarios (with the same total investment budget, $/week). Numerical results are presented in Table 4 and Figure 6.

From Table 4, first, we can see that the construction of the logistics park takes priority over the general logistics node. Second, with the increase in the investment budget in city 2, the number of logistics parks and logistics nodes will increase, and processing capacity will be enhanced gradually. To be specific, when the investment budget is equal to 1,0002,000$/week, one logistics park is constructed; when the investment budget is equal to 3,0008,000$/week, two logistics parks are constructed; and when the investment budget is more than 9,000$/week, two logistics parks and one logistics node are constructed. Finally, the optimal carbon tax gradually declines with the increase in the investment volume. This is because before the construction of the logistics park, there is a single highway with high carbon emissions in the logistics network, and after the construction of the logistics park there will be a single highway and multimodal transportation such as highwayrailway and highwaywaterway in the logistics network, which effectively reduces the carbon tax in the network. As shown in Figure 6, the social welfare increases with the increase in the investment budget for a single city. However, from the perspective of total social welfare, when the investment budget is allocated in the scenario of 5,0006,0004,000 $/week, the total social welfare is the maximum for the urban agglomeration.
5.4. Impact Analysis of Mode Choices on Network Performance
To fully explore the impact of mode choices on network performance, four mode combination strategies are considered: (I) Expressway, (II) ExpresswayRailway, (III) ExpresswayWaterway, and (IV) ExpresswayRailwayWaterway. The first strategy is to simulate the traditional singleroad transport structure. Considering the incomplete accessibility of other transportation modes, this paper reasonably constructs the latter three multimodal strategies based on road transportation. The network performance of different mode strategies under different budget allocation scenarios is shown in Table 5.

As can be seen from Table 5, in all budget allocation scenarios, the total social welfares under multimodal transport strategies II, III, and IV are significantly better than the traditional singleroad strategy I, with an average increase of 42.88%, 41.85%, and 53.97%, respectively. In this context, the total social welfare under the strategy IV is better than the strategies II and III, with an average increase of 7.81% and 8.58%, respectively. Furthermore, strategy III is slightly inferior to strategy II in most scenarios. To reduce the impact on the environment, the logistics authority is also concerned about the average carbon emission rate of the urban agglomeration. Table 5 also shows that the average carbon emission rate under different multimodal transport strategies can effectively reduce CO_{2} emissions. Compared with the strategy I, adopting multimodal transport strategy II, III, and IV results in an average decrease of 61.23%, 55.61%, and 74.41%, respectively. Therefore, the development of multimodal transport network can not only improve social welfare but also effectively reduce CO_{2} emissions. In addition, other realistic conditions of each strategy should be considered when selecting a multimodal strategy. For example, taking strategy I as the benchmark, there is an expected target to reduce the average carbon emission rate in the regional multimodal network design. If the target is set to 50%, all multimodal strategies are feasible. If set to 60%, strategies II and IV are feasible. If set to 70%, only strategy IV is feasible.
6. Conclusion and Future Studies
This paper studied a multimodal green logistics network design problem of urban agglomeration with four marked features: stochastic demand, congestion effects, game behaviors of multiple logistics authorities, and route choice behaviors of logistics users. To describe the game behaviors of logistics authorities in urban agglomeration, two nonlinear bilevel programming models were proposed, decentralized, and centralized decision models. Compared with the centralized decision model, the decentralized decision model is mainly used to analyze the impact of budget allocation of urban agglomeration. At the upper level, each logistics authority aims to maximize the total social welfare of its own city by planning the investment capacity of logistics nodes and determining the carbon tax on transportation service links. The lowerlevel subproblem is a logitbased stochastic user equilibrium (SUE) problem, regarding the logistics demand scenario and carbon tax. Next, to solve the nonlinear bilevel programming models, a quantumbehaved particle swarm optimization (QPSO) embedded with a Method of Successive Averages (MSA) for stochastic user equilibrium flow assignment was proposed.
Numerical examples show that centralized decisions are the optimal decisions for improving the entire logistics network performance, followed by decentralized decisions among multiple logistics authorities under a budget allocation constraint. Meanwhile, compared with the single mode logistics network, the constructed multimode logistics network can reduce the average carbon emission rate by approximately 45%. In addition, by adopting centralized decisions, constructing a multimode logistics network and formulating a reasonable carbon tax, the logistics authority can achieve the overall optimization layout of the green logistics network in urban agglomeration. In short, the proposed models and algorithms can help the logistics authorities make scientific decisions on the multimodal green logistics network design problem of urban agglomeration.
On the basis of this study, some main extensions can be made in the future, such as applying our proposed models and algorithms to a large and realistic logistics network and establishing multiobjective optimization model or robust optimization model.
Appendix
A. Notations
Sets set of all cities in the urban agglomeration set of all cities including the urban agglomeration set of logistics transfer nodes in city , where : set of candidate logistics node in city , where set of arcs of transport service in city , where : set of all logistics transferring service in city , where set of arcs of logistics service in city , where set of all link transport modes in urban logistics service market, where “1”, “2,” and “3” represent the expressway, railway, and waterway, respectively set of all route transport modes in urban logistics service market, where the first three items represent single mode and the last item represents combined mode : set of all OD pairs in the logistics network in city , where : set of all service paths between OD pair in city : set of all logistics demand scenarios.
Indices single city or the urban agglomeration index or logistics node or candidate logistics node index or logistics transport or transferring service index , , or link or route transport mode index or OD pair index service route index logistics demand scenarios index .
General Variables : logistics demand function of city between OD pair in scenario (tons/week) : potential logistics demand between OD pair in city and scenario (tons/week) : transport time function of city on arc by link transport mode in scenario (h) : disutility function on route between OD pair in city and scenario ($) : indicator variable equals 1 if service route uses transport mode : indicator variable equals 1 if link is on route by link transport mode and 0 otherwise : indicator variable equals 1 if a logistics transfer service is made in node along the route and 0 otherwise : expected minimal disutility between OD pair in city and scenario ($) : freight flow on route between OD pair in city and scenario (tons/week) : freight flow on logistics nodes in city and scenario (tons/week) : freight flow on logistics service arc in city and scenario (tons/week).
Continuous Decision Variables (Logistics Authority) : a vector defined as , which means a set of investment capacities for logistics node in city : a decision variable, which means the carbon tax on logistics service arc in city Constants : investment budget for city , where ($) : fare of unit turnover on arc served by link transport mode in city ($/ton km) : transfer cost on arc at logistics node in city ($/ton km) : service capacity on arc in city (tons/km) : unit fare charged at logistics node in city ($/ton) : unit construction cost (fixed cost) at logistics node in city ($/m^{2}) : operator cost of unit turnover on arc served by link transport mode in city ($/ton km) : unit transfer operating cost at logistics node in city ($/ton) : penalty cost for shortage of a unit of flow at each logistics node : average CO_{2} emission per unit turnover by link transport mode (kg/ton km) : expected CO_{2} emission for shipping per unit turnover in city (kg/ton km) : length of arc in city (km) : freeflow transport service time on arc by link transport mode in city (h) : freeflow transfer service time on arc at logistics node in city (h) : average transport time interval for link transport mode (h) : processing capacity of unit construction area (kg/m^{2}) : demand dispersion parameter in the elastic demand function : parameter for representing the perception variation of logistics users on logistics service disutility : parameter for capturing the effects of economies of scale for logistics nodes : service capacity on arc served by link transport mode in city (tons/km) : maximum establish or expansion capacity of logistics node in city : maximum carbon tax on arc in city : minimum subsidy scale of logistics node in city : unit investment subsidy of logistics node in city : probability of scenario in city
B. Proof of Proposition 1
See Proposition 1 and (22)(26).
Proof. Substituting the constraints (25) and (26) directly into the objective function (22), the Lagrangian function for problem (22) can be formulated The KuhnTucker (KT) conditions of problem (B.1) can be given as follows:where
+ and is the corresponding Lagrangian multiplier.
The KT condition (B.4) and condition (B.6) ensures that and ; then condition (B.2) can be rewritten asAs KT condition (B.4), equation (B.7) can easily be transformed to the following logit model for route choice probability:Similarly, combining KT condition (B.4) with (B.7), we can obtainKT condition (B.3) can be rewritten asThen, combining (B.9) and (B.10), we can obtain
C. Supplementary Data



Data Availability
The data used to support the findings of this study are included within the article.
Disclosure
A short version of this manuscript won the first prize in the 2018 Graduate Academic Annual Conference of the School of Transportation Engineering.
Conflicts of Interest
The authors declare no conflicts of interest.
Acknowledgments
The work that is described in this paper was supported by Fundamental Research Funds for the Central Universities of Central South University (No. 2018zzts162), National Natural Science Foundation of China (Nos. 71371181 and 71672193), and Degree Innovation Foundation of Central South University (No. 201831).
References
 R. Babazadeh, J. Razmi, M. Rabbani, and M. S. Pishvaee, “An integrated data envelopment analysis–mathematical programming approach to strategic biodiesel supply chain network design problem,” Journal of Cleaner Production, vol. 147, pp. 694–707, 2017. View at: Publisher Site  Google Scholar
 A. Rezaee, F. Dehghanian, B. Fahimnia, and B. Beamon, “Green supply chain network design with stochastic demand and carbon price,” Annals of Operations Research, pp. 1–23, 2015. View at: Publisher Site  Google Scholar
 Y. Zhao, Q. Xue, Z. Cao, and X. Zhang, “A twostage chance constrained approach with application to stochastic intermodal service network design problems,” Journal of Advanced Transportation, vol. 2018, Article ID 6051029, 18 pages, 2018. View at: Google Scholar
 M. Grote, I. Williams, J. Preston, and S. Kemp, “Including congestion effects in urban road traffic CO2 emissions modelling: Do Local Government Authorities have the right options?” Transportation Research Part D: Transport and Environment, vol. 43, pp. 95–106, 2016. View at: Publisher Site  Google Scholar
 K. Wu, Y. Chen, J. Ma, S. Bai, and X. Tang, “Traffic and emissions impact of congestion charging in the central Beijing urban area: A simulation analysis,” Transportation Research Part D: Transport and Environment, vol. 51, pp. 203–215, 2017. View at: Publisher Site  Google Scholar
 W. Wang, L. Huang, and X. Liang, “On the simulationbased reliability of complex emergency logistics networks in postaccident rescues,” International Journal of Environmental Research and Public Health, vol. 15, no. 1, p. 79, 2018. View at: Publisher Site  Google Scholar
 M. A. Efroymson and T. L. Ray, “A branchbound algorithm for plant location,” Operations Research, vol. 14, no. 3, pp. 361–368, 1966. View at: Publisher Site  Google Scholar
 Q. Sun, F. Tang, and Y. Tang, “An economic tie networkstructure analysis of urban agglomeration in the middle reaches of Changjiang River based on SNA,” Journal of Geographical Sciences, vol. 25, no. 6, pp. 739–755, 2015. View at: Publisher Site  Google Scholar
 M. Beuthe, B. Jourquin, J. Geerts, and C. Koul À Ndjang'Ha, “Freight transportation demand elasticities: a geographic multimodal transportation network analysis,” Transportation Research Part E: Logistics and Transportation Review, vol. 37, no. 4, pp. 253–266, 2001. View at: Publisher Site  Google Scholar
 M. Marufuzzaman and S. D. Ekşioğlu, “Designing a reliable and dynamic multimodal transportation network for biofuel supply chains,” Transportation Science, vol. 51, no. 2, pp. 494–517, 2017. View at: Publisher Site  Google Scholar
 D. Zhang, X. Li, Y. Huang, S. Li, and Q. Qian, “A robust optimization model for green regional logistics network design with uncertainty in future logistics demand,” Advances in Mechanical Engineering, vol. 7, no. 12, 2015. View at: Google Scholar
 G. W. Xiong and Y. Wang, “Best routes selection in multimodal networks using multiobjective genetic algorithm,” Journal of Combinatorial Optimization, vol. 28, no. 3, pp. 655–673, 2014. View at: Publisher Site  Google Scholar  MathSciNet
 Y. Xiao and A. Konak, “The heterogeneous green vehicle routing and scheduling problem with timevarying traffic congestion,” Transportation Research Part E: Logistics and Transportation Review, vol. 88, pp. 146–166, 2016. View at: Publisher Site  Google Scholar
 D. Zhang, X. Wang, S. Li, N. Ni, Z. Zhang, and X. Ma, “Joint optimization of green vehicle scheduling and routing problem with timevarying speeds,” PLoS ONE, vol. 13, no. 2, article e0192000, 2018. View at: Publisher Site  Google Scholar
 C. Mullen and G. Marsden, “Transport, economic competitiveness and competition: A city perspective,” Journal of Transport Geography, vol. 49, pp. 1–8, 2015. View at: Publisher Site  Google Scholar
 J. Li, G. Rong, and Y. Feng, “Request selection and exchange approach for carrier collaboration based on auction of a single request,” Transportation Research Part E: Logistics and Transportation Review, vol. 84, pp. 23–39, 2015. View at: Publisher Site  Google Scholar
 X. Jiang, X. He, L. Zhang, H. Qin, and F. Shao, “Multimodal transportation infrastructure investment and regional economic development: A structural equation modeling empirical analysis in China from 1986 to 2011,” Transport Policy, vol. 54, pp. 43–52, 2017. View at: Publisher Site  Google Scholar
 H. Min, C. S. Ko, and H. J. Ko, “The spatial and temporal consolidation of returned products in a closedloop supply chain network,” Computers & Industrial Engineering, vol. 51, no. 2, pp. 309–320, 2006. View at: Publisher Site  Google Scholar
 M. K. Delivand, A. R. Cammerino, P. Garofalo, and M. Monteleone, “Optimal locations of bioenergy facilities, biomass spatial availability, logistics costs and GHG (greenhouse gas) emissions: a case study on electricity productions in South Italy,” Journal of Cleaner Production, vol. 99, pp. 129–139, 2015. View at: Publisher Site  Google Scholar
 C. Lindsey, H. S. Mahmassani, M. Mullarkey, T. Nash, and S. Rothberg, “Industrial space demand and freight transportation activity: exploring the connection,” Journal of Transport Geography, vol. 37, pp. 93–101, 2014. View at: Publisher Site  Google Scholar
 I. Kumar, A. Zhalnin, A. Kim, and L. J. Beaulieu, “Transportation and logistics cluster competitive advantages in the U.S. regions: A crosssectional and spatiotemporal analysis,” Research in Transportation Economics, vol. 61, pp. 25–36, 2017. View at: Publisher Site  Google Scholar
 R. Dekker, J. Bloemhof, and I. Mallidis, “Operations Research for green logistics  An overview of aspects, issues, contributions and challenges,” European Journal of Operational Research, vol. 219, no. 3, pp. 671–679, 2012. View at: Google Scholar
 M. GhaneEzabadi and H. A. Vergara, “Decomposition approach for integrated intermodal logistics network design,” Transportation Research Part E: Logistics and Transportation Review, vol. 89, pp. 53–69, 2016. View at: Publisher Site  Google Scholar
 N. Turken, J. Carrillo, and V. Verter, “Facility location and capacity acquisition under carbon tax and emissions limits: To centralize or to decentralize?” International Journal of Production Economics, vol. 187, pp. 126–141, 2017. View at: Publisher Site  Google Scholar
 I. Harris, C. L. Mumford, and M. M. Naim, “A hybrid multiobjective approach to capacitated facility location with flexible store allocation for green logistics modeling,” Transportation Research Part E: Logistics and Transportation Review, vol. 66, pp. 1–22, 2014. View at: Publisher Site  Google Scholar
 J. Yang, J. Guo, and S. Ma, “Lowcarbon city logistics distribution network design with resource deployment,” Journal of Cleaner Production, vol. 119, pp. 223–228, 2016. View at: Publisher Site  Google Scholar
 N. Zarrinpoor, M. S. Fallahnezhad, and M. S. Pishvaee, “Design of a reliable hierarchical locationallocation model under disruptions for health service networks: A twostage robust approach,” Computers & Industrial Engineering, vol. 109, pp. 130–150, 2017. View at: Publisher Site  Google Scholar
 J. Gao, Z. Xiao, B. Cao, and Q. Chai, “Green supply chain planning considering consumer’s transportation process,” Transportation Research Part E: Logistics and Transportation Review, vol. 109, pp. 311–330, 2018. View at: Publisher Site  Google Scholar
 T. Ercan, N. C. Onat, O. Tatari, and J. Mathias, “Public transportation adoption requires a paradigm shift in urban development structure,” Journal of Cleaner Production, vol. 142, pp. 1789–1799, 2017. View at: Publisher Site  Google Scholar
 B. Peng, H. Du, S. Ma, Y. Fan, and D. C. Broadstock, “Urban passenger transport energy saving and emission reduction potential: a case study for Tianjin, China,” Energy Conversion and Management, vol. 102, pp. 4–16, 2015. View at: Publisher Site  Google Scholar
 C. Teye, M. G. Bell, and M. C. Bliemer, “Urban intermodal terminals: The entropy maximising facility location problem,” Transportation Research Part B: Methodological, vol. 100, pp. 64–81, 2017. View at: Publisher Site  Google Scholar
 H. Yang, X. Zhang, and Q. Meng, “Stackelberg games and multiple equilibrium behaviors on networks,” Transportation Research Part B: Methodological, vol. 41, no. 8, pp. 841–861, 2007. View at: Publisher Site  Google Scholar
 D. Zhang, F. Zou, S. Li, and L. Zhou, “Green supply chain network design with economies of scale and environmental concerns,” Journal of Advanced Transportation, vol. 2017, Article ID 6350562, 14 pages, 2017. View at: Google Scholar
 C. F. G. Loureiro and B. Ralston, “Investment selection model for multicommodity multimodal transportation networks,” Transportation Research Record, no. 1522, pp. 38–46, 1996. View at: Google Scholar
 D. Zhang, Q. Zhan, Y. Chen, and S. Li, “Joint optimization of logistics infrastructure investments and subsidies in a regional logistics network with CO2 emission reduction targets,” Transportation Research Part D: Transport and Environment, vol. 60, pp. 174–190, 2018. View at: Publisher Site  Google Scholar
 X. Wang, Q. Meng, and L. Miao, “Delimiting port hinterlands based on intermodal network flows: Model and algorithm,” Transportation Research Part E: Logistics and Transportation Review, vol. 88, pp. 32–51, 2016. View at: Publisher Site  Google Scholar
 T. Yamada, B. F. Russ, J. Castro, and E. Taniguchi, “Designing multimodal freight transport networks: A heuristic approach and applications,” Transportation Science, vol. 43, no. 2, pp. 129–143, 2009. View at: Publisher Site  Google Scholar
 Q. Meng and X. Wang, “Intermodal hubandspoke network design: incorporating multiple stakeholders and multitype containers,” Transportation Research Part B: Methodological, vol. 45, no. 4, pp. 724–742, 2011. View at: Publisher Site  Google Scholar
 X. Wang and Q. Meng, “Discrete intermodal freight transportation network design with route choice behavior of intermodal operators,” Transportation Research Part B: Methodological, vol. 95, pp. 76–104, 2017. View at: Publisher Site  Google Scholar
 M. Patriksson, “Robust bilevel optimization models in transportation science,” Philosophical Transactions of the Royal Society A: Mathematical, Physical & Engineering Sciences, vol. 366, no. 1872, pp. 1989–2004, 2008. View at: Publisher Site  Google Scholar  MathSciNet
 R. B. Dial, “A probabilistic multipath traffic assignment model which obviates path enumeration,” Transportation Research, vol. 5, no. 2, pp. 83–111, 1971. View at: Publisher Site  Google Scholar
 H. Yang and M. G. H. Bell, “Models and algorithms for road network design: a review and some new developments,” Transport Reviews, vol. 18, no. 3, pp. 257–278, 1998. View at: Publisher Site  Google Scholar
 J. Sun, W. Xu, and B. Feng, “A global search strategy of quantumbehaved particle swarm optimization,” in Proceedings of the IEEE Conference on Cybernetics and Intelligent Systems, pp. 111–116, Singapore, December 2004. View at: Google Scholar
Copyright
Copyright © 2019 Jiehui Jiang et al. 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.