The Design of a Sustainable Location-Routing-Inventory Model Considering Consumer Environmental Behavior

Abstract: Our aim is to design a sustainable supply chain (SSC) network, which takes into consideration consumer environmental behaviors (CEBs). CEBs not only affect consumers’ demand for products with low carbon emissions, they also affect their willingness to pay premium prices for products with low carbon emissions. We incorporate CEBs into the SSC network model involving location, routing and inventory. Firstly, a multi-objective optimization model comprised of both the costs and the carbon emissions of a joint location-routing-inventory model is proposed and solved, using a multi-objective particle swarm optimization (MOPSO) algorithm. Then, a revenue function including CEBs is presented on the basis of a Pareto set of the trade-off between costs and carbon emissions. A computational experiment and sensitivity analysis are conducted, employing data from the China National Petroleum Corporation (CNPC). The results clearly indicate that our research can be applied to actual supply chain operations. In addition, some practical managerial insights for enterprises are offered.


Introduction
Along with the heightened concerns over the past few decades relating to sustainable supply chains (SSC), governments, enterprises and consumers are becoming increasingly aware of the need to reduce carbon emissions.Governments have introduced a number of regulations, such as carbon taxes, cap-and-trade mechanisms and carbon constraints to mandate carbon emission reductions in SSC management [1].In addition, a few socially responsible enterprises have engaged in voluntary emission reduction programs.Companies such as BP and Nike have taken actions to reduce emissions in order to improve their public image.Wal-Mart and Tesco require their suppliers to reveal their carbon emissions on product labels, where they can be seen by consumers and society.In addition, consumers with higher levels of environmental consciousness are willing to pay a premium price for low carbon products [2][3][4][5].The demand for low carbon products has become greater and greater [6][7][8].It can be safely assumed that low carbon products will become more competitively priced in the future.Clearly, the drive for environmental improvement is increasing.
Traditionally, a supply chain network design problem focuses on minimizing the fixed and operational costs that companies directly incur.Only recently, however, have some studies started taking carbon emissions into account [9][10][11].Many studies indicate that there is a trade-off between the environment and economics in a supply chain [12][13][14].However, it is possible to significantly reduce carbon emissions without greatly increasing costs, using proper supply chain operations [15,16].In general, there is a paradox between cost and carbon emissions in SSC management.
Companies are never going to reduce their carbon emissions until factors such as cost, profits, brand awareness and consumer pressure are involved.Currently, the main drive for carbon emission reduction can be classified into two categories.The first is mandatory emission reductions, which includes features such as carbon taxes and carbon cap policies [17].This approach to carbon emissions is punitive.The alternative method is to encourage enterprises to voluntarily reduce their carbon emissions.This encouragement, in turn, can take the form of two types of motivation.One type is through policies such as carbon allowances and cap-and-trade mechanisms.The second type takes on board market considerations.For example, studies have shown that green products have the marketing potential to endow an enterprise with a good public image, which in turn can improve the relevant products' pricing structure or increase consumer demand [8,18].Looking further into the future, the effects of exploiting the marketing potential of products with low carbon emissions will increase substantially.Creating a SSC network is both a challenge and an opportunity.Presumably, the information already available to society at large has made consumers more environmentally mature, and these mature consumers would like to purchase products with smaller carbon footprints.In this study, we propose the design of a SSC network from a market-driven perspective.Specifically, the purpose of this study is to optimize the profitability of a company through CEB.We decide on the design of the SSC network after considering the number and location of warehouses, the routes from manufacturers to warehouses and from warehouses to retailers, and the inventory polices of the various facilities.Firstly, a multi-objective model is constructed to create a trade-off between cost and carbon emissions.Then, a general revenue objective factoring in CEBs is modeled, based on the relationship between cost and carbon emissions.This study allows us to achieve the best of both worlds, i.e., maximizing the profits of companies, while reducing carbon emissions as much as possible.These achievements also represent the main points of innovation in this paper.
The remainder of the paper is organized as follows: Section 2 reviews the relevant literature.Problem descriptions and assumptions are presented in Section 3. Section 4 describes the multi-objective model that creates the trade-off between cost and carbon emissions, and then constructs the general revenue objective function taking CEBs into consideration.The approach used as a solution for the model is given in Section 5. Results of the computational experiment and a sensitivity analysis are conducted in Section 6; the managerial insights are also illustrated in this section.Finally, our conclusions are presented in Section 7.

