Selection of Cold Chain Logistics Service Providers Based on a Grey AHP and Grey COPRAS Framework: A Case Study in Vietnam

: Choosing the most suitable cold chain logistics service providers (CLPs) is a vital strategic decision for businesses aiming to achieve an effective and sustainable cold supply chain. A sustainable CLP is one that integrates sustainable practices across its whole operation cycle to achieve product quality, on-time deliveries, and satisﬁed customer requirements, while preventing products from going to waste, which is especially important in the context of a developing country. This study aims to evaluate and select the best CLP regarding their sustainability performance. For this evaluation, a multi-criteria decision making (MCDM)-based framework is proposed that integrates the grey analytic hierarchy process (G-AHP) and grey complex proportional assessment (G-COPRAS) methodologies, in which grey numbers are used to express the linguistic evaluation statements of experts. Initially, the evaluation criteria based on service level, economic, environmental, and social dimensions were determined by means of a literature review and experts’ opinions to employ the MCDM approach. The G-AHP was utilized to identify the criteria weights, and then, G-COPRAS was used to select the best CLP among the alternatives. A case illustration in Vietnam is presented to exhibit the presented approach’s applicability. From the G-AHP ﬁndings, product quality, logistics costs, innovation, and effectiveness of cold chain processes, customer experience, and CO emissions of refrigerated vehicle were ranked as the ﬁve most important criteria. From the G-COPRAS analysis, Yoshida Saigon Cold Logistic (CPL-05) is the best CLP. The robustness of the applied integrated MCDM approach was also tested by conducting a comparative analysis, in which the priority rankings of the best CLPs were very similar. The assessment in this study is directed towards enabling managers, practitioners, and stakeholders of cold chain businesses to assess the most efﬁcient CLP in the supply chain in the market and also to devise suitable strategies toward sustainable development.


Introduction
The term "cold supply chain" refers to a process in which all of the operations involved, from raw material storage through to final product distribution, are temperature controlled. An increase in the volume of imported meat and other related products and exported seafood are the main factors has led to the good performance of cold chain logistics service providers. Due to the instability of the COVID-19 pandemic, more and more customers are choosing to use the online market, boosting the demand for cold storage for food and other various criteria and assisting specialists and managers in balancing and weighing different elements in order to simplify and clarify managerial decisions [6]. Using MCDM techniques allows for a clear understanding of the impact of subjective factors on the final ranking of alternatives and provides insight into the priorities and sensitivities of the various stakeholders. Important ideas for using MCDM approaches include handling complicated challenges with few resources and a poor data system [7]. Moreover, in real-world applications and decision making, the selection is further complicated by the existence of uncertainty which is an unavoidable feature due to the vagueness of human judgments and imprecise information. Unquantifiable, incomplete, and non-accessible information, as well as partial ignorance, are examples of imprecise sources [8]. In this direction, the two major methodologies for incorporating uncertainty and ambiguity into the evaluation process are fuzzy sets theory [9] and grey systems theory [10]. Crisp or conventional methodologies are less effective in handling imprecision or vagueness, while fuzzy sets theory and grey systems theory provide a useful paradigm for analyzing systems with imprecise data and successfully dealing with uncertainty.
In view of the above discussion, the following research objectives for this present study are generated:

•
To understand and identify the evaluation criteria responsible for the feasible and sustainable implementation of the CLP in a developing country's perspective.

•
To determine the most sustainable CLP in Vietnam, a CLP that has employed sustainable practices in its procedures.
In order to achieve these objectives, we apply a grey-based multi-criteria method for the evaluation and selection of CLPs. The method is a two-stage MCDM framework that integrates the grey analytic hierarchy process (G-AHP) and the grey complex proportional assessment (G-COPRAS). A case study in Vietnam is deliberated to exhibit the presented framework's applicability. Initially, evaluation criteria are identified through a literature review and experts' opinions to implement the MCDM approach. The G-AHP is used for obtaining the weights of evaluation criteria. Finally, G-COPRAS is applied for ranking the CLP alternatives.
The remainder of the paper is structured as follows: In Section 2, we provide a literature review for the applications of MCDM methods in the area of 3PL selection, selecting CLPs, and research gaps; in Section 3, primarily, we explain the methodologies applied to the case study attempted in this work; in Section 4, the case illustration is described and the results of the method are discussed; in Section 6, the comparative analysis is conducted to check the robustness of the applied integrated method; in Section 6, we further discuss the managerial implications of the present work; in Section 7, we include concluding remarks as well as recommendations for future study.

