1 Introduction

The promotion of environmental sustainability has become crucial in the implementation and day-to-day functioning of complex and global supply networks. A growing number of large multi-national enterprises are implementing more stringent environmental practices; this has also an impact on small and medium-sized enterprises (SMEs) which are often involved in their supply networks (UN Global 2011).

The promotion of tighter environmental standards also implies the adoption of benchmarking approaches for comparing the sustainability performances of supply chains against industry standards. Within this context, Life Cycle Assessment (LCA) methodologies allow the estimation of cumulative environmental impacts across the entire supply chain and against a wide range of indicators, adopting a full product life cycle perspective. Within the current context of growing environmental sustainability, performing an allocation of environmental impacts to the different actors in the supply chain, with the aim of understanding which of them should be deemed responsible for polluting activities and related taxes or abatement costs, is a relevant and timely issue, if proper mitigation measures need to be implemented.

The interest in the topic of environmental pollution responsibility started due to the need to establish mechanisms to perform pollution burden sharing actions across countries (Zhou and Wang 2016); the related literature has then evolved in order to address allocation problems at different levels. According to an analysis of the current literature (see, for instance, Zhou and Wang 2016), most of the available approaches such problems at a macro-level; the application of pollution responsibilities approaches to contemporary multi-tier and multi-stakeholders supply chains is often overlooked.

This paper aims at addressing this research gap, by providing a normative framework based on cooperative game theory aimed at studying the problem of pollution responsibility allocation across multi-tier supply chains.

The paper is organised as follows. After a review of the literature concerned with pollution responsibility allocation, and the identification of relevant gaps (Sect. 2), the paper presents a cooperative game theoretical responsibility model for pollution allocation (Sect. 3). In Sect. 4, the paper details the generic model in the case of a linear supply chain. In Sect. 5, a practical application of the introduced cost allocation rules is provided; finally, some conclusions are drawn.

2 Literature review

The problem of associating pollution responsibilities with economic actors traces its roots in the international trade literature (Leontief and Ford 1970; Proops et al. 1993; Wyckoff and Roop 1994; Imura and Moriguchi 1995; Subak 1995; Bosquet 2000; Korhonen 2002; Bastianoni et al. 2004; Lenzen and Murray 2010). While the interest in the topic started because of the need to establish mechanisms to perform pollution burden sharing actions across countries, the related literature has then evolved in order to address allocation problems at different levels. These include the decomposition of national objectives into regional ones and the distribution of environmental pollution responsibility across partner firms in a given supply chain.

One of the fundamental research questions in the field of environmental pollution responsibility has been represented by the allocation principle to be followed. Scholars have developed different allocation criteria that can be mainly categorized into fairness- and efficiency-inspired principles. The concept of “fairness” is generally linked to equitable distribution issues (Rose 1990); the one of “efficiency” is more related to the minimization of costs related to abatement and mitigation measures. Many emissions allocation methods have been proposed. According to the extensive review performed by Zhou and Wang (2016), these can be classified into multiple categories. One of them is the so-called indicator approach. Through the development of specific indicators, pollution responsibilities, targets or permits are worked out in relation to these measures. Another stream of studies is devoted to the use of optimization approaches, which is based on mathematical programming frameworks for pollution allocation, which try to minimize costs for pollution abatements. For example, Ridgley (1996) combed the use of composite indicators and an optimization method for deriving suitable pollution responsibility allocations at a country level. This rapid scan of the extant literature, coherently with findings from Zhou and Wang (2016), highlights that most of the developed approaches deal with the pollution responsibility problem at a macro-level. Indeed, most of the methodologies tackle the problem at a country or regional level.

While some firm-level approaches can be retrieved, the supply chain perspective has been, so far, largely ignored, with the few available approaches being characterized by very simple and naïve models, mainly dealing with dyadic configurations; the application of pollution responsibilities approaches to contemporary multi-tier and multi-stakeholders supply chains is often overlooked.

Pollution allocation mechanisms often include interactions involving multiple actors. As such, pollution allocation mechanisms may be represented in a very effective way through approaches based on Game Theory. As stated by Zhou and Wang (2016), game theoretic approaches seem to be naturally suited for this purpose. Such allocations result from equilibria among different parties’ interactions or from values, which may satisfy enjoyable properties imposed by a central authority. Compared to the traditional approaches for pollution costs allocation, game theoretic approaches can be deemed as more sophisticated and less immediate. The use of game theoretic approaches is underexploited, especially when dealing with complex and multi-tier supply chains. Concepts from cooperative game theory may assume a normative value, i.e. a value that is imposed by an external authority.

