An Integrated Decision-Making Approach for Green Supplier Selection in an Agri-Food Supply Chain: Threshold of Robustness Worthiness

Along with the increased competition in production and service areas, many organizations attempt to provide their products at a lower price and higher quality. On the other hand, consideration of environmental criteria in the conventional supplier selection methodologies is required for companies trying to promote green supply chain management (GSCM). In this regard, a multi-criteria decision-making (MCDM) technique based on analytic hierarchy process (AHP) and fuzzy technique for order of preference by similarity to ideal solution (TOPSIS) is used to evaluate and rate the suppliers. Then, considering the resource constraint, weight of criteria and a rank of suppliers are taken into account in a multi-objective mixed-integer linear programming (MOMILP) to determine the optimum order quantity of each supplier under uncertain conditions. To deal with the uncertain multi-objectiveness of the proposed model, a robust goal programming (RGP) approach based on Shannon entropy is applied. The offered methodology is applied to a real case study from a green service food manufacturing company in Iran in order to verify its applicability with a sensitivity analysis performed on different uncertainty levels. Furthermore, the threshold of robustness worthiness (TRW) is studied by applying different budgets of uncertainty for the green service food manufacturing company. Finally, a discussion and conclusion on the applicability of the methodology is provided, and an outlook to future research projects is given.


Aim and Originality of the Study
In a competitive environment, supplier selection is one of the most important problems that manufacturing firms encounter. In some industries, the cost of raw materials constitutes a high portion of the final production cost that can be minimized by selecting the appropriate supplier [1]. The supply chain is defined as a network of units for producing and delivering products and services that starts with the supplier or suppliers and ends with the warehouses and buying centers. Supply chain management (SCM) is a set of methods applied for effective integrating of suppliers, producers, warehouses and buying FANP, fuzzy technique for order of preference by similarity to ideal solution (TOPSIS) and a fuzzy MOLP model for supplier selection problem in a green supply chain. A comprehensive eco-friendly technique was introduced by Kumar et al. [21] for supplier selection. They proposed a data envelopment analysis (DEA) with carbon footprint monitoring to consider heterogeneous suppliers and regional emission compliance standards and laws. Bakeshlou et al. [22] addressed the problem of supplier selection using a multi-objective decision-making algorithm. They applied a hybrid fuzzy multi-objective decision-making (MODM) algorithm for a green supplier selection (GSS) problem considering 17 criteria. Rueda et al. [23] proposed a method that was a synthesis of the widely applied instruments applied by agri-food companies for the sustainability promotion of their supply chains. This method also provided an analytical framework to elucidate how the decisions are made. This framework was based on the competitive environment in which firms operate with respect to the location of the raw material, technologies of the suppliers, leverage over upstream suppliers and the characteristics of the end-markets. Shi et al. [24] introduced a novel technique for GSS with interval-valued intuitionistic uncertain linguistic (IVIUL) information. They applied grey relational analysis (GRA) and TOPSIS methods for the assessment and selection of green suppliers. Finally, a case study in an agri-food industry was conducted to demonstrate the applicability of their proposed methodology.
As relevant research carried out with emphasis on operations research (OR) methods for a sustainable agri-food supply chain, Allaoui et al. [25] presented a two-stage hybrid methodology. On the first stage, they implemented AHP and ordered weighted averaging (OWA) aggregation method. Based on the obtained results at the first stage, a multiobjective mathematical model was developed at the second stage to optimize the design of the supply chain network. Banaeian et al. [13] compared the application of three popular multi-criteria supplier selection methods in a fuzzy environment. The incorporation of fuzzy set theory into TOPSIS, VlseKriterijuska Optimizacija I Komoromisno Resenje (VIKOR) and GRA methods was thoroughly discussed. The methods were then utilized to complete a green supplier evaluation and selection study for an actual company from the agri-food industry. Their comparative analysis for the case study indicated that the three fuzzy methods arrive at identical supplier rankings, while fuzzy GRA requires less computational complexity to generate the same results. Gören [26] presented a decision framework for sustainable supplier selection and order allocation in a case study problem. Their framework consists of three integrated components. First, a fuzzy decision-making trial and evaluation laboratory (DEMATEL) helps to identify the weights of sustainable criteria. Second, these weights are considered as inputs in Taguchi Loss Functions for ranking and calculating the rating value of each supplier. Finally, the ranking values were treated as inputs of a bi-objective optimization program to determine the optimal order quantities that should be allocated to each supplier. Liu et al. [27] suggested a novel two-stage integrated model for a GSS problem of fresh product. At the first stage, they took function development (QFD), customer requirement (CR) and company strategy (CS) into account as selection criteria. A fuzzy best-worst method (BWM) was applied to obtain the criteria weights. To rank the suppliers, a fuzzy multi-objective optimization by ratio analysis (MOORA) plus the full multiplicative form, namely, MULTIMOORA was then employed. Recently, Haeri and Rezaei [28] introduced a grey-based GSS model to cover uncertain environments. They applied a novel weight assignment model using the combination of BWM and fuzzy grey cognitive map to obtain the interdependencies among the economic and environmental criteria. They could finally demonstrate the efficiency of the proposed method by analyzing a real case study. Recently, Tirkolaee et al. [29] proposed a hybrid approach based on fuzzy MCDM and weighted goal programming (WGP) to tackle a sustainable-reliable supplier selection problem. They investigated a real case study problem in the lamp supply chain and demonstrated the efficiency of their methodology to determine the optimal policy. Keshavarz-Ghorabaee et al. [30] offered an integrated decision-making approach based on Fermatean fuzzy sets, weighted aggregated sum product assessment (WASPAS) and the simple multi-attribute rating technique (SMART) for evaluating and selecting green construction suppliers under uncertainty. Table 1 provides a comparison of similar supplier selection problems in order to highlight the contributions research has made.

