(Translated by https://www.hiragana.jp/)
On Rider Strategic Behavior in Ride-Sharing Platforms

On Rider Strategic Behavior in Ride-Sharing Platforms

Jay Mulay111Georgia Institute of Technology. Email: jay.mulay@gatech.edu    Diptangshu Sen222Georgia Institute of Technology. Email: dsen30@gatech.edu    Juba Ziani333Georgia Institute of Technology. Email: jziani3@gatech.edu
(August 8, 2024)
Abstract

Over the past decade, ride-sharing services have become increasingly important, with U.S. market leaders such as Uber and Lyft expanding to over 900 cities worldwide and facilitating billions of rides annually. This rise reflects their ability to meet users’ convenience, efficiency, and affordability needs. However, in busy areas and surge zones, the benefits of these platforms can diminish, prompting riders to relocate to cheaper, more convenient locations or seek alternative transportation.

While much research has focused on the strategic behavior of drivers, the strategic actions of riders, especially when it comes to riders walking outside of surge zones, remain under-explored. This paper examines the impact of rider-side strategic behavior on surge dynamics. We investigate how riders’ actions influence market dynamics, including supply, demand, and pricing. We show significant impacts, such as spillover effects where demand increases in areas adjacent to surge zones and prices surge in nearby areas. Our theoretical insights and experimental results highlight that rider strategic behavior helps redistribute demand, reduce surge prices, and clear demand in a more balanced way across zones.

1 Introduction

Over the past decade, ride-sharing services have grown exponentially in importance. U.S. market leaders such as Uber and Lyft have now expanded their operations to over 900 cities worldwide, facilitating billions of rides annually. The rise of ride-sharing platforms reflects how they have adeptly addressed users’ needs for convenience, efficiency, and reliability. One of the major reasons for the popularity of ride-sharing platforms is their ability to provide users with rides from virtually anywhere, at any time, at a relatively affordable price. However, the benefits of ride-sharing platforms can diminish in busy areas and surge zones, incentivizing strategic behavior by riders who may relocate to cheaper and more convenient locations or seek alternate means of transportation entirely.

Much academic work on ride-sharing has recognized the prevalence of strategic behavior among drivers within these platforms. Drivers can strategically decide which locations to serve or when to be online to maximize their earnings during surges; they can turn off the application while transiting to avoid taking rides on their way to more profitable locations; they can selectively accept rides based on profit and convenience, and even cancel rides when needed. While there is a significant focus on driver behavior in the literature, related work has largely neglected rider-side behavior. In particular, riders often strategically relocate to avoid surge zones and to secure cheaper rides during major events such as sports games or concerts.

The goal of this paper is to study the impact of rider-side strategic behavior on surge dynamics. We aim to understand how such rider behavior influences market dynamics, particularly the evolution of demand inside and outside the surge zone. This paper highlights several significant impacts: spillover effects where demand becomes higher than normal and price surges in areas adjacent to surge zones. We provide precise theoretical insights and experimental results to elucidate and quantify these effects. Our analysis reveals that rider strategic behavior helps redistribute and smooth out the demand across different areas, thereby reducing the price gap across the surge boundary.

Summary of Contributions: The main contributions of our work are as follows:

  • We introduce a new model highlighting the phenomenon of strategic riders on ride-sharing platforms in Section 2. In our model, riders can be waiting in one of two zones: a surge and a non-surge zone. Riders faced with high prices during a surge can choose to walk out of the surge zone in search of cheaper rides. We capture, in our initial model, the rate of riders moving across the surge boundary through a simple functional form that depends on the difference in demands between the surge and the non-surge areas. We extend this model to a more realistic agent-based and game-theoretic model in Section 4.3.1.

  • In Section 3, we analyze the dynamics of our initial theoretical model. We first characterize the shape and convergence properties of demand curves in both surge and non-surge zones (Theorems 3.1 and 3.2). We also show the existence of two types of surges: i) a localized type of surge where the surge demand does not accumulate outside the surge boundary (it remains localized as the name suggests) and ii) spill-over type of surge where the surge demand spills over and accumulates in the non-surge zone causing the demand there to increase.

  • Finally, in Section 4, we develop a much more comprehensive game-theoretic, agent-based simulation model of a ride-sharing platform where each rider is modeled as a rational agent. Agents decide to move to the non-surge zone if they get an increased utility from doing so, where this utility depends on prices in each zone and the cost of moving. This model also incorporates additional considerations to take the model one step closer to reality, including strategic riders that may prefer going to the more expensive surge zone, and stochasticity in arrival rates.

    Our agent-based model serves two purposes. First, it shows that our simple initial, theoretical model that solely aimed to capture the impact of strategic rider behavior still provides a good approximation to the more complex and realistic agent-based model, highlighting its value as a simple and tractable alternative. Second, it highlights the role of strategic riders in surge dynamics, in particular quantifying the extent to which this strategic behavior helps smooth out the demand and prices across zones.

1.1 Related Work

There is a rich body of literature on challenges faced by ride-sharing platforms in different settings: Agatz et al. [2012]; Furuhata et al. [2013]; Kooti et al. [2017]; Chan and Shaheen [2012]; Jin et al. [2018] to only name a few.

A prominent line of work involves studying pricing mechanisms for ridesharing platforms and in particular, i) static versus dynamic pricing models Banerjee et al. [2015]; Pandit et al. [2019], and ii) spatial, temporal, and spatio-temporal pricing strategies Bimpikis et al. [2019]; Ma et al. [2022]. Differential pricing strategies have also been explored, with an example being the work of Zhao et al. [2022]. Apart from pricing, there is also a significant amount of work that has studied the general operations management of ride-hailing systems, including empty vehicle rerouting, dispatch, matching and joint optimization  Braverman et al. [2019]; Chen and Cassandras [2020]; Alonso-Mora et al. [2017]; Yan et al. [2020]. Recently, researchers have also considered how ride-sharing platforms can encourage riders in congested areas to walk to nearby locations for pickup as it has been shown to improve the service efficiency by reducing pickup times for drivers  Yao et al. [2021]; Fielbaum et al. [2021].

Our work is closely related to the two lines of research. The first one studies surge pricing, when rider demand vastly outdoes the supply of drivers. Surge pricing aims to better match supply with demand during periods of high demand and lower prices when demand decreases; for example, Castillo et al. [2017] show how low supply relative to demand can lead to market failures and “wild goose chase” phenomena where free drivers are thinly spread and sent to far away locations, hurting both riders’ waiting times and drivers profits. They highlight the role of surge pricing as a solution to this problem. There have also been other works like  Castillo [2023] which study the overall welfare effects of surge pricing. Second, a recent line of work recognizes user incentives and strategic behavior in ride-sharing platforms, particularly driver-side. Drivers employ varied strategies to maximize profits and utility, including choosing which rides to accept, chasing surges, and strategically signing in and out to manipulate supply and prices, to only name a few. There has been increased focus on modeling and understanding this strategic driver behavior, such as Ashkrof et al. [2020]. Garg and Nazerzadeh [2022] and Cashore et al. [2022] study the design of incentive-compatible pricing mechanisms for drivers, while Rheingans-Yoo et al. [2019] study pricing that accounts for drivers’ location preferences. Tripathy et al. [2022]; Yu et al. [2022] focus on how coalitions of drivers can collude to manipulate the system. There have also been attempts to understand how such strategic behavior may affect the marketplace. For example,  Besbes et al. [2021] investigates the supply side’s spatial response to surge pricing when drivers are strategic. They show that optimal surge pricing induces a spatial decomposition of the supply - in some locations supply is exactly matched with demand, while in other locations supply may either be deliberately over-congested or completely diverted away. Unlike the driver side, there has been limited effort to account for riders being strategic.  Chen and Hu [2020] is an exception in this regard, it studies pricing policies for ride-sharing platforms when both drivers and riders are strategic and forward-looking, particularly riders can monitor prices and decide when to accept a ride.

The work which is perhaps most closely related to ours is the recent contribution of  Hu et al. [2022]. Similar to our work combining surge dynamics and rider strategic behavior, Hu et al. [2022] examine two-sided temporal dynamics and show that drivers respond more slowly to surge pricing compared to riders. They propose two types of equilibrium surge pricing policies and characterize the responses of both drivers and riders. In their paper, the main strategic aspect that is modelled for riders is that they can choose to wait for better prices. Our work differs significantly from theirs in the sense that we consider a different form of strategic behavior on the rider-side where they can choose to walk away from the surge zone in search of cheaper rides. We also characterize the effect of this strategic behavior on the surge dynamics itself - something the previous paper does not consider.

2 Model

Refer to caption
Figure 1: The above map schematic represents a geography serviced by a ride-sharing platform. The entire region is divided into hexagonal blocks and the price in each block depends on the market conditions “locally”. Ride-sharing platforms like Uber employ surge pricing algorithms at a “hyperlocal” level Uber [2024], and Uber in particular uses such hexagonal regions for pricing  Uber [2018]. The red block represents a region with abnormally high demand, leading to a surge in prices. We call this zone the “surge zone”. In all other surrounding blocks (marked in green), demand reflects usual levels; we call them the “non-surge zone”.

We study a model of rider dynamics at surge time on ride-sharing platforms where the riders can be strategic and are willing to walk outside of the surge zone to call a ride. This is motivated by behavior observed in practice at large events such as concerts and sports games and backed by empirical research  Selyukh [2015]. The paper takes the point of view of the relatively short time horizon that it takes for a surge to happen and clear (minutes to hours). The main goal of the paper is then to understand ride-sharing dynamics over this type of short time scales, with a focus on understanding how demand clears in surge areas during major events.

To model these dynamics, we study a geography consisting of two main regions: a central region where a major event is happening and where the demand for rides is unusually high, called the surge zone; and the surrounding region, initially facing typical demand levels, called the non-surge zone. The behavior of drivers (supply side) and riders (demand side) on the platform is formalized below.

Driver’s arrivals and behavior:

We consider non-strategic drivers that arrive in each zone at an equal, constant, deterministic rate of μ𝜇\muitalic_μ. We initially focus on non-strategic drivers to isolate and highlight the role of rider-side strategic behavior, which is the main focus of this paper. However, we note that the assumption that rates are equal across zones and that drivers are not strategic, used in the theoretical analysis of Section 3, is later relaxed in the experiments of Section 4.

Rider’s arrivals and behavior:

At time 00, we consider there is an initial demand in both the surge and the non-surge zone. We denote by Ds(t)subscript𝐷𝑠𝑡D_{s}(t)italic_D start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT ( italic_t ) and Dns(t)subscript𝐷𝑛𝑠𝑡D_{ns}(t)italic_D start_POSTSUBSCRIPT italic_n italic_s end_POSTSUBSCRIPT ( italic_t ) the unmet demands at the end of time step t𝑡titalic_t in the surge and non-surge zones respectively, and assume throughout the paper that Ds(0)Dns(0)much-greater-thansubscript𝐷𝑠0subscript𝐷𝑛𝑠0D_{s}(0)\gg D_{ns}(0)italic_D start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT ( 0 ) ≫ italic_D start_POSTSUBSCRIPT italic_n italic_s end_POSTSUBSCRIPT ( 0 ). Starting in the first time step, new riders arrive in each zone at a constant, deterministic rate of λ𝜆\lambdaitalic_λ. We assume444Otherwise, demand can never be cleared and builds up over time, leading to a trivial/uninteresting case. that λ<μ𝜆𝜇\lambda<\muitalic_λ < italic_μ.

Riders can leave each zone via two mechanisms: i) a rider gets matched to a driver and gets a ride; or, ii) a rider can be strategic  Alban [2024] (a departure from much of the related work in ride-sharing), and decide to leave the surge zone in favor of the non-surge zone, for example when the demand in the surge zone is significantly higher than the demand in the non-surge zone. We make the following assumption on rider behavior in case ii):

Assumption 2.1.

The fraction of the unmet rider demand Ds(t)subscript𝐷𝑠𝑡D_{s}(t)italic_D start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT ( italic_t ) that move from the surge to the non-surge zone at every time step is given by f(Ds(t)Dns(t))𝑓subscript𝐷𝑠𝑡subscript𝐷𝑛𝑠𝑡f\left(D_{s}(t)-D_{ns}(t)\right)italic_f ( italic_D start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT ( italic_t ) - italic_D start_POSTSUBSCRIPT italic_n italic_s end_POSTSUBSCRIPT ( italic_t ) ), for a known function f𝑓fitalic_f. Further, f(.)f(.)italic_f ( . ) is non-decreasing, in [0,1]01[0,1][ 0 , 1 ], and is positive only if Ds(t)>Dns(t)subscript𝐷𝑠𝑡subscript𝐷𝑛𝑠𝑡D_{s}(t)>D_{ns}(t)italic_D start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT ( italic_t ) > italic_D start_POSTSUBSCRIPT italic_n italic_s end_POSTSUBSCRIPT ( italic_t ).

The above condition abstracts away riders’ strategic behavior by providing a functional form on the rate at which riders decide to move strategically from the surge to the non-surge zone. First, we assume that f𝑓fitalic_f is non-decreasing; this is a natural assumption, as the higher the difference in demand and how crowded the surge and non-surge zones are, the more riders are willing to move. Second, we note that the rate is always non-negative: this models the fact that riders would sensibly walk from high-density areas like a stadium or a concert venue to lower-density areas, but not the other way around (e.g., a rider in a busy downtown area walking closer to a major venue or stadium experiencing a surge) may be unrealistic. Second, we note that our functional form makes a linearity assumption, where the rate of movement is a function of the difference between demands in the surge and non-surge zones. We provide a motivating explanation of this assumption below:

Example 2.1.

Assume a simple rider behavior where they make strategic decisions based on price multipliers set by the platform. Denote ms(t)subscript𝑚𝑠𝑡m_{s}(t)italic_m start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT ( italic_t ) and mns(t)subscript𝑚𝑛𝑠𝑡m_{ns}(t)italic_m start_POSTSUBSCRIPT italic_n italic_s end_POSTSUBSCRIPT ( italic_t ) respectively the price multipliers in the surge and non-surge zones. Suppose that each agents estimates the price they pay in each zone according to the multipliers, and have a cost c𝑐citalic_c for moving from the surge to non-surge zone. Then the rider’s utility for not moving from the surge zone is ms(t)subscript𝑚𝑠𝑡m_{s}(t)italic_m start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT ( italic_t ), while the utility for moving is mns(t)+csubscript𝑚𝑛𝑠𝑡𝑐m_{ns}(t)+citalic_m start_POSTSUBSCRIPT italic_n italic_s end_POSTSUBSCRIPT ( italic_t ) + italic_c. In this case, a rider moves if and only if cms(t)mns(t)𝑐subscript𝑚𝑠𝑡subscript𝑚𝑛𝑠𝑡c\leq m_{s}(t)-m_{ns}(t)italic_c ≤ italic_m start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT ( italic_t ) - italic_m start_POSTSUBSCRIPT italic_n italic_s end_POSTSUBSCRIPT ( italic_t ). When the costs c𝑐citalic_c are i.i.d. from a distribution with c.d.f. F(.), the rate at which riders move is F(ms(t)mns(t))𝐹subscript𝑚𝑠𝑡subscript𝑚𝑛𝑠𝑡F(m_{s}(t)-m_{ns}(t))italic_F ( italic_m start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT ( italic_t ) - italic_m start_POSTSUBSCRIPT italic_n italic_s end_POSTSUBSCRIPT ( italic_t ) ). Under a simple pricing rule that is proportional to demand, we recover Assumption 2.1: one such example of a pricing rule could be where the price multiplier depends on the supply-demand mismatch, i.e., ms(t)Ds(t)μproportional-tosubscript𝑚𝑠𝑡subscript𝐷𝑠𝑡𝜇m_{s}(t)\propto\frac{D_{s}(t)}{\mu}italic_m start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT ( italic_t ) ∝ divide start_ARG italic_D start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT ( italic_t ) end_ARG start_ARG italic_μ end_ARG and mns(t)Dns(t)μproportional-tosubscript𝑚𝑛𝑠𝑡subscript𝐷𝑛𝑠𝑡𝜇m_{ns}(t)\propto\frac{D_{ns}(t)}{\mu}italic_m start_POSTSUBSCRIPT italic_n italic_s end_POSTSUBSCRIPT ( italic_t ) ∝ divide start_ARG italic_D start_POSTSUBSCRIPT italic_n italic_s end_POSTSUBSCRIPT ( italic_t ) end_ARG start_ARG italic_μ end_ARG.

Another motivating explanation is that riders make decision based on perceived demand around them. A rider that is in a crowded area perceives that they may be a large demand differential between the surge zone and adjacent, outer areas, and may want to switch locations based on that perceived differential.

We note that the assumptions that i) the rate of rider arrivals are deterministic and ii) Assumption 2.1 on the specific functional form of riders’ behavior, made for the theoretical results of Section 3, are relaxed in the experiments of Section 4 in favor of stochastic arrival rates and a game-theoretic model of riders’ decisions on which zone to serve as a function of platform prices.

3 Surge: Dynamics and Convergence

This section presents a theoretical perspective to the dynamics of surge when riders are strategic in nature. The section is organized as follows: we start by deriving the equations of dynamics from first principles in Section 3.1, followed by an analysis of the shapes of surge and non-surge demands and providing bounds on the times to convergence in Section 3.2.

3.1 Equations of Dynamics

In this segment, we explore the dynamics of the unmet demands inside and outside the surge boundary, given by Ds(t)subscript𝐷𝑠𝑡D_{s}(t)italic_D start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT ( italic_t ) and Dns(t)subscript𝐷𝑛𝑠𝑡D_{ns}(t)italic_D start_POSTSUBSCRIPT italic_n italic_s end_POSTSUBSCRIPT ( italic_t ) respectively. We start by formally deriving the one-step equation for the evolution of Ds(t)subscript𝐷𝑠𝑡D_{s}(t)italic_D start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT ( italic_t ) and Dns(t)subscript𝐷𝑛𝑠𝑡D_{ns}(t)italic_D start_POSTSUBSCRIPT italic_n italic_s end_POSTSUBSCRIPT ( italic_t ) from first principles.

Remember that Ds(t)subscript𝐷𝑠𝑡D_{s}(t)italic_D start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT ( italic_t ) is the unmet demand in the surge zone at the end of time period t𝑡titalic_t. We express the unmet demand at the end of period t+1𝑡1t+1italic_t + 1 recursively as a function of the unmet demands at the end of the previous time period t𝑡titalic_t. We do this by accounting for all inflows and outflows of demand from the surge zone.

Inflows:

Since riders can only leave the surge zone by being matched with a ride or choosing to relocate to the non-surge zone at the start of time t𝑡titalic_t, a carryover demand of size Ds(t)subscript𝐷𝑠𝑡D_{s}(t)italic_D start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT ( italic_t ) is already in the surge zone at the start of period t+1𝑡1t+1italic_t + 1. The surge zone also has an exogenous arrival of riders at rate λ𝜆\lambdaitalic_λ which means that λ𝜆\lambdaitalic_λ riders arrive over the course of period t+1𝑡1t+1italic_t + 1. Therefore, the total demand to be cleared in period t+1𝑡1t+1italic_t + 1 is given by Ds(t)+λsubscript𝐷𝑠𝑡𝜆D_{s}(t)+\lambdaitalic_D start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT ( italic_t ) + italic_λ.

Outflows:

Since Ds(t)subscript𝐷𝑠𝑡D_{s}(t)italic_D start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT ( italic_t ) number of riders had already been waiting to be matched since period t𝑡titalic_t, a fraction of them may decide to relocate to the non-surge zone at the start of period t+1𝑡1t+1italic_t + 1. As described earlier, this fraction is captured by f()𝑓f(\cdot)italic_f ( ⋅ ) which is a function of the demand differential Ds(t)Dns(t)subscript𝐷𝑠𝑡subscript𝐷𝑛𝑠𝑡D_{s}(t)-D_{ns}(t)italic_D start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT ( italic_t ) - italic_D start_POSTSUBSCRIPT italic_n italic_s end_POSTSUBSCRIPT ( italic_t ) across the surge boundary. There is also a supply of new drivers at fixed rate μ𝜇\muitalic_μ in each time period, so μ𝜇\muitalic_μ rides arrive over the course of period t+1𝑡1t+1italic_t + 1, get matched with rides and leave the system. Therefore, the total outflow of demand from the surge zone during period t+1𝑡1t+1italic_t + 1 is at most μ+f(Ds(t)Dns(t))Ds(t)𝜇𝑓subscript𝐷𝑠𝑡subscript𝐷𝑛𝑠𝑡subscript𝐷𝑠𝑡\mu+f(D_{s}(t)-D_{ns}(t))\cdot D_{s}(t)italic_μ + italic_f ( italic_D start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT ( italic_t ) - italic_D start_POSTSUBSCRIPT italic_n italic_s end_POSTSUBSCRIPT ( italic_t ) ) ⋅ italic_D start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT ( italic_t ).

Equations of dynamics:

Hence, we can express Ds(t+1)subscript𝐷𝑠𝑡1D_{s}(t+1)italic_D start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT ( italic_t + 1 ) in the following recursive form:

Ds(t+1)=max[0,Ds(t)+(λμ)f(Ds(t)Dns(t))Ds(t)].subscript𝐷𝑠𝑡10subscript𝐷𝑠𝑡𝜆𝜇𝑓subscript𝐷𝑠𝑡subscript𝐷𝑛𝑠𝑡subscript𝐷𝑠𝑡\displaystyle D_{s}(t+1)=\max\left[0,D_{s}(t)+(\lambda-\mu)-f\left(D_{s}(t)-D_% {ns}(t)\right)\cdot D_{s}(t)\right].italic_D start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT ( italic_t + 1 ) = roman_max [ 0 , italic_D start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT ( italic_t ) + ( italic_λ - italic_μ ) - italic_f ( italic_D start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT ( italic_t ) - italic_D start_POSTSUBSCRIPT italic_n italic_s end_POSTSUBSCRIPT ( italic_t ) ) ⋅ italic_D start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT ( italic_t ) ] . (1)

Note that the max(0,)0\max\left(0,\cdot\right)roman_max ( 0 , ⋅ ) ensures that outflow of demand cannot exceed the total demand to be cleared. Similarly, for Dns(t)subscript𝐷𝑛𝑠𝑡D_{ns}(t)italic_D start_POSTSUBSCRIPT italic_n italic_s end_POSTSUBSCRIPT ( italic_t ), we obtain:

Dns(t+1)=max[0,Dns(t)+(λμ)+f(Ds(t)Dns(t))Ds(t)].subscript𝐷𝑛𝑠𝑡10subscript𝐷𝑛𝑠𝑡𝜆𝜇𝑓subscript𝐷𝑠𝑡subscript𝐷𝑛𝑠𝑡subscript𝐷𝑠𝑡\displaystyle D_{ns}(t+1)=\max\left[0,D_{ns}(t)+(\lambda-\mu)+f\left(D_{s}(t)-% D_{ns}(t)\right)\cdot D_{s}(t)\right].italic_D start_POSTSUBSCRIPT italic_n italic_s end_POSTSUBSCRIPT ( italic_t + 1 ) = roman_max [ 0 , italic_D start_POSTSUBSCRIPT italic_n italic_s end_POSTSUBSCRIPT ( italic_t ) + ( italic_λ - italic_μ ) + italic_f ( italic_D start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT ( italic_t ) - italic_D start_POSTSUBSCRIPT italic_n italic_s end_POSTSUBSCRIPT ( italic_t ) ) ⋅ italic_D start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT ( italic_t ) ] . (2)
Boundary Conditions:

In order to solve for Ds(t)subscript𝐷𝑠𝑡D_{s}(t)italic_D start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT ( italic_t ) and Dns(t)subscript𝐷𝑛𝑠𝑡D_{ns}(t)italic_D start_POSTSUBSCRIPT italic_n italic_s end_POSTSUBSCRIPT ( italic_t ), we also need to specify the state of the system at t=0𝑡0t=0italic_t = 0 which constitute our boundary conditions. We set:

Ds(0)=D0;andDns(0)=d0;formulae-sequencesubscript𝐷𝑠0subscript𝐷0andsubscript𝐷𝑛𝑠0subscript𝑑0D_{s}(0)=D_{0};\quad\text{and}\quad D_{ns}(0)=d_{0};italic_D start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT ( 0 ) = italic_D start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ; and italic_D start_POSTSUBSCRIPT italic_n italic_s end_POSTSUBSCRIPT ( 0 ) = italic_d start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ;

D0subscript𝐷0D_{0}italic_D start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT is the large initial demand that arises in the surge zone at t=0𝑡0t=0italic_t = 0 (e.g., people leaving a stadium at the end of a football game), while d0subscript𝑑0d_{0}italic_d start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT is the uncleared demand in the non-surge zone at the same time. We assume that D0d0much-greater-thansubscript𝐷0subscript𝑑0D_{0}\gg d_{0}italic_D start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ≫ italic_d start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT.

The rest of the section studies the evolution of the system dynamics from Equations (1) and (2) above.

3.2 Properties of Surge and Non-surge Demands

Our main contribution in this segment is to characterize the shapes of the surge and non-surge demands over time in Theorems 3.1 and 3.2. Recall that D0subscript𝐷0D_{0}italic_D start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT and d0subscript𝑑0d_{0}italic_d start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT are the initial unmet demands in the surge and non-surge zones, λ𝜆\lambdaitalic_λ is the rate of exogenous demand and μ𝜇\muitalic_μ is the rate of supply in each zone. We start by characterizing properties of Ds(t)subscript𝐷𝑠𝑡D_{s}(t)italic_D start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT ( italic_t ), the unmet demand in the surge area.

Theorem 3.1.

Let Ds(t)subscript𝐷𝑠𝑡D_{s}(t)italic_D start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT ( italic_t ) be the unmet demand inside the surge zone at any time step t{0,1,2,}𝑡012t\in\{0,1,2,...\}italic_t ∈ { 0 , 1 , 2 , … }. Then, Ds(t)subscript𝐷𝑠𝑡D_{s}(t)italic_D start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT ( italic_t ) is always non-increasing in t𝑡titalic_t and must converge to 00 in a finite number of time steps. Further, if τssubscript𝜏𝑠\tau_{s}italic_τ start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT indicates the number of time steps till convergence, then τssubscript𝜏𝑠\tau_{s}italic_τ start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT must satisfy:

D0(μλ)+D0f(D0)τsD0(μλ).subscript𝐷0𝜇𝜆subscript𝐷0𝑓subscript𝐷0subscript𝜏𝑠subscript𝐷0𝜇𝜆\left\lfloor\frac{D_{0}}{(\mu-\lambda)+D_{0}\cdot f(D_{0})}\right\rfloor\leq% \tau_{s}\leq\left\lceil\frac{D_{0}}{(\mu-\lambda)}\right\rceil.⌊ divide start_ARG italic_D start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_ARG start_ARG ( italic_μ - italic_λ ) + italic_D start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ⋅ italic_f ( italic_D start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) end_ARG ⌋ ≤ italic_τ start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT ≤ ⌈ divide start_ARG italic_D start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_ARG start_ARG ( italic_μ - italic_λ ) end_ARG ⌉ .
Proof.

The detailed proof of all parts of the theorem can be found in Appendix A. ∎