Chander and Tulkens (1995) and Filar and Gaertner (1997) developed the first adaptation of Game Theory to the problem in order to investigate the allocation of GHGs emission reductions quota among countries. A similar approach, at an international level, was developed by Eyckmans and Tulkens (2003), while Germain and Steenberghe (2003) utilized a dynamic game framework. Viguier et al. (2006) deployed a two-level game for solving a similar allocation problem. At a firm level, MacKenzie et al. (2008, 2009) utilized, respectively, rank-order contests and incomplete information games for allocating pollution allowances. Chung et al. (2013) utilized dynamic games to evaluate the response of companies to environmental pollution taxes in a spatially distributed supply chain. Ren et al. (2015) proposed a Stackelberg game for studying CO2 reduction targets in a very simple dyadic supply chain.

Surprisingly, the development of normative approaches to pollution cost allocation is less extensive. A practical applications of the Shapley value approach was developed by Liao et al. (2015) for allocating emission allowances across energy producers in Shanghai.

2.1 Contribution of the paper

The next section presents a unified cooperative game theoretical approach to pollution responsibility allocation in a multi-tier supply chain. In Sect. 4, the paper details the generic model in the case of a linear supply chain, by selecting three different responsibility principles (namely: Upstream, Downstream and Local Responsibility). Also, fairness, efficiency and transparency properties are introduced, as desirable characteristics for allocation rules; furthermore, a stability concept for efficient allocations is formulated. In particular, here we present the Shapley value for this unified cooperative game theoretical responsibility model. In doing this, we find interesting connections between the so-called river problem introduced by Dong et al. (2012) in the pollution games literature and our framework. The Shapley allocation satisfies many of the mentioned desirable properties.

In Sect. 5, a practical application of the introduced cost allocation rules is provided, with the reference to the context of dairy supply chains. Finally, some conclusions are drawn.

3 Cooperative responsibility models for environmental supply chains

Consider a supply chain of a given product. Such a supply chain can be seen as a set N of companies engaged, at various stages, in the manufacturing of the product, through sourcing of raw materials, sub-components manufacturing, assembly activities; such companies can be indexed in a given order. A set of processes \( P_{i} \) (with \( |P_{i} | \le \left| N \right| - 1 \)) can be associated with each company \( i \). In this context, a process \( \left( {i, j} \right) \) represents the production of goods by company \( i \) and the supplying of such goods to company j (with i ≠ j) (as shown in Fig. 1).

Fig. 1
figure 1

Structure of a supply chain network

As such, Pi⊆ P, where \( P = \bigcup\nolimits_{i = 1 \ldots n} {P_{i} } \). The set \( P_{i} \) might be empty and this implies that company \( i \) is not supplying goods or services to any supply chain partner. Conversely, the set \( P_{i} \) might be coinciding with the whole set of processes (i.e., \( P \)); this means that company \( i \) is the only supplier in the considered network. Also, each process \( \left( {i,j} \right) \) has an environmental cost that is denoted as \( c_{ij} \). Let \( C \) be the set of all environmental costs. Therefore, we build a mathematical framework as follows:

$$ \left( {N, P, C} \right) $$
(1)

In order to enhance this framework with the concept of environmental responsibility, it can be useful to introduce a \( \left| N \right| \times \left| P \right| \) responsibility matrix \( B = \left( {B_{i, p} } \right) \) where the row index \( i \) represents companies and the column index \( p \) represents processes. Significantly, \( b_{ip} = 1 \) if company \( i \) is responsible (from an environmental point of view) for process \( p \), \( b_{ip} = 0 \) otherwise. Let \( B_{i} \) be the set of the processes for which company \( i \) is responsible from an environmental point of view, that is \( B_{i} = \left\{ {p \in P | b_{ip} = 1} \right\} \). It is important to recall that the environmental responsibility of a process may not depend on the physical location of the process; in other words, \( b_{ip} \) can be equal to 1 even if process does not involve physically company \( i \) (in other words, if the company \( i \) is not involved as a supplier in this process, or, alternatively, if the production of goods involved in the process does not happen at company \( i \) premises). Therefore, the framework (1) can be rewritten as \( \left( {N, P, B, C} \right) \), where \( B \) is the responsibility matrix. The quantity \( \sum\nolimits_{{p \in B_{i} }} {c_{p} } \) can be interpreted as an intuitive cost allocation for company \( i \), summing up all the environmental costs of processes \( p \) for which company \( i \) is deemed responsible. A coalition responsibility set can be defined as \( B_{S} = \cup_{i \in S} B_{i} \); then, the social cost function for each sub-group of companies \( S \subseteq N \) can be defined as follows:

$$ v\left( S \right) = \mathop \sum \limits_{{p \in B_{s} }}^{ } c_{p} $$
(2)

The sum \( v\left( S \right) \) in the formula (2) represents the environmental cost of all the processes for which at least one company belonging to the set S is responsible. Moreover, each cost is counted only once as different companies may be responsible for the same process.Footnote 1 Being \( N \), by definition, a finite set, \( \nu : 2^{N} \to R \) is a function which associates a real value with each subset of \( N \). Following a classical definition, \( \left( {N, \nu } \right) \) represents a cooperative game with a characteristic function. \( N \) also represents the set of players; each subset of the set of players is called a coalition; \( \nu \) is called the characteristic function of the game.

In this case, the characteristic function \( v \) is defined through the above-introduced elements \( P, B \) and \( C \), therefore leading to the following model:

$$ G = \left( {N,P,B,C,v} \right) $$
(3)

Such cooperative games are defined in terms of a characteristic function, which specifies the outcomes that each coalition can achieve for itself. Outcomes are specified in terms of the total utility/cost that a coalition can get and distribute across its members. By assuming the formation of a grand coalition, the main aim of such games is the definition of a solution concept, which allocates costs among each player from \( N \).

In a supply chain context, it can be fair assuming that companies cooperate (by forming a grand coalition) in order to coordinate the activities of the production system. This way, they will incur in different individual production costs because of initial binding agreements among partners; also, they will face costs related to the mitigation of the environmental impacts of their production activities. Cooperative game theory can provide solutions to allocate environmental pollution impacts across each member of the production system, by allowing transferable payments among companies.

In this study, the main aim is to identify a vector \( x = \left( {x_{i} } \right) \in R_{ + }^{\left| N \right|} \) representing a solution to the cooperative model (2), which assigns the cost \( x_{i} \) to each company \( i \).

An allocation is responsible-compatible when, if \( B_{i} \) is empty, then \( {\text{x}}_{\text{i}} \) is null. In the following, we will always refer to allocations which respect the responsible-compatible principle.

An allocation is efficient if the sum of all cost-allocations, i.e. \( \sum\nolimits_{{{\text{i}} = 1}}^{{\left| {\text{N}} \right|}} {x_{i} } \), is equal to the sum of all costs, i.e. \( \sum\nolimits_{i = 1}^{\left| N \right|} {\sum\nolimits_{j = 1}^{\left| N \right|} {c_{ij} } } \).Footnote 2

3.1 Desirable properties for cost allocation rules

Cost allocation rules should respect some desirable properties that have been described in the environmental pollution responsibility literature. In the following, we list and discuss each of these properties; we remark that some of these have been already listed in Gopalakrishnan et al. (2016), while and some of them are novel.

Property 1: Equivalence

If two companies \( i \) and \( j \) are responsible for the same processes, then they must have the same cost allocation. Therefore, if \( B_{i} = B_{j} \) then \( x_{i} = x_{j} \) must be satisfied for any \( i \) and \( j \in N. \)

Property 2: Equal sharing of extra pollution

This property provides fairness to the allocation of extra costs that might arise as the result of additional polluting activities (that could happen over time). If there is an increase in total pollution costs, it is required that companies responsible for this increase should be equally affected by extra burdens.

As a preliminary definition, it can be stated that C′ ≥ C if and only if cij′ ≥cijfor any\( i, j \in N \). We say that C′ > C if C′ ≥ C and there exists a process\( \left( {i, j} \right) \in P \)such that cij′ > cijLet\( G = \left( {N, P, B, C, v} \right) \), \( G^{\prime} = \left( {N, P,B,C^{\prime},v} \right) \)be the two associated cooperative models where\( {\text{C}} \)and C′ are two sets of environmental costs with C′ ≥ C. If\( i, j \in N \)are responsible for any process\( p \), that is\( b_{i, p} = b_{j,p} = 1 \)where\( p \in P \)is characterized by cp′ > cp. Then the following equality\( x_{i} \left( {G^{\prime}} \right) - x_{i} \left( G \right) = x_{j} \left( {G^{\prime}} \right) - x_{j} \left( G \right) \)must be satisfied.

Property 3: No free riding

The equal sharing of extra pollution property (introduced above as Property 2 ) does not guarantee the efficiency of such net cost increase re-allocation. In other words, the property does not require that the sum of extra charges for responsible companies is equal to the net increase of the pollution costs across the supply chain. As such, Property 3 provides an additional efficiency requirement to the above-mentioned situation. This property requires that if the total pollution costs increase, but, at the same time, for some firms the pollution of the processes they are responsible for remains unchanged, the allocation for these firms should remain the same.

