Evaluating Global Container Shipping Companies: A Novel Approach to Investigating Both Qualitative and Quantitative Criteria for Sustainable Development

: The COVID-19 pandemic has implications for the container shipping industry and global supply chains. Measuring the efﬁciency of major international container shipping companies (CSCs) is an important issue that helps them make strategic decisions to improve performance, especially in the context that all businesses and governments are adapting to build back better the post-pandemic world. This paper develops a new integrated approach using both a qualitative assessment tool and a performance assessment tool as a systematic and ﬂexible framework for evaluating the container shipping industry. This new methodology is implemented in two phases to consider both qualitative and quantitative criteria for assessing the performance of CSCs based on efﬁciency. In the ﬁrst phase, qualitative performance evaluation is performed using spherical fuzzy analytical hierarchical process (AHP-SF) to ﬁnd criteria weights and then the grey complex proportional assessment methodology (COPRAS-G) is used to ﬁnd the ranking of CSCs. Qualitative variables are converted into a quantitative variable for use in the data envelopment analysis (DEA) model as an output called an output variable called expert-based qualitative performance (EQP). Then, DEA is performed to identify efﬁcient and inefﬁcient CSCs with the EQP variable and other quantitative parameters (i.e., capacity, lifting, expenses, revenue, and CO 2 emissions). The efﬁciency of 14 major global CSCs is empirically evaluated, and the scores for CSCs’ efﬁciency in all dimensions are measured and examined. The results show that the average cargo efﬁciency of the CSCs is lower than their eco-efﬁciency performance, revealing the operational disruption caused by the pandemic. Moreover, by identifying efﬁcient and inefﬁcient CSCs, our ﬁndings provide practical implications for decision-makers in the maritime ﬁeld and assist in modifying applicable policies and strategies to achieve sustainable performance.


Research Background
Globalization has enhanced the importance of maritime transportation in recent decades as global trade has grown [1]. The container shipping business is an important aspect of the worldwide supply chain's flow of commodities. It is the most efficient and cost-effective method of long-distance transportation of large goods. Thanks to efficient

CSCs' Efficiency Measurements
Running a CSC is difficult due to the lack of significant distinctions in the essential services supplied by CSCs. Shippers may be enticed to move their cargo to another CSC due to a minor difference in freight rates [4]. Global economic volatility, problems in relocating empty containers, and extra capacity brought on by enlarged ships all add to the complexity and challenges of running a container shipping business. Despite these enormous expenses and hurdles, CSCs can survive and improve the operational efficiency of their shipping service production by leveraging their experience and tactics. A production system's function can be defined as transforming input into output [6]. This transformation can be thought of as a process of change-the operand results in extra values to meet the requirements of such changes. Shipping services are a function of both capital and labor [7]. In operational efficiency analysis, CSCs' performance is defined as a process that transforms inputs such as fleet capacity, workforce, and expenses to containers carried. Then, the volume of cargo carried (or so-called lifting) becomes a source of earnings (revenue) as a good output. The efficiency measures for CSCs' performance using different input-output combinations can be different. If lifting is used as an output, named cargo efficiency, this only reflects the relative ability of a company in carrying cargoes. Meanwhile, measuring efficiency using financial indicators of a CSC as outputs, called economic efficiency (or financial efficiency), can give insights into the financial performance of each company. Along with revenue gains, a CSC also causes harmful emissions from operation to the environment and society, which is seen as a bad output in the evaluation system. This leads to the need to consider the ecological efficiency of CSCs.
Due to its operating consequences on the maritime and air environments, the shipping sector has been at the center of attention on a worldwide scale. The impact of shipping activities on environmental pollution is highlighted in various discussions [8][9][10][11][12][13][14]. The global crisis has once again reminded the public of environmental sustainability trends. These trends have been steadily developing over the last decade, but they have accelerated during the epidemic and are continuing to revolutionize maritime transportation and trade. Maritime transportation is under increasing pressure to decarbonize and operate more sustainably; additional challenges have also surfaced in the aftermath of the pandemic. According to the International Maritime Organization (IMO), the share of shipping emissions in global anthropogenic emissions reached nearly 3% in 2018 [15]. Despite maritime transportation being the safest and most energy-efficient mode of transportation compared to other modes, especially when considering transported volumes, the industry's contribution to global CO 2 emissions is expected to increase by 50-250% due to increased international trade until 2050. Global shipping must thus strike the right balance between its role as a private company and its social obligations to protect the environment. CSCs are concentrating on boosting eco-efficiency due to growing worries about climate change and growing knowledge of the detrimental effects of hazardous compounds on human health. The eco-efficiency idea, which is essentially the performance of the economics and the environment combined, has been examined in a number of recent studies and is regarded as the top research for future trends, not least for the shipping sector. To incorporate negative environmental impacts into CSCs' efficiency analysis, the eco-efficiency model investigates how operating costs, workforce, and liftings transform into one desirable output (revenue) and undesirable output (CO 2 emissions) [16].