Theorem 3.1 provides us three key insights about Ds(t)subscript𝐷𝑠𝑡D_{s}(t)italic_D start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT ( italic_t ): i) Ds(t)subscript𝐷𝑠𝑡D_{s}(t)italic_D start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT ( italic_t ) is always non-increasing; ii) Ds(t)subscript𝐷𝑠𝑡D_{s}(t)italic_D start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT ( italic_t ) converges to 00 in finite time; and finally, iii) there are interpretable upper and lower bounds on the number of time steps to convergence. Recall that with μ>λ𝜇𝜆\mu>\lambdaitalic_μ > italic_λ, (μλ)𝜇𝜆(\mu-\lambda)( italic_μ - italic_λ ) represents the excess supply in the zone at each time step after the exogenous demand for that period has been matched and cleared. Therefore, if all riders happened to be non-strategic and chose not to walk, the excess supply in each time step would be used to clear the pending demand from t=0𝑡0t=0italic_t = 0 and it would take exactly D0μλsubscript𝐷0𝜇𝜆\left\lceil\frac{D_{0}}{\mu-\lambda}\right\rceil⌈ divide start_ARG italic_D start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_ARG start_ARG italic_μ - italic_λ end_ARG ⌉ time steps for all the demand to dissipate. This shows that the upper bound is tight. In reality, however, the surge demand will dissipate faster than this because of the effect of unmatched customers from the surge zone gradually relocating to non-surge zones over time. On the other hand, the lower bound is obtained by noting that (μλ)+D0f(D0)𝜇𝜆subscript𝐷0𝑓subscript𝐷0(\mu-\lambda)+D_{0}\cdot f(D_{0})( italic_μ - italic_λ ) + italic_D start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ⋅ italic_f ( italic_D start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) represents the fastest possible rate of clearing demand (D0f(D0)subscript𝐷0𝑓subscript𝐷0D_{0}\cdot f(D_{0})italic_D start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ⋅ italic_f ( italic_D start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) represents the maximum outflow of strategic riders in any time period which can be achieved at t=0𝑡0t=0italic_t = 0 from the surge to the non-surge zone, and μλ𝜇𝜆\mu-\lambdaitalic_μ - italic_λ the rate that is cleared by matching demand to supply). For t>0𝑡0t>0italic_t > 0, the rate of clearing demand can only be smaller because of the smaller value of Ds(t)subscript𝐷𝑠𝑡D_{s}(t)italic_D start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT ( italic_t ) and the smaller demand differential across the surge boundary, leading to a smaller total amount of agents moving. Therefore, the surge demand cannot dissipate earlier than D0(μλ)+D0f(D0)subscript𝐷0𝜇𝜆subscript𝐷0𝑓subscript𝐷0\left\lfloor\frac{D_{0}}{(\mu-\lambda)+D_{0}\cdot f(D_{0})}\right\rfloor⌊ divide start_ARG italic_D start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_ARG start_ARG ( italic_μ - italic_λ ) + italic_D start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ⋅ italic_f ( italic_D start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) end_ARG ⌋ time steps. This bound is exactly tight in the corner case where the rates of clearing are so high that (μλ)+D0f(D0)D0𝜇𝜆subscript𝐷0𝑓subscript𝐷0subscript𝐷0(\mu-\lambda)+D_{0}\cdot f(D_{0})\geq D_{0}( italic_μ - italic_λ ) + italic_D start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ⋅ italic_f ( italic_D start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) ≥ italic_D start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT and the entire demand is cleared in time step 1111.

Now, we investigate the behavior of Dns(t)subscript𝐷𝑛𝑠𝑡D_{ns}(t)italic_D start_POSTSUBSCRIPT italic_n italic_s end_POSTSUBSCRIPT ( italic_t ). We present a complete characterization of the shape properties of Dns(t)subscript𝐷𝑛𝑠𝑡D_{ns}(t)italic_D start_POSTSUBSCRIPT italic_n italic_s end_POSTSUBSCRIPT ( italic_t ) in Theorem 3.2:

Theorem 3.2.

Let Dns(t)subscript𝐷𝑛𝑠𝑡D_{ns}(t)italic_D start_POSTSUBSCRIPT italic_n italic_s end_POSTSUBSCRIPT ( italic_t ) be the unmet demand in the non-surge zone at any time step t{0,1,2,}𝑡012t\in\{0,1,2,...\}italic_t ∈ { 0 , 1 , 2 , … }. Then, there always exists some time step τ0𝜏0\tau\geq 0italic_τ ≥ 0 such that Dns(t)subscript𝐷𝑛𝑠𝑡D_{ns}(t)italic_D start_POSTSUBSCRIPT italic_n italic_s end_POSTSUBSCRIPT ( italic_t ) is non-increasing for all tτ𝑡𝜏t\geq\tauitalic_t ≥ italic_τ and increasing for all t<τ𝑡𝜏t<\tauitalic_t < italic_τ. Further, τ𝜏\tauitalic_τ satisfies:

τ=min{t{0,1,2,.}:f(Ds(t)Dns(t))Ds(t)μλ}.\tau=\min\left\{t\in\{0,1,2,....\}:f\left(D_{s}(t)-D_{ns}(t)\right)\cdot D_{s}% (t)\leq\mu-\lambda\right\}.italic_τ = roman_min { italic_t ∈ { 0 , 1 , 2 , … . } : italic_f ( italic_D start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT ( italic_t ) - italic_D start_POSTSUBSCRIPT italic_n italic_s end_POSTSUBSCRIPT ( italic_t ) ) ⋅ italic_D start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT ( italic_t ) ≤ italic_μ - italic_λ } .

Finally, Dns(t)subscript𝐷𝑛𝑠𝑡D_{ns}(t)italic_D start_POSTSUBSCRIPT italic_n italic_s end_POSTSUBSCRIPT ( italic_t ) must converge to 00 in finite time and if τnsubscript𝜏𝑛\tau_{n}italic_τ start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT indicates the time to convergence, then τnsubscript𝜏𝑛\tau_{n}italic_τ start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT must satisfy:

d0(μλ)τnD0+d0(μλ).subscript𝑑0𝜇𝜆subscript𝜏𝑛subscript𝐷0subscript𝑑0𝜇𝜆\left\lfloor\frac{d_{0}}{(\mu-\lambda)}\right\rfloor\leq\tau_{n}\leq\left% \lceil\frac{D_{0}+d_{0}}{(\mu-\lambda)}\right\rceil.⌊ divide start_ARG italic_d start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_ARG start_ARG ( italic_μ - italic_λ ) end_ARG ⌋ ≤ italic_τ start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT ≤ ⌈ divide start_ARG italic_D start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT + italic_d start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_ARG start_ARG ( italic_μ - italic_λ ) end_ARG ⌉ .
Proof.

The detailed proof of all parts of the theorem can be found in Appendix A. ∎

This theorem provides the following insights: i) if τ=0𝜏0\tau=0italic_τ = 0, Dns(t)subscript𝐷𝑛𝑠𝑡D_{ns}(t)italic_D start_POSTSUBSCRIPT italic_n italic_s end_POSTSUBSCRIPT ( italic_t ) is non-increasing for all t𝑡titalic_t and converges to 00, and ii) if τ>0𝜏0\tau>0italic_τ > 0, Dns(t)subscript𝐷𝑛𝑠𝑡D_{ns}(t)italic_D start_POSTSUBSCRIPT italic_n italic_s end_POSTSUBSCRIPT ( italic_t ) increases initially to attain a peak at τ𝜏\tauitalic_τ, followed by which it decreases and converges to 00. It also provides interpretable upper and lower bounds on the time to convergence. The upper bound corresponds to the extreme case when all riders in the surge zone relocate to the non-surge zone in the first period itself and since the excess supply in each period is (μλ)𝜇𝜆(\mu-\lambda)( italic_μ - italic_λ ) in the non-surge zone, it takes about D0+d0(μλ)subscript𝐷0subscript𝑑0𝜇𝜆\frac{D_{0}+d_{0}}{(\mu-\lambda)}divide start_ARG italic_D start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT + italic_d start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_ARG start_ARG ( italic_μ - italic_λ ) end_ARG time steps to clear a total demand of size D0+d0subscript𝐷0subscript𝑑0D_{0}+d_{0}italic_D start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT + italic_d start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT. Similarly, the lower bound corresponds to the other extreme when riders are non-strategic and choose not to move, in which case it takes about d0(μλ)subscript𝑑0𝜇𝜆\frac{d_{0}}{(\mu-\lambda)}divide start_ARG italic_d start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_ARG start_ARG ( italic_μ - italic_λ ) end_ARG time steps to clear a demand of size d0subscript𝑑0d_{0}italic_d start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT. Thus, both upper and lower bounds are tight.

We have already noted that the value of τ𝜏\tauitalic_τ determines the shape of the demand curve Dns(t)subscript𝐷𝑛𝑠𝑡D_{ns}(t)italic_D start_POSTSUBSCRIPT italic_n italic_s end_POSTSUBSCRIPT ( italic_t ) in the non-surge zone. We will see shortly that the shape of Dns(t)subscript𝐷𝑛𝑠𝑡D_{ns}(t)italic_D start_POSTSUBSCRIPT italic_n italic_s end_POSTSUBSCRIPT ( italic_t ) can be used to characterize the type of surge.

Definition 3.1 (Localized Surge).

When the surge demand does not accumulate outside the surge boundary and the unmet non-surge demand Dns(t)subscript𝐷𝑛𝑠𝑡D_{ns}(t)italic_D start_POSTSUBSCRIPT italic_n italic_s end_POSTSUBSCRIPT ( italic_t ) is non-increasing for all t0𝑡0t\geq 0italic_t ≥ 0, we say that the surge is localized. Formally, τ=0𝜏0\tau=0italic_τ = 0 corresponds to the localized surge scenario.

Definition 3.2 (Spill-over Surge).

When the surge demand accumulates outside the surge zone causing the unmet non-surge demand Dns(t)subscript𝐷𝑛𝑠𝑡D_{ns}(t)italic_D start_POSTSUBSCRIPT italic_n italic_s end_POSTSUBSCRIPT ( italic_t ) to increase (even if for a short period of time), we say that the surge has spilled over. Formally, τ>0𝜏0\tau>0italic_τ > 0 corresponds to the spill-over surge scenario.

τ>0𝜏0\tau>0italic_τ > 0 represents a different scenario where up to time step τ𝜏\tauitalic_τ, Dnssubscript𝐷𝑛𝑠D_{ns}italic_D start_POSTSUBSCRIPT italic_n italic_s end_POSTSUBSCRIPT is increasing because the mass of riders relocating across the surge boundary is large enough that it cannot be cleared instantaneously by the excess supply μλ𝜇𝜆\mu-\lambdaitalic_μ - italic_λ, leading to accumulating demand. Since the high demand inside the surge zone appears to spill-over and affect demand outside the surge boundary, we call this type of surge a spill-over. In contrast, we call the case of τ=0𝜏0\tau=0italic_τ = 0 a localized surge, where the inflow of riders from the surge into the non-surge zone is less than the excess supply, and Dnssubscript𝐷𝑛𝑠D_{ns}italic_D start_POSTSUBSCRIPT italic_n italic_s end_POSTSUBSCRIPT is decreasing.

In the case of a spill-over surge, as Dns(t)subscript𝐷𝑛𝑠𝑡D_{ns}(t)italic_D start_POSTSUBSCRIPT italic_n italic_s end_POSTSUBSCRIPT ( italic_t ) increases, the demand differential Ds(t)Dns(t)subscript𝐷𝑠𝑡subscript𝐷𝑛𝑠𝑡D_{s}(t)-D_{ns}(t)italic_D start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT ( italic_t ) - italic_D start_POSTSUBSCRIPT italic_n italic_s end_POSTSUBSCRIPT ( italic_t ) decreases (recall, Ds(t)subscript𝐷𝑠𝑡D_{s}(t)italic_D start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT ( italic_t ) is steadily decreasing) and so does the incentive for riders to relocate. As a result, beyond τ𝜏\tauitalic_τ, Dns(t)subscript𝐷𝑛𝑠𝑡D_{ns}(t)italic_D start_POSTSUBSCRIPT italic_n italic_s end_POSTSUBSCRIPT ( italic_t ) starts to get cleared by the excess supply μλ𝜇𝜆\mu-\lambdaitalic_μ - italic_λ. This leads to the rapid decrease and eventual convergence of Dns(t)subscript𝐷𝑛𝑠𝑡D_{ns}(t)italic_D start_POSTSUBSCRIPT italic_n italic_s end_POSTSUBSCRIPT ( italic_t ). An interesting highlight of Theorem 3.2 is that Dns(t)subscript𝐷𝑛𝑠𝑡D_{ns}(t)italic_D start_POSTSUBSCRIPT italic_n italic_s end_POSTSUBSCRIPT ( italic_t ) has a relatively simple shape in that it is single peaked.

We now provide a necessary and sufficient condition on the initial demands D0subscript𝐷0D_{0}italic_D start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT and d0subscript𝑑0d_{0}italic_d start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT that guarantees τ=0𝜏0\tau=0italic_τ = 0 or equivalently, the absence of the spill-over effect. This formalizes the intuition expressed above about the relationship between rider relocation and excess supply:

Claim 3.3.

There is no spill-over, i.e., τ=0𝜏0\tau=0italic_τ = 0 if and only if the initial demands D0subscript𝐷0D_{0}italic_D start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT and d0subscript𝑑0d_{0}italic_d start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT satisfy D0f(D0d0)μλsubscript𝐷0𝑓subscript𝐷0subscript𝑑0𝜇𝜆D_{0}\cdot f\left(D_{0}-d_{0}\right)\leq\mu-\lambdaitalic_D start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ⋅ italic_f ( italic_D start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT - italic_d start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) ≤ italic_μ - italic_λ.

Proof.

The proof follows directly from the definition of τ𝜏\tauitalic_τ in Theorem 3.2. ∎

4 Numerical Experiments

In this section, we will start by visualizing the main theoretical insights that we have developed in Section 3 through numerical experiments—specifically the shapes of the demand curves under different surge types, whether the convergence times satisfy our theoretical bounds and how tight the bounds are. We then use simulations to extend our analysis to more general settings with stochastic arrival rates, explicit platform prices, strategic drivers reacting to prices, and a game-theoretic model of rider decision-making.

4.1 Deterministic Setting

First, we present numerical experiments for the simple deterministic setting captured in our model in Section 2 where the exogenous demand λ𝜆\lambdaitalic_λ and the supply of drivers μ𝜇\muitalic_μ in each zone is deterministic and fixed for all time periods. The main objectives for this segment are as follows: i) observing how the unmet demands in the two zones Ds(t)subscript𝐷𝑠𝑡D_{s}(t)italic_D start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT ( italic_t ) and Dns(t)subscript𝐷𝑛𝑠𝑡D_{ns}(t)italic_D start_POSTSUBSCRIPT italic_n italic_s end_POSTSUBSCRIPT ( italic_t ) evolve over time t𝑡titalic_t in localized and spill-over surge scenarios, ii) identifying and highlighting interesting phenomenon like surge inversions that can arise in special cases, and iii) investigating how fast the demands converge compared to the bounds predicted by theory and how the relative size of parameters μ𝜇\muitalic_μ and λ𝜆\lambdaitalic_λ affect this convergence.

