Research Article  Open Access
Carbon Footprint Management of Road Freight Transport under the Carbon Emission Trading Mechanism
Abstract
Growing concern over environmental issues has considerably increased the number of regulations and legislation that aim to curb carbon emissions. Carbon emission trading mechanism, which is one of the most effective means, has been broadly adopted by several countries. This paper presents a road truck routing problem under the carbon emission trading mechanism. By introducing a calculation method of carbon emissions that considers the load and speed of the vehicle among other factors, a road truck routing optimizing model under the cap and trade mechanism based on the Travelling Salesman Problem (TSP) is described. Compared with the classical TSP model that only considers the economic cost, this model suggests that the truck routing decision under the cap and trade mechanism is more effective in reducing carbon emissions. A modified tabu search algorithm is also proposed to obtain solutions within a reasonable amount of computation time. We theoretically and numerically examine the impacts of carbon trading, carbon cap, and carbon price on truck routing decision, carbon emissions, and total cost. From the results of numerical experiments, we derive interesting observations about how to control the total cost and reduce carbon emissions.
1. Introduction
According to an assessment report [1] released by the Intergovernmental Panel on Climate Change, greenhouse gases mainly cause global warming and extreme weather, which gravely harm the ecosystem and human security. Greenhouse gases are principally generated from carbon emissions caused by human activities, such as fossil fuel burning in the production and transport sectors. To reduce the greenhouse gas emissions, the United Nations (UN), European Union (EU), and many other countries and organizations have adopted legislation and designed a variety of mechanisms for regulating the total amount of emissions. Carbon emission trading is one of the most cost effective market mechanisms broadly adopted by the UN, EU, and a number of countries. The European Union Emission Trading System (EU ETS) is by far the largest carbon market in the world, which implements a mandatory cap and trade system in the EU member countries. The number of emissions trading systems around the world is increasing. In addition to the EU ETS, national or subnational systems are already operating in Australia, Japan, New Zealand, and the United States and are planned in Canada, China, South Korea, and Switzerland [2]. To establish a domestic carbon emission trading system, China has launched its first regional market for compulsory carbon trading in the southern city of Shenzhen. Following Shenzhen, other carbon trading pilots have been initiated in six provinces and cities (e.g., Beijing, Tianjin, Shanghai, and Hubei).
The environmental impact of freight transport operations has attracted great attention because it is a major source of carbon dioxide (CO_{2}) emissions. According to a report by Stern [3], global transport is responsible for 14% of the total greenhouse gas emissions, and threequarters of these emissions are from the road sector. In the United Kingdom, road freight transport accounts for 22% of the CO_{2} emissions from the transport sector or 6% of the CO_{2} emissions in the country [4]. The road freight transport in the United States constitutes 78% of the total transport emissions, and the percentage of total greenhouse gas emissions from the transport sector increased from 24.9% to 27.3% between 1900 and 2005 [5]. Similar figures appear in China, where road freight transport accounts for 30% of the CO_{2} emissions from the transport sector, the percentage of total transport carbon emissions has also increased in China in recent years. In the future, governments in the world are expected to implement comprehensive regulations on carbon emissions. First, to reduce carbon emissions, firms can use equipment facilities and vehicles with low energy consumption technologies. Second, they can make optimal operational decisions in the production, inventory, and transport sectors. Compared with the adoption of energy efficient technologies, the optimal operational decisions may be an effective and low cost approach to reduce carbon emissions [6]. Therefore, how does a government design a carbon emission trading mechanism to reduce carbon emissions? How do the firms consider the impact of carbon emissions on freight transport operations? These challenging and interesting research problems are the focus of this paper.
Literature that considers carbon footprint management in supply chain management is sparse. A carbon footprint pertains to the total amount of greenhouse gas emissions produced directly and indirectly by an organization, event, product, or person and is typically expressed in equivalent tons (or kg) of carbon dioxide [7]. Several studies focus on the identification, measurement, and analysis methods of carbon emissions in the supply chain. For instance, Elhedhli and Merrick [8] considered a supply chain network design problem with emission cost and modeled the relationship between CO_{2} emissions and vehicle weight by investigating a concave minimization problem. Pattara et al. [9] introduced a standard and objective methodology and a related tool for calculating carbon footprint. They applied this tool to a wine supply chain by using the life cycle assessment (LCA) methodology. Rizet et al. [10] used a generic methodology to compare the energy consumption and CO_{2} emissions across the supply chains of products, supply chains, and countries. Lee [11] adopted a case study approach to investigate carbon emissions in automobile supply chain management. Sundarakani et al. [12] proposed an analytical and finite difference method to approximate a threedimensional infinite carbon footprint model. Leonardi and Browne [13] focused on the maritime sector and developed a method for the calculation of the carbon footprint of international supply chains. In this paper, we consider the carbon footprint of road freight transport under the carbon emission trading mechanism.
A few studies incorporated carbon emission regulations into the logistics and transport operations decisions. Kim et al. [14] examined the relationship between freight transport costs and carbon dioxide emissions in a given intermodal and truckonly freight networks using a multiobjective optimization technique. Suzuki [15] developed an approach to a timeconstrained, multiplestop, truck routing problem that minimizes the fuel consumption and pollutants emission. Hoen et al. [16] investigated the effect of two regulation mechanisms (i.e., emission cost and emission constraint) on the transport mode selection decision. They suggested that policymakers should implement a constraint on freight transport emissions. Harris et al. [17] examined the impact of the traditional cost optimization approach to strategic modeling on overall logistics costs and the CO_{2} emissions by considering the supply chain structure and different freight vehicle utilization ratios. Tajik et al. [18] addressed a time window pickupdelivery pollution routing problem and presented a new mixed integer programming approach under uncertainty. Different from these studies, we investigate the effect of the carbon emission trading mechanism on road freight transport.
Several scholars attempted to study how to reduce fuel consumption and emissions in maritime, land, and air transport by optimizing speed on the routes [19–21]. For instance, Bektaş and Laporte [22] presented a pollution routing problem (PRP), which accounted not only the travel distance, but also the amount of greenhouse gas emissions, fuel, travel times, and their costs. They also investigated the tradeoffs between various parameters, such as vehicle load, speed, and total cost. Demir et al. [23] presented an adaptive large neighborhood search for the PRP. Demir et al. [24] reported the biobjective pollution routing problem with two conflicting objective functions pertaining to minimization of driving time and fuel. They also proposed an adaptive large neighborhood search algorithm to solve it. Jabali et al. [25] proposed a framework for modeling the CO_{2} emissions in a timedependent vehicle routing context and solved the model via a tabu search procedure. They suggested that limiting vehicle speeds is desired from a total cost perspective. This paper mainly examines the operations decisions of road freight transport when managing the carbon footprints of a firm under the carbon emission trading mechanism.
Studies on the operations decisions that consider the carbon emission trading mechanism are limited. Hua et al. [26] introduced an environmental inventory model under the cap and trade mechanism and derived the optimal order quantity and provided managerial insights. Klingelhöfer [27] examined the effects of emissions trading on investments in endofpipetechnologies by employing duality theory of linear programming and sensitivity analysis. Chaabane et al. [28] presented a mixed integer linear programmingbased framework for sustainable supply chain design that considered life cycle assessment principles in addition to the traditional material balance constraints at each node in the supply chain. They suggested that current legislation and Emission Trading Schemes (ETS) must be strengthened and harmonized at the global level to derive a meaningful environmental strategy. Kwon et al. [29] studied heterogeneous fixed fleet vehicle routing with carbon emission, proposed a tabu search algorithm to solve the problem, and found that we could reduce the amount of carbon emission without sacrificing the cost due to the benefit obtained from carbon trading.
To the best of our knowledge, few related studies on the road freight transport under the carbon emission trading mechanism have been found. Suggesting road transport firms to appropriately respond to carbon emission trading would be very helpful. Based on the classical Traveling Salesman Problem (TSP), we construct an environmental road freight transport model with regard to the cap and trade system. Comparing our model with the classical TSP model allows us to examine the effects of emission trading, carbon cap, and carbon price on freight routing decision, the amount of carbon emissions, and total cost. We employ the combination of theoretical and numerical analyses. We also expect to provide decision makers with managerial insights from the research observations.
The remainder of this paper is organized as follows. This section consists of the introduction and relevant literature. Section 2 describes the carbon footprint management problem of the road freight transport and formulates the environmental truck routing model under the cap and trade mechanism. In Section 3, we theoretically examine the effects of carbon emission trading on the freight operations decisions. Section 4 illustrates the modified tabu search algorithm for solving the proposed model and shows how the different components of the algorithm work together. In Section 5, we use the numerical analysis to investigate the impacts of emission trading, carbon cap, and carbon price on the freight routing decision, the amount of carbon emissions, and the total cost. Finally, Section 5 highlights the main findings and recommends directions for future research.
2. Model Formulation
In this paper, we adopt the TSP model to study a freight truck routing problem under the carbon emission trading mechanism. The carbon emission trading mechanism includes a cap and trade system, which is a marketbased approach to control pollution. In the system, firms or national governments can trade emission allowances under an overall cap (or limit). Governmental authority typically sets the cap based on the total amount of carbon emissions. The cap is allocated or sold to firms in the form of carbon permits or credits, which represents the right to emit a specific carbon. The total number of permits cannot exceed the given cap, thereby limiting the total emissions to the objective level. The firm whose carbon emissions exceed its carbon cap must buy permits in the carbon trading market. Moreover, the firm that requires fewer permits can sell its surplus carbon permits. Thus, the buyer will pay a charge for emissions, whereas the seller will be rewarded for emission reduction. The cap and trade system can reduce overall emissions by rewarding the most efficient companies and providing less efficient companies with incentives to work toward more efficiency over time, while ensuring that nationwide emission limits can be met at the lowest economic cost. The carbon trading mechanism can ensure that the selfinterest firm seeks the economic objective while considering the environmental objective.
In this paper, we employ the TSP model to analyze the routing problem of a supplier under the carbon emissions trading mechanism. We assume that a truck sent by a supplier serves a set of its retailers starting from a distribution center, visits each retailer only once, and returns to the distribution center after serving all retailers. The product demand of each retailer is known and deterministic, and the total demand of all retailers does not exceed the capacity of the truck. The objective function of a supplier is to minimize the total cost consisting of fuel cost, truck usage cost (which depends on the usage time and includes depreciation cost, drivers salary, etc.), and carbon trading costs or gains.
2.1. Notations
In the following, we provide the notations and decision variables used in this paper: , logistics distribution network; , the set of nodes, , where 0 represents the distribution center, and other nodes denote the retailers; , the set of arcs, ; , the travel distance on arc , where ; , the truck’s speed on arc , where ; , the unit cost of fuel; , the cost of truck’s usage per unit time; , the demand of retailer , where ; , the truck’s capacity; , carbon price per unit (kg); , carbon emission quotas during a planning period in the logistics distribution; , a binary decision variable; if a truck travels on arc , then , otherwise , where ; , a nonnegative decision variable; if a truck travels on arc , then the truck’s load is , otherwise ; , a real decision variable; it represents the transfer quantity of carbon emissions. , the set of real numbers.
2.2. Calculation of Carbon Emissions
The calculation of transport carbon emissions is complex and depends on a number of factors, such as modes of transportation, fuel types, total weight of product, and distance traveled [12]. We introduce a calculation method of carbon emissions for road freight transport [16]. We assume that is the amount of carbon emissions of the truck from node to node traveling through the arc , which depends on the amount of fuel consumption on the arc and fuel emissions factor . The arc carbon emissions are calculated with the following equation:The fuel emissions factor is defined as the amount of carbon emissions per unit fuel (liter). It is a truck efficiency measure that converts fuel consumption into carbon emissions and is directly related to the vehicle type and fuel type used. For a practical logistics distribution, the fuel emissions factor is commonly a constant. According to the European carbon emission calculation criteria [16], the fuel emissions factor is set to be kg/L.
A number of related studies adopted the calculation method of fuel consumption by considering only the travel distance, which does not meet the actual condition and is inaccurate as well [14, 16]. In this paper, we comprehensively consider the impacts of distance, load, speed, and many other factors on fuel consumption. The total amount of fuel consumed on the arc is calculated as follows [22]:where is an arcspecific constant associated with road gradient, acceleration, and rolling resistance, is the empty vehicle weight, and is a vehiclespecific constant relying on the air density, drag, and frontal surface area of the vehicle. In most cases, , is related to the vehicle type and set to be for the 20ton truck used in our numerical experiments.
2.3. Decision Model
After considering the fuel and carbon emissions, we obtain the objective function of the supplier . With (1) and (2), we obtain the mixed integer programming model (denoted as MT model) of the supplier for the truck routing problem under the carbon emissions trading mechanism as follows:
The supplier aims to minimize the total cost function (3) that comprises three components, namely, fuel cost, truck usage cost, and carbon trading cost. Constraints (4) and (5) ensure that each retailer would be visited exactly once. Balance of flow is described by constraint (6), in which the flow is increased by the demand of each visited retailer. Constraint (7) guarantees that the total load of the truck would not exceed its capacity. The balance of carbon emissions is imposed by constraint (8), in which the amount of carbon emissions from the road freight transport plus the transfer quantity of carbon emissions equals the carbon emission quotas allocated to the firm. The decision variables are introduced in constraint (9), where and are vectors composed of and , respectively, is binary variable, is nonnegative variable, and is a real variable. if , then the firm sells units of carbon permits to the carbon market, if , then the firm buys units of carbon permits from the carbon market, and if , then the firm does not conduct any carbon trading business.
Now, we define two benchmark models. Without considering carbon emissions, a model denoted by MC with an objective function only minimizing economic costs issubject to constraints (4)–(7).
The decision variables are and , , , .
If the supplier only considers the carbon emission costs, then the decision problem is described as follows (denoted as ME model):
The constraints and decision variables are similar to those in the MC model.
Obviously, the MC and ME models are variations of the MT model which is the most general problem. In the MT, MC, and ME model, the optimal solution is denoted as , , and , the amount of carbon emissions is denoted as , , and , and the total cost is denoted as , , and , respectively. The following conditions evidently hold: , , and .
3. Theoretical Analysis for the Impacts of Carbon Trading Mechanism
In this section, we theoretically examine the impacts of carbon trade, carbon cap, and carbon price on truck routing decisions, carbon emissions, and total cost. We explore these complex parts through numerical experiments in the next section.
Theorem 1. For the MT model, given a fixed carbon price , no matter what the carbon cap changes, the road freight routing decisions and the amount of carbon emissions remain constant.
Proof. From carbon balance constraint (8) of the MT model, we can derive the following conditions:Substitute (12) into (3) and we haveThe constraints are independent of the carbon cap , and is the constant term of the objective function. Thus, the carbon cap merely affects the total cost and exerts no impact on the road freight routing decisions and the amount of carbon emissions .
According to Theorem 1, taking into account the impacts of the variation of carbon cap on the transfer quantity of carbon emissions and total cost, we obtain the following corollaries. From (12) and (13), we can easily derive Corollary 2.
Corollary 2. For the MT model, given a fixed carbon price, then(1)if the carbon cap decreases, then the carbon transfer quantity would decrease and the total cost would increase;(2)if the carbon cap increases, then the carbon transfer quantity would increase and the total cost would decrease.
Corollary 3. For the MT model, given the existence of a threshold of carbon cap and ,(1)if , then the firm should buy units of carbon permits;(2)if , then the firm should sell units of carbon permits;(3)if , then the firm should not conduct the carbon trade.
Proof. From (12), let and . Set , , or , respectively, and the above conditions can be derived easily.
Theorem 1 shows that the carbon cap does not affect the road freight routing arrangement of the firm and is unrelated to the amount of carbon emissions. The reason is that the amount of carbon emissions merely depends on the carbon price, and we disregard the effect of carbon cap. However, this result does not imply that the cap and trade mechanism has no influence on reducing the carbon emissions. Actually, using a market trade system allows the cap and trade mechanism to control the total amount of carbon emissions in a region, country, or even the world; nevertheless, it does not guarantee that all of the firms could reduce their carbon emissions. Thus, several firms may maintain or even increase their carbon emissions. The carbon cap and trade mechanism affects the carbon price by curbing the total amount of carbon emissions and then affects the road freight distribution decision and the amount of carbon emissions of the firm. However, the carbon cap allocated to a single firm cannot affect the carbon price, and also cannot affect the firm’s routing decision and the amount of carbon emissions.
Corollary 2 provides an intuitive result. When the carbon cap decreases, the available carbon permits of the firm will decrease, and the firm has to buy carbon permits, thus undoubtedly increasing the total cost. In contrast, when the carbon cap increases, the available carbon permits of the firm will increase, and the firm will sell its extra carbon permits, thereby reducing the total cost. Corollary 3 indicates that the decision of the firm to conduct carbon trade relies on the carbon cap.
Theorem 4. In the comparison of the MT and MC models, the following results hold:(1);(2)The upper bound of is , and if , then ;(3)If the transfer quantity of carbon emissions , then .
Proof. From (13), we know, for all , there isSpecifically, let , , then (15) can be reformulated as (16) in the following:Since and , we have ; that is, .
From (13), we know that . From (16), there is ; that is, .
Notice that if , then . Thus, .
When , from (13), we obtain , so .
Theorem 4 shows that the cap and trade mechanism can effectively guide the selfinterest firm to reduce carbon emissions, which may fundamentally explain the extensive practical application of this mechanism. Theorem 4 implies that an upper bound exists for the incremental cost arising from the cap and trade mechanism, which is directly related to the carbon price and cap. When the carbon cap is more than a threshold value that equals the amount of carbon emissions of the MC model, the firm will sell its surplus carbon permits and thus reduce the total cost. When the carbon cap is less, our numerical experiments (see Table 2) indicate that the cap and trade mechanism usually will increase the total cost of the firm. Theorem 4 shows that when the firm buys the carbon permits, its total cost is bound to increase. However, our numerical experiments (see Table 2) indicate that when the firm sells carbon permits, its total cost usually decreases, and when the firm neither buys nor sells carbon permits, its total cost may remain constant or increase.
From Theorem 4, the firm would simultaneously reduce carbon emissions and total cost. Corollary 5 provided this result.
Corollary 5. When and , the firm under the cap and trade mechanism can simultaneously reduce carbon emissions and total cost.
Proof. From the proof of Theorem 4, when , there is . From the proof of Theorem 4, when , . So the results hold.
The MT, MC, and ME models are related to each other. After comparing these models, we obtain the following theorem.
Theorem 6. The following conditions hold:(1)If , then , ;(2)If , then , ;(3)If , then , ;(4).
Proof. The MT model has one more carbon balance constraint than the MC model. From the proof of Theorem 1, substituting carbon balance constraint (8) into (13), the MT model will have constraints similar to those in the MC model, and the difference between these two models is simply the objective function. From (13), if , then ; thus , .
From (15), if , then . In particular, let , ; we have , , so , .
From (10), we obtain thatIf , then , so , .
From (13) and (17), we have thatIf , then , where is constant. So , .
We can summarize that , .
From Theorem 4, we can derive the results easily.
Theorem 6( and ) are intuitive. With the carbon cap and trade mechanism, when the carbon permit is free, the firm will disregard the impacts of carbon emissions and adopt the road truck routing decisions that minimize its economic costs. Conversely, when the carbon permit is extremely expensive, to reduce the total costs, the firm will try its best to reduce carbon emissions and adopt the road truck routing decisions that minimize the amount of carbon emissions. Theorem 6 shows that when (i.e., the total economic costs only contain the fuel cost), the economic objective of the firm is consistent with its environmental objective. The reason is that the fuel consumed directly affects carbon emissions. Therefore, the MT, MC, and ME models have the same truck routing solutions. Theorem 6 implies that the ME and MC models provide lower and upper bounds for carbon emissions. The MT model includes both economic and environmental objective and induces medium carbon emissions.
4. A Tabu Search Algorithm
Tabu search (TS) is an iterated local search metaheuristic. It makes use of adaptive memory to guide the search to escape local optima. The algorithm has been successfully used for solving many classes of TSPs and VRPs [29–31]. The TS algorithm starts from an initial solution and tries to make the best possible moves to generate a new solution in the neighborhood of the current solution even if the move may worsen the objective function value. To avoid cycling over a sequence of solutions, a set of moves are temporarily forbidden in a tabu list. The aspiration criterion is used to revoke tabu status of some moves if they create a solution better than the current best solution.
We propose a tabu search (TS) procedure for the MT model and use most of the important TS concepts in the following: tabu list, tabu tenure, aspiration criteria, neighbourhood and moves, and so forth. Next, the main components of the algorithm are described.
4.1. Initial Solution Generation
The initial solution is obtained as follows. First, the tour starts with the depot and each time a customer with the greatest demand from the set of unassigned node is selected and then inserted into the tour using the GENI algorithm presented by Gendreau et al. [32]. Finally, after all nodes are assigned, the solution is improved using the 3opt heuristic [33].
4.2. Neighbourhood Structure and Tabu Rules
Neighborhood operators are applied to explore the neighborhood of the current solution and to produce set of admissible candidate solutions. We consider four different neighborhood operators, that is, Oropt, 2opt, Exchange, and Reverse which are shown in Figure 1. In our TS scheme, the best improving (or least nonimproving) feasible move is chosen after performing each cycle of neighborhood search.
The Oropt scheme involves removing one, two, or three adjacent retailers and inserting them in another position of the route. Figure 1(b) shows the adjacent retailers 2 and 3 are reinserted in another position of the route.
The 2opt is to delete two nonadjacent arcs and add another two to generate a new route. Figure 1(c) shows the arcs and are deleted while the arcs and are added.
The Exchange operator swaps two retailers of the route. It realizes the permutation between two retailers. In Figure 1(d), the retailers 2 and 5 are swapped.
The Reverse movement involves reversing the route direction. Figure 1(e) shows the directions of all the arcs are reversed.
The computational complexity of the neighborhoods Oropt, 2opt, and Exchange is while the computational complexity of Reverse is .
The tabu status of a move is set in the following way: if a node is removed from a position of the route, it cannot return to it during the next iterations, where is the tabu tenure. We adopt the random tabu tenure; that is, is randomly selected between . However, we use an aspiration criterion: the tabu status of a move is revoked if it produces a solution that is better than the ones found in the past.
4.3. Restart Strategy and Perturbation Mechanism
To avoid a solution stagnating in a local optimal area, we consider a restart strategy which chooses a new solution as current solution for proceeding the search. If the bestsofar solution is not improved within iterations, the TS algorithm chooses the solution as the current solution, makes the tabu list empty, and exploits the search experience to lead the search escaping from local optimal.
To change the search direction, we consider a perturbation mechanism to explore more search spaces. Similar to the strategies used in Cordeau and Maischberger [30], our TS algorithm randomly chooses two operators to implement from the four neighborhoods (Oropt, 2opt, Exchange, and Reverse) to diversify the search.
4.4. Termination Condition
The maximum total number of iterations for the TS algorithm is . However, if a best solution is found when the remaining number of iterations is less than , the remaining number of iterations is changed to be .
4.5. Detailed Description of the TS
The detailed procedure of the TS algorithm is listed in the following.
Step 1. Generate the initial solution . Set the initial parameters of the algorithm: limit of the total number of iterations , limit of number of iterations without improving the best solution , and the threshold of updating , set the iteration counter of the whole algorithm , set the iteration counter without improving the best solution , set the best solution , set the current solution , and empty the tabu list.
Step 2. Repeat while the stopping criterion () is not met.
Step 2.1. Generate the neighborhood of by performing the movements Oropt, 2opt, Exchange, and Reverse in sequence, choose the solution that minimizes the total cost and is not tabu or satisfies its aspiration criteria, update the tabu list, and set .
Step 2.2. If is better than the current best solution , then update the best solution, that is, , and set ; else set .
Step 2.3. If , then restart, empty tabu list, and perform two neighborhood operators randomly selected from the neighborhoods Oropt, 2opt, Exchange, and Reverse to change the current search direction.
Step 2.4. If a new best solution is found and , then set , .
Step 2.5. Update to ; that is, .
Step 3. Return the best solution .
5. Numerical Experiments
5.1. Numerical Settings
In this section, we present a series of numerical examples to illustrate the aforementioned analytical results and simultaneously investigate the impacts of carbon cap and carbon price on road truck routing decisions, the amount of carbon emissions, and total cost and provide interesting observations. Our experiments will run with the data generated as realistically as possible. The logistics distribution network is an intercity network consisting of 20 city nodes in five provinces or cities in China (i.e., Shanghai, Zhejiang, Jiangsu, Anhui, and Jiangxi provinces). Shanghai was selected as the distribution center, and the retailers were located in other cities. The relative coordinates of these 20 nodes are shown in Table 1. The distance between two nodes and is defined as the Euclidean distance; namely, .


The demand of each retailer is an integer randomly generated in interval (kg). A truck with capacity tons and empty vehicle weight tons is used to complete the deliveries of all retailers. The parameters related to the measure of the carbon emissions of the truck are set as [22] , , kg/L, CNY/kg, kg, CNY/L, and CNY/h. The speed on each arc is an integer randomly generated in interval (km/h). The delivery trucks seldom visit more than 20 retailers per tour in a practical logistics distribution [15]; hence, we limit the problem sizes to .
Based on the above data, we generate three classes of problems with the number of the retailers in each class being , 15, and 20, respectively. Each class of problem includes 60 examples generated randomly, in which the total demand of the retailers in each instance does not exceed the capacity of the truck, and the test results are calculated from the average value across these 60 examples.
We used the TS algorithm to solve the proposed mixed linear integer programming models. The TS algorithm is coded in C language and compiled on VC++8.0. All the experiments are conducted on a server with CPU P8400 1.58 GHz PC under Windows 8. In terms of preliminary experiments, the parameters chosen for the TS algorithm are summarized in the following: , , , , and .
5.2. Impacts of Carbon Trading on Road Tuck Routing Decisions
In this section, we consider the impacts of carbon trading on truck routing decisions. In the numerical analysis, we compare the results between the MT and MC models in Table 2. To examine the relationships between total cost and carbon emissions of the ME, MC, and MT models, the results of experiments for these models are shown in Table 3. The last column listed in Table 2 reports the average computational time (CPU) over 60 trials. The columns , and of Tables 2 and 3 also show the standard deviation of the total cost across the 60 examples for ME, MC, and MT models, respectively. The data in Tables 2 and 3 verify the theoretical findings.

5.3. Impacts of Carbon Cap and Carbon Price on Carbon Emissions and Total Cost
In this section, we analyze the impacts of carbon cap on the carbon emissions and total cost. To this end, for each class of problems with different numbers of nodes, we select six parameter values of the carbon cap and calculate the amount of carbon emissions and the total cost of various instances based on different carbon caps (see Figure 2). Figure 2 shows that the carbon cap does not affect the carbon emissions of the firm in the three classes of problems with , 15, and 20. The total cost will decrease whereas the carbon cap allocated to the firm will increase. Interestingly, the carbon cap is a key resource in the cap and trade mechanism, which directly affects the cost of the firm. When the firm has more carbon caps, it will save more costs and vice versa.
(a)
(b)
(c)
To examine the impacts of carbon price on the carbon emissions and total cost, we only analyze the data of the 10node problem. The total cost is related to the carbon cap; thus, we discuss the impacts of carbon price on the carbon emissions and total cost considering five scenarios of carbon caps (Figure 3): ; ; ; ; and . The carbon prices are set from 1 to 10.
(a) ;
(b) ;
(c) ;
(d) ;
(e) ;
(f) ;
Figure 3 shows that the carbon emissions follow a “laddershaped” decreasing curve when the carbon price increases. The reason is that the higher the carbon price, the higher the cost paid by the firm to emit carbon or buy carbon permits; hence, the firm will exert efforts to decrease the carbon emissions or buy less carbon permits to reduce the total cost. The appearance of the curve of carbon emissions in the shape of a ladder with a decreasing tendency is attributed to the capacity of the firm to reduce the total cost more than the carbon emissionsminimizing ones in lowcarbon price by adopting the economic costminimizing truck routing strategies. In high carbon price, saving the total cost by adopting the truck routing strategies that keep the carbon emissions as low as possible is a better approach to the firm. Meanwhile, the carbon emissions of the firm will decrease in a stepdown manner (e.g., when in this experiment) only when the carbon price surpasses a threshold; the firm makes the truck routing decisions that minimize carbon emissions, which forms the laddertype curve that declines only after the carbon price reaches a threshold.
Figure 3 illustrates five types of curves of the total cost relative to the carbon price; namely, the total cost initially increases sharply, and then the incremental value decreases; the total cost initially increases and then remains constant; the total cost initially increases and then decreases; the total cost initially remains constant and then decreases; and the total cost initially decreases, and then the decrement value increases. Interestingly, the inflection point of the total cost curve is at the threshold of the carbon price. Therefore, we can conclude that the total cost of the firm depends on the carbon cap.
Figure 3(f) reports the results of the impacts of carbon price on the carbon emissions and total cost when and . The carbon emissions curve is likewise in a decreasing ladder shape and has several thresholds of the carbon price. The curve of the total cost still has a declining trend, and the decrement value increases when the carbon price exceeds several thresholds. The implication is that, for the problems with more nodes, the curve of carbon emissions still has a significant laddertype declining trend, and the curve of the total cost depends on the carbon cap. Another interesting observation is that the conditions of simultaneously reducing the total cost and carbon emissions is and , where is the initial threshold of carbon price.
6. Conclusion
The lowcarbon development of logistics and transport has become an inevitable trend. A number of firms in the world confront an emerging problem of how to manage the carbon footprint in their operations decisions with the subsequent enactment of regulations and legislation on carbon emissions. This paper examined the impacts of the carbon trading mechanism on road truck routing decisions. We initially introduced a calculation method of carbon emissions, including the load and speed of the vehicle. We then constructed an environmental road truck routing model under the cap and trade mechanism based on the classical TSP model. This model explored the impacts of carbon trading, carbon cap, and carbon price on truck routing decisions, the amount of carbon emissions, and the total cost. This model differs from the classical TSP model in the aspect of economic costs. Next, due to the problem complexity, a tabu search algorithm was deployed together with four neighborhoods. We also performed several numerical experiments to illustrate the analytical results and presented interesting implications and managerial insights.
The results of the computational experiments on several realistic instances provided the following important insights.
Under the cap and trade mechanism, the MT model simultaneously considers the economic and environmental aspects. The truck routing optimizing decision actually seeks the balance between minimizing the economic cost and minimizing the carbon emissions. When carbon price is zero, this model aims to minimize the economic cost, and when the carbon price is extremely high, the model aims to minimize the carbon emissions.
Compared with the MC model that minimizes the economic cost, the MT model under the cap and trade mechanism can effectively induce firms to reduce carbon emissions but may increase their total cost. The MT model may also simultaneously reduce carbon emissions and total cost of the firm when the carbon cap is large.
The carbon cap and carbon price significantly affect the road truck routing decision, carbon emissions, and total cost of the firm. The carbon cap allocated to a single firm cannot change the carbon emissions of the firm, but it will affect its total cost. Several conditions exist on whether or not the firm should buy the carbon permits. We find that when the carbon cap is less than a threshold, the firm should buy carbon permits; when the carbon cap is higher than a threshold, the firm would sell carbon permits; and when the carbon cap equals the threshold, the firm would not need to conduct carbon trading. With the increasing carbon price, the carbon emissions demonstrate a “laddertype” declining curve. The carbon emissions would begin to decrease only when the carbon price is higher than a threshold, and the variations of the total costs depend on the carbon cap.
Our study provided several insights into the logistics and transport operations of a firm under the carbon emissions trading mechanism and presented an important reference for governments to develop lowcarbon and sustainable policies. Several extensions are possible for our study. For instance, in this paper, we assume that the carbon cap has no effect on the carbon price, and a single firm deals with its operations decisions under the carbon trading mechanism. Further research may investigate the impacts of carbon pricing decisions and cooperative games among multiple firms on the amount of carbon emissions. Further research may also incorporate demand uncertainty, multiple types of vehicles, and different transport modes into the lowcarbon practice.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgments
The authors are grateful to the referee for his/her constructive comments and suggestions. This work was supported in part by the National Natural Science Foundation of China (Grant nos. 71302035, 71171178, and 71002084), Zhejiang Provincial Natural Science Foundation of China (Grant nos. Z14G020001 and LY14G020001), the Fourth Excellent Yong Talents Foundation of Zhejiang Gongshang University (Grant no. QY1323), and Collaborative Innovation Center of Building of Modern Commercial Circulation System of Zhejiang Gongshang University (14SMXY17ZD).
References
 IPCC, Climate Change 2007: The Physical Science Basis, IPCC, 2007, https://www.ipcc.ch/pdf/assessmentreport/ar4/wg1/ar4wg1spm.pdf.
 International carbon market, 2014, http://ec.europa.eu/clima/policies/ets/index_en.htm.
 N. Stern, Stern Review: The Economics of Climate Change, HM Treasury, London, UK, 2006.
 Department for Transport, Transport Statistics GB:2006 Edition, Department for Transport, London, UK, 2006.
 H. Ohnishi, “Greenhouse gas reduction strategies in the transport sector: preliminary report,” Tech. Rep., OECD/ITF Joint Transport Research Centre Working Group on GHG Reduction Strategies in the Transport Sector, OECD/ITF, Paris, France, 2008, http://www.Internationaltransportforum.org/Pub/pdf/08GHG.pdf. View at: Google Scholar
 S. Benjaafar, Y. Z. Li, and M. Daskin, Carbon Footprint and the Management of Supply Chains: Insights from Simple Models, University of Minnesota, Minneapolis, Minn, USA, 2010.
 Time for change, 2014, http://timeforchange.org/whatisacarbonfootprintdefinition.
 S. Elhedhli and R. Merrick, “Green supply chain network design to reduce carbon emissions,” Transportation Research Part D: Transport and Environment, vol. 17, no. 5, pp. 370–379, 2012. View at: Publisher Site  Google Scholar
 C. Pattara, A. Raggi, and A. Cichelli, “Life cycle assessment and carbon footprint in the wine supplychain,” Environmental Management, vol. 49, no. 6, pp. 1247–1258, 2012. View at: Publisher Site  Google Scholar
 C. Rizet, M. Browne, E. Cornelis, and J. Leonardi, “Assessing carbon footprint and energy efficiency in competing supply chains: review—case studies and benchmarking,” Transportation Research Part D: Transport and Environment, vol. 17, no. 4, pp. 293–300, 2012. View at: Publisher Site  Google Scholar
 K.H. Lee, “Integrating carbon footprint into supply chain management: the case of Hyundai Motor Company (HMC) in the automobile industry,” Journal of Cleaner Production, vol. 19, no. 11, pp. 1216–1223, 2011. View at: Publisher Site  Google Scholar
 B. Sundarakani, R. de Souza, M. Goh, S. M. Wagner, and S. Manikandan, “Modeling carbon footprints across the supply chain,” International Journal of Production Economics, vol. 128, no. 1, pp. 43–50, 2010. View at: Publisher Site  Google Scholar
 J. Leonardi and M. Browne, “A method for assessing the carbon footprint of maritime freight transport: european case study and results,” International Journal of Logistics Research and Applications, vol. 13, no. 5, pp. 349–358, 2010. View at: Publisher Site  Google Scholar
 N. S. Kim, M. Janic, and B. van Wee, “Tradeoff between carbon dioxide emissions and logistics costs based on multiobjective optimization,” Transportation Research Record, vol. 2139, pp. 107–116, 2009. View at: Publisher Site  Google Scholar
 Y. Suzuki, “A new truckrouting approach for reducing fuel consumption and pollutants emission,” Transportation Research D: Transport and Environment, vol. 16, no. 1, pp. 73–77, 2011. View at: Publisher Site  Google Scholar
 K. M. R. Hoen, T. Tan, J. C. Fransoo, and G. J. van Houtum, Effect of Carbon Emission Regulations on Transport Mode Selection in Supply Chains, Eindhoven University of Technology, Eindhoven, The Netherlands, 2010, http://alexandria.tue.nl/repository/books/672727.pdf.
 I. Harris, M. Naim, A. Palmer, A. Potter, and C. Mumford, “Assessing the impact of cost optimization based on infrastructure modelling on CO_{2} emissions,” International Journal of Production Economics, vol. 131, no. 1, pp. 313–321, 2011. View at: Publisher Site  Google Scholar
 N. Tajik, R. TavakkoliMoghaddam, B. Vahdani, and S. Meysam Mousavi, “A robust optimization approach for pollution routing problem with pickup and delivery under uncertainty,” Journal of Manufacturing Systems, vol. 33, no. 2, pp. 277–286, 2014. View at: Publisher Site  Google Scholar
 P. Cariou, “Is slow steaming a sustainable means of reducing CO_{2} emissions from container shipping?” Transportation Research Part D: Transport and Environment, vol. 16, no. 3, pp. 260–264, 2011. View at: Publisher Site  Google Scholar
 J. J. Corbett, H. Wang, and J. J. Winebrake, “The effectiveness and costs of speed reductions on emissions from international shipping,” Transportation Research Part D: Transport and Environment, vol. 14, no. 8, pp. 593–598, 2009. View at: Publisher Site  Google Scholar
 K. Fagerholt, G. Laporte, and I. Norstad, “Reducing fuel emissions by optimizing speed on shipping routes,” Journal of the Operational Research Society, vol. 61, no. 3, pp. 523–529, 2010. View at: Publisher Site  Google Scholar
 T. Bektaş and G. Laporte, “The pollutionrouting problem,” Transportation Research Part B: Methodological, vol. 45, no. 8, pp. 1232–1250, 2011. View at: Publisher Site  Google Scholar
 E. Demir, T. Bektaş, and G. Laporte, “An adaptive large neighborhood search heuristic for the pollutionrouting problem,” European Journal of Operational Research, vol. 223, no. 2, pp. 346–359, 2012. View at: Publisher Site  Google Scholar  MathSciNet
 E. Demir, T. Bektaş, and G. Laporte, “The biobjective pollutionrouting problem,” European Journal of Operational Research, vol. 232, no. 3, pp. 464–478, 2014. View at: Publisher Site  Google Scholar  MathSciNet
 O. Jabali, T. Van Woensel, and A. G. de Kok, “Analysis of travel times and CO_{2} emissions in timedependent vehicle routing,” Production and Operations Management, vol. 21, no. 6, pp. 1060–1074, 2012. View at: Publisher Site  Google Scholar
 G. Hua, T. C. E. Cheng, and S. Wang, “Managing carbon footprints in inventory management,” International Journal of Production Economics, vol. 132, no. 2, pp. 178–185, 2011. View at: Publisher Site  Google Scholar
 H. E. Klingelhöfer, “Investments in EOPtechnologies and emissions tradingresults from a linear programming approach and sensitivity analysis,” European Journal of Operational Research, vol. 196, no. 1, pp. 370–383, 2009. View at: Publisher Site  Google Scholar  MathSciNet
 A. Chaabane, A. Ramudhin, and M. Paquet, “Design of sustainable supply chains under the emission trading scheme,” International Journal of Production Economics, vol. 135, no. 1, pp. 37–49, 2012. View at: Publisher Site  Google Scholar
 Y.J. Kwon, Y.J. Choi, and D.H. Lee, “Heterogeneous fixed fleet vehicle routing considering carbon emission,” Transportation Research Part D: Transport and Environment, vol. 23, pp. 81–89, 2013. View at: Publisher Site  Google Scholar
 J.F. Cordeau and M. Maischberger, “A parallel iterated tabu search heuristic for vehicle routing problems,” Computers & Operations Research, vol. 39, no. 9, pp. 2033–2050, 2012. View at: Publisher Site  Google Scholar
 P. K. Nguyen, T. G. Crainic, and M. Toulouse, “A tabu search for timedependent multizone multitrip vehicle routing problem with time windows,” European Journal of Operational Research, vol. 231, no. 1, pp. 43–56, 2013. View at: Publisher Site  Google Scholar  MathSciNet
 M. Gendreau, A. Hertz, and G. Laporte, “New insertion and postoptimization procedures for the traveling salesman problem,” Operations Research, vol. 40, no. 6, pp. 1086–1094, 1992. View at: Publisher Site  Google Scholar  MathSciNet
 S. Lin and B. W. Kernighan, “An effective heuristic algorithm for the travelingsalesman problem,” Operations Research, vol. 21, pp. 498–516, 1973. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
Copyright
Copyright © 2015 Jin Li et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.