Optimal Inventory Control Strategies for Deteriorating Items with a General Time-Varying Demand under Carbon Emission Regulations

: Climate warming nowadays has caused people to increasingly enhance public awareness about carbon emissions from industries. In the storing industry, inventory management for deteriorating items is crucial in the business competition. To slow down the deterioration and ensure the quality of products, the items are usually stored in certain temperature-controlled environment. However, a lot of carbon emissions of the inventory system are caused by these warehousing activities. In a ﬁnite panning period, this paper studies a continuous review inventory system and proposes inventory models to analyze the impacts of carbon emissions on inventory system for deteriorating items with a general time-varying demand, in which shortages are allowed, and the customer demand during shortage period is partially backlogged till the next replenishment. Under carbon emission regulations, the existence and uniqueness of the optimal solution to each model is explored and comparisons of optimality among the proposed models are given. Numerical examples and robust analysis of the models are presented to illustrate the applicability in practice. Some management insights about inventory policies and emission reduction strategies are obtained.


Introduction
Nowadays, with the globalization of service economy, enterprises often purchase, produce, distribute, maintain, and deliver goods across continents, and there has been a significant increase in the volume of worldwide trade [1]. Facing global customers, it must be quick and timely for enterprises to market products and services, which urges inventory managers to have to make a fine judgement [2]. To enhance the competitiveness of enterprises in the market, an adequate and proper planning of inventory replenishment is required for the effective coordination of the purchase, production, distribution, maintenance, and delivery of products to meet customer needs [3,4]. Thus, inventory service operation has became an important activity for enterprises in the business competition.
Inventory model, as a mean for quantitative analysis in inventory service operation, may help managers determine the optimal inventory replenishment strategy, which can be used to optimize inventory system costs associated with the storage and maintenance of products [5]. It is often assumed that the customer demand rate remains constant throughout the inventory cycle in some economic order quantity (EOQ) inventory models [6][7][8]. However, this basic assumption is often not true in some real situations, and inventory models may be more useful when the customer demand is not a constant, but changes with time [9]. When customer demand follows time distribution, there are some different ways that items are depleted from stock during the storage period. These ways are often referred to as demand patterns. Over the past few decades, many scholars have focused on inventory models with different demand patterns such as Linear demand [10,11], Quadratic demand [12], Exponential demand [13][14][15], and so on. Nevertheless, the impact of environmental factors on the inventory system has not been investigated in the aforementioned inventory problems with the time-varying demand.
In recent years, climate warming and frequent occurrences of natural disasters has made people increasingly aware of environmental issues. The Inter-governmental Panel on Climate Change (IPCC) reports that emissions of carbon dioxide from industrial and economic activities mainly give rise to global warming [16]. To slow down global warming, some countries and regional organizations, including the United Nations (UN), the European Union (EU), and China, have carried out legislation or emission regulations such as carbon cap and carbon market to reduce carbon emissions [17]. Under these emission regulations, enterprises must replan their operations to decrease carbon emissions [18]. For example, Hewlett-Packard indicated that they decreased carbon dioxide release to the air from 26.1 tonnes to 18.3 tonnes in 2010 by changing inventory operations [19]. It is also found that firms can achieve their goal of carbon emission reduction by adjusting their inventory, warehousing, and logistics operations [20]. Up to now, carbon emission issues caused by warehousing and logistic activities have received more and more attention from industry and academia [6].
In the inventory system, another realistic phenomenon is the deterioration of goods. As is well known, some items like refrigerated food, fresh food, vegetable, and fruit have a high deterioration rate during the storage period. To slow down the deterioration, these goods often need to be stored using temperature-controlled equipments [21]. It is found that the carbon footprint mainly stems from the use of cooled/refrigerated storage in the process of storing [8,22]. In the retail industry, it is reported by Walmart (http://www.walmartstores.com/sites/responsibility-report/2012/sustainableFacilities. aspx) that the refrigerants used in grocery stores accounted for a larger percentage of Walmart's greenhouse gas footprint than its truck fleet, and similar case was confirmed in the report of Tesco (http://www.tesco.com/climatechange/carbonFootprint.asp). Therefore, in contrast to nonperishable items, carbon footprint management for the deterioration of inventory is more realistic in modern service operation. However, until now, the total effect of the deterioration and the time-varying demand on the measurable environmental performance is still an unexplored issue, which is exactly the focus of our study.
Toward this end, this paper quantifies the impact of emission regulations on the inventory model for deteriorating items with the time-varying demand. A novel model with sustainability for deterioration items under carbon emission regulations is investigated in the inventory system, in which carbon emissions caused by logistics and warehousing activities are considered, customer demands are assumed to be time-varying, and partial backlogging is allowed when the inventory system is out of stock. Specifically, the following research questions are addressed.

•
How to determine the optimal replenishment strategy and emission reduction plan under carbon emission regulations? • How do enterprises make the best choice when several carbon emission regulations coexist in the market? • What is the effect of some important inventory characteristics such as deterioration, lost sales, carbon price, and carbon cap on the replenishment strategy of the inventory system?
The remainder of this paper is organized as follows. The related literature is reviewed in Section 2. The notations and assumptions related to this study are introduced in Section 3. Section 4 develops basic inventory model. Sections 5 and 6 explore the models considering emission regulations. Then, the existence and uniqueness of the solution to the model are discussed and the comparison of optimality among the proposed models are given in Sections 5-7. In Section 7, the solution procedure of the models is illustrated through numerical examples, and in Section 8, robust analysis of the optimal solution with respect to major parameters is also carried out. The last section summarizes our contributions and discusses the directions of future research.

Literature Review
Our study is related to three streams in literatures for inventory issues. The first is on the inventory operation under the low carbon environment. As global commercial operation becomes increasingly conscious of low carbon concerns, many countries and regional organizations have launched low-carbon regulations to curb carbon emissions. When making business decisions, leading companies are also motivated to reduce carbon emissions by government's low-carbon regulations such as the cap policy and the cap-and-trade policy [23]. Under the cap policy, enterprises are not allowed to exceed mandatory emission limits predetermined by the government. However, under the cap-and-trade policy, carbon emissions can be traded in the carbon market, i.e., this regulation sets a mandatory limit on emissions to all firms engaged in business activities and allows them to buy or sell carbon emission permits. Recently, especially in inventory operation, more and more scholars have focused on the effect of carbon emissions on inventory management. From the perspective of managing carbon footprints in inventory operation, Hua et al. [6] explored how firms optimize carbon footprints under the cap and trade policy, where the optimal order size is given and the effects of key parameters like the carbon cap and the carbon price on the inventory performance are discussed. Using the EOQ framework, Chen et al. [8] investigated a carbon-constrained inventory issue under carbon mechanisms including cap, carbon tax, and cap-and-offset. Toptal et al. [24] extended the EOQ model by taking into account emission reduction investments and analyzed joint decision-making of the inventory order policy and the emission reduction investment in the context of different carbon emission mechanisms. More recently, considering demand linked credit, Dye and Yang [25] developed a sustainable inventory model for deteriorating items with carbon emission constraints. In the manufacturing environment, Tang et al. [26] explored product configuration optimisation with customer satisfaction and low carbon constraints. Xu et al. [27] further discussed inventory models for deteriorating items with the inventory-dependent consumption rate and investigated the effect of carbon emission regulations on inventory performances. From the perspective of supply chain integration, Shen et al. [28] explored a production inventory issue for deteriorating items with collaborative preservation technology investment under a carbon tax policy. Additionally to providing inventory models for studying the replenishment policy of deteriorating items, this paper will also focus on two carbon emission regulations, i.e., the carbon cap policy and the cap and trade policy, as well as the comparison of optimality among the proposed models.
The second stream is on the time-varying demand pattern, as this paper will explore how to make order decisions within a finite planning period when customer demands change with time. Recently, the growth of modern commercial activities have prompted the research on inventory management to meet customer demands. Within the EOQ framework, some extensions to existing inventory models have been explored. Among them, a very interesting extension is the consideration of time-varying demands, which has attracted much attention from academia. For example, Silver and Meal [9] were the first who explored the inventory problems with time varying demands. Donaldson [10] gave an analytical solution to the inventory model with a linear demand trend. Bose et al. [11] also investigated an EOQ model with the linear time-dependent demand and provided an inventory replenishment policy considering allowing shortages and time value of money. Later, EOQ models with the power time-dependent demand was considered by [13][14][15]29]. However, customer demands cannot always increase continuously with time. For instance, the demand rate for some season or fad products increases with time initially, then keeps steady, and finally decreases over time. Considering partial backlogging and the Weibull deterioration rate, Skouri et al. [30] studied inventory models with the general ramp type demand rate under two different inventory systems including beginning with no shortages and beginning with shortages. More recently, Sarkar et al. [12] extended the model of Ghosh and Chaudhuri [31] by investigating an inventory issue for deteriorating products with the time dependent quadratic demand rate and shortages. In the case that shortages are completely backlogged within a fixed order cycle, Cheng and Wang [32] proposed the model with a constant deterioration and trapezoidal type demand. Taking into account the price-sensitive trapezoidal time varying demand and net credit, Shah et al. [33] further developed a joint supplier-buyer inventory model. In a finite time horizon, Wang et al. [34] explored an inventory policy with time-varying demand and deteriorating item under trade credit and inflation. Xu et al. [35] discussed a two-warehouse inventory decision model for a nonperishable product with a generalized trapezoidal-type time-varying demand rate. However, the impact of carbon emission regulations on the inventory system has been ignored in the above researches.
The third stream is on the deterioration inventory issues, in which inventory models for deteriorating items have been continuously extended by more and more scholars so as to reflect more realistic inventory characteristics. As declared in literatures, deterioration phenomena during the storage period are damage, decay, evaporation, spoilage, pilferage, dryness, and so on. Once the deterioration phenomena occurs in the retail inventory system, the loss of commodity utilities may result in decreased usefulness. In real life, many physical goods such as fresh fruits and vegetables, medicine, and blood banks often undergo spoilage or decay over time. In academia, Ghare and Schrader [36] made the first attempt to describe a decaying inventory system, in which the deteriorating rate is assumed to follow the exponential time-varying distribution. Covert and Philip [37] investigated an inventory issue with two-parameter Weibull distributed deterioration. Later, some scholars gave reviews of researches about deteriorating goods in different historical periods. For instance, Nahrnias [38] summarized early works on deteriorating commodities. Based on different inventory environments, such as single and multiple products, deterministic and stochastic demand, quantity discount, and time value of money, Raafat [39] cataloged literatures on perishable commodities in the 1980s. Goyal and Giri [40] provided a detailed review on the deteriorating inventory in the last decade of the 20th century, during which the main concern of scholars is to combine the deteriorating inventory management with supply chain management, information asymmetry, stochastic environment, and other cases. Recently, Bakker et al. [41] gave a up-to-date review in the field of the deteriorating inventory since 2001 by investigating system characteristics such as price discounts, single or multiple items, single or multi-echelon, one or multi-warehouse, and so on. Considering three-parameter Weibull distributed deterioration, Yang [42] discussed a two-warehouse inventory replenishment policy with partial backlogging under inflation. Under two different contracts including cost sharing and two-part tariff, Bai et al. [43] investigated a sustainable supply chain coordination issue with deteriorating items. In contrast to these works, this paper will consider a deteriorating inventory system in which shortages are allowed and the unsatisfied order will be partially backlogged under carbon emission regulations.