Literature Review on 3PL and Selecting CLPs
The selection of third-party logistics (3PL) providers is the most significant issue in any type of business, industry, or supply chain, and numerous studies have been undertaken to develop an effective framework for selecting the best providers in logistics service scopes. Soh [11] employed fuzzy AHP to determine the most potential 3PL provider. The decision making model included four evaluation criteria: investment, service level, communications, and infrastructure, each with its own subcriteria. Ho et al. [12] applied quality function deployment (QFD) and fuzzy AHP to evaluate and select the optimal 3PLs regarding respecting promises, meeting operational and qualitative requirements, creative management, on-time distribution, price, financial stability, error rate, senior manager's availability, and appropriate reaction. To dissect the interrelationships of risks faced by 3PL providers and their customers, Govindan and Chaudhuri [13] employed an MCDM approach called decision making trial and evaluation laboratory (DEMATEL), and concluded that the 3PLs' improvement on internal processes mainly depended on quality management, operational flexibility, and geographic coverage of their services. To achieve well-rounded sustainable 3PLs, Raut et al. [14] indicated the importance of considering ecological, social, and economic dimensions. The authors suggested a hybrid MCDM model that used data envelopment analysis (DEA) and analytic network process (ANP) to assist decision-makers in assessing and obtaining the optimal 3PLs from the perspective of environmental sustainability. Pamucar et al. [15] used the interval rough numbers to handle uncertainty occurring in MCDM problems; the authors presented an integrated interval rough number (IRN) approach based on the best worst method (BWM) and the weighted aggregated sum product assessment (WASPAS) method along multi-attributive border approximation area comparison (MABAC) to evaluate 3PL providers. Jamali and Rasti-Barzoki [16] proposed a game theoretical approach to evaluate 3PLs in a sustainable supply chain by reducing delivery time and carbon emissions. Toward the green logistics outsourcing problem, Vazifehdan and Darestani [17] developed an integrated approach using QFD, fuzzy analytic network process (ANP) method based on fuzzy DEMATEL, and the superiority and inferiority ranking method (SIR) in the petrochemical industry. The 3PL alternatives were evaluated under the most significant factors of flexibility, work experience, quality, green (operation) management of company, logistics green industry, and governments green decision making. Wang et al. [18] developed a hybrid MCDM framework of fuzzy AHP and fuzzy vlsekriterijumska optimizacija i kompromisno resenje (FVIKOR) to select the best 3PL providers in a case study in Vietnam toward sustainability criteria, i.e., service level, economic, environmental, social, and risk dimensions.
In terms of the 3PL selection for cold chain management, Saravanan and Anubama [19] used the AHP to rank CLPs in India's pharmaceutical industry under quality, time delivery, responsiveness, price, flexibility, and technology. Han [20] combined AHP and the technique for order preference by the level of similarity to ideal solution (TOPSIS) method to evaluate CLPs in the fresh food industry, the criteria including logistics level (refrigeration infrastructure, cold chain transportation efficiency, and refrigerated truck utilization), business status (technique level, financial status, and staff level), and service level (customer experience, on-time delivery, and product quality). For managing the value chain of perishable food items or pharmaceutical drugs, known as cold chain management, Singh et al. [21] developed a hybrid model of fuzzy AHP and fuzzy TOPSIS to select an appropriate 3PL to outsource logistics activities of perishable products. Liao et al. [22] implemented the process of evaluating the green levels of CLPs with a hesitant fuzzy linguistic thermodynamic method integrating the cumulative prospect theory and the preference ranking organization method for enrichment evaluation (PROMETHEE) method. The authors used resource, growth ability, performance, and environmental dimensions for the evaluation of the green CLPs. Khan and Ali [23] solved a sustainable supplier selection for the cold supply chain in the context of a developing country, utilizing the interpretive structural modeling (ISM) and fuzzy VIKOR to analyze eight suppliers based on fifteen distinct criteria. Liao et al. [24] conducted a CLP selection in the pharmaceutical industry by a rough set-based gained and lost dominance score method. The authors considered the following criteria for the selection: cost, service performance, infrastructure equipment, and personnel quality.
Based on literature review, the criteria selected for the CLP selection process in this study are summarized in Table 1.