Objectives of Present Study
As for each CSC, its relative performance compared to other firms is measured using several methods. For studying the overall performance of the CSCs, we view the global CSCs as production units and calculate the CSCs' performance in translating their resources into specific forms of outputs. As previously outlined, the efficiency of CSCs concerns the utilization of resources in multiple dimensions and generates different types of outcomes such as lifting (cargo-efficiency model), revenue, and CO 2 emissions (eco-efficiency model). Therefore, assessment should take multiple inputs and outputs into account simultaneously. With this being the case, data envelopment analysis (DEA) is a useful tool for estimating the production frontier and determining the efficiency of different organizations. It has been widely used in studies of relative port and shipping efficiency and, more recently, for efficiency analysis while considering the negative effects of emissions. In existing DEA-related articles, cargo efficiency and economic efficiency of specific CSCs have been investigated, including in studies such as [4,7,[17][18][19][20]. Eco-efficiency has been discussed in the most recent articles [2,16,21].
In addition, the determination of the performance efficiency of a CSC is complicated further by the fact that while considering inputs and outputs that are quantitative (fleet capacity, expenses, workforce, revenue, the amount of CO 2 emissions, etc.) by the DEA models, it is vital to explore the qualitative criteria influencing a shipping enterprise's competitive advantage toward sustainable development. Numerous qualitative factors have recently become increasingly complex as strategic operations, environmental, social, political, counterparty, and customer satisfaction concerns have been discussed in various studies on CSC performance evaluation [22][23][24][25][26][27][28]. In light of this, the multi-criteria decision making (MCDM) techniques are used through experts' judgments or questionnaire surveys. As a result, a method is required to optimize the decision maker's attitude toward each criterion's importance while incorporating multiple aspects [29]. This study develops an integrated AHP-SF (spherical fuzzy analytical hierarchical process), COPRAS-G (grey complex proportional assessment), and DEA methodologies evaluation framework for global CSCs' performances, considering qualitative and quantitative factors. The inputs, outputs, and criteria have been determined by means of reviewing the literature review and experts' discussions. Qualitative assessment results could be integrated as a novel parameter for the quantitative analysis using the DEA model. This framework is comprehensive because, while decision-makers of CSCs can consider expert judgements of qualitative performance, the final efficiency scores from the DEA model suggest that inefficient CSCs make strategic decisions to improve their performance toward sustainable development.
The remaining sections are arranged as follows: Section 3 elaborates on the AHP-SF, COPRAS-G, and slack-based measure DEA (SBM-DEA) techniques. Section 2 gives a review of the literature. The empirical analysis is presented in Section 4. Section 5 concludes by outlining the discussions, limits, and possible directions for further research.

Literature Review on Efficiency Analysis of CSCs
In many industries, for relating efficiency outcomes to features of decision-making organizations, the application of DEA is regarded as an appropriate methodology. DEA was first introduced by Charnes et al. [30,31] to assess the efficiency and productivity of a decision-making unit (DMU) through multiple inputs-outputs. For efficiency analysis in the shipping industry based on quantitative dimensions, DEA models are used so that CSCs can understand their strengths and weaknesses and increase their performance. Two DEA models that are widely used are the Charnes, Cooper, and Rhodes (CCR), and the Banker, Charnes, and Cooper (BCC). Lun and Marlow [17] used two inputs (shipping capacity and operating costs) and two outputs to calculate the efficiency of major global CSCs (profit and revenue). With three inputs (total assets, workforce, and capital expenditures) and sales as an output, Panayides et al. [18] conducted an efficiency analysis of important international shipping corporations. Using a two-stage DEA technique, Bang et al. [7] investigated the operational and financial efficiencies of 14 liner shipping enterprises. In order to assess the effectiveness of international container shipping lines, Gutiérrez et al. [19] developed a bootstrap DEA approach: labor force, the number of ships, and fleet capacity were used as input variables, and containers carried and annual turnover were used as output variables, respectively. Chao et al. [4] used a dynamic network DEA with shared inputs to assess the cargo efficiency and economic efficiency of 13 major global CSCs. Chen et al. [2] considered the environmental efficiency of Taiwanese shipping firms using the network centralized DEA model for resource allocation. The correlations between economic and cargo efficiencies, environmentally adjusted efficiencies, and environmental efficiencies were determined by Gong et al. [16] using the DEA-SBM (slack-based model) for various input-output combinations. Their proposed model was applied for international shipping companies as their business units. Kuo et al. [21] developed two separate DEA models to obtain the cargo and eco-efficiency scores of 10 global CSCs and explore the determinants of the CSCs' efficiencies.
The existing literature reveals that most CSCs' efficiency evaluation approaches have filtered out efficient CSCs on measurable factors (capacity, expenses, workforce, emissions, etc.). In some cases, conducting only quantitative analysis in assessing organizational performance in the shipping industry is insufficient due to the need to investigate a CSC's efficiency from other aspects that are intangible and unmeasurable, such as their impacts on society and their counterparty, or strategy-related efficiency. In this case, experts' involvement is a practical method to provide qualitative judgments for CSCs based on their knowledge. To prioritize the competitive advantages of container shipping liner companies, Bao et al. [28] developed a multi-attribute decision making (MADM) model with intuitionistic fuzzy linguistic variables to evaluate both quantitative and qualitative data with experts' judgments. Each business's efficiency of technological upgrading, workplace safety, gender quality, personnel training, and other quantitative criteria are investigated in this study. The AHP method is most recognized in handling qualitative and subjective measurements of experts through pairwise comparison [32]. Yoon et al. [24] used the AHP method with fuzzy logic for uncertain information to evaluate the performance of CSCs in Vietnam under five main criteria: service, operation, service cost, counterparty, and financial status. The decision-making trial and evaluation laboratory (DEMATEL) method also uses pairwise comparison for collaborative decision-making to evaluate the relative importance of criteria. Based on the experience and knowledge of appropriate shipping industry experts obtained through an effective questionnaire, Hsu and Ho [27] determined key factors toward the performance of CSCs in Taiwan from the perspective of domestic high-tech industry shippers. For evaluating the suitability and relevance of customer demands, customer expenses, customer communication, and customer convenience, the  Table 1 summarizes studies on CSCs' efficiency assessment.  Table 1 summarizes the methods used for global shipping evaluation in papers published from the years 2011-2021. It can be observed that DEA is a common method used for quantitative assessments. Traditional DEA models are CCR (Charnes, Cooper, and Rhodes) [30] and BCC (Banker, Charnes, Cooper) [35] employ radial and orientated techniques. However, the relative efficiency of the radial DEA efficiency measurement will be exaggerated when the input or output has nonzero slacks, and the efficiency result of the orientated DEA measuring technique is equally unreliable. Thus, the slacks-based measure model (SBM) was proposed by Tone [36] to overcome this problem. The SBM is a non-radial DEA model, as opposed to the radial CCR and BCC models. All inputs and outputs change proportionally as a result of the improvement of invalid DMU in the radial DEA model. The main advantage of the SBM model is that it improves inaccurate DMU in a way that directly addresses input or output slacks in addition to equal proportion improvement [37]. Additionally, traditional DEA models usually aim for the highest possible output, but when evaluating an economic system's production or operation efficiency, unexpected outputs are always present alongside expected ones. It is challenging to prevent negative effects, such the release of contaminants. The approach for evaluating efficiency should take into consideration any unfavorable externalities brought on by the production process. As a consequence, Tone suggested expanding the SBM model to examine the connection between production and pollution and to take unexpected output into account. SBM has gained Axioms 2022, 11, 610 6 of 31 popularity as a method of addressing environmentally unfavorable results in efficiency evaluation [38][39][40][41].
However, it is evident that the AHP-SF and COPRAS-G methods have almost been missed in the literature of global shipping evaluation, although they are among the latest MCDM methods. The novel spherical fuzzy sets theory (SFS) combined with AHP for calculating the spherical criteria weights is proposed, while the COPRAS technique under grey theory is utilized for alternative ranking. These methods are among the latest MCDM methods, yet they have been almost missing in the literature of global shipping evaluation. Spherical fuzzy set (SFS) developed in 2019 [42] is a three-dimensional fuzzy set that was created as an extension of the intuitionistic fuzzy set, PFS, and neutrosophic logics, specifically to handle uncertainty during the quantification of expert judgments. SFS improves decision-making over its predecessors by making it more intelligent (comparable to human judgment), which can lead to a decision-making process that is more accurate when evaluating options. Because of this benefit, SFS has lately been used in a variety of applications, including location selection [43,44], renewable energy evaluation [45], supplier selection [46][47][48][49], systems management [50][51][52][53][54][55][56], technology evaluation [57], and strategies selection [58].
The COPRAS-G method ranks alternatives based on their significance and utility degree using a stepwise evaluation procedure [59,60]. The parameters of the alternatives are determined using the grey relational grade and expressed in terms of intervals in this method. The grey systems theory is concerned with the investigation of problems involving small samples and limited information [61]. It is concerned with uncertain systems with only a portion of the information known. Situations in which there is no information are categorized as black, whereas those in which all the information is excellent are categorized as white. On the other hand, these hypothetical situations are rarely present in real-world problems. The areas between these two extremes are described as being grey, ambiguous, or fuzzy. As a result, a grey system is one in which some information is certain and some is unknown. One is constantly in the middle, between the back and white extremes, as there is always some degree of uncertainty (i.e., somewhere in the grey area). Recent applications of COPRAS-G are integrated with other MCDM methods as a more systematic and comprehensive framework. Aghdaie et al. [62] implemented market segment evaluation and selection based on application of fuzzy AHP and COPRAS-G methods. Tavana et al. [63] proposed a novel hybrid social media platform selection model using fuzzy ANP and COPRAS-G. Caisucar et al. [64] used fuzzy-TOPSIS and COPRAS-G approach for validation of portfolio allocation in New Product Development. Kayapinar Kaya and Aycin [65] presented an integrated interval type 2 fuzzy AHP and COPRAS-G methodologies for supplier selection in the era of Industry 4.0.