•
Basic inventory description As shown in Figure 1, a continuous review inventory system running in a finite panning period is considered. Items ordered are perishable in the system, where there is no repair or replenishment of the deteriorated inventory. It is assumed that the items deteriorate at a time-varying rate θ(t), where 0 < θ(t) < 1. The replenishment rate is infinite and the lead time is zero. At the beginning of the inventory cycle, the items from stock are depleted due to the customer demand and the deterioration. When the inventory level reaches zero, shortages occur in the inventory system. As some customers are unwilling to wait backlogged orders during the shortages period, the unsatisfied demand will only be backlogged partially in the next replenishment. Thus, the backlogging rate is adopted to describe this phenomenon as λ(t) = D(t) + δI(t), where D(t) denotes the customer demand, I(t) denotes the inventory level during the shortage period, δ ≥ 0, and I(t) < 0. Namely, the backlogging rate λ(t) depends on the amount of the negative inventory I(t), which implies that the more cumulative customers in the waiting line, the more scale of lost sales. This backlogging rate function is widely adopted in literatures [44].

•
Inventory system's carbon emission assumption As some theoretical and practical studies have observed, delivering and storing items in practice consume significant amounts of energy, resulting in the creation of large amounts of carbon footprints from the transportation and inventory operations [45,46]. Thus, we focus on the carbon footprint caused by logistics and warehousing activities in our model. The following model parameters, model operators, and decision variables are used throughout the paper.