Literature Review
A key driver of any supply chain is its distribution network.This network, however, is generally also the main source of carbon emissions.The operations of a supply chain network consist of three major components, namely location, routing and inventory (LRI).However, most existing literature integrates only any two of the above, i.e., location-routing problems [19], inventory-routing problems [20], and location-inventory problems [21], as their target topics.Ahmadi-Javid and Azad [22] presented for the first time a model to simultaneously optimize location, routing and inventory decisions in a supply chain network.Ahmadi-Javid and Seddighi [23] studied a ternary integration problem that incorporated location, routing and inventory decisions in designing a multi-source distribution network.They then solved the model using a three-phase heuristic.On the whole, very few researchers have studied the ternary integration LRI problem, and fewer still have incorporated carbon emissions into an LRI problem when designing a supply chain network.This is a very important issue, which has unfortunately been largely ignored.
Numerous studies concentrate on the trade-off between the environment and the economy in supply chain management [12,24].According to the most recent papers, three types of research have been conducted and corresponding suggestions made: (i) Translate carbon emissions into cost by introducing carbon regulations, such as a carbon tax, cap-and-trade mechanisms, etc. Kroes et al. [25] investigated the relationship between a firm's environmental performance compliance and their marketing success in the context of stringent cap-and-trade regulations.Benjaafar et al. [12] presented a cost optimization model via translating carbon emissions into unit costs by carbon price.The two studies proved that there is a close relationship between economic costs and carbon emissions.Similar research was conducted by Hua et al. [24], which studied managing carbon footprints in an inventory system under a carbon emission trading mechanism.(ii) A mandatory carbon cap is used to reduce emissions.This policy specifically prohibits companies from emitting any carbon emissions in excess of their carbon cap.Diabat et al. [26] proposed a mixed-integer program model with carbon cap constraints when designing a supply chain network.A carbon-constraint economic order quantity (EOQ) model was provided to reduce emissions by properly adjusting order quantities [16].The effects of carbon-constraint measures are significant.However, it is relatively difficult to implement such policies, as they are currently unacceptable to many companies.Businesses, which are profit-driven, lack the motivation to participate in this non-profitable activity.(iii) Provide a set of Pareto solutions, which shows the trade-off between cost and carbon emissions.The advantage of this method is that it can give a set of non-dominated solutions.In addition, the decision makers can choose their preferred configuration.Wygonik and Goodchild [27] presented trade-offs between cost, service quality and the carbon emissions of an urban delivery system.Wang et al. [28] provided a bi-objective optimization model for a green supply chain network design.One of the two objectives was cost minimization; the other was to minimize carbon emissions.The Pareto results showed that the bi-objective model is an effective tool for solving this kind of problem.However, the terminal decision will be made by managers, and thus, personal preferences will inevitably be involved.
The worldwide reduction framework would involve drawing more companies into carbon reduction activities and also into assuming social responsibilities.In order to determine how to make enterprises voluntarily reduce emissions in the context of an earnings-dominated market, it is first necessary to learn how best to improve the potential motivation for corporations to reduce their carbon emissions.The use of carbon labeling is an effective means to encourage consumers to buy environmentally friendly products.There is, however, a definite need to better understand consumers' responses to eco-labels [28].Consumers' willingness to pay a premium price for products with lower carbon emissions has been shown to be increasing [4,29].Vanclay et al. [30] defined three levels of carbon labeling (from low to high) as green, yellow and black.They then found that after labeling, the black-labeled (highest carbon emission) product sales decreased by six percent, while green-labeled product sales increased by four percent, when all other conditions were basically unchanged.These results imply that the potential effectiveness of carbon labels in emission reductions is significant.However, green products usually cost more than conventional products, which in turn makes green goods more expensive [3].The key issue is whether consumers will be willing to pay a premium price for the green goods.If not, governments may have to subsidize producers who manufacture green products [5].Some studies have shown that the higher the CEBs, the higher the price consumers are willing to pay for environmentally friendly products [2].
Economic globalization and rapid high-tech development have intensified market competition to unprecedented levels.New patterns of product competition will emerge over the next few years, and the manufacture of green products as part of that competition is an irresistible trend.Conrad [3] studied the effects of consumer environmental concerns on price, choice of product and market share in the context of duopoly.Liu et al. [8] proved that, as consumers' environmental awareness increases, retailers and manufacturers with superior eco-friendly operations will benefit in the long run.A model considering the effect of environmental conscious consumers on firms' adoption of cleaner technologies showed that, as pollution intensifies, consumers play a much more positive role in the companies' environmental activities.The consumers' attitudes encourage firms to reduce carbon emissions, even in the absence of emission regulations [7].However, many studies focus on emission reduction through governmental regulations, and rarely through market forces [31].Actually, consumer response and preference for greener products, as well as market competition, combine to strongly encourage companies to adopt environmentally friendly operations.
By reviewing previous studies, we find that very little research has been conducted on LRI optimization as a means to minimize carbon emissions.Fewer still have incorporated CEBs into a revenue model.Indeed, most studies fail to properly integrate market-driven factors-in particular CEBs-and LRI operations and revenue objectives with cost-environment trade off.In this paper, we make the following contributions: (i) The concept of consumer environmental behaviors (CEBs) was proposed and incorporated into a revenue function.CEBs not only affect consumers' demand for low carbon emission products, but also their willingness to pay a premium price for low carbon emission products.(ii) A multi-objective mixed-integer formulation for the trade-off between cost and carbon emissions was presented first.The solution was then found using the multi-objective particle swarm optimization (MOPSO).Hence, a set of distributed Pareto optimal solutions can be obtained.On this basis, revenue function can be maximized.(iii) We conduct a computational experiment based on data from the China National Petroleum Corporation (CNPC) to test the presented models.Then, the Pareto solutions are presented.In addition, a number of sensitivity analyses are implemented on multiple variables.Hence, we obtain interesting managerial insights that may be of use to logistics service firms.

Problem Description
For a supply chain network consisting of manufacturers (M), warehouses (W) and retailers (R), the location of warehouses is potentially significant.In addition, each warehouse has a specific capacity level, which makes the supply chain network more realistic.The goals of our model are to choose and allocate warehouses, schedule vehicle routes and determine an inventory policy to meet retailers' demands taking into consideration CEBs.The framework of the problem is depicted in Figure 1.CEBs-and LRI operations and revenue objectives with cost-environment trade off.In this paper, we make the following contributions: (i) The concept of consumer environmental behaviors (CEBs) was proposed and incorporated into a revenue function.CEBs not only affect consumers' demand for low carbon emission products, but also their willingness to pay a premium price for low carbon emission products.(ii)A multi-objective mixed-integer formulation for the trade-off between cost and carbon emissions was presented first.The solution was then found using the multi-objective particle swarm optimization (MOPSO).Hence, a set of distributed Pareto optimal solutions can be obtained.On this basis, revenue function can be maximized.(iii) We conduct a computational experiment based on data from the China National Petroleum Corporation (CNPC) to test the presented models.Then, the Pareto solutions are presented.In addition, a number of sensitivity analyses are implemented on multiple variables.Hence, we obtain interesting managerial insights that may be of use to logistics service firms.

