1. Introduction and Literature Review
Emissions Trading System (ETS) is a mechanism used by authorities to incentivize the reduction of CO
2-equivalent (CO
2-eq) emissions by companies. The European ETS (EU ETS) was the first large ETS established [
1,
2] and helped to reduce emissions effectively from 2008. Different specific mechanisms can be exploited in an ETS to reach emissions reductions in a given time horizon. While a carbon tax sets the tariff for emitting but does not control the total emission, in a cap-and-trade scenario an upper bound on emissions is set up, and companies can exchange allowances between themselves in a market which contributes to defining the trading price. This mechanism allows the system to progressively reach emission reduction targets. This mechanism can be exploited at an industry level or by companies to implement emissions reduction and sustainability-oriented strategies. However, the setting up of the initial cap is not an easy task. Often it is not directly understandable how the marginal CO
2 price and the cap are related to each other and how they affect the level of investment of companies in green projects. To investigate these relations, in this paper, we describe a mathematical model where the cap is not imposed, and we find the condition under which the emissions are lowered without the need for regulation.
Previous work in the literature addressed cap-and-trade-related problems and applications. Ref. [
3] investigates relations among the lot size policy, the carbon price, and the delivered service level in a 2-echelon supply chain. In [
4], the customer is assumed to prefer green products and the producer has the option to invest in green technology. The problem is modeled as a Stackelberg game on a 2-echelon supply chain. Ref. [
5] considers cap-and-trade in a production and routing problem when companies can sacrifice customer demand when emissions costs increase. The authors account CO
2 by considering truck emissions in the routing process. Ref. [
6] studies whether or not to invest in green technology is determined by balancing the investment cost and its resulting benefit of reduced emissions. They consider a decreasing cost of emissions over time and decreasing emissions with increasing green technology. The authors use the model in realistic instances taken from a company manufacturing printed circuit boards, and the model is solved analytically. Ref. [
7] analyzes the impact of cap-and-trade policies on channel selection for a recycling optimization problem. The problem is formulated as a Mixed Integer Linear Program (MILP) and it considers the maximization of Net Present Value (NPV) and life cycle elements. The case study used for the model analysis is related to an Indian industry for steel drums. The authors find how the carbon emission reduction strategy has implications for optimal recycling channel selection. Ref. [
8] develop a cooperative game model with revenue sharing with cap-and-trade and analyze how carbon prices positively affect emissions reduction. As noted in recent studies [
9], customer interest in supply chain sustainability is important and represents an important pressure to companies for reducing carbon footprint and affordable cost.
Several methodological approaches have been proposed to integrate cap-and-trade in decision support models for strategic and tactical planning of operations. However, decisions involved in a cap-and-trade system are often cross-functional in an organization and affect several or all the companies belonging to an industry. So, it is important to model decisions taken by different actors and/or taken at different time scales. Most of the research work in CO
2 trading scheme considers, as an application field, energy generation. In particular, in [
10], a trading scheme is considered in a unit commitment problem and the resulting problem is solved with ad hoc heuristics. Cap and trade schemes are considered in microgrid-based energy systems in [
11]. The authors propose a MINLP model solved with different non-linear solvers and analyze the impact of carbon emission taxing and cap-and-trade systems on equipment selections. Location routing problems and carbon trading policies are analyzed by [
12]. They focus on the location of green facilities and use a Lagrangian decomposition technique to solve the problem. If we consider the emissions along the supply chain, another interesting contribution is the analysis of the impact of carbon trading schemes on supplier selection strategy. This problem is studied by [
13] with Analytic Hierarchy Process (AHP) and Analytic Network Process (ANP) methods. Moreover, the authors formulate an integer program to optimal allocate orders to suppliers, constrained to the emission trading scheme.
Ref. [
14] study the carbon tax and cap and trade policies using a bilevel programming model. They solve it by using the KKT conditions and assess how the optimal design of the tax system is related to social welfare, inducing cement producers to shift towards more environmentally friendly combustibles. Ref. [
15] provides a framework for evaluating the effects of carbon credit policies on the renewable energy market.
Recently, several studies focused on how to optimize incentives to lower emissions in a given time frame. The optimization of a carbon tax is studied by [
16], where the maximization of GDP is considered while carbon intensity is constrained. Carbon tax and tax allocation are studied by [
17] in a non-linear model where aggregate values are considered for all country industries.
If we consider cap-and-trade optimization, several papers address specific industries. In [
18], a Stackelberg game is used to verify revenue allocation policies for discrete production when a warranty period must be guaranteed and a cap-and-trade mechanism is established. Coordination mechanisms under a cap-and-trade scheme are studied in [
19], by analytical models and Pareto analysis in the cold chain industry. The author finds that cost and revenue-sharing mechanisms may be superior to lower cooperative schemes. Cap-and-trade in the re-manufacturing industry is studied in [
20]. In this case, the authors use a news vendor model to evaluate the profit under the cap-and-trade mechanism.
Table 1 draws a picture of the literature on the topic discussed above. Concerning the solution methods, there are several analytical approaches. However, for strategic and planning/operational decisions involving eco - sustainability, it is important to rely on complex methods, such as branch and price or metaheuristics. As considered by [
21], hybrid meta-heuristics are overtaking pure meta-heuristics for solving sustainability problems. The relevant literature is tabled considering the problem addressed, the CO
2 allocation principle, the objective, and the solution method. The table provides evidence of the research gaps related to the lack of emission account policies where the cap or the budget is not fixed at each period.
In this paper, further pursuing the goal of analyzing the effectiveness of cap-and-trade strategies in reducing emissions in supply chains, a mathematical model for sustainable supply chain management is presented with budgeting and location/allocation decisions, where cap-and-trade is studied in an unregulated scenario to minimize emissions and cost. The model is a bi-objective optimization program that poses the goal of minimizing emissions and supply chain costs in an unregulated scenario w.r.t the cap definition; this means that while trading CO2 is allowed among facilities, no formal limit on the cap of CO2 emissions is imposed. The model encompasses an initial budget for technological investments by the facilities in the considered supply chain which allows plants to reduce their unit production emissions at a different unit production cost.
Theoretical conditions guaranteeing that, if obeyed, the emissions over time have a non-increasing trend are identified. This means that decreasing caps over time can be attained with a self-regulated scenario whether such conditions are attended. Experiments are conducted to analyze what if the model is run without respecting these theoretical conditions and what if one includes additional constraints in the model to take into account the former conditions.
The remainder of the paper is as follows.
Section 2 describes the problem and shows the proposed mathematical formulation.
Section 3 analyzes the experimental behavior of the proposed model.
Section 4 presents some theoretical conditions allowing a correct behavior of the model in producing a non-increasing trend of the emissions over time then guaranteeing the capability of automating the system emission caps over time. Finally,
Section 5 lists some final remarks and poses the basis for future work on this topic.
2. Problem Definition and Mathematical Formulations
A set F of production facilities, each with capacity , with , operates in a single market area with known demand D. There is a planning horizon T discretized into periods indexed in . Given are also a set S of suppliers, each with capacity , with . Let be the graph representing this supply chain, where and . In the following, we will use the terms plant and facility interchangeably. Moving a unit of commodity from to generates a transportation cost and CO2 emissions , at any time .
Each facility may use different technology levels, from 1 to (), to carry out the production over the time horizon, even though only one technology level must be adopted in each period. Moreover, technology levels adopted by each plant must be non-decreasing over time. Based on the level of technology l adopted in plant j, emissions per unit of product in that plant vary; in particular, the higher l, the lower in j. At the same time, the unit production cost in plant j varies with the technology level l adopted; in particular, the higher l, the higher in j.
Each facility , in each time period t can result in the following:
Trade emissions, receiving carbon quotas from or providing carbon quotas to other facilities;
An increase in the technology level owned at to a level in order to reduce emissions for a unit or worked product; this implies an investment by facility j. Trivially, if the technology level of j at time t remains unchanged w.r.t. we have .
To define the mathematical model, we use the following sets:
- T:
the set of time periods indexed by t;
- F:
the set of facilities indexed by j and ;
- S:
the set of suppliers indexed by s;
- L:
the set of technology levels indexed by l and ;
The model parameters are as follows:
- :
CO2 emissions per unit of product manufactured in plant with technology level [ton/unit]; with ;
- :
CO2 emissions per unit of product transported from supplier to plant [ton/unit];
- :
unit production cost in facility using technology level ;
- :
unit transportation cost from supplier to facility ;
- D:
market demand within the time horizon;
- B:
the overall budget that facilities can use to invest in (green) technologies;
: installation cost of technology level in plant assuming a negligible level of technology in j;
- :
capacity of a facility ;
- :
capacity of a supplier ;
- :
weighting value for the emissions term in the objective function.
The decision variables are:
: holds 1 if plant works with a technological level at time and 0 otherwise;
: amount of production in plant at time ;
: amount of products moved from supplier to plant at time ;
: amount of CO2 exchanged from plant to plant at time ;
: supply chain global emissions at time ;
: net CO2 emissions of facility at time (takes into account also carbon quotas traded).
The mathematical model, denoted as
UM (Unregulated Model), is reported below:
The objective function is a convex combination () of two measures, i.e., the overall supply emissions , and the supply chain costs given by the sum of the production and the transportation costs.
Constraint (2) defines the overall emissions at time t, given by the sum of the facility emissions and the transport emissions. Constraint (3) defines the emissions due to facilities: for each facility j at time t, the emissions is given by the sum of the following:
The term which, based on the technological level l chosen for plant j at time t, calculates the product between the unit emission cost associated with technology level l, i.e., , times the flow of products manufactured in j at time t, and
The total CO2 exchanged among plant j and the other facilities.
Constraint (4) says that the sum of manufactured products by all the facilities over the time horizon must be equal to the demand D. Constraint (5) imposes that the technological level of a plant at a certain time t must not be higher than that at a time . Constraint (6) warrants that each facility is associated with exactly one technological level in L at each period. Constraint (7) is the budget constraint, and Constraint (8) calculates the sum of flows entering a plant from the suppliers. Constraints (9) and (10) are capacity constraints while Constraints (11)–(16) define the signs of the variables.
Linearization
Both objective function (1) and Constraint (3) contain non-linear terms. However, they can be linearized by introducing additional variables
,
,
, and
, and replacing the term
so as
. Accordingly, the following constraints must be added:
It can be easily checked that Constraints (17)–(20) connect
to
as detailed in the following:
3. Experimental Analysis of the Model
In this section, we show how the model behaves in trading CO2 among facilities and keeping the emissions under control over time. The goal is to see whether the proposed model produces a non-increasing trend of CO2 emissions over time producing the same trend as the one obtainable by a regulated cap-and-trade model. The model has been coded using the Python library and solved with the solver GUROBITM release 9.5. The machine used for the experiments is equipped with a 12th Gen Intel(R) Core(TM) i7-1260P @ 2.10 GHz processor with 32 GB RAM.
We defined two instances, named
and
, respectively, with different numbers of suppliers and facilities: the former instance has three suppliers and three facilities, and the latter instance has ten suppliers and ten facilities. The complete set of parameters of the two instances are detailed in
Table 2 and
Table 3, respectively, where “baseline values”
bv means that parameters are determined as
bv plus a value chosen uniformly at random between
in
and
in
. In addition, for each baseline instance, different instances are generated by changing
from 0 to 1 with step 0.1. This generation is useful also for drawing the Pareto curve showing the variation of the total emissions compared to the total cost.
In
Figure 1 and
Figure 2, we depict the emissions
over time of the model
UM for instances
and
, respectively. For both instances, it can be easily checked how the total emission does not decrease monotonically over time. For the same instances, in
Table 4 and
Table 5 we report the CO
2 amount sent out from facility
j to facility
, i.e.,
,
.
Also, the emission balance for each facility over time has been computed and reported in
Table 6 for instance
and in
Table 7 for instance
. Here, for each facility
j, we report the CO
2 amount
exchanged in each period with the other facilities: it is positive (negative) if facility
j at time
t sent out more (less) CO
2 to the other facilities than that received from the other facilities. Trivially, zero means that the sum of the CO
2 quotas flowing out
j at time
t equals the overall amount of CO
2 flowing in
j at the same time. Thus, the tables show how much the trading mechanism is used by the facilities to respect the cap constraints. However, being the model unregulated, trading alone is not sufficient to let the total emission decrease. As can be observed by the charts, the trend of
over time is not monotone, which means that the model, while minimizing emissions, cannot guarantee a decreasing cap in an unregulated scenario. In the next section, we prove some conditions in this regard.