A Multimodal Passenger-and-Package Sharing Network for Urban Logistics

-is paper envisions a multimodal passenger-and-package sharing (PPS) network for urban logistics integrating metro, taxi, and truck. A hub-and-spoke structure is designed including hubs located at metro stations and service stores connected to the hubs. Packages are transported by metro on backbone links between the hubs and are carried by taxis or trucks between service stores and hubs, depending on the unit costs of these twomodes and capacity constraint of the taxi. Amixed integer linear programming model for hub location problems—fusing the multiassignment p-hub median problem without capacity constraints and the capacitated multiassignment p-hub covering problem—is formulated to optimize the multimodal PPS network. -e model is implemented based on the real-world data in Shanghai (China) under a series of scenarios to evaluate the network performance from two perspectives: the number of hubs and the proportion of taxi drivers who are willing to carry packages. -e scenarios show that with increased number of hubs, the spatial distribution of hubs disperses from the city center to peripheral areas and more areas can be serviced by taxis. -ere is, however, a trade-off between the operation cost saved by taxis and the establishment cost of an extra hub.-e analysis also presents that if the proportion of taxis willing to carry packages associates with the incentive payments to taxi drivers, an optimal value of incentives exists, by balancing the operation costs of taxis and trucks.


Introduction
e demand for urban parcel delivery booms with the prevalence of e-commerce in the recent years. Take China as an example. Over 40 billion packages were delivered in urban areas in 2017, with an increase of 28% comparing with the previous year [1]. e booming demand for urban parcel deliveries potentially leads to heavy urban freight traffic that inevitably brings negative impacts on traffic congestion and air pollution. To mitigate the negative effects and to improve urban mobility, the European Commission recommends a passenger-andpackage sharing (PPS) system integrating passenger and freight transportation [2], and some typical services was overviewed and introduced by Sampaio et al. [3]. In general, three sharing options are proposed: share of space (e.g., multiuse lanes and night deliveries), share of public transport vehicles and networks, and share of urban delivery facilities such as delivery bays and goods lockers [4][5][6][7]. is study focuses on the PPS network design with passengers and packages sharing the same public transport vehicles and networks.
Current studies on the PPS system mainly focus on vehicle routing for package pickup and delivery, considering various scenarios like taxi sharing or self-owned vehicle routing between public transit stations and customers. For example, Li et al. [8] formulate a share-a-ride problem to optimize taxi routes and schedules to better match the supply of taxis and the demand from passengers and parcel deliveries. e model is evaluated in Tokyo city by Nguyen et al. [9] and extended by Li et al. [10] by integrating stochastic travel times and delivery locations. Masson et al. [11] propose a public transit-based PPS system that firstly uses buses to deliver the packages from a consolidation/distribution center to a set of bus stops and then dispatches a fleet of freighters to deliver packages to customers. Taking into account time windows, schedules, and stochastic demands, Ghilas et al. [12,13] use the public transit as a part of freight journey and optimize the pickup-and-delivery routes of private vehicles between given transfers and demand nodes. Fatnassi et al. [14] assess the potential of integrating personal rapid transit and freight rapid transit vehicles into one network with the objective of minimizing the movements of empty vehicles and the waiting times of passengers and packages. Yildiz and Savelsbergh [15] explore the service and capacity planning problem in store-to-door meal delivery, where the delivery service is provided by individual's willing to participate in such activity. e contribution of this study-different from current studies on vehicle routing in a PPS system-is to design the multimodal PPS network with the hub-and-spoke structure. e idea of the hub-and-spoke network design for multimodal logistics networks is widely supported by past scholars [16]. By integrating multiple modes into one logistics network, the new system is able to achieve economies of scale more effectively, thus improving the transport efficiency [17,18]. Current literatures mainly focus on the intercity multimodal logistics network. Qu et al. [19] design a logistics network integrating ocean and coastal routes, inland waterways, railways, roadways, and airways with nonlinear transfer costs between different modes. Taking into account demand uncertainties and network disruptions, Fotuhi and Huynh [20] identify the rail links to retrofit, existing intermodal hubs to expand, and locations of new hubs in the rail-road network. Similarly, Kim and Ryerson [21] identify a set of unaffected hubs during network disruption in the US ocean-inland freight network.
A conceptual design of the multimodal PPS network for urban parcel delivery is envisioned and drawn in Figure 1, integrating metro carriages, taxis, and trucks. e fusion of metro, taxi, and truck absorbs the advantages of the three modes: metro trains follow predetermined routes and schedules with high reliability, punctuality, and large carrying capacities; taxis and trucks are flexible and are able to provide door-to-door services. e carrying capacity of taxis is constrained by the willingness of drivers to help carrying packages. e capacity of trucks is unlimited since the freighters are fully controlled by the couriers. However, the delivery cost of trucks is higher than the cost of the shared taxis.
To explore the optimal design of a multimodal PPS network, we formulate the network design task into hub location problems (HLPs). e rest of the paper is organized as follows. Section 2 reviews current HLPs and states the novelty of the proposed multimodal PPS network. In Section 3, we propose a variant of the HLP model and develop a mixed integer linear programming model with a modified genetic algorithm to optimize the hub-and-spoke network. Section 4 presents an example of Shanghai for the case study and provides quantitative assessment of the network from the perspectives of the number of hubs and the proportion of taxi drivers who are willing to carry packages. Finally, we conclude the study and offer some discussions.

Literature Review: Hub Location Problems
Hub location problems are firstly proposed by O'Kelly in 1980s [22][23][24]. Typical HLPs can be classified based on their objectives: the p-hub median problems locate p hubs to minimize the total routing cost [25][26][27][28][29][30][31]; the p-hub center problems locate p hubs so that the maximum cost between any OD pair can be minimized [25,32]; and the p-hub covering problems aim to maximize the total flow between all OD pairs covered by p hub facilities [33]. In these studies, the number of hubs to be established, i.e., p, is prespecified. Instead of predetermining the number of hubs, some works assume a setup and operation cost for each hub with the objective of minimizing the sum of hub cost and routing cost [29,34] or maximizing the system profit [35]. Comprehensive reviews can be found in Alumur and Kara [36], Campbell and O'Kelly [37], and Farahani et al. [38]. e network topology of HLPs depends on the structure of backbone links and spoke links. Typically, the hubs are completely connected. However, in some special cases, incomplete interhub network structures, e.g., ring [39], star [40], tree [41,42], and line [43], are adopted. In terms of the design of spoke links for flow assignment, a demand node can be connected to one or multiple hubs, i.e., single or multiple assignment structures. Aykin [44] and Martins de Sá et al. [43] extend the topology of spoke links by considering that the flows between some OD pairs can be shipped directly without being routed through hubs. Besides direct links, Lin and Chen [45] integrate stopovers and feeders to the network. Klincewicz [46] introduce a spoke network in which a demand node is allowed to be linked to a hub through other demand nodes. In addition to the backbone-and-spoke bilevel networks, some researches consider the trilevel network connecting hierarchical hubs and demand nodes [40,47]. e capacity of links is another essential feature of the HLPs. Capacity constraints on flows are usually applied on hubs [25,[48][49][50][51][52], links [53,54], or both [55], as exogenous constraints. Some studies also treat hub capacities as decision variables to select the optimal hub capacity from a list set of predefined capacity levels, each of which is associated with a fixed setup cost [56][57][58]. Consider that the capacity limits make the model more realistic, though it increases the computation time in most cases. e hub location-routing problem, introduced by Nagy and Salhi [59], combines the HLPs with the vehicle routing problems (VRPs), where the vehicle routes are optimized with respect to hub locations and flow assignments simultaneously. e combination of HLPs and VRPs, both of which are NP-hard problems, leads to a great growth in the solution cost. us, the hub locationrouting problems are normally formulated with a hierarchical structure, in which the upper and lower levels deal with the HLPs and VRPs, respectively [59][60][61][62]. In this study, the vehicle route planning is not considered for simplification.
is study applies a complete interhub structure. A demand node is allowed to be connected to multiple hubs to improve the flexibility and capacity of the network, as metro stations are completely connected by metro lines and transfer stations. e spoke links are served by taxis and trucks. e features of the two modes-e.g., the costs to deliver a package and the capacity constraints-vary largely. e study also takes into account the exogenous constraints on spoke links served by taxi due to the drivers' willingness to help carrying packages. e spoke links served by truck and the backbone links served by metro are uncapacitated.
e fusion of these modes enriches the formulation of HLPs. Figure 1 illustrates all potential routes of an origin-destination (OD) pair. e service stores spread out in the city for first-mile package collection and last-mile delivery. e hub facilities are located at metro stations for package storage and transfer, fully connected by the metro network. Between service stores and hubs, packages are transported by taxi with limited carrying capacities and by truck without capacity limits, according to the costs of these two modes. In the hub-and-spoke network, the backbone links and the spoke links are served by metro and taxi/truck, respectively. e envisioned multimodal PPS network also includes the following characteristics. Standardized delivery boxes-fitting into the trunks of taxi and truck-and containers, designed for metro carriages, are used for parcel transfer. e boxes and containers are (un) packed and reassembled at stores and hubs. e requests for package delivery to taxis can be automatically matched by a central operation platform, guiding the taxi drivers if they need to transport packages besides passengers. e platform is also responsible to rent trucks from cooperative on-demand logistic companies to deliver the packages that cannot be covered by taxis and to arrange the transportation of package in the metro network. e following assumptions are made for modeling:

Model Settings and Assumptions.
(i) e service stores are the O/D nodes of the model.
(ii) Packages have to be routed via at least one hub.
(iii) Packages of an OD pair can be routed through multiple hubs.
(iv) On spoke links, the package flows served by taxi are constrained by taxi flows, the proportion of drivers who are willing to carry packages, and the trunk capacity of each taxi. e package flows served by truck are unlimited.
(v) e capacities on hub facilities and backbone links served by metros are unlimited.
Finally, we focus on the low-cost urban courier service without the constraints of the same-day delivery. erefore, the system is evaluated in the dimension of cost rather than time. Table 1 summarizes the inputs, parameters, and variables used in this study. e package flows in this model are treated as multicommodity flows [63]. Each commodity i represents the packages originating from node i. For commodity i, denote z ik and u ik as the flows routed from the origin service store i to hub k by taxi and truck, respectively, y ikl as the flows routed between hubs k and l, and x ilj and v ilj as the flows routed from hub l to destination service store j by taxi and truck, respectively. e multimodal PPS network design problem is formulated as the mixed integer linear programming (MILP) problem, expressed as follows: Objective function (1) is set to minimize the total cost (TC), consisting of 6 terms, which indicate (i) e cost to transport packages (1) From origin stores to hubs by taxi;

Model Formulation.
(2) From hubs to destination stores by taxi; (3) From origin stores to hubs by truck; (4) From hubs to destination stores by truck; (5) Between hubs by metro; (vii) e cost to establish and to operate the hubs, respectively.
Constraint (2) guarantees that there are exactly p hub facilities to be established. Constraint (3) and constraint (4) represent that for commodity i, and the total package flow originated from node i and destined to node j is in accordance with the demand. Constraints (5) and (6) ensure that package flows transported by taxi between demand nodes and hubs are constrained by taxi flows, the WTCP%, and the capacity of the taxi trunk. Constraint (7) is the flow conservation equations for each commodity i at each hub k, of which the left and right sides are outgoing and incoming flows, respectively. Constraint (8) indicates that the package flows related to hub k are greater than zero only if hub k is established, where Q is a large number. Constraints (9) and (10) specify the binary and nonnegative types of decision variables. e proposed model generalizes and is an extension of two classic types of HLPs. If no taxi drivers are willing to carry packages (i.e., r ij � 0), or if the payments for incentive c ij are set to be relatively high, the package flows by taxi become zero. e spoke links are serviced by truck only. In this case, the model becomes a typical multiassignment p-hub median problem without capacity constraints. If the costs to transport packages by truck t ij are set to be relatively high, the model tends to minimize the amount of trucks and to maximize the package flows transported by taxi with capacitated spoke links. In this case, the model becomes similar to a capacitated multiassignment p-hub covering problem.

Modified Genetic Algorithm Method.
e developed MILP model to solve the HLPs is NP-hard due to the interrelationship between locating hub facilities and designing network topology. Past studies focusing on exact algorithms, e.g., branch-and-bound [64], branch-and-cut [65], Lagrangian relaxation [66], and Benders decomposition [67], mainly handle the HLPs with limited scale. To solve the HLPs with a larger scale, heuristic algorithms-e.g., the general variable neighborhood search approach [68], the memetic algorithm mapping two local search heuristics developed by Marić et al. [69], and the genetic algorithm (GA) adopted by Lin et al. [55] for capacitated p-hub median problems-are able to provide the near-optimal solutions with less computation time.
In this study, we propose a modified GA method to solve the large-scale HLP problems. e HLP model consists of binary variables, indicating hub locations, and non-negative continuous variables indicating package flows. e modified GA considers the set of binary variables as an individual. For each individual, the fitness function is the objective function of the model mentioned above, where binary variables are set fixed. e detailed algorithm is presented in Appendix with the parameters empirically set as in Table 2.
To investigate the validity and efficiency of the modified GA, we conduct a series of computational experiments comparing the computation times and optimal solutions generated by the GA with those solved by exact algorithms in CPLEX (with default parameters) [70], listed in Table 3. We use (N, H, p) to denote the problem scale, where N is the number of OD nodes, H is the number of candidate hubs, and p is the number of hubs to establish. Table 3 shows that when N is greater than 110 and/or H is greater than 80, CPLEX fails to provide a feasible solution in some cases within max time limit (i.e., 24 hours).
e results show that the computation time of the exact algorithm increases rapidly with the number of stores and hubs. In terms of the quality of solutions, if optimal solutions can be efficiently provided by the exact algorithms, the GA achieves the same or similar results. If suboptimal solutions Proportion of taxi drivers who are willing to carry packages from node i to j, namely, "WTCP%" θ Average number of packages that could be contained in a taxi trunk g k Cost for establishing and operating hub k c ij Incentive payment per package for taxi drivers from node i to j m kl Transport cost per package from hub k to l by metro t ij Transport cost per package from node i to j by truck p Number of hubs to be established Variables h k Binary variable indicating whether a hub is established at node k. 1 if so, 0 otherwise z ik Package flow routed from node i to hub k by taxi u ik Package flow routed from node i to hub k by tuck y ikl Package flow emanating from node i routed from hub k to hub l by metro x ilj Package flow emanating from node i routed from hub l to node j by taxi v ilj Package flow emanating from node i routed from hub l to node j by truck with gaps can be provided by the exact algorithms, the GA can even achieve better. A feasible solution can still be provided by the GA when the exact algorithm fails to provide a feasible solution in reasonable time.

Data and Study Area.
e study area is encircled by the outer-ring expressway of Shanghai, China, covering 680 square km with great population densities and high economic activities. We tessellate the study area into 171 identical square zones with the length of 2 km. In each zone, one service store is designated to service the area. ere are 106 zones containing at least one metro station, which are considered as the candidate locations to establish hub facilities. In addition, the zones with main land use of agriculture, industry, or nature parks are excluded from the analysis due to scarce demand for delivery service in these areas. us, 127 zones are selected for the optimization of the multimodal PPS network with 94 candidate hub locations. e spatial distribution of zones and candidate hubs is presented in Figure 2.
630 million packages were delivered in Shanghai in 2015, with the spatial distribution presented in Figure 3. e distribution is heterogeneous, and, in general, the demand concentrates in the central areas. In the top ten zones with the highest demands, 31.1% of the total packages are sent out from these areas and 17.2% of the packages are delivered into the areas.
Based on the automatic vehicle location data of about 13,000 taxis-accounting for one fourths of the total taxis operated in Shanghai-in April 2015, we extract the origin and destination of the taxi trips and scale up the number of the trips to synthesize the entire taxi service in Shanghai. e results show that most taxi trips are originated from or destined to metro stations, where the zones with metro stations generate 90.1% of taxi trips. Figure 4(a) presents the spatial distribution of the origins of taxi trips, while Figure 4(b) shows the distribution of the destinations of taxi trips. e two figures suggest that it is reasonable to utilize metro stations as transfer locations between metros and taxis.
We empirically calculate the cost to establish and operate a hub as 8k CNY per day, based on the average cost to lease a shop in metro stations (according to the data from 58.com Inc., a China's Craigslist equivalent), the area of hubs for package storage, and the equipment and labor costs in Shanghai according to the Shanghai Statistic Yearbook in 2016. e cost to carry the packages by truck is estimated according to the pricing plan of renting a mediumsized truck from Lalamove, an on-demand logistics company providing intracity delivery and courier services in Shanghai (for each truck trip: 65 CNY for the first 5 km and 4 CNY/km for more than 5 km (source: Lalamove)). e unit cost for transporting a package by metro is set as 5% of the regular fare (for each customer: 3 CNY for the first 6 km and 1 CNY per 10 km for more than 6 km, at most 15 CNY (source: Shanghai Shentong Metro)). Each taxi and truck is able to carry 10 and 50 packages, respectively. e average distances between the zones are calculated as the Euclidean distances between the centroids of each zone.

Scenario Analysis.
Among a number of factors affecting the system performance, we designed three scenarios to evaluate the proposed multimodal PPS network from the following perspectives: (i) Number of hubs; (ii) Proportion of taxi drivers who are willing to carry packages, named as "WTCP%": (1) WTCP% independent from the incentives (2) WTCP% dependent on the incentives.
Each scenario contains several subscenarios. Table 4 summarizes the attributes of each scenario and subscenario in details. Among the subscenarios, A6, B3, and C3 are the same, while B1 and C1 are the same.
Scenario A consists of 7 subscenarios. In the scenario, the impacts of different numbers of hubs p (i.e., 5, 6, 8, 10, 12, 15, and 20) are tested. e WTCP% is fixed: 50% taxi drivers are willing to carry packages. To compute the delivery cost by taxi on each spoke link, we set that a taxi driver is able to receive an extra 5% of the regular trip fare (for each taxi trip: 14 CNY for the first 3 km, 2.5 CNY/km for 3-15 km, and 3.8 CNY for more than 10 km (source: Shanghai Municipal Transportation Commission)) for every additional package to deliver. With the capacity of 10 packages, a taxi driver may get at most 50% more than the regular fare.
Scenario B consists of 5 subscenarios. Different WTCP% (i.e., 0%, 25%, 50%, 75%, and 100%) are tested with the incentive fixed at 5% of the trip fare. e number of hubs is fixed to 15. e subscenario B1 represents a metro-truck only network without taxi. e formulation becomes a typical multiassignment p-hub median problem without capacity constraints.
Scenario C consists of 5 subscenarios where C1 is the same as B1. Different from scenario B, the WTCP% depends on the incentive. We assume that, with an extra 0%, 2.5%, 5%, 7.5%, and 10% of regular taxi fare for each package, 0%, 25%, 50%, 75%, and 100% of taxi drivers are willing to carry packages. e number of hubs is also fixed to 15. Figure 5(a) illustrates the trends of total cost and transport cost with increased number of hubs, in which the difference of two series of values represents the cost to establish and to operate the hubs. More hubs lead to more packages transported by taxi, which    Journal of Advanced Transportation        Journal of Advanced Transportation reduces the transport costs on spoke links, as shown in Figure 5(b). Nonetheless, the marginal savings of transport cost declines with more hubs established. With 5 additional hubs established, the increase of the hub cost is 40k CNY. From 5 hubs to10 hubs, 18.2% of the packages are shifted from the truck to taxi while the transport cost reduces by 89.5k CNY. From 15 hubs to 20 hubs, the packages' modal shift from the truck to taxi is 7.9% with only 34.2k CNY of transportation cost saved. According to the figures, the total cost reaches minimum with 15 hubs. Beyond 15 hubs, the cost to set up an extra hub becomes greater than the reduction of the transportation cost. Figure 6 presents the proportions of packages transported by taxi and the spatial distribution of hubs by subscenario. With 5 hubs, the system mostly assigns the hubs into the central area with greater demand. As a result, the periphery zones close to the outer-ring expressway are underserved by taxis. Among all 127 zones, there is no zone in which more than 50% of the packages can be served by taxi. e rest of the packages are transported by truck with a higher cost. With 6,8,10,12,15, and 20 hubs, the system covers more periphery zones and the number of zones where over 50% of demand can be served by taxis increases to 2, 11, 31, 68, 110, and 127, respectively. e optimal locations of hubs in each multimodal PPS network are demonstrated in figures. A link connecting to a service store (blue dot) and a hub (red dot) represents the package transported between a service store and a hub by

Journal of Advanced Transportation
taxi. e color of a link is the same as its origin node. e width of the link represents the package flow through the link. In the figures, there are several hub facilities clustered together in the central area. One plausible reason is that the package demand and taxi flow are both higher in the central area. e establishment of hubs in these zones is able to service more demand with limited resources. If more hub facilities are allowed, some hubs spread out near the outerring expressway to expand the coverage of demand. Take the scenario with 15 hubs as an example. e detailed routing procedures for packages from an origin service store to all destination stores are illustrated in Figure 7. e routing process consists of three phases. In the first phase, most of the packages are transported by the taxi to the close hubs. Due to the constrained carrying capacity between the service store and hubs, the rest are transported to hub by trucks. In phase 2, the packages are transported by metro to the hubs connecting to their destination stores. Finally, the packages are routed from the hubs to the respective destination stores in phase 3. On spoke links originated from the hubs in the central area, the packages are mostly transported by taxis, while on the spoke links originated from the hubs in other areas, most packages are transported by trucks due to the insufficiency of taxis.

Scenario B: WTCP% Independent from Incentives.
Under this scenario, the payments for incentive are independent from the WTCP%. If there are no taxi drivers willing to carry packages, the PPS network becomes a metrotruck system, which generates relatively high transport costs on spoke links. As shown in Figures 8(a) and 9(b)  % increases, thanks to lower delivery costs comparing with trucks, the total cost decreases and more packages are carried by taxis on spoke links.
With more taxis willing to carry packages, the optimal hub locations become more dispersed in the city. Figure 9 illustrates the specific hub locations of the five subscenarios under scenario B, with the WTCP% and the average hub distance listed in the caption. If no taxi drivers are willing to carry packages, the hub locations are optimized without considering the distribution of taxi flows. If 25% of taxis are willing to carry packages, the optimal network tends to cluster hubs in the central area of the city where the package demands and taxi supplies are both high. e layout allows more service stores in the central area linking to hubs. With the increase of WTCP%, the central area can be serviced by fewer hubs and some hubs are relocated to peripheral areas to cover more package demands.
With more taxis involved, more trucks are expected to be saved in the logistics system. We divide the study area into three parts-central, middle, and outer parts-according to the locations of Shanghai inner-, middle-, and outer-ring expressways (see Figure 10(a)). With WTCP% increasing from 0% to 100%, 89.32%, 76.43%, and 59.70% of truck flow is saved in the central, middle, and outer parts, respectively (see Figure 10(b)). e central part benefits the most from the PPS network, followed by the middle and outer part. is result implies that the PPS network could mitigate the negative effects of urban logistics on traffic congestion and air pollution in the urban areas.

Scenario C: WTCP% Independent from Incentives.
Under this scenario, the WTCP% depends on incentive payment-the unit cost of taxi delivery. Figures 11(a) and 11(b) present the trends of total cost and proportion of package flows delivered by taxis on spoke links, respectively. Although the increased WTCP% allows more packages to shift from the truck to taxi, the increase of incentive raises the total costs of the system. With greater incentives, the PPS network prefers trucks to transport packages on most spoke links, leading to the reduction of taxi delivery. e results show that under the linear assumption of incentive-WTCP% relationship, the optimal payment for incentive is around 2.5% of regular taxi fare per package, with 25% of WTCP%. In this case, despite of the fact that the proportion taxi usage is not the highest, it generates the minimum total cost due to lower payments for incentive.

Conclusion and Discussion
is study envisions a multimodal PPS network integrating metro, taxi, and truck with a hub-and-spoke network. On spoke links, the PPS network balances the usage of taxis and trucks for package transport according to the unit costs of these two modes and capacity constraints of taxis. e design of the multimodal PPS network is formulated into a HLP with the mixed integer linear programming framework. e developed model is an extension to the two classic types of HLPs: the multiassignment p-hub median problem without capacity constraints if the spoke links are serviced by truck only, and a capacitated multiassignment p-hub covering problem if the model tends to minimize the amount of trucks and to maximize the package flows transported by taxi with capacitated spoke links. A modified GA is developed taking into account the constraints of carrying capacity by taxi.
To demonstrate the potentiality of the multimodal PPS network for urban logistics, the network is adopted based on the real-world data of Shanghai. A series of scenarios are designed to assess the performance of the system from two perspectives, the number of hubs and the proportion of taxi drivers who are willing to carry packages. e scenarios show that with increased number of hubs, the spatial distribution of hubs disperses from the city center to peripheral areas and more areas can be serviced by taxis. More hubs lead to more packages transported by taxi, which reduces the costs on spoke links, but the marginal savings of transport cost declines as well. ere is, nevertheless, a trade-off between the operation cost saving by taxis and the establishment cost of an extra hub. e analysis also presents that with given number of hubs, if the WTCP% is independent from the incentive payments, with more taxis willing to carry packages, the total cost of the system keeps decreasing.   However, if the proportion of taxis willing to carry packages increases with the incentive payments to taxi drivers, an optimal value of incentives exists, by balancing the operation costs of taxis and trucks.
For future works, some hypothetical assumptions made in this study can be further released. For instance, the current network requires the packages to be routed via at least one hub, while a network that allows direct connections between service stores could extend the network coverage. In addition, the influence of various partitioning methods of the study area on the resulting network can be further evaluated. It is also recommended to consider the willingness of taxi drivers to participate in PPS activities as a decision variable in network optimization. It can be an effective means to adjust the carrying capacities since more taxis may be willing to participate in the PPS activities with greater incentives. Finally, the PPS network may lead to the rebalancing of taxis, which may result in different optimized network patterns.

Appendix
Steps and Pseudocodes of the Modified GA e modified GA iterates from Step 3 to 6 and stops if the convergence criterion is satisfied.
Step 1. Coding: a hub location arrangement is represented by an individual consisting of a binary string of length N p . I(k) = 1 if node k is selected to locate a hub facility. Otherwise, I(k) = 0.
Step 2. Initialization: a population of individuals is randomly generated with a size of K. For each individual, p genes are selected from in total n genes without replacement and assigned as 1. e other genes equal 0. Each individual has exactly p number of "1," representing p nodes to locate hub facilities.
Step 3. Fitness evaluation: the fitness function is set as the objective function of the model developed in Section 3.2. Giving selected hub locations, the subproblem becomes a linear programming problem involving continuous variables x, y, z, u, and v that can be solved efficiently to evaluate the fitness of each individual.
Step 4. Selection: the fine-grained tournament selection strategy (FGTS) is adopted to select K individuals from the population [71]. Denote F t as the desired average tournament size. e strategy performs weighted sampling of |F t | + 1 with replacement from the population k 1 times and performs weighted sampling of |F t | with replacement k 2 times. e weight of an individual equals its fitness value. For each sample, the one with the lowest fitness value is selected. e values of k 1 and k 2 are determined by the following equations: (A.1) Step 5. Crossover: two offsprings from a pair of parents are produced. e basic crossover operator randomly exchanges the codes of each gene. e number of "1" in the resulting offspring may not equal p. To address this problem, we simultaneously exchange the codes for two genes to ensure that each offspring has p "1" values. Algorithm 1 describes the crossover.
Step 6. Mutation: the gene code of an individual is randomly shifted from 0 to 1 or from 1 to 0. To ensure that there are p number of "1" in an individual after the mutation operation, we change the codes of two genes I 1 ⟵ First parent; I 2 ⟵ Second parent N p ⟵ Number of potential hub locations i ⟵ 0; j ⟵ N p + 1; while i < j for k ⟵ (i + 1) to (j − 1) if I 1 (k) � 1 and I 2 (k) � 0 i ⟵ k; break end if i ⟵ N p end for for k ⟵ (j − 1) to (i + 1) if I 1 (k) � 0 and I 2 (k) � 1 j ⟵ k; break end if j ⟵ 0 end for if i < j Exchange (I 1 (i), I 2 (i)) and (I 1 (j), I 2 (j)) with probability P c end if end while ALGORITHM 1: Crossover operation.
Journal of Advanced Transportation 13 that have different codes simultaneously. Algorithm 2 describes the mutation. Convergence criterion: the algorithm stops when (1) the average relative change in the fitness function value over stall generation limit G s is less than function tolerance τ or (2) the computation time exceeds time limit T max .

Data Availability
e data used to support the results of this study are available from the corresponding author upon request.

Conflicts of Interest
e authors declare that they have no conflicts of interest.