Research Gaps
Based on the preceding discussion and a thorough analysis of the literature, it is deduced that research into the selection of CLPs in the Vietnamese context is limited, with essentially no studies being conducted to propose and ensure sustainable development. With this in mind, attempts have been made in this study to consider the merits of the grey theory, AHP, and COPRAS approaches for CLP selection, and in a broader sense, for 3PL selection. The AHP, widely accepted in MCDM, has been successfully applied to the ranking process of decision making problems. The method's key benefits include its natural capacity to manage intangibles present in any decision making process and assisting decision-makers in organizing the essential components of an issue in a hierarchical structure [32]. The AHP can be integrated with well-known operation research techniques to solve increasingly challenging situations to measure the proportional significance of all criteria. The characteristics of the AHP are flexibility and its capacity to inspect inconsistencies which enables decision-makers to incorporate subjectivity, experience, and knowledge into the decision process [33]. The AHP provides a proven, efficient means of handling difficult decision making and expediting weights' computations.
In addition to the AHP method, the VIKOR, TOPSIS, and COPRAS methods have been increasingly used to quantitatively evaluate complicated economic, social, and environmental processes. The evaluations aim to choose the best alternatives, ranking the alternatives in the order of their significance. Among these, the COPRAS method was firstly introduced by Zavadskas et al. [34] in 1994. This method compares the alternatives and determines their priorities under the conflicting criteria by considering the criteria weights. Chatterjee et al. [35] conducted a comparative analysis on different methods such as AHP, VIKOR, TOPSIS, and COPRAS regarding computational procedures, effortlessness, probability of visual understanding, and type of the data and concluded that the COPRAS Axioms 2022, 11, 154 6 of 24 strategy outperformed the other methods. The COPRAS method allows calculating the values to be maximized and minimized independently among the criteria, and a typical sample distribution is not required [36]. In the results, the determined utility degree expresses how much better the best choice is, in percentage terms, as compared with other MCDM methods. Based on its successful history in decision making research [37,38] and the fact that there are very few studies on 3PL evaluation in the most recent years (from 2018) [39][40][41], we believe COPRAS is a practical and effective technique for selection and evaluation of alternatives.
In many actual settings, however, the human preference model is ambiguous, and experts may be hesitant or unable to assign accurate numerical values to comparison judgments. According to the literature, most studies have used fuzzy sets theory to handle uncertainty in evaluation. However, the totally grey theory has not been developed with MCDM techniques and used in this research area. One advantage of the grey set over the fuzzy sets is the simplified calculation method and the ability to provide more reliable results [42]. Thus, one of our study's purposes is to demonstrate the applicability of greybased multi-criteria methods other than fuzzy MCDM techniques for choosing the best CLPs, in which the grey approach can also eliminate the vagueness of experts' judgments. The uncertain evaluations were expressed as grey numbers to generate a more accurate and robust ranking for the alternatives [43]. The AHP with the grey theory has been successfully applied for specific problems in [44][45][46], and COPRAS with the grey theory presented in [36,41,[47][48][49][50]. However, the combination of the G-AHP and the G-COPRAS for decision making problems has never been reported, at least for 3PL selection.
Devoted to bridging the gap of the existing literature, one of the study's innovations is its evaluation of logistics providers for the cold supply chain in Vietnam, which has never been done before. When evaluating the alternatives in terms of sustainability, a wide range of aspects are taken into consideration after being identified by means of a literature review and experts' opinions, which is a significant advantage of the proposed study. The combination of AHP, COPRAS, and the grey theory, established as relevant and successful approaches for the cold chain industry evaluation, is another innovation. The proposed research provides decision-makers and practitioners information on business performance to aid in the development of sustainability.

Research Framework
This study applied an integrated MCDM method for evaluating cold chain logistics services providers with the case study in Vietnam. The grey system theory is involved in the presented method to handle the uncertainty and vagueness in the decision making process. First, grey AHP is used to estimate the weights of criteria. The weight scores of the evaluation criteria are conducted based on logistics experts' assessments. Then, the grey COPRAS is applied to rank the alternatives through the utility degree percentage. The comparative analysis is conducted to show the effectiveness and applicability of the integrated method. The findings of the applied grey MCDM method provides reliable and robust results for assessing and improving opportunities for decision-makers and stakeholders in the cold supply chain industry. Figure 1 depicts the research framework of this study.

Grey System Theory
The grey system theory was first introduced by Julong-Deng [51] to solve uncertain decision making problems. In the grey theory, the information degree can be separated into three categories which are "black system", "white system", and "grey system" concerning information which is "unknown", "fully known", and "partially known", respectively. Figure 2 shows the framework of the grey system theory. Step 1: Develop the hierarchy tree and build pair-wise comparisons matrix Step 2: Calculate the normalized grey comparison matrix Step 3: Calculate the grey weight of each criterion Step 4: Conduct a whitenization of the grey weight 0

Grey AHP model
Step 1: Determine the relative grey weight of each criterion Step 2: Identify the decision matrix Step 3: Calculate the normalized decision matrix Step 4: Calculate the weighted normalized decision matrix Step 5: Calculate the relative significance of each alternative 0

Grey COPRAS model
Analyze the cold chain logistics service providers evaluation and selection problems Determine the assessment criteria system and alternatives Define the research methods

Establishing expert groups
Literature review

Comparative analysis of methods and result discussions
Step 6: Determine the utility degree of each alternative

Grey System Theory
The grey system theory was first introduced by Julong-Deng [51] to solve uncertain decision making problems. In the grey theory, the information degree can be separated into three categories which are "black system", "white system", and "grey system" concerning information which is "unknown, "fully known", and "partially known", respectively. Figure 2 shows the framework of the grey system theory.
(1) Step 1: Develop the hierarchy tree and build pair-wise comparisons matrix Step 2: Calculate the normalized grey comparison matrix Step 3: Calculate the grey weight of each criterion Step 4: Conduct a whitenization of the grey weight 0 Grey AHP model Step 1: Determine the relative grey w of each criterion Step 2: Identify the decision matri Step 3: Calculate the normalized deci matrix Step 4: Calculate the weighted norma decision matrix Step 5: Calculate the relative significan each alternative 0

Comparative analysis of methods and result discussions
Step 6: Determine the utility degree of alternative

Grey System Theory
The grey system theory was first introduced by Julong-Deng [51] to so decision making problems. In the grey theory, the information degree can into three categories which are "black system", "white system", and "grey cerning information which is "unknown, "fully known", and "partially kno tively. Figure 2 shows the framework of the grey system theory.  A grey number is called by ⊗x = [x, x] where x and x represent the lower bound and the upper bound of the membership function. The grey number's exact value is unknown, but the interval range within which it falls is known. Let ⊗x 1 = [x 1 , x 1 ] and ⊗x 2 = [x 2 , x 2 ] be two grey numbers, k is a positive real number, and L is the length of grey numbers. The basic arithmetic operations are presented in Equations (1)-(6) as follows: (1)

Grey Analytical Hierarchy Process (G-AHP)
The grey analytical hierarchy process (G-AHP) combines the grey theory with the AHP to reduce subjective judgments in decision making. In this paper, the G-AHP is deployed to estimate the grey weights of criteria, in which grey numbers are used, instead of the classical AHP. The G-AHP used in this study is introduced step by step below [52].
Step 1: Identify the research problem, construct the hierarchical tree and pairwise comparison matrix using experts' evaluations with grey numbers. Table 2 shows the linguistic scale with grey numbers used in the G-AHP. Following that, build the integrated grey comparison matrix using the geometrical aggregation, Equation (7) as follows: Extreme importance EMI [8,10] Step 2: Calculate the normalized grey comparison matrix. The normalization for the grey numbers is given in Equations (8)-(10) as follows: where ⊗x ij is the pairwise comparison concerning the ith criterion over the jth criterion.
Step 3: Calculate the grey weight of each criterion by taking averages of the rows using Equation (11) as follows: where n = {1, 2, . . . , N} is the criterion set.
Step 4: Calculate the whitenization of the grey weight to interpret and analyze them easier and better using Equation (12). The whited value of an interval grey weight is a crisp number whose potential value is between the interval grey weight's upper and lower bounds, as follows: where λ is the whitening coefficient and λ ∈ [0, 1].

Grey Complex Proportional Assessment (G-COPRAS)
Zavadskas et al. [53] first introduced the grey complex proportional assessment (G-COPRAS) to reduce the subjective judgments in the decision making problem. The G-COPRAS method prioritizes the alternative based on the calculation of the utility degree. The G-COPRAS method consists of six steps as follows [54]: Step 1 Determine the important weight of each criterion. In this study, the relative importance of each criterion (i.e., grey weights of each criterion) is estimated by using the G-AHP.
Step 2 Determine the relevant criteria, alternatives, and identify the decision matrix. Suppose that A = {A 1 , A 2 , . . . , A m } is a discrete set of m alternatives, which are ranked by a discrete set C = {C 1 , C 2 , . . . , C n } of n criteria.
The performance ratings of the alternatives with respect to criteria is assessed using the linguistic scale with grey numbers in Table 3. Suppose that there are k experts, and the value of alternative h in the criterion g is calculated using Equation (13). Then, the grey decision matrix is built, Equation (14), as follows: where ⊗G hg represents the importance of alternative h in the criterion g. Step 3 Calculate the normalized decision matrix, Equations (15)-(17) as follows: where ⊗G hg is the pairwise comparison concerning the hth alternative in the gth criterion.
Step 4 Then, the weighted normalized grey decision matrix is calculated, Equation (18), as follows: where ⊗w g is the grey weight of each criterion from the G-AHP.
Step 5 Calculate the relative significance of each alternative. First, we calculate the sums P h of the criterion values, i.e., the larger values are more preferable), Equation (19), as follows: Then, we calculate the sums R h of the criterion value, i.e., the smaller values are more preferable, Equation (20), as follows: Following that, the relative significance of each alternative is computed, as can be seen in Equation (21): Step 6 Determine the utility degree of each alternative. First, the optimality criterion K is calculated using Equation (22). Next, the utility degree of each alternative N h is determined by comparing the alternatives under consideration with the best alternative (i.e., 100% for the best alternative), Equation (23), as follows:

A Case Study in Vietnam
In this study, the two-stage grey MCDM method is applied to conduct a case study of six cold chain logistics companies in Vietnam, which are Panasato Cold Storage and Logistics Co.  Table 4. A panel of 15 experts who had been working in the area of cold chains industry for the past ten years was invited to determine potential alternatives and essential criteria for the evaluation to increase the objectivity of the results as much as possible. They were cold supply managers, supply chain managers in different industries, industrial engineers, and academicians. There was a session in which the committee discussed the factors and potential alternatives in the cold chain sector of Vietnam; many vital considerations were referenced and discussed between experts and specialists to determine critical factors for evaluating and selecting possible options. After discussions, the evaluation indicator system was constructed and finalized as the suitable and comprehensive set of criteria responsible for the feasible implementation of the CLPs from a developing country's perspective. Finally, four evaluation dimensions and 15 criteria were shown; the hierarchical tree for the evaluation of cold chain logistics service providers is presented in Figure 3.

Calculation of Grey Weights with the G-AHP
In this section, G-AHP is applied to compute the grey weight of the evaluation criteria for the selection of cold chain logistics service providers. A total of four dimensions are considered, including environmental (C1), social (C2), economic (C3), and service level (C4), which are decomposed into 15 criteria.
The following procedure shows how to calculate the weight of the eigenvector of four dimensions and the process of the consistency ratio calculation. The initial comparison matrix with grey numbers of G-AHP is shown in Table 5 below.