•
Model Parameters: K is the fixed cost incurred by per order; p is the selling price per unit; c is the purchasing cost per unit, and p > c; h is the holding cost per unit per unit time; c 1 is the lost sales per unit; c 2 is the shortage cost per unit per unit time; T is the fixed inventory cycle, namely, the length of inventory replenishment; S is the maximum shortages per cycle; L is the amount of lost sales per cycle; K is the amount of the carbon emissions per order; h is the amount of the carbon emissions incurred by per unit holding inventory per unit time; c is the amount of the carbon emissions incurred by per unit of delivering; c p is the carbon price per unit carbon emission; E is the capacity of the carbon emissions mandated by the government; θ(t) is the time-varying deteriorating rate when the inventory system is on-hand; D(t) is the time-varying demand rate, which is a general continuous function of time, and is positive in any given time horizon.
• Model Operators: T p(·) is the operator of the total profit for inventory system in [0, T]; E(·) is the operator of the total carbon emissions associated behavior of the on-hand inventory; is the operator of the total profit without considering carbon emissions; is the operator of the total profit under the carbon cap policy; is the operator of the total profit under the carbon cap and trade policy.
• Decision Variables: Thus, considering the inventory operation performance including sale revenue, the purchase cost, the holding cost, the opportunity cost due to lost sales, the shortage cost, and ordering cost, the total profit per cycle can be written as In addition, considering the inventory system's carbon footprint including the carbon emissions incurred by delivering, warehousing, and fixed carbon emissions per order, the amount of carbon emissions per inventory cycle can be calculated as Therefore, combining Equation (6) with the above discussions, the optimization problem in the base inventory system without considering carbon emissions can be described as To solve Problem (8), calculate the first derivative of Π 1 (t 1 ), and we have dΠ 1 (t 1 ) dt. As D(t 1 ) > 0, the optimal behavior of T p(t 1 ) only depends on the auxiliary function is a strictly monotone decreasing function in any given time interval. Therefore, the following theorem can be summarized for determining the optimal solution to M1. Theorem 1. In M1, for the above basic model with total profit Π 1 (t 1 ), there exits a unique optimal solution t * 1 ∈ (0, T) such that F(t 1 ) = 0.
Proof. Through simple derivation, we can get Recalling that F(t 1 ) is a strictly monotone decreasing function with respect to t 1 , Intermediate Value Theorem implies that there exits a unique root t * 1 ∈ (0, T) satisfying F(t * 1 ) = 0. Furthermore, taking the second derivative of T p(t 1 ) with respect to t 1 , we have According to the second order sufficient condition, it can be concluded that t * 1 is the unique maximum value point of T p(t 1 ).
Theorem 1 indicates that the optimal time point t * 1 in the inventory system is unique and satisfies the equation F(t 1 ) = 0, which makes the finding of the optimal solution to M1 easier. Meanwhile, we can easily get Q (t 1 ) > 0, S (t 1 ) < 0, and L (t 1 ) < 0 from Formulas (3)- (5), which indicates that the order quantity increases as t 1 increases, whereas the maximum shortages and the amount of lost sales decrease as t 1 increases. The results are intuitively consistent. That is, in the fixed inventory cycle, the more larger the value of t 1 , the more higher on-hand inventory level. Although the shortage level keeps smaller in the stockout period, it will make the waiting customers become more loyal due to the improvement of service level.