Research Gaps
It is deduced that many previous approaches have suffered from certain shortcomings, such as only considering quantitative parameters or being too subjective in evaluating the shipping businesses. While the problem of evaluating CSCs has been investigated by many scholars, mostly under quantitative variables in the container shipping context (fleet capacity, employees, expenses, revenue, greenhouse gas emissions, etc.), the assessment of the competitive advantage benefits of CSCs is meager. In this direction, evaluation dimensions are often difficult-to-quantify qualitative criteria such as counterparty, strategies, and social and environmental aspects, to name a few. Therefore, it is also essential and instructive to evaluate CSCs' competitive advantages and this can be done using multi-attribute analysis employing expert panels. An integrated approach using both a qualitative assessment tool and a performance assessment tool can be developed as a systematic and flexible framework for evaluating the container shipping industry.
To fulfill the research gap, we used three MCDM methods to find efficient and inefficient international CSCs sensitively. In the first step, the CSCs were evaluated with expert judgments on the sustainability criteria of strategic operation, service level, social and environmental aspects, and counterparty. The spherical fuzzy analytical hierarchical process (AHP-SF) was used to determine the criteria weights. Then, the grey complex proportional assessment (COPRAS-G) method was utilized to find the ranking of CSC. Qualitative variables were transformed into a quantitative variable for use in the DEA model, as an output variable called expert-based qualitative performance (EQP). Then, DEA was performed to identify efficient and inefficient CSCs with the EQP variable and other quantitative parameters (i.e., capacity, lifting, expenses, revenue, emissions, etc.). All steps are explained step by step in Figure 1. and instructive to evaluate CSCs' competitive advantages and this can be done using multi-attribute analysis employing expert panels. An integrated approach using both a qualitative assessment tool and a performance assessment tool can be developed as a systematic and flexible framework for evaluating the container shipping industry.
To fulfill the research gap, we used three MCDM methods to find efficient and inefficient international CSCs sensitively. In the first step, the CSCs were evaluated with expert judgments on the sustainability criteria of strategic operation, service level, social and environmental aspects, and counterparty. The spherical fuzzy analytical hierarchical process (AHP-SF) was used to determine the criteria weights. Then, the grey complex proportional assessment (COPRAS-G) method was utilized to find the ranking of CSC. Qualitative variables were transformed into a quantitative variable for use in the DEA model, as an output variable called expert-based qualitative performance (EQP). Then, DEA was performed to identify efficient and inefficient CSCs with the EQP variable and other quantitative parameters (i.e., capacity, lifting, expenses, revenue, emissions, etc.). All steps are explained step by step in Figure 1. The paper contributes to the existing literature with different concepts of efficiency measures within the shipping industry. To the best of the authors' knowledge, this is the first attempt to utilize the merits of AHP-SF, COPRAS-G, and SBM-DEA for global shipping evaluation. The performance of CSCs is evaluated under both quantitative and qualitative dimensions in terms of sustainability, which is a significant advantage of this study. There is a need for a more systematic and analytical framework for global shipping evaluation; this study is the first integration of AHP-SF, COPRAS-G, and DEA, with the The paper contributes to the existing literature with different concepts of efficiency measures within the shipping industry. To the best of the authors' knowledge, this is the first attempt to utilize the merits of AHP-SF, COPRAS-G, and SBM-DEA for global shipping evaluation. The performance of CSCs is evaluated under both quantitative and qualitative dimensions in terms of sustainability, which is a significant advantage of this study. There is a need for a more systematic and analytical framework for global shipping evaluation; this study is the first integration of AHP-SF, COPRAS-G, and DEA, with the qualitative variables being transformed into a new quantitative variable for the typical CSCs' efficiency model. Thus, the proposed approach and implications are significant materials for decision-makers in the ocean transportation industry to run their CSCs in a sustainable manner, especially in reducing emissions from shipping, along with devising strategies for the industry stakeholders.

