Probabilistic Linguistic VIKOR Method to Evaluate Green Supply Chain Initiatives

To achieve the goals of sustainable economic and environmental protection, more and more firms intend to implement green supply chain (GSC) initiatives in their products. The adoption of GSC initiatives in turn influences the firms’ operations performance. Thus, the firms have to evaluate their performances carefully when implementing GSC initiatives. The performance evaluation of GSC initiatives is a laborious task, which needs to take into account many factors including the inventory level and assurance of supply purchasing-wise and the technical capability and the innovation capability manufacturing-wise, etc. This paper develops a new probabilistic linguistic VIKOR approach to support such an assessment. To do so, a new comparison method of probabilistic linguistic term sets (PLTSs) is first presented to effectively determine the probabilistic linguistic positive ideal solution and the probabilistic linguistic negative ideal solution. Next, a new defuzzification function of PLTSs is proposed to take into account the main-criteria weights and the sub-criteria weights which are represented by PLTSs. Furthermore, several probabilistic linguistic measures are introduced, such as the probabilistic linguistic group utility measure, the probabilistic linguistic individual regret measure and the probabilistic linguistic compromise measure. Finally, the compromise solution is obtained based on these three measures. The desirable advantages of the developed method are summarized as (1) it allows the evaluators to employ PLTSs to express the imprecise performances of the GSC initiatives with respect to various criteria, which greatly improves the elicitation of linguistic information; (2) it presents a probabilistic linguistic compromise solution, which is a maximum probabilistic linguistic group utility for the majority and a minimum probabilistic linguistic individual regret for the opponent. This proposed technique provides a simple and efficient decision making approach to assist the firms to make an appropriate decision in GSC management.


Introduction
Due to increasing public awareness of the need to protect the environment, and the regulatory pressures coming from governments, firms are progressively promoting business practices to assist them in reducing the negative impacts on the environments [1].The implementation of green supply chain (GSC) initiatives is a good choice for firms to obtain profit and market share objectives by lowering environmental impacts and improving ecological efficiency.For example, in the new product development, many GSC initiatives including the green design, the green packaging, the use of green raw materials, etc., can be implemented by various firms.The adoption of new GSC initiatives in turn requires changes in firms' operational capabilities and resources that may have an adverse effect on firms' operations performance [2].Thus, it is necessary for the firms to evaluate their performances showed that the use of the green raw materials and greening the production process can lead to environmental material substitution, waste reduction, and decreasing the consumption of hazardous and toxic materials.In order to enhance cost savings, Zhu et al. [15] proved that green logistics and packaging are good initiatives to cut energy consumption and packaging waste in GSC management.By investigating the relationship between the implementation of GSC management and the economic performance and competitiveness of a sample of Asian firms, Rao and Holt [13] confirmed that GSC management can effectively improve competitiveness.Recently, Govindan et al. [16] investigated the selection problem of green suppliers in the food industry.They showed that the food industry is mainly focused on the quality criterion and the environmental criterion, and the suppliers that are more advantageous in these two aspects are favored.