The Inventory Model under the Carbon Cap Mechanism
In the carbon cap model, the government sets a mandatory quota E on the amount of pollution that the inventory enterprises can emit in a given period. Under this limitation, inventory managers need to determine the optimal ordering quantity for maximizing their total profit at the beginning of the inventory cycle. By combining Equations (6) and (7), the optimization problem in the inventory system can be described as follows, The following theorem is provided for determining the optimal solution to M2.
Theorem 2. In M2, let E( t 1 ) = E, and we have (1) if t 1 ≥ T, then the optimal solution of M2 is equivalent to that of M1; (2) if t 1 < T, then the optimal solution of M2 is the same as that of M1 when F( t 1 ) < 0, and t * 1 = t 1 when F( t 1 ) ≥ 0.

Proof.
(1) If T ≤ t 1 , the mandatory quota E from the government is met automatically and the carbon cap is not binding. Thus, the optimal solution of M2 is equivalent to that of M1; (2) Let t 1 < T. When F( t 1 ) < 0, following the proof of Theorem 4.1, we can get a unique root Under the carbon cap mechanism, the optimal solution to M2 is dependent on the carbon cap E and the sign of F( t 1 ) especially when t 1 < T and F( t 1 ) ≥ 0, which is different from the case for M1. The optimal behavior of the inventory system is absolutely determined by the government's carbon cap, which implies that the government can make inventory managers become more rational by imposing appropriate carbon quota. In what follows, we will discuss how does the carbon cap affect the inventory operation performances.