Let \( G = \left( {N, P, B, C, v} \right) \) and \( G^{\prime} = \left( {N,P,B,C^{\prime},v} \right) \) be the two associated cooperative models where \( {\text{C}} \) , C′ are two sets of environmental costs with C′ ≥ C. Let \( i \in N \) ; suppose that \( c_{p}^{\prime } = c_{p} \) for each process \( p \) satisfying \( b_{i, p} = 1 \) . Then, the equality \( x_{i } \left( {G^{\prime}} \right) = x_{i} \left( G \right) \) must be satisfied.

Property 4: Process independence

Suppose that a single process \( p \) increases its pollution cost \( c_{p} \) ; as a consequence, cost allocations across the supply chain will change. In other words, cost allocations might be affected by the new cost distribution over the set of all processes even if just one single process increases its cost. If this does not happen, it means that the variation of cost allocations solely depends on the process whose cost has increased; then, we say that cost allocations are process independent. As a consequence, this property implies that the firms can transparently relate their cost variations to the processes with increased costs; as such, they may make an official grievance against firms involved in those more costly processes.

Let\( G = \Big( {N, P, B, C, v} \Big) \),\( G^{\prime} = \Big( {N, P, B,C^{\prime},v} \Big) \),\( \tilde{G} = \Big( {N,P,B,\tilde{C},v} \Big) \),\( \tilde{G}^{\prime} = \Big( N,P,B,\tilde{C^{\prime}},v \Big) \)be the cooperative models associated to the cost structures\( C,C^{\prime},\tilde{C},\tilde{C}^{\prime}. \)We assume that\( C^{\prime} \ge C \), \( \tilde{C}^{\prime} \ge \tilde{C} \)and\( c_{p}^{\prime } = c_{p} \), \( \tilde{c}_{p}^{\prime } = \tilde{c}_{p} \)for\( p \in M{ \setminus } \left\{ q \right\} \). We assume that\( \tilde{c}_{q}^{\prime } = c_{q}^{\prime } \)and\( c_{q} = {\tilde{\text{c}}}_{q} \). Then,\( x_{i} \left( {G^{\prime}} \right) - x_{i} \left( G \right) = x_{i} \left( {\tilde{G^{\prime}}} \right) - x_{i} \left( {\tilde{G}} \right) \)must hold for each\( i \in N \).

Property 5: (Weak) unilateral disaggregation stability

Lenzen et al. (2007) and Rodrigues and Domingos (2008) originally discussed the stability of pollution responsibility allocations due to de-merging of companies over supply chains. Let\( G = \left( {N, P, B, C, v} \right) \)be our cooperative model. As a preliminary step, we formalize disaggregated supply chains with reference to the original model.

Let \( P_{ } \left( i \right) \)\( \{ j \in N| \left( {i, j} \right) \in P\} \) and \( R_{ } \left( i \right) \)\( \{ j \in N| \left( {j, i} \right) \in P\} \) be the subsets including firms, who are supplied by \( i \) or supply firm \( i \), respectively. Similarly, let \( P_{i} \) and \( R_{i} \) be the subsets including processes in which firm \( i \) is a supplier or firm \( i \) is supplied by, respectively. Assume that firm \( i \) disaggregates its activities into two companies \( i_{1} \) and \( i_{2} \). After disaggregation, the set of firms is modified as \( N^{\prime} = \left\{ {i_{1} , i_{2} } \right\} \cup N\backslash \left\{ i \right\} \). Then, the new set of processes is composed by \( P^{\prime} = R\left( i \right) \times \left\{ {i_{1} } \right\} \cup \left\{ {i_{2} } \right\} \times P\left( i \right) \cup (i_{1} ,i_{2} ) \cup P\backslash (P_{i} \cup R_{i} ) \). The new set of environmental costs \( C^{\prime} \) is composed as it follows. First of all, we define \( c_{{ji_{1} }}^{\prime } : = c_{ji} \) for \( j \in R\left( i \right) \) and \( c_{{i_{2} j }}^{\prime } ,c_{{i_{1} i_{2} }}^{\prime } \ge 0 \) where \( c_{ij} : = c_{{i_{2} j}}^{\prime } + c_{{i_{1} i_{2} }}^{\prime } \) for each \( j \in P\left( i \right) \). The last equality is named cost-splitting property. For the remaining cases, we simply require that \( c_{hk }^{\prime } : = c_{hk} . \) It is straightforward to think that disaggregated responsibility situations need a new responsibility matrix, that is \( B^{\prime}. \) First of all, such a new responsibility matrix \( B^{\prime} \) must have |N| + 1 rows, that is the new number of firms, and |P| + 1 columns, that is the new number of processes. Such a new responsibility matrix might be defined in several ways. We require some reasonable assumption on it. Let \( h \) be a firm participating in the original game and in the disaggregated game. Then it is immediate to see that \( h \) cannot be firms \( i_{1} ,i_{2} ,i \). For such a firm, we require some consistency property related to the original responsibility matrix \( B \). Firstly, \( B_{hp} = 1 \) with \( p = \left( {{\text{i}}, {\text{j}}} \right) ^{ }_{ } \) implies \( B_{{{\text{h}}p^{\prime}}}^{\prime } = 1 \) with \( p^{\prime} = \left( {i_{2} ,{\text{j}}} \right) \); \( B_{hp} = 1 \) with \( p = \left( {{\text{j}}, {\text{i}}} \right) ^{ }_{ } \) implies \( B_{{{\text{h}}p^{\prime}}}^{\prime } = 1 \) with \( p^{\prime} = \left( {j,i_{1} } \right) \). Finally \( B_{{{\text{h}}p}}^{\prime } = B_{{{\text{h}}p}} \) holds if the process \( p \) does not contain \( i_{1} ,i_{2} \).

