Selecting Suitable, Green Port Crane Equipment for International Commercial Ports

: Responding to the increasing global need for environmental protection, a green port bal-ances economic vibrancy with environmental protection. However, because exhaust emissions (e.g., CO 2 or sulﬁde) are difﬁcult to monitor around ports, data on such emissions are often incomplete, which hinders research on this topic. The present study aimed to ﬁll this gap in this topic. To remedy this problem, this study formulated a new data envelopment analysis (DEA) method for collecting CO 2 emissions data at their source. This method was applied to collect real-world operating data from a large container-handling company in Taiwan. Speciﬁcally, we provide a real example using a novel green energy index to account for undesirable outputs. Our main objective was to formulate two methods that combine: (1) data envelopment analysis based on a modiﬁed slack-based measure, and (2) a multi-choice goal programming approach. The contributions of this paper included the ﬁnding that rubber-tired gantry cranes are the greenest and should be used in ports. Finally, our ﬁndings aid port managers in selecting port equipment that provides the best balance between environmental protection and proﬁtability.


Introduction
Environmental degradation and resource overconsumption are serious global problems, whereas sustainable development benefits both a country and its economy. Human activity is responsible for both environmental protection and environmentally damaging economic growth. Correspondingly, although a country's natural resources (e.g., air, water, soil, and minerals) enable its development, their overexploitation is bound to backfire eventually, leaving future generations with environmental problems, such as wildlife extinction and natural resource depletion.
Human overreliance on fossil fuels has resulted in climate change, which is disruptive at best and destructive at worst. Climate change has and will destroy marine ecosystems, melt glaciers, decimate the Amazon rainforest, and trigger large-scale human migration and conflict [1]. Sea levels will also rise due to climate change, and eroded coastal conditions, the release of inundated land, and the threat of submersion will be disastrous for island nations and low-lying coastal areas. This threat is especially serious given that half of the global population lives within 100 km of a coast [2] and that coastal regions tend to be wealthy.
In response, many coastal governments have begun formulating strategies for sustainable development. Ports are a crucial driver of economic growth, but they are also energy intensive and a source of pollution. To remedy this problem and to ensure sustainable development, the concept of a green port was formulated. The move toward green ports has made considerable progress in many developed countries, as reflected in the San Pedro Bay Clean Air Action Plan (jointly implemented by the Port of Los Angeles, California, and the Port of Long Beach, New York and New Jersey), the Clean Air Initiatives and Harbor Air Management Plan (jointly implemented by port authorities in New York and New Jersey), the Rijnmond Regional Air Quality Action Program (implemented by the Port of Rotterdam, The Netherlands), and the Green Port Guidelines (implemented by the Port of Sydney, Australia).
In the context of these developments, more scholarly attention has been paid to the rational utilization of port resources [3][4][5][6][7][8]. Studies have aimed to assist port managers in formulating feasible policies from a macroscopic perspective that accounts for scaling effects and the balance between economic vibrancy and environmental protection. However, these studies have not considered the sources of environmental damage in and around ports (e.g., sources of CO 2 emissions). In response to this gap in the literature, this study focused on the container-handling system, which is closely related to daily port operations. Even though CO 2 emissions are an important consideration among port authority and terminal operators, the literature on this topic is limited. In general, few studies have focused on evaluating the environmental performance of green ports, likely because port emissions data (pertaining to, for example, CO 2 or sulfide) are difficult to collect. The present study aimed to fill this gap.
Specifically, this real-world example study combined data envelopment analysis (DEA) and multi-choice goal programming (MCGP) to evaluate the green performance of four types of cranes that are commonly used in ports. The findings will aid port managers in making their port greener. The contributions of this paper include: (1) the modified super slack-based measure (SBM-DEA) model, which is presented to evaluate the overall crane equipment efficiency, (2) the MCGP model's use in choosing crane equipment, and (3) the application of the proposed approaches to port operation in Taiwan, alongside the provision of some valuable suggestions.
The remaining parts of the paper are organized as follows: Section 2 reviews the literature on green ports. Section 3 introduces this study's combination of DEA, based on the modified SBM, and the MCGP method that accounts for undesirable outputs. A real-world case example in Taiwan is presented in Section 4. Finally, Section 5 concludes the paper and discusses the managerial implications.