The Effect of the Carbon Cap
Theorem 3. Let E NC be the optimal carbon emissions without considering the carbon emission quota. For M2, both the optimal total profit and the optimal carbon emissions increase first, and then keep unchanged as the carbon cap increases.
Proof. When E ≤ E NC , the optimal solution to M2 happens at the boundary point t * 1 = t 1 , where E( t 1 ) = E. Thus, the optimal total profit and carbon emissions increase as the carbon cap E increases. When E > E NC , the optimal solution to M2 occurs at an interior point in the feasible interval (0, T). Then, the optimality in M2 is independent of the carbon cap E, i.e., the carbon cap is not binding. Thus, the optimal total profit and carbon emissions keep constant as the carbon cap increases.
Theorem 3 shows that, on one hand, the government or policy makers can achieve the effect of reducing emissions by setting an appropriate carbon quota (i.e., E ≤ E NC ), instead of failure in the total control of carbon (i.e., E > E NC ); on the other hand, under the carbon cap mechanism, inventory managers can also maximize their profits by adjusting order and storage (i.e., E ≤ E NC ). Next, a comparison between M1 and M2 is given. Theorem 4. Let Q NC and TP NC be the optimal order quantity and the optimal total profit without considering any carbon emission limits, and Q CC , E CC , and TP CC be the optimal order quantity, the optimal carbon emissions, and the optimal total profit under the carbon cap, respectively. Compared with M1, we have (1) Q NC ≥ Q CC ; (2) TP NC ≥ TP CC ; and (3) E NC ≥ E CC .
Proof. By Theorem 2, when t 1 ≥ T, or when t 1 < T and F( t 1 ) < 0, the cap remains at a higher level, and this limitation is not binding. Then, the optimal ordering policy with the cap is the same as that without considering carbon emissions, that is, Q NC = Q CC , TP NC = TP CC , and E NC = E CC < E. However, when t 1 < T and F( t 1 ) ≥ 0, Theorem 1 shows that the cap is low. Under this restriction, the optimal time point only happens at t * 1 = t 1 . Noting the monotonicity of Q(t 1 ) and TP(t 1 ), we get Q CC < Q NC , TP NC > TP CC , and E NC > E CC = E. In summary, we obtain that Q NC ≥ Q CC , TP NC ≥ TP CC , and E NC ≥ E CC .
Theorem 4 indicates that the carbon emission reduction can be achieved under the cap policy, but the total profit in M2 is less than or equal to that in M1. Thus, the inventory decision-maker may have little incentive to accept this amount of emission reduction in a mandatory manner adopted in M2.

The Inventory Model under the Cap and Trade Mechanism
The cap and trade mechanism, as one of the most effective market-dependent mechanisms, has been adopted by many countries and regions in the world. This mechanism also sets a mandatory quota E on the amount of pollution to related industries in a given period, but different from the carbon cap mechanism in Section 5, it allows the firms to buy or sell any carbon credit within this quota from the external carbon market. In this section, let X = E(t 1 ) − E, where X > 0 (or < 0) means that the inventory firm buys (or sells) the corresponding carbon credit from the market, and X = 0 means the inventory firm neither buy or sell any credit in the carbon market. Therefore, the model with the cap and trade in the inventory system can be formulated as Combining with Equations (6) and (7), the objective function in M3 can be written as To solve M3, calculate the first derivative of Π 3 (t 1 ), and then we have dΠ 3 (t 1 ) dt. Thus, the following conclusion can be derived.
Proof. As the proof is similar to that of Theorem 1, we leave the proof to the reader. Remark 1. Theorems 1 and 5 also show that the optimal inventory decision is independent of any time-varying demand type. These findings are consistent with [47], who investigated inventory models with completely backlogging from the perspective of operation cost analysis. As carbon emission constraints are not considered in [47], our result can be regarded as a generalization of his work in this perspective. Next, the effect of the cap and the carbon price on the inventory operation performance is analyzed under the cap and trade.