Application of MCDM Techniques in Green Supply Chain
Usually, the evaluation of GSC management or the selection of green suppliers is involved with multiple conflicting criteria, which is a kind of typical MCDM problem.Over the last few decades, many MCDM techniques have been applied effectively to the evaluation of GSC management and the selection of green suppliers in GSC management.For example, Azadi [17] introduced two-stage targetsetting DEA (data envelopment analysis) methods to the GSC management of public transportation service providers.Mathiyazhagan et al. [18] employed an AHP (Analytical Hierarchy Process) technique to investigate the pressures for GSC adoption in the mining and mineral industry context.Govindan et al. [16] applied PROMETHEE (a preference ranking organisation method for enrichment evaluations) to select the green suppliers in the food supply chain.On the other hand, owing to the fact that the crisp values are not adequate or sufficient for the decision makers (DMs) to express their preferences, fuzzy set and/or its extensions are used in GSC management.Many fuzzy MCDM methods [19][20][21][22] have been applied in GSC management.For instance, Govindan et al. [23] presented an intuitionistic fuzzy DEMATEL (decision-making trial and evaluation laboratory) technique for developing green practices and performances in process of GSC management.Kusi-Sarpong et al. [24] introduced a joint rough sets and fuzzy TOPSIS (technique for order preference by similarity to ideal solutions) to evaluate the GSC practices in the mining industry.Rostamzadeh et al. [25] also proposed a fuzzy VIKOR technique for assessing GSC practices.Based on linguistic preferences contexts, Shen et al. [26] introduced a fuzzy multiple criteria approach for evaluating green supplier's performance in GSC.Tseng et al. [27] introduced a convergence of interval-valued triangular fuzzy numbers-grey relation analysis technique to enhance performance of the GSC management in a Taiwanese electronic focal firm.Moreover, Zhang et al. [28] developed a hesitant trapezoidal fuzzy QUALIFLEX (qualitative flexible multiple criteria method) method to evaluate GSC initiatives in order to achieve sustainable economic and environmental performance under the contexts of hesitant trapezoidal fuzzy numbers-based comparative linguistic expressions.Wang and Chan [2] developed a hierarchical fuzzy TOPSIS approach to assess improvement areas when implementing GSC initiatives.In addition, a series of hybrid MCDM methods [6,[29][30][31][32][33][34] have recently been developed and applied in green supply chain.On the other hand, Mardani et al. [35] presented a systematic review of MCDM techniques in sustainable and renewable energy systems problems.Mardani et al. [36] provided a review on the application and use of MCDM approaches in regard to energy management problems from 1995 to 2015.

The Classical VIKOR Approach
The classical VIKOR approach is one of the well-known MCDM methods, which focuses on the multi-criteria ranking index based on the particular measure of "closeness" to the ideal solution.Consider a MCDM problem which includes a set of feasible alternatives {a 1 , a 2 , • • • , a m } and a set of the predefined evaluation criteria {c 1 , c 2 , • • • , c n }. x ij represents the criteria values of the alternative a i over the criterion c j and w j is used to denote the weight of the criterion c j .In the classical VIKOR method, all criteria values are the crisp values.The decision steps of the classical VIKOR method for dealing with this MCDM problem are summarized as below: Step 1. Determine the positive ideal solution (x j ) + (x j ) + = max m i=1 x ij (1) and the negative ideal solution (x j ) − Step 2. Calculate the value S i and the value R i : Step 3. Determine the value Q i : Step 4. Identify the best solution or the set of compromise solutions Using Equations ( 3)-( 5), three ranking lists of alternatives are obtained according to the decreasing order of the values of S i , R i and Q i , respectively.The alternative a σ(1) (it is the alternative with the first position in the ranking list derived by Q i ) is the best solution if the following two conditions satisfied: Acceptable stability in decision making: The alternative a σ(1) is also the alternative with the first position in the ranking lists derived by S i and/or R i .
If the condition C 1 is not satisfied, then the set of alternatives {a σ(1) , a σ(2) , • • • , a σ( f ) } is the compromise solution in which the maximum value of f is determined by the following equation: . On the other hand, if the condition C 2 is not satisfied, then the set of alternatives {a σ(1) , a σ(2) } is the set of compromise solutions.
The VIKOR approach is widely used to solve complicated decision making problems with crisp values, such as the selection of a renewable energy project [37], post-earthquake sustainable reconstruction [38], financial performance evaluation [39], and selection of the alternative-fuel modes [40], etc.With the rapid development of human social-economic activities, however, it is more and more difficult for the DMs to employ crisp values to express their preferences [41].Instead, various fuzzy forms of decision data have been developed to model the imprecise preferences [42][43][44][45][46].To effectively solve the MCDM problems with various fuzzy environments, many extension forms of the VIKOR approach have recently been developed, such as the interval-valued VIKOR [47], the triangular fuzzy VIKOR [48][49][50], the trapezoidal fuzzy VIKOR [51], the intuitionistic fuzzy VIKOR [52], the interval-valued intuitionistic fuzzy VIKOR [53], the hesitant fuzzy VIKOR method [54,55] and hesitant fuzzy linguistic VIKOR method [56], etc.Recently, Mardani et al. [57] presented a systematic review of the VIKOR technique in the application areas of sustainability and renewable energy.Although useful, these aforementioned VIKOR-based methods fail to solve the real-world MCDM problems with PLTSs.