Calculation of Grey Weights with the G-AHP
In this section, G-AHP is applied to compute the grey weight of the evaluation criteria for the selection of cold chain logistics service providers. A total of four dimensions are considered, including environmental (C1), social (C2), economic (C3), and service level (C4), which are decomposed into 15 criteria.
The following procedure shows how to calculate the weight of the eigenvector of four dimensions and the process of the consistency ratio calculation. The initial comparison matrix with grey numbers of G-AHP is shown in Table 5 below. The grey number is transformed to the crisp number to check the consistency ratio (CR) of the performance rating from experts' opinions. The crisp matrix of four dimensions is shown in Table 6. The normalized matrix of G-AHP is, then, obtained by dividing each value in a column of the matrix by its column sum to obtain the priority vector. As shown in Table 7, the priority vector is constructed by averaging the row elements in the normalized matrix. The largest eigenvector (λ max ) is computed to determine the consistency index (CI), the random index (RI), and the consistency ratio (CR), as follows: such that n = 4, we have RI = 0.9, and the consistency ratio (CR) is calculated below.
As shown in CR = 0.0332 < 0.1, the result was satisfactory, and the pairwise comparison matrix was consistent. The same procedure is used to determine 15 criteria. Table A1 shows the integrated grey comparison matrix of the 15 criteria (Appendix A).
As the results from the G-AHP calculation, the grey weights and their transformed crisps are shown in Table 8. From the results, for example, the grey weight of criteria "the green design of cold storages" (C11) has the lowest weight (lower bound) at 0.0563 and the highest weight (upper bound) at 0.0981. Similar to the procedure, the grey weight of criteria "utilization of refrigerant" (C12) has the lowest weight at 0.0367 and the highest weight of 0.0648. The significance levels of 15 criteria of the G-AHP are displayed in Figure 4. The results indicate that the five most significant criteria are C42 (product quality), C31 (logistics costs), C32 (innovation and effectiveness of cold chain processes), C24 (customer experience), and C13 (CO emissions of refrigerated vehicle) with significance levels of 9.53%, 9.18%, 9.18%, 8.52%, and 8.08%, respectively. Meanwhile, C35 (cold chain network management) is specified as the least significant criterion, with a value of 4.10% as compared with other considered criteria. The findings suggest that decision-makers focus on "C42", "C31", "C32", "C24", and "C13" for improving the performance of the cold chain logistics service.