The Effect of the Cap and Trade
Theorem 6. Under the cap and trade mechanism, we have (1) if X < 0, then the optimal total profit Π 3 (t 1 ) increases as the unit carbon price increases; if X > 0, then the optimal total profit Π 3 (t 1 ) decreases as the unit carbon price increases; (2) the optimal carbon emissions decrease as the unit carbon price increases; (3) as the carbon cap increases, the optimal carbon emissions keep unchanged, but the optimal total profit π 3 (t 1 ) increases.
Proof. By Theorem 5, the first-order condition for maximizing the total profit Π 3 (t 1 ) is G(t 1 ) = 0. Taking the implicit derivative of two sides of the equation G(t 1 ) = 0 with respect to c p , we can obtain which implies that the monotonicity of Π 3 (t 1 ) with respect to c p lies on the sign of X. When X < 0, the optimal total profit increases as the unit carbon price increases, and vice versa.
(2) From Equation (7), we can easily get E (t 1 ) > 0. Since dt 1 dc p < 0, we have dE(t 1 ) dc p = E (t 1 ) dt 1 dc p < 0, which shows that the optimal carbon emissions decrease as the unit carbon price increases.
(3) Noting that the parameter E is not in the implicit equation G(t 1 ) = 0, we obtain dE(t 1 ) d E = 0, and similarly, we can get Thus, the optimal carbon emissions keep unchanged as the carbon cap increases, but the optimal total profit Π 3 (t 1 ) increases as the carbon cap increases.
It can be obtained from Theorem 6 that, under the cap and trade, once carbon emissions caused by the inventory system do not exceed the carbon quota, inventory decision-makers can sell these surplus carbon emissions for improving their total profit performances. In contrast with this, excessive carbon emissions beyond the quota will reduce their profit performances. Furthermore, when the carbon price is more higher in the carbon market, the inventory decision-maker can also properly cut down carbon emissions caused by warehousing and logistics, and then earn more profit by adjusting the amount of carbon emissions. In addition, note that the carbon cap does not affect the optimal decision result in this time-varying inventory system, and this is the same as those in literatures [6,20,25]. However, when the government's initial carbon cap is higher, inventory decision-makers can gain more earnings by selling more surplus quota (i.e., . Through the comparison between M1 and M3, we can obtain the following theorem. Theorem 7. Let Q CT and E CT be the optimal order quantity and the optimal carbon emissions under the cap and trade policy. Compared with M1, we have (1) Q NC > Q CT ; and (2) E NC > E CT .

Proof. First, an auxiliary function is defined as
Obviously, when λ = c and γ = h, M(t 1 ) = F(t 1 ); when λ = c + c p c and γ = h + c p h, M(t 1 ) = G(t 1 ). Furthermore, taking the implicit derivative of two sides of the equation M(t 1 ) = 0 with respect to λ, we can get dt 1 dλ < 0 and dt 1 dγ < 0, which indicates that the optimal time point t 1 is a strictly monotone decreasing function on λ and γ, respectively. Let t NC 1 and t CT 1 be the optimal time point without considering any carbon emission limits and the optimal time point under the cap and trade, respectively. By Theorem 1 and Theorem 5, we can obtain t NC 1 > t CT 1 , as c + c p c > c and h + c p h > h. Recalling that Q (t 1 ) > 0 and E (t 1 ) > 0, we have Q NC > Q CT and E NC > E CT .
It can be concluded from Theorem 7 that, different from the compulsory mechanism (i.e., the carbon cap policy), the cap and trade policy appears more modest, and it can really induce inventory decision-makers to reduce carbon emissions by lowing the on-hand inventory level or allowing more shortages. Owing to the complexity of the problem, the comparison between M2 and M3 will be discussed through numerical examples in the next section.

Numerical Examples
In this section, the proposed models are illustrated through three numerical examples and the corresponding managerial insights are obtained. The programs are run using Matlab R2013a. Specially, two demand types are considered for illustrating the proposed models in Example 1 and 2.

Example 1
Suppose that the customer demand follows exponential time-varying distribution D(t) = 250e −0.49t . The other input parameters in the inventory system are set as follows, p = 60, c = 25, In this case, the following optimal results are obtained.