Basic Concept of PLTSs
The mathematical form of PLTS is introduced as follows.
Definition 1 [4].Given an LTS L = { l i |i = 0, 1, • • • , τ }, a PLTS based on L can be defined as: where L (k) (p (k) ) is the linguistic term L (k) associated with the probability p (k) , and #L(p) is the number of all different linguistic terms in L(p).
To facilitate its understanding and application, Zhang [58] provided a new mathematical form for expressing the PLTS.
τ} be a set of the subscripts of linguistic terms in L, and let Λ ⊆ Γ be a subset of Γ; then the new mathematical expression form of PLTS is provided as below: It is easily observed that the commonly used expression forms of linguistic information, i.e., the ordinal scales-based linguistic variables [59] and the hesitant fuzzy linguistic term sets (HFLTSs) [60], are the special cases of PLTSs because (1) if #Λ = 1 and ∑ σ∈Λ p σ = 1 the PLTS is mathematically reduced to an ordinal scales-based linguistic variable [59] then the PLTS is mathematically reduced to the HFLTS [60] which is considered to have a possibility distribution across the entire LTSs [61].
In practice, the PLTS 2), l 2 (0.6), l 3 (0.1)}, L 3 (p) = {l 4 (0.9)}.It is easy to see that the PLTS L 1 (p) is a complete PLTS, L 2 (p) and L 3 (p) are the partial PLTSs.On the other hand, it is easily observed that #Λ 1 = 2, #Λ 2 = 3, #Λ 3 = 1.That is to say, the numbers of linguistic terms in these three PLTSs are different.In practice, these situations in which the numbers of linguistic terms in different PLTSs assessments are different and some PLTSs assessments are not complete PLTSs but are partial PLTSs usually happens in the decision making process.To ensure the rationality of the operational results among PLTSs, Zhang [58] further introduced the corresponding normalization method for normalizing the PLTSs before conducting the operation of PLTSs.
then the normalization process of the set of PLTSs ℘ can be divided into the following two steps [58]: (1) If the PLTS L i (p) is a partial PLTS, i.e., ∑ σ∈Λ i p i σ < 1, then the PLTS L i (p) should be normalized into the complete PLTS Li (p) as below: , and the probabilities of all the added linguistic terms are zero.
To effectively measure the deviations and closeness of different PLTSs assessments, Zhang [58] presented the probabilistic linguistic distance measure which is introduced as follows.
The PLTS is a useful tool to model the uncertainty in qualitative MCDM.Benefiting from the development of information techniques, the decision made by a large number of DMs or participants has become more frequent in recent years.The decision making problems in which the decision data are represented by PLTSs are more and more common nowadays, which have recently received increasing attention from many researchers and practitioners [62][63][64][65].For example, Zhang et al. [64] developed a consistency-based risk assessment approach to handle the group decision making problem with probabilistic linguistic preference relation.Lin et al. [63] developed a probabilistic uncertain linguistic TOPSIS approach to solve the selection problem of cloud storage services.

Proposed Probabilistic Linguistic VIKOR Approach
To solve the GSC initiatives evaluation problem in which the evaluation values are expressed in PLTSs, we next develop a probabilistic linguistic VIKOR multi-criteria analysis approach, which is motived by the classical VIKOR method.Firstly, the evaluating framework of the green supply chain initiatives is presented in Section 4.1.According to the idea of the classical VIKOR approach, there are two key issues to be addressed in the developed method, i.e., identifying the probabilistic linguistic ideal solutions which are conducted in Section 4.2 and calculating the closeness of each solution to the ideal one by a series of probabilistic linguistic measures which is implemented in Section 4.3.Finally, the procedures of the developed method are introduced in Section 4.4.