The above conditions ensure that, if a firm is responsible for processes involving the firm \( i \) in the original game, is also responsible for associated processes involving \( i_{1} ,i_{2} \) after the disaggregation. The latter means that the disaggregation does not create externalities on others’ responsibilities for processes. In addition, we require that \( \left( {i_{1} ,i_{2} } \right) \notin B^{\prime}_{{i_{1} }} \cap B^{\prime}_{{i_{2} }} \); this means that, at most, one of the new companies is responsible for the internal process \( \left( {i_{1} ,i_{2} } \right), \) due to the disaggregation. If \( B^{\prime} \) satisfies both conditions we say that a responsibility matrix \( B^{\prime} \) is \( B^{ } \)-compatible.

We have defined disaggregated models \( G^{\prime} = \left( {N^{\prime}, P^{\prime}, B^{\prime}, C^{\prime}, v} \right) \) from the original model \( {\text{G}} = \left( {N^{ } , P^{ } , B^{ } , C^{ } ,v} \right) \). The multiplicity of disaggregated models depends on the multiple ways to split costs and/or on the multiple ways to implement new pollution responsibilities over supply chains.

Within this context, we say that \( {\text{G}} = \left( {N^{ } , P^{ } , B^{ } , C^{ } , v} \right)^{ } \) satisfies the unilateral disaggregation stability property if there exists a \( B \)-compatible responsibility matrix \( B^{\prime} \) and, then, there exists a firm \( i \) such that \( x_{{i_{1} }} \left( {G^{\prime}} \right) + x_{{i_{2} }} \left( {G^{\prime}} \right) > x_{i} \left( G \right) \) for any disaggregated model \( G^{\prime} = \left( {N^{\prime}, P^{\prime}, B^{\prime}, C^{\prime}, v} \right) \) (Unilateral Disaggregation Stability). We say that firm \( i \) is the aggregating company and \( B^{\prime} \) is the aggregating responsibility matrix. We say that \( G^{ } \) satisfies weak unilateral disaggregation stability if \( x_{j} \left( G \right) \ge x_{j} \left( {G^{\prime}} \right) \) (Weak Unilateral Disaggregation Stability) for each company \( j \) different from i.

Assume that cost allocations satisfy unilateral disaggregation stability. In virtue of the cost-splitting property, it is straightforward to prove that:

$$ \mathop \sum \limits_{h = 1}^{n} \mathop \sum \limits_{k = 1}^{n} c_{hk} \ge \mathop \sum \limits_{h = 1}^{n + 1} \mathop \sum \limits_{k = 1}^{n + 1} c_{hk}^{\prime } . $$

Since cost allocations are efficient, we have

$$ \mathop \sum \limits_{h = 1}^{n} x_{h} \left( G \right) = \mathop \sum \limits_{h = 1}^{n} \mathop \sum \limits_{k = 1}^{n} c_{hk} = \mathop \sum \limits_{h = 1}^{n + 1} x_{h} \left( {G^{\prime}} \right) = \mathop \sum \limits_{h = 1}^{n + 1} \mathop \sum \limits_{k = 1}^{n + 1} c_{hk}^{\prime } . $$

We say that firm \( i \) is the aggregating company. It follows that

$$ \sum\limits_{h = 1}^{n} {x_{h} } \left( G \right) - x_{\text{i}} \left( G \right) > \sum\limits_{h = 1}^{n + 1} {x_{h} } \left( {G^{\prime } } \right) - x_{{i_{1} }} \left( {G^{\prime } } \right) - x_{{i_{2} }} \left( {G^{\prime } } \right) $$

