Next Article in Journal
Analysis of Psychological Factors Influencing Mathematical Achievement and Machine Learning Classification
Previous Article in Journal
Spatial Demo-Genetic Predator–Prey Model for Studying Natural Selection of Traits Enhancing Consumer Motility
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Joint Optimal Production Planning in Supply Chains under Diverse Low-Carbon Policies

1
College of Intelligent Manufacturing, Tianjin Sino-German University of Applied Sciences, Tianjin 300350, China
2
Department of Supply Chain Management and Information Systems, Rennes School of Business, 2 Rue Robert d’Arbrissel, 35065 Rennes, France
3
Industrial and Manufacturing Systems Engineering Department, Hong Kong University, Hong Kong
*
Authors to whom correspondence should be addressed.
Mathematics 2023, 11(15), 3379; https://doi.org/10.3390/math11153379
Submission received: 4 March 2023 / Revised: 21 June 2023 / Accepted: 28 June 2023 / Published: 2 August 2023

Abstract

:
This paper studies the impact of diverse carbon emission regulations on optimal operation decisions by exploiting the economic production lot model of a multi-manufacturer system in a low-carbon environment. The optimal production planning for multiple enterprises in terms of cost optimization and carbon emission optimization are presented in various scenarios of the fully decentralized system under carbon taxation, cap and trade, and mandatory carbon cap. We prove a unified framework for modeling the scenarios enabled by carbon tax and cap-and-trade policies. Furthermore, the decision vector of the optimal production lot is obtained by analyzing the joint decision scenario based on shared carbon emission caps. We give a comparative analysis of two settings of fully decentralized and carbon quota sharing. The results show that the decision of joint production and stock preparation can reduce the total operational costs of the multi-enterprise system, but the realized carbon emissions may not change or only increase slightly. Our study provides a useful reference for government carbon emission regulation and enterprise operation decision optimization under carbon neutrality and a carbon peak environment.

1. Introduction

In a low-carbon economy, firms can reduce carbon emissions by either adopting direct emission reduction technology in response to the pressure of carbon emissions or adjusting their operations strategies to meet various low-carbon constraints. In this paper, the typical economic production quantity (EPQ) model is used to examine production tactics under the constraint of diverse low-carbon policies, such as mandatory carbon cap, cap and trade, and carbon taxation. Although a couple of basic carbon emission regulations have been explored in the literature, these policies, including mandatory carbon emission policy, are among the most common and easy to operate. This study shows the feasibility of adjusting operation strategy to realize carbon emission control by taking the cost optimization problem under low-carbon regulations.
Supply chain multi-links like procurement, production, and inventory usually generate carbon emissions in different forms, and the trade-off among the links also has an important impact on carbon emissions control, so the joint decision-making among these functions has attracted much attention. For example, production quantity, equipment installation frequency, and inventory level are interrelated, and affect the cost and carbon emission level simultaneously. To solve these problems, the economic order quantity model under carbon emission constraints is used to discuss the optimal order or production quantity and emission control, which is also extended to the newsvendor environment. Furthermore, there is research on joint decision-making of different supply chain functions based on cost and emission reduction considerations. Similarly, research concerning the optimal production-inventory decision of closed-loop supply chains under carbon tax policy or considering the comprehensive decision of inventory path and operation adjustment emerges. Hence, researchers have refined from general goods to special goods inventory and extended from a simple secondary to a closed-loop supply chain. However, the issues related to carbon emission quota sharing under the mandatory carbon cap policy that the present study pays attention to have not been reflected yet in the existing literature.
Driven by practice and theory, researchers have realized the importance of integrating carbon emissions and sustainability into the supply chain from the relationship between supply chain management and environmental protection. From the angle of measuring the carbon footprint of the supply chain, carbon emission reduction investment often brings great cost pressure to enterprises and affects supply chain efficiency and sustainable development. Hence, there are studies on carbon emission reduction technology investment and carbon emission reduction initiative choice at the supply chain level, such as discussing the joint investment behavior of supply chain members in emission reduction technology. Although these researches pay attention to the interaction of carbon emission reduction strategies at the supply chain level, which provides new insights for inter-enterprise emission reduction, it seldom involves the influence of carbon emission quota sharing on inter-enterprise decision-making interaction.
Low carbon policy can also influence strategic supply chain decisions like supply chain network design. The existing studies usually relate carbon emissions to transportation links for investigating how to choose and optimize transportation modes to reduce carbon emissions as much as possible. Furthermore, research exists in response to the tactical level through operations integration, such as order quantity and inventory when considering carbon emission regulation and carbon trading, which points out the necessity and importance of considering carbon emissions in supply chain operations. Along this vein, this paper integrates carbon emissions into the traditional EPQ model to explore optimal production decisions under diverse scenarios.
The distinct effects of carbon policy heterogeneity on supply chain operations are also stressed in the literature. In view of the effect of emission reduction achieved in practice, examining the influence of different carbon regulations on a firm’s decision-making and social welfare not only helps enterprises to cope with various policy constraints flexibly but also inspires the government to formulate optimal policies. In this paper, we use nonlinear optimization theory to build a mathematical programming model to study the impact of various low-carbon policies on multi-enterprise operations. In this sense, our models show the power of simple models to describe low-carbon operation problems and the significance of describing low-carbon joint decision-making problems among enterprises in some typical scenarios. Particularly, the influence of different strategies of using carbon caps on the operations in this paper has not been covered previously.
The remainder of the paper is organized as follows. A literature review is given in Section 2. In Section 3 presents the multi-firm decentralized system under various regulative low-carbon policies including carbon tax, cap and trade and mandatory cap policies. Section 4 proceeds to discuss the joint economic production decisions under caps sharing. A complete comparison and numerical analysis of these two scenarios come in Section 5. Concluding remarks are in Section 6.

2. Literature Review

This study is related to several streams of literature on low-carbon environment, which will be reviewed as follows.
First, our work is relevant to the literature on considering factors of carbon emissions and sustainability in EOQ/inventory models. Ping He et al. found the optimal production and emissions of a firm under the regulations of total control and trading and carbon tax [1]. Vincent Hovelaque et al. proposed an EOQ model with carbon constraints. The model considered the link between inventory policy, carbon emissions, price, and environment-related demand to maximize the firm’s profits while reducing emissions [2]. Ali Bozorgi et al. proposed an inventory model applicable to items that must be stored at a specific temperature [3]. Taleizadeh et al. constructed four new sustainable EPQ models to address different inventory shortage situations in the production system [4]. Liao et al. incorporated demand uncertainty into the EOQ model and solved the multi-objective optimization problem by deducing the optimal acquisition strategies for remanufactured products under a cap-and-trade policy [5]. S. Wang et al. compared the EOQ and JIT inventory models considering carbon emissions generated during inventory and transportation processes and obtained the cost difference point [6]. Rabta constructed an EOQ model under the circular economy and represented the circularity level of product production in exponential form [7]. Gharaei et al. developed a sustainable EOQ model for growing products under environmental costs and constraints. They found that planning projects considering sustainable development factors not only increased order quantity but also significantly shortened the growth cycle [8].
Our study is related to the research comparing inventory management under different carbon policies. Gökce Palak et al. analyzed the impact of different carbon regulation mechanisms on the costs and emissions of replenishment programs and supply chains in the context of biofuels [9]. Chung-Yuan Dye and Chih-Te Yang investigated the effects of trade credits and different carbon policies on inventory management models [10]. A bi-objective continuous review inventory model considering cost and carbon emission was developed by Konur et al. [11]. Tang et al. studied the (R, Q) inventory control models under three carbon policies: restricted percentage reduction target, carbon tax, and carbon offset [12]. Yu et al. innovatively considered supply strategy and investment in preservation technology to reduce the deterioration rate [13]. Mishra et al. constructed three inventory models under different shortage situations in the context of manufacturers investing in green technology and found that the sustainable EPQ carbon tax and cap partial backlogging model has higher profit and lower green investment cost [14]. Ruidas et al. established incomplete production inventory models under four popular carbon policies. They compensated for the shortcomings of previous studies by setting carbon emission parameters as interval numbers rather than fixed values. In contrast to the studies considering multiple carbon policies [15], As’ad et al. investigated the impact of different implementation methods of carbon quota policy (i.e., implementing in the entire planning scope or each period) on the sustainable, dynamic batch model for frozen products [16]. Turken et al. constructed a single-vendor-multi-buyer-multi-product model in which the buyer and the vendor are subject to various environmental regulations [17].
Our study is also relevant to the literature on the investment behavior of companies in emissions reduction and inventory management issues. Ayşegül Toptal et al. analyzed the inventory replenishment and carbon reduction investment decisions of retailers under three different carbon regulatory policies [18]. Li et al. developed a single-warehouse multi-retailer inventory management model and found that considering carbon emission costs increases the reorder interval for both the warehouse and retailers [19]. Lu et al. investigated the joint investment behavior of supply chain companies in reducing emissions under cap-and-trade and carbon offset policies [20]. J.-Y. Lee studied the optimal order quantity and carbon reduction investment in the EOQ model under the cap-and-price policy and found that investing in carbon reduction reduces carbon emissions per replenishment and per production unit [21]. Mishra et al. studied the replenishment problem of non-instantaneous deteriorating seasonal products based on joint pricing [22]. Mashud et al. investigated sustainable inventory issues for retailers selling non-instantaneous deteriorating products, considering green technology investment in both owned and leased warehouses, and found that investing in green technology can reduce emission costs and increase total profits [23]. Qi et al. developed a joint decision-making model for emission reduction investment and ordering under the cap-and-trade policy and investigated the impact of risk aversion behavior and investment coefficients on the optimal decision for risk-averse companies [24]. Halat et al. studied the inventory cooperation game problem of multiple supply chains under carbon tax policy, considering four structures: decentralized, downward vertical, upward vertical, and horizontal cooperation [25]. Fu et al. innovatively examined the impact of emission asymmetry and carbon tax on the production and green technology investment of competing companies [26].
The literature on production strategies and pricing issues in the context of carbon emissions is also close to our work. Du et al. considered a low carbon premium for consumers in an aggregate control and trading mechanism to promote manufacturers’ participation in low-carbon production [27]. Arda Yenipazarli used the Stackelberg game model to investigate the effect of an emissions tax on manufacturers’ optimal production and pricing strategies [28]. Xu et al. investigated the production and pricing problems of a supply chain comprising manufacturers and retailers of two made-to-order (MTO) products [29]. Liu et al. proposed three optimization models to analyze the impact of three regulatory models, mandatory carbon capacity, carbon tax, and cap-and-trade, on remanufacturers’ production quantity decisions [30]. Wang et al. investigated the manufacturing and remanufacturing problems of manufacturers considering capital and carbon emission constraints [31].
Other literature close to our study is about integrated inventory control and transportation planning issues in the context of carbon emissions. Konur and Schaefer obtained optimal order quantities for retailers under different carrier regulations and illustrated the impact of regulatory parameters on carrier preferences [32]. Hariga et al. developed three models to evaluate the impact of carbon emissions in the transportation and storage of cold chain products in a multi-stage supply chain [33]. Rout et al. constructed a sustainable supply chain inventory management model for a single supplier and multiple buyers, considering product deterioration and imperfect production under carbon emission constraints [34]. Konur et al. constructed a mixed-integer nonlinear programming model considering zoning to study integrated regional partitioning and fleet composition [35].
Several other relevant studies have explored existing emissions reduction methods or policies. Park et al. discussed the impact of levying carbon costs on supply chain structure and social welfare from the perspective of central policymakers in three different monopolistic competitive environments [36]. Zhao et al. studied the long-term joint emission reduction and low carbon propaganda in the supply chain by constructing three differential games [37]. Zhou et al. developed a three-stage game model involving consumers, firms, and the government [38]. Xu et al. considered a two-tier supply chain under aggregate control and transaction regulation and compared the profits of decentralized and centralized systems [39]. Chen et al. proposed that punishing the carbon emissions of individual companies may increase the carbon emissions of the entire supply chain due to negative externalities [40]. Lee et al. suggested that one company in the supply chain supporting another company’s carbon emissions reduction efforts benefits all participants [41]. Gao et al. compared two methods for companies to reduce emissions, namely direct investment in green technology and indirect carbon offsetting [42]. Guo et al. constructed a dynamic game theory model and compared the business models of credit wholesalers and voluntary brokers for public emissions reduction projects [43]. Kaur et al. proposed a sustainable disaster-resistant supply chain management model for procurement and logistics under a cap-and-trade policy, establishing an effective and optimal balance between a firm’s economic benefits and environmental responsibility [44]. Yenipazarli built a two-stage duopoly model to investigate the impact of consumer demand with brand preferences, competitor actions, and carbon tax policy on firms’ environmental R&D under uncertainty [45]. Y. Wang et al. studied a low-carbon supply chain consisting of a small and medium-sized manufacturer and a leading retailer. They found that altruism can improve the efficiency of the supply chain [46]. Feng et al. investigated the profit distribution rules for joint replenishment of retailers under cap-and-trade and divided retailers into two categories, efficient with altruistic behavior and inefficient. They found that the altruism parameter plays a critical role in carbon quota transfers [47]. More research articles related to this field and green emissions reduction are found in the literature review papers [48,49,50,51].
Though closely related to our study is the Chen et al. (2019) [52], their work only considers single firm EOQ production problem. Moreover, they never consider the joint production optimization and carbon caps sharing situation, which is exactly what we examine in the EPL setting in this paper. Other close work in some studies discusses integrating carbon emissions and sustainability into the supply chain from the relationship between supply chain management and environmental protection. For example, Carbon-Trust (2006) focused on developing carbon footprint measurement tools at the supply chain level [53]. After obtaining the energy consumption data of each link in the supply chain of different products, the carbon footprint of each product can be obtained by this tool. Other researchers focus on and measure the carbon footprint of the supply chain, such as Cholette and Venkat (2009) [54], Sundarakani et al. (2010) [55], and Chaabane et al. (2012) [56]. Furthermore, considering that carbon emission reduction investment often brings great cost pressure to enterprises and affects supply chain efficiency and sustainable development, there are studies on carbon emission reduction technology investment and carbon emission reduction initiative choice at the supply chain level. Chen et al. (2019) studied the EOQ problem of a single supplier and a single buyer. They found that punishing each enterprise for emissions alone might lead to higher overall emissions of the supply chain [40].