Problem Description
For a supply chain network consisting of manufacturers (M), warehouses (W) and retailers (R), the location of warehouses is potentially significant.In addition, each warehouse has a specific capacity level, which makes the supply chain network more realistic.The goals of our model are to choose and allocate warehouses, schedule vehicle routes and determine an inventory policy to meet retailers' demands taking into consideration CEBs.The framework of the problem is depicted in Figure 1.In Figure 1, the operations generate cost and CO2 in a supply chain network involving location, routing and inventory.θ is the green level coefficient of products, which is decided by the CO2 emissions from the LRI operations, which can in turn be calculated by Equation (14);  is the consumer environmental behaviors (CEBs), and a larger indicates that consumers are willing to pay a higher premium for greener products.CEBs can be calculated as  In Figure 1, the operations generate cost and CO 2 in a supply chain network involving location, routing and inventory.θ is the green level coefficient of products, which is decided by the CO 2 emissions from the LRI operations, which can in turn be calculated by Equation ( 14); τ is the consumer environmental behaviors (CEBs), and a larger τ indicates that consumers are willing to pay a higher premium for greener products.CEBs can be calculated as τ " ş g g τpgqβ, where τpgq is the CEBs of consumer group g, β is a correction factor of CEBs over time.We assume β ě 1, because CEBs would not decrease over time; g is the consumer group with the worst CEBs, and g is the consumer group with the best CEBs.p is the price of the product, which is decided by θ and the CEBs τ.The market demand of a product depends on p, θ and τ.Conversely, operation-induced emissions and cost will be influenced by the market, which is important, especially in a situation of oversupply.To maximize profits, supply chain enterprises will certainly endeavor to meet consumers' preferences, so as to improve their businesses' performance.

Assumptions
We assume that the consumers are under symmetric information regarding products' carbon emissions.With the preferences displayed by CEBs, we aim to find the optimal supply chain network design and operational strategy.
(i) In this paper, the CEB choices focus on the carbon emissions from the LRI, including sourcing, production and/or recovery.It is reasonable to choose supply chain services as the study object, as they represent a major source of carbon emissions.(ii) There is no difference among delivery routes, and the road conditions are nearly the same.In other words, the carbon emissions and costs are only affected by the distance travelled.(iii) Each warehouse is assumed to follow a pQ, Rq inventory policy.That is, when the inventory of a warehouse reaches the reorder point R, a fixed quantity Q is ordered from the upper stream plant.(iv) The discussed products/services are in an oversupplied market.CEBs are in positive correlation with market demand.We assume the consumer demand function is expressed as: where D 0 is the initial demand without considering CEBs or a premium for greener products, λ 1 is the market inverse demand coefficient, λ 2 is the attraction coefficient with the environmentally friendly level of products, and τ is the consumers' environmental preference for low carbon products.Obviously, the market demand is a decreasing function of price, and an increasing function of θ and τ.

A Multi-Objective Model for Cost and Carbon Emissions
There is a trade-off between cost and carbon emissions in supply chain operations.Generally, a set of optimal Pareto solutions pc x , e x q can be obtained, and particularly, the extreme values on the Pareto curve are pc, eq and pc, eq.The aim of this paper is to find the optimal solution pc x ˚, e x ˚q in the supply chain; one which will maximize profits while taking CEBs into consideration.For these operations, we should make the following decisions: (i) Location decisions-how many warehouses should be opened, and where to locate the opened warehouses.(ii) Routing decisions-how to assign the vehicle routes from manufacturers to warehouses (M-to-W) and from warehouses to Retailers (W-to-R).(iii) Inventory decisions-what is the order quantity, and how many safety stocks should be maintained?(iv) What is the most appropriate level of green to choose?Thus, the decision variables can be denoted as Q j is the order quantity of warehouse j.
The multi-objective function includes cost and carbon emissions from location, routing and inventory.First, the cost is composed using the following terms: , and safety stock is , where h o is the ordering cost, µ i is the demand by retailer i, h j is the hold cost per unit; L j is the lead time of DC j; z α is left α-percentile of standard normal random variable Z, i.e., PpZ ď z α q " α (α is the desired percentage of retailers' orders that should be satisfied); σ 2 i is the variance of demand from retailer i.
The carbon emissions are composed of the following terms: (i) Carbon emissions from facilities.The carbon emissions of a warehouse location can be denoted as ř jPJ fj y j , where fj is the carbon emissions of building warehouse j. is the M-to-W carbon emissions per distance, and t2 is the carbon emissions per distance from warehouse j (or retailer k) to retailer k 1 .(iii) Carbon emissions from inventory.The inventory emissions come from the working inventory and safety stock, which are ÿ mPpIYJq z miv ´ÿ mPpIYJq z imv " 0, @i P I, @v P V, @m P pI Y Jq R iv ´Rmv `pn r ˆzmiv q ď n r ´1, @i P I, @m P pI Y Jq, @v P V (9) y j " t0, 1u , @j P J, @n P N j (10) x p kj " t0, 1u, @k P K, @j P J x r ji " t0, 1u, @j P J, @i P I z miv " t0, 1u , @m P pI Y Jq, i P I, @v P V.
Equation ( 2) minimizes the cost of the CLRIP, where the first three terms are the fixed location cost, inventory cost, and routing cost, respectively.Equation (3) minimizes the carbon emissions.Equation (4) restricts the inventory in warehouse j to remain within its capacity N j .Equation ( 5) restricts the load of each vehicle to within its capacity V c .Equation (6) ensures one and only one vehicle serves any retailer.Equation (7) requires that each vehicle serves no more than one warehouse.The flow conservation Equation (8) states that a vehicle entering a node must also leave the node, so as to ensure the route is circular.The sub-tour elimination Equation (9) guarantees that each tour contains a warehouse, from which the tour originates and some retailers [32], where R iv is an auxiliary variable defined for retailer i for sub-tour elimination in the route of vehicle v, n r is the number of retailers.Equations ( 10)-( 13) enforce the decision variables to remain within their respective domains.