For these experiments, we choose a simple functional form of f()𝑓f(\cdot)italic_f ( ⋅ ) which captures the fraction of waiting customers in the surge zone who choose to walk across the surge boundary in each time period. We use f(x)=max(0,kxμ)𝑓𝑥0𝑘𝑥𝜇f(x)=\max\left(0,\frac{kx}{\mu}\right)italic_f ( italic_x ) = roman_max ( 0 , divide start_ARG italic_k italic_x end_ARG start_ARG italic_μ end_ARG ) where k𝑘kitalic_k is a scalar multiplier. Therefore the rate of movement from the surge to the non-surge zone is given by

max(0,k(DsDns)μ).0𝑘subscript𝐷𝑠subscript𝐷𝑛𝑠𝜇\displaystyle\max\left(0,\frac{k(D_{s}-D_{ns})}{\mu}\right).roman_max ( 0 , divide start_ARG italic_k ( italic_D start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT - italic_D start_POSTSUBSCRIPT italic_n italic_s end_POSTSUBSCRIPT ) end_ARG start_ARG italic_μ end_ARG ) . (3)

The value of k𝑘kitalic_k is chosen appropriately to ensure that f()[0,1]𝑓01f(\cdot)\in[0,1]italic_f ( ⋅ ) ∈ [ 0 , 1 ]. A higher value of k𝑘kitalic_k indicates that a large fraction of customers choose to move. Recall that the input x𝑥xitalic_x to f()𝑓f(\cdot)italic_f ( ⋅ ) is the demand difference between the surge and non-surge zones (Assumption 2.1). We additionally use the supply μ𝜇\muitalic_μ as a normalizing factor, referring back to our pricing rule example in Example 2.1, and noting that intuitively, the higher the supply, the lesser the waiting times and the likelihood of a surge, and the less likely agents are to move. Finally, our justification for using a simple linear form for f()𝑓f(\cdot)italic_f ( ⋅ ) is that it is easy to implement and interpret, but we will show later in Section 4.3 that despite being simple, our linear model of f()𝑓f(\cdot)italic_f ( ⋅ ) has very good explaining power and can accurately capture the effect of moving riders due to high prices in the surge zone.

Refer to caption
Figure 2: Spill-over surge: Demand curves over time for D0=1000subscript𝐷01000D_{0}=1000italic_D start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = 1000, d0=200subscript𝑑0200d_{0}=200italic_d start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = 200, λ=30𝜆30\lambda=30italic_λ = 30 and μ=50𝜇50\mu=50italic_μ = 50. We set k=0.005𝑘0.005k=0.005italic_k = 0.005 and simulate over T=65𝑇65T=65italic_T = 65 time steps. Ds(t)subscript𝐷𝑠𝑡D_{s}(t)italic_D start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT ( italic_t ) converges in τs=33subscript𝜏𝑠33\tau_{s}=33italic_τ start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT = 33 time steps while Dns(t)subscript𝐷𝑛𝑠𝑡D_{ns}(t)italic_D start_POSTSUBSCRIPT italic_n italic_s end_POSTSUBSCRIPT ( italic_t ) converges in τn=30subscript𝜏𝑛30\tau_{n}=30italic_τ start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT = 30 time steps. Observe that they lie within our theoretically computed bounds in Theorems 3.1 and 3.2 (indicated by dotted lines).

In Figure 2, we demonstrate a spill-over surge scenario. Note that for the chosen values of the parameters D0subscript𝐷0D_{0}italic_D start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT, λ𝜆\lambdaitalic_λ, μ𝜇\muitalic_μ, we have D0f(D0d0)>μλsubscript𝐷0𝑓subscript𝐷0subscript𝑑0𝜇𝜆D_{0}\cdot f\left(D_{0}-d_{0}\right)>\mu-\lambdaitalic_D start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ⋅ italic_f ( italic_D start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT - italic_d start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) > italic_μ - italic_λ which we showed in Claim 3.3 to be a sufficient condition for surge to spill-over. This results in Dns(t)subscript𝐷𝑛𝑠𝑡D_{ns}(t)italic_D start_POSTSUBSCRIPT italic_n italic_s end_POSTSUBSCRIPT ( italic_t ) having a single peak at some t>0𝑡0t>0italic_t > 0 and then decreasing to 00. We also observe that Ds(t)subscript𝐷𝑠𝑡D_{s}(t)italic_D start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT ( italic_t ) decreases monotonically to 00, as expected.

Refer to caption
Figure 3: Localized surge: Demand curves over time for D0=1000subscript𝐷01000D_{0}=1000italic_D start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = 1000, d0=200subscript𝑑0200d_{0}=200italic_d start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = 200, λ=30𝜆30\lambda=30italic_λ = 30 and μ=50𝜇50\mu=50italic_μ = 50. We set k=0.001𝑘0.001k=0.001italic_k = 0.001 and simulate over T=65𝑇65T=65italic_T = 65 time steps. Ds(t)subscript𝐷𝑠𝑡D_{s}(t)italic_D start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT ( italic_t ) converges in τs=41subscript𝜏𝑠41\tau_{s}=41italic_τ start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT = 41 time steps while Dns(t)subscript𝐷𝑛𝑠𝑡D_{ns}(t)italic_D start_POSTSUBSCRIPT italic_n italic_s end_POSTSUBSCRIPT ( italic_t ) converges in τn=20subscript𝜏𝑛20\tau_{n}=20italic_τ start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT = 20 time steps (within theoretically computed bounds). Note the the distinctive difference in shape of Dns(t)subscript𝐷𝑛𝑠𝑡D_{ns}(t)italic_D start_POSTSUBSCRIPT italic_n italic_s end_POSTSUBSCRIPT ( italic_t ) from Figure 2.

Now, Figure 3 demonstrates a scenario with a low value of k𝑘kitalic_k—here the mass of riders relocating in each time period is always smaller than μλ𝜇𝜆\mu-\lambdaitalic_μ - italic_λ, so it can be matched instantaneously by the excess supply in the non-surge zone. This leads to a localized surge scenario (τ=0𝜏0\tau=0italic_τ = 0) where the unmet demands in both zones are found to decrease monotonically as predicted by theory. This also causes Dns(t)subscript𝐷𝑛𝑠𝑡D_{ns}(t)italic_D start_POSTSUBSCRIPT italic_n italic_s end_POSTSUBSCRIPT ( italic_t ) to converge much faster than Ds(t)subscript𝐷𝑠𝑡D_{s}(t)italic_D start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT ( italic_t ).

Surge Inversion:

We highlight the existence of cases where the non-surge demand Dns(t)subscript𝐷𝑛𝑠𝑡D_{ns}(t)italic_D start_POSTSUBSCRIPT italic_n italic_s end_POSTSUBSCRIPT ( italic_t ) exceeds the surge demand Ds(t)subscript𝐷𝑠𝑡D_{s}(t)italic_D start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT ( italic_t ). We call this phenomenon a surge inversion. An example of surge inversion is presented below in Figure 4. We have already seen that when the surge does not remain localized, Dns(t)subscript𝐷𝑛𝑠𝑡D_{ns}(t)italic_D start_POSTSUBSCRIPT italic_n italic_s end_POSTSUBSCRIPT ( italic_t ) increases initially and attains a peak before decreasing monotonically.

Refer to caption
Figure 4: Surge Inversion: Demand curves over time for D0=1000subscript𝐷01000D_{0}=1000italic_D start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = 1000, d0=300subscript𝑑0300d_{0}=300italic_d start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = 300 λ=30𝜆30\lambda=30italic_λ = 30, μ=50𝜇50\mu=50italic_μ = 50. We set k=0.05𝑘0.05k=0.05italic_k = 0.05. Observe that for this set of parameters, Dns(t)subscript𝐷𝑛𝑠𝑡D_{ns}(t)italic_D start_POSTSUBSCRIPT italic_n italic_s end_POSTSUBSCRIPT ( italic_t ) exceeds Ds(t)subscript𝐷𝑠𝑡D_{s}(t)italic_D start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT ( italic_t ) at t=1𝑡1t=1italic_t = 1.

However, in some cases (at high values of k𝑘kitalic_k), the mass of riders who choose to walk over in the first few time periods is so significant that the surge demand Ds(t)subscript𝐷𝑠𝑡D_{s}(t)italic_D start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT ( italic_t ) dips below the peak of Dns(t)subscript𝐷𝑛𝑠𝑡D_{ns}(t)italic_D start_POSTSUBSCRIPT italic_n italic_s end_POSTSUBSCRIPT ( italic_t ) leading to a surge inversion. As soon as inversion happens, the status quo is maintained for all future time steps because as per Assumption 2.1, riders will never choose to walk back inside the surge zone again even if the demand outside is higher555When riders are allowed to walk back, we expect to see cycles where the zone with the highest demand and highest price switches over time; we do not aim to characterize such cycles in this paper.. Henceforth, both zones will clear demand at a fixed rate of (μλ)𝜇𝜆(\mu-\lambda)( italic_μ - italic_λ ) until they converge to 00, implying that Ds(t)subscript𝐷𝑠𝑡D_{s}(t)italic_D start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT ( italic_t ) converges faster than Dns(t)subscript𝐷𝑛𝑠𝑡D_{ns}(t)italic_D start_POSTSUBSCRIPT italic_n italic_s end_POSTSUBSCRIPT ( italic_t ) during surge inversion. It is important to note that surge inversion does not violate any of our theoretical bounds on Ds(t)subscript𝐷𝑠𝑡D_{s}(t)italic_D start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT ( italic_t ) or Dns(t)subscript𝐷𝑛𝑠𝑡D_{ns}(t)italic_D start_POSTSUBSCRIPT italic_n italic_s end_POSTSUBSCRIPT ( italic_t ).

Variation of results over the parameter space:

Finally, we conclude the discussion on the deterministic setting with an analysis of how changing the value of the different parameters like D0subscript𝐷0D_{0}italic_D start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT, k𝑘kitalic_k, μ𝜇\muitalic_μ and λ𝜆\lambdaitalic_λ may affect the results: we are primarily interested in seeing i) how the system transitions between the two types of surge as a function of the parameters, and ii) how different parameters affect the time to convergence of the demands relative to each other.

In Figure 5, we fix the initial demands D0subscript𝐷0D_{0}italic_D start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT, d0subscript𝑑0d_{0}italic_d start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT and the parameter k𝑘kitalic_k associated with f()𝑓f(\cdot)italic_f ( ⋅ ) and only vary the exogenous demand λ𝜆\lambdaitalic_λ and supply rate μ𝜇\muitalic_μ. Firstly, as the magnitude of the excess supply (μλ𝜇𝜆\mu-\lambdaitalic_μ - italic_λ) increases, the gap between the convergence times τssubscript𝜏𝑠\tau_{s}italic_τ start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT and τnsubscript𝜏𝑛\tau_{n}italic_τ start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT increases — Dns(t)subscript𝐷𝑛𝑠𝑡D_{ns}(t)italic_D start_POSTSUBSCRIPT italic_n italic_s end_POSTSUBSCRIPT ( italic_t ) is found to converge much faster than Ds(t)subscript𝐷𝑠𝑡D_{s}(t)italic_D start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT ( italic_t ). The excess supply also determines the type of surge—with a smaller value of (μλ𝜇𝜆\mu-\lambdaitalic_μ - italic_λ), spill-over surges are more likely as long as the supply rate μ𝜇\muitalic_μ is not too large; this is in accordance with Claim 3.3. The size of the supply μ𝜇\muitalic_μ is also important independently of λ𝜆\lambdaitalic_λ, as it plays a role in the rate at which agents move, given in Equation 3. Because the initial demands D0subscript𝐷0D_{0}italic_D start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT and d0subscript𝑑0d_{0}italic_d start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT are fixed, a smaller value of μ𝜇\muitalic_μ indicates that the demand gap across the surge boundary relative to the supply rate is higher, which increases the incentive for riders to walk. This causes the Ds(t)subscript𝐷𝑠𝑡D_{s}(t)italic_D start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT ( italic_t ) and Dns(t)subscript𝐷𝑛𝑠𝑡D_{ns}(t)italic_D start_POSTSUBSCRIPT italic_n italic_s end_POSTSUBSCRIPT ( italic_t ) curves to become steeper and approach each other, as the demand spontaneously redistributes.

[][μ𝜇\muitalic_μ = 30, λ𝜆\lambdaitalic_λ = 5]Refer to caption [][μ𝜇\muitalic_μ = 30, λ𝜆\lambdaitalic_λ = 15]Refer to caption [][μ𝜇\muitalic_μ = 30, λ𝜆\lambdaitalic_λ = 25]Refer to caption

[][μ𝜇\muitalic_μ = 85, λ𝜆\lambdaitalic_λ = 10]Refer to caption [][μ𝜇\muitalic_μ = 85, λ𝜆\lambdaitalic_λ = 50]Refer to caption [][μ𝜇\muitalic_μ = 85, λ𝜆\lambdaitalic_λ = 70]Refer to caption

[][ λ𝜆\lambdaitalic_λ = 15, μ𝜇\muitalic_μ = 20]Refer to caption [][λ𝜆\lambdaitalic_λ = 15, μ𝜇\muitalic_μ = 40]Refer to caption [][λ𝜆\lambdaitalic_λ = 15, μ𝜇\muitalic_μ = 60]Refer to caption

[][λ𝜆\lambdaitalic_λ = 60, μ𝜇\muitalic_μ = 70]Refer to caption [][λ𝜆\lambdaitalic_λ = 60, μ𝜇\muitalic_μ = 90]Refer to caption [][λ𝜆\lambdaitalic_λ = 60, μ𝜇\muitalic_μ = 110]Refer to caption