Ranking Alternatives with the G-COPRAS
In this step, the G-COPRAS ranks the alternatives. The preference grey weight of each criterion is obtained from the G-AHP. According to the process of the G-COPRAS, the integrated grey decision matrix of alternative concerning criteria is presented in Table A2 (Appendix A).
Through a numerical example, we demonstrated the calculation process of the G-COPRAS for C1 (environmental) including C11 (the green design of cold storages), C12 (utilization of refrigerant), and C13 (CO emissions of refrigerated vehicle). The calculation for C2 (social), (C3) (economic), and (C4) (service level) is conducted as the same procedures.
The important weight of each criterion is determined. The grey weights of all criteria are obtained from the G-AHP result, as can be seen in Table 9. The integrated grey decision matrix is identified. There are 15 experts used for CLP performance evaluation (k = 15). Table 10 shows the initial evaluation matrix of C1 (environmental). The integrated grey decision matrix of C1 (environmental) is calculated in Table 11 using Equations (13) and (14) in Section 3.4, example for CLP-01 as follows:   Example for CLP-01: The normalized decision matrix is calcualted using Equations (15)- (17) in Section 3.4. The results are presented in Table 12, example for CLP-01 as follows: Example for CLP-01:  Then, the weighted normalized grey decision matrix is calculated using Equation (18) in Section 3.4, as can be seen in Table 13, example for CLP-01 as follows: Example for CLP-01: Next the relative significance of each alternative is calculated according to benefit criteria and cost criteria using Equations (19)- (21) in Section 3.4. Finally, the utility degree of each alternative is determined using Equations (22) and (23) in Section 3.4.
The evaluation of the utility degree of the G-COPRAS is shown in Table 14. From the result, Yoshida Saigon Cold Logistic Co., Ltd. (CLP-05) achieves the highest performance with a utility degree of 100%. ABA Cooltrans Business Solutions Corporation (CLP-06) ranks second, with a utility degree of 94.96%, and Panasato Cold Storage and Logistics Co., Ltd. (CLP-01) ranks third with a utility degree of 93.69%. Meanwhile, Tan Bao An Logistics Joint Stock Company (CLP-03) has the lowest performance with a utility degree of 66.38%. The ranking performance is visualized in Figure 5.