Example 2
Consider the time-quadratic demand D(t) = 2t 2 + 5t + 10, and modify the values of the other parameters in Example 1. The following input parameters are adopted, p = 200, c = 20, h = 10, (1) Figure 2 illustrates the impacts of the deterioration coefficient on the optimal inventory performance. It can be seen from Figure 2 that, no matter the inventory model is M1, M2, or M3, both the optimal time point t * 1 and the optimal total profit T p * always decrease, but the optimal order quantity Q * increases as deterioration coefficient θ 0 increases. In addition, the optimal total carbon emissions E * in M1 and M3 decrease as deterioration coefficient θ 0 increases, whereas E * in M2 keeps constant first and then decreases as deterioration coefficient θ 0 increases. These results indicate that, under the deterioration inventory environment, the above optimal behaviors can reduce the carbon emissions of the inventory system to some extent, but they also lead to lower profit of the inventory system. Therefore, when the items in the inventory system are higher deteriorating, inventory decision-makers can adopt lower on-hand inventory level and more shortages to reduce the adverse impact of deterioration on the inventory system and increase inventory performance. (2) Figure 3 illustrates the impacts of the coefficient of lost sales on the optimal inventory performance. In Figure 3, for M1, M2, and M3, as the coefficient of lost sales δ increases, the optimal order quantity Q * and the optimal profit T p * decrease, but the optimal time point t * 1 increases. Meanwhile, as the coefficient of lost sale δ increases, the optimal total carbon emissions E * increase for M1 and M3, but keep constant for M2. The increase in the coefficient of lost sales means that, when shortages happen in the inventory system, more customers are unwilling to wait backlogged shortages and turn to other retailers. Thus, to please end customers, the retailer should make order adjustment by adopting a higher on-hand inventory level instead of more shortages. However, this will also lead to higher carbon emissions and higher operation cost of the inventory system.
(3) Figure 4 illustrates the impacts of carbon cap on the optimal inventory performance. It can be observed from Figure 4 that the optimal t * 1 , Q * , T p * , and E * for M2 all increase first, and then keep unchanged as the carbon cap E increases. For M3, the optimal t * 1 , Q * , and E * keep constant, but the optimal T p * increases as the carbon cap E increases, which also verifies the conclusion in Theorem 6(3) from the perspective of numerical analysis. That is, under the cap and trade, the quantity of carbon quota does not affect the optimal inventory replenishment policy and the optimal emission reduction strategy of the inventory system. However, the inventory operator may increase the profit performance of the inventory system by selling more carbon emission permits during the increase process of carbon quota. (4) Figure 5 illustrates the impacts of carbon price on the optimal inventory performance. As shown in Figure 5, as carbon price per unit c p increases, the optimal t * 1 , Q * , and E * for M3 all increase, but the optimal T p * decreases first and then increases gradually, which is consistent with Theorem 6(1). That is, when the real carbon emissions of the inventory system exceed the carbon quota set by the government, the optimal total profit decreases with the increase of carbon price. However, once the real carbon emission level is lower than the carbon quota, selling the surplus carbon cap will increase the profit of the system during the period of rising carbon prices in the carbon market. (5) In addition, from the comprehensive view, Figures 2-5 show that the optimal total profit in M3 is always higher than that in M1 or M2, and the optimal carbon emissions in M3 are lower than that in M1 or M2. This implies that, compared with M1 and M2 cases, the M3 case (namely, the cap and trade mechanism) can achieve lower carbon emissions and higher profit performances, which enriches the conclusion of Theorem 6.3 from the perspective of numerical analysis.

Sensitivity Analysis and Managerial Insights
To illustrate the application of the model, the robustness of the model and managerial insights will be discussed in this section. As the environment of the inventory system is often dynamic, errors in the model assumptions and parameters usually happen due to the uncertainties in making inventory operation decisions. Thus, it is important for inventory managers to explore the robustness of the model. Here, we examine the impact of errors in estimating model parameters through using the same data in Example 1. The measurements of impacts are Q i  Table 1 lists the robustness results of the models obtained by varying the value of each inventory parameter while maintaining the same initial values of other parameters. From Table 1, we can reach the following conclusions.

