A Multi-commodity Network Flow Model for Sustainable Performance Evaluation in City Logistics: Application to the Distribution of Multi-tenant Buildings in Tokyo

: The distribution of goods in crowded city centers is a major challenge. In this paper, we propose a methodology for evaluating the performance of a parcel distribution network in city logistics. This methodology encompasses the main entities of a two-tier distribution system made up of carriers, huge shopping centers (multi-tenant buildings) and intermediate depots (urban consolidation centers), as well as the parcel ﬂows between them. This methodology aims to optimize the transport ﬂows (distance traveled) of a given distribution network while also quantifying the impact in terms of sustainable development by measuring gas emissions. Two di ﬀ erent states of the network with di ﬀ erent connectivity degrees are evaluated and compared: the current state of the network as well as its future state. The transport network modeling is based on a network ﬂow, which is expressed in linear programming and implemented with an optimization solver. The validation of this methodology is based on the parcel distribution of the Multi-tenant Buildings of the city of Tokyo. The ﬁndings are that the network with greater connectivity between the entities brings signiﬁcant traveled distance reduction as well as a reduction of emissions of CO2. Another ﬁnding is that the grouping of the parcels (i.e., pooling) brings a reduction of the distance traveled compared to the transport organization without grouping and contributes to a reduction in the number of trucks.


Introduction
City logistics aims at optimizing logistics activities in urban areas with the support of advanced information and communication technologies towards sustainable development [1,2]. Our study deals with the transport flow analysis in the context of the last-mile delivery. Even if the last mile represents a small part of the whole transportation chain, it is known to represent a significant part (about one-third) of the total transportation costs [3]. This problem generates traffic jams as well as an increase in air and noise pollution; it has a strong impact on the quality of life in cities and consequently undermines the objectives of sustainable development.
In this article, we focus on the distribution of goods from the outskirts of the city to the main buildings located in the city center of a megalopolis. In this context, two-tier distribution systems have emerged for many years; therefore, the main components of the distribution networks studied in this paper are as follows: • The multi-tenant buildings (MTB) are shopping complexes containing hundreds of shops, restaurants, offices, open spaces, and even sometimes a theater. These buildings with high population densities are the place of complex relationships between the owners, the tenants, and the management companies from which multiple problems can emerge. Let us mention, for instance, sustainable development problems such as energy consumption [4], and organizational problems such as regulating flows between different stakeholders. In this article, we focus precisely on the distribution flows of parcels to these buildings.

•
The urban consolidation center (UCC) is a way to consolidate the flow of deliveries in the last mile [5]. UCCs are central elements in vertical collaboration that make it possible to group parcels that have been addressed for certain customers in a single flow. These centers facilitate last-mile delivery by providing equipment that allows large vehicles to enter these spaces and consolidating the flows that can be formatted to suit the delivery mode more adapted to the last mile.

•
The carriers (CRRs) are independent companies that receive transport orders from the suppliers of the final customers located in the MTBs. They deliver parcels from their depot to a given final destination (i.e., MTB) and use intermediate logistics platforms such as UCC to optimize their own transport performance.
This paper aims at modeling and optimizing freight flows in order to compare the performances of different distribution networks. The measured performances cover traditional indicators such as distance and number of vehicles, but they also include indicators related to greenhouse gases and particulate emissions that are important in terms of sustainable development. The main objective of this paper is to highlight the interest of an open distribution network where the infrastructure can be shared by different transport operators. The following sections present a literature review, define the problem and working hypothesis, and present a linear programming model to optimize flows. Lastly, we report the experiments carried out for the freight distribution in Tokyo. A conclusion ends the paper.

Literature Review
Most of the work in transport planning involves the development of demand forecasting models coupled with network models to analyze distribution flows. For instance, [6] proposed an Origin-Destination synthesis model for a multi-commodity transport system. A commodity corresponds to a good (or a parcel) that has to be transported from a source node to a destination node. Our contribution also concerns the multi-commodity context, focusing on analyzing the distribution flows.
Research in the field of distribution in urban logistics is mainly oriented towards the design of distribution networks. The so-called Service Network design term is extensively used to refer to many design planning problems in transportation and distribution systems [7,8] proposed a low carbon network planning model for a city logistic distribution system; their model aims to decide opening distribution centers, and the objective function integrates operational costs with or without tax carbon for comparison purposes. [9] proposed to select among four kinds of typical distribution modes: self-distribution mode, mutual distribution mode, third party distribution mode, and common distribution mode while also reducing the carbon emissions. [10] presented a mathematical model to build a green distribution system minimizing the total cost made up of economic cost, environmental cost, and socioeconomic cost. The focus of our study is similar to the two-tier city logistics network proposed by [11]; however, ours is more close to the real-life case as the base network of freight carriers' depots, urban consolidation centers, and multi-tenant buildings already exists (for example, the Tokyo case study presented later in the paper). We propose a complete integration and analyze the performance of these existing (real-life) networks with and without complete integration.
At the core of the studied network are the urban consolidation centers. There are several studies [12,13] addressing advantages of collaboration in the form of reducing the negative impact of urban transport. For example, [14] used a multi-agent simulation in cases of urban goods delivery using urban consolidation centers and evaluated the system based on the multiple criteria including reducing delivery costs and CO2, NOx, and suspended particulate matter (SPM) emissions. [15] also reported similar advantages (reduced number of tours needed to meet the customer's demands and decreased the required distance) in a simulation to evaluate a UCC for distribution in Athens. This challenging topic has been reviewed in depth in [16]. The authors distinguished between vertical collaboration, which relies on breaking down transportation routes into multiple tiers, and horizontal collaboration, focusing on the question of centralized versus decentralized transportation planning. In this paper, we focus on vertical collaboration using UCCs in the last-mile delivery.
Apart from any studies which have reported the role and the potential impact of UCC on reducing nuisance [17][18][19][20][21][22][23], some real-life pilot projects implementing UCCs have also been reported [24,25]. For example, [26] reported a new type of urban consolidation center (Binnenstadservice.nl) in Nijmegen, which became successful by focusing on receivers rather than on carriers. On the other hand, the mixed success of the UCCs is reported due to several obstacles such as an increase of overall cost or space limitation in the city center. The analysis of the reasons for the failures or successes of UCCs is beyond the scope of this article and can be found in [27]. The recent review of UCC literature proposed by [28] highlights important findings, including the following three points: (i) many studies relate to the optimization of potential initiatives, but very few focus on existing initiatives; (ii) very few studies are oriented towards the overall impact of the UCCs on the distribution system; (iii) few studies measure or assess their environmental impact. Our contribution addresses some of these shortcomings. Indeed, we focus on the one hand on the whole network that includes the UCCs and that covers several entities of the distribution chain, and on the other hand, we carry out an analysis at the sustainable level of an existing and operating network in the city of Tokyo.
In order to assess the impact of distribution networks implementing UCCs in the centers of large cities, we chose to model the distribution of parcels from carriers to MTBs as a multi-commodity network flow (MCNF) problem.
An exhaustive review of the MCNF is beyond the scope of this article, and readers can refer to [29] for a comprehensive study of commodity flow models. The MCNF can be solved in polynomial time using linear programming if the flows are real values. However, the MCNF is an NP-hard problem if the flows are integers so it remains a difficult problem to solve [30]. Many formulations and many solution methods have been proposed to solve this problem using advanced mathematical programming approaches [31]. For instance, [32] presented a branch and price and cut algorithm to solve integers MCNF. The MCNF has been used in many domains such as telecommunication, production planning, logistics, and transportation. Regarding transportation, let us mention for instance a less than truckload (LTL) transportation problem which has been formulated as an MCNF in [33] and [34]. In our study, the overall objective is to ship the commodities through the network at the lowest possible cost while satisfying the arc capacities; furthermore, the carriers have to consolidate shipments in order to minimize the total time. To our knowledge, there is no work using a multi-commodity network flow model to assess the transportation flows in city logistics in the context of distribution of multi-tenant buildings. The model described in the next section is an extension of the integral MCNF, oriented towards the calculation of the flow of trucks into the network.

Materials and Methods
Our objective is to evaluate the performance of a parcel distribution network in terms of emissions and distance traveled. Our methodology can be described in the following three steps, corresponding to the next three sub-sections: • First we set out the working hypotheses in order to define the studied problem (Section 3.1). In this section, two key definitions are introduced: one characterizes the degree of connectivity of a distribution network (labeled as "integration"); the second concerns the grouping level of the flow of transported parcels (labeled as "pooling"). • Next, we propose a model of the distribution network in linear programming (Section 3.2). • Last, we address the present the experimental framework in order to validate our solution approach based on an example case in Tokyo (Section 3.3).

Problem Definition and Hypothesis
The first key definition relates to the degree of connectivity of the network used by different partners. We have considered two distribution networks in this study, depending on the degree of accessibility of the urban consolidation centers and the level of connectivity of transportation infrastructure: • The first distribution network is made up of UCCs, which are owned by transportation companies. Each facility is used by a unique owner; the service and space are not shared with other companies. In this case, UCCs are operated independently of each other, without any integration between different transport partners. Each UCC is dedicated to the deliveries to one or more consumers' locations in a specific multi-tenant building. In this configuration, only a single path exists from a carrier to an MTB with no alternative solution. Figure 1 represents a simplified case containing three MTBs, three UCCs, and three carriers with their own distribution chain. As there is no connectivity (integration) between UCCs, this network configuration is labeled as NI (i.e., No Integration).

•
The second distribution network is with an open transport infrastructure, where various transport operators share space, handling resources and other operational services. In this configuration, the UCCs are interconnected; therefore, the carriers can choose among alternative solutions from their depot to an MTB in order to optimize their own performance. These alternative solutions generate multiple potential paths for the parcels. For instance, many UCCs can be used in a chain to deliver parcels to the final destination. In this network, the transport operations performed by the transport companies are considered as fully integrated. Hence, this network is labeled as FI (Full Integration). Figure 2 represents this distribution network with integration of the urban consolidation centers.
The second key definition is about the pooling of parcels. From a general point of view, pooling consists in sharing of business processes and resources to provide scalable services for customers to gain optimized performance. Pooling requires companies to collaborate; for instance, by using joint resources with finite capacities. In the field of transport, sharing the urban consolidation centers among transport operators is a way to improve the overall performance of the distribution. Another way to further improve the transport performance could be the pooling of the vehicles.
Therefore, in the remainder of this study, we assume that the distribution network can be operated in different modes, whether the pooling of parcels from different carriers in the same truck is allowed or not. These modes are described as follows:

•
In the first mode, labeled as NoPooling (NP), the carriers are independent and the parcels coming from different source nodes (i.e., carriers' depot) cannot share the same truck.

•
In the second mode, labeled as WithPooling (WP), parcels from different source nodes to the same destination can be pooled in the same truck. Notice that the pooling potentially applies to all the arcs of the network, but its impact on the parcels flows is different depending on the degree of connectivity of the network.
In order to evaluate the transport flows of these various distribution networks and to carry out some comparisons between them, we made the following assumptions: • We focused on an overall (i.e., macroscopic) evaluation of the flows. Consequently, we focused on a mid-term evaluation without considering the routes optimization of the individual trucks.
• We focused on a one-day period of operation for the performance evaluation. Therefore, we did not address the specific time aspects (such as arrival time of the parcel) related to the entities (UCCs, MTBs) and to the flows between them.

•
We carried out performance evaluations based on a cumulative assessment of the costs (i.e., distance, time, etc.) occurring on a full path from source to destination (which consists of many arcs shown in Figures 1 and 2).

•
We focused on the forward distribution flows from the carriers (i.e., sources) to the final destinations (i.e., MTB). The number of parcels returned from the MTB is usually very low; therefore, the return flows are not as critical as the forward flows. Consequently, no return flows from the MTBs back to the carrier are taken into account in this study.

•
We assumed that the local administrative authorities can impose some limitations on the truck flows on certain arcs (i.e., maximum number of trucks). This limitation is in line with the aim of being able to reduce the negative impacts of freight transport in certain parts of the city. One example is the work of [35], which proposes an approach to reduce the environmental footprint of freight deliveries near sensitive urban facilities such as hospitals, schools, and retirement homes.

•
Lastly, we assumed that the goal of UCCs located near the city center is to ban direct deliveries from a source to a final destination in order to reduce congestion or other environmental factors in densely populated areas. Therefore, in our model, the parcels flows have to cross UCCs in their journey from a carrier to an MTB, instead of being directly delivered to MTB.

Modeling of the Distribution Networks
The studied multi-commodity network flow (MCNF) problem is modeled by a directed graph G = (V, A) where V is the set of nodes and A the set of arcs. The nodes are associated to each entity of the distribution network: UCCs, MTBs, and the carriers' depots. For each node i ∈ V, the set labeled suc(i) (respectively, pre(i)) characterizes the set of successors (respectively, predecessors) of node i in graph G.
Set A of arcs of the graph is partitioned into two subsets A 1 and A 2 , corresponding to the two levels (L1 and L2) of the distribution network (see Figures 1 and 2). These two separate levels are defined as follows: • The flows between carriers and UCCs (i.e., lower level) make up the first level of the two-tiers distribution problem. This level, labeled as L1, characterizes the part of the transportation network used by larger vehicles, which use traffic lanes outside the city center.

•
The flows between UCCs and MTBs (i.e., upper level), make up the second level of the two-tiers distribution problem. This level is labeled as L2. The flows (parcels) exchanged between two UCCs are also considered as a part of this level.
Set L defines the set of levels (i.e., L = {1,2}) and, set A l contains the level l arcs. A non-negative capacity cap i,j is assigned to each arc (i, j) ∈ A and characterizes the maximum number of parcels going through this arc. The arc capacity is expressed as a multiple of the capacities of trucks TruckCap l used in level l ∈ L. A travel cost (cost i,j ) is also defined for each arc, representing the distance between nodes i and j. In addition, we introduce set K of the delivery demands. Each delivery demand k (k ∈ K) has to be delivered from the depot of the carrier i (i.e., source node) to MTB j (i.e., destination node). This demand is modeled by assigning a positive quantity of parcels (supply k,i ) at the source node i and by assigning the corresponding negative value (−supply k,j ) at the destination node j.
Variables nbt i,j,k and nbtUT i,j,k correspond to, respectively, the number of full and partial (incomplete) trucks needed to satisfy the demand k on arc (i,j). Variable nbpUT i,j,k is the quantity of parcels corresponding to an incomplete truck (less than a full truckload) of the total demand k assigned to arc (i,j). Integer variable x i,j,k quantifies the number of parcels of demand k on arc (i,j).
We propose two formulations of the studied problems based on whether the pooling of parcels from different carriers in the same truck is allowed (Formulation WithPooling) or not (Formulation NoPooling). Both formulations use the same objective function (see Equation (1)) which minimizes the cumulated distances of all trucks over all the arcs of the graph. Equation (2) represents the capacity constraint on each arc (i,j). Equation (3) is the flow balance equation between input and output on arcs, which has to be verified for each node and for each demand. Equation (4) represents the truckload capacity constraint, which is used to calculate the number of full trucks and the number of parcels in the incomplete trucks. Equation (5) limits the number of parcels in an incomplete truck to the capacity of this truck. Equation (6) calculates the number of incomplete trucks, which is equal to zero or one. Equations (7)-(9) are corresponding equations in WithPooling formulation to Equations (4)-(6) of NoPooling formulation. Equation (6) and Equation (9) are linearized in the model implemented in the solver.
(NoPooling) min x i,j,k = nbt i,j,k · TruckCap l + nbpUT i,j,k ∀ l ∈ L, ∀ (i, j) ∈ A l , ∀k ∈ K (4) nbtUT i,j,k = nbpUT i,j,k / TruckCap l ∀l ∈ L, ∀ (i, j) ∈ A l , ∀k ∈ K The formulation WithPooling is a variant of the previous model, where the variables nbt i,j , nbtUT i,j and nbpUT i,j are no longer dependent on index k due to pooling of all demands.

Experimental Setup
In order to evaluate the distribution networks (Figures 1 and 2), we created some test scenarios based on the freight distribution in Tokyo. It should be recalled that the capital of Japan is one of the largest metropolitan areas in the world (about 39 million inhabitants) and one of the largest economies in the world. We firstly introduce the context of the distribution of Tokyo city as to why it has been chosen to define our test scenarios. Then, we describe the input parameters common to all test scenarios followed by the input factors, which vary in each test scenario. Lastly, we describe the performance indicators used to evaluate the distribution networks.

Distribution Context of Tokyo
The approach proposed in this paper for modeling and performance evaluation is generic and has no strong limitation regarding the number of MTBs and UCCs. However, our test scenarios are based on the logistics practices currently implemented in some parts of Tokyo, which is comparable to the situation depicted in Figure 1 (i.e., No Integration). It can be described as follows:

•
There are three main MTBs of Tokyo: Ginza six (G6) in the Chuo ward, Tokyo Sky Tree Town (TT) in the Sumida ward and Tokyo Midtown (MT) in the Minato ward as shown in Figure 3. • Four UCCs are in charge to distribute the parcels to these MTBs. The UCCs of Midtown and Ginza six are located in the basement of the MTB itself and serve to perform consolidation of deliveries addressed to the different floors. The MTB "Tokyo Sky Tree Town" is supplied by two external UCCs (Ariake and Shinsuna) as shown in Figure 2. All of these UCCs are dedicated to a specific MTB and there is no integration between them.

•
We have also considered many depots, owned by the carriers, mainly located in the area of the Tokyo bay where the goods are arriving by linehaul trucks. Each carrier can serve any of the MTBs.  Figure 3 gives the locations of MTBs, UCCs, and carriers on the city map. It also shows a possible path from a carrier to Shinsuna UCC and from this UCC to Tokyo Sky Tree Town as an example.

Common Input Parameters
There are two main input parameters, which are fixed for all scenarios: • The demands (supply k,i ), which are based on a combination of real and estimated delivery data. The real data comes from a major transportation company and logistics provider, operating the UCC at the Ginza six in Tokyo. The demands of the other MTBs are estimated using a proportional rule according to the surface area ratio of the considered MTB and Ginza six. Based on the real data, 18 carriers have been used in the distribution of parcels to the MTB terminals; one of these carriers has a dominant position and covers 43% of deliveries, whereas the top four carriers cover 88% of the demand.

•
The cost (cost i,j ) is the distance from node i to node j. The distances are average of the three distances obtained from Google Maps at 07:00 a.m., 11:00 a.m., and 2:00 p.m. in order to integrate the different traffic congestion situations. UCCs in Ginza six and Mid Town are located within these buildings; therefore, in these two cases, the distances (cost i,j ) between UCC and MTB equals to zero.

Tests Scenarios
The tests scenarios are created using different input factors, which are as follows: • Two possible values (45 or 180 parcels) for the capacity of the trucks (TruckCapl) have been considered. The lower capacity (i.e., 45 parcels) corresponds to a scenario where small trucks are used in the last-mile delivery starting from the UCC. These could be Light Electric Freight Vehicles as they are quiet, agile, emission-free, and take less space than conventional trucks and vans. The higher capacity (i.e., 180 parcels) corresponds to the average number of parcels for a three-ton truck, which is commonly used by the carriers. Consequently, we created two sets of possible scenarios based on types of vehicle in the network: either a homogenous fleet (HM), in which all vehicles have the same capacity over the entire network (180-180), or a heterogeneous fleet (HT) (45-180), with lower-capacity vehicles in the upper level (i.e., for (i, j) ∈ A 2 ) and larger vehicles in the lower level of the network (i.e., for (i, j) ∈ A 1 ). Table 1 shows the full list of test scenarios resulting from these input factors. The arcs capacities cap i,j of the network are unlimited (UL) or limited (LI). In the latter case, the values are assigned to the arcs in order to limit the emissions on a sub-area of the network.

Performance Indicators
In addition to the cumulated total distance of the trucks (TD) provided by the objective function (Equation (1)), some other performance indicators, including some sustainable development indicators (i.e., emissions) have also been calculated in a post-optimization stage. The list of indicators is defined below: • The cumulated total distance of the trucks TD (i.e., Equation (1)).

•
The number of full trucks (NFT), the number of incomplete trucks (NIT), and the total number of trucks (NT) are calculated according to Equations (10)- (12) in the NoPooling formulation of the problem.
NT= NFT + NIT Similar indicators are also calculated for the WithPooling formulation as well.

•
In addition, we define the percent gap of a given indicator "X", labeled as gapX (z1, z2). It is the relative variation of "X", generated by switching the two-levels input factor z, and it is defined as: gapX (z1, z2) = (Xz1-Xz2)/Xz2, keeping all other input factors unchanged. For example, the percent gap in distance generated by the capacities of the network is reported as gapTD (LI, UL) = (TD LI -TD UL )/ TD UL .

Results
The experiments (or scenarios) that have been conducted, aimed to assess the impact of the input factors on the performance indicators mentioned in the previous section. The models have been implemented with the IBM ILOG CPLEX solver. Table 2 shows some information about the size of the models: number of nodes and arcs and also number of variables and constraints. The computation times to get the optimal solutions are less than 30 minutes in every experiment.

Total Cumulative Distance
Firstly, we focus on a network where the capacities on the arcs are limited (i.e., scenarios numbered 1-8). The results presented in Figures 4 and 5 are intended to compare the total cumulated distance of the two distribution networks (NI and FI) depending on whether the capacities of trucks are homogeneous (HM) (180-180 parcels) or not (HT) (45-180 parcels) and depending also on whether the pooling of the parcels is allowed (WP) or not (NP). Figures 4 and 5 show that the distance is significantly reduced when integration between UCCs is considered.   Tables 3 and 4 give the impact (i.e., reduction gaps) of various input parameters. The "scenarios" column of these tables omits the comparison variable and replaces it with "xx"; for example, the label "NP-xx-LI-180" means that the impact of integration has been compared by reporting gapTD (NI, FI). Table 3. Gap (percent) of total cumulated distance-Limited arcs capacities.

Scenarios gapTD (NI,FI) Scenarios gapTD (NP,WP)
NP-xx-LI-HM 11.7 xx-NI-LI-HM 18.4 WP-xx-LI-HM 38.9 xx-FI-LI-HM 43.5 NP-xx-LI-HT 9.9 xx-NI-LI-HT 11.9 WP-xx-LI-HT 22.3 xx-FI-LI-HT 23.9 Table 4. Gap (percent) of total cumulated distance-Full integration.  Table 3a shows that the gap varies between 10% and 39% when integration is considered. On the other hand, the distance reduction gap in Table 3b varies between 12% and 44% when the effect of pooling (WP or NP) is evaluated; proving a high impact of pooling on the distance reduction as compared to the integration. Furthermore, Table 3 shows that the reduction in distance is greater when the homogeneous fleet is used as compared to the case when heterogeneous fleet is used; for instance, this reduction drops from 43.5% to 23.9% in the case FI. Indeed, the limitation of the truck capacities at the upper level of the network (45 parcels) generates an increase in the number of trucks and therefore in distance.
Secondly, we analyze the impact of the arcs capacities by comparing the limited capacities scenarios (number 1 to 8) with unlimited capacities scenarios (number 9 to 16). The impact of arcs capacities on the total distance is provided by the gap indicator gap (LI, UL) in Table 4a,b. This indicator, which is characteristic of the dataset used, can be interpreted as a deviation from a minimum distance boundary that can only be reached if there is no capacity constraint on the arcs (i.e., gap = 0 means that the smallest minimum distance is reached).
These distance gaps in Table 4a,b are only for the full integration scenarios as the capacities of the arcs play no role (i.e., gap = 0) when there is no integration (NI), whether they are pooling or not. However, the impact of the arcs capacities of the network is significant with the integrated distribution network (FI) as seen in Table 4.  Figure 6 considers homogeneous fleet while Figure 7 is for heterogeneous fleet. Both of these figures have the same trend. The integrated version of the network allows a reduction of CO2 compared to the nonintegrated version; this reduction varies from 5% to 23% in the heterogeneous fleet case and varies from 5% to 36% in the homogeneous fleet case.  The reductions in all emissions (CO 2 , NOx, SPM) are better for the FI-WP scenario, which exceeds the other three scenarios. Hence, the combination of pooling (WP) and network integration (FI) gives the best results in terms of emissions. This corresponds to a significant reduction in the cumulative total distance observed in the previous analysis (i.e., the TD indicator). The NI-WP scenario is intermediate in terms of emissions, demonstrating a very positive impact of pooling, even in non-integration (NI) scenarios. In contrast, FI-NP and NI-NP scenarios without pooling provide almost identical results that are the worst in terms of emission reductions. Figure 8 focuses on the distribution of CO 2 emissions over the two levels of the distribution network (i.e., L 1 and L 2 ) in the different scenarios. This figure shows that the introduction of zero-emission vehicles such as Light Electric Freight Vehicles in the last few kilometers of delivery would significantly reduce CO 2 emissions by keeping only low-level emissions (grey bars). A comparison of FI-WP scenarios with homogeneous and heterogeneous fleets shows that the increase in emissions generated by the use of smaller trucks at the high level (FI-WP-LI-HT) can be offset by the use of zero-emission vehicles, reducing CO 2 emissions to get the same orders of magnitude (101 and 109 kg/day) for both scenarios FI-WP-LI-HM and FI-WP-LI-HT. The average reduction of the CO2 emissions for all scenarios is 35%.

Number of Trucks
In this last part, we evaluate the impact of the pooling on the number of trucks used in the whole network. We focus on the variation of the number of full trucks (∆ NFT) and the variation of the number of incomplete trucks (∆ NIT), depending on whether pooling is allowed or not. Figure 9 shows that in the four scenarios shown, pooling makes it possible, on the one hand, to increase very significantly the number of full trucks (white part of bars) and simultaneously reduces the number of incomplete trucks (grey bars). The most significant variations come from the scenarios with full network integration (FI) compared to the scenarios without integration (NI). For example, in the full network integration scenarios, heterogeneous fleet scenario reduces the number of incomplete trucks by 63 units (NIT decreased from 84 to 21) and increases the number of full trucks by 26 units (NFT increased from 21 to 46). The net reduction of the total number of trucks is significant for all scenarios, as shown in Table 5 ("Total" column). The average reduction in the number of trucks for all scenarios is 29%. We can conclude that pooling has a very strong impact on reducing the number of trucks. The breakdown of this variation by network level shows that scenarios without integration (NI) are unable to reduce the number of trucks, which is in line with the fact that these scenarios are very constrained in terms of flow because they do not offer transport alternatives from the depot to the UCC. Table 5 also shows that scenarios with integration (FI) generate a higher average reduction on the lower level. This reduction is less important when heterogeneous fleet is used because the increased number of lower capacity trucks in the upper level balances the decreasing effect of pooling.

Conclusions
In this article, we propose a methodology for evaluating the performance of a parcels distribution network in urban logistics. The originality of this methodology is as follows: • The representation of the distribution problem in the form of a multi-commodity network flow which is, to our knowledge, a novelty in the field of urban logistics. The linear programming model makes it possible to optimize the total distance traveled and makes it possible to evaluate the impact in terms of emissions of the studied network. This model encompasses the UCCs, which are central elements in vertical collaboration; it also offers the possibility to limit the impact of the distribution flows in specific areas by using the arc capacities of the network.

•
We reported several comparisons. Firstly, the comparison of two types of distribution networks with different degrees of connectivity (i.e., integrated (FI) vs. nonintegrated (NI) networks); these two types of networks correspond to the current state (NI) and a proposed future state of the distribution system (FI). Secondly, we reported the comparison of two types of organization regarding the parcel distribution (With Pooling (WP) or No Pooling (NP)).

•
Lastly, this approach is generic and can be applied to optimize parcel distribution in city centers of large urban areas.
The experimental study of the distribution of parcels to the MTBs (large shopping centers) in the city of Tokyo highlights the benefits of using a fully integrated distribution network compared to a network without integration. A very important reduction of the traveled distance (up to 39%) and a significant reduction in emissions (up to 36% of CO2 emissions) can be obtained in fully integrated distribution networks. It also highlights the impact of pooling, which brings a reduction of up to 29% of the distance traveled compared to no pooling organization. Our study also makes it possible to quantify the reduction in emissions due to the introduction of zero-emission vehicles, such as Light Electric Freight Vehicles at the upper level (L2) of the network; this reduction is about 35% on average in all scenarios. Overall, the combination of an integrated network and the pooling of parcels improves the reduction of emissions, travel distances, and number of trucks, leading to the implementation of sustainable urban freight transport systems.
Future studies will focus on the following issues: (i) introducing time in the flow modeling, which is a key component for enhancing the model in order to carry out multiperiod optimization on a given planning horizon and also to be able to represent the behavior of the UCC over time (i.e., load over time); (ii) carrying out financial assessments of the different networks, taking into account the administrative costs of the UCCs; (iii) study the impact of new multi-tenant buildings resulting from the implementation of transportation infrastructure dedicated to the development of a new generation of ultra-high-speed trains in Tokyo; (iv) applying our distribution modeling approach to other cities in order to carry out comparative studies between the practices used in different countries in order to identify best practices.