Figure 5: Illustration of surge dynamics in the deterministic setting with variation in parameters μ𝜇\muitalic_μ and λ𝜆\lambdaitalic_λ. For all sets of results here, the following parameter values are fixed: D0=1000subscript𝐷01000D_{0}=1000italic_D start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = 1000, d0=200subscript𝑑0200d_{0}=200italic_d start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = 200, and k=0.005𝑘0.005k=0.005italic_k = 0.005.

In Figure 6, we explore how the demand curves evolve when i) the volume of the initial surge demand D0subscript𝐷0D_{0}italic_D start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT changes (all else fixed), and ii) the parameter k𝑘kitalic_k in f()𝑓f(\cdot)italic_f ( ⋅ ) changes (all else fixed). From sub-figures (a)-(b)-(c), it is clear that as D0subscript𝐷0D_{0}italic_D start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT increases, the surge type gradually transitions from the localized to the spill-over type which is highly intuitive. We make a similar observation in sub-figures (d)-(e)-(f) where D0subscript𝐷0D_{0}italic_D start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT is fixed, but k𝑘kitalic_k increases gradually. A higher k𝑘kitalic_k indicates that a higher fraction of riders choose to walk in each time step. Again, we see a gradual evolution in the shape of Dns(t)subscript𝐷𝑛𝑠𝑡D_{ns}(t)italic_D start_POSTSUBSCRIPT italic_n italic_s end_POSTSUBSCRIPT ( italic_t ) as it transitions from being monotonically decreasing (localized surge) to increasing first and then decreasing (spill-over surge).

Across both Figures 5 and 6, it is important to note how close the convergence times τssubscript𝜏𝑠\tau_{s}italic_τ start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT and τnsubscript𝜏𝑛\tau_{n}italic_τ start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT are to the bounds. In localized surge settings, τssubscript𝜏𝑠\tau_{s}italic_τ start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT is very close to its upper bound and τnsubscript𝜏𝑛\tau_{n}italic_τ start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT is very close to its lower bound. These observations can be explained by noting that in localized surge, the mass of riders walking across the surge boundary is not too high—as a result, the surge and non-surge demands only change mildly (surge demand decreases, non-surge demand increases) from the scenario where riders are non-strategic and do not walk. However, as we move to spill-over surge setting, our bounds become looser which is again along expected lines (in the case of surge demand, the upper bound completely ignores the effect of walking while the lower bound overestimates it by using the rate of walking in the first time step which is the highest).

[][D0subscript𝐷0D_{0}italic_D start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = 300, k=0.005𝑘0.005k=0.005italic_k = 0.005 ]Refer to caption [][D0subscript𝐷0D_{0}italic_D start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = 500, k=0.005𝑘0.005k=0.005italic_k = 0.005]Refer to caption [][D0subscript𝐷0D_{0}italic_D start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = 1500, k=0.005𝑘0.005k=0.005italic_k = 0.005]Refer to caption

[][k=0.0005𝑘0.0005k=0.0005italic_k = 0.0005, D0subscript𝐷0D_{0}italic_D start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = 1000]Refer to caption [][k=0.001𝑘0.001k=0.001italic_k = 0.001, D0subscript𝐷0D_{0}italic_D start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = 1000]Refer to caption [][k=0.01𝑘0.01k=0.01italic_k = 0.01, D0subscript𝐷0D_{0}italic_D start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = 1000]Refer to caption

Figure 6: Demonstration of how the shapes of Ds(t)subscript𝐷𝑠𝑡D_{s}(t)italic_D start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT ( italic_t ) and Dns(t)subscript𝐷𝑛𝑠𝑡D_{ns}(t)italic_D start_POSTSUBSCRIPT italic_n italic_s end_POSTSUBSCRIPT ( italic_t ) change as we traverse the parameter space in k𝑘kitalic_k and D0subscript𝐷0D_{0}italic_D start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT. All other parameters remain fixed: d0=200subscript𝑑0200d_{0}=200italic_d start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = 200, λ=30𝜆30\lambda=30italic_λ = 30 and μ=50𝜇50\mu=50italic_μ = 50.

4.2 Extension 1: Stochastic Demands & Supplies

One immediate extension is to consider stochastic exogenous demands and supplies for both zones at all time steps. Instead of having a fixed exogenous demand λ𝜆\lambdaitalic_λ and a fixed supply μ𝜇\muitalic_μ at each time step, we now assume that at any time t{0,1,2,.T}t\in\{0,1,2,....T\}italic_t ∈ { 0 , 1 , 2 , … . italic_T }, the exogenous demand is drawn from a distribution F()𝐹F(\cdot)italic_F ( ⋅ ) with mean number of arrivals per unit time λ𝜆\lambdaitalic_λ. Similarly, the supply is drawn from another distribution G()𝐺G(\cdot)italic_G ( ⋅ ) with mean number of arriving vehicles per unit time μ𝜇\muitalic_μ. We choose a Poisson distribution for F𝐹Fitalic_F and G𝐺Gitalic_G, a standard process for modeling stochastic arrival Kingman [1993].

Refer to caption
Figure 7: Spill-over surge under stochastic demands and supplies: Demand curves over time for D0=1000subscript𝐷01000D_{0}=1000italic_D start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = 1000, d0=200subscript𝑑0200d_{0}=200italic_d start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = 200, λ=30𝜆30\lambda=30italic_λ = 30 and μ=50𝜇50\mu=50italic_μ = 50. We set k=0.005𝑘0.005k=0.005italic_k = 0.005 and simulate over T=65𝑇65T=65italic_T = 65 time steps. For one realization of supplies and demands, Ds(t)subscript𝐷𝑠𝑡D_{s}(t)italic_D start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT ( italic_t ) converges in τs=31subscript𝜏𝑠31\tau_{s}=31italic_τ start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT = 31 time steps while Dns(t)subscript𝐷𝑛𝑠𝑡D_{ns}(t)italic_D start_POSTSUBSCRIPT italic_n italic_s end_POSTSUBSCRIPT ( italic_t ) converges in τn=29subscript𝜏𝑛29\tau_{n}=29italic_τ start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT = 29 time steps (compared to τs=33subscript𝜏𝑠33\tau_{s}=33italic_τ start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT = 33 and τn=30subscript𝜏𝑛30\tau_{n}=30italic_τ start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT = 30 in the deterministic case).
Refer to caption
Figure 8: Localized surge under stochastic demands and supplies: Demand curves over time for D0=1000subscript𝐷01000D_{0}=1000italic_D start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = 1000, d0=200subscript𝑑0200d_{0}=200italic_d start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = 200, λ=30𝜆30\lambda=30italic_λ = 30 and μ=50𝜇50\mu=50italic_μ = 50. We set k=0.001𝑘0.001k=0.001italic_k = 0.001 and simulate over T=65𝑇65T=65italic_T = 65 time steps. For one realization of supplies and demands, Ds(t)subscript𝐷𝑠𝑡D_{s}(t)italic_D start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT ( italic_t ) converges in τs=42subscript𝜏𝑠42\tau_{s}=42italic_τ start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT = 42 time steps while Dns(t)subscript𝐷𝑛𝑠𝑡D_{ns}(t)italic_D start_POSTSUBSCRIPT italic_n italic_s end_POSTSUBSCRIPT ( italic_t ) converges in τn=20subscript𝜏𝑛20\tau_{n}=20italic_τ start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT = 20 time steps (compared to τs=41subscript𝜏𝑠41\tau_{s}=41italic_τ start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT = 41 and τn=20subscript𝜏𝑛20\tau_{n}=20italic_τ start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT = 20 in the deterministic case).

We replicate our experiments with stochastic demands and supplies for the same parameter combinations as in the deterministic case (Figures 7 and 8). We can see that the trends broadly generalize from the deterministic to the stochastic case, in terms of convergence and shape properties of Ds(t)subscript𝐷𝑠𝑡D_{s}(t)italic_D start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT ( italic_t ) and Dns(t)subscript𝐷𝑛𝑠𝑡D_{ns}(t)italic_D start_POSTSUBSCRIPT italic_n italic_s end_POSTSUBSCRIPT ( italic_t ).

4.3 Extension 2: Agent-based, Game-theoretic simulation of ridesharing system with real-time market matching and passenger movement

Our theoretical results, as well as our experimental results of Section 4.2, rely on several strong assumptions. The first one is that, given the unmet demand Ds(t)subscript𝐷𝑠𝑡D_{s}(t)italic_D start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT ( italic_t ) in the surge zone and the unmet demand Dns(t)subscript𝐷𝑛𝑠𝑡D_{ns}(t)italic_D start_POSTSUBSCRIPT italic_n italic_s end_POSTSUBSCRIPT ( italic_t ) in the outside zone, the number of riders moving every day is a function of DsDnssubscript𝐷𝑠subscript𝐷𝑛𝑠D_{s}-D_{ns}italic_D start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT - italic_D start_POSTSUBSCRIPT italic_n italic_s end_POSTSUBSCRIPT. The second one is that drivers are not strategic and arrive in each zone at a rate that is independent of demands and prices, while drivers may chase surges in real life. In this section, we provide experiments showing that our theoretical highlights are robust to our theoretical modeling choice, highlighting how our model of Section 2 provides a reasonable, first-order approximation of the main effect we are modeling. Our experiments rely on an agent-based, game-theoretic model of how riders decide to walk out of the surge zone and how drivers decide between the surge and non-surge zones.

4.3.1 Game-theoretic Model.

Below, we describe our game-theoretic, game-theoretic model of interactions between riders, drivers, and the platform.

Driver Model:

We consider a finite number of strategic drivers. Drivers arrive in each time step at a rate of 2μ2𝜇2\mu2 italic_μ, according to a Poisson distribution666In this sense, we do not model drivers who strategically decide when to come online.. However, drivers can strategically choose which of the surge and non-surge zones to serve. Given price multipliers Pssubscript𝑃𝑠P_{s}italic_P start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT and Pnssubscript𝑃𝑛𝑠P_{ns}italic_P start_POSTSUBSCRIPT italic_n italic_s end_POSTSUBSCRIPT in the surge and non-surge zones respectively, each driver decides which zone to go to according to a multinomial logit model, a seminal and practically well-motivated model of decision-making Hausman and McFadden [1984]: i.e., a driver goes to the surge zone with probability

γs=exp(τPs)exp(τPs)+exp(τPns),subscript𝛾𝑠𝜏subscript𝑃𝑠𝜏subscript𝑃𝑠𝜏subscript𝑃𝑛𝑠\gamma_{s}=\frac{\exp(\tau P_{s})}{\exp(\tau P_{s})+\exp(\tau P_{ns})},italic_γ start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT = divide start_ARG roman_exp ( italic_τ italic_P start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT ) end_ARG start_ARG roman_exp ( italic_τ italic_P start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT ) + roman_exp ( italic_τ italic_P start_POSTSUBSCRIPT italic_n italic_s end_POSTSUBSCRIPT ) end_ARG ,

and to the non-surge zone with probability

γns=exp(τPns)exp(τPs)+exp(τPns),subscript𝛾𝑛𝑠𝜏subscript𝑃𝑛𝑠𝜏subscript𝑃𝑠𝜏subscript𝑃𝑛𝑠\gamma_{ns}=\frac{\exp(\tau P_{ns})}{\exp(\tau P_{s})+\exp(\tau P_{ns})},italic_γ start_POSTSUBSCRIPT italic_n italic_s end_POSTSUBSCRIPT = divide start_ARG roman_exp ( italic_τ italic_P start_POSTSUBSCRIPT italic_n italic_s end_POSTSUBSCRIPT ) end_ARG start_ARG roman_exp ( italic_τ italic_P start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT ) + roman_exp ( italic_τ italic_P start_POSTSUBSCRIPT italic_n italic_s end_POSTSUBSCRIPT ) end_ARG ,

for some positive parameter τ>0𝜏0\tau>0italic_τ > 0.

Rider Model:

We consider a finite number of riders, arriving at a rate λ𝜆\lambdaitalic_λ in each zone, according to a Poisson distribution. There is also an initial number of riders D0subscript𝐷0D_{0}italic_D start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT in the surge zone and d0subscript𝑑0d_{0}italic_d start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT in the non-surge zone.

Each rider has a willingness to pay, denoted by v𝑣vitalic_v, which determines whether they choose to accept a ride, given a price p𝑝pitalic_p. If vp𝑣𝑝v\geq pitalic_v ≥ italic_p, then the rider can be matched with a ride at the current price, otherwise they continue to wait. Each rider is also associated with a cost c𝑐citalic_c for moving from the surge zone to the non-surge zone. The costs c𝑐citalic_c are drawn i.i.d. from a distribution 𝒟𝒟\mathcal{D}caligraphic_D with cumulative distribution function F(.)F(.)italic_F ( . ). In this case, if the prices set by the platform are Pnssubscript𝑃𝑛𝑠P_{ns}italic_P start_POSTSUBSCRIPT italic_n italic_s end_POSTSUBSCRIPT and Pssubscript𝑃𝑠P_{s}italic_P start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT, a user with cost c𝑐citalic_c in the surge zone walks to the non-surge zone if and only if

cPsPns.𝑐subscript𝑃𝑠subscript𝑃𝑛𝑠c\leq P_{s}-P_{ns}.italic_c ≤ italic_P start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT - italic_P start_POSTSUBSCRIPT italic_n italic_s end_POSTSUBSCRIPT .

In particular, the probability that a given user decides to walk to the non-surge zone is then given by F(PsPns)𝐹subscript𝑃𝑠subscript𝑃𝑛𝑠F(P_{s}-P_{ns})italic_F ( italic_P start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT - italic_P start_POSTSUBSCRIPT italic_n italic_s end_POSTSUBSCRIPT ). In our experiments, we choose 𝒟𝒟\mathcal{D}caligraphic_D to be a symmetric truncated Gaussian distribution (left-truncated at 00) for which the target mean 𝒟meansubscript𝒟𝑚𝑒𝑎𝑛\mathcal{D}_{mean}caligraphic_D start_POSTSUBSCRIPT italic_m italic_e italic_a italic_n end_POSTSUBSCRIPT and standard deviation 𝒟stdsubscript𝒟𝑠𝑡𝑑\mathcal{D}_{std}caligraphic_D start_POSTSUBSCRIPT italic_s italic_t italic_d end_POSTSUBSCRIPT are specified. 𝒟meansubscript𝒟𝑚𝑒𝑎𝑛\mathcal{D}_{mean}caligraphic_D start_POSTSUBSCRIPT italic_m italic_e italic_a italic_n end_POSTSUBSCRIPT and 𝒟stdsubscript𝒟𝑠𝑡𝑑\mathcal{D}_{std}caligraphic_D start_POSTSUBSCRIPT italic_s italic_t italic_d end_POSTSUBSCRIPT directly influence what fraction of the current surge demand Ds(t)subscript𝐷𝑠𝑡D_{s}(t)italic_D start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT ( italic_t ) chooses to relocate across the surge boundary.

Platform Prices:

We consider a platform using a uniform price multiplier Pssubscript𝑃𝑠P_{s}italic_P start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT in the surge zone and Pnssubscript𝑃𝑛𝑠P_{ns}italic_P start_POSTSUBSCRIPT italic_n italic_s end_POSTSUBSCRIPT in the non-surge zone. We are in particular interested in studying the price gap ΔP=PsPnsΔ𝑃subscript𝑃𝑠subscript𝑃𝑛𝑠\Delta P=P_{s}-P_{ns}roman_Δ italic_P = italic_P start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT - italic_P start_POSTSUBSCRIPT italic_n italic_s end_POSTSUBSCRIPT. Note that ΔPΔ𝑃\Delta Proman_Δ italic_P represents an additive surge premium over the non-surge price; Garg and Nazerzadeh [2022] have shown that additive surge pricing has desirable properties, such as being incentive-compatible for drivers. Additionally,  Ma et al. [2022] demonstrate that spatial mispricing (i.e., prices in adjacent areas being vastly different) during price surges incentivizes drivers to “cherry-pick” rides and is a leading cause of market failure. Therefore, guaranteeing low ΔPΔ𝑃\Delta Proman_Δ italic_P across adjacent zones is highly desirable, and we will show in this section how rider strategic behavior contributes to reducing this gap ΔPΔ𝑃\Delta Proman_Δ italic_P.

We derive ΔPΔ𝑃\Delta Proman_Δ italic_P following the approach from the previous work of Yu et al. [2023], the platform set prices to satisfy a simple equilibrium condition: namely, the wait times in both zones should be roughly the same. Here, we consider a slight simplification of this condition, where we use the ratio of demand to supply as a proxy for wait times.777If there are D𝐷Ditalic_D riders total to clear, S𝑆Sitalic_S drivers come at every time step, and DS𝐷𝑆D\geq Sitalic_D ≥ italic_S, then the average wait time is D/S1𝐷𝑆1D/S-1italic_D / italic_S - 1. Letting Rs(t)subscript𝑅𝑠𝑡R_{s}(t)italic_R start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT ( italic_t ) and Rns(t)subscript𝑅𝑛𝑠𝑡R_{ns}(t)italic_R start_POSTSUBSCRIPT italic_n italic_s end_POSTSUBSCRIPT ( italic_t ) be the average number of drivers going to the surge and non-surge zones (respectively) at time t𝑡titalic_t, the platform’s goal is to set equilibrium prices to induce

Ds(t)Rs(t)=Dns(t)Rns(t).subscript𝐷𝑠𝑡subscript𝑅𝑠𝑡subscript𝐷𝑛𝑠𝑡subscript𝑅𝑛𝑠𝑡\frac{D_{s}(t)}{R_{s}(t)}=\frac{D_{ns}(t)}{R_{ns}(t)}.divide start_ARG italic_D start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT ( italic_t ) end_ARG start_ARG italic_R start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT ( italic_t ) end_ARG = divide start_ARG italic_D start_POSTSUBSCRIPT italic_n italic_s end_POSTSUBSCRIPT ( italic_t ) end_ARG start_ARG italic_R start_POSTSUBSCRIPT italic_n italic_s end_POSTSUBSCRIPT ( italic_t ) end_ARG .

Remember that Rs(t)subscript𝑅𝑠𝑡R_{s}(t)italic_R start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT ( italic_t ) and Rns(t)subscript𝑅𝑛𝑠𝑡R_{ns}(t)italic_R start_POSTSUBSCRIPT italic_n italic_s end_POSTSUBSCRIPT ( italic_t ) are a function of the price multipliers, following a multinomial-logit-based decision process; given an arrival rate of 2μ2𝜇2\mu2 italic_μ, we have that in expectation,

Rs(t)=2μexp(τPs)exp(τPs)+exp(τPns),Rns(t)=2μexp(τPns)exp(τPs)+exp(τPns).formulae-sequencesubscript𝑅𝑠𝑡2𝜇𝜏subscript𝑃𝑠𝜏subscript𝑃𝑠𝜏subscript𝑃𝑛𝑠subscript𝑅𝑛𝑠𝑡2𝜇𝜏subscript𝑃𝑛𝑠𝜏subscript𝑃𝑠𝜏subscript𝑃𝑛𝑠R_{s}(t)=2\mu\cdot\frac{\exp(\tau P_{s})}{\exp(\tau P_{s})+\exp(\tau P_{ns})},% \leavevmode\nobreak\ \leavevmode\nobreak\ R_{ns}(t)=2\mu\cdot\frac{\exp(\tau P% _{ns})}{\exp(\tau P_{s})+\exp(\tau P_{ns})}.italic_R start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT ( italic_t ) = 2 italic_μ ⋅ divide start_ARG roman_exp ( italic_τ italic_P start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT ) end_ARG start_ARG roman_exp ( italic_τ italic_P start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT ) + roman_exp ( italic_τ italic_P start_POSTSUBSCRIPT italic_n italic_s end_POSTSUBSCRIPT ) end_ARG , italic_R start_POSTSUBSCRIPT italic_n italic_s end_POSTSUBSCRIPT ( italic_t ) = 2 italic_μ ⋅ divide start_ARG roman_exp ( italic_τ italic_P start_POSTSUBSCRIPT italic_n italic_s end_POSTSUBSCRIPT ) end_ARG start_ARG roman_exp ( italic_τ italic_P start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT ) + roman_exp ( italic_τ italic_P start_POSTSUBSCRIPT italic_n italic_s end_POSTSUBSCRIPT ) end_ARG .

The equilibrium condition for the price multipliers then becomes

exp(τPns)Ds(t)Dns(t)=exp(τPs),𝜏subscript𝑃𝑛𝑠subscript𝐷𝑠𝑡subscript𝐷𝑛𝑠𝑡𝜏subscript𝑃𝑠\exp(\tau P_{ns})\cdot\frac{D_{s}(t)}{D_{ns}(t)}=\exp(\tau P_{s}),roman_exp ( italic_τ italic_P start_POSTSUBSCRIPT italic_n italic_s end_POSTSUBSCRIPT ) ⋅ divide start_ARG italic_D start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT ( italic_t ) end_ARG start_ARG italic_D start_POSTSUBSCRIPT italic_n italic_s end_POSTSUBSCRIPT ( italic_t ) end_ARG = roman_exp ( italic_τ italic_P start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT ) ,

yielding the following closed-form expression for the price gap across the surge boundary:

ΔP=PsPns=1τlog(Ds(t)Dns(t)).Δ𝑃subscript𝑃𝑠subscript𝑃𝑛𝑠1𝜏subscript𝐷𝑠𝑡subscript𝐷𝑛𝑠𝑡\displaystyle\Delta P=P_{s}-P_{ns}=\frac{1}{\tau}\log\left(\frac{D_{s}(t)}{D_{% ns}(t)}\right).roman_Δ italic_P = italic_P start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT - italic_P start_POSTSUBSCRIPT italic_n italic_s end_POSTSUBSCRIPT = divide start_ARG 1 end_ARG start_ARG italic_τ end_ARG roman_log ( divide start_ARG italic_D start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT ( italic_t ) end_ARG start_ARG italic_D start_POSTSUBSCRIPT italic_n italic_s end_POSTSUBSCRIPT ( italic_t ) end_ARG ) . (4)

When the total unmet demand Ds(t)+Dns(t)subscript𝐷𝑠𝑡subscript𝐷𝑛𝑠𝑡D_{s}(t)+D_{ns}(t)italic_D start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT ( italic_t ) + italic_D start_POSTSUBSCRIPT italic_n italic_s end_POSTSUBSCRIPT ( italic_t ) falls below the total supply, we note that the surge is cleared, and in this case ΔP=0Δ𝑃0\Delta P=0roman_Δ italic_P = 0 (both zones return to normal pricing).

Finally, platforms usually also need to impose a cap M𝑀Mitalic_M on the surge price either due to regulatory interventions  Lien [2018] or to make sure that prices do not grow unbounded during emergencies  News [2015]. We incorporate the same into our pricing model.

4.3.2 Experimental results.

In this segment, we provide simulation results using the agent model described above where we capture how the unmet demands Ds(t)subscript𝐷𝑠𝑡D_{s}(t)italic_D start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT ( italic_t ) and Dns(t)subscript𝐷𝑛𝑠𝑡D_{ns}(t)italic_D start_POSTSUBSCRIPT italic_n italic_s end_POSTSUBSCRIPT ( italic_t ) evolve over time. We compare these results directly against:

  • Results from our theoretical model (Theory); and

  • A baseline where riders are not strategic, called Non-Strategic Benchmark (NSB). This baseline uses the game-theoretic model described in Section 4.3.1, but with the caveat that riders never move from the surge to the non-surge zone. Equivalently, riders have a large or infinite cost c𝑐citalic_c for moving across zones.

Henceforth, for brevity, we will refer to the strategic agent model as (SA), the theoretical model as (Theory) and the non-strategic agent benchmark model as (NSB).

Comparison to Theory model:

We explore the explaining power of our theory model with respect to the strategic agent simulation model (SA). In our theory model, the effect of agents relocating across the surge boundary is modelled using a function f()𝑓f(\cdot)italic_f ( ⋅ ) which depends on the demand differential between the surge and the non-surge zones. Recall that we have assumed a simple linear form of f()𝑓f(\cdot)italic_f ( ⋅ ) given by f(x)=kxμ𝑓𝑥𝑘𝑥𝜇f(x)=\frac{kx}{\mu}italic_f ( italic_x ) = divide start_ARG italic_k italic_x end_ARG start_ARG italic_μ end_ARG. In Figure 9, we present the demand dynamics of our theory model for the following set of parameters (D0=2000subscript𝐷02000D_{0}=2000italic_D start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = 2000, d0=250subscript𝑑0250d_{0}=250italic_d start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = 250, λ=30𝜆30\lambda=30italic_λ = 30, μ=50𝜇50\mu=50italic_μ = 50, 𝒟mean=15subscript𝒟𝑚𝑒𝑎𝑛15\mathcal{D}_{mean}=15caligraphic_D start_POSTSUBSCRIPT italic_m italic_e italic_a italic_n end_POSTSUBSCRIPT = 15, 𝒟std=8subscript𝒟𝑠𝑡𝑑8\mathcal{D}_{std}=8caligraphic_D start_POSTSUBSCRIPT italic_s italic_t italic_d end_POSTSUBSCRIPT = 8) for different values of k𝑘kitalic_k. Observe that k𝑘kitalic_k can be tuned to emulate the SA model very closely, which we show in Figure 10.

[][k=0.0005𝑘0.0005k=0.0005italic_k = 0.0005]Refer to caption [][k=0.0010𝑘0.0010k=0.0010italic_k = 0.0010]Refer to caption

[][k=0.0025𝑘0.0025k=0.0025italic_k = 0.0025]Refer to caption [][k=0.0040𝑘0.0040k=0.0040italic_k = 0.0040]Refer to caption

Figure 9: Comparison of SA model with different instances of Theory model, each with a different value of k𝑘kitalic_k.
Refer to caption
Figure 10: At k=0.0016𝑘0.0016k=0.0016italic_k = 0.0016, we observe that the Theory model curves closely match the SA simulation model curves. This is for the following set of parameters: D0=2000subscript𝐷02000D_{0}=2000italic_D start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = 2000, d0=250subscript𝑑0250d_{0}=250italic_d start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = 250, λ=30𝜆30\lambda=30italic_λ = 30, μ=50𝜇50\mu=50italic_μ = 50, 𝒟mean=15subscript𝒟𝑚𝑒𝑎𝑛15\mathcal{D}_{mean}=15caligraphic_D start_POSTSUBSCRIPT italic_m italic_e italic_a italic_n end_POSTSUBSCRIPT = 15 and 𝒟std=8subscript𝒟𝑠𝑡𝑑8\mathcal{D}_{std}=8caligraphic_D start_POSTSUBSCRIPT italic_s italic_t italic_d end_POSTSUBSCRIPT = 8. For any set of parameters, the optimal k𝑘kitalic_k can be tuned.
Comparison to Non-Strategic Benchmark Model:

We also compare the price curves and demand curves of our game-theoretic simulation model with strategic agents (SA) with its non-strategic counterpart (NSB) in Figure 11.

[][Demand Curves]Refer to caption [][ΔPΔ𝑃\Delta Proman_Δ italic_P curves]Refer to caption

Figure 11: Comparison of demand curves (a) and price curves (b) between the Strategic Agent Model (SA) and the Non-Strategic Agent Benchmark Model (NSB) with the following parameters: D0=2000subscript𝐷02000D_{0}=2000italic_D start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = 2000, d0=250subscript𝑑0250d_{0}=250italic_d start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = 250, λ=30𝜆30\lambda=30italic_λ = 30, μ=45𝜇45\mu=45italic_μ = 45, 𝒟mean=8subscript𝒟𝑚𝑒𝑎𝑛8\mathcal{D}_{mean}=8caligraphic_D start_POSTSUBSCRIPT italic_m italic_e italic_a italic_n end_POSTSUBSCRIPT = 8, 𝒟std=8subscript𝒟𝑠𝑡𝑑8\mathcal{D}_{std}=8caligraphic_D start_POSTSUBSCRIPT italic_s italic_t italic_d end_POSTSUBSCRIPT = 8. Note that 𝒟𝒟\mathcal{D}caligraphic_D represents the cost distribution of agents. The willingness to pay values of agents are drawn i.i.d. from a symmetric truncated normal distribution with mean =7absent7=7= 7 and standard deviation =2absent2=2= 2. Time to convergence τ=82𝜏82\tau=82italic_τ = 82 time steps.

