## Analysis of Transshipment in Three-Sided Meal Delivery Services via Microhubs

Linxuan Shi 1 and Zhengtian Xu * 1

1 Department of Civil and Environmental Engineering, The George Washington University, Washington DC, United States

April 8, 2025

## Abstract

This paper introduces and analyzes a novel transshipment strategy for meal delivery. In this approach, the service area is partitioned into smaller sub-areas, with deliverers assigned to operate exclusively within these sub-areas. Meanwhile, a centrally located microhub functions as a logistic depot to facilitate the batching and transfer of meal packages toward different sub-areas. We model the meal delivery system with transshipment using networked G/G/m queues and analytically approximate two critical system performance metrics-customer waiting time and vehicle miles traveled-to evaluate the effectiveness of the proposed strategy. The performance achieved by transshipment is benchmarked against that of the classic pickup-and-delivery strategy without transshipment, both predicted using continuous approximations. For the latter, we enhance the existing modeling by incorporating the delivery distance profiles of individual orders to better match the meal delivery context. Our comparisons indicate that meal delivery via transshipment outperforms the non-transshipping counterpart across both metrics under either high-demand or low-supply conditions, with particular advantages in servicing larger areas or handling long-distance orders. This conclusion is corroborated by a numerical experiment using empirical meal delivery data from Meituan, which suggests that an optimally configured transshipment strategy can significantly improve service performance for both customers and deliverers during peak lunch hours and in the busiest districts. While transshipment continues to reduce vehicle miles traveled by deliverers during non-peak hours, it results in longer customer waiting times compared to the benchmark without transshipment as demand decreases.

Keywords - meal delivery; transshipment; pickup-and-delivery problems; continuous approximation; discrete-event simulation

* Corresponding author. E-mail address: zhengtian@gwu.edu (Z. Xu).

## 1 Introduction

Meal delivery services connect customers and restaurants through mobile applications, allowing users to order meals online with just a few taps on their smartphones. The platform assigns meal orders to a team of affiliated deliverers, who complete the fulfillment process by picking up the meal package from the restaurant and delivering it to the customers. The lasting COVID19 pandemic prompted an irreversible shift among consumers from indoor dining to online meal delivery (Kaczmarski, 2024), under which takeout channels have become increasingly vital to restaurant operations in the post-pandemic era (Shi and Xu, 2024). The significant surge in customer demand has heightened the need for meal delivery platforms to innovate supply management to ensure timely delivery. To prioritize customer satisfaction, platforms have been continuously seeking faster last-mile delivery methods, such as strategic order batching and assignment (Liu et al., 2021; Simoni and Winkenbach, 2023; Ye et al., 2024), advanced vehicle routing strategies (Steever et al., 2019; Kohar and Jakhar, 2021), and new delivery modes utilizing autonomous vehicles and drones (Liu, 2019; Shi et al., 2022). As intermediaries in a three-sided service market involving customers, deliverers, and restaurants, meal delivery platforms must carefully consider and manage the differing interests of these stakeholders when introducing new delivery operations (Agatz et al., 2024).

This study conceptualizes a new operational model for meal delivery by enabling transshipment through a microhub, aiming to enhance overall efficiency and provide long-term benefits to all stakeholders. The concept of microhubs mirrors logistic depots in line-haul logistics, allowing for sorting, storing, and transferring packages, but with a focus on urban neighborhoods. As illustrated in Figure 1c, the availability of these depots re-structures freight distribution into a transshipment-based operation, under which packages are transferred from line-haul trucks to local distribution networks for final delivery. Compared to the delivery model without transshipment, shown in Figure 1a, this approach brings services closer to end customers and enables more flexible operations within targeted service areas. In some use cases, it alleviates congestion by consolidating freight vehicle trips and lowers emissions by facilitating transfers to low- or zero-emission fleets for last-mile deliveries (Gunes et al., 2024).

Unlike the traditional line-haul distribution system shown in Figure 1a, which typically involves one cluster of origins and another of destinations, the meal delivery system features orders with pickup and drop-off points dispersed across the entire region, as illustrated in Figure 1b. In our proposed design, illustrated in Figure 1d, the service region is divided into pie-shaped sub-areas, with one microhub located at the junction, to facilitate localized batching for order collection and distribution. The microhub serves as an intermediate node for the temporal storage and transfer of meal packages, while also coordinating pickup and drop-off tasks both within and across sub-areas. Individual meal packages are picked up from the restaurants where they originated and transported to the microhub. Subsequently, one deliverer collects the package from the microhub and delivers it to the customer. Delivery vehicles are dispatched within each area to perform pickup, drop-off, and transfer tasks, with each vehicle dedicated to serving a specific sub-area. In practical implementation, deliverers possess better knowledge in routing for regions they are familiar with and prefer to operate exclusively within these areas for greater efficiency. In contrast to transshipment in line-haul logistics, which exploits transportation economies of scale by moving goods in large batches over shared long-haul segments, our design of transshipment through a microhub, tailored for intra-city many-to-many distribution contexts, aims to gain efficiency through localized order batching. To evaluate the comparative efficiency gains provided by introducing transshipment, this study employs continuous approximations (CA) to

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.01047v2-figure-f0e46089-31ac-4fd3-9393-4ff5547170a3.jpg)

(a) Line-haul logistics without transshipment

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.01047v2-figure-7eebce2e-f77b-44c9-8116-9bc06af6f0e4.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.01047v2-figure-2c99b05e-fc41-4fbf-b0b3-b2c08a629cd8.jpg)

(b) Meal delivery without transshipment

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.01047v2-figure-a92de1fb-2cd4-40dc-9645-1a2490e45bf2.jpg)

- (c) Line-haul logistics with transshipment

(d) Meal delivery with transshipment

Figure 1: Comparison between transshipment in line-haul logistics and meal delivery.

model meal delivery operations. We estimate two key metrics-deliverers' Vehicle Miles Traveled (VMT) and customers' Average Waiting Time-both with and without transshipment.

Compared to traditional optimization-based models used to address the vehicle routing problem with pickup and delivery (VRPPD) at the operational level (Toth and Vigo, 2002), CA models are extensively applied to tackle big-picture questions about logistics and transportation systems (Ansari et al., 2018). These applications include evaluating new designs in public transit (Chen and Nie, 2017a,b) and comparing new strategies with existing ones for dial-a-ride services (Xu et al., 2020; Ouyang and Yang, 2023). While optimization-based approaches have significantly improved in solution quality over the past decade, their heavy formulations may fall short in providing a comprehensive understanding of the problem's inherent characteristics, which could be critical for managerial insights. Additionally, concerns about data quality, particularly regarding precision and uncertainty, cast doubt on the reliability of optimal solutions in practical applications (Daganzo, 1987b). CA models, serving as a complement to optimization models, have proven effective in addressing these challenges. According to Daganzo et al. (2012), effective CA models can yield five major benefits: they require fewer data inputs, exhibit reduced computational complexity, provide enhanced system representation, improve transparency, and facilitate deeper insights.

For adopting the CA approach, the challenge lies in describing the system of interest as simply as possible while maintaining accurate estimations, a principle known as creating a parsimonious model. For a meal delivery system, its market conditions can be characterized primarily by three exogenous components: the coverage of service areas, the level of customer demand, and the service capacity provided by deliverers. Given specific market conditions, the platform

will make strategic decisions on how to divide the area into smaller sub-areas as well as the strategy for batching and fulfilling pickup and drop-off requests within each sub-area. The system performances are influenced by both the exogenous market conditions and the operating strategy taken by the platform in fulfilling meal deliveries. With the help of CA models, we can address strategic-level questions to identify the market conditions under which introducing transshipment can be advantageous for meal delivery compared to the standard VRPPD strategy.

The remainder of this paper is organized as follows. Section 2 first reviews related work to position the contributions of our study to the literature. Section 3 introduces the design of our proposed strategy for meal delivery with transshipment, followed by a detailed explanation of the derivations to estimate two key system performance metrics: vehicle miles traveled and customer waiting time. Then, Section 4 derives an approximation for VRPPD, tailored to the specific characteristics of meal delivery, which serves as a benchmark for comparison. Comprehensive numerical experiments are conducted in Section 5 to validate the accuracy of our approximations and to compare the performance of the proposed transshipment strategy against the prevalent VRPPD approach. Finally, Section 6 summarizes the findings and concludes the paper.

## 2 Literature Review

In the design of freight distribution systems, there are two main classes based on the presence of transshipment: systems with transshipment and those without. The proposed integrated meal delivery system with a microhub falls into the former category. On the other hand, within each partitioned sub-area, the operations of deliverers align with pickup-and-delivery problems (PDPs) without transshipment.

## 2.1 Distribution system without transshipment

Most studies on distribution systems without transshipment focus on identifying the most costeffective routes to service demand nodes within a network, such as PDPs. PDPs constitute an important family of urban logistics, where goods are transported from various origins to corresponding destinations. The PDP can be categorized into three primary groups based on the nature of the origin-destination pairs: one-to-one problems, where requests start and end at one origin and one destination; one-to-many problems, involving the transport of items from one origin to multiple destinations; and many-to-many problems, where multiple entities (or commodities) are transported between numerous origins and destinations (Daganzo, 2005). The one-to-one distribution problem primarily aims to optimize the one-dimensional dispatching headways along the timeline in response to dynamic demand, which is traditionally addressed through dynamic programming (Newell, 1971). The one-to-many and many-to-many problems are grounded in classic routing problems, such as the Traveling Salesman Problem (TSP) and the Vehicle Routing Problem (VRP), and they vary based on the number of vehicles deployed.