Green Ports
In general, green port construction involves aspects such as improving water quality, supervising air quality, ensuring noise control, managing waste and hazardous cargo, conducting environmental education and training, and maintaining biodiversity in the port area. Scholars have researched these aspects.
In analyzing the water circulation patterns in the port of Ensenada (one of Mexico's most important ports), Espino et al. [9] suggested the use of a wave energy pumping system to gradually dilute the concentration of pollutants in the port area. Otene and Nnadi [10] focused on water quality indices and water quality conditions in the Port of Harcourt, Nigeria. Their study collected water samples from four key locations in the port and analyzed the water quality parameters using standard methods. Their findings indicated the poor state of environmental monitoring, thus aiding the port's managers. Lee et al. [11] analyzed a comprehensive 2010-2011 data set on marine environmental trends, including those of water quality, along the coast of Busan New Port. Their findings aided port managers in monitoring the impact of projects on the offshore marine environment around the port. Bolognese et al. [12] noted that, in contrast to the many studies that have investigated the management of noise from transportation, few studies have investigated the management of noise from port operations. The authors investigated the North Tyrrhenian Sea Port by collecting data from monitoring systems, noise measurements, and citizen complaints. Their findings indicated a neglect of noise levels by port managers. Reviewing the regulations and literature on environmental issues in port management systems, Vaio et al. [13] conducted semi-structured interviews with users of an Italian port to explore how port management control systems assist port authorities in the decision-making process. To help port managers improve management efficiency during ship mooring, their study also assessed the efficiency of port waste management.
Focusing on official regulations, Prati et al. [14] investigated the air quality in the Port of Naples through two experiments. Measurements were made at 15 points within the port. In addition, a laboratory was established within the port area to take continuous measurements of pollutant concentrations, ambient parameters, particulate matter (PM) levels, and wind direction and intensity. Their findings indicated that ship emissions contributed the most to SO 2 concentrations as compared to the concentrations of other pollutants. Kontos et al. [15] focused on the impact of gas emissions from cruise ships and passenger vessels on air quality and human health risks in the area around the Port of Thessaloniki. They estimated the surface concentration of pollutants caused by passenger ship traffic using the CALPUFF dispersion models for 2013, and their study also forecast trends for future environmental conditions within the port area. Casazza et al. [16] used 3D modeling to achieve the effective regulation of air quality within a port area. Their study not only enabled air pollution monitoring in ports, but also provided a new methodology in support of local environmental management systems. Progiou et al. [17] demonstrated that navigation emissions from ships are an important component of the total emissions, whether of a port, port city, or country. Their study used atmospheric models to simulate the dispersion of air pollutants, and their findings indicated a significant increase in activity in the Port of Piraeus over the last decade, especially from merchant ships.
As evident in the preceding literature review, studies have typically monitored the environment around port areas through monitoring stations, thus gaining a macro-level understanding [18][19][20][21][22][23]. In this study, we classified research topics on green performance evaluation methods, and several representative publications are listed in Table 1. Table 1. Variables selection of the recent studies for the green performance evaluation. Few studies have monitored greenhouse gas emissions at their source. The cranes in a port are one such source, as they emit greenhouse gases when continually loading and unloading cargo. Therefore, the construction of an effective evaluation approach for selecting environmentally friendly cranes is a research problem of practical importance, and it is this problem (and gap in the literature) that this study aimed to address.

DEA Applied in Green Ports
Among the many existing methods for evaluating performance, DEA is well known by many managers and researchers alike because of its unique advantages in processing multiple inputs and outputs. The conventional DEA model was first proposed by Charnes et al. [27] in 1978. It was based on linear programming, which is a quantitative method of evaluating the relative effectiveness of comparable units of the same type. As DEA became methodologically more sophisticated with time, it has developed into a new field that integrates operations research, management science, and mathematical economics. Subsequently, Banker et al. [28] extended the DEA model to cover variable returns to scale (VRS). Since then, DEA models have been extended to other practical domains in the form of super-efficiency models [29][30][31], cross-efficiency models [32,33], SBM models [34,35], super-SBM models [36,37], and network DEA models [38][39][40].
Although DEA methods have often been used to evaluate performance with respect to CO 2 emissions [41][42][43][44][45][46], few have applied DEA to green ports specifically. Using an inseparable input-output SBM-DEA model, Na et al. [35] analyzed how environmentally friendly eight major container ports in China were by using 2005-2014 environmental monitoring data. Their results indicated that the eight ports significantly differed in their CO 2 emission levels and that their pure technical environmental efficiency was low. Li et al. [45] noted that the rapid development of China's port industry has led to serious problems with CO 2 emissions. Specifically, the authors analyzed 2013-2018 data on 16 Chinese port companies; the ports were segmented by size and complexity criteria in the analysis. Using an improved non-radial directional distance function, they determined the performance of these ports with respect to CO 2 emissions. Wang et al. [46] constructed three DEA models to evaluate the environmental efficiency gained by cooperation between ports under the conditions of environmental control, non-environmental control, and PM emissions. They collected and analyzed data from 11 major Chinese ports and found that ports in the eastern region of China performed the best with respect to environmental friendliness.

Methodology
In the present study, a two method data envelopment analysis, based on a modified slack-based measure and a multi-choice goal programming approach, was used to evaluate the performance of green port crane equipment.

Modified SBM-DEA Model
The modified SBM-DEA model is described via the sets of indices, variables, and parameters presented in Table 2.  respectively, of the jth DMU(j = 1, . . . , n). The production possible set (PPS) given by the DMUs is as follows: where v i and u r are nonnegative intensity vectors, indicating that the preceding definition corresponds to a situation of constant returns to scale (CRS). The original DEA-CCR model proposed by Charnes et al. [27] is a nonlinear programming model, which traditionally analyzes all positive data. Through the Charnes-Cooper transformation [47], the efficiency of DMU-k can be formulated as follows: Equation (1) is the basic DEA-CCR model in multiplier form. The dual model presented in the envelopment form is as follows: Subsequently, Banker et al. [28] extended Equation (2) to cover VRS. However, the two radial approaches may be limited by some of the inefficient components not reflected in the measurement results (such as the mix inefficiencies). To address this problem, Tone [34] proposed the following SBM model: where s − i , s + r denote the inefficient components. In Equation (3), Tone [34] defined the evaluated DMU to be efficient if, and only if, the optimal solution of s − * i = s + * r = 0 for all i and r (or equivalently, the efficiency ρ * = 1). To further enhance the discrimination In Equation (4), the new PPS can be defined as: Note that for the inefficient DMUs, the efficiency evaluated by Model (4) is necessarily one. That is, Model (4) is only effective for distinguishing between efficient DMUs. Thus, applications typically use Model (3) and Model (4) in combination.
Fang et al. [49] noted that Model (4) does not incorporate slacks explicitly, and they suggested adding two slack variables (w − i , w + r ) to account for the incorporated slacks of the first two constraints of Model (4). Furthermore, because our variable of CO 2 emissions was considered an undesirable output in this study, referencing Fang et al. [49], we supposed that n DMUs obtain m inputs, s outputs, and g undesirable outputs. Let three vectors x i ∈ R m , y r ∈ R s and u h ∈ R g (h = 1, . . . , g) denote m, s, and g, respectively. Correspondingly, we can obtain the matrices X, Y and U as follows: Note that, because all the research data are nonnegative, we obtained X > 0, Y > 0 and U > 0. The new PPS can be defined as follows: where the intensity vector λ ∈ R n , and the preceding definition of PPS, corresponds to the CRS in envelopment form.
In fact, the original SBM-DEA model involved calculating the ratio of the average input reduction to the average output growth when evaluating the efficiency. In other words, the purpose of the objective function of the SBM-DEA model is to determine the most appropriate extent of improvement between inputs and outputs. Thus, the SBM-DEA model can be referred to as a non-radial model or non-oriented model. One advantage of this model is that it allows the analyst to evaluate the efficiency by analyzing the maximum adjustable quantity of each vector instead of only analyzing the improvement of one dimension (inputs or outputs) alone. In this study, we aimed to minimize both the inputs and undesired outputs. Therefore, we propose the following model to evaluate the super-efficiency: where w − i , w + r , and w − h denote the incorporated slacks (or super-efficient components) of inputs, good outputs, and undesirable outputs, respectively. In Equation (6), the constraints w + r ≤ y rk (r = 1, . . . , s) and w − h ≤ u hk (h = 1, . . . , g) ensure that the computed superefficiency value is always nonnegative.
Similar to Equation (4), Equation (6) is such that when DMU-k is located outside of the new PPS (5), the efficiency value of DMU-k is greater than 1; this DMU is then evaluated as an efficient unit. In other words, Equation (6) can determine the minimum distance (w − i , w + r , and w − h ) between the efficient frontier and the evaluated DMU. However, for any evaluated DMU-k that falls within the region of the new PPS (5), the minimum distance (w − i , w + r , and w − h ) is necessarily zero; that is, Equation (6) cannot determine the gap between the evaluated DMU and its target. Thus, in this study, we propose the following model to calculate the efficiency of inefficient DMUs: where w − * i , w + * r , and w − * h are the optimal solutions that are calculated using Equation (6), and the optimal solution of the new variables s + * i , s − * r , and s + * h denote the inefficient components of the evaluated DMU. Therefore, we formulate efficiency as follows: In this study, to determine the optimal loading tool that has satisfactory green performance, we further define a new green energy index (GI j ), which is obtained by first calculating the super-efficiency value, DMU-j(j = 1, . . . , n), before calculating the maximum value, E max = n max j=1 φ * j . Finally, the green energy index GI j can be calculated as follows: In order to ensure the possibility of the proposed approach, this study provides an unambiguous flow chart in Figure 1. The flow chart for this DEA application is as follows: Sustainability 2021, 13, x FOR PEER REVIEW 9 of 21 In order to ensure the possibility of the proposed approach, this study provides an unambiguous flow chart in Figure 1. The flow chart for this DEA application is as follows:

Modified SBM-DEA Model for Evaluating Green Energy Index Process
Step 1: First, determine the research subjects of the study; Step 2: Select appropriate variables (green energy index) through the literature review; Step 3: Construct the suitable modified SBM-DEA model for the performance evaluation; Step 4: Based on the research results of previous steps, this study should interview some relevant enterprises to investigate the practical operational data of the selected variables; Step 5: Use the new proposed model (such as the new Equation (6) and (7) in this study) to conduct performance evaluation processing on the collected data. If some infeasibility problem occurs, it is necessary to go back to Step 3 and make appropriate adjustments to the constructed modified SBM-DEA model. After doing so, repeat Step 5; Step 6: Each evaluated DMU can obtain its super-efficiency ) ,..., 1 ( * n j j =  with Formula-

Modified SBM-DEA Model for Evaluating Green Energy Index Process
Step 1: First, determine the research subjects of the study; Step 2: Select appropriate variables (green energy index) through the literature review; Step 3: Construct the suitable modified SBM-DEA model for the performance evaluation; Step 4: Based on the research results of previous steps, this study should interview some relevant enterprises to investigate the practical operational data of the selected variables; Step 5: Use the new proposed model (such as the new Equations (6) and (7) in this study) to conduct performance evaluation processing on the collected data. If some infeasibility problem occurs, it is necessary to go back to Step 3 and make appropriate adjustments to the constructed modified SBM-DEA model. After doing so, repeat Step 5; Step 6: Each evaluated DMU can obtain its super-efficiency φ * j (j = 1, . . . , n) with Formulation (8); Step 7: Then, determine the maximum value E max = n max j=1 φ * j to calculate the green energy index GI j (j = 1, . . . , n) for each evaluated DMU via Formulation (9). Based on the results of the green energy index, obtain performance ranking results and screen out the efficient benchmarking DMUs; Step 8: Conduct an in-depth analysis of the most suitable adjustment for the each DMU; Step 9: Finally, check the modified SBM-DEA model evaluation results through the MCGP method (Equations (10)- (14)).

MCGP Model for Choosing Suitable Crane Equipment
The MCGP approach encompasses the many modified GP methods in the literature. Chang (2008) developed a multi-choice aspiration level model for solving multi-objective decision-making problems [50]. A typical MCGP problem has the following structure.
In the real-world decision-making problem of choosing crane equipment, the goals are often related. This problem is represented in the following MCGP equations: Minimize subject to As illustrated in Equations (11)-(13), selection restrictions are absent for any single goal, but some goals are dependent on another. For example, we can add the auxiliary constraint b i ≤ b i+1 + b i+2 to the MCGP model, where b i , b i+1 and b i+2 are binary variables. Thus, b i+1 or b i+2 must be equal to 1 if b i = 1. This means that if goal one was achieved, then either goal two or goal three was also achieved.

Data Collection and Evaluating Green Performance of Various Cranes
This real-world case study aimed to evaluate the green performance of various cranes used to load and unload cargo in port operations. The four most common types of cranes used in international commercial ports, both in general and by the prominent container-handling company in Taiwan in particular, are as follows: gantry cranes (GCs), rail-mounted gantry (RMG) cranes, rubber-tired gantry (RTG) cranes, and empty container handlers (ECHs). This study collected and analyzed the 2018-2020 data on these cranes (Tables 3-5).

Modified SBM-DEA Model Analysis Result
In general, the selection of input and output variables is critical in the application of DEA. This is because the evaluation results become highly variable when the set of research variables change. The Taiwanese company investigated in this study was large and operated many cranes (including nine RMG cranes). The data for all cranes of each type also differed little. Thus, the data used in this study were the average values for each crane type.
In Tables 3-5, the basic information on each crane is presented from the second to sixth columns from the left, and the performance values, as computed using Equations (6) and (7) jointly, are presented in the seventh column from the left. The penultimate and final columns present the value of the green energy index (GI j ) and the ranking for all four crane types, respectively.
The results indicated that the green performance ranking among the cranes differed little from 2018 to 2019, and that the efficiency value of three crane types (RTG, GCs, and ECHs) exceeded 1. Thus, these three crane types operated efficiently throughout the years, with RTG cranes having the best green performance and being the most efficient. In 2020 (Table 4), in contrast to previous years, RMGs and ECHs switched rankings and the green performance of ECHs was determined to be inefficient. RTG cranes still had the best, and thus most stable, green performance and were found to be optimal for use in global commercial ports.

Suitable Adjustment for Each Variable Analysis
Simply providing port managers with the green performance of the four types of evaluated cranes can judge whether their operations are efficient or not, but it fails to inform them of the quantitative analysis of the advantages and disadvantages for each evaluated crane. In order to explore the advantages of efficient cranes and the disadvantages of inefficient ones, this study conducted an in-depth analysis of the most suitable adjustment for the four evaluated cranes. Then, port managers can effectively address or improve operational weaknesses, which can highlight the practical management significance and value of this study. The corresponding analysis results are presented in Table 6. In this table, the positive values of the change rate denote the advantage of each variable, while the negative values denote the disadvantage. For example, GCs produced an average of 278,928 kg of CO 2 emissions in 2019, but its relative benchmark value should have been 147,983 kg for that year and, as such, the decrease change rate was −46.95%. In other words, the poor performance of the GC was due to the large waste of resources in terms of carbon dioxide emissions, which is also an important reason for its relative inefficient performance. Therefore, managers should strengthen the scientific supervision at this level in future operations.  Table 6 presents the quantitative results for the tradeoff among the variables for each crane type. The results indicated the target that should be learned for each variable in a given year and the extent of adjustment (expressed in terms of an adjustment ratio) for each variable in the optimal tradeoff. For the input and undesired outputs, the adjustment ratio was calculated by subtracting the original resource value from the target value, and then dividing this difference by the original values. A positive adjustment ratio represented the performance of the learning benchmark in that direction as being not yet as good as that of the evaluated unit. In other words, a positive adjustment ratio can be interpreted as representing the advantage for a given crane type.
Conversely, if the value of the adjustment ratio for an item is negative, it represents a disadvantage for a given crane type. For good output variables, this study used reverse processing, in which the original data value was subtracted from the target value and this difference was divided by the target value. This was done to allow positive numbers to also represent advantages. Table 6 presents the adjustment ratios for all crane types. The RTG was the best crane type with respect to all variables, especially in energy consumption and total energy cost, with average three-year advantages of 14.74% and 11.81%, respectively. The GC was the second-best crane type, and it was superior primarily in operational duration. Thus, the GC is especially advantageous when used to load and unload the same type of containers. Finally, RMG cranes and ECHs were disadvantageous due to their high energy consumption and high total energy cost. Between the two, ECHs emitted less CO 2 and had a better operational duration. These results are visualized in  type with respect to all variables, especially in energy consumption and total energy cost, with average three-year advantages of 14.74% and 11.81%, respectively. The GC was the second-best crane type, and it was superior primarily in operational duration. Thus, the GC is especially advantageous when used to load and unload the same type of containers. Finally, RMG cranes and ECHs were disadvantageous due to their high energy consumption and high total energy cost. Between the two, ECHs emitted less CO2 and had a better operational duration. These results are visualized in Figures 2-5.         In Figures 2-5, which each present the adjustment ratios for a given crane type for all variables, the solid line segment indicates the average value of the adjustment ratio for each year. As mentioned previously, positive and negative values indicate advantages and disadvantages, respectively. The characteristic patterns presented in these four figures remain largely consistent with those highlighted by the average evaluation results. From , the most obvious advantages and disadvantages of each evaluated DMU are clear. For example, the biggest advantages of RTG cranes are shown in the X2 and X3 levels, while the other three evaluated cranes expose the potential problem of resource waste in these two investments. Meanwhile, only RTG cranes and ECHs had certain advantages for U1 (CO 2 emission volume (kg)).

Using MCGP to Solve the Case Example of Choosing Between Crane Equipment
To solve the crane equipment selection problem of choosing between crane types, the analyst must define the MCGP model according to the following goals. In this case, suppose that the decision maker had the following set of priorities derived from the DMU results for RTG cranes in Table 5 s collected data and evaluation results for 2020 and Table 6's suitable adjustment for each variable RTG crane value.
To solve the problem of choosing between crane equipment problems, it is necessary to define the MCGP model according to the following goals: 1.
The first goal is Y1: the working capacity is the RTG benchmark. The DMU of the RTG was 71,341 and 77,989 in the Table 6 RTG 2020 results; 71,341 was the benchmark value and 77,989 was the AVE value; 2.
The second goal is U1: the emission volume is the RTG benchmark. The DMU of the RTG was 88,143 and 92,889 in the Table 6 RTG 2020 results; 3.
The third goal is X1: the operational duration is the DMU of the RTG crane's input. The DMU of the RTG was 3397 and 4014 in the Table 6 RTG 2020 results; 4.
The fourth goal is X2: the energy consumption is the DMU of the RTG crane's input. The DMU of the RTG was 153,429 and 186,553 in the Table 6 RTG 2020 results; 5.
The fifth goal is X3: the total energy cost is the DMU of the RTG crane's input. The DMU of RTG was 504,154 and 610,634 in the Table 6 RTG 2020 results. Table 7 provides the coefficient and goal values for solving the MCGP crane equipment selection problem. The MCGP model notation is introduced as follows: Indices: i 1,2, . . . , n index of crane equipment type; j 1,2, . . . , j index of deviation corresponding to the goals. Parameters: U1 j is the CO 2 emission volume of crane equipment type j, where j =1, 2, 3, 4; Y1 j is the working capacity of crane equipment type j, where j =1, 2, 3, 4; X1 j is the energy consumption of crane equipment type j, where j =1, 2, 3, 4; X2 j is the energy consumption of crane equipment type j, where j =1, 2, 3, 4; X3 j is the total energy cost of crane equipment type j, where j =1, 2, 3, 4; d + j , d − j are the maximum and minimum deviation of goal j; e + t , e − t are the maximum and minimum deviation of y t − g i,max/min . Decision variables: X i is the order quantity of vendor i; y i is the binary integer 1 i f the order is o f f ered by crane euipment i 0 otherwise We provide the following MCGP model programming in Table 8: Table 8. MCGP model solution programming for the crane equipment section.

Conclusions
Green ports are becoming increasingly prominent with the increased need for environmental protection globally. However, few studies have monitored exhaust or PM emissions (such as CO 2 or sulfide) around ports due to the difficulty of doing so, and the data obtained are incomplete.
To fill this gap in the literature, this study measured CO 2 emissions at their source, specifically container-handling cranes (which are indispensable to port operations). Five key variables, including CO 2 emissions, were identified based on consultations with experts. Subsequently, (1) we applied a method that combined a modified SBM-DEA model with the MCGP method to account for undesirable outputs, and (2) we defined a novel green energy index to evaluate green performance. Our findings determined (1) the crane type with the best green performance, and (2) how advantages and disadvantages are balanced in the use of each crane type. These findings can help port managers select the best crane equipment that makes their port greener, smarter, and more profitable.

Managerial Implications
We present the following managerial prescriptions based on our findings. First, we recommend RTG cranes because they are the most environmentally friendly when used in international commercial ports and strike the best tradeoff between environmental protection and profitability. Second, RMG cranes and ECHs consume considerable energy, which constitutes a point of concern that port managers must pay attention to. Third, to mitigate environmental harm and commercial loss, port managers should replace outdated equipment or, if they are unable to do so, supervise outdated equipment more intensely. Fourth, port managers can invest more in researching and developing smarter port equipment, which incorporates, for example, big data or Internet of Things technology. Smart port equipment minimizes operational waste and, thus, mitigates their environmental impact and enhances profitability.

Limitations
(i) To mitigate the disadvantages of the DEA method, we used the MCGP method to verify the DEA results. To better cope with uncertainty, decision makers can use the novel MCGP method in conjunction with the multi-criteria decision-making approach; (ii) Our study, a real-world case example, used data collected in Taiwan.

Future Directions
Future studies can use other new DEA methods to solve crane equipment selection problems. Additionally, other mathematical models, such as new MCGP models, can be combined with our study's model.