The last formula implies that \( x_{h} \left( G \right) > x_{h} \left( {G^{\prime}} \right) \) for some companies \( h. \) This means that there exists, at least, one firm who might consider the disaggregation option as a profitable one. If each company can get a benefit from disaggregation, then the whole supply chain can be considered as weakly stable, even if the aggregating firm prefers not to disaggregate its activities. In addition, if the cost of original game is strictly superior to the cost of the disaggregated game, there are internal incentives along supply chains to cause disaggregation. Even if cost allocations satisfy the unilateral disaggregation stability property, internal resources might still provide incentives to companies, which do not directly demerge, to collude and cause the disaggregation of the aggregating firm.

3.2 Additional remarks

The first two properties (Equivalence and Equal Sharing of Extra Pollution) respect fairness by definition. It must be highlighted that, according to an egalitarianism principle, each company across a supply chain is granted the same right to pollute (Zhou and Wang 2016). As such, the total pollution across a supply chain might be increased due to unilateral actions from individual partners. However, in this case, an egalitarian perspective does not necessarily imply distributive justice concerning pollution responsibility allocation deriving from extra pollution (Rose 1990). Within this context, the equal sharing of extra pollution, no free riding and process independence properties provide additional commitment to substantive distributive justice.

4 Linear supply chains and pollution responsibility

Here, we specify our model \( \left( {N, P, B, C, v} \right) \) in (2) where \( \upsilon \) has been already defined in (1). The specifications deal with the structure of the process set \( P \) and of the responsibility matrices \( B \). A generic process \( \left( {i, j} \right) \) represents the production of goods by company \( i \) and the supply of such goods to company \( j \). In the linear supply chain, the set \( P_{i} \) is reduced to a singleton for each company, at most. This implies that company \( i \) is supplying goods or services to a single supply chain partner. Then a process \( \left( {i, j} \right) \) can be identified as an element of the set of companies. In addition, each process has an environmental cost \( c_{i,i + 1}^{ } \), which can be simply expressed by \( c_{i } \). In the context, we can assume that existence of a final company |N| + 1, which is the company such that \( P_{\left| N \right| + 1} = \emptyset . \) From company |N|, a process starts and supplies company |N| + 1. We assume that the final company has no responsibility, that is characterized by the following assumption \( B_{\left| N \right| + 1} = \emptyset . \) Because allocations are responsible compatible by our initial choice, we know that \( x_{\left| N \right| + 1} = 0 \). Then, we disregard the last company \( \left| N \right| + 1 _{ } \) because its cost allocation is null. Therefore the number of processes \( \left| P \right|_{ } \) is equal to the number of pertinent firms is \( \left| N \right| = n_{ } \).

Here, we define how the responsibility matrix \( B_{ } \) is characterized according to environmental supply chain philosophies. In this context, three principles can be employed to solve the allocation problem. A Local Responsibility principle (LR), according to which each company \( i \) is strictly responsible for the pollution costs, related to the production activities strictly happening at its premises. An Upstream Responsibility principle (UR), stating that upstream suppliers (dealing with raw material extraction, sub-component manufacturing and other energy intensive activities) are responsible not only for pollution happening at their premises, but can also influence the environmental performance of downstream partners. A Downstream Responsibility principle (DR), stating that downstream partners in the supply chain are responsible for the polluting activities happening at upstream suppliers’ premises.

The three above principles can be exemplified by introducing the following responsibility matrices \( B \). The following diagonal unitary \( n \times n \) matrix characterizes LR:

$$ B_{ } = \left( {\begin{array}{*{20}c} 1 & 0 & \ldots & 0 \\ 0 & 1 & \ldots & 0 \\ \ldots & \ldots & \ldots & \ldots \\ 0 & 0 & \ldots & 1 \\ \end{array} } \right) $$

The UR principle says that player \( j \) has the responsibility to participate in pollution mitigation expenses of each downstream supply chain partner. Then, the following unitary upper triangular \( n \times n \) responsibility matrix characterizes the UR principle.

$$ B_{ } = \left( {\begin{array}{*{20}c} 1 & 1 & \ldots & 1 \\ 0 & 1 & \ldots & 1 \\ \ldots & \ldots & \ldots & \ldots \\ 0 & 0 & \ldots & 1 \\ \end{array} } \right) $$

Similarly, the following unitary lower triangular \( n \times n \) responsibility matrix characterizes DR principle.

$$ B_{ } = \left( {\begin{array}{*{20}c} 1 & 0 & \ldots & 0 \\ 1 & 1 & \ldots & 0 \\ \ldots & \ldots & \ldots & \ldots \\ 1 & 1 & \ldots & 1 \\ \end{array} } \right) $$

After having defined the structure of process set P and the responsibility matrices \( B \), the characteristic function, which has been defined in (1), depends just on the set of environmental costs, that is \( C \). It will assume different forms for each of the three principles. By adopting the LR principle, the characteristic function (1) turns out to be:

$$ \nu_{1} \left( S \right) = \mathop \sum \limits_{i \in S} c_{i} . $$

Conversely, it will assume the following form if a UR principle is adopted:

$$ \nu_{2} \left( S \right) = \mathop \sum \limits_{{i = j_{min} }}^{n} c_{i} ,\quad j_{ \hbox{min} } = \hbox{min} S. $$

Finally, by adopting a DR principle, the characteristic function will take the following form:

$$ \nu_{3} \left( S \right) = \mathop \sum \limits_{i = 1}^{{j_{max} }} c_{i} ,\quad j_{max} = \hbox{max} S. $$

Therefore the model \( \left( {N, P, B, C, v} \right) \) in (2) can be simplified as \( \left( {N, C, \nu } \right) \), according to the shape of the above-mentioned characteristic functions. In the following section, an allocation rule for the environmental responsibility costs is illustrated. This allocation rule is efficient, in such a way to ensure that \( \sum\nolimits_{i = 1}^{n} {x_{i} } = \sum\nolimits_{i = 1}^{n} {c_{i} } \) among all the supply chain partners for all the three games introduced above \( \left( {N,C,\nu_{i} } \right). \)

5 The Shapley allocation for linear supply chains

Among the possible allocation rules that can be derived from cooperative game theory, the Shapley value enjoys desirable properties such as efficiency and fairness (Owen 1995). The game \( \left( {N, C,\nu_{1} } \right) \) admits a Shapley value \( {\text{x}}_{\text{i}} = c_{i} \). The formula can be explicitly expressed by the following equations:

$$ \begin{array}{*{20}l} {x_{1} = c_{1} } \hfill \\ {x_{2} = c_{2} } \hfill \\ {x_{3} = c_{3} } \hfill \\ { \cdots = \cdots } \hfill \\ {x_{n} = c_{n} = 0} \hfill \\ \end{array} $$

This particular allocation can be called Local Responsibility Sharing (LRS). For an arbitrary company \( i \), if \( c_{i} = 0 \) then \( x_{i} = \) 0. In this case, LRS allocation is null for firm \( i \), even if \( B_{i} \ne \emptyset . \) Mathematical results provided in Dong et al. (2012) are helpful to have a closed- formula for the Shapley value for the game \( \left( {N, C, \nu_{2} } \right) \). According to such results, \( \left( {N, C, \nu_{2} } \right) \) admits a Shapley value which is equivalent to \( x_{i} = \sum\nolimits_{s = i}^{n} { {\text{c}}_{\text{s}} } /s \). The formula can be explicitly expressed by the following equations:

$$ \begin{array}{*{20}l} {x_{1} = c_{1} + \frac{{c_{2} }}{2} + \frac{{c_{3} }}{3} + \cdots + \frac{{c_{n} }}{n} = c_{1} + \frac{{c_{2} }}{2} + \frac{{c_{3} }}{3} + \cdots + \frac{{c_{n - 1} }}{n - 1}} \hfill \\ {x_{2} = \frac{{c_{2} }}{2} + \frac{{c_{3} }}{3} + \cdots + \frac{{c_{n} }}{n} = \frac{{c_{2} }}{2} + \frac{{c_{3} }}{3} + \cdots + \frac{{c_{n - 1} }}{n - 1}} \hfill \\ {x_{3} = \frac{{c_{3} }}{3} + \cdots + \frac{{c_{n} }}{n} = \frac{{c_{3} }}{3} + \cdots + \frac{{c_{n - 1} }}{n - 1}} \hfill \\ { \cdots = \cdots } \hfill \\ {x_{n} = \frac{{c_{n} }}{n} = 0} \hfill \\ \end{array} $$

This particular allocation can be called Upstream Equal Sharing (UES) cost allocation. Similarly, thanks to the result provided in Ni and Wang (2007) and Gómez-Rúa (2013), the game \( \left( {N, C, \nu_{3} } \right) \) admits a Shapley value, i.e. \( {\text{x}}_{{\text{i}}} = \mathop \sum \nolimits_{{{\text{s}} = 1}}^{\text{i}} {{\text{c}}}_{{\text{s}}} /\left( {{\text{n}} - {{\text{s}}} + 1} \right) \). The formula can be explicitly expressed by the following equations:

$$ \begin{array}{*{20}l} {x_{1} = \frac{{c_{1} }}{n}} \hfill \\ {x_{2} = \frac{{c_{1} }}{n} + \frac{{c_{2} }}{n - 1}} \hfill \\ {x_{3} = \frac{{c_{1} }}{n} + \frac{{c_{2} }}{n - 1} + \frac{{c_{3} }}{n - 2}} \hfill \\ { \cdots = \cdots } \hfill \\ {x_{n} = \frac{{c_{1} }}{n} + \frac{{c_{2} }}{n - 1} + \frac{{c_{3} }}{n - 2} + \cdots + c_{n} = \frac{{c_{1} }}{n} + \frac{{c_{2} }}{n - 1} + \frac{{c_{3} }}{n - 2} + \cdots \frac{{c_{n - 1} }}{n - 1}} \hfill \\ \end{array} $$