Spherical Fuzzy Analytical Hierarchy Process (AHP-SF)
Spherical fuzzy set (SFS) was created by Kutlu Gündoudu and Kahraman [45] as the most current emergence of the fuzzy sets, which can better manage uncertainties and ambiguities in the decision-making process. As shown in Figure 2, each spherical fuzzy number contains the membership, non-membership, and hesitancy functions associated with the interval [0, 1]. qualitative variables being transformed into a new quantitative variable for the typical CSCs' efficiency model. Thus, the proposed approach and implications are significant materials for decision-makers in the ocean transportation industry to run their CSCs in a sustainable manner, especially in reducing emissions from shipping, along with devising strategies for the industry stakeholders.

Spherical Fuzzy Analytical Hierarchy Process (AHP-SF)
Spherical fuzzy set (SFS) was created by Kutlu Gündoudu and Kahraman [45] as the most current emergence of the fuzzy sets, which can better manage uncertainties and ambiguities in the decision-making process. As shown in Figure 2, each spherical fuzzy number contains the membership, non-membership, and hesitancy functions associated with the interval [0,1].
Comparison of IFS, PFS, NS, and SFS in 3D space.

Definition 1:
Singer value SFS F S of the universe of discourse X is presented by Equations (1)- (3).
and hesitancy levels of x to F S , respectively.
Compared to the traditional AHP model, the AHP-SF model offers numerous benefits. While the method collects data from experts, it may not accurately reflect the expressed opinions. As a result, AHP-SF may readily eliminate the uncertainty caused by expert opinion in the comparison matrix. The AHP-SF model was used to calculate the weights of the criteria in this paper. The AHP-SF process has six steps, as follows [57].
Step 1: Define the hierarchical structure of the problem.
The hierarchical structure is organized with the research goal (level 1) and the list of criteria C = {C 1 , C 2 , . . . C n } (level 2) within n ≥ 2.
Step 2: Construct pairwise comparison matrices. As demonstrated in Table 2, the pairwise comparison matrices are built in consideration of spherical fuzzy linguistic scales. The score indices (SI) are determined by Equations (17) and (18).
for the AMI, VHI, HI, SMI, and EI.

Linguistics Terms Spherical Fuzzy Number Score Index
Absolutely more importance (AMI) (0.9, 0.1, 0.0) 9 Very high importance (VHI) Step 3: Check consistency. The corresponding SI is transformed from the linguistics scales. The consistency ratio (CR), which must be less than 10%, is next evaluated for the pairwise comparison matrices.
Step 4: Determine the spherical fuzzy local weights of the criteria. Using the SWAM operator to determine the weight of each condition by Equation (19).
Step 5: Aggregate spherical fuzzy weights. Equation (20) is used to defuzzify the criterion weights. They are then normalized using Equation (21). The final ranking scores are aggregated using the multiplication operator in Equation (22).
Using Equation (23), spherical fuzzy arithmetic addition over global weights is used to determine the final AHP-SF score ( F).
Step 6: Defuzify the final score of each criterion. Sort the criteria list by their defuzzified final rating: the greater the number, the better. In the following step, the criterion weights are employed for the COPRAS-G model.

Grey Complex Proportional Assessment (COPRAS-G)
Julong [66] created the grey theory to investigate uncertainty with ambiguous information. The grey hypothesis is separated into three categories, referred to as the "white system", "black system", and "grey system" depending on how much of the knowledge is "completely known", "unknown", and "partially known", respectively. The grey theory notion is shown in Figure 3. Using Equation (23), spherical fuzzy arithmetic addition over global weights is used to determine the final AHP-SF score (̃).
Step 6: Defuzify the final score of each criterion. Sort the criteria list by their defuzzified final rating: the greater the number, the better. In the following step, the criterion weights are employed for the COPRAS-G model.