The Evaluating Framework of the Green Supply Chain Initiatives
In practice, the performance evaluation of GSC initiatives usually requires the firms to take many conflicting criteria into account, such as the inventory level, the suppliers and the materials, etc., which should be considered purchasing-wise; the technical capability, the innovation capability and the others should be considered manufacturing-wise; the inbound and outbound logistics, and the packaging, etc., should be considered logistics-wise; and the saleability, marketability and the others should be considered marketing-wise.The framework of the evaluation problem of the GSC initiatives introduced by the studies [2,28] is shown in Figure 1.
It is easily observed from Figure 1 that this evaluation problem is a hierarchical MCDM problem with a two-layered structure.The first layer includes four main-criteria {c 1 , c 2 , c 3 , c 4 }.The weights of the main-criteria are denoted by {w(c 1 ), w(c 2 ), w(c 3 ), w(c 4 )} which are represented by PLTSs in this study.Every main-criterion c j (j ∈ {1, 2, 3, 4}) includes #c j sub-criteria {c j(1) , c j(2) , • • • , c j(#c j ) } in the second layer, where #c j denotes the number of sub-criteria in the main-criterion c j .The weights of the sub-criteria {c j(1) , c j(2) , • • • , c j(#c j ) } are represented by {w(c j(1) ), w(c j(2) ), • • • , w(c j(#c j ) )} which also take the form of PLTSs.Let x ij(k) be the criterion value of the alternative a Then, this evaluation problem is mathematically expressed as a probabilistic linguistic decision matrix R which is shown in Table 1.
others should be considered marketing-wise.The framework of the evaluation problem of the GSC initiatives introduced by the studies [2,28] is shown in Figure 1.   .Then, this evaluation problem is mathematically expressed as a probabilistic linguistic decision matrix R which is shown in Table 1.

To Determine the Probabilistic Linguistic Ideal Solutions
To effectively solve the previous decision making problem, we first need to identify the probabilistic linguistic positive ideal solution and the probabilistic linguistic negative ideal solution among all alternatives over each sub-criterion, respectively.To do so, a new comparison method is next developed for comparing the magnitude of PLTSs.

To Determine the Probabilistic Linguistic Ideal Solutions
To effectively solve the previous decision making problem, we first need to identify the probabilistic linguistic positive ideal solution and the probabilistic linguistic negative ideal solution among all alternatives over each sub-criterion, respectively.To do so, a new comparison method is next developed for comparing the magnitude of PLTSs.Proof.
Based on the expectation values and variance values of PLTSs, the comparison law for PLSTs is introduced as follow.
Theorem 1.Given a PLTS L(p), let L(p) be the normalized form of L(p); then the following statements hold: The proof of Theorem 1 is straightforward.According to the definitions of the expectation and variance of PLTSs, it is easy to show that statements (1) and ( 2) in Theorem 1 hold.Based on statements ( 1) and ( 2) in Theorem 1, statement (3) in Theorem 1 holds apparently.
Using the above comparison method of PLTSs, the positive ideal solution (L j(k) (p)) + among all alternatives over the sub-criterion c j(k and the negative ideal solution (L j(k) (p)) − is identified as follows: Therefore, the probabilistic linguistic positive ideal solution for this kind of probabilistic linguistic MCDM problem is denoted by: and the probabilistic linguistic negative ideal solution is represented by:

Measuring the Closeness Indices of Alternatives to the Ideal solutions
Next, the distance measure of PLTSs developed by Zhang [58] is used to measure the closeness index to the ideal solution.Before that, we need to defuzzy the weights of the main-criteria and the weights of their sub-criteria because they are represented by PLTSs in our work.To this end, we first develop a defuzzification function of PLTSs.

Definition 8. Given a set of PLTSs
where Exp(L i (p)) is the expectation values of L i (p and the parameter ∂ is to balance the expectation values and the variance values in decision making process, which should be bigger than 0.5 according to the developed comparison law of PLTSs. Furthermore, the normalized defuzzification values of the weights of the main-criteria w(c j ) and the sub-criteria w(c j(k) ) are obtained as follows: As a result, the overall weight of the sub-criterion c j(k) , which is denoted by ij , is the product of the normalized defuzzification values of the weight of this sub-criterion c j(k) and the normalized defuzzification values of the weight of the corresponding main-criterion c j .
Then, the probabilistic linguistic group utility measure of the alternative ) is defined by the following equation: Furthermore, the probabilistic linguistic individual regret measure of the alternative a i (i = 1, 2, • • • , m) is introduced by the following equation: According to the idea of the classical VIKOR method, the compromise solution is usually a maximum group utility for the majority and a minimum individual regret for the opponent, simultaneously.To this end, the compromise measure of the alternative a i (i = 1, 2, • • • , m) is defined as below: where the parameter η is the weight of the strategy of the majority of criteria or the maximum overall utility, and is usually taken as 0.5.According to the decreasing order of the values of PLS i , PLR i , PLQ i (i = 1, 2, • • • , m), three ranking lists of alternatives are obtained.The compromise solution for the MCDM problem is the alternative a σ(1) (it is the alternative with the first position in the ranking list derived by The alternative a σ( 1) is also the alternative with the first position in the ranking lists derived by PLS It is common that these two conditions are not usually satisfied simultaneously in the practical MCDM process.If Condition C 1 is not satisfied, the set of alternatives {a σ(1) , a σ(2) , • • • , a σ(µ) } is the compromise solution in which the maximum value of µ is determined by the following equation: PLQ(a σ(µ) ) − PLQ(a σ(1) ) < 1/(m − 1); and if Condition C 2 is not satisfied, then the set of alternatives {a σ(1) , a σ(2) } is the set of compromise solutions.

Procedures for the Developed Method
Based on the above analysis, the schematic diagram of the proposed method for solving the MCDM problem with PLTSs is provided in Figure 2, and the decision procedures of the proposed method can be summarized as follows: Step 1. Formulate the MCDM problem, identify the probabilistic linguistic decision matrix, and determine the probabilistic linguistic weights of the main-criteria and sub-criteria.
Step 2. Determine the probabilistic linguistic positive ideal solution and the probabilistic linguistic negative ideal solution by employing Equation ( 16) and Equation ( 17), respectively.Step 3. Calculate the overall weight of each sub-criterion by using Equation (20).

A Case Study for the Evaluation of Green Supply Chain Initiatives
In this section, a real-world MCDM problem involved with the selection of an appropriate time window for a fashion company to implement green raw material is introduced to demonstrate the decision making process and usefulness of the developed probabilistic linguistic VIKOR approach.

Description
To improve saleability and secure future growth in the wide market, more and more fashion retail firms intend to employ new green materials in their products.The use of new green materials in turn requires changes in firms' operational capabilities and resources that may have an adverse effect on firms' operations performance [2].Thus, the firms need to evaluate their performances carefully when implementing GSC initiatives.In this section, we consider a real-life MCDM problem involved with the selection of an appropriate time window for a fashion company to implement

A Case Study for the Evaluation of Green Supply Chain Initiatives
In this section, a real-world MCDM problem involved with the selection of an appropriate time window for a fashion company to implement green raw material is introduced to demonstrate the decision making process and usefulness of the developed probabilistic linguistic VIKOR approach.
Then, the normalized defuzzification values of the weights of the main-criteria and the corresponding sub-criteria are obtained by using Equation (19).These results are illustrated in Table 5.Furthermore, the overall weights of the sub-criteria under each main-criterion are obtained by employing Equation (20) and are also shown in Table 5. Afterwards, the probabilistic linguistic group utility measures of four alternative implementation time windows are calculated by using Equation (21) as below: PLS 1 = 0.7440, PLS 2 = 0.5204, PLS 3 = 0.2626, PLS 4 = 0.6742, and the probabilistic linguistic individual regret measures of these four alternative implementation time windows are computed by using Equation ( 22) as follows: PLR 1 = 0.0921, PLR 2 = 0.1535, PLR 3 = 0.0898, PLR 4 = 0.1228.
Using Equation ( 23), the probabilistic linguistic compromise measures of these four alternative implementation time windows are calculated as follows: According to the decreasing order of the values of PLS i , PLR i , PLQ i (i = 1, 2, 3, 4), three ranking lists of these four alternative implementation time windows are obtained as follows: Obviously, the compromise solution for the above assessment problem is alternative a 3 (i.e., Implement in 8 months) since it satisfies the following two conditions: (1) condition C 1 , i.e., PLQ(a 3 ) − PLQ(a 1 ) = 0.5180 ≥ 1 4−1 (= 0.3333); (2) condition C 2 , i.e., alternative a 3 is also the alternative with the first position in the ranking lists derived by PLR i (i = 1, 2, 3, 4).

Discussion
The results derived for the developed probabilistic linguistic VIKOR technique show that implementing the green raw material in 8 months (a 3 ) is the appropriate time window for the fashion company and should be recommended among the four possible time windows.The main reason is that there are potential gaps in capability and resources in its supply chain in order to successfully implement the green raw material now [2].It is easily noticed that compared with implementing now (a 1 ), the company has made full preparations in terms of manufacturing processes, production capacity, and technical and innovation capabilities when implementing green raw material in eight months, and based on a comparison with implementing in twelve months (a 4 ), the company has marketing advantages over its rival competitors because few competitors have already launched similar green initiatives.Therefore, the result obtained by the probabilistic linguistic VIKOR method is reasonable.Comparing the previous studies [2,12] on the evaluation of GSC initiatives, the developed probabilistic linguistic VIKOR technique allows the evaluators to employ PLTSs to express the imprecise performances of the GSC initiatives with respect to various criteria, which greatly improves the elicitation of linguistic information.Moreover, it presents a probabilistic linguistic compromise solution which is a maximum probabilistic linguistic group utility for the majority and a minimum probabilistic linguistic individual regret for the opponent.Although the proposed probabilistic linguistic VIKOR technique is used in this paper to assist a fashion retail company for selecting an appropriate time window to implement green raw material, this proposed technique can also be used by different companies in other industry sections, as it is extremely flexible and can be reduced into various special models according to the different actual needs.

Conclusions
The real-world GSC management problems in which the decision data are represented by PLTSs are more and more common nowadays.To adequately deal with this kind of GSC management problem with PLTSs, in this paper, we have developed a new probabilistic linguistic VIKOR method.Using this developed method, the managers of a fashion retail firm have successfully selected a suitable time window to implement green raw material in order to achieve the goals of sustainable economic and environmental protection.The main contributions of this paper are summarized as follows: (1) a new comparison method of PLTSs has been presented to effectively identify the ideal solutions for the probabilistic linguistic MCDM problems; (2) a new defuzzification function of PLTSs has been proposed to fully take into account the weights of criteria under probabilistic linguistic contexts; (3) a new probabilistic linguistic VIKOR approach has been developed to effectively evaluate various challenges that the firms usually face when implementing GSC initiatives.
Nevertheless, this study also contains several limitations which may serve as suggestions for further research.First, we assumed that all criteria are independent in the developed approach; further study can consider integrating the Choquet integral into the proposed method to take the dependency between criteria into account.Meanwhile, with the rapid development of human social-economic activities, some new and important main-criteria and sub-criteria should be discussed and added in the developed approach because the evaluation main-criteria and sub-criteria used in this paper were adopted in [2].Last but not least, further research should consider developing a decision support system based on the proposed approach to facilitate the managers of the firms to make decisions quickly and rationally.

Figure 1 .
Figure 1.The framework of the performance evaluation of GSC initiatives.
the second layer, where # j c denotes the number of sub-criteria in the main-criterion j c .The weights of the sub-criteria

Figure 1 .
Figure 1.The framework of the performance evaluation of GSC initiatives.

Table 1 .
The probabilistic linguistic decision matrix.

Step 4 . 19 Step 6 .
Compute the probabilistic linguistic group utility measures of alternatives by using Equation (21).Step 5. Calculate the probabilistic linguistic individual regret measures of all alternatives by using Equation (22).Step 6. Calculate the compromise measures of alternatives by employing Equation (23).Step 7. Identify the compromise solution(s) according to the values of PLS i , PLR i , PLQ i .Sustainability 2017, 9, 1231 12 of Calculate the compromise measures of alternatives by employing Equation (23).Step 7. Identify the compromise solution(s) according to the values of , ,

Figure 2 .
Figure 2. Schematic diagram of the proposed probabilistic linguistic VIKOR technique.

Figure 2 .
Figure 2. Schematic diagram of the proposed probabilistic linguistic VIKOR technique.

Table 1 .
The probabilistic linguistic decision matrix.

Table 5 .
The defuzzification values of the main-criteria weights and the sub-criteria weights.