The Revenue Model Considering CEBs
This study focuses on the effects of CEBs on the task of designing a supply chain network, which includes making LRI decisions.CEBs not only affect consumers' willingness to pay premium prices for greener products, but they also affect the market demand for such products.This willingness to pay varies greatly across industries and consumer groups and also changes in intensity over time [4].If anything, carbon emissions due to logistics operations have been a concern for a considerable length of time, as these operations are a major source of emissions.We are interested in determining how to maximize earnings, as well as how to improve competition, through the influence of CEBs in three supply chain network structures which include location, routing and inventory considerations.
As non-green products have already been in circulation for many years, the general optimal decision is based on cost minimization.In this study, however, in addition to cost, we also consider carbon emissions as a benchmark.There is a terminal consumer group with an average CEB in the market.The green level θ is closely related to the carbon emissions from the supply chain.It has been proven that the carbon emissions and cost are in negative correlation, and thus, we assume that the optimal cost corresponds to a poor performance in relation to carbon emissions, and vice versa.Specifically, with an operation map, we can connect inputs rc, cs to corresponding outputs re, es.Then θ can be denoted as θ " e{e x (14) where e x is the actual carbon emission, and thus (θ ´1) is the carbon abatement ratio.p is the price of non-green products, pc, eq represents the cost and carbon emissions, and p x is the price with pc x , e x q.
As we know, the marginal cost of carbon reduction increases by degrees.The "low hanging fruit" effect also indicates that initial basic improvement is easier, but the cleanup is harder.Thus, the above situation is considered, and the price of a product with green level θ is It is worth noting that product price is a quadratic function of θ, since the environmental improvement has an increasing marginal cost, and production price is worked out to the costing.The quadratic function is commonly used to describe the cost related to the product's environmental improvement.That is, each additional increment of emissions reduction is more difficult, and hence costlier to achieve [8].Also, from the market's perspective, consumers with CEBs are willing to pay a premium price for green products.The greener the product, the more expensive it will be.In addition, for an advanced green product, too, even a small improvement will result in a significant price increase.This increase is deemed to be reasonable.
The aim is to find an optimal portfolio pc x , e x q under this context.The basic profit function can be defined as where constant c 0 is the unit cost of raw materials, and ε is other expenditure, which can be ignored in most cases.Substituting Equations ( 1)- (15) into Equation ( 16), then: Based on Equation (17), if there is no CEB, the enterprise loses the motivation to reduce carbon emissions, which is consistent with the traditional model.We assume the traditional model has revenue of Π C , with the only measure being the cost, and we mark it as model PC.In this condition, θ " 1, τ " 0, thus: Enterprises have an incentive to join in carbon reduction practices only when Π ´Πc ą 0. c 0 is the same constant in Π and Π c , and thus can be ignored.Then, enterprises will participate in carbon emissions when ppθ 2 ´cx qpD 0 ´λ1 pθ 2 `1 2 λ 2 τθq ą pp ´cqpD 0 ´λ1 pq, and thus . Actually, θ is a function of e x ; the relationship between c x and e x is important, and it will be solved in the next section.

Particle Swarm Optimization Algorithm
The particle swarm optimization (PSO) algorithm was first proposed by Kennedy and Eberhart [32].It is a population-based optimization technique and is becoming very popular, due mainly to its simplicity of implementation and ability to quickly converge to a reasonably good solution [33].It has been extensively applied to many complex network optimization problems.In the PSO heuristic procedure, a swarm of particles is retained in the search process.Each particle follows Sustainability 2016, 8, 211 9 of 20 a specific trajectory in the search space, and each step of the particle determines a trial solution.Each particle has knowledge of its previous best experience, as well as the best global experience of the entire swarm.The current best fitness of, i.e., the best solution found so far by particle p is represented by x pbest p , while the global best fitness among all particles is represented by x gbest .The velocity and position of particle p at iteration (time) t in dimension d are represented by v pd ptq and x pd ptq, respectively.Each particle updates its direction at time t according to Equation (19) in the following [32]: v pd ptq " ωv pd pt ´1q `c1 r 1 ptqpx pbest pd ´xpd ptqq `c2 r 2 ptqpx gbest ´xpd ptqq (19) where ω is the inertia influencing the local and global ability of the particle; usually a value between 0.2 to 0.6 is recommended; c 1 and c 2 are cognitive and social learning rates, respectively, and r 1 ptq and r 2 ptq are two uniform random numbers such that r 1 , r 2 P r0, 1s.The position of particle p is then updated according to Equation (20) in the following The update of velocity and the position process is repeated for every dimension and for all particles in the swarm.Eventually the swarm as a whole, like a flock of birds collectively foraging for food, is likely to move close to an optimum of the fitness function [33].

The Hybrid PSO
The multi-objective model contains location and routing assignments involving binary decisions.Multi-objective programming problems with binary variables cannot be directly processed using the Multi-objective Particle Swarm Optimization (MOPSO) heuristic procedure.Following Shankar et al. [33], the velocity of a particle should be modified if x d is binary.The modified velocity can be updated as: v pd ptq " v pd pt ´1q `r1 ptqpx pbest pd ´xpd pt ´1qq `r2 ptqpx gbest ´xpd pt ´1qq where r 1 ptq and r 2 ptq are two random numbers.The position of particle p can be updated as: where ρ pd is a uniformly distributed random number such that ρ pd P r0, 1s and spv pd ptqq is the probability threshold given by spv pd q " 1 1 `expp´v pd ptqq .In the MOPSO heuristic procedure, the velocity and positions of the continuous particles are updated according to Equations ( 19) and (20), respectively, while those of the binary variables are updated according to Equations ( 21) and ( 22), respectively.

An Improved Constraint of the MOPSO
In order to improve the ability of the heuristic procedure to search the edges crossing unconnected parts of the feasible region, and also to obtain global non-dominated solutions, some infeasible solutions that are near the feasible solutions are retained in the swarm at the beginning of the search process.A constraint that restricts the infeasibility degree of the constraints is used.At the end of the solution process, all particles retained in the swarm must be feasible.Any infeasible particles will be deleted from the external file gradually, throughout the progress of the search process.A dynamic self-adapting process is needed to control the infeasibility degree in the heuristic procedure.In the multi-objective programming model, the th inequality constraint can be written as g pxq ď 0 and the ' th equality constraint can be written as h 1 pxq ´δ " 0. The infeasibility of a trail solution x can be quantized as follows: Cpxq " ÿ maxpg pxq, 0q `ÿ In Equation ( 23), δ is a permissible deviation, such that δ ą 0 and is very small.If x P X, Cpxq " 0. A dynamic infeasibility threshold ε is used that guarantees the final solutions are all feasible.This threshold is defined as: where ε 0 is the initial allowable deviation of all the constraints.Obviously, ε decreases with the increase in the number of evolutionary generations.In the searching process, the solution x is retained if Cpxq ď ε; otherwise it is discarded.

Selecting the Optimal Particles
The solution of the MOPSO optimization problem is different from a single objective optimization problem.With a single objective problem, it is easy to know which particles are the personal best (pbest) and global best (gbest).With the MOPSO, however, it is difficult to judge which particles are pbest and gbest, because the particles are often non-dominate solutions.However, it is important to pick suitable pbest and gbest particles, since each particle must change its position, as guided by pbest and gbest.Each particle moves toward the non-dominated frontier during the search process [34].
The selection for pbest is relatively simple compared to gbest.A method called Prandom is used in this study, according to which a single pbest is maintained.Pbest is replaced if a new value < pbest, or else, if the new value is found to be mutually non-dominating with pbest, one of the two is randomly selected to be the new best [35].Before the selection for gbest, there are still some works to illustrate.In the MOPSO algorithm, we usually store the non-dominate solutions in archive, and the archive has a limited capacity.Thus, in order to maintain the archive, the crowding distance should be measured as a base for reserving or discarding non-dominate solutions.The crowding distance dt ij can be calculated as: where f l pXq denotes the objective functions in the dimension l.According to Equation ( 25), the crowding distance matrix can be indicated as:  1 2 S S S   .In the same way, archive A can also be divided into three types.Set A1 is the non-dominate solutions, which dominate at least one of the particles in S. Set A2 is the non-dominate solutions which have the same position with the particles in S. Set A3 is comprised of the other non-dominate solutions.2 shows the mapping relations of S and A. For the MOPSO algorithm, the diversity and convergence of population are contradictory issues.One contradiction is the diversity, which guarantees the global best while avoiding the local optimal.The other is the convergence, which promotes particles approaching the Pareto frontier as far as possible.Hence, for particles in S1, if we select non-dominate solutions in A1, which dominate the particles as gbest, the search engines would speed up.However, this can lead to a premature problem.For particles in S2, the global best selection strategy would lead those particles moving to less crowded regions to improve the capability of a global search.
Regardless, each non-dominate solution in the archive has its unique feature.To maintain the diversity of an algorithm, each should have a chance to become a global guide.When paired with these factors, two properties, ri f and pi f , are given for non-dominate solutions in the archive.ri f denotes how often the non-dominate solution is selected as gbest, and pi f denotes how many particles in the current population select the non-dominate solution as gbest.Generally, the size of For the MOPSO algorithm, the diversity and convergence of population are contradictory issues.One contradiction is the diversity, which guarantees the global best while avoiding the local optimal.The other is the convergence, which promotes particles approaching the Pareto frontier as far as possible.Hence, for particles in S1, if we select non-dominate solutions in A1, which dominate the particles as gbest, the search engines would speed up.However, this can lead to a premature problem.For particles in S2, the global best selection strategy would lead those particles moving to less crowded regions to improve the capability of a global search.
Regardless, each non-dominate solution in the archive has its unique feature.To maintain the diversity of an algorithm, each should have a chance to become a global guide.When paired with these factors, two properties, f ri and f pi , are given for non-dominate solutions in the archive.f ri denotes how often the non-dominate solution is selected as gbest, and f pi denotes how many particles in the current population select the non-dominate solution as gbest.Generally, the size of f pi should be restricted.If one gbest is selected by too many particles, the result would be particles converging to a limited region.Based on our experience, we use f pi ď 0.05N, where N is the number of particles [36].
Putting the above pieces together, the global best can be selected as follows: (i) For each particle in S1, we select a non-dominate solution that randomly dominates the particle from A1 as gbest, but f pi ď 0.05N is necessary.If no solution is found, the gbest should be selected from the A1 with greater crowding distance and smaller f ri .(ii) For each particle in S2, a random probability model is employed to select gbest from the A2 with greater crowding distance and smaller f ri .
The pseudo-code of MOPSO algorithm depicting the entire process is given as follows: (1) Initialize positions and velocities of all particles.
(2) Set the current particle position as Pbest.

Computational Experiment
In this section, we evaluate the presented model using a set of numerical data from a real case.The problem is solved by the MOPSO method with Matlab 7.01 on a PC with Intel core i5 and 2.4 GHz.Then, the effects of CEBs and green levels on the decision process are comprehensively analyzed.Finally, some managerial insights are presented.