Comparative Analysis
In this section, a comparison among MCDM methods is performed to illustrate th applicability and rationale of the provided G-AHP and G-COPRAS methodologies. COP RAS [55], G-EDAS [56], and G-TOPSIS [57] methodologies are used to rank the alterna tives in this work. Table 15 displays the ranking results of the compared methods, whic are also depicted in Figure 6. According to the comparative findings, the rating of the col chain logistics providers with a case study in Vietnam produced the same results as th methods provided in this work, which is the integration of the G-AHP and G-TOPSIS, G AHP and G-EDAS, and G-AHP and G-COPRAS methodologies. The rankings of the G AHP and COPRAS methods deviate somewhat from the suggested method. There is slight distinction between Freight Mark Vietnam Co., Ltd. (CLP-04), Yoshida Saigon Col

Comparative Analysis
In this section, a comparison among MCDM methods is performed to illustrate the applicability and rationale of the provided G-AHP and G-COPRAS methodologies. CO-PRAS [55], G-EDAS [56], and G-TOPSIS [57] methodologies are used to rank the alternatives in this work. Table 15 displays the ranking results of the compared methods, which are also depicted in Figure 6. According to the comparative findings, the rating of the cold chain logistics providers with a case study in Vietnam produced the same results as the methods provided in this work, which is the integration of the G-AHP and G-TOPSIS, G-AHP and G-EDAS, and G-AHP and G-COPRAS methodologies. The rankings of the G-AHP and COPRAS methods deviate somewhat from the suggested method. There is a slight distinction between Freight Mark Vietnam Co., Ltd. (CLP-04), Yoshida Saigon Cold Logistic Co., Ltd. (CLP-05), and ABA Cooltrans Business Solutions Corporation (CLP-06). These results show the effectiveness of the applied techniques for evaluating and selecting the cold chain supplier selection with a case study in Vietnam. As a consequence of the above comparison research, it is reasonable to conclude that the ranking orders of the considered alternatives are credible and may serve as a valuable reference for users or decision-makers in picking the appropriate cold chain suppliers, particularly during the COVID-19 pandemic.   Figure 6. Comparative analysis of methods.

Research Implications
The research findings were discussed with experts to analyze the issues related to sustainability-focused CLP initiatives and there was good agreement with the results. In this sense, this study has several beneficial implications for managers who both outsource and offer logistics activities, stakeholders, as well as academicians, as follows: For managers who intend to outsource their cold logistics activities, this study reveals some challenges in evaluating and selecting the right logistics providers for sustainability in the cold supply chain context, which has managerial implications for managers to support their businesses while remaining competitive in the market. Especially for a developing country, sustainable 3PL decisions have become necessary as organizations are aware of society and environment protection and increasingly compete on environmental, economic, and social supply chain capabilities. Important criteria and dimensions of criteria related to CLP selection have been recognized. The proper understanding of these criteria and dimensions can help in managing the challenges or barriers in developing

Research Implications
The research findings were discussed with experts to analyze the issues related to sustainability-focused CLP initiatives and there was good agreement with the results. In this sense, this study has several beneficial implications for managers who both outsource and offer logistics activities, stakeholders, as well as academicians, as follows: For managers who intend to outsource their cold logistics activities, this study reveals some challenges in evaluating and selecting the right logistics providers for sustainability in the cold supply chain context, which has managerial implications for managers to support their businesses while remaining competitive in the market. Especially for a developing country, sustainable 3PL decisions have become necessary as organizations are aware of society and environment protection and increasingly compete on environmental, economic, and social supply chain capabilities. Important criteria and dimensions of criteria related to CLP selection have been recognized. The proper understanding of these criteria and dimensions can help in managing the challenges or barriers in developing sustainabilityfocused CLP selection criteria and evaluation decisions. All criteria considered in this research are helpful in addressing various challenges among organizations in developing sustainability aspects in 3PL-related decisions.
From the academic point of view, this study presents an integrated approach based on the G-AHP and G-COPRAS methods for CLP selection under limited or no quantitative information. Expert committees and grey theory are particularly suitable to address the challenge of assessing sustainability practices while minimizing the effects of imprecise or missing data. The G-AHP successfully provided consistent criteria ratings, whereas the G-COPRAS generates alternative rankings based on utility degree in the form of percentages. The comparative analysis will allow practitioners to test the observation stability.
Further, from the presented case study, this study provides a scientific means for managers who provide their logistics services in the cold supply chain to determine their strengths and weaknesses and to also improve their sustainable processes. The important benefit that this study offers is the development of evaluation criteria using literature and feedback from experts. For the whole society in the supply chain context, applying a comprehensive selection method (as presented in this paper based on G-AHP and G-COPRAS methods) can make the supply chain more transparent for internal and external stakeholders, therefore, helping businesses gain more accountability and decrease their vulnerability.