From the demand curves in (a), it is clear that the SA model leads to a redistribution of demand due to the effect of strategic agents moving from the surge zone to the non-surge zone. This also results in a more desirable price curve for the SA model with the price gap ΔPΔ𝑃\Delta Proman_Δ italic_P across the surge boundary diminishing and returning to normal levels faster than the NSB model.

The heatmap in Figure 12 illustrates how much improvement the SA model achieves in terms of price gap ΔPΔ𝑃\Delta Proman_Δ italic_P relative to the NSB model for different combinations of parameters 𝒟meansubscript𝒟𝑚𝑒𝑎𝑛\mathcal{D}_{mean}caligraphic_D start_POSTSUBSCRIPT italic_m italic_e italic_a italic_n end_POSTSUBSCRIPT and 𝒟stdsubscript𝒟𝑠𝑡𝑑\mathcal{D}_{std}caligraphic_D start_POSTSUBSCRIPT italic_s italic_t italic_d end_POSTSUBSCRIPT. The blue regions in the heatmap represent high improvement regimes while the purple regions represent low improvement regimes. For example, there is an average 75% improvement in SA price gaps compared to NSB for the parameter combinations used in Figure 11. Observe that the highest improvements are achieved at low values of 𝒟meansubscript𝒟𝑚𝑒𝑎𝑛\mathcal{D}_{mean}caligraphic_D start_POSTSUBSCRIPT italic_m italic_e italic_a italic_n end_POSTSUBSCRIPT. This is because low mean values of the waiting cost distribution indicate that a large fraction of the riders have an incentive to cross the surge boundary. This leads to demand redistribution and lowering of the price gap. On the other hand, at high values of 𝒟meansubscript𝒟𝑚𝑒𝑎𝑛\mathcal{D}_{mean}caligraphic_D start_POSTSUBSCRIPT italic_m italic_e italic_a italic_n end_POSTSUBSCRIPT with low 𝒟stdsubscript𝒟𝑠𝑡𝑑\mathcal{D}_{std}caligraphic_D start_POSTSUBSCRIPT italic_s italic_t italic_d end_POSTSUBSCRIPT, we see the smallest levels of improvement. This is again intuitive as a high mean waiting cost coupled with a small standard deviation implies that most riders have high waiting costs and are unwilling to move. Therefore, the demand remains significantly higher in the surge zone for longer periods of time, leading to higher price gaps on average, and marginal improvements compared to the NSB model. Finally, at high values of 𝒟stdsubscript𝒟𝑠𝑡𝑑\mathcal{D}_{std}caligraphic_D start_POSTSUBSCRIPT italic_s italic_t italic_d end_POSTSUBSCRIPT, there is a mix of riders with low and high waiting costs due to the high variability in the distribution of costs. As a result, the price gaps become invariant of 𝒟meansubscript𝒟𝑚𝑒𝑎𝑛\mathcal{D}_{mean}caligraphic_D start_POSTSUBSCRIPT italic_m italic_e italic_a italic_n end_POSTSUBSCRIPT and stabilize (between 50%percent50-50\leavevmode\nobreak\ \%- 50 % and 70%percent70-70\leavevmode\nobreak\ \%- 70 %) as we move towards the right of the map.

Refer to caption
Figure 12: Relative Difference in ΔPΔ𝑃\Delta Proman_Δ italic_P between the SA and NSB models, given by ΔPSAΔPNSBΔPNSB×100Δsubscript𝑃𝑆𝐴Δsubscript𝑃𝑁𝑆𝐵Δsubscript𝑃𝑁𝑆𝐵100\frac{\Delta P_{SA}-\Delta P_{NSB}}{\Delta P_{NSB}}\times 100divide start_ARG roman_Δ italic_P start_POSTSUBSCRIPT italic_S italic_A end_POSTSUBSCRIPT - roman_Δ italic_P start_POSTSUBSCRIPT italic_N italic_S italic_B end_POSTSUBSCRIPT end_ARG start_ARG roman_Δ italic_P start_POSTSUBSCRIPT italic_N italic_S italic_B end_POSTSUBSCRIPT end_ARG × 100, for different combinations of 𝒟meansubscript𝒟𝑚𝑒𝑎𝑛\mathcal{D}_{mean}caligraphic_D start_POSTSUBSCRIPT italic_m italic_e italic_a italic_n end_POSTSUBSCRIPT and 𝒟stdsubscript𝒟𝑠𝑡𝑑\mathcal{D}_{std}caligraphic_D start_POSTSUBSCRIPT italic_s italic_t italic_d end_POSTSUBSCRIPT. Reported values are time averages over T=500𝑇500T=500italic_T = 500 time steps. A more negative ratio indicates a higher degree of improvement. Demand and supply parameters are given by D0=2000subscript𝐷02000D_{0}=2000italic_D start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = 2000, d0=250subscript𝑑0250d_{0}=250italic_d start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = 250, μ=45𝜇45\mu=45italic_μ = 45 and λ=30𝜆30\lambda=30italic_λ = 30.
Effects of different parameters on the demand and price curves:

We now evaluate our game-theoretic simulation models for different parameter combinations where the parameters of interest are the exogenous demand and supply in each zone, the mean and scale parameters of the agents’ cost distribution 𝒟𝒟\mathcal{D}caligraphic_D and the initial demands D0subscript𝐷0D_{0}italic_D start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT and d0subscript𝑑0d_{0}italic_d start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT.

  • Ratio of Supply to Exogenous Demand (μλ)𝜇𝜆\left(\frac{\mu}{\lambda}\right)( divide start_ARG italic_μ end_ARG start_ARG italic_λ end_ARG ): For this set of experiments, we fix D0=2000subscript𝐷02000D_{0}=2000italic_D start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = 2000, d0=250subscript𝑑0250d_{0}=250italic_d start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = 250, 𝒟mean=12subscript𝒟𝑚𝑒𝑎𝑛12\mathcal{D}_{mean}=12caligraphic_D start_POSTSUBSCRIPT italic_m italic_e italic_a italic_n end_POSTSUBSCRIPT = 12, 𝒟std=8subscript𝒟𝑠𝑡𝑑8\mathcal{D}_{std}=8caligraphic_D start_POSTSUBSCRIPT italic_s italic_t italic_d end_POSTSUBSCRIPT = 8 and λ=30𝜆30\lambda=30italic_λ = 30. We vary the supply rate μ𝜇\muitalic_μ to evaluate the performance of the system for different supply-to-demand ratios.

    [][μλ=1.25𝜇𝜆1.25\frac{\mu}{\lambda}=1.25divide start_ARG italic_μ end_ARG start_ARG italic_λ end_ARG = 1.25, τ=170𝜏170\tau=170italic_τ = 170]Refer to caption [][μλ=1.8𝜇𝜆1.8\frac{\mu}{\lambda}=1.8divide start_ARG italic_μ end_ARG start_ARG italic_λ end_ARG = 1.8, τ=45𝜏45\tau=45italic_τ = 45]Refer to caption [][μλ=2.5𝜇𝜆2.5\frac{\mu}{\lambda}=2.5divide start_ARG italic_μ end_ARG start_ARG italic_λ end_ARG = 2.5, τ=24𝜏24\tau=24italic_τ = 24]Refer to caption

    [][μλ=1.25𝜇𝜆1.25\frac{\mu}{\lambda}=1.25divide start_ARG italic_μ end_ARG start_ARG italic_λ end_ARG = 1.25]Refer to caption [][μλ=1.8𝜇𝜆1.8\frac{\mu}{\lambda}=1.8divide start_ARG italic_μ end_ARG start_ARG italic_λ end_ARG = 1.8]Refer to caption [][μλ=2.5𝜇𝜆2.5\frac{\mu}{\lambda}=2.5divide start_ARG italic_μ end_ARG start_ARG italic_λ end_ARG = 2.5]Refer to caption

    Figure 13: Price and demand curves for SA and NSB models for different supply-demand ratios.

    As we increase the supply μ𝜇\muitalic_μ relative to the exogenous demand λ𝜆\lambdaitalic_λ, the excess supply in each time period increases, leading to faster dissipation of surge and restoration of prices and unmet demands to normal levels. This can be observed by the decreasing values of τ𝜏\tauitalic_τ in Figure 13 as we go from left to right.

  • Ratio of Initial Demands D0subscript𝐷0D_{0}italic_D start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT and d0subscript𝑑0d_{0}italic_d start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT: For this set of experiments, we fix d0=250subscript𝑑0250d_{0}=250italic_d start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = 250, 𝒟mean=12subscript𝒟𝑚𝑒𝑎𝑛12\mathcal{D}_{mean}=12caligraphic_D start_POSTSUBSCRIPT italic_m italic_e italic_a italic_n end_POSTSUBSCRIPT = 12, 𝒟std=8subscript𝒟𝑠𝑡𝑑8\mathcal{D}_{std}=8caligraphic_D start_POSTSUBSCRIPT italic_s italic_t italic_d end_POSTSUBSCRIPT = 8 and λ=30𝜆30\lambda=30italic_λ = 30 and μ=60𝜇60\mu=60italic_μ = 60. D0subscript𝐷0D_{0}italic_D start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT is varied to capture different relative sizes of uncleared demands in the two zones. The ratio of D0subscript𝐷0D_{0}italic_D start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT to d0subscript𝑑0d_{0}italic_d start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT represents the magnitude of the demand disparity caused by the surge. As D0subscript𝐷0D_{0}italic_D start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT increases, the maximum ΔPΔ𝑃\Delta Proman_Δ italic_P is found to increase. This is because the large demand disparity forces the platform to set much higher prices in the surge zone to capture a larger proportion of the total supply and still equalize waiting time across the two zones. Additionally, since μ𝜇\muitalic_μ and λ𝜆\lambdaitalic_λ remain unchanged, the total demand clearing rate in the system is still the same but there is larger demand to be cleared. So, the system takes much longer to dissipate the surge, i.e., τ𝜏\tauitalic_τ increases from left to right. Refer to Figure 14.

    [][D0d0=5subscript𝐷0subscript𝑑05\frac{D_{0}}{d_{0}}=5divide start_ARG italic_D start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_ARG start_ARG italic_d start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_ARG = 5, τ=24𝜏24\tau=24italic_τ = 24]Refer to caption [][D0d0=10subscript𝐷0subscript𝑑010\frac{D_{0}}{d_{0}}=10divide start_ARG italic_D start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_ARG start_ARG italic_d start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_ARG = 10, τ=47𝜏47\tau=47italic_τ = 47]Refer to caption [][D0d0=20subscript𝐷0subscript𝑑020\frac{D_{0}}{d_{0}}=20divide start_ARG italic_D start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_ARG start_ARG italic_d start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_ARG = 20, τ=91𝜏91\tau=91italic_τ = 91]Refer to caption

    [][D0d0=5subscript𝐷0subscript𝑑05\frac{D_{0}}{d_{0}}=5divide start_ARG italic_D start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_ARG start_ARG italic_d start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_ARG = 5]Refer to caption [][D0d0=10subscript𝐷0subscript𝑑010\frac{D_{0}}{d_{0}}=10divide start_ARG italic_D start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_ARG start_ARG italic_d start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_ARG = 10]Refer to caption [][D0d0=20subscript𝐷0subscript𝑑020\frac{D_{0}}{d_{0}}=20divide start_ARG italic_D start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_ARG start_ARG italic_d start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_ARG = 20]Refer to caption

    Figure 14: Price and demand curves for SA and NSB models for different ratios of initial demands.
  • Effect of the mean parameter of distribution 𝒟𝒟\mathcal{D}caligraphic_D: We can vary parameters of agents’ cost distribution 𝒟𝒟\mathcal{D}caligraphic_D keeping all other system parameters fixed with values: D0=2000subscript𝐷02000D_{0}=2000italic_D start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = 2000, d0=250subscript𝑑0250d_{0}=250italic_d start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = 250, μ=60𝜇60\mu=60italic_μ = 60, λ=30𝜆30\lambda=30italic_λ = 30. Recall that the parameters of 𝒟𝒟\mathcal{D}caligraphic_D regulate what fraction of agents would choose to relocate in each time period as a function of the price gap ΔPΔ𝑃\Delta Proman_Δ italic_P. For the results in Figure 15, we fix 𝒟stdsubscript𝒟𝑠𝑡𝑑\mathcal{D}_{std}caligraphic_D start_POSTSUBSCRIPT italic_s italic_t italic_d end_POSTSUBSCRIPT = 5 and vary the 𝒟meansubscript𝒟𝑚𝑒𝑎𝑛\mathcal{D}_{mean}caligraphic_D start_POSTSUBSCRIPT italic_m italic_e italic_a italic_n end_POSTSUBSCRIPT. As 𝒟meansubscript𝒟𝑚𝑒𝑎𝑛\mathcal{D}_{mean}caligraphic_D start_POSTSUBSCRIPT italic_m italic_e italic_a italic_n end_POSTSUBSCRIPT increases, we see a lesser number of strategic agents leaving the surge zone, causing the SA and NSB curves to approach each other as we move further to the right. Note that changing the parameters of 𝒟𝒟\mathcal{D}caligraphic_D has no impact on the convergence times τ𝜏\tauitalic_τ, as expected.

    [][𝒟mean=5subscript𝒟𝑚𝑒𝑎𝑛5\mathcal{D}_{mean}=5caligraphic_D start_POSTSUBSCRIPT italic_m italic_e italic_a italic_n end_POSTSUBSCRIPT = 5, τ=36𝜏36\tau=36italic_τ = 36]Refer to caption [][𝒟mean=10subscript𝒟𝑚𝑒𝑎𝑛10\mathcal{D}_{mean}=10caligraphic_D start_POSTSUBSCRIPT italic_m italic_e italic_a italic_n end_POSTSUBSCRIPT = 10, τ=37𝜏37\tau=37italic_τ = 37]Refer to caption [][𝒟mean=15subscript𝒟𝑚𝑒𝑎𝑛15\mathcal{D}_{mean}=15caligraphic_D start_POSTSUBSCRIPT italic_m italic_e italic_a italic_n end_POSTSUBSCRIPT = 15, τ=36𝜏36\tau=36italic_τ = 36]Refer to caption

    [][𝒟mean=5subscript𝒟𝑚𝑒𝑎𝑛5\mathcal{D}_{mean}=5caligraphic_D start_POSTSUBSCRIPT italic_m italic_e italic_a italic_n end_POSTSUBSCRIPT = 5]Refer to caption [][𝒟mean=10subscript𝒟𝑚𝑒𝑎𝑛10\mathcal{D}_{mean}=10caligraphic_D start_POSTSUBSCRIPT italic_m italic_e italic_a italic_n end_POSTSUBSCRIPT = 10]Refer to caption [][𝒟mean=15subscript𝒟𝑚𝑒𝑎𝑛15\mathcal{D}_{mean}=15caligraphic_D start_POSTSUBSCRIPT italic_m italic_e italic_a italic_n end_POSTSUBSCRIPT = 15]Refer to caption

    Figure 15: Price and demand curves for SA and NSB models for varying 𝒟meansubscript𝒟𝑚𝑒𝑎𝑛\mathcal{D}_{mean}caligraphic_D start_POSTSUBSCRIPT italic_m italic_e italic_a italic_n end_POSTSUBSCRIPT

5 Discussion & Future Work

In this work, we have explored an important aspect of strategic behavior of riders during price surges on ride-sharing platforms—riders monitor prices in real-time and may choose to walk away from the surge zone in search of cheaper rides. We have shown analytically how such strategic behavior has the potential to redistribute demand spatially and through extensive numerical experiments, demonstrated how it can equalize prices across the surge boundary faster. However, there are many interesting avenues of future work.

In recent years, we have seen an uptick in ride-sharing platforms related to riders deciding to walk when faced with high prices and long waiting times. Platforms have been trying to directly and explicitly incentivize riders to walk to and from pick-up and drop-off points by providing them with cheaper rides if they do so. Uber’s ExpressPool [Forbes, 2018] and Lyft’s Shared Saver [Lyft, 2019] are key examples of such a “walk and save” feature. In doing so, the platform can reduce driver downtime while improving operational efficiency. Over the years, these features have been modified and re-introduced multiple times which demonstrates an ongoing rider-side demand for such an option. This raises an important mechanism design question—how should platforms design rider incentives in an optimal way so that they can be encouraged to avail them while still protecting platform profits?

In addition to leaving the surge zone, another important aspect of strategic behavior on the rider side is waiting for prices to drop. Earlier works have modeled this effect in isolation; in the real world, however, strategic rider behavior is often a combination of both aspects: riders often initially wait at the surge zone for prices to drop, and then consider walking to other locations after too long of a wait. One immediate extension of our work is to model this more complex rider-side behavior with both waiting and walking effects simultaneously.

References

  • Agatz et al. [2012] Niels Agatz, Alan Erera, Martin Savelsbergh, and Xing Wang. Optimization for dynamic ride-sharing: A review. European Journal of Operational Research, 223(2):295–303, 2012.
  • Alban [2024] Jack Alban. ‘thank you b/c i be paying the $70 like a fool’: Traveler shares hack to getting uber for cheap from the airport. https://www.dailydot.com/news/uber-airport-hack/, 2024.
  • Alonso-Mora et al. [2017] Javier Alonso-Mora, Samitha Samaranayake, Alex Wallar, Emilio Frazzoli, and Daniela Rus. On-demand high-capacity ride-sharing via dynamic trip-vehicle assignment. Proceedings of the National Academy of Sciences, 114(3):462–467, 2017.
  • Ashkrof et al. [2020] Peyman Ashkrof, Gonçalo Homem de Almeida Correia, Oded Cats, and Bart Van Arem. Understanding ride-sourcing drivers’ behaviour and preferences: Insights from focus groups analysis. Research in Transportation Business & Management, 37:100516, 2020.
  • Banerjee et al. [2015] Siddhartha Banerjee, Carlos Riquelme, and Ramesh Johari. Pricing in ride-share platforms: A queueing-theoretic approach. Available at SSRN 2568258, 2015.
  • Besbes et al. [2021] Omar Besbes, Francisco Castro, and Ilan Lobel. Surge pricing and its spatial supply response. Management Science, 67(3):1350–1367, 2021.
  • Bimpikis et al. [2019] Kostas Bimpikis, Ozan Candogan, and Daniela Saban. Spatial pricing in ride-sharing networks. Operations Research, 67(3):744–769, 2019.
  • Braverman et al. [2019] Anton Braverman, Jim G Dai, Xin Liu, and Lei Ying. Empty-car routing in ridesharing systems. Operations Research, 67(5):1437–1452, 2019.
  • Cashore et al. [2022] J Massey Cashore, Peter I Frazier, and Eva Tardos. Dynamic pricing provides robust equilibria in stochastic ride-sharing networks. In Proceedings of the 23rd ACM Conference on Economics and Computation, pages 301–302, 2022.
  • Castillo [2023] Juan Camilo Castillo. Who benefits from surge pricing? Available at SSRN 3245533, 2023.
  • Castillo et al. [2017] Juan Camilo Castillo, Dan Knoepfle, and Glen Weyl. Surge pricing solves the wild goose chase. In Proceedings of the 2017 ACM Conference on Economics and Computation, pages 241–242, 2017.
  • Chan and Shaheen [2012] Nelson D Chan and Susan A Shaheen. Ridesharing in north america: Past, present, and future. Transport reviews, 32(1):93–112, 2012.
  • Chen and Cassandras [2020] Rui Chen and Christos G Cassandras. Optimization of ride sharing systems using event-driven receding horizon control. IFAC-PapersOnLine, 53(4):411–416, 2020.
  • Chen and Hu [2020] Yiwei Chen and Ming Hu. Pricing and matching with forward-looking buyers and sellers. Manufacturing & Service Operations Management, 22(4):717–734, 2020.
  • Fielbaum et al. [2021] Andres Fielbaum, Xiaoshan Bai, and Javier Alonso-Mora. On-demand ridesharing with optimized pick-up and drop-off walking locations. Transportation research part C: emerging technologies, 126:103061, 2021.
  • Forbes [2018] Forbes. Uber launches new cheaper rides, but there’s one big catch. https://www.forbes.com/sites/bizcarson/2018/02/21/uber-cheap-ride-expresspool-how-to-use/, 2018.
  • Furuhata et al. [2013] Masabumi Furuhata, Maged Dessouky, Fernando Ordóñez, Marc-Etienne Brunet, Xiaoqing Wang, and Sven Koenig. Ridesharing: The state-of-the-art and future directions. Transportation Research Part B: Methodological, 57:28–46, 2013.
  • Garg and Nazerzadeh [2022] Nikhil Garg and Hamid Nazerzadeh. Driver surge pricing. Management Science, 68(5):3219–3235, 2022.
  • Hausman and McFadden [1984] Jerry Hausman and Daniel McFadden. Specification tests for the multinomial logit model. Econometrica: Journal of the econometric society, pages 1219–1240, 1984.
  • Hu et al. [2022] Bin Hu, Ming Hu, and Han Zhu. Surge pricing and two-sided temporal responses in ride hailing. Manufacturing & Service Operations Management, 24(1):91–109, 2022.
  • Jin et al. [2018] Scarlett T Jin, Hui Kong, Rachel Wu, and Daniel Z Sui. Ridesourcing, the sharing economy, and the future of cities. Cities, 76:96–104, 2018.
  • Kingman [1993] John Kingman. Poisson processes. https://www.google.com/books/edition/Poisson_Processes/VEiM-OtwDHkC?hl=en&gbpv=0, 1993.
  • Kooti et al. [2017] Farshad Kooti, Mihajlo Grbovic, Luca Maria Aiello, Nemanja Djuric, Vladan Radosavljevic, and Kristina Lerman. Analyzing uber’s ride-sharing economy. In Proceedings of the 26th international conference on world wide web companion, pages 574–582, 2017.
  • Lien [2018] Tracey Lien. Honolulu lawmakers agree to cap surge pricing for uber and lyft. https://www.latimes.com/business/technology/la-fi-tn-uber-lyft-price-cap-20180607-story.html, 2018.
  • Lyft [2019] Lyft. Introducing shared saver, our most affordable ride. https://www.lyft.com/blog/posts/introducing-shared-saver, 2019.
  • Ma et al. [2022] Hongyao Ma, Fei Fang, and David C Parkes. Spatio-temporal pricing for ridesharing platforms. Operations Research, 70(2):1025–1041, 2022.
  • News [2015] ABC News. How uber chooses its surge price cap in emergencies. https://abcnews.go.com/Business/uber-chooses-surge-price-cap-emergencies/story?id=28494303, 2015.
  • Pandit et al. [2019] Vivek Nagraj Pandit, Datar Mandar, Manjesh K Hanawal, and Sharayu Moharir. Pricing in ride sharing platforms: static vs dynamic strategies. In 2019 11th International Conference on Communication Systems & Networks (COMSNETS), pages 208–215. IEEE, 2019.
  • Rheingans-Yoo et al. [2019] Duncan Rheingans-Yoo, Scott Duke Kominers, Hongyao Ma, and David C Parkes. Ridesharing with driver location preferences. arXiv preprint arXiv:1905.13191, 2019.
  • Selyukh [2015] Alina Selyukh. Uber surge price? research says walk a few blocks, wait a few minutes. https://news.wgcu.org/2015-10-29/uber-surge-price-research-says-walk-a-few-blocks-wait-a-few-minutes, 2015.
  • Tripathy et al. [2022] Manish Tripathy, Jiaru Bai, and H Sebastian Heese. Driver collusion in ride-hailing platforms. Decision Sciences, 2022.
  • Uber [2018] Uber. H3: Uber’s hexagonal hierarchical spatial index. https://www.uber.com/en-IN/blog/h3/, 2018.
  • Uber [2024] Uber. What’s happening when prices surge? https://www.uber.com/us/en/marketplace/pricing/surge-pricing/, 2024.
  • Yan et al. [2020] Chiwei Yan, Helin Zhu, Nikita Korolko, and Dawn Woodard. Dynamic pricing and matching in ride-hailing platforms. Naval Research Logistics (NRL), 67(8):705–724, 2020.
  • Yao et al. [2021] Danli Yao, Simai He, and Zhen Wang. A new ride-sharing model incorporating the passengers’ efforts. Naval Research Logistics (NRL), 68(4):397–411, 2021.
  • Yu et al. [2022] Chenkai Yu, Hongyao Ma, and Adam Wierman. Price cycles in ridesharing platforms. arXiv preprint arXiv:2202.07086, 2022.
  • Yu et al. [2023] Chenkai Yu, Hongyao Ma, and Adam Wierman. Price cycles in ridesharing platforms. In International Conference on Web and Internet Economics, pages 618–636. Springer, 2023.
  • Zhao et al. [2022] Daozhi Zhao, Ziwei Yuan, Mingyang Chen, and Shuang Yang. Differential pricing strategies of ride-sharing platforms: choosing customers or drivers? International Transactions in Operational Research, 29(2):1089–1131, 2022.

Appendix A Auxillary Results & Missing Proofs

A.1 Convergence

Theorem A.1.

The total unmet demand in the system, given by Ds(t)+Dns(t)subscript𝐷𝑠𝑡subscript𝐷𝑛𝑠𝑡D_{s}(t)+D_{ns}(t)italic_D start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT ( italic_t ) + italic_D start_POSTSUBSCRIPT italic_n italic_s end_POSTSUBSCRIPT ( italic_t ), must converge to 00 in finite time.

Proof.

Ds(t)0subscript𝐷𝑠𝑡0D_{s}(t)\geq 0italic_D start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT ( italic_t ) ≥ 0 and Dns(t)0subscript𝐷𝑛𝑠𝑡0D_{ns}(t)\geq 0italic_D start_POSTSUBSCRIPT italic_n italic_s end_POSTSUBSCRIPT ( italic_t ) ≥ 0 for all t𝑡titalic_t. Therefore, Ds(t)+Dns(t)>0subscript𝐷𝑠𝑡subscript𝐷𝑛𝑠𝑡0D_{s}(t)+D_{ns}(t)>0italic_D start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT ( italic_t ) + italic_D start_POSTSUBSCRIPT italic_n italic_s end_POSTSUBSCRIPT ( italic_t ) > 0 if and only if at least one of Ds(t)subscript𝐷𝑠𝑡D_{s}(t)italic_D start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT ( italic_t ) and Dns(t)subscript𝐷𝑛𝑠𝑡D_{ns}(t)italic_D start_POSTSUBSCRIPT italic_n italic_s end_POSTSUBSCRIPT ( italic_t ) is strictly greater than 00. Therefore, in each time period t𝑡titalic_t, at least a μλ>0𝜇𝜆0\mu-\lambda>0italic_μ - italic_λ > 0 amount (this amount is in fact 2(μλ)>μλ>02𝜇𝜆𝜇𝜆02(\mu-\lambda)>\mu-\lambda>02 ( italic_μ - italic_λ ) > italic_μ - italic_λ > 0 if both demands are positive) of unmet demand can be met every time step. This implies that Ds(t)+Dns(t)subscript𝐷𝑠𝑡subscript𝐷𝑛𝑠𝑡D_{s}(t)+D_{ns}(t)italic_D start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT ( italic_t ) + italic_D start_POSTSUBSCRIPT italic_n italic_s end_POSTSUBSCRIPT ( italic_t ) must reach 00 in at most D0+d0μλsubscript𝐷0subscript𝑑0𝜇𝜆\left\lceil\frac{D_{0}+d_{0}}{\mu-\lambda}\right\rceil⌈ divide start_ARG italic_D start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT + italic_d start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_ARG start_ARG italic_μ - italic_λ end_ARG ⌉ time steps. ∎

Corollary A.2.

Each of Ds(t)subscript𝐷𝑠𝑡D_{s}(t)italic_D start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT ( italic_t ) and Dns(t)subscript𝐷𝑛𝑠𝑡D_{ns}(t)ita