3. Decentralized System under Regulative Low-Carbon Policies

This study starts from the classic economic production quantity (EPQ) model to focus on multiple firms constrained by regulative low-carbon policies. Consider N manufacturers make production decisions independently in a low-carbon environment. As stated in the EPQ model, manufacturer i faces his own deterministic market demand d i . And his production rate is a constant p i , which can be regarded as production capacity. To let the following formulation make sense, we need condition d i / p i < 1 . Manufacturers organize production and seek optimal output or supply levels by considering demand characteristics and production capacity comprehensively. In the framework of the joint inventory optimization model based on the extended model of economic production quantity (EPQ) model, we consider that each manufacturer produces one kind of product with the demand d i > 0 and variable production cost c i . Similarly, we assume the supply rate a constant p i 0 . The ordering cost and the storage cost rates are a i 0 and h i 0 , respectively. We eliminate the impact of environmental change on inventory costs as well as the limit of order quantity and storage capacity.
In the low-carbon era, the regulative carbon policies will substantially affect the business and economic activities relative to the situation without carbon constraints. We, therefore, consider three common policies, i.e., mandatory cap, cap and trade, and carbon taxation, and examine their impacts on supply chain operations. To reflect the acting effect of carbon policies, we assume that all ordering, inventory, and production will generate carbon emissions, the corresponding rate coefficient of which is shown in Table 1. Furthermore, assume that these coefficients do not change over time.
Note that the time dimension of the above notations is measured in years for all i N . Throughout the paper, we attach an embellishment hat “ ” to indicate their relationship to carbon emission.
In this part, we first examine the scenario of completely decentralized decision-making for a dyadic supply chain and then consider the situation with a single supplier and multiple manufacturers. As mentioned above, the order quantity Q i implies that the order period and frequency satisfy L i = Q i / d i and m i = d i / Q i , respectively. The retailer’s costs are made up of the annual ordering cost a i d i / Q i , the annual inventory cost h i Q i ( p i d i ) / ( 2 p i ) , and the annual production cost c i d i , where Q i ( p i d i ) / ( 2 p i ) represents the annual average inventory level. The carbon emission associated with the supply chain operations includes ordering-incurred emission a ^ i d i / Q i , inventory-incurred emission h ^ i Q i ( p i d i ) / ( 2 p i ) , and production-incurred emission c ^ i d i .

3.1. Carbon Tax or Cap-and-Tarde Policy-Based Optimal Production Quantity

In a low-carbon economy, it is natural to think about the impact of low-carbon policies on economic activities and operations management. Therefore, we consider the impact of two typical decarbonization schemes, i.e., carbon taxation and subsidy, on supply chain operations. We refer to carbon taxation and cap-and-trade policies as carbon pricing policies since both charge costs according to the units of emission permits. The pricing of carbon taxation reflects the emission penalty, while the pricing in cap and trade regulation represents the opportunity cost of consuming carbon rights. Unlike prior literature, we mathematically formulate these two schemes in a unified framework as follows, except that the term C G H G i vanishes under the cap taxation policy. We give the n-firm system total cost A C ( Q i , , Q n ) as follows:
min Q i 0   A C ( Q i , , Q n ) = i = 1 N [ a i d i Q i + h i p i d i 2 p i Q i + c i d i + τ ( a ^ i d i Q i + h ^ i p i d i 2 p i Q i + c ^ i d i C G H G i ) ]
where τ ( τ > 0 ) represents the tax rate under carbon tax policy or emission permit price under cap and trade policy, which can be regarded as the governmental penalty rate on carbon emission. Here we denote the term E ^ i = a ^ i d i Q i + h ^ i p i d i 2 p i Q i + c ^ i d i C G H G i as the net carbon emission generated under each policy, and C G H G i disappears as cap tax regulation applies. Under the cap-and-trade policy, the positive (negative) E ^ i means permits purchasing cost (sales revenue). Analyzing the general total cost incorporating net emission under taxation or cap and trade gives the proposition as follows:
Corollary 1. 
The optimal solutions of the economic production lot model for a single firm under carbon tax policy or cap and trade policy gives optimal production quantity  Q i , τ = a i + τ a ^ i h i + τ h ^ i 2 p i d i p i d i ,  frequency  m i , τ = h i + τ h ^ i a i + τ a ^ i d i ( p i d i ) 2 p i ,  emission  E ^ i , τ = ( a ^ i h i + τ h ^ i a i + τ a ^ i + h ^ i a i + τ a ^ i h i + τ h ^ i ) d i ( p i d i ) 2 p i + c ^ i d i  and total cost  A C = 2 ( a i + τ a ^ i ) ( h i + τ h ^ i ) d i ( p i d i ) 2 p i + c i d i + τ ( c ^ i d i C G H G i )  with the following features: (1)  A C τ 0  holds if  C G H G i 2 a ^ i h ^ i d i ( p i d i ) p i + c ^ i d i ; (2)  A C τ τ < 0  if  C G H G i > ( h ^ i a i h i + a ^ i h i a i ) d i ( p i d i ) 2 p i + c ^ i d i ; (3)  A C τ τ > 0  for  0 < τ < τ 0  and  A C τ τ 0  for  τ τ 0  if  2 a ^ i h ^ i d i ( p i d i ) p i + c ^ i d i + c ^ i d i < C G H G i < ( h ^ i a i h i + a ^ i h i a i ) d i ( p i d i ) 2 p i + c ^ i d i ;  (4) Systemwide emission  E ^ τ  always decreases in emission penalty rate  τ , where  τ 0  is the unique solution of equation  ( h ^ i a i + τ a ^ i h i + τ h ^ i + a ^ i h i + τ h ^ i a i + τ a ^ i ) d i ( p i d i ) 2 p i = C G H G i c ^ i d i  if the following condition  2 a ^ i h ^ i d i ( p i d i ) p i < C G H G i c ^ i d i < ( h ^ i a i h i + a ^ i h i a i ) d i ( p i d i ) 2 p i  holds.
Proposition 1. 
The optimal solutions of economic production lot model for multiple firms under carbon tax policy or cap and trade policy gives optimal production quantity  Q i , τ = a i + τ a ^ i h i + τ h ^ i 2 p i d i p i d i ,  total emission  E ^ τ = i N [ ( a ^ i h i + τ h ^ i a i + τ a ^ i + h ^ i a i + τ a ^ i h i + τ h ^ i ) d i ( p i d i ) 2 p i + c ^ i d i ] = E ˜ τ ( τ ) + i N c ^ i d i  and systemwide cost  A C = i = 1 n [ 2 ( a i + τ a ^ i ) ( h i + τ h ^ i ) d i ( p i d i ) 2 p i + c i d i + τ ( c ^ i d i C G H G i ) ]  with following features: (1)  A C τ 0  holds if  i N C G H G i inf { E ˜ τ } + i N c ^ i d i ; (2)  A C τ < 0  if  i N C G H G i > sup { E ˜ τ } + i N c ^ i d i ; (3)  A C τ τ > 0  for  0 < τ < τ 0  and  A C τ τ 0  for  τ τ 0  if  inf { E ˜ τ } + i N c ^ i d i < i N C G H G i < sup { E ˜ τ } + i N c ^ i d i ;  (4) Systemwide emission  E ^ τ  always decreases in emission penalty rate  τ , where  τ 0  is the unique solution of equation  i N [ ( h ^ i a i + τ a ^ i h i + τ h ^ i + a ^ i h i + τ h ^ i a i + τ a ^ i ) d i ( p i d i ) 2 p i ] = i N ( C G H G i c ^ i d i )  if the following condition  inf { E ˜ τ }   +   i N c ^ i d i   <   i N C G H G i   <   sup { E ˜ τ }   +   i N c ^ i d i  holds, where we denote  inf { E ˜ τ }   =   i = 1 N [ 2 a ^ i h ^ i d i ( p i d i ) p i ]  and  sup { E ˜ τ } = i = 1 N [ ( h ^ i a i h i + a ^ i h i a i ) d i ( p i d i ) 2 p i ] .
The proofs of Corollary 1 and Proposition 1 are shown in Appendix A. The above summary shows that the joint production decision among multiple enterprises can overcome the problem caused by the insufficient carbon emission of a single enterprise. Under this situation, enterprises have a broader decision-making space. How the carbon penalty rate affects the total cost depends on the comparison between the total emission quota and several thresholds. Specifically, increasing the penalty rate will reduce (resp. increase) the total cost when the total carbon quota is greater (resp. less) than a certain threshold, while increasing the rate will make the total cost increase and decrease when the total carbon quota is between certain thresholds. However, the associated optimal emission is always decreasing with respect to the carbon price. This is intuitive since the higher carbon price can better curb the total emission. Nevertheless, it does not always benefit the control of overall cost since a balance exists between the operational and emission costs.
Proposition 2. 
The impact of the emission penalty rate  τ  on the optimal production for multiple firms exhibits the following features:  Q i , τ τ 0  if  a i h i a ^ i h ^ i , and  Q i , τ τ < 0  otherwise.
For manufacturer i, the problem of solely optimizing operational cost is equivalent to the situation without carbon emission constraints. Therefore, from Equation (1) can be directly derived the associated order quantity Q i * = lim τ 0 Q i , τ = 2 a i d i p i h i ( p i d i ) , optimal ordering frequency m i * = d i Q i * = h i ( p i d i ) d i 2 a i p i and minimum average cost A C ( Q i * ) =   2 a i m i * + c i d i . Further analysis shows the equal optimal ordering cost and inventory cost a i m i * . The carbon emission of solely operational cost optimization can be expressed as follows:
E M ( Q i * ) = lim τ 0 E ^ i , τ = ( a ^ i h i a i + h ^ i a i h i ) d i ( p i d i ) 2 p i
Accordingly, we can get the optimal inventory cost A C o ( Q i * ) and optimal ordering cost A C h ( Q i * ) as shown in Table 2.
Optimizing the emission cost solely is the other angle to observe the objective function of total cost in Formula (1). Taking the same manner as in the scenario of solely optimizing operational cost to examine the obtained Q i , τ , E ^ i , τ and A C i , τ yields the emission cost solely optimization-based optimal values, such as production quantity Q ^ min , i = lim a i , h i 0 Q i , τ = 2 a ^ i p i d i h ^ i ( p i d i ) , ordering frequency m ^ min , i = h ^ i d i ( p i d i ) 2 a ^ i p i , average cost A C ( Q ^ min , i ) = lim a i , h i 0 A C i , τ = τ [ 2 a ^ i h ^ i d i ( p i d i ) / p i + ( c ^ i d i C G H G i ) ] , and carbon emission E M ( Q ^ min , i ) = lim a i , h i 0 E ^ i , τ = 2 a ^ i h ^ i d i ( 1 d i / p i ) + c ^ i d i .   E M ( Q ^ min , i ) is increasing in production rate p i , which is consistent with intuition. Furthermore, A C ( Q ^ min , i ) and E M ( Q ^ min , i ) degenerate into economic order quantities when d i / p i 0 , while the cost and carbon emission tend to be infinite when d i / p i 1 . Therefore, our model under EPL situation is substantially general to incorporate the EOQ-based settings, which can be considered a special case of ours. The pairs of system decision variables from the perspective of cost optimization and carbon emission are shown in Table 3.