Case Study
We consider the experiment based on a case study from the petrochemical industry.Specifically, data from the Northeast Chemical Sales Company of the China National Petroleum Corporation (CNPC) (Beijing, China) was studied.CNPC is a large group, and its supply chain network is responsible for transporting petrochemicals from plants, via warehouses, to retailers.This transportation operation involves location, routing and inventory decisions, as well as the creation of considerable carbon emissions.In this paper, a section of the operational data was analyzed.Specifically, this case study involves two plants, five potential warehouses with retailing functions and eight retailers.Each trajectory is relative to a routing cost and the amount of carbon emission.The routings and distances are shown in Figure 3.The parameters of the warehouses and retailers are listed in Tables 1-3.In addition, the market inverse demand coefficient is set as 5500, and the attraction coefficient with the environment level of products is set as 5000.CEB is 1, and the routing cost per distance of M-to-W and W-to-R are all equal to q.The order cost is 500, and the capacity of each vehicle is 1500.In addition, f A j " 29.3N j , the carbon emissions per distance are 0.17, and the carbon emissions per inventory are 0.00276.In addition, ˆ29.3A j j f N  , the carbon emissions per distance are 0.17, and the carbon emissions per inventory are 0.00276.

Numerical Analysis
According to the above data and the approach used to solve the question in Section 5, the trade-off between cost and carbon emissions can be shown as Figure 4.The result provides decision makers with decidedly indifferent choices.In conclusion, all the points on the Pareto line are the solutions, but the managers themselves cannot directly decide.If CEBs are incorporated, the optimal solution is unique (Figure 5).Clearly, revenue first increases and then decreases with increasing carbon emissions.The increasing gradient is greater than the decreasing gradient.This result illustrates that a proper carbon reduction policy can improve corporate revenue, but excessive carbon reduction activities would have a negative impact.Figure 5 shows that the maximum attainable revenue is ¥601,230,000.The optimal configuration can be shown as follows: The location decision is to open Cangzhou, Jinan, and Zhengzhou.The routing decision is divided into two parts.(i) As regards the routing from plants to warehouses, the first decision is that Daqing serves Cangzhou and Jinan, while Fushun serves Zhengzhou.(ii) Considering transportation from warehouses to retailers, the routing schedule of the Cangzhou warehouse is Cangzhou-Tianjin-Beijing-Baoding-Cangzhou-Shijiazhuang-Taiyuan-Cangzhou.The routing schedule for the Jinan warehouse is Jinan-Liaocheng-Linyi-Qingdao-Jinan, and the routing schedule for the Zhengzhou warehouse is Zhengzhou-Xinyang-Yuncheng-Zhengzhou.The order quantities of the three warehouses are 4792, 3156 and 2834 tons, respectively.
We are interested in how CEBs affect companies' decision making.As we know, CEBs mainly affect demand.The effect of consumers' environmental preference on demand for products with different carbon emissions is shown in Figure 6.We vary the CEBs from 1 to 1.8 and obtain a series of demand vs. carbon emissions.Clearly, the curves move from left to right as the coefficient increases from 1 to 1.8, which implies that with the same carbon emission levels, larger CEBs lead to greater demand.This is due to the fact that when consumers pay closer attention to environmental protection, enterprises are more likely to take actions that will improve their environmental protection levels.Then, consumers with greater environmental awareness will buy more products from those enterprises with superior eco-friendly operations.This is a virtuous cycle.However, the marginal cost of implementing environmental improvements increases by degrees.As we know, the ultimate goal of enterprise management is to maximize benefits.Similarly, we adjust the carbon emissions variable to analyze the effect of CEBs on revenue (Figure 7).Clearly, revenue increases as the CEBs move from 1 to 1.8.However, the degree of revenue growth is clearly slower than the increasing CEBs.That is to say, the initial improvement brought about by the CEBs greatly affects the operation of supply chain enterprises, but this effect will weaken because of the high costs associated with further reductions of carbon emissions.
trade-off between cost and carbon emissions can be shown as Figure 4.The result provides decision makers with decidedly indifferent choices.In conclusion, all the points on the Pareto line are the solutions, but the managers themselves cannot directly decide.If CEBs are incorporated, the optimal solution is unique (Figure 5).Clearly, revenue first increases and then decreases with increasing carbon emissions.The increasing gradient is greater than the decreasing gradient.This result illustrates that a proper carbon reduction policy can improve corporate revenue, but excessive carbon reduction activities would have a negative impact.Figure 5 shows that the maximum attainable revenue is ¥601,230,000.The optimal configuration can be shown as follows: The location decision is to open Cangzhou, Jinan, and Zhengzhou.The routing decision is divided into two parts.(i) As regards the routing from plants to warehouses, the first decision is that Daqing serves Cangzhou and Jinan, while Fushun serves Zhengzhou.(ii) Considering transportation from warehouses to retailers, the routing schedule of the Cangzhou warehouse is Cangzhou-Tianjin-Beijing-Baoding-Cangzhou-Shijiazhuang-Taiyuan-Cangzhou.The routing schedule for the Jinan warehouse is Jinan-Liaocheng-Linyi-Qingdao-Jinan, and the routing schedule for the Zhengzhou warehouse is Zhengzhou-Xinyang-Yuncheng-Zhengzhou.The order quantities of the three warehouses are 4792, 3156 and 2834 tons, respectively.trade-off between cost and carbon emissions can be shown as Figure 4.The result provides decision makers with decidedly indifferent choices.In conclusion, all the points on the Pareto line are the solutions, but the managers themselves cannot directly decide.If CEBs are incorporated, the optimal solution is unique (Figure 5).Clearly, revenue first increases and then decreases with increasing carbon emissions.The increasing gradient is greater than the decreasing gradient.This result illustrates that a proper carbon reduction policy can improve corporate revenue, but excessive carbon reduction activities would have a negative impact.Figure 5 shows that the maximum attainable revenue is ¥601,230,000.The optimal configuration can be shown as follows: The location decision is to open Cangzhou, Jinan, and Zhengzhou.The routing decision is divided into two parts.(i) As regards the routing from plants to warehouses, the first decision is that Daqing serves Cangzhou and Jinan, while Fushun serves Zhengzhou.(ii) Considering transportation from warehouses to retailers, the routing schedule of the Cangzhou warehouse is Cangzhou-Tianjin-Beijing-Baoding-Cangzhou-Shijiazhuang-Taiyuan-Cangzhou.The routing schedule for the Jinan warehouse is Jinan-Liaocheng-Linyi-Qingdao-Jinan, and the routing schedule for the Zhengzhou warehouse is Zhengzhou-Xinyang-Yuncheng-Zhengzhou.The order quantities of the three warehouses are 4792, 3156 and 2834 tons, respectively.emissions variable to analyze the effect of CEBs on revenue (Figure 7).Clearly, revenue increases as the CEBs move from 1 to 1.8.However, the degree of revenue growth is clearly slower than the increasing CEBs.That is to say, the initial improvement brought about by the CEBs greatly affects the operation of supply chain enterprises, but this effect will weaken because of the high costs associated with further reductions of carbon emissions.We assume that product pricing is a function of carbon emissions.However, it is not a hard and absolute fact that pricing is the single, key factor.Our analysis (Figure 8) shows that higher prices emissions variable to analyze the effect of CEBs on revenue (Figure 7).Clearly, revenue increases as the CEBs move from 1 to 1.8.However, the degree of revenue growth is clearly slower than the increasing CEBs.That is to say, the initial improvement brought about by the CEBs greatly affects the operation of supply chain enterprises, but this effect will weaken because of the high costs associated with further reductions of carbon emissions.We assume that product pricing is a function of carbon emissions.However, it is not a hard and absolute fact that pricing is the single, key factor.Our analysis (Figure 8) shows that higher prices  We assume that product pricing is a function of carbon emissions.However, it is not a hard and absolute fact that pricing is the single, key factor.Our analysis (Figure 8) shows that higher prices generate greater revenue.What is important is that the higher the price of a product, the bigger the revenue will be obtained with lower carbon emissions.This illustrates that a higher price for green products can stimulate a reduction in carbon emissions, but that higher price can also curb product demand (Figure 9).Clearly, product pricing is increasing with a reduction in carbon emissions.When the price increases, the product demand decreases.In this case, the consumer's willingness to pay is the most important factor.Hence, the enterprise should encourage consumers to improve their CEBs, and pay closer attention to purchasing green products.
revenue will be obtained with lower carbon emissions.This illustrates that a higher price for green products can stimulate a reduction in carbon emissions, but that higher price can also curb product demand (Figure 9).Clearly, product pricing is increasing with a reduction in carbon emissions.When the price increases, the product demand decreases.In this case, the consumer's willingness to pay is the most important factor.Hence, the enterprise should encourage consumers to improve their CEBs, and pay closer attention to purchasing green products.