Grey Complex Proportional Assessment (COPRAS-G)
Julong [66] created the grey theory to investigate uncertainty with ambiguous information. The grey hypothesis is separated into three categories, referred to as the "white system," "black system," and "grey system" depending on how much of the knowledge is "completely known," "unknown," and "partially known," respectively. The grey theory notion is shown in Figure 3. Let ⨂ = [ , ] represent a grey number with denoting the lower limit and denoting the upper limit of the membership function.
The concept of grey complex proportional assessment (COPRAS-G), which aims to lessen subjective judgements in the decision-making process, was initially created by Zavadskas et al. [59]. There are many benefits to using COPRAS-G in the decision-making process, including: (1) fewer calculation steps than TOPSIS or WASPAS methods; (2) the ability to calculate the values to be maximized and minimized separately among the criteria; (3) usually, distribution of samples is not necessary; and (4) compared to other MCDM methods, the estimated utility degree of COPRAS-G reveals how much better the optimal alternative is than the other in percentage terms [65]. Let ⊗x = [x, x] represent a grey number with x denoting the lower limit and x denoting the upper limit of the membership function.
Let ⊗x 1 = [x 1 , x 1 ] and ⊗x 2 = [x 2 , x 2 ] be two grey numbers; ε denotes a positive real number and L denotes the length of grey number. The basic grey number arithmetic operations are shown in Equations (24)- (29).
The concept of grey complex proportional assessment (COPRAS-G), which aims to lessen subjective judgements in the decision-making process, was initially created by Zavadskas et al. [59]. There are many benefits to using COPRAS-G in the decision-making process, including: (1) fewer calculation steps than TOPSIS or WASPAS methods; (2) the ability to calculate the values to be maximized and minimized separately among the criteria; (3) usually, distribution of samples is not necessary; and (4) compared to other MCDM methods, the estimated utility degree of COPRAS-G reveals how much better the optimal alternative is than the other in percentage terms [65].
All of the evaluation criteria in this work are subjective, and it is impossible to precisely describe their values. Grey numbers are therefore a useful way for the COPRAS-G technique to communicate the actual situation during the decision-making process. According to the utility degree computation, the COPRAS-G technique prioritizes the alternative. There are six steps in the COPRAS-G method [68].
Step 1: Identify the relevant criteria and alternatives. 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.
Step 2: Construct the decision matrix.
Utilize the linguistic scale with grey values in Table 3 to evaluate how well the alternatives performed in relation to the criteria. Suppose that there are k experts, and the value of alternative h in the criterion g is calculated using Equation (30). Following that, the grey decision matrix is built, as can be seen in Equation (31).
where ⊗G hg is the importance of alternative h in the criterion g.  [6,9] Very Good (VG) [9,10] Step 3: Determine the important weight of each criterion. The AHP-SF method is employed to calculate the significance level of criteria.
Step 4: Determine the weighted normalized decision matrix. First, use Equations (32)- (34) to develop the normalized grey decision matrix.
where ⊗G hg represents the pairwise comparison from a group of decision-makers with respect to the h th alternative in the g th criterion.
Weighted normalized grey decision matrix is then created, as shown in Equation (35).
where w g is the important weight of each criterion.
Step 5: Determine the relative significance of each alternative. First, we compute the sums P h of the criterion values (the larger values are better) using Equation (36).
Next, we compute the sums R h of the criterion value (the smaller values are better) using Equation (37).
After that, Equation (38) is used to calculate the relative relevance of each alternative.
Step 6: Calculate the utility degree of each alternative. First, Equation (39) is used to derive the optimality criterion K. The utility degree of each alternative N h is then determined by contrasting them with the best alternative (which has a utility degree of 100%), as shown in Equation (40).

Data Envelopment Analysis (DEA)
Data envelopment analysis (DEA), which does not rely on the premise that the data are normal, is an effective non-parametric methodology. In the DEA model, homogeneity and isotonicity are two essential premises. Before using the DEA model, the correlation between the inputs and outputs must be confirmed, and it must be a complete positive linear correlation. The correlation of Pearson's (r) of two variables (x) and (y) is calculated, as can be seen in Equation (41) [69].
where n is the sample size; x i and y i denote the individual sample points indexed with i; and x = 1 n ∑ n i=1 x i is the mean of the sample which is analogous for y. Slacks-based measure (SBM) is one of the many widely utilized DEA models. It can be used for efficiency analysis either without taking into account the negative environmental impact [36] or by adjusting for the negative outputs, a so-called DEA undesirable output model [70]. This paper analyzes the economic efficiency and environmental efficiency of the shipping companies at the same time. Hence, the adjusted DEA model ensures consistency and enables comparability. The list of symbols and notations used in the adjusted DEA model is presented as follows [16].
• n: number of decision-making units (DMUs), as shipping companies in this paper • DMU j : the j th DMU, j = 1, 2, . . . , n, has m inputs, s 1 desirable outputs, s 2 undesirable outputs • x ij (i = 1, 2, . . . , m): the i th input of the j th DMU, the matrix as X ∈ R m×n • y grj (r = 1, 2, . . . , s 1 ): the r th desirable output, the matrix as y g ∈ R s 1 ×n • y bqj (q = 1, 2, . . . , s 2 ): the q th undesirable output, the matrix as y b ∈ R s 2 ×n • δ = 1: the model adjusts the desirable output by the negative impact of shipping emissions (i.e., undesirable output model) • δ = 0: otherwise The adjusted DEA model with or without considering the negative impact of shipping emissions can be written as: . . , s b s 2 o ) denote slack variables of inputs, desirable outputs, and undesirable outputs, respectively. λ = (λ 1 , . . . , λ n ) denotes a positive weight vector. The subscript o represents the DMU evaluated.
Through the Charnes-Cooper transition [70], Equation (42) can be transformed into Equation (43) with an integer t (linear programming model) as follows.
The adjusted DEA model can be applied for situations either without considering the environmental impact or by adjusting for the bad outputs and can be solved as a linear programming model. The optimal value of the objective function is at 1, or the efficiency score is 100%. Otherwise, the DMU is inefficient.