To estimate the minimum total traveling distance for one-to-many TSP, CA models suggest that the optimal distance can be approximated as dTSP ≈ kTSP √ AN , where A represents the area of service region, and N denotes the number of destinations to visit within this area. Through a swath heuristic, which splits the service area into several strips and then determines the optimal swath width, Daganzo (1984b) estimated the constant factor kTSP is 0.9 for Euclidean metric and 1.15 for L 1 metric if the service zone is not too narrow. Different from the swath-split strategy, del Castillo (1999) constructs the suboptimal tour via a ring-radial separation. Their

numerical experiments demonstrate that the average Euclidean length of the tours produced by this approach scales with √ AN , showing performance essentially on par with Daganzo's swath heuristic. Building on the TSP approximation, the multiple-vehicle delivery problem, or capacitated vehicle routing problem, typically assumes each vehicle can only deliver to a limited number of points. The expected all-vehicle travel distance is expressed similarly to the TSP as dVRP ≈ 2 r ( s / S N ) + kVRP √ AN , where r represents the distance from the depot to the service area's center, s and S are demand expectation and vehicle capacity, respectively (Daganzo, 1984a). This formula, incorporating a 'line-haul distance' for initial travel from the depot to the customers and a 'local distance' for detours within each region, is also extended to a ring-radial network and other rectangular zones with various orientations (Daganzo and Newell, 1986a; Newell and Daganzo, 1986).

The many-to-many distribution problem, usually referred to as the multi-commodity problem, assumes each demand has a specific origin and destination. The many-to-many distribution system can be divided into two categories: VRPPD and dial-a-ride problem (DARP), where the former considers the transportation of goods, and the latter focuses on passenger transportation. The key difference is that DARP usually incorporates an extra hard or soft constraint to specifically address and minimize customers' inconvenience (Toth and Vigo, 2002). Daganzo (1978) initially formulated a heuristic strategy designed to estimate customers' waiting and riding times, which can be represented as functions of the number of vehicles, request density, service area, and vehicle traveling speed. The optimal approach advocates for alternating pickups and deliveries, consistently opting for the nearest pick-up or delivery point. The effectiveness of this approach was illustrated through a comparison against two alternative strategies: (i) the bus routes to the nearest feasible point, and (ii) the bus collects a fixed number of passengers before proceeding with deliveries. To estimate the total service time, several simulation models (Wilson et al., 1969, 1976) and empirical models (Flushberg Sr and Wilson, 1976) have also been developed for the many-to-many DARP service. Stein (1978) proved that the many-to-many TSP with pickup and delivery constraints (TSPPD), adheres to a similar formula as the traditional TSP with dTSPPD ≈ kTSPPD √ AN , where √ 2 kTSP ≤ kTSPPD ≤ 2 kTSP . Through a probabilistic analysis, Stein split the optimal tour into several segments and demonstrated that the optimal kTSPPD asymptotically converges to 4 √ 2/3. In this study, we adopt the methodology from Daganzo (1978) by incorporating additional practical considerations specific to meal delivery to estimate deliverers' VMT and customers' waiting time under the VRPPD strategy. These estimations are then used to benchmark the performance of a distribution system that implements transshipment.

## 2.2 Distribution system with transshipment

In one-to-many or many-to-many distribution, transshipment can be viewed as a terminal sorting, storing, and transferring items from one vehicle to another. Incorporating transshipment can significantly reduce the total miles traveled by vehicles. But at the same time, it may introduce new handling and holding costs at the terminal, attributable to the increased batch sizes. Designing the distribution system with transshipment typically involves determining the number of terminals to be utilized, their strategic locations, the planning of routes and schedules for various types of vehicles, and the assignment of customers to specific terminals and routes (Daganzo, 2005). Daganzo and Newell (1986b) first employed a CA method to find near-optimal solutions for the design guidelines in a one-to-many distribution system, such as the location of transshipping points, their areas of influence, vehicle frequency, and fleet size. Smilowitz and Daganzo (2007) then advanced the design strategy for an integrated package distribution system consid-

ering service level, i.e., overnight and longer deadlines. They separated the complex system into a series of subproblems to determine the densities of various terminal levels (e.g., consolidation terminal, breakbulk terminal, airport) and their corresponding dispatching headways through a CA method. Extensive studies have employed the CA method to design distribution systems or identify the optimal location for transshipment across various application contexts, such as delivery consolidation in freight transportation (Campbell, 2013; Xie and Ouyang, 2015; Ghaffarinasab et al., 2018), inventory management within supply chain (Naseraldin and Herer, 2011; Tsao et al., 2012), and agricultural collection during harvest season (Wiles and van Brunt, 2001).

To evaluate the benefit of the transshipment design, Daganzo (1987a) first delved into the 'break-bulk' role of the transshipment in many-to-many logistic networks using CA. The investigation highlighted the cost-saving benefits of allowing one and then two transshipment stops, attributing these savings to the reduction in line-haul distance traveled. They offset the additional costs associated with transshipment, including delivery delays and handling costs at the terminal. In the modern logistics network with multiple suppliers and customers, numerous studies have employed CA to evaluate the effects of delivery consolidation. To reduce transportation costs in urban delivery, companies collaborate to consolidate shipments at a central terminal rather than independently dispatching goods from the distribution centers to customers. Kawamura and Lu (2007) analyzed the economic efficiency of cooperative delivery consolidation based on a logistics cost analysis and concluded that, without accounting for societal benefits such as reduced illegal parking and congestion, the appeal of consolidation to the U.S. industry is marginal. Chen et al. (2012) extended this work by considering a different network setting for small businesses and using an additional dataset for parameter estimation. Their study theoretically demonstrated that delivery consolidation becomes cost-effective when the number of customers and suppliers is large, and the terminal operation cost is low. Lin et al. (2016) then expanded their CA model by incorporating energy consumption and PM2.5 emissions as additional costs, along with policy factors like commercial vehicle size restriction in urban areas. Their findings indicate that when there is an economy of scale or high customer density, the consolidation design yields both logistics and environmental benefits. Similarly, Pahwa and Jaller (2022) developed a multi-echelon last-mile distribution model and suggested that consolidation strategies are well-suited for areas with dense demand. In contrast, for regions with lower population density, outsourcing solutions such as crowdsourced delivery and customer self-collection emerge as more cost-effective alternatives.

In urban logistics, existing research on transshipment has mostly focused on the design of consolidation centers for logistic companies, with an emphasis on the operational costs, such as handling and fixed pipeline inventory expenses associated with additional transshipment. However, these costs are relatively minor when applied to the temporary storage of meal packages awaiting transshipment. The effectiveness of transshipment strategies for meal delivery services, as well as the optimal system design, has not yet been thoroughly investigated and remains unclear. We thus propose to examine the cost-effectiveness of transshipment specifically for meal delivery from the perspectives of both deliverers and customers. In this study, we develop stylized models for meal delivery operations, both with and without transshipment, to compare the efficiency of the two operational models.

## 3 Integrated Meal Delivery System with Transshipment

This section introduces our proposed system design for transshipment and key performance metrics related to different stakeholders. It begins with an overview of the system setup for the transshipment operations, followed by a breakdown of the lifecycle for meal orders from creation to package fulfillment, from their creation to final package delivery. Given the system inputs, continuous approximations are applied to estimate the average vehicle miles traveled by deliverers, whilst a queuing network is conceptualized to analyze the time that meal packages spend at various delivery stages.

## 3.1 Meal delivery with transshipment

Consider a meal delivery service offered in a region of area A , with an average order arrival flux λ (i.e., order arrivals per hour per square mile), a total of m deliverers. For the transshipment operations illustrated in Figure 1d, the region is divided into K pie-shaped sub-areas, and each deliverer is assigned to service within a specific sub-area. Each sub-area k ∈ K is characterized by a coverage area Ak , an order arrival rate λ k (= λ Ak ) , and a dedicated fleet of mk deliverers. Any regional partition defined by { Ak , mk } k ∈ K must satisfy the constraints ∑ k mk = m and ∑ k Ak = A . Given a partition of sub-areas, individual delivery orders can be categorized as either intra-zonal or inter-zonal. Intra-zonal orders have both their pickup and drop-off locations within the same sub-area, whereas inter-zonal orders involve pickup and drop-off locations spanning different sub-areas.

Deliverers operate cyclically and exclusively within their designated sub-areas to perform pickup and drop-off tasks. In each sub-area k , a deliverer is dispatched once there are nk points accumulated to visit. These nk points comprise both the meals waiting for drop-off at the microhub and the meals awaiting pickup at restaurants. Upon dispatch, the deliverer completes a round trip that begins and ends at the microhub, touring through the nk pending pickup and drop-off points. Regardless of whether the orders are intra-zonal or inter-zonal, all meal packages are first picked up from their originating restaurants and transported to the microhub during the current delivery cycle. At the microhub, they are sorted according to their destinations and await drop-off. Then, in a subsequent cycle, one deliverer assigned to the sub-area containing an order's drop-off point retrieves the package from the microhub and delivers it to its final destination. As a result, the distribution tasks in each delivery cycle within a sub-area closely resemble a one-to-many-to-one (1-M-1) pickup and delivery problem, involving two distinct sets of commodities: some are transported from the microhub to customers, while others are picked up at restaurants and returned to the microhub (Toth and Vigo, 2002).

## 3.2 Lifetime of meal delivery orders

Figure 2 outlines the lifecycle of a delivery order for meal packages under the proposed transshipment strategies. Each delivery order progresses through four sequential stages, defined as follows:

- 1. Meal Preparation : Interval from when an order is placed until it is ready for pickup;
- 2. Pickup : Interval from the meal package being ready to its arrival at the microhub;

Figure 2: Lifetime of individual meal delivery orders

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.01047v2-figure-28ea3098-1d2b-4634-8818-16f3d9be985b.jpg)

- 3. Transfer : Interval during which the package is held at the microhub for the next delivery tour;
- 4. Drop-off : Interval for transporting the package from the microhub to the customer.

Of these, the meal preparation stage is unaffected by the introduction of transshipment and will therefore be left out in the following analysis. The pickup stage can be further broken down into three phases: the period from when the package is ready for pickup until a deliverer is dispatched, the interval between the dispatch of the deliverer and the actual pickup of the order, and the duration required for the package to be transported from the restaurant to the microhub. The combined duration of the latter two phases indeed constitutes the total time for a complete delivery tour executed by a deliverer. Additionally, it should be noted that when a meal is ready for pickup, it does not guarantee the immediate dispatch of a deliverer. Two conditions must be met within each sub-area: first, nk pending pickup and drop-off tasks must accumulate in sub-area k ; second, after nk tasks have been accumulated, the batch must wait until a deliverer assigned to sub-area k becomes available to handle the service. Due to the presence of interzonal orders, the operations across different sub-areas are closely interdependent. To capture these cross-zonal interactions, a queuing network with K dependent queues is proposed to pin down the waiting time for pickup as well as the transfer time, which can also be interpreted as the waiting time for drop-off.

## 3.3 Queuing network model

The K sub-areas collectively form a Jackson queuing network in Figure 3, where each sub-area operates as a separate queue with its own servers. The following assumptions are made for this network of K single-station queues:

Figure 3: Cross-zonal meal delivery process as a queuing network

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.01047v2-figure-8b7e1215-a45e-4735-af34-6e5da661c3a3.jpg)

- · New orders arrive at sub-area k following a Poisson process with an arrival rate of λ k . The pickup point for each order is randomly located within sub-area k .
- · A package picked up from sub-area i , after being transferred to the microhub, will join the queue at sub-area k for drop-off with probability pi k , . The drop-off points of orders are also random within each sub-area k . Once a package is dropped off, it exits the queuing network completely.
- · The k -th sub-area is served by mk delivery vehicles. When the number of orders in the queue reaches nk or more and a deliverer becomes idle, a batch of exactly nk orders is taken for service.
- · The service time Sk for meal packages at sub-area k are independent and identically distributed, with a mean E [ Sk ] and a variance Var [ Sk ] .
- · The microhub has sufficient capacity to accommodate all incoming packages, and the queue at each sub-area has unlimited waiting spots.

The arrival of meal orders in the queue for each sub-area can be viewed as a Poisson process with an arrival rate δ k , which sums the arrival rates of newly generated orders in sub-area k for pickup and the collected orders from all sub-areas that need to be dropped off to sub-area k , i.e.,

$$\delta _ { k } = \lambda _ { k } + \sum _ { i = 1 } ^ { K } \lambda _ { i } p _ { i, k }.$$

The queued orders are serviced in batches of size nk . For the i -th order arriving before reaching the nk threshold in a single batch, its waiting time is given by ( nk -i ) / δ k . Consequently, the average waiting time W a k of an arbitrary pickup or drop-off order in sub-area k , from its arrival

until the accumulation of nk orders for a batch, is calculated as follows,

$$W _ { k } ^ { a } = \frac { 1 } { n _ { k } } \sum _ { i = 1 } ^ { n _ { k } } \frac { n _ { k } - i } { \delta _ { k } } = \frac { n _ { k } - 1 } { 2 \delta _ { k } }.$$

After accumulating a sufficient number of orders for a batch, the batched orders may also experience additional waiting time if a deliverer is not immediately available. For each sub-area, this process of forming and delivering batch orders by deliverers can be modeled as a G G m / / queue. The inter-arrival time of consecutive order batches follows an Erlang distribution with the shape parameter nk and the rate parameter δ k . The service rate by a deliverer is determined by the touring time required to complete a Traveling Salesman Problem (TSP) for fulfilling the batched orders. The number of servers in each queue equals the number of deliverers mk . Thus, the additional queuing time for batched orders W q k can be estimated with a widely used approximation for a G G m / / queue by Larson and Odoni (1981), presented as follows,

$$W _ { k } ^ { q } \approx \frac { \left [ \frac { n } { \partial _ { k } } + \frac { 1 } { m _ { k } } \, \text{Var} [ S _ { k } ] \right ] \cdot \frac { \delta _ { k } } { n _ { k } } } { 2 \left ( 1 - \frac { \delta _ { k } \, \mathbb { E } [ S _ { k } ] } { n _ { k } \, \cdot \, m _ { k } } \right ) }.$$

## 3.4 Service time approximations

The service process for a single queue follows a TSP tour undertaken by deliverers. Daganzo (1984b) showed that, in a regularly shaped area of size A with N randomly distributed points to visit, the expected TSP tour length is proportional to √ AN . To adapt this formulation for the pie-shaped sub-areas, we introduce an extra distance R ′ to account for radial forward and return travel between the microhub and the farthest points within each pie-shaped region. Proposition 1 derives the expected value of R ′ .

Proposition 1. For a circular sector of radius R containing N nodes, the expected farthest distance, E [ R , ′ ] is given by E [ R ′ ] = 2 N 2 N + 1 R, where E [ R ′ ] is independent of the sector's angle.

This formula suggests that when there are fewer points to visit, the radial traveling distance often remains short, as deliverers likely do not need to venture far from the microhub. However, as the number of points increases, deliverers tend to travel closer to the outer perimeter of the pie-shaped area. By incorporating this radial distance, the expected service time in sub-area k under the modified TSP approximation is expressed as follows:

$$\mathbb { E } [ S _ { k } ] = \frac { 1 } { v } \underset { \wedge \,. } { \cdot } ( a \cdot \sqrt { A _ { k } n _ { k } } + b \cdot \mathbb { E } [ R ^ { \prime } _ { k } ] ),$$

$$\mathbb { E } [ R ^ { \prime } _ { k } ] = \frac { \mu } { 2 n _ { k } + 1 } R _ { k },$$

where v denotes the speed of delivery vehicles, and Rk is the radius of sub-area k , which approximates the average radial distance from the vertex to the perimeter border of the pie-shaped area. It is worth noting that the origin-destination pairs of individual meal delivery orders may be correlated according to certain distance profiles, implying that pickup and drop-off locations are not entirely random within the service region. However, in the transshipment system, each order's pickup and drop-off points are separated into two distinct delivery tours. As a result, all pickup and drop-off points within a single tour correspond to different orders. Consequently,

it is reasonable to assume that the access points visited by individual deliverers within each sub-area are randomly distributed in space.

For the variance of TSP touring distances Var [ Sk ] , it is intuitive that its value is linearly proportional to the size Ak of the sub-area. Regarding its relationship with the number of points nk , we observe that when nk is small, the spatial dispersion of points can be highly uncertain, with points either clustering together or spreading across the entire region. Conversely, as nk increases, the points tend to distribute more uniformly throughout the region, leading to a more stable variance. Thus, Var [ Sk ] could be inversely proportional to nk and converge to a certain positive value as nk approaches infinity. Based on the above deduction, we propose the following hypothetical form for describing the service time variance under the TSP:

$$\text{Var} [ S _ { k } ] = \frac { A _ { k } } { v ^ { 2 } } \cdot \left ( \frac { \alpha } { n _ { k } } + \beta \right ).$$

The hyperparameters a b , and α β , in Eqs. (4) and (6) will be calibrated later through simulation, where the effectiveness of the proposed approximations will also be demonstrated.

## 3.5 System performance metrics

Given a regional partition { Ak , mk } k ∈ K and the order batching sizes { nk } k ∈ K , the queuing network model predicts the stationary states of the meal delivery system under transshipment and subsequently yields the two key system performance metrics.

Deliverers' vehicle miles traveled. For each sub-area, when nk orders are collected, one deliverer is dispatched to visit these nk points, while other idle deliverers, if any, wait near the microhub waiting for the next delivery. Thus, the expected VMT per hour by all deliverers in sub-area k , denoted as Qk , is given by the product of the number of delivery cycles performed per hour and the expected traveling distance for visiting nk points. Specifically,

$$Q _ { k } = \frac { \delta _ { k } } { n _ { k } } \cdot \mathbb { E } [ S _ { k } ] \cdot v.$$

Customers' waiting time. As previously noted, each customer's waiting time can be divided into four components, with the meal preparation time excluded from this analysis. Combining the rest three components, the expected total waiting time for meal orders delivered from subarea i to k , denoted as Wi k , , can be expressed as:

$$W _ { i, k } = W _ { i } ^ { a } + W _ { i } ^ { q } + \mathbb { E } [ S _ { i } ] + W _ { k } ^ { a } + W _ { k } ^ { a } + \frac { \text{Var} [ S _ { k } ] + \mathbb { E } ^ { 2 } [ S _ { k } ] } { 2 \, \mathbb { E } [ S _ { k } ] },$$

where W a i + W q i + E [ Si ] accounts for the pickup time, covering the period from order generation to arrival at the microhub; W a k + W q k represents the transfer time at the microhub; and the last term in Eq. (8) stands for the drop-off time. Specifically, the time for package drop-off is derived based on the occurrence of a random event along a TSP tour, where the touring time follows a given distribution with a known mean and variance (Larson and Odoni, 1981).

## 4 A VRPPD Benchmark without Transshipment

For the benchmark operations without transshipment, individual deliverers continuously traverse the service area, picking up meal packages and delivering them directly to their destina-

tions. In this section, we adopt a heuristic for the VRPPD, originally proposed by Daganzo (1978), to estimate system performance without transshipment. The heuristic is further enhanced by incorporating practical considerations specific to the delivery distances of individual meal orders.

## 4.1 A heuristic for meal delivery without transshipment

The heuristic devised by Daganzo (1978) directs deliverers to continuously route to the nearest point pending a visit, which can be either the origin of a pickup order or the destination of a drop-off order already on the vehicle. At equilibrium, a deliverer is equally likely to visit an origin/pickup point or a destination/drop-off point. By assuming that both origins and destinations are independently and randomly distributed across the service area, Daganzo (1978) showed that the trip chains continue with equal probabilities, and the number of packages on each delivery vehicle simply equals the number of packages waiting for pickup. Unlike the setting in Daganzo (1978), distance-based delivery fees and platform-imposed limitations on ordering distances create strong bundling between origins and destinations of meal orders, shaping delivery distances into specific profiles. The empirical data from Meituan calibrates the distance profile as a Rayleigh distribution, with the goodness-of-fit detailed in Appendix B.

Accounting for the distance bundling between origin-destination pairs of meal delivery results in more sophisticated state transition probabilities for deliverers' trip chains, as illustrated in Figure 4. Let NO denote the total number of packages waiting for pickup and ND the number of packages en-route on a delivery vehicle. Then, immediately following a visit to a drop-off point Dpre , shown in Figure 4a, the deliverer has a probability of NO NO + ND to visit a pickup point Opost ′ and ND NO + ND to visit another drop-off point Dpost ′ . These probabilities arise from the assumption that all meal orders are generated i.i.d. across the region according to a spatial Poisson process. On the other hand, following a visit to a pickup point Opre , Figure 4b shows that the deliverer's next visit falls into either of the two cases: (1) the corresponding drop-off point Dpost , with probability γ , or (2) any other access point within the service region, with probability 1 -γ . The latter case can be further decomposed into visits to either another pickup point Opost ′ or a dropoff point Dpost ′ that is different from the destination of the order just picked up. The probability γ represents the likelihood that, after visiting a pickup point, the distance to its corresponding

Figure 4: Scenario enumeration of state transitions between the types of preceding and succeeding points visited consecutively by individual deliverers.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.01047v2-figure-08122787-bcec-41bc-b729-0035a9b5db23.jpg)

drop-off point is shorter than to any other point to visit. This probability can be derived from Proposition 2.

Proposition 2. Let X be a random variable representing the distance to the nearest point that adheres to a spatial Poisson process with density NO + ND . Let Y be a random variable following a Rayleigh distribution π R 2 with a scale parameter σ &gt; 0 . Then, the probability that X greater than Y is given by:

$$\gamma = I P ( X > Y ) = \frac { 1 } { 1 + \frac { 2 ( N _ { O } + N _ { D } ) \sigma ^ { 2 } } { R ^ { 2 } } }.$$

As suggested by Eq. (9), if the delivery orders are shorter, with a smaller σ , the probability of directly visiting the drop-off point immediately after picking up an order without detours increases. Conversely, with denser visiting points in space, the likelihood of immediately dropping off an order following its pickup decreases.

The state transitions spelled out in Figure 4 lead to the following relationship based on the Law of total probability:

$$\mathbb { P } ( O _ { p o s t } ) = \mathbb { P } ( D _ { p r e } ) \cdot \mathbb { P } ( O _ { p o s t } | D _ { p r e } ) + \mathbb { P } ( O _ { p r e } ) \cdot \mathbb { P } ( O _ { p o s t } | O _ { p r e } ),$$

$$\Leftrightarrow \quad \mathbb { P } ( O _ { p o s t } ) = \mathbb { P } ( D _ { p r e } ) \cdot \frac { \dot { N } _ { O } } { N _ { O } + N _ { D } } + \mathbb { P } ( O _ { p r e } ) \cdot ( 1 - \gamma ) \frac { \dot { N } _ { O } } { N _ { O } + N _ { D } }.$$

Moreover, the equilibrium condition implies that I P ( Opre ) = I P ( Dpre ) = I P ( Opost ) = I P ( Dpost ) = 50%, which then gives rise to the following equation:

$$\gamma = \frac { N _ { O } - N _ { D } } { N _ { O } }.$$

## 4.2 Service rate estimations

The service rate µ provided by individual deliverers is given by:

$$\mu = \left ( \frac { 2 \, \mathbb { E } [ d ] } { \upsilon } \right ) ^ { - 1 },$$

where d is the distance from an arbitrary visit point to its nearest point. To derive this, we analyze the two scenarios separately to determine the distances dD and dO , which respectively denote the distance to the nearest point after visiting a drop-off and a pickup point, given NO + ND points to visit within the service area.

First, with the intensity of the spatial Poisson process given by δ = NO + ND π R 2 , the cumulative distribution function (CDF) and probability density function (PDF) of dD are FdD ( r ) = 1 - - · e δ π r 2 and f dD ( r ) = 2 πδ r · e - · δ π r 2 . The expected value of dD is thus,

$$\alpha _ { \ } J d _ { D } \vee \nu - \nu \nu \cdot \mathfrak { c } \quad \cdot \, \min \, \epsilon \, \cos \, \nu \, \cos \, \nu \, \cos \, \nu \, \cos \, \nu \, \cos \, \nu \, \cos \, \nu \, \cos \, \nu \, \cos \, \nu \, \cos \, \nu \, \cos \, \nu \, \cos \, \nu \, \cos \, \nu \, \cos \, \nu \, \cos \, \nu \, \cos \, \nu \, \cos \, \nu \, \cos \, \nu \, \cos \, \nu \, \cos \, \nu \, \cos \, \nu \, \cos \, \nu \, \cos \, \nu \, \cos \, \nu \, \cos \, \nu \, \cos \, \nu \, \cos \, \varphi \, \cos \, \varphi \, \cos \, \varphi \, \cos \, \varphi \, \cos \, \varphi \, \cos \, \varphi \, \cos \, \varphi \, \cos \, \varphi \, \cos \, \varphi \, \cos \, \varphi \, \cos \, \varphi \, \cos \, \varphi \, \cos \, \varphi \, \cos \, \varphi \, \cos \, \varphi \, \cos \, \varphi \, \cos \, \varphi \, \cos \, \varphi \, \cos \, \varphi \, \cos \, \varphi \, \cos \, \varphi \, \cos \, \varphi \, \cos \, \varphi \colon \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \pi \, \cos \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \alpha \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \gamma \, \cos \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \beta \, \gamma \,.$$

$$= \frac { \sqrt { \pi } R } { 2 \sqrt { N _ { O } + N _ { D } } }.$$

Specifically, the first equality follows from the definition, the second applies the integration by parts, the third uses the Gaussian integral formula, and the final equality substitutes the corresponding expression for δ .

Second, immediately after visiting a pickup point, the minimum distance to the next point is determined by the closest point between its corresponding drop-off point and other existing visit points, as derived in Proposition 3.

Proposition 3. As a continuation of Proposition 2, the expected value of the minimum of X and Y is given by E [ min ( X Y , )] = 1 2 √ π N O + N D R 2 + 1 2 σ 2 .

Combining the two cases, the average distance d from an arbitrary visit point to its nearest point is obtained as follows,

$$\cdot \quad \mathbb { E } [ d ] = \mathbb { P } ( D _ { p r e } ) \cdot \mathbb { E } [ d _ { D } ] + \mathbb { P } ( O _ { p r e } ) \cdot \mathbb { E } [ d _ { O } ] = \frac { \sqrt { \pi } R } { 4 \sqrt { N _ { O } + N _ { D } } } + \frac { 1 } { 4 } \sqrt { \frac { \pi } { N _ { O } + N _ { D } } + \frac { 1 } { 2 \sigma ^ { 2 } } }.$$

## 4.3 System performance metrics

When the demand for meal delivery is low, it is unnecessary to have all deliverers actively in service. Therefore, we define an additional variable m ′ as the total number of deliverers in service, which should not exceed the total number of available deliverers m . Additionally, system equilibrium requires that the arrival rate ( λ A ) equals the total service rate ( µ m ′ ), i.e.,

$$\lambda A = \mu m ^ { \prime }.$$

Combining all the components above, given the demand flux λ and the number of deliverers in service m ′ , the equilibrium state of the VRPPD system without transshipment can be determined by solving the system of equations (9)-(14).

Deliverers' vehicle miles traveled . Similar to Eq. (7), the total VMT by deliverers is given by the product of the number of points visited by a deliverer per hour and the expected distance between two adjacent visits:

$$Q _ { V R P P D } = 2 \lambda A \cdot \mathbb { E } [ d ] = m ^ { \prime } \cdot v.$$

Customers' waiting time . The average waiting time for customers, including both the waiting time for package pickup and the riding time on a vehicle, is given by:

$$W _ { V R P P D } = \frac { N _ { O } } { \lambda A } + \frac { N _ { D } } { \mu }.$$

## 5 Numerical Analysis

We then evaluate the potential of the proposed meal delivery operations with transshipment by benchmarking them against the traditional pickup-and-delivery model without transshipment.

First, we set the stage for numerical experiments by formulating two optimal system designs to establish the performance envelope under different market conditions. Next, we run simulations to validate the accuracy of our analytical model in predicting average customer waiting time and VMT by deliverers. A comprehensive performance evaluation is then conducted to compare meal delivery systems with and without implementing transshipment. Finally, we introduce a dataset from the Meituan meal delivery platform to compare the performance of these strategies in a real-world scenario.

## 5.1 Experimental setup

To ensure a fair head-to-head comparison of operations using the proposed CA models, this evaluation considers the context of an isotropic circular area of radius R with an average order arrival flux λ and a total of m deliverers. For each market condition, we optimize the system design for both transshipment and non-transshipment operations and compare the resulting system performance.