Managerial Insights
In the current business climate, enterprises and consumers have gradually come to recognize the importance of environmental protection.Both businesses and consumers are more inclined to make an effort to reduce carbon emissions.In particular, consumers with greater environmental awareness are happy to pay a premium price for low carbon emission products.This willingness, which is based on increased environmental awareness, provides an opportunity for logistics enterprises.Our study is consistent with the work of Liu et al. [8], which found that, with consumers' greater environmental awareness, more of them are willing to pay higher prices for low carbon emission products.In turn, the enterprises that produce those products can earn greater revenue.revenue will be obtained with lower carbon emissions.This illustrates that a higher price for green products can stimulate a reduction in carbon emissions, but that higher price can also curb product demand (Figure 9).Clearly, product pricing is increasing with a reduction in carbon emissions.When the price increases, the product demand decreases.In this case, the consumer's willingness to pay is the most important factor.Hence, the enterprise should encourage consumers to improve their CEBs, and pay closer attention to purchasing green products.The relationship between product pricing -carbon emissions, and product pricing-product demand.

Managerial Insights
In the current business climate, enterprises and consumers have gradually come to recognize the importance of environmental protection.Both businesses and consumers are more inclined to make an effort to reduce carbon emissions.In particular, consumers with greater environmental awareness are happy to pay a premium price for low carbon emission products.This willingness, which is based on increased environmental awareness, provides an opportunity for logistics enterprises.Our study is consistent with the work of Liu et al. [8], which found that, with consumers' greater environmental awareness, more of them are willing to pay higher prices for low carbon emission products.In turn, the enterprises that produce those products can earn greater revenue.The relationship between product pricing -carbon emissions, and product pricing-product demand.

Managerial Insights
In the current business climate, enterprises and consumers have gradually come to recognize the importance of environmental protection.Both businesses and consumers are more inclined to make an effort to reduce carbon emissions.In particular, consumers with greater environmental awareness are happy to pay a premium price for low carbon emission products.This willingness, which is based on increased environmental awareness, provides an opportunity for logistics enterprises.Our study is consistent with the work of Liu et al. [8], which found that, with consumers' greater environmental awareness, more of them are willing to pay higher prices for low carbon emission products.In turn, the enterprises that produce those products can earn greater revenue.
In addition, the companies that produce low carbon emissions should also make a concerted marketing effort to shift consumers' traditional purchasing decision criteria and transform those buyers into a group with a preference for low carbon emission products and services.This study indicates that the returns can be substantial if consumers who are currently not interested in purchasing environmentally friendly products make even a little progress.Moreover, the results show that low carbon emission operations cost more than the operations that do not consider carbon emissions.However, when the CEBs are positive, an optimal degree of carbon reduction will maximize revenue.Sadly, the unavoidable fact is that most consumers loathe paying to pay premium prices for low carbon emission products.If enterprises are going to implement sustainable decisions, they must be certain of CEBs.