Empirical Analysis
This study analyzes 14 publicly traded shipping companies, including Maersk, CMA-CGM, COSCO, and Hapag-Lloyd. Data from 30 CSCs were obtained because they now own 80% of the worldwide container vessel fleet [16,71], but 16 CSCs were excluded from the analysis due to missing data. Their primary operations are container transportation and dry bulk freight. Most of them also offer logistics and shipping-related services, and others own or operate ports. As displayed in Table 4, each CSC is regarded as a DMU. During the evaluation period, the selected CSCs in this analysis controlled more than 70% of the overall fleet capacity in the world liner shipping industry. Toward the comprehensive evaluation implementation for CSCs' performances, the aggregation of qualitative and quantitative data for efficiency analysis is proposed. The concepts of quantitative and qualitative variables are described in the forthcoming subsections. Steps of the defined method are applied in the upcoming two phases.

Phase 1: Qualitative Efficiency Analysis
It is vital to investigate the elements that influence CSCs' competitive advantage and deepen understanding of the mechanisms that underpin their influence because CSCs account for a large amount of the shipping industry's market share. When CSCs are faced with globalized competition, the following dimensions are their standard potential advantages: service quality requirements, logistics capabilities, technology upgrading, maintenance of business relationships, flexible operations, and market orientation, to name a few. Fanam et al. [72] reported that freight costs, service quality, scheduling, equipment handling, and information technology (IT) would influence container liner shipping advantages. Magnus et al. [73] claimed that economies of scale, pricing, network coordination, cost reduction, regulation, and shipper relationships were all key factors affecting the competitiveness of a CSC. When evaluating CSCs discussed by Hsu and Ho [27], other influential factors are transport reliability and security, corporate reputation, professional expertise, and integrated logistics operations. Notably, in rapidly evolving market circumstances where the focus is on green and sustainable development, social and environmental factors in the container shipping context have recently received increasing attention. For container shipping operators to improve their green performance, Yang [25] indicated that environmental conventions/directives/regulations, internal green practices, and external green collaborations are noteworthy dimensions.
According to previous research, many factors influence the competitive advantage of a shipping liner. Based on the literature review and expert's opinions, a list of criteria was selected in this study to evaluate CSCs' competitive advantage benefits (see Table 5). These factors are frequently described using difficult-to-quantify qualitative descriptions. When evaluating organizational performance in the shipping business, it can be challenging to undertake a quantitative study in some circumstances. As a result, experts are frequently asked to provide qualitative descriptions based on their knowledge and expertise. However, one of the major drawbacks of using expert panels is that their recommendations are frequently accompanied by ambiguity and hesitancy. Thus, the use of proposed methodologies AHP-SF and COPRAS-G is an effective way to handle this problem, presented in the upcoming sub-sections. For this analysis, a panel of 15 experts with at least ten years of professional experience in the maritime industry was invited to evaluate criteria and DMUs, as shown in Table 6. Some of the experts are well-experienced port state control officers, master mariners, senior lecturers, shipping administration professionals, and freight forwarders. The others are managers and executives in the shipping industry. Refers to the assurance of a safe and equitable workplace for all employees. Bao et al. [28] Internal green practices (SE2) Defined as many internal green shipping practices and operations that a CSC can implement and manage independently to reduce the environmental impacts of daily activities.
Yang [25], Di Vaio et al. [26], Lirn et al. [74] Environmental institutional pressures (SE3) The adoption and implementation of conventions, directives, regulations, and strategies on container transport to protect the environment.

Reliability (SL1)
Refers to on-time performance, responsibility display to customers, accuracy of transshipment, ability to handle cargo at the destination in safe and sound condition, and lower probability of shutting out or roll-over of containers at transshipment port.
Iqbal and Siddiqui [23], Hsu and Ho [27] Flexibility and responsiveness (SL2) Defined as how fast a shipping line is to cater and adapt to the changing needs and requirements.
Iqbal and Siddiqui [23],Čirjevskis [78] Quality of service (SL3) Refers to quality control and inspection for a variety of available and value-added services of a CSC can provide, commitment to continuous improvement.

The Use of AHP-SF for Determination Criteria Weights and Results
The AHP-SF model is employed to weight the qualitative performance criteria of CSCs. In this stage, four selected dimensions including Counterparty (C), Social and Environmental aspects (SE), Service Level (SL), Operation (O), and their sub-criteria for qualitative performance evaluation were determined as shown in Table 5. Figure 4 depicts the decision problem's hierarchical structure. Figure 4 shows a hierarchy of the significance of the criteria as determined by the decision-makers using the linguistic weighting factors.
The computation of the four primary dimensions shown below serves as an example of the AHP-SF methodology. The pairwise comparison matrix utilizing linguistic terms, the non-fuzzy comparison matrix, and the normalized comparison matrix of the four primary dimensions are all displayed in Tables 7-9. The consistency ratio of the pairwise comparison for experts was calculated as follows. Note that WSV denotes weighted sum value, CV denotes consistency vector, D denotes a considered dimension, SI denotes score index.    Dimension  AMI  VHI  HI  SMI  EI  SLI  LI  VLI  ALI   C  4  3  3  2  2  1  SE  C  3  2  2  4  4  SL  C  1  2  3  2  4  3  O  SE  3  3  3  2  1  3  SL  SE  2  4  3  2  1  3  O  SL  1  4  3  3  1 3 O With the four main dimensions (n = 4), the largest eigenvector (λ max ) was calculated to identify the consistency index (CI), the random index (RI), and consistency ratio (CR) as follows:

Left Criteria Is Greater Right Criteria Is Greater
λ max = 4.085 + 4.079 + 4.127 + 4.124 4 = 4.104 Such that n = 4, RI = 0.9, and the CR value is calculated as follows: As shown in CR = 0.038 < 0.1, the pairwise comparison matrix was consistent, and the result was satisfactory.
Following that, Table 10 calculates the integrated spherical fuzzy comparison matrix. The resulting spherical fuzzy weights for each dimension were then determined and are displayed in Table 11 as a result. For explanation, the following calculation was presented for the spherical fuzzy weights of criteria D1, which is Counterparty (C), with spherical fuzzy weights (α, β, γ) = (0.432, 0.542, 0.304), as follows:  The AHP-SF weights of the four main dimensions consist of three parameters: the membership function (α), non-membership function (β), and hesitancy function (γ) of the element x ∈ X. The crisp weights of the four main dimensions were calculated based on the abovementioned calculation. The most significantly correlated dimension to qualitative performance is Service Level (SL) with a value of 0.279, followed by Operation (O) with a value of 0.258, Social and Environmental aspects (SE) with a value of 0.237, and Counterparty (C) is the last significantly dimension with a value of 0.226. Therefore, the same procedures are used to determine the importance level for 15 criteria. Table A1 presents the integrated spherical fuzzy comparison matrix with 15 criteria (Appendix A). Table 12 shows the spherical fuzzy weights and crisp weights of the AHP-SF model. The significant level of each criterion is calculated based on the geometrical mean, defuzzification, and normalization procedures. For example, the spherical fuzzy weights of the criteria External green collaborations (C1) have a membership function (α) at 0.496, non-membership function (β) at 0.478, and hesitancy function (γ) at 0.342. Similar to the procedure, the spherical fuzzy weights of the criteria Relationship (C2) have membership function (α), non-membership function (β), and hesitancy function (γ) of 0.458, 0.529, and 0.325, respectively. The significance levels of 15 criteria of the AHP-SF model are displayed in Figure 5. The results indicate that the five most significant criteria for qualitative performance evaluation of CSCs are Quality of service (SL3), Equipment system and IT application (O4), Integrated logistics operations (O3), Flexibility and responsiveness (SL2), and Environmental institutional pressures (SE3), with significance levels of 7.89%, 7.72%, 7.58%, 7.33%, and 7.05%, respectively. Meanwhile, Professionalism (O5) is specified as the least significant criterion, with a value of 5.12% compared to other considered criteria. The findings suggest that decision-makers focus on "SL3", "O4", "O3", "SL2", and "SE3" for improving the qualitative performance of CSCs.  Next, the COPRAS-G approach will employ the crisp numbers discovered at the conclusion of the AHP-SF stage as weight values for qualitative criteria.

The Use of COPRAS-G and Results
In this step, COPRAS-G will convert qualitative inputs into a single quantitative output variable termed "expert-based qualitative performance" (EQP). Using COPRAS-G methodology, CSCs are graded according to how well they handle qualitative data. The AHP-SF model is used to determine the preference weight for each criterion. The integrated grey decision matrix of alternatives regarding criteria, in accordance with the COP-RAS-G method, is shown in Table A2 (Appendix A). Table 13 provides an evaluation of the COPRAS-G model's utility level. Maersk (CSC-05) attains the best qualitative performance with a utility degree of 100% consequently. COSCO (CSC-01) comes in second with a utility degree of 96.97%, while Evergreen (CSC-11) comes in third with a utility degree of 94.41%. Meanwhile, with a utility degree of 53.75%, Sinotrans (CSC-03) has the lowest qualitative performance. Next, the COPRAS-G approach will employ the crisp numbers discovered at the conclusion of the AHP-SF stage as weight values for qualitative criteria.

The Use of COPRAS-G and Results
In this step, COPRAS-G will convert qualitative inputs into a single quantitative output variable termed "expert-based qualitative performance" (EQP). Using COPRAS-G methodology, CSCs are graded according to how well they handle qualitative data. The AHP-SF model is used to determine the preference weight for each criterion. The integrated grey decision matrix of alternatives regarding criteria, in accordance with the COPRAS-G method, is shown in Table A2 (Appendix A). Table 13 provides an evaluation of the COPRAS-G model's utility level. Maersk (CSC-05) attains the best qualitative performance with a utility degree of 100% consequently. COSCO (CSC-01) comes in second with a utility degree of 96.97%, while Evergreen (CSC-11) comes in third with a utility degree of 94.41%. Meanwhile, with a utility degree of 53.75%, Sinotrans (CSC-03) has the lowest qualitative performance. At the end of the second step, firstly, criteria weights based on 15 qualitative criteria are calculated with AHP-SF methodology as fuzzy numbers and then calculated into crisp numbers. Secondly, according to the COPRAS-G method, qualitative variables are transformed into only one quantitative variable as an output called expert-based qualitative performance (EQP) and will be used in DEA as an output.

