Designing a Sustainable Green Closed-Loop Supply Chain under Uncertainty and Various Capacity Levels

: The ever-increasing concerns of the growth in the volume of waste tires and new strict government legislations to reduce the environmental impact of the end-of-life (EOL) tires have increased interest among companies to design a sustainable and efficient closed-loop supply-chain (CLSC) network. In the real world, the CLSC network design is subject to a variety of uncertainties, such as random and fuzzy (epistemic) uncertainties. Designing a reliable and environmentally cautious CLSC with consideration of risks and the uncertainty of the parameters in the network is nec-essary for a successful supply-chain network. This study proposes a sustainable and environmentally cautious closed-loop supply-chain network for the tire industry, by considering several recovery options, including retreading, recycling, and energy recovery. This study aims to design and develop a robust multi-objective, multi-product, multi-echelon, multi-cycle, multi-capacity, green closed-loop supply-chain network under hybrid uncertainty. There are two types of uncertainties associated with the parameters in the network. There is an uncertainty associated with the demand, which is expressed in some future scenarios according to the probability of their occurrences, as well as fuzzy-based uncertainty associated with return rates, retreading rates, recycling rates, procurement, and production costs, which are expressed with possibilistic distributions. In order to deal with this hybrid uncertainty, a robust fuzzy stochastic programming approach has been proposed, and the proposed mixed integer programming model is applied to a case study in the tire industry to validate the model. The result indicates the applicability of the proposed model and its efficiency to control the hybrid uncertainties and the risk level in the network.


Introduction
In today's business world, companies need to pay close attention to their supplychain management systems and logistics network, in order to effectively compete in a highly competitive market. A great supply-chain system is a necessity for any manufacturing organization. Most manufacturers put their focus primarily on the forward supply chain; however, the key to success is to have a holistic view of the entire forward and reverse logistic network. Reverse logistics is a network where selected products transfer up the stream from the end users for repair, remanufacturing, or recycling [1]. Environmental concerns, government regulations, and the potential economic benefits are some of the motives for manufacturers to move toward this direction [2,3]. Besides environmental regulations, consumer's awareness and their responsibility against the environment are some of the motives. The tire industry is one of the industry sectors moving toward this direction. Worldwide, more than a billion tires are manufactured every year, and a significant number of them are thrown away after being used. According to the US Environmental Protection Agency (EPA), in the United States alone, 280 million tires are being scraped every year. The average passenger tire weighs about 22 pounds, and the average commercial semi-truck tire weighs about 110 pounds, and these tires contain oil and other combustible carbon compounds, which create a fire hazard. Figure 1. Tire fires are difficult to extinguish and can break the tire down into hazardous compounds, including gases, heavy metals, and oil. Tire piles can also create a breeding ground for pests such as rodents and mosquitos that contribute to the spread of diseases [4,5]. Although the forward tire supply chain is well developed, there is no systemic and comprehensive reverse logistics system in place. Most collection, retreading, and recycling centers are small individual shops with a limited capacity. Government regulations and requirements vary from state to state. This lack of a comprehensive plan can be seen by a large volume of tires still ending up in the landfills. There are many factors that are considered by the decision maker (DM) when designing a closed-loop supply-chain (CLSC) network, such as costs, environmental impacts, social dimensions, uncertainty in the parameters, governmental regulations, disruption and operational risks, reliability of the network, etc. Uncertainty with the parameters is probably the biggest challenge facing the supply-chain network (SCN) design. Random and fuzzy uncertainties are generally the two types of uncertainty in the design of a supply-chain network [6]. Possibilistic, stochastic, and robust programming are typically the methods utilized to control these uncertainties. Stochastic programming can be used when there are sufficient data available to estimate the probability distribution of the uncertain parameters that is used for modeling the uncertain parameters under different scenarios in the network [7]. For instance, several scenarios for demand may be related to economic conditions, occurrence of natural disasters, changes in customer preferences, development of new technologies, etc. However, the value of specific parameters in each scenario could be imprecise and can be determined by possibilistic distribution.
Possibilistic programming is used when there are no sufficient historical data available, but the possibility distribution of uncertain parameters can be estimated by using the expert's knowledge in the field. Robust programming is a method that allows the DMs to adjust the risk level of output results with respect to the uncertainty of parameters [8]. Combining the stochastic, possibilistic, and robust programming can supplement the weaknesses of each individual method and provide more reliable results for the DMs in the supply-chain network.
In this study, robust fuzzy stochastic programming was used to cope with hybrid uncertainties in the model. Moreover, in this paper, we provide a holistic and realistic modeling approach for the coordination of a forward and reverse tire management system that extends the life span of re-manufacturable tires, recovers the value of scrap tires, and eliminates the disposal in the landfills. In that respect, we examine different recovery options simultaneously, such as retreading, recycling, and energy recovery. Our model aims to maximize the total CLSC profit, minimize the total greenhouse gas (GHG) generated along the CLSC and maximize the social responsibilities which is presented in the form of the total number of jobs created by expanding the network in the forward direction and establishing a new reverse logistics network.
As the body of literature about CLSC network design shows, most models ignore the fact that products can be remanufactured multiple times and remain in the supply-chain network to satisfy the customer demand, rather than being recycled prematurely. We refer to this as a multi-cycle supply-chain network. This can impact model assumptions and outputs with respect to overall market demand, the production rate of manufacturing centers, the return rate of end-of-life (EOL) products, and the capacity of the facilities in the network. Also, a genuinely green supply-chain network must consider and plan for the return of EOL remanufactured products for final processing. This has also been neglected in most previous studies.
This study aims to design a closed-loop supply-chain network considering the economic, environmental, and social factors with the following contributions:


Proposing a new sustainable, robust, multi-objective, multi-product, multi-echelon, multi-cycle, multi-capacity, green closed-loop supply-chain network by considering the real-world situation where used tires are retreaded multiple times, and the retread tires are being returned to the collection center for final processing.  Presenting a novel Robust Fuzzy/Possibilistic Stochastics Programming (RFSP) model that controls scenario-based uncertainty (disruption risks) and fuzzy-based uncertainty (operational risks) in the network simultaneously and provide a more reliable output results for the managers and the DMs.  Evaluating the possibility and feasibility of expanding the existing forward logistic network at different capacity levels to cope with market demand under various scenarios, and also the feasibility of establishing a new reverse logistic network at different capacity levels to cope with the return rate of EOL tires.  Validation of the proposed model by sensitivity analysis through a case study.
The paper is structured as follows. In Section 2, the relevant literature on closed-loop supply-chains network design and tire industry is reviewed; some works concerning uncertainty in the supply chain are also analyzed. In Section 3, the main problem characteristics, assumptions, and mathematical formulation are presented. In Section 4, a proposed Robust Fuzzy/Possibilistic Stochastics Programming (RFSP) approach and augmented epsilon constraint method are described in detail. In Section 5, the proposed model is applied to an example problem inspired by a case in the U.S. The computational results and the sensitivity analysis are also presented in this section. Finally, in Section 6, some conclusions and future work are discussed.

Literature Review
Designing a robust and flexible supply-chain network (SCN) that includes optimizing the number, capacity, and location of network facilities has received a lot of attention these days [9]. A closed-loop supply-chain network (CLSC) considers both forward and reverse supply-chain networks in an integrated manner. There have been many papers published in this subject area and especially green closed-loop supply-chain network. Govindan et al. (2014) and Devika et al. (2015) provided a comprehensive survey of the CLSC network [10,11]. Moreover, a comprehensive review on environmentally conscious manufacturing and product recovery was provided by Ilgin and Gupta (2010) [9]. Jindal and Sangwan (2016) presented the design and optimization of a multi-objective closed-loop supply chain, considering the economic and environmental factors under uncertainty using the interactive ε-constraint method [12]. Sarkar et al. (2017) studied the effect of transportation cost and carbon emission cost on remanufacturing. They proposed a multi-echelon closed-loop supply-chain model with third-party logistics (3PL) [13]. In another study, they developed a multi-attribute closed-loop supply-chain model for self-healing polymers-based returnable transport packaging. Their proposed multi-objective supply-chain model aimed to maximize profit while minimizing carbon emission of the system [14]. Dulman and Gupta (2018) studied the financial impact of sensors on closed-loop supply-chain systems [15]. Zhou and Gupta (2019) provided a pricing strategy for new and remanufactured high-tech products. Their model aimed to maximize the profit associated with customer demand based on the manufacturer, retailer, and joint supply chain [16]. In another study, they proposed partial least square method to explore the factors that affect value depreciation rate and price differentiation between new and remanufactured iPhone and iPad [17]. Aldoukhi and Gupta (2019) proposed a new model for designing a closed-loop supply-chain network by considering downward product substitution policy under four carbon-emission regulation policies. They used robust optimization to deal with uncertainty related to product demand and number of returned products [18]. In another study, they proposed a multi-objective model to design a closed-loop supplychain (CLSC) network, with the aim of minimizing the total cost, minimizing the carbon emission, and maximizing the service level of the retailers. They used physical programing approach to model the problem [19]. Fadhel and Gupta (2019) evaluated the food waste valorization alternatives from a sustainability point of view. They estimated energy utilization and GHG emission reduction for each potential food waste processing technique [20]. Gupta et al. (2019) studied the optimal control of production and maintenance operations in smart custom manufacturing systems, using a heuristic algorithm [21]. De and Giri (2020) studied a closed-loop supply chain (CLSC), focusing on managing, scheduling, and routing problems, to achieve economic and environmental sustainability [22]. Mogale et al. (2018) designed a dual market facility network by considering manufacturing the common components centrally and locally [23]. Oliveira and Machado (2021) provided a literature review on the application of optimization methods in a closed-loop supply chain [24]. Zahedi et al. (2021) proposed a closed-loop supply-chain network by considering multi-task sales agencies and multi-mode transportation, and they solved their model by using metaheuristic and hybrid algorithm [25]. Guarnieri et al. (2020) proposed a theoretical model to assess reverse logistics using multi-criteria decision-making for e-waste [26]. Vieira et al. (2020) studied the barriers in implementing reverse logistics for e-waste and investigated the widely applied multi-criteria decision aid (MCDA) methods used to prioritize these barriers [27].
There have been many researches in the area of sustainability [28]. Mogale et al. (2019) designed a sustainable food-grain supply-chain distribution system with the aim of minimizing the cost and carbon dioxide emission simultaneously [29]. Ullah et al. (2020) investigated the optimal remanufacturing strategy and reusable packaging capacity under stochastic demand and return rate [30]. Ghadge et al. (2020) studied obstacles in sustainability implementation in cold-food supply-chain networks by considering a case of UK artisan cheese producers [31]. Gupta et al. (2020), presented a compilation of six recent papers on a variety of topics to demonstrate the pioneer research activity within responsible and sustainable manufacturing [32]. Sazvar et al. (2021) designed a sustainable CLSC for pharmaceutical by considering manufacturer's brand and waste management [33].
Inconsistent, unreliable, and uncertain data are one of the biggest challenges facing the SCN [34][35][36]. In the real world, there are many parameters in SCN that are uncertain, such as customer demands, costs of opening new facilities at a specific location, manufacturing cost, cost of raw material, transportation cost, etc. [37]. According to Mula et al. (2006Mula et al. ( /2007, the two main types of uncertainty are randomness and epistemic [38,39]. Stochastic programming and possibilistic programming are the most common methods to deal with randomness and epistemic (fuzzy) uncertainty, respectively [40,41]. A flexible robust optimization approach for scenario-based stochastic programming models was proposed by Mulvey et al. (1995) [42]. Their approach was later developed by Yu and Li, 2000 [7].  also developed a new approach, called robust possibilistic programming [43]. Mousazadeh et al. used the robust possibilistic approach for a health-service network design [44]. Tozanli and Gupta et al. (2017) provided a literature review paper for environmentally concerned logistics operation (ECLO) in fuzzy environment. They utilized over 800 papers between 1994 and 2017 and performed a detailed analysis of ECLO with an emphasis on fuzzy application [45]. Vahdat et al. (2017) provided a twostage stochastic programming modeling to design a multi-period, multistage, and singlecommodity CLSC network under uncertainty [46]. Parakash et al. (2018) proposed a CLSC network under risks and demand uncertainty [47]. Yolmeh and Saif (2020) provided a supply-chain network design with assembly and disassembly line balancing under uncertainty [48]. Wang et al. (2021) studied various incentive mechanism in green supply chain under demand uncertainty [49]. Tordecilla et al. (2021) provided a review of simulation optimization methods for designing a resilient supply chain under uncertainty [50].
There has been some research on the tire CLSC network. In 1997, Ferrer published a paper and estimated the optimal number of times each tire can be retreaded [51]. Sasikumar et al. (2010) provided an optimization model for a tire remanufacturing case [52]. According to Chopra and Meindl (2015), between 2009 and 2011, only 3 percent of total sales by US firms within the tire industry were retread tires [53]. One of the challenges facing the tire industry is the coordination of the entire forward and reverse logistics network and managing the flow of the new and used tires [54].  [64]. Mehrjerdi at al. (2020) studied sustainability and resiliency of a closed-loop supply chain, using multiple sourcing strategy [65]. Lokesh et al. (2020) also designed a sustainable tire supply chain with a fuzzy goal planning approach [66].

Tire Industry and Recovery Options
According to Modern Tire Dealer (MTD), in 2019, 267.8 million new passenger tires, 37.9 new light truck tires, and 25.4 million new medium and heavy truck tires were shipped within the United States. EOL tires due to the large quantities being manufactured and the durability of the material being used, pose a significant challenge for the government and the supply chain. Disposing scrap tires in landfills consumes valuable space and can also result in the penetration of pollutants and metal contents into the underground water. With recent research and advancement in the composition of tire, tires have become even more durable. When the tire treads are worn out, the majority of the time, the tire casing is still in good condition and can be reused many times. Disposing of in the landfill or incineration for energy recovery only shortens the life span of the tire casing. There are few recovery options available for EOL tires, some of which are as follows: Reuse or repurposing the tire is the most environmentally friendly option that extends the life cycle of the tires. Tires with sufficient thread depth can be sold in the used tire market. Whole tires can also be used for a number of applications, including erosion control, playground equipment, crash barriers on the side of the highway, and boat bumpers at the docks.
Retread or recap is a process in which the worn thread of the tire is removed from the tire casing and replaced with new treads. Study shows 25 gallons of oil needed to produce a new truck tire compared to only 7 gallons of oil to retread a used one. That results in a significant reduction in greenhouse gas emission, preserving raw material, and minimizing the amount of waste disposed of in the landfills.
Waste tires are typically shredded before they are recycled. Recycled tires can be used as crumb rubber in asphalt pavement or as an aggregate in Portland cement concrete. Tires are also cut up and used in garden beds to hold in the water and to prevent weeds from growing [67].
A large percentage of scrap tires are used as fuel. Tire-derived fuel (TDF) is a fuel derived from scrap tires. Scrap tires, either shredded or as a whole, is used to supplement coal in power plants, cement kilns, paper mills, etc. Research has shown that atmospheric contamination dramatically increases when tire rubber is used as fuel [68][69][70][71][72].
Disposal in the landfill is the most harmful option due to the large volume and high percentage of void space inside the tire. Tires can trap gas that causes them to become buoyant and damage the landfill covers.
Designing an economically optimized and environmentally friendly CLSC network is a prerequisite for tire manufacturers not only to earn a profit but also to decrease waste, preserve natural resources, and save landfills and the environment, with the ultimate goal of sustainable development. This has motivated this study to take a closer look at the challenges facing the tire industry and come up with a solution that creates a win-win situation for companies in the tire industry network and the environment. Based on the reviewed papers reported in Table 1, in this study, we aimed to identify the uncertain parameters in the logistics network and find the solution for dealing with these uncertainties in order to reach a robust modeling in the logistics network. We introduce a robust multiobjective, multi-product, multi-echelon, multi-cycle, multi-capacity, green closed-loop supply-chain network under hybrid uncertainty with the aim of reducing costs, reducing environmental impacts, and increasing the social responsibilities and the reliability of the network.

Problem Description and Formulation
As shown in Figure 2, we are considering a multi-product and multi-echelon CLSC network consisting of suppliers, production centers, distribution centers, collection centers, retreading centers, recycling, and energy-recovery centers. In the forward flow, various raw materials are procured from the available choice of suppliers. After tires are manufactured at the production facilities, they are distributed through the distribution centers and the retailers to the end users. The location of retailers is fixed, and there might be unsatisfied demands due to uncertainty in demand. In the reverse flow, used tires are transferred to collection centers, where the initial inspection is performed, and based on the condition of the tires, they are either routed to retreading center or the recycling centers.
The main objective in the CLSC network design problem is to maximize the total profit of the supply chain by choosing the optimal number and location of facilities, their capacities, and the flow of products between the facilities, while reducing the environmental impacts and maximizing the social benefits. In real life, the product demand could be a function of economic condition; therefore, in this study, the future demand was expressed under a few scenarios. The return rates, retreading rates, recycling rates, procurement, and production costs are fuzzy; therefore, the possibilistic approach was utilized to deal with these parameters.

Problem Assumption and Notation
The main characteristics and assumptions are as follows:  The potential location of new manufacturing centers, distribution centers, retreading centers, and recycling centers is fixed and predefined.  Production, distribution, retreading, and recycling capacities are limited, and all are subject to capacity level U.  Capacity of the energy recovery centers is unlimited, but it has the most negative environmental impact.  The return, recycling, and retreading rates of new and retread tires are fuzzy.  Raw material procurement costs and production costs are fuzzy, and their fuzzy expected value are utilized.  Customer demand is stochastic and is expressed in different scenarios.  Ground transportation is the only transportation method used in the network.  There is cost associated with the collection of used tires. This provides an incentive for the end users to return the EOL tires.  Tires can be retreaded two times, retread-1 and retread-2, and each tire will be identified how many times it has been retreaded.  The value of the tires that were retreaded once (retread-1) is greater than the value of the tires that were retreaded two times (retread-2), and both are less than the value of the new tires.  The demand for new tires must be satisfied, and there is a penalty for unsatisfied demand. There is no obligation to meet the demand for retread tires. Table 2 shows the list of indices, parameters and variables related to the network under study.  Return rate of EOL retread-1 tire from center c (fuzzy) 2 Return rate of EOL retread-2 tire from center c (fuzzy) 1 Retread rate of EOL new tires (fuzzy) 2 Retread rate of EOL retread-1 tires (fuzzy) Recycling rate (fuzzy) Unit price of the new tire p Unit price of the retread-1 tire p Unit price of the retread-2 tire p Value of the recycled material for unit tire p

Problem Formulation
In the following, we introduce three objective functions , , and , under scenario , representing the total profit of the network, the total environmental impacts, and the social benefits, respectively.
Equation (1) maximizes the total profit of the supply chain. The components of the Equation (1) include the following: the overall cost of opening new centers, cost of procurement of raw material, overall manufacturing cost, collection cost, overall retreading cost, overall recycling cost, overall energy recovery cost, total transportation cost, revenue from sales of new tires, revenue from sales of tires that are retreaded once, revenue from sales of tires that are retreaded twice, and revenue from sales of recycled material.
Equation (2) minimizes the overall environmental impact, which is expressed as the amount of greenhouse gases (GHGs) generated at production centers, retreading centers, recycling centers, and energy recovery centers, as well as during the transportation of raw material and finished tires between centers. Among greenhouse gases, this research is focusing on carbon dioxide (CO2) emission.
Constraint (4) ensures that the total quantity of raw material sent to tire manufacturing centers is less than or equal to the maximum capacity of each supplier for that raw material. Constraints (5) and (6) ensure that the total quantity of tires produced at each manufacturing center is less than or equal to the maximum capacity of each existing and potential manufacturing center. Constraint (7) ensures that only one production facility at capacity level u can be open at location f. Constraint (8) guarantees that the raw material needed at production centers are satisfied. Constraint (9) is an equilibrium constraint indicating the total quantity of tires transferred from each manufacturing center f to distribution center d is equal to the quantity of tires manufactured at that center. Constraints (10) and (11) ensure that the maximum capacity of the existing and potential distribution centers are not violated. Constraint (12) guarantees that only one distribution center at capacity level u can be open at any potential location d. Constraints (13)-(15) are equilibrium constraints to ensure the total quantity of new and retread (retread-1 and retread-2) tires transferred out of each distribution center is equal to the total quantity of tires enter that center. Constraint (16) ensures that the market demand for new tires is satisfied. Constraints (17) and (18) indicate that the demand for retread-1 and retread-2 tires does not have to be entirely satisfied. It should be noted that the demand parameters are stochasticscenario-based. Constraints (19)-(21) calculate the total quantity of tires returned to each collection center. It should be noted , 1 , 2 are fuzzy parameters that will be discussed in the next section. Constraint (22) controls the maximum capacity of the collection centers. Constraints (23)- (25) are equilibrium constraints which create flow balance at each collection center. Constraint (26) ensures the flow of products into and out of each retreading center do not exceed the capacity of that center. Constraints (27) and (28) calculate the quantity of the tires that are being successfully retreaded at each retreading center. Moreover, ( 1 ) and ( 2 ) are fuzzy parameters. Constraints (29) and (30) calculate the total quantity of new and retread tires that cannot be retreaded at retreading centers and are being transferred to recycling centers. Constraint (31) controls the capacity of the recycling centers. Constraint (32) calculates the total quantity of products transfers to each recycling center. It should be noted that is a fuzzy parameter. Constraint (33) creates equilibrium at recycling centers. Constraint (34) calculates the total number of tires being incinerated at each energy recovery center. Ultimately, constraint (35) enforces the binary and non-negative constraints on the corresponding decision variables.
The proposed tire CLSC model is a Mixed Integer Linear Programming (MILP) problem under hybrid uncertainty. Fuzzy uncertainty related to return rate, retreading rate, recycling rate, procurement cost, production cost, and scenario-based uncertainty related to demand parameters. In the following, we discuss a method that can simultaneously deal with these hybrid uncertainties.

The Proposed Approach
Robust Fuzzy/Possibilistic Stochastics Programming (RFSP) is a novel approach in mathematical programming that deals with hybrid uncertain parameters [6]. It is an expansion of possibilistic programing [43] in which the possible scenarios of the disturbance in the system are additionally being considered and addressed using stochastic scenariobased programming method [42]. An example of a disturbance in the system would be a change in customer behavior that creates a noticeable change in demand in comparison to the initial forecast. These changes can create a new scenario in the system with respect to demand parameters.
Since the proposed tire CLSC problem faces both scenario-based uncertainty related to demand parameters and fuzzy uncertainty associated with return rates, retreading rates, recycling rates, procurement, and production costs, Robust Fuzzy/Possibilistic Stochastics Programming method is proposed in this research to deal with hybrid uncertain parameters.
In the following, we first describe the basics of robust stochastic programming and possibilistic programming. Then, the proposed RFSP approach is explained; and at the end, this approach is used to control the uncertainties in the proposed tire CLSC problem.

Stochastic Programming
In general, the scenario-based stochastic programming approach can be described as follows: where is the objective function under scenario ∈ , are the scenario dependent variables, are the variables independent of the scenario (i.e., binary variables), , , , , are the parameters under scenario ∈ , ( , ) are crisp parameters and is the probability of occurrence of scenario ∈ [7,73].
Since scenario-based programming models are not sensitive to deviations of the value of the objective function, in order to prevent excessive deviation of the of the objective function under different scenarios, the variance of the objective function is also considered. Therefore, the resulting robust stochastic programming model would be as follows: where denotes the importance or weighing scale to measure the tradeoff between optimality and the cost under different scenarios. denotes the absolute deviation from the mean under each scenario.
represent the deviation for violation of the control constraint, and ω is the robustness coefficient or feasibility robustness multiplier which its value can be equal to the penalty cost of the shortage. It should be noted that the above model is known as the Mulvey Robust stochastic model [42], although the above model is the linear version of the model based on the research of Yu and Li [7].
In the above model, the following linearization was utilized:

Possibilistic Programming
If some of the data in an optimization problem are imprecise due to insufficient data or unavailability of the information, but their approximate values can be expressed in fuzzy numbers, fuzzy mathematical programming can be used [74]. If fuzzy data in the optimization problem are expressed as Trapezoidal Fuzzy Number (TFN) and the possibility of occurrence of the possibilistic parameters are measured by Pos. and Nec., then possibilistic programming can be used, where Pos. is the maximum level of occurrence of parameters and Nec. is the minimum possibility level.
The general form of possibilistic programming can be expressed as follows: Where Me(.) is a fuzzy measure; ℱ = ( ( ) , ( ) , ( ) , ( ) ) and ̃= ( ( ) , ( ) , ( ) , ( ) ) are the fuzzy parameters in the form of TFN, Figure 3; is the minimum confidence level; and ( ̃ ) is the fuzzy expected value of the objective function, which can be expressed as follows: And the following relationship can exist [70]: One of the disadvantages of classical programming is that the value is determined subjectively by the decision-maker (DM) and is usually set inactively or ultimately interactively and may not be optimal. This issue was first addressed in the research conducted by Pishvaee (2012) [43].
In their proposed model, taking into account the two concepts of feasibility robustness and optimality robustness, the optimal levels of value are obtained proactively and objectively. This approach has been used in several studies [74][75][76][77][78][79][80].

Robust Fuzzy/Possibilistic Stochastic Programming
The RFSP approach is an extension of the combination of scenario-based stochastic programming and possibilistic programming [43] in which random scenarios and fuzzy data are controlled simultaneously [6]. The RFSP approach can be summarized as follows: where is the coefficient for estimating fuzzy numbers, ( ) is a non-negative descending function that calculates the violation of fuzzy constraints with confidence level . Moreover, ( ) can be expressed as the difference between the substituted crisp numbers, ( ( ) − (2 − 1)( ( ) − ( ) )), and the worst-case value of the imprecise parameters (TFNs) can be expressed as follows: calculates shortage under each scenario. It should be noted that denotes the weighing scale to measure the tradeoff between optimality and the cost under different scenarios, denotes the absolute deviation from the mean under each scenario. represents the deviation for violation of the control constraint, and is the model robustness coefficient or feasibility robustness multiplier.

Coping with Multi-Objective Functions
In order to solve multi-objective decision-making (MODM) optimization problems, several approaches have been proposed, such as Weighted Sum Method (WSM), Epsilon Constraint (EC), Augmented Epsilon Constraint (AEC), Goal Programing (GP), Lexicographic (Lex), and so on. The general form of a MODM problem is as follows:

Epsilon Constraint
Let us assume the first objective function is the main objective, and the remaining objectives will have an upper bound of epsilon and will be presented as constraints of the model. In this case, the EC method is being used, and the resulting single objective model can be described as follows: In the above model, the first objective is considered as the main objective, and the remaining objectives are considered as constraints and are bound by maximum . By changing the , different results can be obtained, which may be weakly efficient or not efficient.
This issue can be resolved with minor modification, which is known as the AEC method [81]. To better implement the AEC method, the appropriate range of epsilons can be obtained by using the Lex method [82]. In the AEC method, the appropriate range of is first determined and then the Pareto front obtained for different values of .

Appropriate Interval for with Lex Method
In order to find the appropriate interval for related to the objective function (i = 2, …,n), the following optimization problem is first being solved: where = ( , * ) is the optimum value of and , * is the vector of decision variables that optimizes .
Then, with the solution that optimizes the objective function , the value of the remaining objective functions is obtained as follows: PayOff = Min f (x) f (x) = PayOff x ∈ X j = 1,2, … , n ; j ≠ i (48) where = ( , , * ) is the optimum value of and , , * is the vector of decision variables that optimizes . Using the Lex method, the following payoff matrix is obtained: After calculating the payoff table for i=1,…..n, the following can be defined: Using the above definition, the appropriate rage for each e can be obtained by using the Lex method: The value R(f ) is used to normalize the objectives in the AEC objective functions.

Augmented Epsilon Constraint
The AEC model can be expressed as follows: where s are non-negative variables, ϕ is a parameter for normalizing the first objective with respect to f and ϕ = ( ) ( ) . In the proposed AEC method in this study, first e ∈ [Min(f ), Max(f )] is obtained using Lex method, then Equation (45) is solved, which provides an optimum solution.

Computational Experiment
To check the validity, practicality, and performance of the model, this model was applied for the design and optimization of a network in Northeast USA, as illustrated in Figure 4. The relevant data presented in Tables 3-6  In this research, there are ten types of medium/heavy truck tires considered. There are five groups of raw materials (synthetic polymers, natural rubber, steel, fillers, and curing compounds) required to manufacture a tire and five suppliers to purchase the raw material from, which are all located out of the state. There are also two existing production facilities and five distribution centers at medium-capacity level, two energy recovery centers, and forty commercial retail locations. We are considering the feasibility of expanding the network by adding two additional manufacturing facilities, five additional distribution centers, ten collection centers, four retreading centers, and five recycling centers. We believe that no tire should end up in the landfills, and by providing appropriate incentives, we can encourage customers to return the EOL tires for collection and processing. These facilities are selected such that they are in the proximity of the facilities they are linked to in the network flow and also to provide maximum coverage to the entire region. The goal is to optimize the number, location, and capacity of the facilities, increase the efficiency and profitability of the network and reduce the amount of greenhouse gases generated in the network while maximizing the total number of jobs that can be created. The demand for truck tires depends on the performance of the logistics and transportation industry in the country, which is directly correlated to the condition of the economy. We assume three scenarios for future demand, corresponding to demand at the time of poor economy (recession), during a good steady economy, and when the economy is growing. Ground transportation is the only mean of transferring the tires and raw material between the facilities. Assuming a 22-foot full truckload on each transfer, the cost of transferring one unit of tire per kilometer is $0.02 and $0.002 for one unit of raw material. The amount of Co2 generated per unit of tire transferred between facilities, using ground transportation is 0.002 kilogram per kilometer and is 0.0002 kg per kilometer for one unit of raw material being transferred.

Computational Results and Sensitivity Analysis
In order to evaluate the performance of the proposed model, first, the bi-objective version of the model is solved by considering the maximization of the profit and minimization of the environmental impact for different values of λ. The model is solved using GAMS software and a PC with the following configuration: Intel ® Core™ i7-1165G7 Processor (2.80 GHz) and 16 GB DDR4 3200MHz. Figures 5 and 6 illustrate the changes in objective functions with respect to changes to the risk factor (λ). It can be observed that by increasing λ, meaning, by penalizing the deviation of the objective function from the mean value (reducing the risk taken by DM), both total profit and scenario variability are decreased. The optimality robustness of the model can be controlled by adjusting the risk factor (λ) to DM's preferences. Reducing the risk factor to the lowest (λ = 0) results in the highest profit and highest variability.  Table 7 shows the optimum configuration of the network, using the mean approach (λ = 0). In this configuration, two additional production centers will be open at high capacity and five additional distribution centers will be open, namely two at low, two at medium, and one at high capacity, to cope with the demand. In the reverse logistics, six collection centers will be open at low capacity; three retreading centers will be open, one at medium and two at high capacity; and two recycling centers open at low capacity. The network flow for this configuration is illustrated in Figure 7. The penalty cost for unsatisfied demand or shortage (ω) is another tool at DM's disposal to manage shortage in the system. An increase in the value of ω results in a decrease in mean profit and an increase in feasibility robustness. It can be seen in the optimum configuration of the network that the appropriate capacity levels of the open centers were determined in order to eliminate the shortage. The optimum confidence level of the fuzzy constraints ( ) was also determined based on the objective functions and the constraints which resulted in zero violation of fuzzy constraints (ϑ( )).   164.926 Δ (shortage) 0 ( ) 0 (confidence level of fuzzy constraints) = 1/ =0.5/ = 0.5/ = 1/ = 1/ = 1 Further sensitivity analysis is performed to evaluate the impact of changes to the problem parameters on the optimal solution and the model behavior. In that respect, we evaluate the impact of increasing the demand by 20% without making any changes to the number of potential facilities. As illustrated in Table 8, an increase in demand has resulted in an increase in total profit and GHG production which is the result of an increase in production. The capacity of one of the distribution centers has also increased from low to high level. It can be noted that there is some unsatisfied demand due to the capacity constraint at production centers.
We further evaluate the model by increasing the demand by 20% and also considering opening an additional production facility in order to eliminate the bottle neck at production centers. As it is shown in Table 9, the total profit further increased, the new production facility is open at a high-capacity level, and the capacity of the other two production facilities decreased to medium, in order to adjust with the overall demand, which resulted in zero unsatisfied demand. There was also an increase in the capacity of distribution centers, collection centers, recycling centers, and retreading centers. This is reflected in the total number of the job opportunities created along the network. 107,151 ( ) 0 (confidence level of fuzzy constraints) = 1/ =1/ = 1/ = 1/ = 1/ = 1 Table 9. Summary of optimum result (demand increase by 20% and considering opening an additional production center). 190.680 Δ (shortage) 0 ( ) 0 (confidence level of fuzzy constraints) = 1/ =0.5/ = 0.5/ = 1/ = 1/ = 1

Value of Coefficients
In this section, the model is evaluated by optimizing all three objective functions, which are (1) maximizing the total profit, (2) minimizing the environmental impacts, and (3) maximizing the total number of jobs created along the network. As illustrated in Table  10, the total number of job opportunities has been increased by opening more centers at a high-capacity level. This has resulted in a reduction in the total profit due to the opening cost of the new facilities. The total amount of GHG remains unchanged, as it is directly correlated with the demand and the quantity of products being processed in each center. The overall unsatisfied demand remains at zero. The optimum confidence level of fuzzy constraints is obtained, which results in zero violation of fuzzy constraints ( ( )). 168.635 Δ (shortage) 0 ( ) 0 (confidence level of fuzzy constraints) = 1/ =1/ = 1/ = 0.5/ = 0.5/ = 0.5

The Managerial Implication
There are many decision factors facing the DMs in designing a complex closed-loop supply-chain network. Consideration of disruption risk in addition to operational risk is one of the biggest challenges facing the DMs and the managers in supply-chain network. An example of disruption risk could be changes to economic condition, introduction of new technologies, changes in customer habit and preferences. Operational risks can be related to uncertainty in demand, procurement, production and transportation cost. RFSP is a powerful tool that provides flexibility to the DM to control these risks in the design of CLSC network simultaneously. By adjusting the risk factors (λ) to the DM's preferences and their degree of conservatism, the optimality robustness of the network can be controlled. Feasibility robustness and the acceptable degree of violation of the constraints of the network can also be controlled by adjusting the value of Ψ and ω.
It is evident that this robust optimization method is capable of controlling scenario variability, possibilistic variability, and unsatisfied demand by adjusting the weighting scale based on the DM's preferences.

Conclusions
Recovery of EOL tires has been a challenge due to the large volume being produced and the material composition of the tires that makes them difficult to decompose. With all recent advancements in recycling options for EOL tires, a significant volume of them still ends up in the landfill. This indicates the importance of a need for a systemic, sustainable CLSC network for the tire industry. In this study, we developed a mixed-integer linear programming model for a CLSC network for the tire industry by considering several recovery options and taking into account the hybrid uncertainty associated with some parameters in the network, with the aim of optimizing the total profit, while minimizing the environmental impact and maximizing the social responsibility of the network. A robust fuzzy stochastic programming method was proposed to cope with the scenario-based uncertainty associated with the demand and fuzzy-based uncertainty associated with return rate, retreading rate, recycling rate, procurement cost, and production cost simultaneously.
To illustrate the performance and practicality of the proposed model, this model was applied for the design and development of a tire CLSC network in the Northeast USA. The result indicates the efficiency and applicability of the model in real-world problems and the feasibility of designing a profitable and environmentally friendly CLSC network for the tire industry. This model and proposed solution can be further extended and applied to the supply-chain network design of other industries.
One of the limitations of exact solution methods is their inability to solve large-scale problems; therefore, for a future study, a larger-scale problem using metaheuristics methods, with a focus on the impact of low-quality, low-cost imported tires and governmental intervention, regulations, and incentives on the performance and profitability of the network, can be considered.