Conclusions and Future Studies
Outsourcing cold chain logistics services is a growing trend that brings businesses a number of undeniable benefits. Specifically, CLPs assist businesses in the cold supply chain in Vietnam and many other countries in thriving and ensuring a more efficient goods distribution process and a complete supply chain that can meet the criteria of sustainable development (service and product quality, technology, environmentally friendly innovations, etc.) to satisfy customers' requirements. Thus, choosing a suitable CLP is a vital strategic decision. In this paper, a grey-based MCDM method for CLP selection in a developing country is presented. The suggested method allows the G-AHP to determine the weights of the assessment criteria and subsequently the G-COPRAS to rank the alternatives. A case study in Vietnam is used to demonstrate the reliability of the suggested integrated framework. Fifteen evaluation criteria were chosen from the literature after consulting with industry experts. Product quality, logistics costs, innovation and effectiveness of cold chain processes, customer experience, and CO emissions of refrigerated vehicle are the five evaluation criteria that obtained the highest weight priority in the analysis. The final ranking indicates that Yoshida Saigon Cold Logistic company (CPL-05) is the best CLP among alternatives. A comparative analysis was conducted to test the robustness of the approach, with the results illustrating that the applied methods reach common CLP rankings. This shows that the applied approach is practical in nature.
The contributions of the present study are three-fold. First, our study aims to evaluate cold supply chain logistics providers in the context of Vietnam, which has never been investigated before in the literature. While evaluating the alternatives, a set of sustainability criteria is examined, which is a significant advantage of this study. Second, in terms of method, the combination of AHP, COPRAS, and the grey theory is proposed to solve the problem that has been identified as appropriate and effective methodologies for the cold chain industry evaluation. The application of grey-based MCDM methods for the evaluation of 3PL in general and selecting CLPs, in particular, is limited. Thus, this study shows the applicability of other grey-based multi-criteria methods for the logistics providers evaluation in an integrated manner. The presented approach was applied for Vietnam to show its effectiveness and the applicability. On the one hand, compared to the classical AHP, the introduced G-AHP method can solve intricate problems in uncertain conditions, which makes it a very effective technique and provides more information on specific elements of experts' preferences. In addition, very few studies have used the integration of the G-COPRAS tool for 3PL evaluation. The grey relationship in G-COPRAS makes it more practical when dealing with uncertain factors and the method creates results with higher accuracy, fewer computation steps, and a shorter calculation time. Third, for managerial implications, the assessments presented in this research offer valuable materials for managers, practitioners, and stakeholders of cold chain businesses not only to assess the most efficient CLP in the supply chain in the market but also to devise suitable strategies toward sustainable development.
Although the methodology adopted in this study has been conducted successfully in prioritizing different alternatives and factors, it was not without some limitations. One limitation can be the use of the AHP method. Although the consistency check in the present study has been fulfilled, it is inconceivable to neglect the inconsistency in the pairwise comparison matrix that might occur in practice for other problems. The best worst method (BWM) can overcome this drawback as it unburdens decision-makers by requiring fewer pairwise comparisons than the conventional AHP procedure [58] or the multi-level parsimonious analytic hierarchy process (PAHP) method, which can significantly reduce survey [59]. The analytic network process (ANP) method can also be a better option to avoid the interrelationship of factors. Hence, these methods are recommended for future studies. Another limitation is that the evaluation process of CLPs depends on the involvement of experts; thus, the results are based on personal opinions, knowledge, and judgment. To avoid this limitation, 15 experts were utilized to provide different preferences. Thus, different multi-criteria evaluation techniques such as VIKOR, TOPSIS, WASPAS, PROMETHEE, and MACBETH (measuring attractiveness by a categorical based evaluation technique) could be employed to achieve the same goal, and findings could be compared.
The methodology used in this study can also be extended within the dynamic and uncertain environment in future research by integrating novel criteria factors responsible for coordination in the supply chain, especially those regarding the current crisis. In addition, the present study could be applied to specific supply chain scenarios in different industries and countries with other multi-criteria methods as outlined earlier to see if the findings are generalizable.