For transshipment operations, the experiments in this section focus on symmetric partitioning strategies, where the area is evenly divided into K sub-areas, with deliverers uniformly distributed across them. For brevity, we retain the notations of variables previously defined for individual sub-areas with the subscript k . However, since all sub-areas are now identical, these notations will now refer to a representative sub-area rather than a specific one. The optimal design under symmetric partitioning reduces to the following optimization program, which determines the optimal number of partition K and the batch size nk per delivery cycle in sub-areas:

$$& \min _ { K, n _ { k } \geq 0 } K \cdot ( \pi _ { Q } \cdot Q _ { k } + \pi _ { W } \cdot \lambda _ { k } W _ { k, k } ) \\ s. t. & \quad ( 2 ) - ( 8 ) \\ & A _ { k } = \frac { A } { K }, \ \lambda _ { k } = \lambda A _ { k }, \ \delta _ { k } = 2 \lambda _ { k }, \ m _ { k } = \frac { m } { K } \\ & \delta _ { k } \leq \frac { m _ { k } } { \mathbb { E } [ S _ { k } ] / n _ { k } }$$

$$\delta _ { k } \leq \frac { m _ { k } } { \mathbb { E } [ S _ { k } ] / n _ { k } }$$

The objective (17a) represents the generalized variable costs in the meal delivery system, encompassing both the operational costs incurred by deliverers and the time costs experienced by customers, weighted by the coefficients π Q and π W . The equalities in (17b) specify the allocation of demand and supply within the sub-areas according to the partition. Finally, the condition enforced by (17c) ensures that the meal delivery system with transshipment remains in steady operations.

For the non-transshipment benchmark modeled as a VRPPD, the optimal system design seeks to determine the number of deliverers actively participating in the service m ′ , along with the total number of packages waiting for pickup NO and the number of packages en route on a delivery vehicle ND . The corresponding optimization problem is formulated as:

$$\min _ { m ^ { \prime }, N _ { O }, N _ { D } \geq 0 } \pi _ { Q } \cdot Q _ { V R P P D } + \pi _ { W } \cdot \lambda A \cdot W _ { V R P P D }$$

$$\begin{matrix} \text{s.t.} & ( 9 ) - ( 1 6 ) \\ & m ^ { \prime } \leq m \end{matrix}$$

where constraint (18b) ensures that the number of active deliverers in service does not exceed the total available deliverers.

## 5.2 Simulation and measurement

Simulations are conducted to demonstrate the effectiveness of proposed approximations for TSP tours-Eqs. (4), (6)-and for customer waiting time-Eqs. (2), (3), (8). The simulated results are compared against the predicted approximations across a range of scenarios, varying the values of { λ , R m n K , , , } , to measure the approximation accuracy. Using the trajectory data from Meituan deliverers, the average touring speed for meal delivery is calibrated to v = 4.15 miles per hour, accounting for delays caused by traffic lights and the time required to access restaurants and customer drop-off points.

## Measurement of approximations for deliverers' TSP tours

We vary the values of R , n , and K and simulate 1,000 trips for each case, where n access points are randomly distributed across a pie-shaped circular sector with radius R and angle θ = 2 π / K . The simulated results are obtained through the built-in solvers provided by Google OR-Tool's vehicle routing package, and we then calculate the average and variance of travel

(a) Mean w.r.t. R

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.01047v2-figure-fd2d0816-beaf-4375-9268-c68c09bedc3c.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.01047v2-figure-6bec58e4-5cb0-4a5b-ab7f-038916b1feaf.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.01047v2-figure-59c13440-f0e9-4cd0-9d02-42d119fb8262.jpg)

(b) Mean w.r.t. n

(c) Mean w.r.t. K

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.01047v2-figure-df19ae1b-8d0a-47cd-9c0e-61fe0ab2d274.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.01047v2-figure-63fb42b8-6b7b-4564-b789-223ea3961f45.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.01047v2-figure-5b403012-5419-4808-9894-f3c31d27910a.jpg)

(f) Var w.r.t.

K

Figure 5: Comparison between measured and predicted mean and variance of the traveling distance in TSP tours w.r.t. 1) radius size R , 2) number of points n , and 3) number of divisions K .

distance across the 1,000 simulations for each combination of R , n , and K . We calibrate Eq. (4) and Eq. (6) for E [ Sk ] and Var [ Sk ] using the simulation results and obtain the parameter values a = 0.64, b = 1.28, α = 0.42, and β = 0.07. Figure 5 illustrates the fitting curves for the expected value and variance of the traveling distance in TSP tours. The results indicate that the approximations using Eqs. (4) and (6) achieve R 2 values of 98.64% and 85.98%, respectively, when compared to the optimal solutions provided by Google OR-Tool.

## Measurement of approximations for customers' waiting time

The 'ground truth' arrivals of delivery requests and customers' waiting time are simulated over a 6-hour period to experiment with the order fulfillment process. The arrival of requests follows a spatio-temporal Poisson process with a given flux λ . Once a sub-area accumulates n orders and there is a deliverer available, one deliverer is dispatched to visit the n access points following the routing provided by the Google OR-Tool. The simulation platform records the critical time intervals for each order, including the total waiting time from order generation to delivery, along with the three decomposed components described in Section 3.2. To ensure steady-state conditions, the first hour is treated as a warm-up period and excluded from the analysis. The average waiting time of all orders is calculated over the remaining five hours. The baseline parameters are set as follows: order arrival flux λ = 50 orders/hour mi 2 , fleet size · m = 100, area radius R = 1.5 miles, batch size n = 10 orders, and number of sub-area partition K = 4. For each plot in Figure 6, four out of the five parameters are held constant, while the total and three decomposed components of customer waiting time are plotted as functions of the varying parameter.

Figures 6a and 6c exhibit a non-monotonic trend in customer waiting time as a function of demand intensity and service area size, respectively. As the demand intensity λ or service radius R increases, customer waiting time initially decreases and then starts to rise. The initial decrease occurs because, when the order arrival rate is low or the service area is small, it takes longer to accumulate n orders for batch service, resulting in excessively long customer wait times. Increasing the demand intensity or service area size raises the order arrival rate, thereby reducing the holding time in delivery. However, in the increasing segment, when the arrival rate becomes sufficiently high and approaches the service capacity of the available fleet, the system becomes congested, leading to longer waiting times. In addition to demand-side influences, Figure 6b illustrates that a shortage of deliverers also results in system congestion and prolonged customer waiting times. However, once the fleet size reaches a sufficient level, adding more deliverers does not yield noticeable improvements, as many will remain idle.

Interestingly, Figure 6d demonstrates that a higher batch size of orders n consistently results in longer customer waiting times due to the increased touring time per delivery cycle. However, in terms of VMT, as shown in Eq. (7), batching more orders per delivery tour reduces the overall VMT by deliverers, highlighting a clear trade-off in system design. Regarding the number of partitioned sub-areas, Eqs. (2) and (3) indicate that finer partitioning of sub-areas (i.e., higher K ) increases the time required to accumulate n orders within each sub-area but reduces the travel distance per delivery cycle. Consequently, the total waiting time exhibits a non-monotonic trend, as depicted in Figure 6e. These trade-offs over the choices of n and K justify the existence of optimal system designs with moderate configurations, rather than extreme values.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.01047v2-figure-8ce7ae6c-755d-4853-b626-b1d7f3c21417.jpg)

(a)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.01047v2-figure-c1db7c23-dd79-4a46-a5c0-7d891fdd4fef.jpg)

(c)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.01047v2-figure-799f33a6-3393-4024-b46f-08bba2533c66.jpg)

(b)

Figure 6: Comparison between measured and predicted customers' waiting time w.r.t. (a) order arrival flux λ , (b) fleet size m ; (c) area radius R , (d) batch size of orders n , (e) number of sub-area partitions K .

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.01047v2-figure-58e561fb-532e-4e96-bc28-4072c5ae42ba.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.01047v2-figure-51297fc4-b7a5-4ddc-ac70-8fb56e2ce5b6.jpg)

## 5.3 Effectiveness of transshipment vs. non-transshipment

We next compare the performance of meal delivery systems with and without implementing transshipment strategies. Following the values recommended by Nourbakhsh and Ouyang (2012), we set the operational cost per vehicle mile at π Q = $2/veh-mile and customers' value of time at π W = $20/hour. We analyze three scenarios where the radius of the service region is set at R = { 1.2, 1.5, 1.8 } miles, and the order arrival flux λ ranges from 5 to 170 orders/hour mi 2 , re-· spectively. For each scenario, the total number of available deliverers is set to m = { 64, 100, 144 } , maintaining a consistent ratio between the number of orders and deliverers. For all cases, given the area A = π R 2 , fleet size m , and order arrival flux λ , we solve problems (17) and (18) to determine the optimal designs for transshipment and non-transshipment operations, respectively. We then compare the average customer waiting time (minutes/customer) and the average deliverer VMT (miles/hour-vehicle) achieved under transshipment with the corresponding metrics in the non-transshipping VRPPD counterpart.

Figure 7 presents a comparison of results and highlights the conditions under which transshipment demonstrates relative strength over non-transshipment strategies. From Figure 7a, it is evident that while both strategies involve order consolidation onto individual vehicles, the localized batching facilitated by transshipment consistently reduces VMT compared to the VRPPD approach. Without transshipment, when demand exceeds a certain threshold-around λ ≈ 20 in Figure 7a-all deliverers are fully deployed and travel continuously between pickup and drop-off points, causing the average VMT per hour to converge toward the vehicle's operating speed. In contrast, our transshipment design allows deliverers to wait at the microhub until the number of visiting points accumulates to the batching size n , thereby resulting in significant VMT savings.

Furthermore, as shown in Figure 7b, for a given service range R with a fixed supply, the transshipment strategy demonstrates its advantages by achieving lower customer waiting times compared to the VRPPD counterparts as demand increases within the examined region. In Figure 8, we demonstrate the optimal partitioning and batching results obtained from solving the optimal design problem (17). It suggests that, with sufficient service capacity, system efficiency can potentially improve as demand grows by raising the number of sub-areas and the batching