Contributions and Structure of the Manuscript
According to the literature, only a few studies have been practically developed in which there were hybrid methods of economy and GSS criteria, and attention to the order allocation and the supplier selection problem in a green agri-food supply chain. Moreover, the robust optimization technique was not investigated in the literature on our subject. Thus, the innovation of our paper is summarized as follows:

1.
Applying a AHP-fuzzy TOPSIS method for rating the suppliers in a green agri-food supply chain network; 2.
Developing a robust goal programming (RGP) model to determine the order allocation under uncertain condition; 3.
Introducing and investigating the threshold of robustness worthiness (TRW); 4.
Investigating the applicability of the proposed methodology through an illustrative case study.
The rest of the paper is organized as follows. Sections 2 and 3 present the proposed methodology to design a green agri-food supply chain design. An AHP-fuzzy TOPSIS method is introduced in Section 2, followed by the RGP-based mathematical model in Section 3 to deal with uncertainty and multi-objectiveness. The computational results of the case study problem are explained in Section 4. A discussion of the results is provided in Section 5 and, finally, Section 6 concludes the work with a note about future research directions.

Methodology
In this section, the proposed integrated approach for GSS is presented based on AHPfuzzy TOPSIS and RGP. The analytic hierarchical process (AHP) is used for the calculation of relative weights of the supplier selection criteria; the fuzzy TOPSIS is used for rating the suppliers according to the selection criteria; finally, the weight of each supplier is identified. Then, the problem is formulated as a bi-objective model based on real-world assumptions and limitations, in order to determine the optimal values of the decision variables, which are described in Section 3. To cope with two objectives of the problem under uncertainty of the demand parameter, an RGP technique is applied to the model. This technique guarantees feasibility, robustness and optimality of the solutions [37]. Figure 1 depicts the proposed method.

Analytic Hierarchical Process (AHP)
The AHP technique helps decision-makers to identify the alternative that best befits their goals. It was first introduced by Saaty [38] as a measurement theory to provide pairwise comparisons, which depend on expert opinions.
AHP is implemented in four decision-making phases: 1. designing the hierarchy, 2. performing the pairwise comparisons, 3. calculating the weights, and system consistency.

Phase 1: The Hierarchy Tree Construction
The phase 1, important criteria of decision making must be hierarchically expressed in the form of a decision tree where each criterion can include several sub-criteria.

Phase 2: Pairwise Comparisons
The criteria of each level of the AHP are compared with the ones of the upper level in a pairwise manner, and their weights are calculated. They are called relative weights. These weights are then combined and identify the final weight of each solution. These comparisons are given by decision makers as linguistic variables; i.e., if criterion i is compared to criterion j, the decision maker expresses the importance of i over j as one of the following choices based on Saaty's nine-point scale (cf. Table 2): Table 2. Comparison criteria according to the relative importance degree.

