key: cord-0952665-40itos99 authors: Nagurney, Anna title: Supply chain game theory network modeling under labor constraints: Applications to the Covid-19 pandemic date: 2021-01-03 journal: Eur J Oper Res DOI: 10.1016/j.ejor.2020.12.054 sha: 48c5522a5fd8f619ebcf8b4ce32c1622abd3218c doc_id: 952665 cord_uid: 40itos99 The Covid-19 pandemic has brought attention to supply chain networks due to disruptions for many reasons, including that of labor shortages as a consequences of illnesses, death, risk mitigation, as well as travel restrictions. Many sectors of the economy from food to healthcare have been competing for workers, as a consequence. In this paper, we construct a supply chain game theory network framework that captures labor constraints under three different scenarios. The appropriate equilibrium constructs are defined, along with their variational inequality formulations. Computed solutions to numerical examples inspired by shortages of migrant labor to harvest fresh produce; specifically, blueberries, in the United States, reveal the impacts of a spectrum of disruptions to labor on the product flows and the profits of the firms in the supply chain network economy. This research adds to the literature in both economics and operations research. In this paper, we take up the challenge of constructing a supply chain network game theory framework for the modeling of competition among firms that produce a differentiated but substitutable product. The novelty of the framework lies in that it explicitly includes labor availability. Economists have emphasized the use of labor, as well as capital, as essential components of production functions (see Mishra, 2007 and the references therein). The inclusion of labor into a competitive supply chain network for differentiated products, however, has not been addressed. This is an important area of research since only through a full supply chain network perspective can one identify the impacts of labor availability and possible capacity disruptions in the pandemic on profits, costs, and consumer prices. Three different scenarios of labor availability are considered in our modeling framework, with accompanying constraints. In the first scenario, there is a bound on labor availability associated with each link of the supply chain network of each firm in the supply chain network economy. This is the most restrictive scenario in that labor capacities are imposed on individual links and, hence, there is not the freedom of movement from firm to firm and across tiers as in the other scenarios. In the second scenario, there is a bound on the labor availability associated with an activity tier of the supply chain networks, such as production, transportation, storage, and distribution. This scenario is relevant to the farming sector, since farmers compete for seasonal migrant workers for harvesting of the products, with the Covid-19 creating shortfalls of such labor in many parts of the globe (see, e.g., Corbishley, 2020 ) . In the third scenario, there is a single bound on labor availability in the supply chain network economy and labor is free to move across a tier or between tiers. Barrero et al. (2020) argue for the reallocation of labor because of the economic impacts of the Covid-19 crisis, which this scenario enables the evaluation and quantification of. The governing equilibrium conditions for the first scenario correspond to a Nash (1950 Nash ( , 1951 Equilibrium, whereas those for the second and third scenarios correspond to a Generalized Nash Equilibrium (GNE) (cf. Debreu, 1952 and Arrow & Debreu, 1954 ) , since the feasible sets associated with the firms' strategies are common, that is, shared, because of the respective labor constraints. Although game theory supply chain network models are now fairly well-established (cf. Nagurney, 2006; Nagurney, Alvarez Flores, & Soylu, 2016; Nagurney, Dong, & Zhang, 2002; Qiang, Ke, Anderson, & Dong, 2013; Saberi, 2018; Saberi, Cruz, Sarkis, & Nagurney, 2018; Toyasaki, Daniele, & Wakolbinger, 2014; Yu, Cruz, & Li, 2019 , and the references therein) and advances continue to be made (see, e.g., Choi, 2020 ) , Generalized Nash Equilibrium models have only recently been applied to supply chains. For example, Nagurney, Yu, and Besik (2017) introduced a competitive supply chain network equilibrium model with outsourcing in which firms competed for limited capacity at shared distribution facilities. However, labor was not explicitly considered. Also related to our theme here, in part, is the use of GNE supply chain models in disaster relief (cf. Nagurney et al., 2016 Nagurney et al., , 2018 Nagurney & Dutta, 2019; Nagurney, 2020a; Nagurney, 2020b ) and in healthcare ( Nagurney & Dutta, 2019 ) . We emphasize that both Nash Equilibrium and Generalized Nash Equilibrium are applied in this paper, for the first time, for supply chain network competition with labor. Nagurney (2020b) introduced labor into a supply chain network optimization model for perishable food in the Covid-19 pandemic. Therein, arc multipliers captured perishability, akin to the work of Yu and Nagurney (2013) , but only labor bounds on links were considered. Nagurney (2020c) considered a nonperishable product and also provided a system-optimization perspective for the supply chain network in the case of labor under different scenarios and in the case of elastic or fixed demands for the firm's product. In contrast to the framework in our paper, there was no competition among firms and no game theory concepts were utilized. The Covid-19 pandemic was declared only several months ago and already, over this relatively short period of time, has completely disrupted both lives and the economy. In terms of additional related literature, which is nascent, but, nevertheless, growing, we note that Queiroz, Ivanov, Dolgui, and Wamba (2020) detailed a research agenda through a structured literature review of Covid-19 related work and supply chain research on earlier epidemics. Ivanov (2020a) , in turn, described simulation-based research concentrating on the potential impacts on global supply chains of the Covid-19 pandemic. Ivanov and Das (2020) argued for the necessity of a new perspective due to the coronavirus Covid-19 outbreak through the use of intertwined supply networks (ISNs). An ISN is an entirety of interconnected supply chains which, in their integrity, yield the provision of society and markets with goods and services. Currie et al. (2020) in their paper identified multiple, complex challenges due to the COVID-19 pandemic and elaborated on how simulation modelling could help to support informed decision-making. Ivanov (2020b) , in turn, introduced a novel notion -that of a viable supply chain (VSC), in which viability is considered as an underlying supply chain property spanning three perspectives of agility, resilience, and sustainability. The contributions in the paper can assist firms in guiding their decisions on the recovery and re-building of their supply chains after crises of long duration such as the COVID-19 pandemic. Ivanov and Das (2020) captured the ripple effect of an epidemic outbreak in global supply chains in their model, with consideration of the velocity of pandemic propagation, the duration of production, distribution and market disruption, and a demand decline. In addition, the authors analyzed pandemic supply risk mitigation measures and possible recovery paths, accompanied by a relevant discussion of prospective global supply chain (re)-designs. The paper is organized as follows. In Section 2 , we construct the supply chain network game theory framework for differentiated products under three distinct labor scenarios. We identify the supply chain network structure, describe the behavior of the profitmaximizing firms, and the underlying constraints. For each scenario, we present two alternative variational inequality formulations of the governing equilibrium conditions. For Scenarios 2 and 3, which are Generalized Nash Equilibrium problems, we use the concept of a Variational Equilibrium to enable the variational inequality formulations. We also discuss existence results. The alternative variational inequality formulations allow for the implementation of an effective algorithm, which we present in Section 3 , along with computes solutions to a series of numerical examples. The numerical examples are motivated by shortages of migrant labor to harvest fresh produce, in particular, blueberries in the United States. The results of the paper are summarized in Section 4 , along with our conclusions, and suggestions for future research. There are I firms in the supply chain network economy that produce a substitutable product and compete noncooperatively in the production, transportation, storage, and distribution of their products to demand markets. The firms also compete with one another for labor, since labor is essential to the above network economic activities. Each firm is represented as a network of its economic activities as drawn in Fig. 1 . Observe that, according to Fig. 1 , the networks of the individual firms do not have any links in common. Table 1 contains the basic notation for the model. All vectors are column vectors. JID: EOR [m5G; 14:17 ] Fig. 1. The supply chain network topology of the model with labor. Notation for the supply chain game theory modeling framework with labor. The set of paths in firm i 's supply chain network terminating in demand market k ; i = 1 , . . . , I; k = 1 , . . . , n R . The set of all n P i paths of firm i ; i = 1 , . . . , I. P The set of all n P paths in the supply chain network economy. x p ; p ∈ P i k The nonnegative flow on path p originating at firm node i and terminating at demand market k ; i = 1 , . . . , I; k = 1 , . . . , n R . We group firm i 's product path flows into the vector x i ∈ R n P i + . We emphasize that x i is the vector of strategic variables of firm i . We then group all the firms' product path flows into the vector x ∈ R n P + . f a The nonnegative flow of the product on link a, ∀ a ∈ L . We group all the link flows into the vector f ∈ R n L + . l a The labor on link a (usually denoted in person hours). Positive factor relating input of labor to output of product flow on link a, ∀ a ∈ L . l a The upper bound on the availability of labor on link a under Scenario 1, ∀ a ∈ L l t The upper bound on labor availability for tier t activities under Scenario 2, with tier t = 1 being production; tier t = 2 refers to transportation, and, so on, until t = T, which corresponds to distribution. Here, T + 1 corresponds to the electronic commerce tier. l The upper bound on labor availability under Scenario 3. d ik The demand for the product of firm i at demand market k ; i = 1 , . . . , I; k = 1 , . . . , n R . We group the { d ik } elements for firm i into the vector d i ∈ R n R + and all the demands into the vector The total operational cost associated with link a, ∀ a ∈ L . πa Cost of a unit of labor on link a ρ ik (d) The demand price function for the product of firm i at demand market k ; i = 1 , . . . , I; k = 1 , . . . , n R . Each firm i ; i = 1 , . . . , I, owns n i M production facilities; can make use of n i D distribution centers, and can provide its product to the n R demand markets. Let L i denote the links comprising the supply chain network of firm i ; i = 1 , . . . , I, that it owns/controls, with a total of n L i elements. The links of L i include firm i 's links to its production nodes; the links from production nodes to the distribution centers, the storage links, and the links from the distribution centers to the demand markets. L then denotes the full set of links in the supply chain network economy with L = ∪ I i =1 L i with a total of n L elements. Let now G = [ N, L ] denote the graph consisting of the set of nodes N and the set of links L in Fig. 1 . Each firm seeks to determine its optimal product quantities that maximize its profits by using Fig. 1 as a schematic, coupled with the labor volumes, which we discuss the constraints on, further below. Specifically, production links from the top-tiered nodes i ; i = 1 , . . . , I, representing firm i, in Fig. 1 are connected to the production nodes of firm i, which are denoted, respectively, by: [m5G; 14:17 ] shipments to the demand markets. For example, in the case of demand markets corresponding to homes, such links can capture electronic commerce. They also can denote direct deliveries to consumers at demand markets where the producers are farms and such distribution channels have been initiated because of the pandemic (see, for example, Shea, 2020 ) . Of course, the supply chain network topology in Fig. 1 can be adapted for the specific application under consideration, with the appropriate addition/deletion of nodes, links, and/or supply chain network tiers. The demand for each firm's product at each demand market must be satisfied by the product flows from the firm to that demand market. Hence, the following conservation of flow equations must hold for each firm i : i = 1 , . . . , I: (1) Moreover, the path flows must be nonnegative; that is, for each firm i ; i = 1 , . . . , I: The link flows of each firm i ; i = 1 , . . . , I, are related to the path flows by the expression: where δ ap = 1 , if link a is contained in path p, and 0, otherwise. According to (3) , the flow of a firm's product on a link is equal to the sum of that product's flows on paths that contain that link. We now discuss how labor is related to product flow. In particular, we assume that the product output on each link is a linear function of the labor input. This corresponds to what is known as a linear production function in economics. Hence, we have that The greater the value of α a , the more productive the labor on the link. The utility function of firm i, U i ; i = 1 , . . . , I, is the profit, given by the difference between its revenue and its total costs: The first expression after the equal sign in (5a) is the revenue of firm i . The second expression in (5a) is the total operational costs for the supply chain network L i of firm i and the third expression captures the total labor costs of firm i . The functions U i ; i = 1 , . . . , I, are assumed to be concave, with the demand price functions being monotone decreasing and continuously differentiable and the total link cost functions being convex and also continuously differentiable. The optimization problem of each firm i ; i = 1 , . . . , I, is, hence, for firm i as follows: subject to: (1) -(4) . We now consider three distinct scenarios in terms of the labor availability in the supply chain network economy and the associated bound(s) on specific links. In Scenario 1, the additional constraints on the fundamental model described above are: l a ≤l a , ∀ a ∈ L. (6) According to (6) , there is an upper bound on labor associated with each link in the supply chain network of Fig. 1 . In this scenario, unlike the subsequent two scenarios, the feasible sets of the individual firms will depend only on their specific strategies and not on the strategies of the other firms, as we will show below. In Scenario 2, we consider that there is a bound on labor associated with each activity tier in the supply chain network economy, that is, in addition to the original constraints (1) through (4) , the firms are now faced with the following constraints: (7 , 2) and so on, until Observe that, in Scenario 2, unlike in Scenario 1, the firms now have shared, that is, common labor constraints. Hence, their underlying feasible sets will no longer be disjoint. This will result in a Generalized Nash Equilibrium, rather than a Nash Equilibrium, and we will elaborate further below when we present the variational inequality formulations of the model under the different scenarios. Hence, in this scenario, labor is movable across a tier within a firm or across firms. Scenario 2 allows for the movement of labor across the firms' different production facilities; across the distinct storage facilities, etc. Since one can expect that the skill set of a worker will be transferable across production facilities of a differentiated, but substitutable product; similarly, across distribution facilities or freight provision services, this scenario is quite reasonable in the pandemic. Scenario 3 may be interpreted as being the least restrictive of the scenarios considered here in that labor can be transferable across different activities of production, transportation, storage, and distribution. In the pandemic, we are seeing that some employees are assuming different tasks in supply chain networks than they had been doing previously. This may enhance agility and flexibility, provided that labor has the requisite skills or the skills can be acquired fairly readily. For example, farmers now in the US have been innovating in terms of direct sales to consumers and taking on different tasks, such as deliveries of their farmed products to homes (see Woolever, 2020 ) . This can benefit both producers and consumers. In Scenario 3, in addition to constraints (1) through (4) , the firms are now faced with the following single constraint: We now reformulate the objective function of each firm i ; i = 1 , . . . , I, given by (5b) in path flow variables exclusively. We are able to do this because of expressions (1), (2) , and (4) , which, recall, relates labor to product flow. Specifically, we can redefine the total operational cost link functions as: ˜ c a (x ) ≡ˆ c a ( f ) , ∀ a ∈ L, and the demand price functions as: ˜ In addition, using a result of Nagurney (2020b) , we know that, in view of (3) and (4) : l a = p∈ P x p δap αa , for all a ∈ L . Recall also that, according to Table 1 , x i denotes the vector of strategies, which are the path flows, for each firm i ; i = 1 , . . . , I. We can redefine the utility/profit functions ˜ JID: EOR [m5G; 14:17 ] and group the profits of all the firms into an I-dimensional vector ˜ U , such that Objective function (5b) , in lieu of the above, can now be expressed as: Furthermore, it readily follows that constraint (6) for Scenario 1 can be reexpressed exclusively in path flows; the same holds for constraints (7, 1) through (7,T+1) for Scenario 2, and for constraint (8) for Scenario 3. We now state the governing equilibrium conditions for the different scenarios and provide alternative variational inequality formulations for each scenario. We define the feasible set K i for firm i thus: In Scenario 1, each firm competes noncooperatively until the following equilibrium is achieved. A path flow pattern x * ∈ K is a supply chain network Nash Equilibrium if for each firm i ; i = 1 , . . . , I: where ˆ According to (11) , a supply chain Nash Equilibrium is established if no firm can improve upon its profits unilaterally. We know that K is a convex set. Applying the classical theory of Nash equilibria and variational inequalities, under our imposed assumptions on the underlying functions, it follows that (cf. Gabay & Moulin, 1980; Nagurney, 1999 ) the solution to the above Nash Equilibrium problem (see Nash, 1950; Nash, 1951 ) coincides with the solution of the variational inequality problem: determine x * ∈ K, such that where ·, · represents the inner product in the corresponding Euclidean space, which here is of dimension n P , and ∇ Existence of a solution to variational inequality (12) is guaranteed since the feasible set K is compact and the utility functions are continuously differentiable, under our imposed assumptions (cf. Kinderlehrer and Stampacchia, 1980 ) . We now provide an alternative variational inequality to (12) over a simpler feasible set. We introduce Lagrange multipliers λ a associated with the constraint (6) for each link a ∈ L and we group the Lagrange multipliers for each firm i 's network L i into the vector λ i . We then group all such vectors for the firms into the vector λ ∈ R n L + . Also, we define the feasible sets: Then, using similar arguments as in Theorem 1 in Nagurney et al. (2017) , the following result is immediate. The supply chain network Nash Equilibrium satisfying Definition 1 is equivalent to the solution of the variational inequality: determine the vector of equilibrium path flows and the vector of optimal Lagrange multipliers, (x * , λ * ) ∈ K 1 , such that: The above feasible set K 1 is the nonnegative orthant. This feature enables the implementation of an algorithm, which we describe in the next section, which is an iterative procedure that yields closed form expressions at an iteration for the path flows and the link Lagrange multipliers. In Scenarios 2 and 3, on the other hand, not only do the utility functions of the firms depend on their own strategies and those of the other firms, but the feasible sets do, as well. This happens because of the corresponding labor constraints, which capture that labor is "shared." Therefore, in Scenarios 2 and 3, the governing concept is no longer that of Nash Equilibrium, but, rather, it is that of Generalized Nash Equilibrium (see Debreu, 1952 and Arrow & Debreu, 1954 ) . This presents additional challenges, since, as noted in Nagurney et al. (2017) , Generalized Nash Equilibrium problems cannot be directly formulated as variational inequality problems, but, instead, quasi-variational inequalities are often used as the formulation (see, e.g., Facchinei & Kanzow, 2010 ) . It is well-known (cf. Luna, 2013 and the references therein) that quasi-variational inequality problems are, nevertheless, much harder to solve than finite-dimensional variational inequality problems. We can utilize a refinement of the Generalized Nash Equilibrium (GNE) in order to secure a variational inequality formulation. The refinement is a Variational Equilibrium and it is a specific type of GNE (see Kulkarni & Shanbhag, 2012 ) . In a Generalized Nash Equilibrium defined by a Variational Equilibrium, the Lagrange multipliers associated with the common/shared constraints are all the same. As noted in Nagurney et al. (2017) , this provides a fairness interpretation and is reasonable from an economic standpoint. We denote the set of shared constraints under Scenario 2 by We have the following definitions: A product flow vector x * ∈ K ∩ S 1 is a Generalized Nash Equilibrium under Scenario 2, if for each firm i ; i = 1 , . . . , I: where, as defined previously, ˆ x i * ≡ (x 1 * , . . . , x i −1 * , x i +1 * , . . . , x I * ) . A strategy vector x * is said to be a variational equilibrium of the above Generalized Nash Equilibrium according to Definition 2 if x * ∈ K ∩ S 1 is a solution of the variational inequality: Clearly, a solution x * to variational inequality (16) exists. For Scenario 2, we now associate the nonnegative Lagrange multiplier μ t with labor constraint (7 , t ) , for t = 1 , . . . , T + 1 . We define the vector of Lagrange multipliers μ ∈ R T +1 + and the feasi- } . An alternative variational inequality formulation to that of (16) , whose solution corresponds to the Variational Equilibrium under Scenario 2 is: determine (x * , μ * ) ∈ K 2 , such that: The variational inequality (17) will also allow for the application of an effective and efficient computational procedure. Recall that in Scenario 3 there is a single labor constraint, and labor is allowed to move freely among the supply chain network economic activities across the firms. Proceeding in a similar manner as for Scenario 2, we define S 2 , where S 2 ≡ { x | (8) holds } . We then can state the following definitions: A product flow vector x * ∈ K ∩ S 2 is a Generalized Nash Equilibrium Under Scenario 3, if for each firm i ; i = 1 , . . . , I: Of course a solution to variational inequality (19) is guaranteed to exist since the underlying feasible set K ∩ S 2 is compact. As for Scenarios 1 and 2, we now provide an alternative variational inequality formulation. We denote the Lagrange multiplier associated with the labor constraint (8) by γ . Also, we define the feasible set K 3 ≡ { (x, γ ) | (x, γ ) ∈ R n P +1 + } . We then have the following result. An alternative variational inequality to that of (16) is: Observe that the feasible set K 3 is also the nonnegative orthant and this will allow effective solution of variational inequality (20) . Also, it is important to emphasize that the Lagrange multipliers associated with the respective constraint sets in Scenarios 1, 2, and 3, provide valuable information, upon problem solution. In particular, they quantify the value of an additional unit of the labor resource associated with the respective constraint. Having a computational procedure that can compute both the equilibrium product path flows (from which the labor values are then obtained from: l a = p∈ P x p δap αa , for all links a ∈ L ) and the equilibrium Lagrange multipliers enriches the toolbox for managerial insights and decision-making. We present the algorithm in Section 3 . All of the above six variational inequalities can be put into standard variational inequality form, where (cf. Nagurney, 1999 ) the finite-dimensional variational inequality problem VI (F , K) , is to de- where F is a given continuous function from K to R N , K is a given closed, convex set, and ·, · denotes the inner product in Ndimensional Euclidean space. We, specifically, recommend and, then, apply the modified projection method of Korpelevich (1977) to compute solutions to the supply chain network game theory model under labor different scenarios. This algorithm was recently also applied to solve problem of human migration with policy interventions in Nagurney and Daniele (2021) . For the alternative variational inequality formulations of Scenarios 1, 2, and 3, given, respectively, by (13), (17) , and (20) , each iteration of the modified projection method yields closed form expressions for the path flows and for the associated Lagrange multipliers. The modified projection method is guaranteed to converge if the function F (X ) that enters the variational inequality (cf. (21) ) is Lipschitz continuous and monotone. These are reasonable conditions for the supply chain network game theory model with labor under the scenarios considered here. The definitions of monotonicity and Lipschitz continuity of F (X ) are as follow. The function F (X ) is said to be monotone, if and the function F (X ) is Lipschitz continuous, if there exists a constant L > 0 , known as the Lipschitz constant, such that The steps of the modified projection method are given below, with τ denoting an iteration counter: Step 0: Initialization Initialize with X 0 ∈ K. Set the iteration counter τ := 1 and let β be a scalar such that 0 < β ≤ 1 L , where L is the Lipschitz constant. Step 1: Computation Compute X τ by solving the variational inequality subproblem: Step 2: Adaptation JID: EOR [m5G; 14:17 ] Compute X τ by solving the variational inequality subproblem: Step 3: Convergence verification If | X τ − X τ −1 | ≤ , with > 0 , a pre-specified tolerance, then stop; otherwise, set τ := τ + 1 and go to Step 1. We now present the numerical examples. The modified projection method was implemented in FORTRAN and a Linux system at the University of Massachusetts Amherst used for the computation of solutions to the subsequent numerical examples. The algorithm was initialized as follows. We initialized the elastic demand for each demand market at 40 and equally distributed the demand among the paths connecting each demand market from each origin node (firm). The Lagrange multipliers were initialized to 0. The modified projection method was deemed to have converged if the absolute difference of the path flows differed by no more than 10 −7 and the same for the Lagrange multipliers. Our numerical examples are inspired by the summer harvest season for fresh produce and, specifically, the need to harvest berries, including blueberries. Much of the seasonal fresh produce harvesting in the United States is done by migrant workers. This is also the case in other parts of the world, including Europe. The harvesting of blueberries in the northeastern United States, with major growing areas being New Jersey and Maine, has been especially challenged this summer because of the Covid-19 pandemic (see Tully, 2020 and Russell, 2020 ) . Workers have become ill. Some have died, whereas others have not been able to secure permission to travel in time to harvest this perishable product. Interestingly, although operations researchers have been tackling human migration networks for several decades now, it is only recently that regulations have been incorporated in such network-based models (see, e.g., Daniele, 2021 and Cappello, 2021 ). Here we focus on seasonal migrants, rather than those who wish to relocate permanently to new locations. Our numerical examples are stylized but we did utilize Internet available resources to obtain blueberry price and picking data in the United States (see Galinato & Hong, 2016 and howmuchisit.org, 2018 ) . The flow variables are in pounds of blueberries, the prices are in dollars per pound, and labor is in person hours. In Section 3.1 , the numerical examples are of Scenario 1 and we solve the governing variational inequality (13) . In Section 3.2 , the examples are of Scenario 3, and we solve the governing variational inequality (20) . Examples 1, 2, and 3 have the supply chain network topology given in Fig. 2 . In these examples there is a single firm, with Example 1 serving as the baseline, which is a blueberry farm. It has two planting sites (production locations), access to a single distribution center, and the blueberries are distributed to two demand markets. The total operational cost functions on the links are: The costs for labor are: π a = 10 , π b = 10 , π c = 15 , π d = 15 , π e = 20 , π f = 17 , π g = 18 . The modified projection method yields the following solution: The equilibrium path flows are: The equilibrium link labor values are: As expected, all the Lagrange multipliers are equal to 0: The demand price at demand market 1 is: 5.98 and at demand market 2: 7.93 and the computed respective demands are: 177.23 and: 365.16. The profit of the firm is: 1,684.47. Example 2 has the same data as Example 1, except that now we consider the following situation. Workers have been getting sick and are not as available to pick the blueberries on the first harvesting site. Specifically, now the bound on the labor on link a has been drastically reduced to: l a = 10 . The modified projection method yields the following equilibrium solution: The equilibrium path flow pattern is: The equilibrium link labor values are: JID: EOR [m5G; 14:17 ] l * a = 10 . 00 , l * b = 11 . 70 , l * c = 2 . 40 , l * d = 2 . 92 , l * e = 10 . 65 , l * f = 1 . 72 , l * g = 3 . 60 . The Lagrange multipliers are all equal to 0, except for: The demand price at demand market 1 is: 5.98 and at demand market 2: 7.93 with the computed respective demands being, respectively: 172.18 and: 360.22 . The profit of the firm is: 1,680.61. Since the labor amount is at the bound on link a, the associated Lagrange multiplier is positive. Due to the labor disruption because of the pandemic, the profit of the firm decreases. The demands also decrease. Example 3 is a variant of Example 2 and has the same data except that now, due to illnesses, there is no labor available on link a ; hence, l a = 0 . 00 . The modified projection method converges to the following equilibrium solution. The equilibrium path flow pattern is: The equilibrium link labor values are: The Lagrange multipliers are all equal to 0, except for: We can see from this example the value of increasing labor availability on link a, as revealed through the significantly higher value of the Lagrange multiplier λ * a . The profit of the firm now drops to 1466.78. The demand price at demand market 1 is: 5.99 and at demand market 2 it is: 7.93. The computed respective demands are: 139.38 and: 327.99. The demand for blueberries drops at each demand market. Clearly, the farmer should do everything possible to secure the health of the workers at his production/harvesting facilities, so that the blueberries can be harvested in a timely manner and so that profits do not suffer. In Examples 4 through 6, we evaluate the impact of a competitor, that is, there is now another farm growing blueberries in the general area of the first farm. The underlying supply chain network topology is now as in Fig. 3 . The new firm also has two production sites, a single distribution center, and serves the same two locations as demand markets as Firm 1 does. The data for the first firm are identical to the data in Example 2 except with a modification of the firm's demand price functions in order to capture competition with the new firm. In particular, the demand price functions of Firm 1 are now: The Lagrange multipliers are all equal to 0.00 except to l * a = 4 . 93 . The product prices at equilibrium are: The profit for Firm 1 is: 1671.86 and the profit for Firm 2 is: 1134.61. The profit for firm 1 rises ever so slightly, whereas that for Firm 2 decreases. Example 6 has the same data as Example 5 except that we now consider a major disruption in terms of the spread of Covid-19 at the distribution centers of both firms with the bounds on labor corresponding to the associated respective links now being: l e = 5 , l l = 5 . The modified projection method computes the following equilibrium path flow pattern: The equilibrium link labor values are: All computed equilibrium Lagrange multipliers are now equal to 0 except for those associated with the distribution center link labor bounds which are: The product prices at equilibrium are now: The profit for Firm 1 is now drastically reduced to: 1218.74 and the profit for Firm 2 also declines but by a much smaller amount to 1126.73. In this subsection, we compute solutions via the modified projection method to numerical examples of Scenario 3. Specifically, the supply chain network topology remains as in Fig. 3 . The data for Examples 7 through 11 in this set are identical to that in Example 6, except there are no longer any link bounds but, rather, only a single labor bound l . Example 7 has l = 80 ; Example 8 has l = 70 ; Example 9 has l = 60 ; Example 10 has l = 30 , and Example 11 has l = 10 . The computed equilibrium path flows for these examples are reported in Table 2 with the equilibrium labor values given in Table 3 . The computed associated equilibrium Lagrange multiplier is as follows. For Example 7, the labor bound is not tight and, therefore, γ * = 0 . 00 . For Example 8, γ * = 3 . 99 ; for Example 9: γ * = 23 . 79 ; for Example 10: γ * = 43 . 59 and for Example 11: γ * = 68 . 00 . As for the profits, it is interesting to see how the disruptions in labor availability impact the profits negatively. In Example 7, Firm 1 earns a profit of 1,675.58, whereas Firm 2 earns a profit of 1,134.42. In Example 8, with labor being at the capacity, Firm 1 earns a profit of 1,671.23, whereas Firm 2 earns a profit of 1,131.76, which is only a small decrease. However, there is a big decrease in profit for both firms in Example 9. In Example 9, Firm 1 earns a profit of 1506.94, whereas Firm 2 earns a profit of 1022.56. In Example 10, with the labor capacity further reduced, Firm 1 earns only a profit of 1105.06, whereas Firm 2 earns a profit of only 753.56. And, in our final example, Example 11, with the labor bound l = 10 , Firm1's profit drops further to 523.19, whereas Firm2's profit drops to 228.91. Again, we emphasize that, although the above examples are stylized, and, clearly, examples of Scenario 2 can also be effectively computed using the proposed algorithm, the supply chain network game theory framework provides decision-makers with the flexibility of evaluating impacts of labor reductions due to Covid-19 and, of course, also investments in terms maintaining labor bound values at high levels. The latter can be achieved, in part, through social distancing, proper sanitation, and providing uncrowded living facilities for migrant workers. Furthermore, our numerical examples were inspired by crises faced by farmers due to labor shortages as a consequence of the Covid-19 pandemic. The supply chain network game theory frame-work can be adapted to different industrial sectors where labor plays an important role, such as, for example, in medical supply production from PPEs to ventilators and, in the future, if all goes well, also for medical treatments for Covid-19 patients and vaccines for everyone. The computed solutions demonstrate that the proposed algorithm is effective for the game theory supply chain network model with labor constraints and that the model is illuminating. Furthermore, having such an algorithm allows for the calculation of the quantitative impacts on equilibrium product flows, product prices, and firm profits of labor availability and possible disruptions under different scenarios. Finally, the computed equilibrium Lagrange multipliers associated with the various labor constraints provide the value in relaxing a specific labor constraint. In this paper, we introduced a supply chain network game theory framework consisting of multiple profit-maximizing firms competing noncooperatively in the production, storage, and ultimate distribution of their differentiated products in the presence of labor constraints. Specifically, we considered three different scenarios of labor constraints consisting of bounds on labor availability on supply chain links in the first scenario; bounds on labor availability across a tier of the supply chain networks of the first in the second scenario, and, finally, in the third scenario, a single bound on the labor in the supply chain network economy. The framework is a contribution to the operations research and to the economics literatures and is inspired by the challenges associated with labor and accompanying shortages in the Covid-19 pandemic. The governing concept in the case of labor Scenario 1 is that of Nash Equilibrium, whereas, because of the shared constraints in terms of labor, the governing concept in the case of labor Scenarios 2 and 3 is that of a Generalized Nash Equilibrium. For each labor scenario, we provide two variational inequality formulations, with the variational inequality that includes the Lagrange multipliers associated with the labor constraints being especially amenable to solution via our proposed algorithm, since both the equilibrium product flows and the equilibrium Lagrange multipliers can be computed iteratively using explicit formulae. Qualitative results are presented in terms of existence of the equilibrium patterns, along with solutions to a series of numerical examples. The numerical examples are grounded in disruptions taking place now in terms of shortages of migrant workers to harvest fresh produce, in particular, berries, and, specifically, blueberries in parts of the United States. Such labor shortfalls are not only a critical issue now in the United States but also in other parts of the globe, as a consequence of the Covid-19 pandemic. We emphasize that the framework constructed in his paper is a partial equilibrium one since the focus was on oligopolistic competition in supply chain networks of relevance to the Covid-19 pandemic. In future research, we expect that the supply chain network game theory model(s) herein can be adapted to specific industrial sectors and the data fitted empirically. We also expect that different labor constraints may be identified and perhaps even different functions than the linear ones that we utilized here to relate labor to product flow. Also, interesting avenues for future research are the use of cooperative game theory as well as the incorporation of costs associated with training/education to enable the reallocation of labor with appropriate skills where they may be shortfalls. This is happening, for example, in practice (cf. Reuters, 2020 ) , where we are seeing former airline workers being retrained as healthcare workers. Finally, it would be interesting to construct a general equilibrium model. JID: EOR [m5G Existence of an equilibrium for a competitive economy Coronavirus pandemic leads to idaho potato market woes COVID-19 and labour reallocation: Evidence from the US. voxEU.org US factories in China are open, but with 'severe' worker shortage Global shortage of personal protective equipment Farm-labor crisis under Covid-19 sends countries scrambling. Wolf Street The food chain weakest link: Slaughterhouses. The New York Times How simulation modelling can help reduce the impact of COVID-19 Fix it': Amazon workers demand protections as COVID-19 cases grow in Southern California facilities . California: Desert Sun A social equilibrium existence theorem Available at Generalized Nash equilibrium problems On the uniqueness and stability of Nash equilibria in noncooperative games 2015 cost estimates of establishing and producing conventional high bush bluberries in western Washington Why disinfectant wipes aren Dual sourcing under supply disruption with risk-averse suppliers in the sharing economy Competitive pricing of substitute products under supply disruption Will the pandemic accelerate automation in supply chains? Part 1: Labour shortages in a time of high demand . Reuters Coronavirus triggers acute farm labour shortages in europe Predicting the impacts of epidemic outbreaks on global supply chains: A simulation-based analysis on the coronavirus outbreak (COVID-19/SARS-cov-2) case Viable supply chain model: integrating agility, resilience and sustainability perspectives -Lessons from and thinking beyond the COVID-19 pandemic Coronavirus (COVID-19/SARS-cov-2) and supply chain resilience: A research note COVID-19 dashboard by the center for systems science and engineering (CSSE) at Johns Hopkins University (JHU) An introduction to variational inequalities and their applications Without federal protections, farm workers risk coronavirus infection to harvest crops The extragradient method for finding saddle points and other problems On the variational equilibrium as a refinement of the generalized Nash equilibrium COVID-19 risks to global food security Decomposition and approximation methods for variational inequalities, with applications to deterministic and stochastic energy markets A brief history of production functions Network economics: A variational inequality approach, second and revised edition Supply chain network economics: Dynamics of prices, flows, and profits A variational equilibrium network framework for humanitarian organizations in disaster relief: Effective product delivery under competition for financial funds How coronavirus is upsetting the blood supply chain . The Conversation Perishable food supply chain networks with labor in the Covid-19 pandemic Optimization of supply chain networks with inclusion of labor: Applications to Covid-19 pandemic disruptions A generalized Nash equilibrium model for post-disaster humanitarian relief Human migration networks and policy interventions: Bringing population distributions in line with system-optimization A supply chain network equilibrium model Competition for blood donations Supply chain network capacity competition with outsourcing: A variational equilibrium framework Equilibrium points in n-person games Noncooperative games Available at: COVID-19 and the food and agriculture sector: Issues and policy responses WHO director-general's opening remarks at the media briefing on COVID-19 -11 As COVID-19 cases rise, truckers share their stories of diagnosis and recovery. Land line Piglets aborted, chickens gassed as pandemic slams meat sector The closed-loop supply chain network with competition, distribution channel investment, and uncertainties Impacts of epidemic outbreaks on supply chains: Mapping a research agenda amid the COVID-19 pandemic through a structured literature review Critical supply shortages -The need for ventilators and personal protective equipment during the Covid-19 pandemic Laid-off SAS airline staff offered fast-track healthcare training COVID-19 cases among blueberry workers worsen farm labor shortage JID: EOR [m5G Sustainable, multiperiod supply chain network model with freight carrier through reduction in pollution stock A competitive multiperiod supply chain network model with freight carriers and green technology investment option From dirt to doorstep, during coronavirus crisis massachusetts farms find ways to come to you A variational inequality formulation of equilibrium models for end-of-life products with nonlinear constraints How you get your berries: Migrant workers who fear virus, but toil on . The New York Times Local farms embrace change in the face of coronavirus . Baltimore Magazine The sustainable supply chain network competition with environmental tax policies Competitive food supply chain networks with application to fresh produce The author is very grateful to Jose Figueroa and Stephen Cumberbatch, the excellent systems administrators in Engineering Computer Services at the University of Massachusetts Amherst, who, in the pandemic, assisted her and kept the computer server that she utilizes operational. She also acknowledges the remarkable essential workers -from healthcare workers to farmers and grocery store workers as well as freight service providers who enabled our functionality during this very challenging period in history.The author thanks the three anonymous reviewers plus the Editor for helpful comments and suggestions on two earlier versions of this paper. Also, the total operational costs associated with Firm2's supply chain network L 2 are:The costs for labor for Firm 2 are: π h = 11 , π i = 22 , π j = 15 , π k = 15 , π l = 18 , π m = 18 , π n = 18 .The link labor productivity factors for Firm 2 on its supply chain network are: α h = 23 , α i = 24 , α j = 100 , α k = 100 , α l = 70 , α m = 100 , α n = 100 .The bounds on labor, in turn, are:The four new paths associated with Firm 2 are:The modified projection method yields the following equilibrium solution: the equilibrium path flow pattern is:The equilibrium link labor values are:The Lagrange multipliers are all equal to 0.00 except for l * a = 4 . 9295 .The product prices at equilibrium are: The profit for Firm 1 is: 1671.80 and the profit for Firm 2 is: 1145.06. With the addition of a competitor, the prices of Firm 1 for its blueberries drop at the demand markets and its profit is also negatively impacted. Example 5 has the same data as Example 4 except that we modified the demand price functions ρ 21 (d) and ρ 22 (d) to include a cross term, so that: The new computed equilibrium path flow and Lagrange multiplier patterns are given below, as well as the equilibrium labor values on the links. The equilibrium path flow is now: l * a = 10 . 00 , l * b = 11 . 65 , l * c = 2 . 40 , l * d = 2 . 91 , l * e = 10 . 62 , l * f = 1 . 72 , l * g = 3 . 59 ,