M1 M2 M3
Parameters Chang (%) (1) Although deterioration parameter θ 0 , lost sales coefficient δ, lost sales cost c 1 , and shortage cost c 2 are not easy to be estimated accurately, the effects of these parameters on the order quantity, the emissions reduction, and the total profit are very slight. For instance, in M1, for 20% under or overestimation of θ 0 , Q 1 Q 1 changes from −0.092% to 0.09%, E 1 E 1 changes from 0.104% to −0.103%, and Π 1 Π 1 changes from 0.118% to −0.117%. This shows that estimating errors in θ 0 , δ, c 1 , and c 2 may not lead to large fluctuations in the optimal solutions. Therefore, the models are robust to the parameters θ 0 , δ, c 1 , and c 2 .
(2) The total profit in the inventory system is more sensitive to the price parameter p and the cost parameter c than other parameters; the order quantity in the inventory system is more sensitive to inventory cycle T than other parameters; and the emission reduction is more sensitive to the parameters T, c, and h.
(3) It can be observed that K and K have no influence on the optimal solutions of the models, but the total profit in the inventory system is scarcely sensitive to these two parameters. It shows that the impacts of changing K and K on the optimal replenishment policy are very slight.
(4) As discussed above, it can be concluded that the proposed inventory models are basically robust.
From the above discussion on the inventory parameters affecting the inventory performance, the following managerial insights can be obtained.
Considering the high sensitivity of total annual profit to the selling price p and the purchasing cost c, it is recommended that the inventory manager should pay particular attention to the values of those inventory parameters when making inventory decisions. In addition, the emission reduction is more sensitive to the parameters T, c, and h, but those parameters can be exactly estimated by the actual data from the market. Thus, these parameters will not have much influence on the inventory replenishment policy.
Observing that the order quantity, the emissions reduction, and the total profit are more sensitive to the underestimation of T and E than the overestimation of these two parameters, it is recommended that inventory managers should adopt the optimistic estimation rather than the conservative estimation in evaluating parameters T and E when making decisions in practice.
Noting the utilization of a generalized time-varying factor makes the application scope of the models broader, it is convenient for the decision maker to handle the problems with particular circumstances such as linear demand, time-quadratic demand, exponential demand, and so on.

Conclusions and Future Research
In this paper, we develop inventory models with partial backlogging for deteriorating items, and investigate the impacts of carbon emission regulations on the inventory system. In contrast to other related literatures considering environmental regulations, the proposed models impose a general time-varying demand in the system, as the customer demands often change with time in the real deteriorating inventory system. More importantly, under emission regulations, the existence and uniqueness of the optimal solution to each model is discussed and the comparison of optimality among the proposed models is carried out. Numerical examples and robust analysis of the models are presented to illustrate the applicability in practice. Key contributions of this paper include the following.
(1) Under environmental regulations, this article incorporates the deterioration into the inventory system, and the utilization of a generalized time-varying demand makes the application scope of the models broader. For instance, if the demand is set as a specific function (e.g., exponential demand in Example 1 and time-quadratic demand in Example 2), the proposed models can be simplified as the inventory models under certain demand pattern.
(2) Following the above discussion, it can be found that inventory managers can gain more benefits from the cap and trade regulation compared with other cases. Both the theoretical and numerical results show that, under the cap and trade, the quantity of carbon quota does not affect the optimal inventory replenishment policy and the optimal emissions reduction strategy, but the inventory manager can increase the profit performance of the inventory system by selling more carbon emission permits. Furthermore, when the items in the inventory system are higher deteriorating, lower carbon emission reduction can be achieved by lower on-hand inventory level and more shortages, but this will also lead to lower profit. Once lost sales increase in the inventory system, the retailer can adopt a higher on-hand inventory level instead of more shortages. However, this will also lead to higher carbon emissions and higher operation cost of the system.
(3) Following the robust analysis, it can be observed that the impacts of the deterioration, lost sales, and shortages on the optimal solutions are very slight. The total profit in the inventory system is more sensitive to the selling price and cost per unit than other inventory parameters; the order quantity is more sensitive to the inventory cycle than other parameters; and the emission reduction is more sensitive to the inventory cycle, the carbon cap, and the carbon emissions incurred by per unit holding. However, these parameters do not have much influence on the inventory policies, and thus the proposed models are basically robust.
In view of the importance of the inferences for inventory policies and emissions reduction strategies, there are also some limits in this paper. For example, this paper only considers the deterministic inventory issues within a fixed inventory cycle. In the future, more practical situations, such as the stochastic demand and the varying inventory cycle, can be incorporated into the models. In addition, this paper only explore the impacts of the carbon cap policy and the cap-and-trade policy on the inventory system, future study may focus on the low-carbon investment and more carbon regulations like carbon tax, carbon offset policy, and so on.