Figure 7: Comparison between the transshipment and non-transshipment strategies on (a) VMT and (b) customer waiting time w.r.t. order arrival flux.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.01047v2-figure-82328e6f-f14e-4128-9107-84bc6344c097.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.01047v2-figure-0631b523-bd3a-422a-bc37-6f1565563175.jpg)

Figure 8: Optimal partitioning and batching strategic decisions for meal delivery system with transshipment w.r.t. order arrival flux.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.01047v2-figure-f904b1bc-a935-41cd-a9e2-14880f343b34.jpg)

size of orders. However, at lower demand levels, transshipment through the microhub becomes less efficient than the direct pickup-and-delivery strategy due to the additional time required for routing and transferring. Notably, Figure 7b also illustrates that the efficiency gains enabled by transshipment are influenced by the size of the service region. While keeping the ratio of total demand to total supply constant, larger service areas tend to exacerbate the challenges faced by the VRPPD model due to the increased frequencies of long-distance cross-regional delivery tours. In contrast, the localization strategy enabled by transshipment mitigates the impact of regional size on meal delivery and maintains efficiency by dynamically adjusting the number of sub-areas to accommodate varying sizes of the service region. In summary, transshipment through a microhub offers significant advantages over non-transshipment strategies in scenarios with high demand or insufficient supply, particularly in servicing areas of larger sizes.

## 5.4 A case study of Meituan meal delivery

We conduct a case study comparing transshipment and non-transshipment models using Meituan meal delivery data. The original dataset comprises 654,343 order requests recorded between October 17 and 24, 2022, in an undisclosed city, publicly provided by Meituan. After filtering out requests located far from the central cluster, the master dataset retains 484,650 order records, 4,571 deliverers, and 206,748 trajectories. In Figure 9a, we define a hexagonal service area with a radius of R = 2.5 miles, intended for transshipment operations through a single microhub, and assume that orders are uniformly and randomly distributed within the hexagon. For Figure 9b, demand is defined as the number of orders occurring per hour within the service area, while supply is measured as the total active labor hours of deliverers derived from their trajectory data. Both demand and supply values are averaged over the week to represent a typical day.

For each pair of demand and supply data from 10:00 to 21:00, we independently solve the optimal design problem with transshipment (17) and without (18). The relative differences in two key metrics-customers' average waiting time and deliverers' vehicle miles traveled-are shown in Figure 10, where greener shades indicate that the transshipment model is more efficient than

Figure 9: Visualization of Meituan meal delivery data: (a) a map showing customer locations and the targeted area with a radius of R = 2.5 miles for analysis; (b) hourly variations in demand and supply within the targeted service area throughout the day.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.01047v2-figure-fa3bfb99-282f-460c-8f83-80d646b7df47.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.01047v2-figure-29388321-4094-49f4-bf0b-49f21fd93c8c.jpg)

Figure 10: Savings in deliverers' vehicle miles traveled and customers' average waiting time achieved by transshipment compared to non-transshipment.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.01047v2-figure-277bbdb9-9076-40ac-a1df-ff0a8506208f.jpg)

the non-transshipment counterpart. As anticipated from previous analysis, the transshipment strategy constantly reduces VMT across all time periods compared to the non-transshipment approach. For customer waiting time, during lunch and dinner peak hours, when demand significantly exceeds supply, the transshipment strategy achieves up to a 45% reduction, with the greatest saving occurring between 11:00 and 12:00, coinciding with the highest demand period of the day. However, during non-peak hours, the VRPPD model without transshipment outperforms the transshipment strategy in providing shorter customer waiting times. These findings suggest an alternating approach for meal delivery, where transshipment is adopted during peak hours, while the system reverts to the classic pickup-and-delivery model during non-peak hours.

In Appendix C, we reduce the service range covered by individual microhubs for transshipment and deploy multiple microhubs to cover the entire region. The results further confirm that transshipment is particularly advantageous in high-demand areas during peak hours.

## 6 Conclusion

This paper investigates a novel meal delivery strategy that leverages a microhub for transshipment. An analytical model was developed using a CA method to evaluate the system's expected performance under the transshipment strategy. To estimate customer waiting time, we decomposed the lifecycle of a meal order and adapted a queuing network to approximate the time spent by orders at different service stages. A comprehensive set of simulations was conducted to validate the accuracy of these approximations. The analytical approximations were then incorporated into an optimal design problem to solve for the performance envelope under various exogenous market conditions. The numerical results indicate that when the system is optimally configured, the proposed meal delivery strategy through transshipment significantly reduces vehicle miles traveled compared to the conventional VRPPD strategy without transshipment. Additionally, the transshipment strategy effectively decreases expected customer waiting time, particularly under high demand or low supply conditions, and proves especially advantageous for servicing larger areas. However, in scenarios with sufficient supply, excessive reliance on transshipment may lead to longer customer waiting times. These insights into the applicability of the transshipment strategy were finally validated using Meituan meal delivery data.

Future research can explore several promising directions. First, this study examines the simplest scenario of a single, centrally placed microhub for transshipment in an isotropic area. A natural extension would be to analyze heterogeneous settings where both demand and supply are unevenly distributed and determine the optimal location of the microhub to accommodate spatial variations in delivery patterns. Second, future research can also explore scenarios involving multiple microhubs optimally positioned to facilitate transshipment within a larger area. In such cases, a two-echelon delivery system could be investigated, where micromobility options such as e-bikes or e-scooters are dedicated to servicing smaller areas covered by individual microhubs, while heavy-duty delivery vehicles periodically transport packages between microhubs. Third, toward the end of the paper, we sketchily rationalize an alternating strategy for meal delivery to adopt transshipment during peak periods and switch back to non-transshipment strategies during non-peak periods. However, the current modeling framework treats different time periods independently, without factoring in transition costs, start-up delays, or relocation complexities. Future work could incorporate these dynamic transition effects to provide a more realistic evaluation of the dynamic system design. Lastly, it would be valuable to endogenize both customer demand and deliverer supply in the model. Customer demand might fluctuate in response to expected waiting time, while deliverer availability could be influenced by overall travel distances and order completion rates. Integrating these behavioral and economic factors would yield more accurate predictions of system performance and help platform operators design more robust and adaptive delivery strategies.

## Acknowledgement

This research utilized data from the TSL Data-Driven Research Challenge, provided by Meituan. We extend our heartfelt gratitude to Meituan for sharing their service data. This work is partially supported by the National Science Foundation (CMMI-2443338).

## Declaration of generative Al and Al-assisted technologies in the writing process

During the preparation of this work the authors used ChatGPT in order to enhance the readibility of the paper. After using this tool/service, the authors reviewed and edited the content as needed and took full responsibility for the content of the published article.

## Nomenclature

Table 1: Notation list of parameters and variables

| Notation           | Description                                                             | Unit          |
|--------------------|-------------------------------------------------------------------------|---------------|
| Input Parameters   | Input Parameters                                                        |               |
| λ                  | Expected number of orders generated per hour per unit area              | 1/hour · mi 2 |
| m                  | Total size of delivery fleet                                            |               |
| v                  | Speed of delivery vehicles                                              | mph           |
| A                  | Total area of the region                                                | mi 2          |
| a , b , α , β      | Parameters related to TSP approximations in a pie-shaped area           |               |
| π Q                | Operating cost per vehicle mile traveled                                | $             |
| π W                | Value of time                                                           | $/hour        |
| Decision Variables | Decision Variables                                                      |               |
| K                  | Number of sub-area partitions                                           |               |
| A k                | Area of each sub-area k                                                 | mi 2          |
| λ k                | Expected number of new orders generated per hour in sub-area k          | 1/hour        |
| δ k                | Number of pickup and drop-off locations to visit in sub-area k          | 1/hour        |
| p i , k            | Probability of an order picked up in sub-area i being dropped off in k  |               |
| m k                | Number of deliverers assigned to sub-area k                             |               |
| n k                | Batch size of packages per delivery tour in sub-area k                  |               |
| S k                | Traveling time of individual deliverer in sub-area k per delivery cycle | hour          |
| Q                  | Expected vehicle miles traveled by deliverers per hour of operation     | mph           |
| W a k              | Average waiting time for accumulating n k packages in sub-area k        | hour          |
| W q k              | Average package holding time for available deliverers in sub-area k     | hour          |
| W                  | Average total waiting time per customer                                 | hour          |

## Appendix A Proof of Propositions

This appendix presents the proofs for the propositions introduced in Sections 3 and 4.

## A.1 Proof of Proposition 1

Consider a pie-shaped circular sector with radius R and central angle θ , containing N nodes randomly distributed across the entire region. Let Di denote the distance from point i within the sector to the center of the circle. Then, the distance from the origin to the farthest node is defined as R ′ = max { Di } ∈ i N , and its CDF and PDF are given by the following expressions:

$$F _ { R ^ { \prime } } ( r ) = \mathbb { P } ( \max \{ D _ { i } \} _ { i \in N } \leq r ) = \mathbb { P } ( D _ { 1 } \leq r, \cdots, D _ { N } \leq r ) = \left ( \frac { \frac { \theta } { 2 \pi } \pi r ^ { 2 } } { \frac { \theta } { 2 \pi } \pi R ^ { 2 } } \right ) ^ { N } = \frac { r ^ { 2 N } } { R ^ { 2 N } },$$

$$f _ { R ^ { \prime } } ( r ) = 2 N \frac { r ^ { 2 N - 1 } } { R ^ { 2 N } }.$$

The expected value of R ′ is thus given by

$$\mathbb { E } [ R ^ { \prime } ] = \int _ { 0 } ^ { R } r \cdot f _ { R ^ { \prime } } ( r ) \, \mathrm d r = \frac { 2 N } { 2 N + 1 } R.$$

## A.2 Proof of Proposition 2

For the Spatial Poisson process, let δ = NO + ND π R 2 denote the density of access points across the entire area. The survival function and PDF of the distance X from a point to its nearest neighboring point are given by

$$& \mathbb { P } ( X > x ) = \mathbb { P } ( \text{no points within distance } x ) = e ^ { - \delta \cdot \pi x ^ { 2 } }, \\ & f _ { X } ( x ) = \delta \cdot 2 \pi x \cdot e ^ { - \delta \cdot \pi x ^ { 2 } }, \ x \geq 0.$$

Similarly, for the variable Y , which follows the Rayleigh distribution, the survival function and PDF are given by:

$$& \mathbb { P } ( Y > y ) = e ^ { - \frac { y ^ { 2 } } { 2 \sigma ^ { 2 } } }, \\ & f _ { Y } ( y ) = \frac { y } { \sigma ^ { 2 } } \cdot e ^ { - \frac { y ^ { 2 } } { 2 \sigma ^ { 2 } } }, \, y \geq 0.$$

Thus, the probability of X &gt; Y can be derived as:

$$\text{$line, the probability of $X > Y$ can be aerive as} \colon \\ \mathbb { P } ( X > Y ) & = \int _ { 0 } ^ { \infty } \int _ { y } ^ { \infty } f _ { X } ( x ) \ f _ { Y } ( y ) \, \mathrm d \, x \, \mathrm d \, y \\ & = \int _ { 0 } ^ { \infty } \int _ { y } ^ { \infty } 2 \delta \pi x e ^ { - \delta \pi x ^ { 2 } } \cdot \frac { y } { \sigma ^ { 2 } } e ^ { - \frac { y ^ { 2 } } { 2 \sigma ^ { 2 } } } \, \mathrm d \, x \, \mathrm d \, y \\ & = \int _ { 0 } ^ { \infty } \left ( \frac { y } { \sigma ^ { 2 } } e ^ { - \frac { y ^ { 2 } } { 2 \sigma ^ { 2 } } } \right ) \left ( \int _ { y } ^ { \infty } 2 \delta \pi x e ^ { - \delta \pi x ^ { 2 } } \, \mathrm d \, x \right ) \mathrm d \, y.$$

First, considering the integral over x , we have IP ( X &gt; y ) = ∫ ∞ y 2 δπ xe -δπ x 2 d x = e -δπ y 2 . Substituting it back yields I P ( X &gt; Y ) = ∫ ∞ 0 y σ 2 · e -( δπ + 1 2 σ 2 ) y 2 d y . By letting a = δπ + 1 2 σ 2 , u = ay 2 , and d u = 2 ay d y , we obtain

$$\mathfrak { 1 } \, u = 2 a y d \, y, \, \text{we obtain} \\ \mathbb { P } ( X > Y ) & = \frac { 1 } { \sigma ^ { 2 } } \cdot \frac { 1 } { 2 a } \int _ { 0 } ^ { \infty } e ^ { - u } \mathrm d \, u \\ & = \frac { 1 } { 1 + 2 \delta \pi \sigma ^ { 2 } } \\ & = \frac { 1 } { 1 + \frac { 2 ( N _ { O } + N _ { D } ) \sigma ^ { 2 } } { R ^ { 2 } } }.$$

## A.3 Proof of Proposition 3

For variable X , the survival function is given by:

$$\mathbb { P } ( X > z ) = \mathbb { P } ( \text{no points within distance} \, z ) = e ^ { - \delta \cdot \pi z ^ { 2 } } = e ^ { - \frac { N _ { O } + N _ { D } } { R ^ { 2 } } z ^ { 2 } }.$$

For variable Y , by substituting σ = 1 √ 2 µ , its survival function becomes:

$$\mathbb { P } ( Y > z ) = e ^ { - \frac { z ^ { 2 } } { 2 \sigma ^ { 2 } } } = e ^ { - \mu z ^ { 2 } }.$$

The expectation of min ( X Y , ) is then given by:

$$\text{The expectation of min} ( X, Y ) \text{ is then given by} \text{:} \\ \mathbb { E } [ \min ( X, Y ) ] & = \int _ { 0 } ^ { + \infty } \mathbb { I } P ( X > z, Y > z ) \, d z \\ & = \int _ { 0 } ^ { + \infty } \mathbb { I } P ( X > z ) \cdot \mathbb { I } P ( Y > z ) \, d z \\ & = \int _ { 0 } ^ { + \infty } e ^ { - \left ( \frac { N _ { O } + N _ { D } } { R ^ { 2 } } + \mu \right ) z ^ { 2 } } \, d z \\ & \text{Gaussian integral} \left [ \frac { 1 } { 2 } \sqrt { \frac { \pi } { N _ { O } + N _ { D } } { R ^ { 2 } } + \frac { 1 } { 2 \sigma ^ { 2 } } }.$$

## Appendix B Validation of Rayleigh Distribution Assumption for Order Distances

This appendix details the calibration and validation of using Rayleigh distribution to model the delivery distance d between the origin and destination of meal orders. The PDF of the Rayleigh distribution is given by:

$$f ( d | \sigma ) = \frac { d } { \sigma ^ { 2 } } e ^ { - \frac { d ^ { 2 } } { 2 \sigma ^ { 2 } } }, \ \ d > 0, \ \sigma > 0,$$

where σ is the scale parameter. Given the observed distances { d 1, d 2, · · · , dn } for n orders, the log-likelihood function is:

$$\mathcal { L } ( \sigma ) = \sum _ { i = 1 } ^ { n } \left ( \ln d _ { i } - 2 \ln \sigma - \frac { d _ { i } ^ { 2 } } { 2 \sigma ^ { 2 } } \right ).$$

To obtain the maximum likelihood estimate (MLE) of σ , we differentiate L ( σ ) with respect to σ , set the derivative to zero, and solve:

$$\frac { \partial \mathcal { L } ( \sigma ) } { \partial \sigma } = - \frac { 2 n } { \sigma } + \frac { \sum _ { i = 1 } ^ { n } d _ { i } ^ { 2 } } { \sigma ^ { 3 } } = 0.$$

From the Meituan meal delivery order data, the MLE of σ is calculated as:

$$\hat { \sigma } = \sqrt { \frac { \sum _ { i = 1 } ^ { n } d _ { i } ^ { 2 } } { 2 n } } = 0. 8 3.$$

Figure B1 compares the fitted Rayleigh distribution with the histogram of 423,449 observed delivery distances for meal orders recorded during the week of October 17-24, 2022, between 10:00 AMand 10:00 PM each day. The average delivery distance is approximately 1 mile, with a higher concentration of orders occurring below this distance. This skewed distribution is well captured by the Rayleigh model, which achieves a mean squared error of 0.001249 miles.

Figure B1: Fitting of the Rayleigh distribution to meal delivery order distances.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.01047v2-figure-8eef7cbd-7697-4e03-991f-fd1749ef5c36.jpg)

## Appendix C Case study with smaller transshipment range

In this appendix, we reduce the transshipment service range through a single microhub from R = 2.5 in Section 5.4 to R = 1.5, dividing the entire region is divided into 10 hexagons, each with distinct demand and supply patterns. Each microhub is used exclusively for transshipment within its respective hexagon. Figure B1 illustrates the zonal divisions, while Figure B2 presents the demand and supply for each hexagon at different times of the day. The two numbers in each cell indicate the number of orders and the labor hours available, respectively. Those hexagons with higher demand are highlighted in redder shades.

Figures B3 and B4 illustrate the relative savings achieved by comparing transshipment against non-transshipment in terms of average customer waiting time and vehicle miles traveled by deliverers across different areas and time periods. Among the ten divisions, Hexagons 2, 3, 5, and 6 are identified in Figure B2 as the busiest areas. In Figure B3, Hexagons 2 and 3 demonstrate improved customer waiting times when transitioning to transshipment during the lunch peak hour, which corresponds to the highest demand-to-supply ratio. However, in Hexagons 5 and 6, despite having the highest order volumes among all zones, the available supply is also sufficient, making the transshipment strategy less effective compared to the non-transshipping counterpart. Comparing the results with those in Section 5.4 suggests that transshipment becomes less efficient on smaller service ranges, which aligns with the findings in Section 5.3. While VMT savings are consistently observed in all scenarios when using transshipment, the benefits in customer waiting time remain highly sensitive to the balance between demand and supply.

Figure B1: A map showing customer locations and the service range covered by each microhub, with a radius of R = 1.5 miles.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.01047v2-figure-aa0dee99-323d-4cb9-ba58-a7cae1dc2c3e.jpg)

Figure B2: Hourly variations in demand and supply across the ten hexagons throughout the day.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.01047v2-figure-900bbb74-805c-4f23-b7aa-6b5bc23d72ba.jpg)

Figure B3: Savings in customer waiting time achieved through transshipment across the ten hexagonal areas throughout the day.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.01047v2-figure-414f5daa-0df8-42e7-9039-ebfbdcb7cbe7.jpg)

Figure B4: Savings in vehicle miles traveled by deliverers through transshipment across the ten hexagonal areas throughout the day.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.01047v2-figure-ce680876-dddf-4758-80eb-85c30fd826bf.jpg)

## References

| Agatz, N., Cho, S.-H., Sun, H., and Wang, H. (2024). Transportation-enabled services: Concept, framework, and research opportunities. Service Science , 16(1):1-21.                                                                                   |
|-------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------|
| Ansari, S., Basdere, ¸ M., Li, X., Ouyang, Y., and Smilowitz, K. (2018). Advancements in continu- ous approximation models for logistics and transportation systems: 1996-2016. Transportation Research Part B: Methodological , 107:229-252.         |
| Campbell, J. F. (2013). A continuous approximation model for time definite many-to-many trans- portation. Transportation Research Part B: Methodological , 54:100-112.                                                                                |
| Chen, P. W. and Nie, Y. M. (2017a). Analysis of an idealized system of demand adaptive paired- line hybrid transit. Transportation Research Part B: Methodological , 102:38-54.                                                                       |
| Chen, P. W. and Nie, Y. M. (2017b). Connecting e-hailing to mass transit platform: Analysis of relative spatial position. Transportation Research Part C: Emerging Technologies , 77:444-461.                                                         |
| Chen, Q., Lin, J., and Kawamura, K. (2012). Comparison of urban cooperative delivery and direct delivery strategies. Transportation research record , 2288(1):28-39.                                                                                  |
| Daganzo, C. (2005). Logistics systems analysis . Springer Science & Business Media.                                                                                                                                                                   |
| Daganzo, C. and Newell, G. (1986a). Design of multiple-vehicle delivery tours-i: a ring-radial network. Transportation Research Part B , pages 345-364.                                                                                               |
| Daganzo, C. F. (1978). An approximate analytic model of many-to-many demand responsive transportation systems. Transportation Research , 12(5):325-333.                                                                                               |
| Daganzo, C. F. (1984a). The distance traveled to visit n points with a maximum of c stops per vehicle: An analytic model and an application. Transportation science , 18(4):331-350.                                                                  |
| Daganzo, C. F. (1984b). The length of tours in zones of different shapes. Transportation Research Part B: Methodological , 18(2):135-145.                                                                                                             |
| Daganzo, C. F. (1987a). The break-bulk role of terminals in many-to-many logistic networks. Operations Research , 35(4):543-555.                                                                                                                      |
| Daganzo, C. F. (1987b). Increasing model precision can reduce accuracy. Transportation Science , 21(2):100-105.                                                                                                                                       |
| Daganzo, C. F., Gayah, V. V., and Gonzales, E. J. (2012). The potential of parsimonious models for understanding large scale transportation systems and answering big picture questions. EURO Journal on Transportation and Logistics , 1(1-2):47-65. |
| Daganzo, C. F. and Newell, G. F. (1986b). Configuration of physical distribution networks. Net- works , 16(2):113-132.                                                                                                                                |
| del Castillo, J. M. (1999). A heuristic for the traveling salesman problem based on a continuous approximation. Transportation Research Part B: Methodological , 33(2):123-152.                                                                       |
| Flushberg Sr, M. and Wilson, N. H. (1976). A descriptive supply model for demand-responsive transportation system planning.                                                                                                                           |

| Ghaffarinasab, N., Van Woensel, T., and Minner, S. (2018). Acontinuous approximation approach to the planar hub location-routing problem: Modeling and solution algorithms. Computers & Operations Research , 100:140-154.   |
|------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------|
| Gunes, S., Fried, T., and Goodchild, A. (2024). Seattle microhub delivery pilot: Evaluating emis- sion impacts and stakeholder engagement. Case Studies on Transport Policy , 15:101119.                                     |
| Kaczmarski, M. (2024). Which company is winning the restaurant food delivery war? https://secondmeasure.com/datapoints/food-delivery-services-grubhub-uber-eats- doordash-postmates/.                                        |
| Kawamura, K. and Lu, Y. (2007). Evaluation of delivery consolidation in us urban areas with logistics cost analysis. Transportation research record , 2008(1):34-42.                                                         |
| Kohar, A. and Jakhar, S. K. (2021). A capacitated multi pickup online food delivery problem with time windows: a branch-and-cut algorithm. Annals of Operations Research , pages 1-22.                                       |
| Larson, R. C. and Odoni, A. R. (1981). Urban operations research . Number Monograph.                                                                                                                                         |
| Lin, J., Chen, Q., and Kawamura, K. (2016). Sustainability si: logistics cost and environmental impact analyses of urban delivery consolidation strategies. Networks and Spatial Economics , 16:227-253.                     |
| Liu, S., He, L., and Max Shen, Z.-J. (2021). On-time last-mile delivery: Order assignment with travel-time predictors. Management Science , 67(7):4095-4119.                                                                 |
| Liu, Y. (2019). An optimization-driven dynamic vehicle routing algorithm for on-demand meal delivery using drones. Computers & Operations Research , 111:1-20.                                                               |
| Naseraldin, H. and Herer, Y. T. (2011). A location-inventory model with lateral transshipments. Naval Research Logistics (NRL) , 58(5):437-456.                                                                              |
| Newell, G. F. (1971). Dispatching policies for a transportation route. Transportation Science , 5(1):91- 105.                                                                                                                |
| Newell, G. F. and Daganzo, C. F. (1986). Design of multiple vehicle delivery tours-ii other metrics. Transportation Research Part B: Methodological , 20(5):365-376.                                                         |
| Nourbakhsh, S. M. and Ouyang, Y. (2012). A structured flexible transit system for low demand areas. Transportation Research Part B: Methodological , 46(1):204-216.                                                          |
| Ouyang, Y. and Yang, H. (2023). Measurement and mitigation of the 'wild goose chase' phe- nomenon in taxi services. Transportation Research Part B: Methodological , 167:217-234.                                            |
| Pahwa, A. and Jaller, M. (2022). A cost-based comparative analysis of different last-mile strate- gies for e-commerce delivery. Transportation Research Part E: Logistics and Transportation Review , 164:102783.            |
| Shi, L. and Xu, Z. (2024). Dine in or takeout? trends on restaurant service demand amid the covid-19 pandemic. Service Science .                                                                                             |
| Shi, L., Xu, Z., Lejeune, M., and Luo, Q. (2022). An integer l-shaped method for dynamic order fulfillment in autonomous last-mile delivery with demand uncertainty. arXiv preprint arXiv:2208.09067 .                       |

| Simoni, M. D. and Winkenbach, M. (2023). Crowdsourced on-demand food delivery: An or- der batching and assignment algorithm. Transportation Research Part C: Emerging Technologies , 149:104055.                                             |
|----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------|
| Smilowitz, K. R. and Daganzo, C. F. (2007). Continuum approximation techniques for the design of integrated package distribution systems. Networks: An International Journal , 50(3):183-196.                                                |
| Steever, Z., Karwan, M., and Murray, C. (2019). Dynamic courier routing for a food delivery service. Computers & Operations Research , 107:173-188.                                                                                          |
| Stein, D. M. (1978). An asymptotic, probabilistic analysis of a routing problem. Mathematics of Operations Research , 3(2):89-101.                                                                                                           |
| Toth, P. and Vigo, D. (2002). The vehicle routing problem . SIAM.                                                                                                                                                                            |
| Tsao, Y.-C., Mangotra, D., Lu, J.-C., and Dong, M. (2012). A continuous approximation approach for the integrated facility-inventory allocation problem. European journal of operational research , 222(2):216-228.                          |
| Wiles, P. G. and van Brunt, B. (2001). Optimal location of transshipment depots. Transportation Research Part A: Policy and Practice , 35(8):745-771.                                                                                        |
| Wilson, N., Sussman, J., Goodman, L., and Higonnet, B. (1969). Simulation of a computer aided routing system (cars). Technical report, Institute of Electrical and Electronics Engineers (IEEE).                                             |
| Wilson, N. H. M., Weissberg, R. W., and Hauser, J. (1976). Advanced dial-a-ride algorithms research project. Technical report.                                                                                                               |
| Xie, W. and Ouyang, Y. (2015). Optimal layout of transshipment facility locations on an infinite homogeneous plane. Transportation Research Part B: Methodological , 75:74-88.                                                               |
| Xu, Z., Yin, Y., and Ye, J. (2020). On the supply curve of ride-hailing systems. Transportation Research Part B: Methodological , 132:29-43.                                                                                                 |
| Ye, A., Zhang, K., Chen, X. M., Bell, M. G., Lee, D.-H., and Hu, S. (2024). Modeling and manag- ing an on-demand meal delivery system with order bundling. Transportation Research Part E: Logistics and Transportation Review , 187:103597. |
Linxuan Shi Zhengtian Xu
Analysis of Transshipment in Three-Sided Meal Delivery Services via Microhubs
This paper introduces and analyzes a novel transshipment strategy for meal
delivery. In this approach, the service area is partitioned into smaller
sub-areas, with deliverers assigned to operate exclusively within these
sub-areas. Meanwhile, a centrally located microhub functions as a logistic
depot to facilitate the batching and transfer of meal packages toward different
sub-areas. We model the meal delivery system with transshipment using networked
G/G/m queues and analytically approximate two critical system performance
metrics -- customer waiting time and vehicle miles traveled -- to evaluate the
effectiveness of the proposed strategy. The performance achieved by
transshipment is benchmarked against that of the classic pickup-and-delivery
strategy without transshipment, both predicted using continuous approximations.
For the latter, we enhance the existing modeling by incorporating the delivery
distance profiles of individual orders to better match the meal delivery
context. Our comparisons indicate that meal delivery via transshipment
outperforms the non-transshipping counterpart across both metrics under either
high-demand or low-supply conditions, with particular advantages in servicing
larger areas or handling long-distance orders. This conclusion is corroborated
by a numerical experiment using empirical meal delivery data from Meituan,
which suggests that an optimally configured transshipment strategy can
significantly improve service performance for both customers and deliverers
during peak lunch hours and in the busiest districts. While transshipment
continues to reduce vehicle miles traveled by deliverers during non-peak hours,
it results in longer customer waiting times compared to the benchmark without
transshipment as demand decreases.