This particular allocation can be called Downstream Equal Sharing (DES) allocation method. Here, we verify if LRS, DES and UES allocation rules satisfy the desirable properties introduced in Sect. 3.

Proposition 1

LRS, DES and UES allocation rules satisfy the equivalence property.

Proposition 2

LRS, DES and UES allocation rules satisfy the equal sharing of extra pollution property.

Proposition 3

LRS, DES and UES allocation rules satisfy the no free riding property.

Proposition 4

LRS, DES, UES cost allocations satisfy the process independence property.

Because the Shapley value is an efficient allocation, then UES cost allocations are efficient.

Proposition 5

LRS allocation rule does not satisfy the unilateral disaggregation stability property. We assume that \( c_{n - 1 } \ne 0 \) and \( n > 2. \) Then, UES allocation rule satisfies unilateral disaggregation stability property. UES allocation rule satisfies weak unilateral disaggregation stability.

6 Applications

In order to demonstrate the applicability of the Shapley allocation (and related allocation rules) to real-world problems, the developed approach has been tested on an example related to the dairy supply chain. Dairy supply chains enjoy a very linear structure, in which the following stakeholders can be identified (Vergé et al. 2013):

  • Feed producer, producing animal feed from raw agricultural products;

  • Milk Producer, managing lactating animals (generally cows) and obtaining raw milk;

  • Processor, performing important operations for converting milk into a product fit for human consumption (heat treatments and pasteurization) and transforming it into a final product to be sold onto the market (also taking care of the packaging element);

  • Distributor, taking care of the logistical processes throughout the supply chain, involving cold chain operations;

  • Retailer, making the product available to the final consumer.

  • Consumers, buying and consuming the final product.

In a typical dairy supply chain, environmental pollution can be characterized through carbon equivalent emissions (expressed in Kg CO2-eq per Litre). Combining findings from several sources (Vergé et al. 2013; Thoma et al. 2013; Ormond and Goodman 2015), environmental impacts (measured in terms of Kg CO2-eq per litre of product) happening at each stage of the supply chain can be reported as shown in Table 1. For sake of simplicity, pollution costs might be deemed proportional to these environmental impacts.

Table 1 Environmental impacts across a typical dairy supply chain

The results of the three allocation principles shown in Sects. 3 and 4 (LRS, UES, DES) to a typical dairy supply chain are shown in the following Table 2. It can be seen how the three allocation principles provide very different results, allocating different quota of environmental impacts to different actors. Under the LRS rule, the highest proportion of environmental impacts (and, therefore, of associated mitigation costs) is assigned to the Milk Producer, respecting a simple proportionality mechanism. UES and DES rules develop more complex mechanisms; notably, the DES rule strongly penalizes the Retailer and the Consumer.

Table 2 Environmental impact allocation according to LRS, UES, DES rules

It must be pointed out that the proposed allocation rules might be considered not mutually exclusive; indeed, convex combinations of these rules might be developed, as shown in Table 3, where LRS and DES are combined together. This could be done, for instance, in order to introduce, within a LRS framework, elements of downstream responsibility.

Table 3 Hybrid Allocation Rules

7 Conclusions and future works

This study shed light on normative frameworks for the problem of pollution responsibility allocation across multi-tier supply chains through cooperative game theory concepts. After a thorough literature review, the paper has presented a cooperative game theoretical responsibility model for pollution allocation; the paper has detailed the generic model in the case of a linear supply chain, by selecting three different responsibility principles (namely: Upstream, Downstream and Local Responsibility) and developing some associated pollution responsibility allocation rules. These rules have been tested against some desirable properties that have been proposed. Furthermore, a practical application of the introduced cost allocation rules has been provided in the context of dairy supply chains.

The work proposed in this paper can be further extended in future researches, in order to overcome some of the limitations that characterize the current approaches. First of all, the model has been implemented on a simple linear supply chain, in which each company supplies and receives goods from a single partner; future work could extend the framework to more complex supply networks, in which multiple companies operate at each tier. Also, different pollution responsibility schemes and different supply chains structures (including, for instance, reverse elements, such as theorized by Jacobs and Subramanian 2012) might be employed in our model. A new set of normative properties could be defined, with the specific aim to assess the stability of international supply chains, especially against pollution costs differently perceived and taxed in nonhomogeneous countries.