Relative Comparison of Criterion i Against j (Regarding the Intended Goal)
Relative Importance Degree (Score) Equal importance 1 Weak priority of i against j 3 Strong priority of i against j 5 Very strong priority of i against j 7 Absolute priority of i against j 9 Now, a pairwise comparison matrix is generated employing Saaty's nine-point scale. Considering n criteria, we can construct a square matrix (A applying the pairwise comparison of criterion i with criterion j), which is shown in Equation (1), where a ij represents the importance of ith criterion against jth criterion. We have a ij = 1, when i = j, and a ji = 1

Phase 3: The Extraction of Weights from the Pairwise Comparison Matrix
The weight extraction method depends on the consistency or inconsistency of the decision matrix. Since decision makers are not always consistent, AHP technique permits some small opinion inconsistency. If the condition represented in Equation (2) holds, then the matrix is called consistent: a ij × a jk = a ik ; ∀i, j, k. ( For further information, please see Saaty and Vargas [39].

Similarity to the Ideal Solution
In the real world, due to incomplete or inaccessible information, data are not deterministic; they are, as we say, fuzzy data. Actually, in many issues, unclear and unrealistic concepts are observed in the decision-making data; deterministic data are not sufficiently appropriate in real conditions. Therefore, we applied TOPSIS with fuzzy data to prioritize the potential suppliers. In this method, the ranking of solutions and the criteria relative weights have been given by linguistic variables whose corresponding fuzzy numbers are presented in Table 3.  (5,7,9) High (H) (7,9,10) Very High (VH) (9,10,10) In this technique, the decision matrix contains m rows and n columns. Subsequently, six steps of the fuzzy TOPSIS technique are outlined.

Step 0: Creating the Decision Matrix
Considering the number of criteria and solutions, the decision matrix is created as follows: x ij is the performance of the alternative i (i = 1, 2, . . . , m) in relation to the criterion j (j = 1, 2, . . . , n). In this matrix, shown in Equation (3), the criterion that has a positive utility is the profit criterion, and the criterion that has a negative utility is the cost criterion. Furthermore, the importance degree of each criterion in the decision making has to be identified in this step, which is represented by W = (w j1 , w j2 , w j3 )-a fuzzy factor. Now, the fuzzy decision matrix ( D) is defined as follows, which comes from the potential use of eigenproblem [40,41]: where X j denotes the fuzzy performance vector of all solutions for criterion j, and A i represents the fuzzy performance vector of all criteria for alternative i.

Step 1: Normalizing the Decision Matrix
In this step, the values of the decision matrix become unscaled by normalization. Consequently, the normalized fuzzy decision matrix ( R) is defined: with respect to r ij = and a − j = min(a ij ) for negative criteria, where a ij , b ij , c ij are three elements of the triangular fuzzy number for rating alternative i with respect to criterion j.

Step 3: Finding the Fuzzy Positive Ideal Solution and Fuzzy Negative Ideal Solution
First, M( v ij ) is used to compare fuzzy numbers-the mean of removals is calculated in the following way: The maximum value of M( v ij ) is selected as the fuzzy positive ideal solution (FPIS) for the jth criterion V + j ; it is represented as

Step 4: Calculating the Distance of Solutions from the FPIS and FNIS
The distances of solutions from the FPIS and FNIS or fuzzy positive ideal reference point (D + ij ) and fuzzy negative ideal reference point (D − ij ) are calculated according to Equations (7) and (8): Here,(D + ij ) and (D − ij ) are deterministic numbers. The distances of the solution i from the positive ideal solution (S + i ) and negative ideal solution (S − i ) are determined by Equations (9) and (10), respectively: 2.2.6.
Step 5: Calculating the Similarity Index The similarity index or closeness coefficient of each alternative (C i ) is defined by: 2.2.7.
Step 6: Ranking the Solutions The calculated similarity indices are ranked in descending order.

Robust Goal Programming (RGP)
The RGP approach is implemented to deal with the uncertainty and bi-objectiveness of the model. To the best of our knowledge, RGP has not been investigated enough in the literature while it is regarded as an efficient exact solution approach [37,42]. In other words, it consists of two efficient methods; i.e., GP that causes the constraints to be less strict and the solution space to be larger [43], and robust optimization that protect the solutions against worst-case realizations of the uncertain parameters [44]. Hence, we apply RGP to the model in order to generate robust solutions and guarantee the feasibility of the solutions under uncertain conditions.

The Robust Goal Programming Mathematical Model
The aim of this stage is to design a novel mathematical model based on the real-world assumptions in order to determine the optimal allocation of raw material to the suppliers within a green agri-food supply chain design so that total cost (cf. Equation (21)) and total weighted purchase (cf. Equation (22)) are minimized and maximized, respectively. This model is applied to allocate orders to the selected suppliers that have been prioritized according to the proposed technique in Section 2.
Here, a two-objective linear programming model is proposed for the supplier selection problem with multiple sourcing that minimizes the final purchasing costs and maximizes the total purchasing amount from the high-ranked suppliers considering their priority values. Of course, this model includes a set of constraints such as purchasing order, supplier capacity and assumptions applied in formulating the multi-criteria linear programming model. The number of items allocated to a supplier is also taken into account. The order quantity of each item from a supplier is constrained by the supplier capacity, too.

Interval-Based Robust Optimization
There are different methods to deal with the uncertainty. In this paper, the intervalbased robust optimization method based on the Bertsimas and Sim model is used, because it offers a linear model that outperforms the others [45,46]. Accordingly, the demand parameter (D r ) is considered uncertain; i.e., it is defined as the interval [D r − D r , D r + D r ] in which D r is the average demand of item r, and D r is the deviation rate of demand of item r. Here, D r is defined according to the uncertainty level (ψ), determined by D r = ψD r . This parameter leads to uncertainty in the maximum acceptable defect rate and constraint and the calculation equation of the order quantity. The corresponding equations with uncertainty are: where q rs and Q r represent the average defect rate of item r from supplier s and maximum acceptable defect rate of item r, respectively. Then, the parameter Γ (uncertainty budget) is defined for the rth item such that Γ r ∈ [0, |J r |], where J r represents the set of coefficients in row r of the constraints; Equations (12) and (13) [46]. In fact, it ensures the robustness adjustment against the protection level of the solution. Bertsimas and Sim [46] demonstrated that it is unlikely that all of the coefficients change concurrently. In fact, a subset of the coefficients is allowed to have an impact on the solution adversely. On the other hand, due to the symmetric distribution of variables, even if more than the parameters change, the obtained robust solution will remain feasible with a very high probability. Therefore, this is called the conservation level for each item. The parameter may be an integer or not. If so, then, the problem will be deterministic. On the other hand, if Bertsimas' and Sim's approach is transformed to Soyster's method, then the most degree of conservation [46] is obtained. Thus, the fluctuation of the robustness is adjusted against the conservatism level of the solution. Finally, the proposed linear robust form of Equation (12) is [47]: where Γ 1 r denotes the uncertainty budget associated with the demand uncertainty of rth item in Equation (12), which can take value in [0, 1].
Similarly, Equation (13) would be robust as follows [48]: where Γ 2 r denotes the uncertainty budget associated with the demand uncertainty of rth item in Equation (13), which can take value in [0, 1]. Now, a robust mixed-integer linear goal programming (RMILGP) model is proposed to formulate the problem, based on the suggested robust model by Bertsimas and Sim [46].

Objective Functions
Equations (16) and (17) represent the 1st and 2nd objectives of our model, respectively. As is clear, the 1st objective is to minimize the total cost of the proposed green supply chain, and the 2nd objective is to maximize the total purchasing value obtained by different suppliers and food items.

RGP Model
x rs ≤ c rs y rs , (r = 1, 2, · · · , R; s = 1, 2, · · · , S) x rs ≥ c rs y rs , (r = 1, 2, · · · , R; s = 1, 2, · · · , S) x rs , l r ≥ 0, y rs ∈ {0, 1}, (r = 1, 2, · · · , R; s = 1, 2, · · · , S) Equation (18) denotes the main objective of the RGP model with the aim of minimizing the weighted sum of positive deviation from the 1st objective (cf. Equation (15)) and negative deviation from the 2nd objective (cf. Equation (16)). Since the 1st objective is a minimization function and the 2nd objective is a maximization one with different units, each related part in the RGP objective is divided by its ideal value in order to be normalized. Equation (19) considers the total cost that consists of the purchasing cost, ordering cost, transportation cost, holding cost and inventory shortage cost, which is written as a constraint using RGP. Objective (Equation (20)) focuses on the total weighted purchasing value, which is written as another constraint using RGP. Note that TCP and TVP are calculated by solving the single-objective model by 1st and 2nd objectives separately. Equation (21) is a robust constraint that limits the average defect rate of raw materials to the maximum acceptable defect rate. Equation (22) is another robust constraint with the same notations of robustness defines that the order quantity of item r must be equal its total demand minus its quantity that is not fulfilled. Equation (23) guarantees that the total order quantity of item r from supplier s has to not exceed the capacity of item r. Equation (24) reflects that the order quantity of item r from supplier s must not be less than the minimum capacity of supplier s (or, due to economic considerations, the order quantity of item r from supplier s must be higher than the minimum capacity of supplier s). Equation (25) indicates the type of the decision variables.

Computational Analysis
This section presents the results of a real case study in order to evaluate the applicability of the offered methodology; through supplier evaluation and selection process, the aim is to provide the Green Service Food Company in Iran with a clear image of the structure of activities performed by the purchasing and supply manager, with the hope that it can improve the economic performance and efficiency of the company. Investigating the list of effective factors on the supplier evaluation and selection process and interviewing the relevant experts of raw food material supply and food industry specialist, 10 major criteria of price, transportation cost, flexibility, technology, quality, on-time delivery, failure, order fulfillment, green competencies and environmental management system have been identified. It should be noted that the first eight criteria are regarded as conventional criteria [13], and the last two criteria are taken into account as the most important environmental criteria [49]. Moreover, we consider five potential suppliers for each item. The obtained values for TVP and TCP are 1,106,279.309 and 201,057.224; both are considered as the main determinative input parameters to the model.

Supplier Ranking
The main criteria to supply each item were extracted through interview with the Quality Assurance (QA) expert of the company to determine the importance degree of each criterion and the performance of each potential supplier in relation to each criterion. His opinion was collected in linguistic terms, and then was converted into triangular fuzzy numbers (see Table 2). The hierarchy of the supplier selection is represented in Figure 2. Furthermore, the final ranking of the suppliers is given in Table 4.

Raw Food Supplier Ranking Overall Weight of Suppliers
The details of the supplier ranking calculations are provided in Appendix A.

Model Results
After prioritizing the suppliers and calculating the importance of each supplier (W r ), the proposed mathematical model has been solved using the CPLEX solver of GAMS software. This model determines the optimum allocation order of each item and finds out the best solution for the allocation order problem. It should be noted that the uncertainty level (ψ) is determined by the market expert to be 0.1, and the conservatism levels (Γ 1 r and Γ 2 r ) are considered to be 0.5 for each item. The advantages of the 1st objective (W 1 ) and the 2nd objective (W 2 ) are determined using the entropy approach proposed by Shannon [50]. Shannon entropy is only a particular, limit case of Renyi entropy [51] and has important applications regarding hierarchical tree structures [52]. It is an appropriate scale to be implemented on different cases of assessment in various decision-making processes [53]. Moreover, it can be employed to calculate the quantity of useful information presented by data itself. We take into account three combinations based on the market, production and QA experts as alternatives. The values of the 1st and 2nd objectives are considered as two evaluation criteria, thus the decision-making matrix is defined in Table 5. Now we need to normalize the above 3 × 2 matrix, the result is shown in Table 6. In the following, the entropy of each criterion (E j ) should be calculated according to Equation (26).
where i and j denote the criterion and alternative indexes, k is a fixed value and p ij represents the values in the normalized matrix (cf. Table 6). Now, the deviation degree of criteria (dev j ) is calculated: In the final step, the assigned weight to each criterion (wg j ) is obtained: The obtained values for E j , dev j and wg j are given in Table 7. Now, we should choose the best alternative considering the assigned weights of criteria. Since the 1st criterion (1st objective function) is a negative criterion (the lower value it has, the more desirable it is) and the 2nd criterion (2nd objective function) is a positive criterion (the larger the value, the more desirable), the following equation is defined to determine the score of each alternative: Table 8 presents the score of each alternative. As can be seen, the 3rd alternative is chosen. By implementing the mathematical model, the allocation values have been determined; they represent the demand and the purchasing quantity of each item. Furthermore, the purchasing quantity of each item and the quantity of not fulfilled items are reported in Table 9. The obtained values for the RGP, 1st and 2nd objective functions with deviations values are shown in Table 10. In the following, a sensitivity analysis is conducted on the uncertainty levels of 0, 0.05, 0.1, 0.15 and 0.2 to study the instability of the uncertainty and address specificity and accuracy [54,55]. It should be noted that the uncertainty level of 0 leads to a certain (deterministic) model that is not robust anymore. Under these different conditions, the obtained results are compared with each other, which would help managers to find the optimal policy at each level of uncertainty. To this end, the obtained values for the 1st and 2nd objective functions are illustrated in Figures 3 and 4. Note that the trade-off between objectives is dependent on the purchasing quantity. It is expected that, with the increase in purchasing quantity, the 1st and 2nd objectives also increase. However, a larger uncertainty level may violate this trade-off. This would result in more shortage and, consequently, higher penalty costs, which increase the 1st objective and decrease the 2nd objective.
According to Figures 3 and 4, the obtained results for ψ = 0 show the status of the problem in a certain condition, which has the best values for both objectives. Under uncertain conditions, which range from ψ = 0.05 to ψ = 0.2, the status is becoming worse and worse. These deviations from the best values are the imposed costs to remain feasible in different uncertain conditions. As is obvious, the maximum value of the 1st objective and the minimum value for the 2nd objective are obtained for ψ = 0.2. Therefore, managers should pay more attention, consider more resource and have a secondary plan in order to overcome these increases in total cost and shortages in fulfilling the required demands. These efforts may be made after analyzing and evaluating the behavior of the objectives under different conditions.

Robustness Acceptance Threshold
One of the most effective factors in robust decision-making is to determine to what extent it is valuable to make the problem robust. The main factor is to consider the budget constraint of the organization so that the budget value would specify how much the robustness is possible and valuable or not. To this end, according to the market expert's opinion, we consider a fixed uncertainty level of 0.1 and analyze the problem for different values of conservatism levels to find the threshold of robustness worthiness (TRW). Moreover, the budget limitation of USD 1, 600, 000 is applied to the proposed robust model. Table 11 shows the total cost against different combinations of conservatism levels, and Figure 5 highlights the threshold. As is clear from Table 11 and Figure 5, the approximate value of 0.785 is the recommended threshold value for both conservatism levels. The budget limitation would be violated for a greater value than 0.785, and the robustness becomes unworthy due to problem's infeasibility.

Discussion
The numerical results provided an optimal plan to the management of Green Service Food Company under uncertain demands. First, the priority of the potential suppliers for each raw food material was determined to identify the superior suppliers. Based on real data, the optimization process was then performed to determine the total cost considering the maximization of the total purchasing value. Accordingly, the optimal supply plan was found including the purchasing quantity of each item from the best possible suppliers. This plan can be regarded as an efficient tool to conduct each supply for raw food materials.
To study the instability of the real world and robustness of the solutions, a sensitivity analysis was performed on the parameter of uncertainty level. It was demonstrated that the total cost increases by 7.85% for an increase of 20% in the uncertainty level. Moreover, the total purchasing value increases by 2.51% for a decrease of 20% in the uncertainty level. On the other hand, an analytical experiment was designed to study the robustness acceptance threshold as an applicable option under uncertain conditions. To this end, a budget limitation was first proposed by the management, and then different tests were performed for different values of conservatism levels (budgets of uncertainty). It was revealed that the 18th breakpoint (Γ 1 r = Γ 2 r = 0.785) exceeds the budget limitation of USD 1, 600, 000. This point can be considered as the critical stage of robustness.
The main contributions of this study were to design a framework addressing the robustness, sustainability and applicability of the green supplier selection problem in an agri-food supply chain. To do so, an integrated methodology was designed based on the gaps in the literature to provide the required requirements. It can inferred that the results obtained by this research can be employed in the future as an applicable managerial tool in decision-making processes of the company under uncertainty.

Conclusions and Future Research Directions
SCM, which includes different tasks, is one of the most important competitive factors for modern companies. At the initial stages, the supplier selection process is a critical step, aiming at selecting the best supplier. In this study, a GSS problem was addressed in an area where there are several raw agri-food materials and suppliers under uncertain condition. The proposed method consists of a novel integration of AHP and fuzzy TOPSIS (AHP-fuzzy TOPSIS) and RGP approach. At the first stage, AHP was employed to identify the criteria weights; then, the fuzzy TOPSIS method was applied to identify the importance degree of each supplier for supplying each item at the second stage. At the third stage, the allocation quantity to each supplier was determined by solving the proposed RMILGP model executed by CPLEX solver and using Shannon entropy method to determine the weights of the objective functions. To validate the proposed methodology, it was implemented by the Green Service Food Company in Iran, and the optimal policy was found out. Practically, this method was developed as a hybrid method under economic and GSS criteria and addressed the allocation of the orders for a multi-resource supplier selection problem in a GSCM context. Finally, a sensitivity analysis was performed on different uncertainty levels to support the managers optimally elaborate on the behavior of the objectives and optimal policies. Moreover, an analysis of different values of the conservatism levels was conducted to investigate TRW considering the budget limitation of the company. Finally, our proposed methodology could provide a robust and sustainable solution to the problem through integrating MCDM and robust optimization. This can be also considered for different variants of the problem under different real-world conditions and it is expected that the same performance will be represented.
As a future suggestion, the proposed methodology is capable of being implemented as a practical framework in similar supply chains such as flower supply chains. To extend the current study based on the main limitations, other uncertainty study methods such as grey systems [56], fuzzy programming [57] and stochastic optimal control [58] can be employed and compared to the current robust optimization approach. Moreover, Paretobased solution methods, e.g., well-known multi-objective algorithms, may be implemented and compared to the proposed RGP. Considering a planning horizon and supply chain reliability maximization as a 3rd objective in the model would make it even closer to real-world status and impact. Eventually, studying the application of electric vehicles in agri-food freight transportation [7] and blockchain technology within the agri-food supply chain [59] can be worthy of investigation to meet the requirements of this era. Acknowledgments: The author would like to thank students Nosek K. and Novak J. for their cooperation on the topic.

Conflicts of Interest:
The authors declare no conflict of interest.

Abbreviations
The following abbreviations are used in this manuscript: The Sets, Indices and Parameters r = 1, 2, · · · , R Index of items; s = 1, 2, · · · , S Index of suppliers; P rs Purchasing cost of item r from supplier s; O rs Order cost of item i from supplier s; t i Transportation cost of item i; d s Distance from the supplier s; H r Holding cost of item r;

W rs
Overall weight (priority value) of the supplier s to supply item r (obtained by AHP-fuzzy TOPSIS method); q rs Average defect rate of item r from supplier s; Q r Maximum acceptable defect rate of item r; D r Uncertain demand of item r; D r Average demand of item r at each period; D r Deviation quantity in the demand of item r; c rs Maximum acceptable capacity of supplier s for item r; c rs Minimum acceptable capacity of supplier s for item r; ψ Uncertainty level; Γ 1 and Γ 2 Conservatism levels of constraints from Equations (12) and (13), respectively; α r Transformation coefficient of weight supplier s to cost; B r Penalty cost for a unit shortage of item r; W 1 Advantage of the 1st objective based on the decision-maker idea; Advantage of the 2nd objective based on the decision-maker idea; TCP Ideal value of the 1st objective; TVP Ideal value of the 2nd objective.

Variables x rs
Purchasing quantity of item r from supplier s; y rs Binary variable; if the item r is supplied from supplier s is 1, otherwise 0; l r Shortage quantity of item r; dl + 1 Positive deviation from the ideal value of 1st objective; dl − 1 Negative deviation from the ideal value of 1st objective; dl + 2 Positive deviation from the ideal value of 2nd objective; dl − 2 Negative deviation from the ideal value of 2nd objective.

Appendix A. The Calculations Related to AHP-Fuzzy TOPSIS
This section provides the calculations related to the proposed AHP-fuzzy TOPSIS step by step. First, different criteria are numbered and then the decision matrix for each item is constructed. Now, the decision matrix for meat is as follows: Table A5. Decision matrix considering suppliers and similarity index for meat.