Phase 2: Finding the Ranking of Efficient and Inefficient CSCs with DEA and Final Results
The qualitative performance values obtained from COPRAS-G analysis are considered an output variable called EQP. In this phase, we select quantitative variables for efficiency analysis after reviewing the existing literature in the container shipping industry. Two efficiency measurement models are established for this analysis, as shown in Figure 6. First is the cargo efficiency model of CSCs, which is measured by using four inputs (owned-in and chartered-in fleet capacity, employee, and operating costs) to produce one output (lifting). Then, the eco-efficiency model with the EQP output variable is developed by using three inputs (employee, operating costs, and lifting) and three outputs (revenue, CO 2 emissions, and EQP). Data collection of input and output variables is shown in Table A3 (Appendix A). Table 14 shows the description and data sources of inputs and outputs (data values of 2020) used in the DEA analysis. Data on chartered-in fleet capacity, owned fleet capacity, employees, operating costs, lifting, revenue, and CO 2 emissions were collected from two main sources: the Alphaliner website [78] and related reports of each CSC. The statistical analysis of input and output variables is shown in Table 15. The variables have a positive connection as in Table A4 (Appendix A), with all correlation coefficients between input/output items being significant, indicating a sound isotonicity of our DEA model.
The efficiency scores and rankings of 14 CSCs in both models were obtained, as shown in Table 16, after the DEA analysis of this system was performed in the DEA Solver software package. In total, five out of 14 CSCs (Sinotrans, SITC (data values of 2020) used in the DEA analysis. Data on chartered-in fleet capacity, owned fleet capacity, employees, operating costs, lifting, revenue, and CO2 emissions were collected from two main sources: the Alphaliner website [78] and related reports of each CSC. The statistical analysis of input and output variables is shown in Table 15. The variables have a positive connection as in Table A4 (Appendix A), with all correlation coefficients between input/output items being significant, indicating a sound isotonicity of our DEA model.  Annual reports, website, and related reports of each company

Discussions and Conclusions
Measuring the efficiency of major global CSCs is a contentious issue for both academics and practitioners. More significantly, in the case of maritime transportation, businesses and governments realize that adapting to the post-pandemic world and rebuilding better entails providing economic, social, and environmental value and developing new business opportunities. Compared to most of the previous studies, mainly quantitative-oriented approaches, this study enables the inclusion of both quantitative and qualitative factors in the evaluation process for the container shipping sector. A more systematic and analytical framework for global shipping evaluation that integrates AHP-SF, COPRAS-G, and DEA has been proposed in this paper for the first time in the literature. In terms of qualitative performance analysis, the criteria weights from AHP-SF results indicated that quality of service, equipment system and IT application, integrated logistics operations, flexibility and responsiveness, and environmental institutional pressures are the most im-

Discussions and Conclusions
Measuring the efficiency of major global CSCs is a contentious issue for both academics and practitioners. More significantly, in the case of maritime transportation, businesses and governments realize that adapting to the post-pandemic world and rebuilding better entails providing economic, social, and environmental value and developing new business opportunities. Compared to most of the previous studies, mainly quantitative-oriented approaches, this study enables the inclusion of both quantitative and qualitative factors in the evaluation process for the container shipping sector. A more systematic and analytical framework for global shipping evaluation that integrates AHP-SF, COPRAS-G, and DEA has been proposed in this paper for the first time in the literature. In terms of qualitative performance analysis, the criteria weights from AHP-SF results indicated that quality of service, equipment system and IT application, integrated logistics operations, flexibility and responsiveness, and environmental institutional pressures are the most important criteria that influence CSCs' competitive advantage. COPRAS-G indicated efficient and inefficient CSCs under experts' qualitative evaluation. Expert-based qualitative performance (EQP) is considered an output variable for the DEA performance analysis. The DEA analysis indicated that only a few CSCs were efficient concerning all dimensions.
From the final results, it is found that the average cargo efficiency score of all CSCs (0.6407) is lower than the average score in the eco-efficiency model with EQP (0.7899). Our result implies that, during the COVID-19 pandemic, CSCs have had a positive impact on the environment. This can be explained by the fact that the COVID-19 pandemic has crippled the global shipping industry in its initial stages. Although preliminary calculations for 2020 show a significant fall in transportation emissions due to lower activity during the COVID-19 outbreak, there is no doubt that emissions will rebound after 2020. As a result, CSC decision-makers should focus on strengthening environmental strategies as soon as possible. Rather, it is indicated that Maersk, CMA-CGM, and Hapag-Lloyd, three of the largest shipping companies in the world in terms of capacity, have the worst efficiencies during the research period, indicating that they have not been maximizing their resources to obtain the best possible outcomes in both operational and environmental performance. Thus, there is a lot of opportunities for them to enhance efficiency performance. The improvement can be achieved by making better use of their current resources rather than intensifying the resources. Moreover, CSCs that are efficient in only one dimension (COSCO, ONE, HMM, and Matson) should dissect the weaknesses and redevelop their strategies for improving both efficiencies.
The integrated method used in this study to evaluate the container shipping sector aids in the synthesis of several criteria in various CSCs to provide sustainable performances that take operational, economic, environmental, and social factors into account. This study suggests that inefficient CSCs make strategic steps to improve their performance, which typically helps to evaluate relative efficiency in the literature on the shipping industry. Our findings can assist the managers of CSCs in identifying their respective strengths and weaknesses in running businesses. On the other hand, identifying the efficient CSCs is very important from the perspectives of shippers and logistics service providers to select appropriate shipping suppliers. Due to the incompleteness of data, our study's limitation is that we only considered 14 CSCs as a sample, despite the fact they owned more than 70% fleet capacity of the container shipping market during the evaluation period. Therefore, it is advised that researchers examine as many DMUs as they can. Additionally, due to lacking data, this study only takes CO 2 emissions into account as an unwanted output in the eco-efficiency model. Future research may incorporate other emissions (SOx, NOx, etc.) as unwanted outputs in the eco-efficiency model. Furthermore, it is possible to look at external variables affecting CO 2 and other emissions. In terms of methodology, sensitivity analysis and/or comparative analysis with other MCDM techniques are suggested for researchers to validate the results.

Acknowledgments:
The authors appreciate the support from the National Kaohsiung University of Science and Technology, Taiwan; and Hong Bang International University, Vietnam.

Conflicts of Interest:
The authors have declared that no competing interests exist.
Appendix A   Table A1. The integrated spherical fuzzy comparison matrix of the SF-AHP model.   Note: ** denotes correlation is significant at the 0.01 level (2-tailed), * denotes correlation is significant at the 0.05 level (2-tailed).