Conclusions
This paper discussed the effects of consumer environmental behaviors (CEBs) on the design of a sustainable supply chain.CEBs not only affect consumers' willingness to pay premium prices for low carbon emission products, but also the overall demand for low carbon emission products.We introduced a sustainable supply chain network model based on the joint optimization of location, routing and inventory, taking carbon emissions into consideration.The distinguishing feature of our model is its consideration of the CEBs, which affect both carbon emission decisions and product demand.
First, a multi-objective model is constructed, which provides a trade-off between costs and carbon emissions.The MOPSO algorithm is used to solve the model, and then a Pareto optimal set can be obtained.After that, we model the revenue function based on the Pareto solutions.In the computational experiments, we test the model by the data from the Northeast Chemical Sales Company of CNPC.We first obtain the Pareto optimal curve, which provides a portfolio of configurations for decision makers.Then, we can use the same technique to obtain the revenue curves from different carbon emissions.Hence, the unique optimal revenue levels and the relevant decisions can be acquired.Finally, the sensitivity of the case study was analyzed.We are interested in the effects of CEBs on the demand and revenue in a three-level supply chain.The results show that more positive CEBs result in greater demand and higher revenue.We also observe that the pricing of low carbon operations is critical.Therefore, enterprises should make marketing efforts to strengthen consumers' environmental preferences.Companies should support their claims to consumers and ensure the degree of CEBs before implementing their carbon emission reduction policies.
Further research is required to determine more specific factors pertaining to CEBs in a supply chain (e.g., the decision makers' appetite for risk, the expectations of market development and the effects of government intervention via carbon emission reduction policies and legislation), so that the model will be more adaptive to real-life scenarios.

Figure 1 .
Figure 1.The framework of supply chain network.

Figure 1 .
Figure 1.The framework of supply chain network.

2
(i) Location cost.The cost of warehouse location is ř jPJ f j y j , where f j is the single cost of opening warehouse j. (ii) Routing cost occurs in the distribution from M-to-W and from warehouse to retailer (W-tod mi z miv , respectively, where t 1 is the M-to-W routing cost per distance; d kj is the distance from manufacturer k to warehouse j; t 2 is W-to-R routing cost per distance; and d mi is the distance from warehouse j (or retailer k k) to retailer k 1

(
ii) Carbon emissions from routing.The routing emissions from the M-to-W and W-to-R transportations are denoted as ř z miv , respectively, where t1 divided into two types.One set (S1) is comprised of particles that are dominated by at least one of the non-dominate solutions in A. The other set (S2) is comprised of particles that are not dominated by any one solution in A.

Figure 2 .
Figure 2. The mapping relationship of archive A and population S.

Figure 2 .
Figure 2. The mapping relationship of archive A and population S.

Figure 3 .
Figure 3.The network of two plants, five potential warehouses and 14 retailers.

Figure 3 .
Figure 3.The network of two plants, five potential warehouses and 14 retailers.

Figure 4 .
Figure 4. Pareto optimal curve between cost and carbon emissions.

Figure 5 .
Figure 5.The relationship between carbon emissions and revenue.

Figure 4 .
Figure 4. Pareto optimal curve between cost and carbon emissions.

Figure 4 .
Figure 4. Pareto optimal curve between cost and carbon emissions.

Figure 5 .Figure 5 .
Figure 5.The relationship between carbon emissions and revenue.

Figure 6 .
Figure 6.The demand in different consumer environmental preferences varying with carbon emissions.

Figure 7 .
Figure 7.The revenue in different consumer environmental preferences varying with carbon emissions.

Figure 6 .
Figure 6.The demand in different consumer environmental preferences varying with carbon emissions.

Figure 7 .
Figure 7.The revenue in different consumer environmental preferences varying with carbon emissions.

Figure 7 .
Figure 7.The revenue in different consumer environmental preferences varying with carbon emissions.

Figure 8 .
Figure 8.The revenue vs. carbon emissions in different product pricing.

Figure 9 .
Figure 9.The relationship between product pricing -carbon emissions, and product pricing-product demand.

Figure 8 .
Figure 8.The revenue vs. carbon emissions in different product pricing.

Figure 8 .
Figure 8.The revenue vs. carbon emissions in different product pricing.

Figure 9 .
Figure 9.The relationship between product pricing -carbon emissions, and product pricing-product demand.
Figure 9.The relationship between product pricing -carbon emissions, and product pricing-product demand.
where τ( ) g is the CEBs of consumer group g , β is a correction factor of CEBs over time.We assume β 1  , because CEBs would not decrease over time; g is the consumer group with the worst CEBs, and g is the consumer group with the best CEBs.p is the price of the product, which is decided by θ and the CEBs  .The market demand of a product depends on p , θ and  .Conversely, operation- respectively, where ĥj is the carbon emissions per holding inventory.It is worth mentioning that carbon emissions from inventory mainly refer to the energy consumption and product emissions during storage.(iv) Other emissions, including emissions from purchasing, production and recovery.The purchasing emission is Pur 1 .ř iPI µ i , where Pur 1 is carbon emissions from purchase per unit.The production emission is Pn 1 .ř iPI µ i , where Pn 1 is carbon emissions from production per unit.The recovery emission is Rcy 1 .ř iPI µ i , where Rcy 1 is carbon emissions from recovery per unit.iPI ř mPpIYJq t 2 d mi z miv `ř jPJ ř iPI ph o

Table 1 .
The parameters of potential warehouses.

Table 2 .
The area and fixed location cost of potential DC.

Table 3 .
The parameters of retailers.

Table 1 .
The parameters of potential warehouses.

Table 2 .
The area and fixed location cost of potential DC.

Area of Location (m 2 ) Fixed Location Cost (¥) Hold Cost (¥/ton Day)
Figure 6.The demand in different consumer environmental preferences varying with carbon emissions.