Abstract

In this study, the author proposes a new carbon taxing policy. This proposed carbon tax has two tax components. The first component is constant, and the second component depends on the green efficiency of production. The green efficiency of production is measured by the average amount of emissions per unit production in an assessment year. The green efficiency-based tax component can be reset every year. Lesser average emission rate indicates better green efficiency. The second component of this proposed carbon tax forces the firm to improve the green efficiency of production, which results in cleaner production. The author incorporates this new carbon tax policy in a production-inventory system with a price-sensitive demand rate. A rule is provided for the implementation of this new tax. Emissions during setup, production, and storage are considered as independent random variables. The firm has the opportunity of investing in green technologies to improve green efficiency. A profit maximization policy is adopted to solve the developed model. A solution algorithm is also provided. The model is illustrated by numerical examples with randomly generated model parameters. The results of numerical examples show the environmental benefits of the proposed carbon tax.

1. Introduction

Over the past decades, global warming has received growing attention worldwide. It poses a serious threat to our planet because its consequences are devastating, far-reaching, and long-lasting. Global warming gradually damages our planet’s ecosystem, which results in the extinction of many plant and animal species. Rise of sea level, depletion of the ozone layer, increase in average, earth’s temperature, extreme weather conditions, frequent and intense storm, drought, flood, and increased clear-air turbulence are some of the catastrophic effects of global warming (https://www.carbonpricingleadership.org/who/ [1, 2]). Industrialization, urbanization, transportation, and electricity generation are some human activities, which emit a huge amount of GHGs because these activities need to burn fossil fuels. The Intergovernmental Panel on Climate Change (IPCC), which includes more than 1,300 scientists from the United States and other countries, forecasts a temperature rise of 2.5 to 10 degrees Fahrenheit over the next century (https://climate.nasa.gov/effects/). Tollefson [3] predicted, “two degrees of warming could destroy ecosystems on around 13% of the world’s land area, increasing the risk of extinction for many insects, plants, and animals. Holding warming to 1.5°C would reduce that risk by half.” There is a growing consensus that carbon emissions are a leading contributor to global warming [4]. The industrial sector is one of the major contributors of GHGs emissions. In 2014, the United States (US) Environmental Protection Agency (EPA) reported that industry was the third major contributor to GHG emissions, contributing 21% of the total GHG emissions in the US. The industrial process consists of production/manufacturing, transportation, storage, waste disposal, etc. Each of these activities emits a huge amount of GHGs. The amount of emissions varies from industry to industry. The growing environmental consciousness of the users of the end products compels the manufacturers to tend towards green products. In a manufacturing system, three main phases of emissions are: production setup/order placing, production process, and storage [1, 4]. Emission rates in these three sources are different. Efficient production and inventory control, wise managerial decisions, green technology implementation, and environmental awareness can considerably reduce carbon emissions during manufacturing activities [1]. Strict environmental legislation is also essential to mitigate emissions. Regulatory bodies of many developed and developing countries are enacting carbon policies to reduce emissions. Carbon tax and emission trading systems (also known as cap-and-trade) are the most common mechanisms implemented in many countries to mitigate emissions. Many countries adopted a carbon tax policy to reduce carbon emissions and promote energy saving. Chile, Colombia, Finland, Argentina, Sweden, Switzerland, Netherland, Norway, Alberta, UK, and Italy are some countries/states who implemented a carbon tax. Environmental legislation affects the optimal inventory decision and must be incorporated in designing a sustainable inventory model [5]. There is a growing interest among many researchers to include environmental parameters in designing their models. A number of research papers are available in the existing literature which includes sustainability in the production-inventory system and supply chain [528]. García-Alvarado et al. [18] analyzed a stochastic closed-loop inventory model with environmental constraints. They adopted a discrete-in-time periodic review policy with a finite horizon. Battini et al. [16] incorporated sustainability in a single product replenishment problem under the traditional economic order quantity (EOQ) framework. Glock et al. [11] developed an inventory model with price- and quality-sensitive demand. They considered a single product with sustainability as a quality attribute and measured in terms of the level of scraps and emission in the supply chain. Zanoni et al. [26] published a research article on the joint economic lot size model with price and environmentally sensitive demand in a vendor-buyer supply chain. Absi et al. [6, 12] developed single item, multisourcing, and carbon-constrained lot size models. Hovelaque and Bironneau [21] investigated an inventory system with price and emission-dependent demand under environmental conditions. They showed the environmental benefit of a carbon tax. He et al. [22] investigated the production lot-sizing issues of a firm under a carbon tax and cap-and-trade regulations. Dye and Yang [5] and Tsao et al. [25] developed sustainable inventory models with trade credit under various environmental regulations. It is important to note that none of these studies included green investment as one of the model parameters. A firm should invest money on green technologies for cleaner production that also reduces the burden of the carbon tax. The green investment amount is an important decision parameter, and it should be included as a model parameter, while developing a sustainable model. Some researchers realized this fact and considered the green investment as a model parameter in their models. Jiang and Klabjan [29] investigated a joint production capacity and green investment decision policy under command-and-control and market-based regulations for an emission intensive company with stochastic demand. Lou et al. [30] incorporated green technology investment as a model parameter in a supply chain model to analyze the consequences of emission trading policy. Toptal et al. [31] analyzed a retailer’s joint decision policy on inventory replenishment and green investment with emission regulation policies. Bi et al. [32] studied on governmental subsidy policy to motivate firms’ adoption of green investment when consumers are environmentally discerning. Datta [1] published a research article on a production-inventory system with defective product under a carbon tax and green investment. Zheng et al. [33] considered green technology investment in a duopoly manufacturers’ game model under the cap-and-trade policy. Recently, Datta et al. [34] published a research article with a hybrid carbon tax and emission source-based green investment. None of the above studies considered random emissions. In the real world, the amount of GHGs emitted by a firm to produce each unit of the product is not exactly the same, varies randomly. Hence, the inclusion of random emissions in a sustainable model will make the model more realistic, particularly when the emissions rate is more volatile.

In the present article, the author proposes a novel carbon taxing policy. This proposed carbon tax consists of a constant tax component and another component based on green efficiency of production. We analyze the effects of this taxing policy on the optimal decisions in the production and inventory system of a manufacturing firm. The green efficiency of production is measured by the average emission per unit production in an assessment year. This new carbon tax can be highly effective if applied to large firms, which are the main polluters. It guarantees cleaner production. In this paper, the proposed carbon tax will be termed as green efficiency-based (GE-based) carbon tax. The model considers independent random emissions released from three main emission sources (setup, production process, and warehouse/stock-holding) of the production system. The model also incorporates green investment as a decision parameter. The demand rate is assumed to be price sensitive. To solve the model, the author adopts a profit maximization policy with the selling price, green investment, and production run-time as decision parameters. To the best of the author’s knowledge, none of the existing published articles focuses on random emissions and green efficiency-based carbon tax system. The author provides a solution algorithm. Five randomly generated numerical problems are solved to illustrate the model. To justify the contribution of this study, a comparison of this study with some related published works is shown in Table 1. The rest of the paper is organized as follows. Section 2 describes the assumptions and notations used in the model. Sections 3 and 4, respectively, present the model development and the solution algorithm. In Section 5, the author proposes how this proposed carbon tax can be implemented. Numerical illustrations are presented in Section 6. Section 7 shows the environmental benefits of the proposed carbon tax. Section 8 performs an analysis to observe the effects of the new green efficiency-based tax component on the optimum solutions. In Section 9, the author presents concluding remarks and future scopes.

2. Assumptions and Notations

(i)Time horizon is infinite.(ii)Shortages are not permitted.(iii)The product is nondeteriorating.(iv)Emissions during setup, production, and holding (storage) are independent random variables.(v)Demand rate is a decreasing linear function of selling price.(vi)Carbon tax per unit emissions has two tax components: one tax component is constant and the other component depends on the green efficiency of production.: production rate, which is known and constant (in units/year).: production run-time (in year).: length of a production cycle (in year).: selling price of each unit (in $), : the demand rate (in units/year). The demand rate is a decreasing function of s. We take it as , where and .: unit cost (in $).: carrying cost per unit per year (in $).: setup cost per production cycle (in $).: total emissions for setup activities, a random variable (in ton).: emissions per unit production run-time due to machining operations, a random variable (in ton).: emissions per unit per year for holding inventory, a random variable (in ton).: probability density function of ei defined in [0, Ri]. The joint density function of is because are independent random variables.: the expected value of defined by .: the variance of defined by .: total emissions during a production cycle (in ton).: average emissions per year during a production cycle (in ton).: index representing the green efficiency of production (in ton/unit). This index is the average amount of emissions per unit production. A smaller value of index indicates cleaner production.: carbon tax per unit emission (in $). This carbon tax is the proposed carbon tax, an increasing function of green efficiency index “.” We take it in the form where .APR: average profit per year in a production cycle (in $).EAP: EAP = E (APR), expected average profit per year in a production cycle (in $).: total $ investment in green technology. This investment is a one-off investment. The green investment $I reduce the total emissions to a proportion of .: cap/budget on green investment (in $).: is at optimum.

3. The Model Development

We assume that the firm starts its production after all necessary green installations. The production starts at time t = 0 and continues until time During production run-time , the inventory level accumulates at the rate of per year and reaches the maximum level at time During the period , the inventory level depletes at the rate of D (s) per year and becomes zero at time t = T. The cycle length is T years. Let us now calculate the average emissions per year in a production cycle.

Total emission in a production cycle without investment is

Using the relation , we obtain

With the green investment, it becomes .

The average emission per unit time without investment is , and with the green investment, it becomes .

The index representing the green efficiency of production is defined by

Costs, revenue, and the average profit per year during the cycle [0, T]:

Total revenue by selling the goods = .

The average profit per year in a production cycle is

Using the relation , we obtain

Because of the presence of three random variables in APR, we require to maximize the expected average profit , defined by

On simplifying, we obtain

The following optimization problem represents the proposed system:

Let us express EAP in powers of τ as follows:where

One can observe that

In the above equation, are functions of the decision variables

The following are two important properties derived from the structural characteristics of EAP.

Property 1. For fixed s and I, the expected average profit EAP is concave in τ, and it attains its maximum value in the interval:

Proof. To prove the concavity of EAP in τ, we need to show that .
For fixed s and I, the expressions become constant, and so EAP depends on τ only. Now, differentiating EAP twice successively with respect to τ, we obtainWe obtain because and τ > 0. This proves the concavity of EAP in τ.
The value of τ, which maximizes EAP, is the positive root of .
Now, where
For any fixed values of equation (14) is a polynomial equation in τ of degree 4 with constant coefficients. We observe that there is only one change in sign in the polynomial . So, by Descartes’ rule of sign, the polynomial equation (14) has at most one positive root. We further observe that and as . Therefore, by the property of polynomial equations, the equation (14) has an odd number of positive roots. Combining these results, we conclude that the equation (14) has a unique positive root. This positive root maximizes EAP because EAP is concave in τ.
Let us now prove that the positive root of the polynomial equation (14) lies in the open intervalWe obtainSo, for But, .
Hence, the positive root lies in the open interval .
We complete the proof of Property 1.

Property 2. For fixed s and τ, the function EAP is concave in green investment I.
We obtainwhereWe observe that are functions of For any fixed values of s and τ, the values of are constant.
Now, because Hence, EAP is concave in I. This completes the proof of Property 2.
Let us now find a decision rule to find the optimum value of the green investment I for any fixed values of s and τ. We have already proved the concavity of EAP in I for fixed values of s and τ. Therefore, the optimum value of I can be obtained by solving the equation .
Nowwhere .
By the property, the above quadratic equation has exactly one positive root.
Its positive root isWe knowand . Therefore, the following rule can be used to select the optimum value of I, for fixed s and τ:It is a tedious task to prove the joint concavity of EAP in τ, s, and I because the Hessian matrix is too complex. Therefore, we develop a solution algorithm for finding the optimum solution of the developed model based on the above properties.

4. Solution Algorithm

We develop the algorithm based on two concavity properties that we discussed in Section 3. In the algorithm, we consider the selling price“s” as an input variable. For each value of “s,” we find the optimal solution by using equations (14) and (20) and applying rule (22). For this purpose, the variable “s” is discretized by taking the initial value as “0” and incrementing it by a specific step size in each step within its domain. The solution algorithm is given below.Solution algorithm:Step 1: input model parameters Also, set the value of z, the step size of s.Step 2: set , where L is a large positive number. This, -L is the initial guess of the optimum expected average profit EAP.Step 3: set .Step 4: calculate Also, set green investment I = 0.Step 5: solve equation (14) for with most recent values of and I. We can use bisection method to solve this equation taking starting interval as mentioned in Property 1.Step 6: use equations (20) and (22) to find I with most recent values of and .Step 7: repeat steps 5 and 6 until and I become stable. It converges very rapidly.Step 8: calculate the expected average profit (EAP) by equation (8).If  > , reset  = EAP, , and Step 9: where z is an increment in . For one decimal place accuracy in , the value of z can be taken as 0.1.Step 10: if , go to step 4. Else next step.Step 11: display/print the values of .Step 12: stop.

In the above algorithm, steps 4 to 8 calculate the optimum values of the decision variables , I and, EAP for most recent value of . Step 11 gives the optimum solution of the developed model.

5. Implementation of the Proposed Carbon Tax

Many countries have currently implemented a constant carbon tax system that might be increased or decreased in subsequent years. The constant carbon tax component “a” in this proposed model can be easily implemented according to the existing carbon taxing system adopted by many countries. Every carbon emitter must pay this tax irrespective of whether the emitter is a small firm or a large firm. The second tax component “b” (GE-based) would only be imposed on the large firms. This component represents an additional tax burden for these polluters since they make a significant contribution to climate change. This new carbon tax ensures cleaner production since the firm will try to reduce by investing in green technology and waste reduction. Ideally, the tax rate “b” should be variant for different industry-based firms since the energy usage and greenhouse gases (GHGs) released by a firm for producing one unit of their product depends on the type of the industry. While estimating the tax rate “b” (in $) for the firms belonging to a particular industry, we can consider the best performer (with the least ) among these firms in the previous year as the benchmark. It means that, if the government decides to impose an additional carbon tax of $x per unit emission for the best performing firm with the least in the previous year, then “b” can be taken as $(x/G). It is mandatory for every firm on which this additional tax is applicable to disclose its annual production and emission data to the regulatory body/government. The second example in the next section provides some ideas about how the government can estimate the tax component “b.”

6. Numerical Examples

In this section, we present two numerical examples. The purpose of Example 1 is to validate the model by checking the correctness of the algorithm provided in Section 4. Example 2 demonstrates how the GE-based tax component “b” can be estimated and the impact of the proposed GE-based carbon tax on the firms’ optimal solutions.

Example 1. We consider five numerical problems by randomly generated model parameters from a specified range for each parameter using uniform distributions. Table 2 displays the ranges of the parameters and their generated values. For the sake of simplicity, except , all generated numbers are rounded to the nearest integers. These problems are solved by using the proposed algorithm, and the results are displayed in Table 3. In Table 3, a “” indicates the value at optimum.

Example 2. In this example, we consider the same five problems as in Example 1 and assume that these problems represent five firms belonging to the same industry sector. They have the same carbon tax structure. To demonstrate the estimation of tax component “b,” and the impact of the GE-based carbon tax on the firms’ optimal solutions, we consider two consecutive years. In the first year, all firms will pay the carbon tax at a constant rate “a.” Each firm’s optimal policy is to maximize its expected average profit (EAP) per year. In the second year, the government will implement the GE-based carbon tax to all the firms. At the beginning of the second year, the government will announce the values of the tax components “a” and “b” for the second year to all five firms. Accordingly, the firms will adjust their optimal decisions. We consider the first year’s constant carbon tax rate as $20 per ton. Using the proposed algorithm, we find the optimal solution of all five firms and presented in Table 4.
We can observe that the least value of the green index is 0.01835 of Firm 1. In the second year, suppose the government wants to impose an additional tax of $50 per ton on this best-performing firm, if it maintains the same . We set b = 50. This gives, b = 50/ = 2724.79, which is rounded to 2725. The government announces the second year’s tax components as a = 20 and b = 2725. Now, all the firms recalculate their optimal solutions, presented in Table 5.
We observe that the green efficiency is improved in the second year due to the imposition of GE-based carbon tax.

7. Environmental Benefits of GE-Based Carbon Tax

In this section, we perform two types of analyses to justify the proposed GE-based carbon tax’s environmental benefits. First, we perform a comparative analysis between the GE-based carbon tax and the equivalent constant carbon tax on all five problems of Example 1. Here, the equivalent constant carbon tax means the effective carbon tax per unit emission of GE-based carbon tax at optimum, i.e., . Second, we consider Example 2 and compare the results between the first year and the second year. Table 6 shows the results of GE-based carbon tax and the equivalent constant carbon tax of all five problems.

From these results, the following facts are observed in all five problems.At optimum, the average emission per unit time () is lesser for GE-based carbon tax compared to an equivalent constant carbon tax.At optimum, the green index () is lesser for GE-based carbon tax compared to an equivalent constant carbon tax.At optimum, the green investment amount () is higher for GE-based carbon tax compared to an equivalent constant carbon tax.

In the second analysis, we calculate the percentage changes of the , , and due to the imposition of the GE-based tax component in the second year (Example 2). The results are presented in Table 7.

We can observe in Table 7 that, in the second year, there is a massive percentage reduction in average emissions (column 3) and significant improvement in green efficiency (column 2). Also, there is a considerable increase in green investments (column 4). The percentage reduction in expected average profit for each firm (except Firm 2) is significantly less. Due to an insufficient green budget, Firm 2 could not achieve the actual optimum solution. Without budget restriction, in the second year, Firm 2 will generate an optimum solution:  = 0.01252,  = 3.3167,  = 32307.59, and  = 3688.14. Therefore, Firm 2 will increase its profit in the second year. We further observe in Tables 3 and 4 that, the optimum selling prices () are the same in both the years, indicating that no loss of demand after the imposition of GE-based carbon tax. These analyses justify the environmental benefits of the proposed GE-based carbon tax.

8. Effects of the Green-Based Tax Component “b” on the Optimum Solutions

In this section, we perform an analysis on the five problems defined in Example 1 to observe the effects of the new tax component “b” on the optimum values , , and . We take the range of “b” as 10–100 and present the results graphically in Figures 13. Figure 1 shows the changes in optimum green efficiency () for changes in “b.” Figure 2 shows the changes in optimum average emission () for the changes in “b.” Figure 3 represents the changes in optimum green investment () for the changes in “b.”

We observe the following properties from this analysis: is strictly decreasing functions of “b.” This result indicates the better green efficiency for higher “b.” is strictly decreasing functions of “b.” This result shows the lower emissions for higher “b.” is an increasing function of “b.” This result indicates that the firm will invest more in green technologies if “b” increases.

These observations further justify the positive environmental benefits of this proposed GE-based carbon tax.

9. Concluding Remarks

A carbon tax is a powerful way to combat carbon emissions. Many countries have enacted a carbon tax to mitigate carbon emissions. This carbon tax is constant and can be modified in subsequent years. In the present study, we propose a novel green efficiency-based carbon taxing policy. This proposed carbon tax has two tax components, a constant tax component and a green efficiency-based tax component. This proposed carbon taxing policy is different from all existing carbon policies. We also suggest how this new tax can be implemented. A solution algorithm is provided. The following are some insights on the proposed GE-based carbon tax derived from this study.This carbon tax has a very significant positive effect on the environment.It has almost no effect on annual demand.It has a less negative effect on the optimal profit. However, in some situations, it can have a positive impact on the optimal profit.The proposed carbon tax ensures cleaner production and has better environmental benefits than the equivalent constant carbon tax at optimum.For fixed values of the selling price (s) and green investment (I), the profit function is concave in production run-time (τ). This property confirms the existence of a unique optimum value of production run-time for every set of values of selling price and green investment amount.For fixed values of the selling price (s) and production run-time (τ), the profit function is concave in green investment (I). This property confirms the existence of a unique optimum value of the green investment for every set of selling price and production run-time.

The results of some related models can be derived as special cases of this model. Some of the special cases are mentioned below.(a)A constant carbon tax model with tax rate “a” can be obtained by substituting b = 0. This gives, , by equation (14) and .(b)The model with constant carbon tax “a” and without green investment can be obtained by substituting b = I = 0. In this case, .(c)The model with constant demand rate D and without carbon tax can be obtained by substituting D(s) = D and a = b = I = 0. This gives, , which agrees with the corresponding basic EOQ formula with a finite replenishment rate.

The corresponding cap-and-trade policy model with carbon price $a and annual carbon quota Q units can be derived from this developed model by substituting b = 0 and adding aQ to EAP. One can extend this model by incorporating imperfect products and variable production rates. A hybrid carbon policy model would also be of great interest. This model can also be extended by incorporating inflations.

Data Availability

No data were used to support this study.

Conflicts of Interest

The authors declare that they have no conflicts of interest.