3.2. Mandatory Cap Policy-Based Optimal Production Quantity

In this section, we check the same system under the situation with another regulative policy, the mandatory cap scheme, which differs from the previous carbon pricing policy. Actually, the mandatory cap means a very strict constraint on the total emission, which has a distinct acting principle and mechanism compared with the carbon pricing regime. Therefore, the economic production quantity optimization problem constrained by the mandatory carbon emission cap can be expressed as the following nonlinear optimization model:
min Q i 0   A C ( Q i ) = a i d i Q i + h i p i d i 2 p i Q i + c i d i s . t .   E M ( Q i ) = a ^ i d i Q i + h ^ i p i d i 2 p i Q i + c ^ i d i C G H G i   for   i N
Note that A C ( Q i )  is the annual total operating cost when the production quantity is Q i , and E M ( Q i )  represents the corresponding annual total carbon emissions. The constraint in the above model means that the annual total carbon emissions cannot exceed the total constraint C G H G i . The left-hand side expression of constraint contains emissions incurred by delivering and storing goods, respectively.
From the above description, it can be seen that to solve the inequality, the limit of mandatory carbon emission constraint should satisfy the inequality below:
C G H G i E M ( Q min , i ) = 2 a ^ i h ^ i d i ( p i d i ) / p i + c ^ i d i
Otherwise, enterprises will be unable to meet the carbon emission constraint no matter how they organize production and order, and this situation will also make the mandatory carbon emission reduction regulation meaningless. Moreover, the mandatory carbon cap requires the actual production Q to meet the condition E M ( Q i ) C G H G i .
Therefore, solving the above model yields the optimal economic production quantity Q i * for minimizing annual operational costs under the constraint of a mandatory cap. A simple analysis of the constrained nonlinear optimization model shows that the goal of enterprise operation is to minimize the operational cost and meet the carbon emission constraints. To better reflect these two goals, we need to analyze them separately.
Based on the above perspectives, we investigate the nonlinear optimization problem considering both carbon emission and cost optimization constraints. First, observing the constraint, we can rewrite the inequality into the following expression:
h ^ i p i d i 2 p i Q i 2 ( C G H G i c ^ i d i ) Q i + a ^ i d i 0
Two roots of the equation are Q i , 1 = p i [ C ^ G H G i C ^ 2 G H G i 2 a ^ i h ^ i ( p i d i ) d i / p i ] h ^ i ( p i d i ) and Q i , 2 = p i [ C ^ G H G i + C ^ 2 G H G i 2 a ^ i h ^ i ( p i d i ) d i / p i ] h ^ i ( p i d i ) , where we denote C ^ G H G i = C G H G i c ^ i d i . The solution of equation (5) shows that the value range of the production quantity that meets the constraint is:
Q i , 1 Q i Q i , 2
In this way, the requirements for the value of economic production quantity are very clear.
Considering the objective function A C ( Q i ) , a convex function of production quantity Q i , and the evaluation of decision variables from the perspective of solely cost optimization in Table 2, the optimal solution of above constrained nonlinear programming can be obtained as follows:
Q ^ i = { Q i ,   if   Q i , 1 Q i Q i , 2 ; Q i , 1 ,   if   Q i Q i , 1 ; Q i , 2 ,   if   Q i Q i , 2 ;   for   i { 1 , 2 , , n }
Equation (7) shows that the solely cost optimization-based production quantity and the production boundary determined by pure carbon emission constraints together determine the optimal economic production quantity in the supply chain subject to mandatory cap policy. The above formula shows that the optimal economic production quantity problem in a low-carbon environment can be equivalent to the problem without carbon constraint only when the feasible range of emission constraint meets certain conditions, that is, the emission constraint fails to act in affecting the production quantity from the perspective of solely cost optimization. From the opposite point of view, the governmental mandatory carbon emission restriction policy can play its due role only if it is set at a reasonable level. The inequalities A C ( Q i ) A C ( Q ^ i , min ) and E M i ( Q i ) E M i ( Q ^ i , min ) show a contradictory relationship between the optimal operational cost of solely cost optimization and optimal carbon emission of solely emission optimization.
Considering system decision from the perspective of utility, we can incorporate the emission constraint as a penalty term in the objective function as follows, which has a similar format to Formula (1).
min Q i 0   A C s ( Q i ) = a i d i Q i + h i p i d i 2 p i Q i + c i d i + s i [ E M ( Q i ) C G H G i ] = ( a i + s i a ^ i ) d i Q i + ( h i + s i h ^ i ) p i d i 2 p i Q i + c i d i + s i ( c ^ i d i C G H G i )
In the above formula, si can be regarded as the unit penalty rate of carbon emission rights exceeding the quota or the equivalence benefit of quota surplus. When there are no other constraints in Equation (8), the optimal value of the corresponding decision variable is Q i * * = 2 p i d i p i d i a i + s i a ^ i h i + s i h ^ i , having a similar mathematical structure with the constraint-free optimal decision Q i * = 2 p i d i p i d i a i h i and the solely emission optimization-based decision Q E M i * = 2 p i d i p i d i a ^ i h ^ i . Comparing these optimal decisions gives the following Proposition 3 from a utility perspective.
Proposition 3. 
Under mandatory carbon cap, investigating the relationship between the ratio of ordering and inventory holding associated carbon emission rates and the ratio of cost rates of the same kind yields: (1) If  a ^ i h ^ i a i h i , according to the above analysis,  C G H G i  can enable  Q E M i *  available, and accordingly  Q ^ i * = Q i * *  will be the systemwide optimal value; (2) If  a ^ i h ^ i < a i h i ,the systemwide optimal solution in this situation is:
Q ^ i * = { Q i * * , if   Q i , 1 Q i * * Q i , 2 ; Q i , 1 , if   Q i * * Q i , 1 ; Q i , 2 , if   Q i * * Q i , 2 ;
Proof. 
According to a i + s i a ^ i h i + s i h ^ i a i h i a ^ i h ^ i a i h i and a i + s i a ^ i h i + s i h ^ i a ^ i h ^ i a i h i a ^ i h ^ i , the relationship between the optimal order decision in the three scenarios is Q E M i * Q i * * Q i *   if   a ^ i h ^ i a i h i and Q i * Q i * * Q E M i * if a i h i a ^ i h ^ i . Therefore, the above proposition can be obtained. □

4. Joint Production Stocking Decision under Carbon Caps Sharing

This Section still considers the same supply chain structure as Section 2 but investigates different scenarios. Because retailers have certain differences in various operating parameters and carbon emission constraint parameters, they greatly differ in considering carbon emission constraints in simultaneously optimizing operations. This difference allows enterprises to take advantage of the differences in emission reduction endowments to better organize and make operational decisions. Therefore, this section examines the situation when retailers make completely independent decentralized decisions among operators but share carbon emissions. This situation enables enterprises to maintain good independence in operation but form an alliance to deal with carbon emissions and jointly face carbon emission regulations.
Based on the results in the above section, we can obtain the optimal decentralized operational decision problem for joint emission reduction under the scenario of decentralized operational decision-making but shared carbon emission allowances, as shown in the constrained nonlinear programming below:
min A C ( Q 1 , Q 2 , , Q n ) = i N [ a i d i Q i + h i p i d i 2 p i Q i + c i d i ]
s . t . E M ( Q 1 , Q 2 , , Q n ) = i N [ a ^ i d i Q i + h ^ i p i d i 2 p i Q i + c ^ i d i ] C ˜ G H G
Including C ˜ G H G = i N C G H G i , i N [ 2 a ^ i h ^ i d i ( p i d i ) / p i + c ^ i d i ] C ˜ G H G is satisfied.
The Hessian matrix of the above formula is:
[ 2 A C / Q 1 2 , 2 A C / Q 1 Q 2 , , 2 A C / Q 1 Q n 2 A C / Q 1 Q 2 , 2 A C / Q 2 2 , , 2 A C / Q 2 Q n 2 A C / Q n Q 1 , 2 A C / Q n Q 2 , , 2 A C / Q n 2 ] = [ 2 a 1 d 1 Q 1 3 , 0 , , 0 0 , 2 a 2 d 2 Q 2 3 , , 0 0 , 0 , , 2 a n d n Q n 3 ]
Therefore, the function A C ( Q 1 , Q 2 , , Q n ) is convex since the Hessian matrix H A C is positively definite, and there is a global minimum. Similarly, the Hessian matrix of the function E M ( Q 1 , Q 2 , , Q n ) is also positively definite with a global minimum.
The solving the nonlinear programming problem under this situation yields Propositions 4–6, and their related proofs.
Proposition 4. 
The solution of Equations (10) and (11) is:
Q ˜ i = { 2 a i p i d i h i ( p i d i ) = Q i ,   when   Q Γ a i + γ a ^ i h i + γ h ^ i 2 p i d i p i d i ,   when   Q Γ ( + n \ Γ )
where  Γ = { ( Q 1 , Q 2 , , Q n ) | i N [ a ^ i d i Q i + h ^ i p i d i 2 p i Q i + c ^ i d i ] C ˜ G H G } ,  Q = ( Q 1 , Q 2 , , Q n ) ,  Γ  is the boundary of  Γ  and  γ  the Lagrange multiplier determined by Equation (15).
Proof. 
The Lagrange function is L A C ( Q 1 , Q 2 , , Q n ) = i N [ a i d i Q i + h i p i d i 2 p i Q i + c i d i ] γ { C ˜ G H G i N [ a ^ i d i Q i + h ^ i p i d i 2 p i Q i + c ^ i d i ] } , the corresponding K-T condition includes γ { C ˜ G H G i N [ a ^ i d i Q i + h ^ i p i d i 2 p i Q i + c ^ i d i ] } = 0 ,   L A C = A C γ E M = 0 and γ 0 . Considering L A C = [ , ( a i + γ a ^ i ) d i Q i 2 + ( h i + γ h ^ i ) p i d i 2 p i , ] , that is Q ˜ i = a i + γ a ^ i h i + γ h ^ i 2 p i d i p i d i , which is derived from the function ( a i + γ a ^ i ) d i Q i 2 + ( h i + γ h ^ i ) p i d i 2 p i = 0 .
If γ = 0 , Q ˜ i = 2 a i p i d i h i ( p i d i ) = Q i , then Equation (13) is established.
( Q 1 , Q 2 , , Q n ) Γ = { ( Q 1 , Q 2 , , Q n ) | i N [ a ^ i d i Q i + h ^ i p i d i 2 p i Q i + c ^ i d i ] C ˜ G H G }
and there is A C ( Q 1 , Q 2 , , Q n ) = i N A C i ( Q i ) = i N [ 2 a i h i d i ( p i d i ) / p i + c i d i ] .
E M ( Q 1 , Q 2 , , Q n ) = i N E M i ( Q i ) = i N [ a ^ i h i d i ( p i d i ) 2 a i p i + h ^ i a i d i ( p i d i ) 2 h i p i + c ^ i d i ]
If γ > 0 , then
Q ˜ i = Q i , γ = a i + γ a ^ i h i + γ h ^ i 2 p i d i p i d i
where ( Q ˜ 1 , Q ˜ 2 , , Q ˜ n ) { ( Q ˜ 1 , Q ˜ 2 , , Q ˜ n ) | C ˜ G H G = i N [ a ^ i d i Q ˜ i + h ^ i p i d i 2 p i Q ˜ i + c ^ i d i ] } , that is
C ˜ G H G = i N [ ( a ^ i h i + γ h ^ i a i + γ a ^ i + h ^ i a i + γ a ^ i h i + γ h ^ i ) ( p i d i ) d i 2 p i + c ^ i d i ]
The optimal yield and Lagrange multiplier values for decentralized decision-making but shared carbon emission constraints are determined by Equations (14) and (15) simultaneously, that is, the optimal yield can be obtained by substituting γ into (14). □
Proposition 5. 
The prerequisite for Equation (15) to have a positive real solution is:
i N ( 2 a ^ i h ^ i ( p i d i ) d i / p i + c ^ i d i ) C ˜ G H G i N [ ( a ^ i h i a i + h ^ i a i h i ) ( p i d i ) d i 2 p i + c ^ i d i ]
Proof. 
Let the right side of the Equation (15) be expressed as f ( γ ) = i N [ ( a ^ i h i + γ h ^ i a i + γ a ^ i + h ^ i a i + γ a ^ i h i + γ h ^ i ) ( p i d i ) d i 2 p i + c ^ i d i ] , then get
f ( γ ) = i N [ ( h ^ i a i h i a ^ i ) 2 2 ( a i + γ a ^ i ) 3 / 2 ( h i + γ h ^ i ) 3 / 2 ( p i d i ) d i 2 p i ] 0
f ( γ ) = i N [ ( a ^ i a i + γ a ^ i + h ^ i h i + γ h ^ i ) 3 ( h ^ i a i h i a ^ i ) 2 4 ( a i + γ a ^ i ) 3 / 2 ( h i + γ h ^ i ) 3 / 2 ( p i d i ) d i 2 p i ] 0
From Equations (17) and (18), the function f ( γ ) decreases monotonically. There is a unique real root γ 0 for Equation (15) if it exists.
lim γ f ( γ ) f ( γ ) f ( 0 ) = i N [ ( a ^ i h i a i + h ^ i a i h i ) ( p i d i ) d i 2 p i + c ^ i d i ]
where lim γ f ( γ ) = i N ( 2 a ^ i h ^ i ( p i d i ) d i / p i + c ^ i d i ) , therefore, Formula (16) is established. □
To make a more detailed explanation of the above propositions, without losing generality, we give the following proposition reflecting the optimal decision of homogeneous enterprises.
Proposition 6. 
For N homogeneous enterprises in the situation of joint carbon emission reduction, the following expression can be established:
a i + γ a ^ i h i + γ h ^ i = { 1 2 h ^ i ( 1 n C ˜ G H G c ^ i d i ) 2 p i ( p i d i ) d i     2 p i 4 h ^ i 2 ( p i d i ) d i ( 1 n C ˜ G H G c ^ i d i ) 2 a ^ i h ^ i ,   if   a ^ h ^ a h 1 2 h ^ i ( 1 n C ˜ G H G c ^ i d i ) 2 p i ( p i d i ) d i   +   2 p i 4 h ^ i 2 ( p i d i ) d i ( 1 n C ˜ G H G c ^ i d i ) 2 a ^ i h ^ i ,   if   a h a ^ h ^
γ = { h i [ 2 p i 4 h ^ i 2 ( p i d i ) d i ( 1 n C ˜ G H G c ^ i d i ) 2     2 p i 4 h ^ i 2 ( p i d i ) d i ( 1 n C ˜ G H G c ^ i d i ) 2 a ^ i h ^ i ] 2 a i a ^ i h ^ i [ 2 p i 4 h ^ i 2 ( p i d i ) d i ( 1 n C ˜ G H G c ^ i d i ) 2     2 p i 4 h ^ i 2 ( p i d i ) d i ( 1 n C ˜ G H G c ^ i d i ) 2 a ^ i h ^ i ] 2 , if   a ^ h ^ a h h i [ 2 p i 4 h ^ i 2 ( p i d i ) d i ( 1 n C ˜ G H G c ^ i d i ) 2   +   2 p i 4 h ^ i 2 ( p i d i ) d i ( 1 n C ˜ G H G c ^ i d i ) 2 a ^ i h ^ i ] 2 a i a ^ i h ^ i [ 2 p i 4 h ^ i 2 ( p i d i ) d i ( 1 n C ˜ G H G c ^ i d i ) 2   +   2 p i 4 h ^ i 2 ( p i d i ) d i ( 1 n C ˜ G H G c ^ i d i ) 2 a ^ i h ^ i ] 2 , if   a h a ^ h ^
Proof. 
Make x = ( a i + γ a ^ i ) / ( h i + γ h ^ i ) , the roots of Equation (15) are x 1 and x 2 , where we have x 1 = 1 2 h ^ i ( 1 n C ˜ G H G c ^ i d i ) 2 p i ( p i d i ) d i 2 p i 4 h ^ i 2 ( p i d i ) d i ( 1 n C ˜ G H G c ^ i d i ) 2 a ^ i h ^ i and x 2 = 1 2 h ^ i ( 1 n C ˜ G H G c ^ i d i ) 2 p i ( p i d i ) d i + 2 p i 4 h ^ i 2 ( p i d i ) d i ( 1 n C ˜ G H G c ^ i d i ) 2 a ^ i h ^ i .
As mentioned previously, a h a + s a ^ h + s h ^ a ^ h ^   if   a ^ h ^ a h   and   a ^ h ^ a + s a ^ h + s h ^ a h   if   a h a ^ h ^ were established. The inequality - 2 p i 4 h ^ i 2 ( p i d i ) d i ( 1 n C ˜ G H G c ^ i d i ) 2 a ^ i h ^ i < 2 p i 4 h ^ i 2 ( p i d i ) d i ( 1 n C ˜ G H G c ^ i d i ) 2 a ^ i h ^ i 2 p i 4 h ^ i 2 ( p i d i ) d i ( 1 n C ˜ G H G c ^ i d i ) 2 a ^ i h ^ i and a ^ i h ^ i 2 p i 4 h ^ i 2 ( p i d i ) d i ( 1 n C ˜ G H G c ^ i d i ) 2 also holds. There are x = x 1   if   a ^ h ^ a h and x = x 2   if   a h a ^ h ^ . Then solving equation x 1 = a i + γ a ^ i h i + γ h ^ i and x 2 = a i + γ a ^ i h i + γ h ^ i gives the Equation (21). Therefore, the above proposition can be proved.
Substituting Q ˜ i into constraint (11) yields inequality E M ( Q ˜ 1 , Q ˜ 2 , , Q ˜ n ) i N E M i ( Q ^ i , min ) . Thus, we have A C ( Q 1 , Q 2 , , Q n ) = i N [ ( a i h i + γ h ^ i a i + γ a ^ i + h i a i + γ a ^ i h i + γ h ^ i ) ( p i d i ) d i 2 p i + c i d i ] and total emissions C ˜ G H G . Retailer i’s emissions is E M i ( Q ˜ i )   =   ( a ^ i h i + γ h ^ i a i + γ a ^ i   +   h ^ i a i + γ a ^ i h i + γ h ^ i ) ( p i d i ) d i 2 p i   +   c ^ i d i     E M i ( Q ^ i , min ) . □
The comparison of centralized and decentralized situations is shown in Table 4.
Observation 1. 
From previous Propositions, we can have  C ˜ G H G i N [ ( a ^ i h i a i + h ^ i a i h i ) ( p i d i ) d i 2 p i + c ^ i d i ] .  The optimal solution satisfying (10) and (11) is  Q ˜ = ( Q 1 , Q 2 , , Q n ) ,where the constraint (11) is inactive and retailers decide according to the state without carbon emission constraint. Hence, to reduce the carbon emission of the whole supply chain system, the total carbon emission limit of the system should meet the condition Formula (16). Particularly, Formula (16) does not imply that the constraint in Formula (3) will be satisfied for every retailer (e.g.,  i ). On the contrary, it is possible that the retailer does not satisfy Formula (3) but satisfies the overall carbon emission constraint (16) when considering systemwide optimization. That is,  Q ^ i = Q i  while  Q ˜ i = Q i , γ  or  Q ^ i = Q i , 1 ( or   Q i , 2 )  while  Q ˜ i = Q i , γ  for retailer i.
The following discussion mainly focuses on the situations with condition (16) established.
Proposition 7. 
In the situation of carbon quota sharing, as follows is the relationship between carbon emissions-oriented optimal decision  Q ˜ i , total cost-oriented optimal decision  Q i , and optimal systemwide decision  Q ^ i  under different parameters configurations:
(1) 
If  a ^ h ^ a h  for retailer i, then  Q i < Q ˜ i = a i + γ a ^ i h i + γ h ^ i 2 p i d i p i d i = a i ( 1 + γ a ^ i / a i ) h i ( 1 + γ h ^ i / h i ) 2 p i d i p i d i < Q ^ i , min .
(2) 
If  a ^ h ^ < a h  for retailer i, then  Q ^ i , min < Q ˜ i = a i + γ a ^ i h i + γ h ^ i 2 p i d i p i d i = a i ( 1 + γ a ^ i / a i ) h i ( 1 + γ h ^ i / h i ) 2 p i d i p i d i < Q i .
(3) 
If  a ^ h ^ = a h , then  Q ˜ i = Q i = Q ^ i , min .
The above Proposition 7 shows the influence of different parameter Settings on the optimal decisions of three scenarios; that is, comparing the ratio between carbon emission rates caused by production lot and inventory holding and the ratio between cost rates caused by them will determine the size relationship between optimal decisions.
The above formulation mainly discusses and compares the optimal production plan and corresponding cost under the two situations of individual decision-making and joint decision-making. Joint decision-making means that all enterprises implement centralized decision-making based on sharing carbon emission quota, that is, to make joint optimal production arrangements from the perspective of maximizing system benefits. The purpose of exploring the significance of centralized decision-making is to demonstrate that the joint optimal decision-making in this situation can achieve better operational performance than the traditional decentralized decision-making situation to demonstrate that the government’s relaxation of the carbon quota use standard has important social and economic benefits. This analysis lays a good foundation for designing a reasonable and feasible guarantee mechanism for implementing carbon quota sharing and joint optimal production decision-making among enterprises in the case of low-carbon policy improvement.

5. The Comparative Analysis of Two Scenarios

5.1. Comparison of Two Decision-Making Situations under Carbon Emission Constraints

In this part, we mainly analyze the optimal operational cost and carbon emission under two situations, i.e., completely decentralized decision and joint emission reduction.
(1) If Q ^ i = Q i for i N , that is, constraint (16) is established as constraint (2) is satisfied, thus Q ˜ i = Q i , γ .
(2) Q ^ i = Q i , 1   or   Q i , 2 holds for i W , where W   = { i | a ^ i h i d i ( p i d i ) 2 a i p i   +   h ^ i a i d i ( p i d i ) 2 h i p i   +   c ^ i d i   C G H G i } , and Q ^ j = Q j for j N \ W ,  then. In the case of joint production decisions with carbon emission caps sharing for retailers i N , there are two situations concerning whether Formula (16) is established or not:
(2-i) If the following inequality holds i N ( 2 a ^ i h ^ i ( p i d i ) d i / p i + c ^ i d i ) C ˜ G H G i N [ ( a ^ i h i a i + h ^ i a i h i ) ( p i d i ) d i 2 p i + c ^ i d i ] , i.e., Formula (16) still established. And then the optimal production economic batch of retailers is: Q ˜ i = Q i , γ .
(2-ii) If C ˜ G H G > i N [ ( a ^ i h i a i + h ^ i a i h i ) ( p i d i ) d i 2 p i + c ^ i d i ] , then Formula (16) is not established. At this point,  in Formula (15) have no positive real solution according to prior propositions. In this situation there is Q ˜ i = 2 a i p i d i h i ( p i d i ) = Q i , that is, all retailers can make decisions based solely on cost optimization.
In situation (2-ii), the carbon emission constraint does not affect the cost-optimal decision, and there is no difference between scenarios of solely cost-optimizing and solely emission decision. When some enterprises belong to set W and others belong to set N \ W , sharing carbon emissions will realize the reallocation of carbon emissions in the retailer alliance, regardless of the case (2-i) or (2-ii).
Proposition 8. 
Let the left-hand side of inequality (11) be function  E M ( Q ) = i N E M i ( Q i ) , including vector  Q = [ Q 1 , Q 2 , , Q n ] ,  E M i ( Q i )  is expressed as in (2). The lower bound of total system cost under completely decentralized and joint emission reduction decision makings is related as follows:
inf { A C ( Q ) : E M ( Q ) C ˜ G H G } inf { i N A C i ( Q i ) : E M i ( Q i ) C G H G i }
Proof. 
First, we examine the domain of the objective function defined by the constraints of the two situations. Obviously, there is
{ ( Q 1 , Q 2 , , Q n ) | a ^ d Q i + h ^ p d 2 p Q i + c ^ d C G H G i ,   for   all   i N } { ( Q 1 , Q 2 , , Q n ) | i N [ a ^ i d i Q i + h ^ i p i d i 2 p i Q i + c ^ i d i ] i N C G H G i }
That is, { ( Q 1 , Q 2 , , Q n ) | E M i ( Q i ) C G H G i for   all   i N } { ( Q 1 , Q 2 , , Q n ) | E M ( Q ) C ˜ G H G } holds.
Although A C ( Q ) and i N A C i ( Q i ) are equal in value, using A C ( Q ) instead of the latter expression aims to distinguish two decision scenarios. Therefore, when the function expressions are equal, the lower bound of the function with a large value space of independent variables must not be greater than the lower bound of the function with a small value space, so Equation (21) is established and Proposition 8 is proved. □
The discussion above and Proposition 8 show that enterprises strictly constrained by carbon emission allowances in the case of independent decision-making can relax the carbon emission constraints by sharing emission allowances so that all retailers in their system make decisions according to solely cost optimization. Essentially, the enterprises with more carbon emission limits (in the relaxed state of constraints) will transfer part of their allowances to the enterprises with tight emission limits, making the cost of the whole system lower. In essence, firms with relatively more carbon allowances (the constraints are relaxed) will transfer some of their allowances to firms with tight allowances, making the systemwide cost less.
A C ( Q ˜ 1 , Q ˜ 2 , , Q ˜ n ) = inf { A C ( Q ) : E M ( Q ) C ˜ G H G } inf { i N A C i ( Q i ) : E M i ( Q i ) C G H G i } i N A C i ( Q ^ i )
Therefore, the difference between optimal total operational costs under two situations is:
Δ A C = i N A C i ( Q ^ i ) A C ( Q ˜ 1 , Q ˜ 2 , , Q ˜ n ) = { i W [ a i d i Q i , k + h i p i d i 2 p i Q i , k + c i d i ] + i N \ W [ 2 a i h i d i ( p i d i ) / p i + c i d i ] i N [ ( a i h i + γ h ^ i a i + γ a ^ i + h i a i + γ a ^ i h i + γ h ^ i ) ( p i d i ) d i 2 p i + c i d i ] ,   if   W ; 0 ,   if   W = .
where k { 1 , 2 } .
It can be seen that the carbon quotas sharing scenario can meet the overall carbon emission constraint while reducing the systemwide operational cost of the downstream retailer’s alliance.

5.2. Numerical Analysis

This subsection takes a one-supplier triple-retailer supply chain as an example for computational study. Assume that only one product is produced and sold and that the market demand and supply characteristics are consistent with the assumptions made at the beginning. The setting of parameters is shown as follows:
NotationDescription
p [2.5, 5.6, 4.8]
d [1.2, 4.1, 2.9]
a [12.3, 13.4, 15.7]
h [0.61, 0.38, 0.50] Annual inventory carrying cost per unit product
c [5.8, 4.2, 4.6] Variable cost of production per unit of product
a ^ [2.3, 4.7, 3.6] Carbon emissions produced by a single order
h ^ [0.017, 0.023, 0.033] Carbon emissions caused by annual inventory holding per unit of product
c ^ [0.25, 0.18, 0.22] Carbon emissions per unit product
C G H G [4.7, 6.3, 5.8] Annual total inventory emissions
Q Economic production batch, decision variables
Q * Optimal economic production batch without carbon emission constraint
Q min Optimal economic production batch with carbon emission minimized
Q ^ * The optimal economic production batch with carbon emission constraint and cost minimization is considered simultaneously
Q ˜ * Optimal economic production batch for joint operation decisions
The numerical results in Table 5 verify the propositions and related conclusions in the preceding section. The results of completely decentralized decision-making and carbon quota-sharing decision-making show that sharing carbon emission allowances can reduce the total operating cost of the overall system of the alliance but, at the same time, meet the total carbon emission quota constraints of the system. At the same time, the actual carbon emission may be increased, but the total social utility may still be improved.

6. Concluding Remarks

Inspired by emission reduction practices and governmental carbon regulations, we study in this paper a two-layer supply chain consisting of a single supplier and multi-retailer under carbon emission constraints. The optimal output under completely decentralized decision and joint emission reduction decision are given, respectively. At the same time, the carbon emission reduction effects under the two decision-making situations are compared and analyzed. The results show that sharing a carbon emission quota can effectively make full use of carbon emission resources, minimize the overall cost of the supply chain, and satisfy carbon emission constraints. This paper also refers to the condition of joint decision, that is, the effective range of the carbon emission constraint coefficient. The research of this paper has certain enlightening significance for enterprise operation decisions in a low-carbon environment. Enterprises can realize effective utilization of carbon emission resources using horizontal alliance. One limitation is that we formulate the focused problem in a deterministic setting. Therefore, studying joint production planning with caps sharing facing stochastic demand can be a future research direction.

Author Contributions

Writing—original draft, C.H.; Writing—review & editing, J.B. and K.-K.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no potential conflict of interest.

Appendix A

Proof of Proposition 1. 
  • Consider
min Q i 0   A C ( Q i , , Q n )   =   i = 1 N [ a i d i Q i   +   h i p i d i 2 p i Q i   +   c i d i   +   τ ( a ^ i d i Q i + h ^ i p i d i 2 p i Q i + c ^ i d i C G H G i ) ] .
Its Hessian matrix can be shown as follows:
A C = i = 1 N [ a i d i Q i + h i p i d i 2 p i Q i + c i d i + τ ( a ^ i d i Q i + h ^ i p i d i 2 p i Q i + c ^ i d i C G H G i ) ] = i = 1 N [ ( a i + τ a ^ i ) d i Q i + ( h i + τ h ^ i ) p i d i 2 p i Q i + c i d i + τ ( c ^ i d i C G H G i ) ] H A C = [ 2 A C / Q 1 2 , 2 A C / Q 1 Q 2 , , 2 A C / Q 1 Q n 2 A C / Q 1 Q 2 , 2 A C / Q 2 2 , , 2 A C / Q 2 Q n 2 A C / Q n Q 1 , 2 A C / Q n Q 2 , , 2 A C / Q n 2 ] = [ 2 ( a 1 + τ a ^ 1 ) d 1 Q 1 3 , 0 , , 0 0 , 2 ( a 2 + τ a ^ 2 ) d 2 Q 2 3 , , 0 0 , 0 , , 2 ( a n + τ a ^ n ) d n Q n 3 ]
which implies the existence of a global optimizer. Hence, we can obtain the optimal solution as below from the first order condition of the objective function.
A C Q i = ( a i + τ a ^ i ) d i Q i 2 + ( h i + τ h ^ i ) p i d i 2 p i Q i * = a i + τ a ^ i h i + τ h ^ i 2 p i d i p i d i
Proof of Corollary 1. 
A C i , τ = 2 ( a i + τ a ^ i ) m i , τ = 2 ( a i + τ a ^ i ) ( h i + τ h ^ i ) d i ( p i d i ) 2 p i + c i d i + τ ( c ^ i d i C G H G i ) , A C i , τ τ = ( h ^ i a i + τ a ^ i h i + τ h ^ i + a ^ i h i + τ h ^ i a i + τ a ^ i ) d i ( p i d i ) 2 p i + ( c ^ i d i C G H G i ) { A C i , τ τ > 0 ,   if   c ^ i d i > C G H G i ; A C i , τ τ > 0 ,   if   c ^ i d i < C G H G i   and   [ a ^ i h i + a i h ^ i + τ ( a ^ i h ^ i + a ^ i h ^ i ) ] 2 d i ( p i d i ) 2 p i > ( C G H G i c ^ i d i ) 2 ( a i + τ a ^ i ) ( h i + τ h ^ i ) A C i , τ τ 0 ,   if   c ^ i d i < C G H G i   and   [ a ^ i h i + a i h ^ i + τ ( a ^ i h ^ i + a ^ i h ^ i ) ] 2 d i ( p i d i ) 2 p i < ( C G H G i c ^ i d i ) 2 ( a i + τ a ^ i ) ( h i + τ h ^ i ) [ a ^ i h i + a i h ^ i + τ ( a ^ i h ^ i + a ^ i h ^ i ) ] 2 d i ( p i d i ) 2 p i > ( C G H G i c ^ i d i ) 2 ( a i + τ a ^ i ) ( h i + τ h ^ i ) [ h ^ i a i + τ a ^ i h i + τ h ^ i + a ^ i h i + τ h ^ i a i + τ a ^ i ] > ( C G H G i c ^ i d i ) 2 p i d i ( p i d i )
Thus, the above inequality holds always, if
h ^ i a i + τ a ^ i h i + τ h ^ i + a ^ i h i + τ h ^ i a i + τ a ^ i 2 a ^ i h ^ i ( C G H G i c ^ i d i ) 2 p i d i ( p i d i )
and there exists two solutions if
h ^ i a i + τ a ^ i h i + τ h ^ i + a ^ i h i + τ h ^ i a i + τ a ^ i ( C G H G i c ^ i d i ) 2 p i d i ( p i d i ) 2 a ^ i h ^ i
i.e.,
a i + τ a ^ i h i + τ h ^ i [ p i ( C G H G i c ^ i d i ) 2 2 h ^ i 2 d i ( p i d i ) + p i ( C G H G i c ^ i d i ) 2 2 h ^ i 2 d i ( p i d i ) a ^ i h ^ i ] 2
and
0 < a i + τ a ^ i h i + τ h ^ i [ p i ( C G H G i c ^ i d i ) 2 2 h ^ i 2 d i ( p i d i ) p i ( C G H G i c ^ i d i ) 2 2 h ^ i 2 d i ( p i d i ) a ^ i h ^ i ] 2
Therefore, we have
0 τ a i h i [ p i ( C G H G i c ^ i d i ) 2 2 h ^ i 2 d i ( p i d i ) + p i ( C G H G i c ^ i d i ) 2 2 h ^ i 2 d i ( p i d i ) a ^ i h ^ i ] 2 h ^ i [ p i ( C G H G i c ^ i d i ) 2 2 h ^ i 2 d i ( p i d i ) + p i ( C G H G i c ^ i d i ) 2 2 h ^ i 2 d i ( p i d i ) a ^ i h ^ i ] 2 a ^ i
and
τ a i h i [ p i ( C G H G i c ^ i d i ) 2 2 h ^ i 2 d i ( p i d i ) p i ( C G H G i c ^ i d i ) 2 2 h ^ i 2 d i ( p i d i ) a ^ i h ^ i ] 2 h ^ i [ p i ( C G H G i c ^ i d i ) 2 2 h ^ i 2 d i ( p i d i ) p i ( C G H G i c ^ i d i ) 2 2 h ^ i 2 d i ( p i d i ) a ^ i h ^ i ] 2 a ^ i .
Proof of Proposition 1. 
A C i , τ = 2 ( a i + τ a ^ i ) m i , τ = 2 ( a i + τ a ^ i ) ( h i + τ h ^ i ) d i ( p i d i ) 2 p i + c i d i + τ ( c ^ i d i C G H G i ) , A C τ = i N [ ( h ^ i a i + τ a ^ i h i + τ h ^ i + a ^ i h i + τ h ^ i a i + τ a ^ i ) d i ( p i d i ) 2 p i ] i N ( C G H G i c ^ i d i ) = E ˜ ( τ ) + i N ( c ^ i d i C G H G i )   { A C τ > 0 ,   if   i N c ^ i d i > i N C G H G i ; A C τ > 0 ,   if i N c ^ i d i < i N C G H G i   and   i N E ˜ ( τ ) > i N C G H G i i N c ^ i d i ; A C τ 0 ,   if   i N c ^ i d i < i N C G H G i   and   i N E ˜ ( τ ) < i N C G H G i i N c ^ i d i . 2 A C τ 2 = E ˜ ( τ ) τ < 0 i N E ˜ ( τ ) > i N C G H G i i N c ^ i d i
Therefore, we have: (1) A C τ 0 holds if i N C G H G i i N [ 2 a ^ i h ^ i d i ( p i d i ) p i ] + i N c ^ i d i ; (2) A C τ τ < 0   if   i N C G H G i > i N [ ( h ^ i a i h i + a ^ i h i a i ) d i ( p i d i ) 2 p i ] + i N c ^ i d i ; (3) A C τ τ > 0 for 0 < τ < τ 0 and A C τ τ 0 for τ τ 0 if inf { E ˜ τ } + i N c ^ i d i < i N C G H G i < sup { E ˜ τ } + i N c ^ i d i ; (4) Systemwide emission E ^ τ always decreases in emission penalty rate τ . For emission
i N E ^ i , τ = i N [ ( a ^ i h i + τ h ^ i a i + τ a ^ i + h ^ i a i + τ a ^ i h i + τ h ^ i ) d i ( p i d i ) 2 p i + c ^ i d i ] = E ˜ τ + i N c ^ i d i a ^ i h i + τ h ^ i a i + τ a ^ i h ^ i a i + τ a ^ i h i + τ h ^ i a ^ i h ^ i ( h i + τ h ^ i ) a i + τ a ^ i a ^ i h ^ i a i h i E ^ i , τ τ = ( a i h ^ i a ^ i h i ) 2 [ a i + τ a ^ i h i + τ h ^ i a ^ i ( a i + τ a ^ i ) 2 h i + τ h ^ i a i + τ a ^ i h ^ i ( h i + τ h ^ i ) 2 ] a i + τ a ^ i h i + τ h ^ i a ^ i ( a i + τ a ^ i ) 2 > h i + τ h ^ i a i + τ a ^ i h ^ i ( h i + τ h ^ i ) 2 a ^ i ( a i + τ a ^ i ) > h ^ i ( h i + τ h ^ i ) a ^ i h i > a i h ^ i
E ^ i , τ τ < 0 always holds. □
Proof of Proposition 2. 
min Q i 0   A C ( Q i ) = a i d i Q i + h i p i d i 2 p i Q i + c i d i + τ ( a ^ i d i Q i + h ^ i p i d i 2 p i Q i + c ^ i d i C G H G i ) = ( a i + τ a ^ i ) d i Q i + ( h i + τ h ^ i ) p i d i 2 p i Q i + ( c i + τ c ^ i ) d i τ C G H G i Q i , τ = a i + τ a ^ i h i + τ h ^ i 2 p i d i p i d i = a ^ i h ^ i τ + a i a ^ i τ + h i h ^ i 2 p i d i p i d i = a ^ i h ^ i [ 1 + a i / a ^ i h i / h ^ i τ + h i / h ^ i ] 2 p i d i p i d i a i h i < a ^ i h ^ i a i a ^ i < h i h ^ i Q i , τ τ > 0 a i h i > a ^ i h ^ i a i a ^ i > h i h ^ i Q i , τ τ < 0

References

  1. He, P.; Zhang, W.; Xu, X.Y.; Bian, Y. Production lot-sizing and carbon emissions under cap-and-trade and carbon tax regulations. J. Clean. Prod. 2015, 103, 241–248. [Google Scholar] [CrossRef]
  2. Hovelaque, V.; Bironneau, L. The carbon-constrained EOQ model with carbon emission dependent demand. Int. J. Prod. Econ. 2015, 164, 285–291. [Google Scholar] [CrossRef]
  3. Bozorgi, A.; Pazour, J.; Nazzal, D. A new inventory model for cold items that considers costs and emissions. Int. J. Prod. Econ. 2014, 155, 114–125. [Google Scholar] [CrossRef]
  4. Taleizadeh, A.A.; Soleymanfar, V.R.; Govindan, K. Sustainable economic production quantity models for inventory systems with shortage. J. Clean. Prod. 2018, 174, 1011–1020. [Google Scholar] [CrossRef]
  5. Liao, H.; Deng, Q. A carbon-constrained EOQ model with uncertain demand for remanufactured products. J. Clean. Prod. 2018, 199, 334–347. [Google Scholar] [CrossRef]
  6. Wang, S.; Ye, B. A comparison between just-in-time and economic order quantity models with carbon emissions. J. Clean. Prod. 2018, 187, 662–671. [Google Scholar] [CrossRef] [Green Version]
  7. Rabta, B. An Economic Order Quantity inventory model for a product with a circular economy indicator. Comput. Ind. Eng. 2020, 140, 106215. [Google Scholar] [CrossRef]
  8. Gharaei, A.; Almehdawe, E. Optimal sustainable order quantities for growing items. J. Clean. Prod. 2021, 307, 127216. [Google Scholar] [CrossRef]
  9. Palak, G.; Eksioglu, S.D.; Geunes, J. Analyzing the impacts of carbon regulatory mechanisms on supplier and mode selection decisions: An application to a biofuel supply chain. Int. J. Prod. Econ. 2014, 154, 198–216. [Google Scholar] [CrossRef]
  10. Dye, C.-Y.; Yang, C.-T. Sustainable trade credit and replenishment decisions with credit-linked demand under carbon emission constraints. Eur. J. Oper. Res. 2015, 244, 187–200. [Google Scholar] [CrossRef]
  11. Konur, D.; Campbell, J.F.; Monfared, S.A. Economic and environmental considerations in a stochastic inventory control model with order splitting under different delivery schedules among suppliers. Omega Int. J. Manag. Sci. 2017, 71, 46–65. [Google Scholar] [CrossRef]
  12. Tang, S.; Wang, W.; Cho, S.; Yan, H. Reducing emissions in transportation and inventory management: (R, Q) Policy with considerations of carbon reduction. Eur. J. Oper. Res. 2018, 269, 327–340. [Google Scholar] [CrossRef]
  13. Yu, C.; Qu, Z.; Archibald, T.W.; Luan, Z. An inventory model of a deteriorating product considering carbon emissions. Comput. Ind. Eng. 2020, 148, 106694. [Google Scholar] [CrossRef]
  14. Mishra, U.; Wu, J.-Z.; Sarkar, B. A sustainable production-inventory model for a controllable carbon emissions rate under shortages. J. Clean. Prod. 2020, 256, 120268. [Google Scholar] [CrossRef]
  15. Ruidas, S.; Seikh, M.R.; Nayak, P.K. A production inventory model with interval-valued carbon emission parameters under price-sensitive demand. Comput. Ind. Eng. 2021, 154, 107154. [Google Scholar] [CrossRef]
  16. As’ad, R.; Hariga, M.; Shamayleh, A. Sustainable dynamic lot sizing models for cold products under carbon cap policy. Comput. Ind. Eng. 2020, 149, 106800. [Google Scholar] [CrossRef]
  17. Turken, N.; Geda, A.; Takasi, V.D.G. The impact of co-location in emissions regulation clusters on traditional and vendor managed supply chain inventory decisions. Ann. Oper. Res. 2021, 1–50. [Google Scholar] [CrossRef]
  18. Toptal, A.; Ozlu, H.; Konur, D. Joint decisions on inventory replenishment and emission reduction investment under different emission regulations. Int. J. Prod. Res. 2014, 52, 243–269. [Google Scholar] [CrossRef] [Green Version]
  19. Li, Z.; Hai, J. Inventory management for one warehouse multi-retailer systems with carbon emission costs. Comput. Ind. Eng. 2019, 130, 565–574. [Google Scholar] [CrossRef]
  20. Lu, C.-J.; Yang, C.-T.; Yen, H.-F. Stackelberg game approach for sustainable production-inventory model with collaborative investment in technology for reducing carbon emissions. J. Clean. Prod. 2020, 270, 121963. [Google Scholar] [CrossRef]
  21. Lee, J.-Y. Investing in carbon emissions reduction in the EOQ model. J. Oper. Res. Soc. 2020, 71, 1289–1300. [Google Scholar] [CrossRef]
  22. Mishra, U.; Wu, J.-Z.; Tsao, Y.-C.; Tseng, M.L. Sustainable inventory system with controllable non-instantaneous deterioration and environmental emission rates. J. Clean. Prod. 2020, 244, 118807. [Google Scholar] [CrossRef]
  23. Mashud, A.H.M.; Pervin, M.; Mishra, U.; Daryanto, Y.; Tseng, M.L.; Lim, M.K. A sustainable inventory model with controllable carbon emissions in green-warehouse farms. J. Clean. Prod. 2021, 298, 126777. [Google Scholar] [CrossRef]
  24. Qi, Q.; Zhang, R.-Q.; Bai, Q. Joint decisions on emission reduction and order quantity by a risk-averse firm under cap-and-trade regulation. Comput. Ind. Eng. 2021, 162, 107783. [Google Scholar] [CrossRef]
  25. Halat, K.; Hafezalkotob, A.; Sayadi, M.K. Cooperative inventory games in multi-echelon supply chains under carbon tax policy: Vertical or horizontal? Appl. Math. Model. 2021, 99, 166–203. [Google Scholar] [CrossRef]
  26. Fu, K.; Li, Y.; Mao, H.; Miao, Z. Firms’ production and green technology strategies: The role of emission asymmetry and carbon taxes. Eur. J. Oper. Res. 2023, 305, 1100–1112. [Google Scholar] [CrossRef]
  27. Du, S.; Tang, W.; Song, M. Low-carbon production with low-carbon premium in cap-and-trade regulation. J. Clean. Prod. 2016, 134, 652–662. [Google Scholar] [CrossRef]
  28. Yenipazarli, A. Managing new and remanufactured products to mitigate environmental damage under emissions regulation. Eur. J. Oper. Res. 2016, 249, 117–130. [Google Scholar] [CrossRef]
  29. Xu, X.; Zhang, W.; He, P.; Xu, X. Production and pricing problems in make-to-order supply chain with cap-and-trade regulation. Omega-Int. J. Manag. Sci. 2017, 66, 248–257. [Google Scholar] [CrossRef]
  30. Liu, B.; Holmbom, M.; Segerstedt, A.; Chen, W. Effects of carbon emission regulations on remanufacturing decisions with limited information of demand distribution. Int. J. Prod. Res. 2015, 53, 532–548. [Google Scholar] [CrossRef]
  31. Wang, Y.; Chen, W.; Liu, B. Manufacturing/remanufacturing decisions for a capital-constrained manufacturer considering carbon emission cap and trade. J. Clean. Prod. 2017, 140, 1118–1128. [Google Scholar] [CrossRef]
  32. Konur, D.; Schaefer, B. Integrated inventory control and transportation decisions under carbon emissions regulations: LTL vs. TL carriers. Transp. Res. Part E Logist. Transp. Rev. 2014, 68, 14–38. [Google Scholar] [CrossRef]
  33. Hariga, M.; As’ad, R.; Shamayleh, A. Integrated economic and environmental models for a multi stage cold supply chain under carbon tax regulation. J. Clean. Prod. 2017, 166, 1357–1371. [Google Scholar] [CrossRef]
  34. Rout, C.; Paul, A.; Kumar, R.S.; Chakraborty, D.; Goswami, A. Integrated optimization of inventory, replenishment and vehicle routing for a sustainable supply chain under carbon emission regulations. J. Clean. Prod. 2021, 316, 128256. [Google Scholar] [CrossRef]
  35. Konur, D.; Geunes, J. Integrated districting, fleet composition, and inventory planning for a multi-retailer distribution system. Ann. Oper. Res. 2019, 273, 527–559. [Google Scholar] [CrossRef]
  36. Park, S.J.; Cachon, G.P.; Lai, G.; Seshadri, S. Supply Chain Design and Carbon Penalty: Monopoly vs. Monopolistic Competition. Prod. Oper. Manag. 2015, 24, 1494–1508. [Google Scholar]
  37. Zhao, D.; Xu, C.; Wang, Q. Differential strategies of joint emission reductions and low-carbon promotion considering competing retailers. Control. Decis. 2014, 29, 1809–1815. [Google Scholar]
  38. Zhou, W.; Huang, W. Contract designs for energy-saving product development in a monopoly. Eur. J. Oper. Res. 2016, 250, 902–913. [Google Scholar] [CrossRef]
  39. Xu, J.; Chen, Y.; Bai, Q. A two-echelon sustainable supply chain coordination under cap-and-trade regulation. J. Clean. Prod. 2016, 135, 42–56. [Google Scholar] [CrossRef]
  40. Chen, X.; Benjaafar, S.; Elomri, A. On the effectiveness of emission penalties in decentralized supply chains. Eur. J. Oper. Res. 2019, 274, 1155–1167. [Google Scholar] [CrossRef]
  41. Lee, S.; Park, S.J. Who should lead carbon emissions reductions? Upstream vs. downstream firms. Int. J. Prod. Econ. 2020, 230, 107790. [Google Scholar] [CrossRef]
  42. Gao, F.; Souza, G.C. Carbon Offsetting with Eco-Conscious Consumers. Manag. Sci. 2022, 68, 7879–7897. [Google Scholar] [CrossRef]
  43. Guo, S.; Han, J.; Wang, W. Business models for public emission reduction projects. Oper. Res. Lett. 2023, 51, 190–196. [Google Scholar] [CrossRef]
  44. Kaur, H.; Singh, S.P. Sustainable procurement and logistics for disaster resilient supply chain. Ann. Oper. Res. 2019, 283, 309–354. [Google Scholar] [CrossRef]
  45. Yenipazarli, A. Incentives for environmental research and development: Consumer preferences, competitive pressure and emissions taxation. Eur. J. Oper. Res. 2019, 276, 757–769. [Google Scholar] [CrossRef]
  46. Wang, Y.; Yu, Z.; Jin, M.; Mao, J. Decisions and coordination of retailer-led low-carbon supply chain under altruistic preference. Eur. J. Oper. Res. 2021, 293, 910–925. [Google Scholar] [CrossRef]
  47. Feng, H.; Zeng, Y.; Cai, X.; Qian, Q.; Zhou, Y. Altruistic profit allocation rules for joint replenishment with carbon cap-and-trade policy. Eur. J. Oper. Res. 2021, 290, 956–967. [Google Scholar] [CrossRef]
  48. Andriolo, A.; Battini, D.; Grubbstrom, R.W.; Persona, A.; Sgarbossa, F. A century of evolution from Harris’s basic lot size model: Survey and research agenda. Int. J. Prod. Econ. 2014, 155, 16–38. [Google Scholar] [CrossRef] [Green Version]
  49. Jaehn, F. Sustainable Operations. Eur. J. Oper. Res. 2016, 253, 243–264. [Google Scholar] [CrossRef]
  50. Chelly, A.; Nouira, I.; Frein, Y.; Hadj-Alouane, A.B. On The consideration of carbon emissions in modelling-based supply chain literature: The state of the art, relevant features and research gaps. Int. J. Prod. Res. 2019, 57, 4977–5004. [Google Scholar] [CrossRef]
  51. Zhou, P.; Wen, W. Carbon-constrained firm decisions: From business strategies to operations modeling. Eur. J. Oper. Res. 2020, 281, 1–15. [Google Scholar] [CrossRef]
  52. Chen, X.; Benjaafar, S.; Elomri, A. The carbon-constrained EOQ. Oper. Res. Lett. 2013, 41, 172–179. [Google Scholar] [CrossRef]
  53. Carbon-Trust. Carbon Footprint in Supply Chain: The Next Step for Business. 2006. Available online: https://www.carbontrust.com/our-work-and-impact/guides-reports-and-tools/carbon-footprints-in-the-supply-chain-the-next-step-for-business (accessed on 1 June 2023).
  54. Cholette, S.; Venkat, K. The Energy and Carbon Intensity of Wine Distribution: A Study of Logistical Options for Delivering Wine to Consumers. J. Clean. Prod. 2009, 17, 1401–1413. [Google Scholar] [CrossRef]
  55. Sundarakani, B.; de Souza, R.; Goh, M.; Wagner, S.M.; Manikandan, S. Modeling Carbon Footprints across the Supply Chain. Int. J. Prod. Econ. 2010, 128, 43–50. [Google Scholar] [CrossRef]
  56. Chaabane, A.; Ramudhin, A.; Paquet, M. Design of Sustainable Supply Chains under the Emission Trading Scheme. Int. J. Prod. Econ. 2012, 135, 37–49. [Google Scholar] [CrossRef]
Table 1. Classification of relevant literature.
Table 1. Classification of relevant literature.
Problem Presented in the LiteratureSources
Optimization of EOQ/inventory model considering carbon emissions[1,2,3,4,5,6,7,8,52]
Comparison of Inventory Management under different carbon policies[9,17]
The investment behavior of companies in emissions reduction and inventory management issues[18,26]
Production strategies and pricing issues in the context of carbon emissions[27,28,29,30,31]
Integrated inventory control and transportation planning issues in the context of carbon emissions[32,33,34,35]
Research exploring existing emission reduction methods or policies[36,37,38,39,40,41,42,43]
Sustainable disaster resilient supply chain management under carbon policy[44]
Impact of carbon tax policy and other policies on corporate environmental R&D[45]
The role of altruism in low carbon supply chains[46,47]
Review of literature in the field of carbon constraint[48,49,50,51]
Table 2. Summary of notation.
Table 2. Summary of notation.
NotationExplanation
p i Annual supply (production) rate of product i per year
d i Annual demand rate of product i per year
a i Single order cost per batch
h i Annual inventory holding cost per unit product
c i Variable production cost per unit product
a ^ i Carbon emissions from a single batch order
h ^ i Carbon emission caused by inventory holding product per unit
c ^ i Carbon emissions generated by production of unit product
C G H G i Firm i’s emission quota/cap
Q i Economic production quantity, decision variable
Q i * Optimal economic production quantity without carbon emission reduction constraints
Q min * Optimal economic production quantity size considering only the minimization of carbon emissions
Q ^ i * Optimal economic production quantity with consideration of carbon emission reduction constraints and cost minimization
Q ˜ i * Optimal economic production quantity size in joint operation decision
Table 3. Comparison of system decision variables from two perspectives.
Table 3. Comparison of system decision variables from two perspectives.
ItemOptimizing Operational Cost SolelyOptimizing Carbon Cost Solely
Optimal order quantity Q i , Q ^ i , min 2 a i d i p i h i ( p i d i ) 2 a ^ i d i p i h ^ i ( p i d i )
Optimal ordering frequency m i , m ^ i , min h i d i ( p i d i ) 2 a i p i h ^ i d i ( p i d i ) 2 a ^ i p i
Optimal total cost A C 2 a i h i d i ( p i d i ) p i + c i d i τ [ 2 a ^ i h ^ i d i ( p i d i ) p i   +   ( c ^ i d i C G H G i ) ]
Optimal ordering cost A C o a i h i ( p i d i ) d i 2 p i a i h ^ i d i ( p i d i ) 2 a ^ i p i
Optimal inventory cost A C h a i h i ( p i d i ) d i 2 p i h i a ^ i d i ( p i d i ) 2 h ^ i p i
Carbon emission E M i ( a ^ i h i a i   +   h ^ i a i h i ) d i ( p i d i ) 2 p i   +   c ^ i d i 2 a ^ i h ^ i d i ( 1 d i p i ) + c ^ i d i
Table 4. Comparison of decisions under completely decentralized and quota-sharing strategies.
Table 4. Comparison of decisions under completely decentralized and quota-sharing strategies.
Completely Decentralized SituationDecentralized But Shared Carbon Emission Allowance
Optimal order quantity Q ^ i = { Q i ,   if   Q i , 1 Q i Q i , 2 ; Q i , 1 , if   Q i Q i , 1 ; Q i , 2 , if   Q i Q i , 2 ; i = 1 , 2 , , n , Q ˜ i = { 2 a i p i d i h i ( p i d i ) = Q i , if   Q Γ a i + γ a ^ i h i + γ h ^ i 2 p i d i p i d i = Q i , γ , if   Q Γ ( + n \ Γ )
Optimal order frequency m ^ i = d i Q ^ i m ˜ i = d i Q ˜ i
Optimal total cost A C i ( Q i ) = { 2 a i h i d i ( p i d i ) p i + c i d i ,   if   Q i , 1 Q i Q i , 2 ; A C i ( Q i , 1 ) , if   Q i Q i , 1 ; A C i ( Q i , 2 ) , if   Q i Q i , 2 ; A C i ( Q ˜ i ) = { 2 a i h i d i ( p i d i ) p i + c i d i ,   if   Q Γ ; ( a i h i + γ h ^ i a i + γ a ^ i + h i a i + γ a ^ i h i + γ h ^ i ) ( p i d i ) d i 2 p i + c i d i ,   if   Q Γ ( + n \ Γ ) .
Optimal order cost A C o ( Q i ) = a i m i = a i h i ( p i d i ) d i 2 p i A C o ( Q ˜ i ) = { 2 a i h i d i ( p i d i ) p i + c i d i ,   if   Q Γ ; a i h i + γ h ^ i a i + γ a ^ i ( p i d i ) d i 2 p i ,   if   Q Γ ( + n \ Γ ) .
Optimal inventory cost A C h ( Q i ) = a i m i = a i h i ( p i d i ) d i 2 p i A C h ( Q ˜ i ) = { a i h i d i ( p i d i ) p i + c i d i ,   if   Q Γ ; h i a i + γ a ^ i h i + γ h ^ i ( p i d i ) d i 2 p i ,   if   Q Γ ( + n \ Γ ) .
Total carbon emission E M i ( Q i ) = a ^ i h i d i ( p i d i ) 2 a i p i + h ^ i a i d i ( p i d i ) 2 h i p i + c ^ i d i E M i ( Q ˜ i ) = { 2 a ^ i h ^ i d i ( p i d i ) p i + c ^ i d i ( a ^ i h i + γ h ^ i a i + γ a ^ i + h ^ i a i + γ a ^ i h i + γ h ^ i ) ( p i d i ) d i 2 p i + c ^ i d i
Table 5. Numerical analysis of cost and carbon emission in settings of decentralized decision making and joint emission reduction.
Table 5. Numerical analysis of cost and carbon emission in settings of decentralized decision making and joint emission reduction.
Cap vector C G H G [2.2, 3.0, 4.5][1.5, 1.8, 2.0][1.3, 1.5, 1.6][0.83, 1.27, 1.17][0.77, 1.25, 1.20][0.74, 1.24, 1.18][0.72, 1.25, 1.17]
Critical decentralized production Q , 1 = [ Q 1 , 1   Q 2 , 1   Q 3 , 1 ] [1.46, 8.62, 2.72][2.32, 19.22, 7.97][2.79, 28.59, 11.80][5.46, 51.70, 32.97][6.24, 57.59, 27.13][6.73, 61.89, 30.39][7.10, 57.59, 32.97]
Q , 2 = [ Q 1 , 2   Q 2 , 2   Q 3 , 2 ] [428.4, 725.7, 588.6][269.2, 325.55, 200.6][223.5, 218.8, 135.5][114.45, 121, 48.49][100.1, 108.62, 58.92][92.82, 101.08, 52.59][87.92, 108.62, 48.49]
Production for minimum emission Q min = [ Q min , 1 Q min , 2 Q min , 3 ] [24.99, 79.09, 39.98][24.99, 79.09, 39.98][24.99, 79.09, 39.98][24.99, 79.09, 39.98][24.99, 79.09, 39.98][24.99, 79.09, 39.98][24.99, 79.09, 39.98]
Production for minimum total cost Q * = [ Q 1   Q 2   Q 3 ] [9.65, 32.86, 21.45][9.65, 32.86, 21.45][9.65, 32.86, 21.45][9.65, 32.86, 21.45][9.65, 32.86, 21.45][9.65, 32.86, 21.45][9.65, 32.86, 21.45]
Optimal decentalized production Q ^ * = [ Q ^ 1   Q ^ 2   Q ^ 3 ] [9.65, 32.86, 21.45][9.65, 32.86, 21.45][9.65, 32.86, 21.45][9.65, 51.70, 32.97][9.65, 57.59, 27.13][9.65, 61.89, 30.39][9.65, 57.59, 32.97]
Optimal production with caps sharing Q ˜ * = [ Q ˜ 1   Q ˜ 2   Q ˜ 3 ] [9.65, 32.86, 21.45][9.65, 32.86, 21.45][9.65, 32.86, 21.45][10.00, 34.96, 22.23][10.42, 37.37, 23.14][11.05, 40.72, 24.44][11.30, 41.98, 24.94]
Total cost with minimum emission A C ( Q min ) [11.51, 21.94, 18.44][11.51, 21.94, 18.44][11.51, 21.94, 18.44][11.51, 21.94, 18.44][11.51, 21.94, 18.44][11.51, 21.94, 18.44][11.51, 21.94, 18.44]
Total cost with minimum cost A C ( Q * ) [10.02, 20.56, 17.59][10.02, 20.56, 17.59][10.02, 20.56, 17.59][10.02, 20.56, 17.59][10.02, 20.56, 17.59][10.02, 20.56, 17.59][10.02, 20.56, 17.59]
Optimal decentralized production cost A C ( Q ^ * ) [10.02, 20.56, 17.59][10.02, 20.56, 17.59][10.02, 20.56, 17.59][10.02, 20.91, 17.98][10.02, 21.10, 17.70][10.02, 21.26, 17.85][10.02, 21.10, 17.98]
Optimal caps-sharing production cost A C ( Q ˜ * ) [10.02, 20.56, 17.59][10.02, 20.56, 17.59][10.02, 20.56, 17.59][10.02, 20.57, 17.59][10.03, 20.59, 17.60][10.05, 20.64, 17.62][10.06, 20.67, 17.63]
Optimal costs comparison analysis i N [ A C i ( Q ^ * ) A C i ( Q ˜ * ) ] 0 (→)0 (→)0 (→)0.73 (↓)0.60 (↓)0.82 (↓)0.74 (↓)
Targeted minimum emission E M ( Q min ) [0.52, 1.23, 1.16][0.52, 1.23, 1.16][0.52, 1.23, 1.16][0.52, 1.23, 1.16][0.52, 1.23, 1.16][0.52, 1.23, 1.16][0.52, 1.23, 1.16]
Emission with minimum total cost E M ( Q * ) [0.63, 1.43, 1.26][0.63, 1.43, 1.26][0.63, 1.43, 1.26][0.63, 1.43, 1.26][0.63, 1.43, 1.26][0.63, 1.43, 1.26][0.63, 1.43, 1.26]
Emission of decentralized production E M ( Q ^ * ) [0.63, 1.43, 1.26][0.63, 1.43, 1.26][0.63, 1.43, 1.26][0.63, 1.27, 1.17][0.63, 1.25, 1.20][0.63, 1.24, 1.18][0.63, 1.25, 1.17]
Emission of caps-sharing production E M ( Q ˜ * ) [0.63, 1.43, 1.26][0.63, 1.43, 1.26][0.63, 1.43, 1.26][0.62, 1.40, 1.25][0.61, 1.37, 1.24][0.60, 1.34, 1.22][0.59, 1.33, 1.22]
Comparison analysis of two- scenario emissions i N [ E M i ( Q ^ * ) E M i ( Q ˜ * ) ] 0 (→)0 (→)0 (→)−0.20 (↑)−0.14 (↑)−0.11 (↑)−0.09 (↑)
Feasibility test of caps sharing production C G H G i E M i ( Q min ) [1.68, 1.77, 3.34][0.98, 0.57, 0.84][0.78, 0.27, 0.44][0.31, 0.04, 0.01][0.25, 0.02, 0.04][0.22, 0.01, 0.02][0.20, 0.02, 0.01]
i N C G H G i i N E M i ( Q min ) 6.792.391.490.360.310.250.23
i N C G H G i i N E M i ( Q i * ) 6.381.981.08−0.05−0.10−0.16−0.18
Parameters for KKT condition γ NaNNaNNaN0.471.082.082.51
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Hu, C.; Bian, J.; Lai, K.-K. Joint Optimal Production Planning in Supply Chains under Diverse Low-Carbon Policies. Mathematics 2023, 11, 3379. https://doi.org/10.3390/math11153379

AMA Style

Hu C, Bian J, Lai K-K. Joint Optimal Production Planning in Supply Chains under Diverse Low-Carbon Policies. Mathematics. 2023; 11(15):3379. https://doi.org/10.3390/math11153379

Chicago/Turabian Style

Hu, Chenglin, Junsong Bian, and Kin-Keung Lai. 2023. "Joint Optimal Production Planning in Supply Chains under Diverse Low-Carbon Policies" Mathematics 11, no. 15: 3379. https://doi.org/10.3390/math11153379

APA Style

Hu, C., Bian, J., & Lai, K. -K. (2023). Joint Optimal Production Planning in Supply Chains under Diverse Low-Carbon Policies. Mathematics, 11(15), 3379. https://doi.org/10.3390/math11153379

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop