Basic Ship-Planning Support System Using Big Data in Maritime Logistics for Simulating Demand Generation

Dynamic changes in the global market demand affect ship development. Correspondingly, big data have provided the ability to comprehend the current and future conditions in numerous sectors and understand the dynamic circumstances of the maritime industry. Therefore, we have developed a basic ship-planning support system utilizing big data in maritime logistics. Previous studies have used a ship allocation algorithm, which only considered the ship cost (COST) along limited target routes; by contrast, in this study, a basic ship-planning support system is reinforced with particularized COST attributes and greenhouse gas (GHG) features incorporated into a ship allocation algorithm related to the International Maritime Organization GHG reduction strategy. Additionally, this system is expanded to a worldwide shipping area. Thus, we optimize the operationlevel ship allocation using the existing ships by considering the COST and GHG emissions. Finally, the ship specifications demanded worldwide are ascertained by inputting the new ships instance.


Introduction
Economic growth creates dynamic changes in the global market, including the maritime industry. On a worldwide scale, huge demand for industrial cargo transport promotes intense ship movement and replacement [1]. As the demand for newbuilding specifications also changes according to the circumstances, developing a ship with adequate specifications to satisfy this demand is essential.
Meanwhile, the use of digital infrastructure to alter a business model and implement value-producing opportunities naturally generates a large data stream called big data. With the reduction in cost of data collection tools, a large amount of data can be obtained from various sources and formats [2,3]. This is significant for a broader understanding of the current and future conditions of various industries. Therefore, big data analytics will be a critical advantage in the future [4].
In the maritime industry, big data are being generated through advancements in navigation systems [5,6]. Together with the voyage data recorder, an automatic identification system (AIS) is required by the International Convention for the Safety of Life at Sea to aid navigation and avoid collisions of ocean-going ships. Towards its development, the deployment of satellite-based AIS receiver enables an accurate ship's geospatial monitoring on a worldwide scale [7]. From the ship side, an AIS transponder also transmits a ship's identification number, position, course, speed, and destination. These systems uphold the digitalization of previously analog-stored data, such as ship specifications, port limitations, and sailing routes [8].
The collection of static and dynamic big data in maritime logistics has allowed various studies to be conducted. Safety improvement, energy efficiency, logistics optimization, 2 of 23 and predictive analysis have been broadly discussed [9][10][11][12]. Similarly, many data-driven studies have introduced a demand forecasting application for both regional and global scales. Forecasting analyses, such as the cargo throughput and shipbuilding market, have also been presented in some studies [13][14][15][16][17][18][19][20][21][22][23][24][25]; however, the demand for new ship specifications is unlikely to have been covered.
With big data in maritime logistics, ship operation monitoring is becoming relevant. Besides, shipbuilding tends to apply risk-based design rather than rules-based design, with aims of compatibility of design and performance [26,27]. Likewise, the International Maritime Organization (IMO) greenhouse gas (GHG) reduction strategy [28,29] set out a new future guideline. Therefore, it is crucial to examine the actual ship operation characteristics in the particular route to optimize its future ship design, both costand GHG-effectively.
Our prior studies examined the demand for new ship specifications by proposing a basic ship-planning support system using big data in maritime logistics [30]. We built the system and conducted the simulations in the form of a basic study. The proposed system was applied to the target ship of Capesize dry bulk carrier, which operated on relatively fixed routes, such as Australia and Brazil to East Asia routes [31]. Assuming the target ship operated in a time-charter contract manner, we built an algorithm to replicate its ship bidding scheme. However, the scope of that study was limited by the target routes and ship allocation considerations. The previous system delivered the simulations by only considering the ships' fuel costs; by contrast, this study has proposed an enhancement to our basic ship-planning support system. Additionally, we have broadened the scope of the geographical area to a global scale to understand the ship specifications in demand. Finally, we suggest two attributes to be considered in the ship allocation algorithm: ship cost (COST) aspects and GHG emissions.

Related Studies Utilizing Big Data in Maritime Logistics
In recent years, studies utilizing big data in maritime logistics have emerged as a tool for predictive analytics. These have mainly discussed energy efficiency improvement, accident avoidance, and logistics optimization. Focusing on both regional and global scales, some studies have examined the forecasting of the cargo throughput, shipbuilding market, and new ship specifications [9][10][11][12]. Table 1 lists several studies on these topics.
Regarding container throughput forecasting, Yang and Chang [13], Jugović et al. [14], Akar and Esmer [15], and Gökkuş et al. [16] constructed a long-term model predicting cargo demand based on regression analysis and neural networks using a combination of ship and port data. Moreover, Jia et al. [17] evaluated a regression model by inputting the AIS data of bulk carriers to estimate its payload and measure the future global carbon footprint per transport mode. Similarly, Zhou and Hu [18] examined an AIS-based estimation of iron ore volume. The actual shipment was obtained by inputting the static and dynamic data into the constructed back-propagation neural network. In addition, Kanamoto et al. [19] discussed the applicability of the AIS data of dry bulk carriers to forecast future shipping demand using a logit model.
For the case of shipping and newbuilding markets, Prochazka et al. [20] analyzed the effectiveness of AIS data to comprehend the behavior of the crude oil tanker spot-market charterer using quantile regression. Sharma and Sha [21] formulated a newbuilding index forecast based on neural networks and genetic algorithms using statistical data. Similarly, Wada et al. [22][23][24] suggested a system dynamics model to forecast the shipbuilding market demand, together with long-term GHG emissions reduction measures. Finally, Lee and Jung [25] built a platform to collect big data in shipbuilding and forecast the ship order quantity of container ships and bulk carriers, adopting an autoregressive model. The abovementioned studies have discussed a macro-level forecasting model for fluctuations in the shipping and shipbuilding markets. However, they have not examined the operation-level demand forecasting to understand competitive ship specifications following market changes.

Optimization Studies in Ship Deployment and Contracts
Several studies proposed an optimization method in ship deployment and contracts. Zhang et al. [32] analyzed the cold chain shipping mode selection between containerized and bulk reefers. Adopting the value-based management, each ship is deployed considering the economic and environmental objectives. Similarly, Venturini et al. [33] discussed the berthing times and positioning optimization in container terminals. A novel mathematical formulation is introduced to integrate the berth allocation problem and sailing speed optimization. It aims to investigate the trade-off between fuel consumption reductions and dwelling time extension. Up to 40% potential GHG emissions reduction is reported compared with when the ship sails at its design speed. These studies reported the correlation between the ship's sailing speed and GHG emissions.
Arslan and Papageorgiou [34] adopted multi-stage stochastic programming to deal with the bulk carrier renewal problem, focusing on ship-sizing and deployment. It was reported to provide the total cost reduction. Lin and Liu [35] developed a genetic algorithm to solve the tramp shipping routing problem of Handymax dry bulk carrier. The combined mathematical model was proposed to overcome ship allocation and cargo movement problems. These studies formulated a novel method to optimize the operation-level ship allocation. However, a particular case study of an actual ship allocation was not considered.
Yang et al. [36] discussed the diverse applications of AIS data, including the ship performance evaluation. Dynamic route planning in the Baltic Sea Region to reduce ship owners' COST and GHG emissions was examined [37]. The reported benefit can surpass the cost in the proposed cost-benefit analysis by reducing sailed distances and GHG emissions. It proved the applicability of AIS data to evaluating the ship allocation in the specific region. Dynamic optimization was specified as a future research opportunity to attain the optimized ship allocation considering various routes and external factors.
Bai et al. [38] examined the potency of financial hedging and operational risk management strategies to overcome significant earning risks in the tramp shipping scheme. Ordinary least square regression and the Bayesian belief network were proposed to take into account the data of 31 tramp shipping companies, such as ship age, sailing distances, and AIS data. Operational risk management strategies were reported to mitigate bunker price and freight rate fluctuation risks effectively. It confirmed that ship allocation optimization is attainable utilizing big data in maritime logistics.

Originalities in This Study
According to our literature review, no study has been carried out on demand forecasting at the particular operational level to recognize the benefit of new ship specifications in the future. Therefore, Arifin et al. [30] proposed a basic ship-planning support system by integrating big data into maritime logistics (see Table 1). We have explored the possibility of simulating operation-level ship allocation and successfully built a ship allocation algorithm as the core of our system. The following improvements have been suggested to the aforementioned study: • Target route and cargo type: Prior simulations [30] were executed on limited target routes and cargo types, such as the Australia, Brazil, Japan, and South Korea routes of iron ore. In contrast, this study has broadened the scope of simulations to accommodate iron ore and coal routes worldwide. • Ship allocation algorithm: The key performance index (KPI) in the prior ship allocation algorithm [30] considered only the fuel cost as the COST variable. Therefore, the previously generated ship allocation simulation directly reflected the estimated fuel cost, which was a fraction of the actual COST. In this context, this study has proposed more realistic COST attributes [1], such as the costs of fuel, operation, depreciation, and stockpiling.

Ship Movement Data
This dataset presents dynamic data of the ship movement and its attributes. In this study, the port calling data of ships obtained from the Market Intelligence Network (MINT) database of IHS Markit Ltd. [39] are defined as the ship movement data, which provide the historical ship position in a port-level manner for various ship types. The data acquisition between the port calling data of MINT and the operation reports of the ship's operator have been evaluated by Wada et al. [40]. More than 95% of the ship operation reports from 2015 to 2018 are covered in the ship movement data.
We then extracted the ship movement data of the Capesize dry bulk carriers, which are dry bulk carriers with a deadweight tonnage (DWT) of more than 100,000 [31]. The time range of the data was from January 2015 to May 2019. Initially, the ship movement data contained 328,670 entries, including calling types of anchorage, port, and terminal, which enclose the ships' IMO and name, country of origin and destination, port of origin and destination, dates of arrival and departure, and draughts of arrival and departure. In this study, we aggregated the data into voyage records. They were defined as the set of a ship origin port's departure data and destination port's arrival data. Thus, any operation between these entries was considered an outlier by referring to Arifin et al. [30]. Later, the entry numbers were reduced to 36,849 entries of laden and ballast.

Ship Data
This dataset presents static data of the ships' specifications. The ship data from the Sea-web Ships database of IHS Markit Ltd. (London, UK) were used in this study [41]. We downloaded the ship data of the Capesize dry bulk carriers (DWT of 100,000 or more) [31], for ships built from 1962 to 2021. This dataset contained 3253 entries, including the technical and non-technical attributes of ships, such as the IMO, name, built year, displacement, DWT, length, breadth, draught, design speed, main engine power, operator, technical manager, and registered owner.

Port Data
This dataset presents static data of the port and its specifications. The ship data from the Sea-web Ports of IHS Markit Ltd. and the IHS Fairplay Ports and Terminals Guide were used in this study [42,43]. We downloaded the worldwide port data served by dry bulk carriers containing 2529 entries, which included the variables of the ports' latitude-longitude, size limitations of served ships, and handled cargo type.

Route Data
This dataset presents static data of the routes and its sailing distances. The distance table of the IHS Fairplay Ports and Terminals Guide was used in this study [43]. We extracted the worldwide sailing distances served by dry bulk carriers containing 17,030 entries.

System Configuration
Three distinct models and a ship allocation algorithm were assembled in our basic ship-planning support system, as shown in Figure 1. We modeled a global network model, cargo movement model, and ship model using the available big data in maritime logistics. The global network model defines the ports worldwide and their attributes in its sailing network. The cargo movement model covers an estimation of cargo movement between ports. The ship model includes each ship and its specifications and predicts the ships' operation conditions.

Port Data
This dataset presents static data of the port and its specifications. The ship data from the Sea-web Ports of IHS Markit Ltd. and the IHS Fairplay Ports and Terminals Guide were used in this study [42,43]. We downloaded the worldwide port data served by dry bulk carriers containing 2529 entries, which included the variables of the ports' latitudelongitude, size limitations of served ships, and handled cargo type.

Route Data
This dataset presents static data of the routes and its sailing distances. The distance table of the IHS Fairplay Ports and Terminals Guide was used in this study [43]. We extracted the worldwide sailing distances served by dry bulk carriers containing 17,030 entries.

System Configuration
Three distinct models and a ship allocation algorithm were assembled in our basic ship-planning support system, as shown in Figure 1. We modeled a global network model, cargo movement model, and ship model using the available big data in maritime logistics. The global network model defines the ports worldwide and their attributes in its sailing network. The cargo movement model covers an estimation of cargo movement between ports. The ship model includes each ship and its specifications and predicts the ships' operation conditions. The ship allocation algorithm reconstructs the ship allocation, which is the core of our basic ship-planning support system. An input of the cargo movement demand, route list, and existing ship list is needed to observe the changes in ship allocation by adding future scenarios, such as presenting the new ships in the ship list. Additionally, we can identify new ships that are replacing existing ships from their results, which makes it possible to assess competitive ships in demand.

Global Network Model
This model contains the ports' specifications, sailing routes, and sailing distances. We extended its scope to include a worldwide sailing network served by Capesize dry bulk carriers. Thus, the model contained 564 ports and 1042 sailing routes. Additionally, we characterized the ports into main ports and port clusters.

Cargo Movement Model
This model provides the cargo movement and type between the ports. We used the arrival draught in the ship movement data to estimate the cargo amount. Later, this model results in various cargo movements arranged in origin-destination (OD) tables and is used as the cargo amount to be transported in the ship allocation algorithm.

Cargo Movement Model
Cargo movement between ports

Ship Model
Ships' specification and operations

Ship Allocation Algorithm
Allocation Algorithms COST Calculation GHG Calculation The ship allocation algorithm reconstructs the ship allocation, which is the core of our basic ship-planning support system. An input of the cargo movement demand, route list, and existing ship list is needed to observe the changes in ship allocation by adding future scenarios, such as presenting the new ships in the ship list. Additionally, we can identify new ships that are replacing existing ships from their results, which makes it possible to assess competitive ships in demand.

Global Network Model
This model contains the ports' specifications, sailing routes, and sailing distances. We extended its scope to include a worldwide sailing network served by Capesize dry bulk carriers. Thus, the model contained 564 ports and 1042 sailing routes. Additionally, we characterized the ports into main ports and port clusters.

Cargo Movement Model
This model provides the cargo movement and type between the ports. We used the arrival draught in the ship movement data to estimate the cargo amount. Later, this model results in various cargo movements arranged in origin-destination (OD) tables and is used as the cargo amount to be transported in the ship allocation algorithm.

Ship Model
This model depicts the ships' specifications, which covers the technical and nontechnical variables, such as the principal dimensions and ownership information. In addition, this model generates the operation conditions, draught rate, average sailing speed, and port staying time when a ship is allotted to a particular route. This model also predicts the operation conditions for each ship, which were later inputted into the GHG and COST calculations.

Ship Allocation Algorithm
This model computes the COST and GHG emissions when a ship is designated to a specific route. Then, this model allocates a ship onto a route with the highest merit to transport a certain cargo amount in a time-charter contract manner by adopting a greedy algorithm that considers the COST-and GHG-optimized ship allocation. Thus, this model optimizes the ship allocation and the global COST and GHG emissions, which are visualized in a region-level manner to understand the changes in the new ship allocation.

Global Network Model
This model defines each port and route served by the Capesize dry bulk carriers. We extracted a list of ports visited by such carriers from the ship movement data on a worldwide scale. The various ports were divided into main ports and port clusters; the former represented each port with significant port calling numbers, and the latter included all the other ports served by the Capesize dry bulk carriers. The port and route network contains 52 main ports and 541 other ports as 79 port clusters, as shown in Figure 2.

Ship Model
This model depicts the ships' specifications, which covers the technical and non-technical variables, such as the principal dimensions and ownership information. In addition, this model generates the operation conditions, draught rate, average sailing speed, and port staying time when a ship is allotted to a particular route. This model also predicts the operation conditions for each ship, which were later inputted into the GHG and COST calculations.

Ship Allocation Algorithm
This model computes the COST and GHG emissions when a ship is designated to a specific route. Then, this model allocates a ship onto a route with the highest merit to transport a certain cargo amount in a time-charter contract manner by adopting a greedy algorithm that considers the COST-and GHG-optimized ship allocation. Thus, this model optimizes the ship allocation and the global COST and GHG emissions, which are visualized in a region-level manner to understand the changes in the new ship allocation.

Global Network Model
This model defines each port and route served by the Capesize dry bulk carriers. We extracted a list of ports visited by such carriers from the ship movement data on a worldwide scale. The various ports were divided into main ports and port clusters; the former represented each port with significant port calling numbers, and the latter included all the other ports served by the Capesize dry bulk carriers. The port and route network contains 52 main ports and 541 other ports as 79 port clusters, as shown in Figure 2. The main ports are mainly located in China (21 ports), Japan (7 ports), Australia (6 ports), and Brazil (4 ports) and accounted for 77.2% of the iron ore and coal trade in 2018, as seen from the worldwide port calling numbers of the Capesize dry bulk carriers [39], showing that most of the port calling numbers were served toward these countries. This tendency matches the worldwide top exporters and importers ranking of iron ore and coal in 2018 [31]. Additionally, this model also contains each port limitation and the sailing distance between the ports. The port limitations are set by port data and the largest possible ship visiting a port. Finally, the sailing distances were extracted from the IHS Fairplay Ports and Terminals Guide [43].

Cargo Movement Model
This model defines the worldwide port-to-port cargo movement demand served by the Capesize dry bulk carriers. This study defined the cargo movement demand between ports based on the estimated cargo movement (CARGO), which was calculated by considering the ship's DWT and draught when it arrived at the destination port [30]. The main ports are mainly located in China (21 ports), Japan (7 ports), Australia (6 ports), and Brazil (4 ports) and accounted for 77.2% of the iron ore and coal trade in 2018, as seen from the worldwide port calling numbers of the Capesize dry bulk carriers [39], showing that most of the port calling numbers were served toward these countries. This tendency matches the worldwide top exporters and importers ranking of iron ore and coal in 2018 [31]. Additionally, this model also contains each port limitation and the sailing distance between the ports. The port limitations are set by port data and the largest possible ship visiting a port. Finally, the sailing distances were extracted from the IHS Fairplay Ports and Terminals Guide [43].

Cargo Movement Model
This model defines the worldwide port-to-port cargo movement demand served by the Capesize dry bulk carriers. This study defined the cargo movement demand between ports based on the estimated cargo movement (CARGO), which was calculated by considering the ship's DWT and draught when it arrived at the destination port [30].
The set of available ships and routes of certain cargo types are represented as I and J, respectively. The cargo movement (CM ij ) per trip for each ship i serving the route j (t), and the annual cargo movement (CA j ) of the route j (t), are calculated as follows: where DWT i is the ship i DWT, d ij is the ship i arrival draught (m), and D i is the ship i maximum draught (m). The CA j for each cargo type was collected later in the OD tables. Hence, these tables contained the annual cargo movement in a port-level manner. This study classified the ship movement data into two voyage modes by referring to Arifin et al. [30]; laden and ballast. We applied Equation (1) to calculate the cargo movement (CM ij ) carried by ship i only for its laden voyage. Then, the transported annual cargo movement (CA j ) is calculated by using Equation (2). It summarized all the laden instances in the route j towards the target year. Finally, Table 2 shows a sample of the iron ore OD table. These data were used as the cargo movement demand in the ship allocation algorithm.

Ship Specifications
This model contains each ship and its attributes. A ship is individually attached to several static variables as follows:

Operation Condition Prediction
In addition to what has already been discussed, each ship is tied to its operation condition whenever a ship is assigned to a certain route. Therefore, it is necessary to know this information for all available routes. A ship's operation conditions throughout a particular route are assumed to be its typical performance, despite of the actual contract. This study defines the operation conditions as follows: • Draught rate: The ratio between the arrival draught and the ship's design draught at the destination port; • Average sailing speed: The ship's average sailing speed starting from its departure time at the origin port to its arrival time at the destination port; • Port staying time: The ship's port staying time starting from its arrival time until its departure time at the destination port. Figure 3 shows an overview of the proposed deep learning (DL) model. We constructed a DL model distinctively to predict the operation conditions in each voyage mode: laden and ballast. With each operation condition as the target feature, the model takes into account the following Capesize features as an input: • Ship movement data [39,43]: 36,596 entries of worldwide operations from January 2015 to May 2019, which consist of dates of arrival and departure, ports of origin and destination, and sailing distances; • Ship specifications [41]: 1658 entries, which consist of DWT, length, breadth, draught, design speed, main engine power, and built year; • Port data [42]: 564 entries, which consist of port limitations of length, breadth, and draught.
2015 to May 2019, which consist of dates of arrival and departure, ports of origin and destination, and sailing distances; • Ship specifications [41]: 1658 entries, which consist of DWT, length, breadth, draught design speed, main engine power, and built year; • Port data [42]: 564 entries, which consist of port limitations of length, breadth, and draught. The same architecture is used to predict the abovementioned operation conditions Its input features are represented as binary vectors by using the one-hot encoding. These are normalized using Min-Max scaling before the training process.
The proposed DL model is comprised of three hidden layers of 128 nodes equipped with Rectified Linear Unit activation functions. We applied the optimizer of the Roo Mean Square prop and the loss function of mean absolute error loss. Finally, we applied the holdout method to train our model, in which 70% was for the training data, 20% for the validation data, and 10% for the test data.
The average operation conditions of the worldwide routes (AVG) are the baselines that affirm the benefits of adopting the DL model. The DL prediction results, indicating a higher accuracy than the AVG, are shown in Table 3. Later, we will use the prediction results of the DL model to calculate the CARGO, GHG, and COST variables for each ship

GHG Calculation
In this study, we calculated the total GHG, which has been referred to in prior IMO GHG studies [28,29]. It comprised the GHG emissions from the main engine and auxiliary machinery (an auxiliary engine and a boiler). This calculation used ship data and port data [41,42], as well as the results of the global network model and ship model. First, we calcu lated the actual power at which the main engine operated on average in sailing mode. The actual power of the main engine ( ) of a ship operating in the route (kW) is calcu lated as follows [28,29,32,33]:  The same architecture is used to predict the abovementioned operation conditions. Its input features are represented as binary vectors by using the one-hot encoding. These are normalized using Min-Max scaling before the training process.
The proposed DL model is comprised of three hidden layers of 128 nodes equipped with Rectified Linear Unit activation functions. We applied the optimizer of the Root Mean Square prop and the loss function of mean absolute error loss. Finally, we applied the holdout method to train our model, in which 70% was for the training data, 20% for the validation data, and 10% for the test data.
The average operation conditions of the worldwide routes (AVG) are the baselines that affirm the benefits of adopting the DL model. The DL prediction results, indicating a higher accuracy than the AVG, are shown in Table 3. Later, we will use the prediction results of the DL model to calculate the CARGO, GHG, and COST variables for each ship.

GHG Calculation
In this study, we calculated the total GHG, which has been referred to in prior IMO GHG studies [28,29]. It comprised the GHG emissions from the main engine and auxiliary machinery (an auxiliary engine and a boiler). This calculation used ship data and port data [41,42], as well as the results of the global network model and ship model. First, we calculated the actual power at which the main engine operated on average in sailing mode. The actual power of the main engine (p ij ) of a ship i operating in the route j (kW) is calculated as follows [28,29,32,33]: where P i is the ship i maximum main engine power (kW), d ij is the ship i actual draught (m), D i is the ship i maximum draught (m), v ij is the ship i average sailing speed (kn), V i is the ship i design speed (kn), and ηw ij and η f ij are the correction factors of weather and fouling, respectively. The P i , D i , and V i of each ship were acquired from the ship data [41], d ij and v ij were generated by the ship model, and ηw ij and η f ij were considered as constants with values 0.909 and 0.917, respectively [28]. After calculating p ij , the main engine fuel consumption (A ij ) for a ship i serving the route j (t) is calculated as follows: where p ij is the actual power of the main engine (kW), s f c i is the ship i specific fuel oil consumption of the main engine (g/kWh), and st ij is the round-trip sailing time (d). In this study, the main engine of ship i was assumed to be a slow-speed diesel engine fueled by heavy fuel oil (HFO). Thus, s f c i for ships built before and after 2001 were assumed to be 185 and 175 g/kWh, respectively [28].
In addition to the main engine, we estimated the auxiliary machinery fuel consumption (an auxiliary engine and a boiler). The daily average fuel consumption of the auxiliary engine (a i ) and boiler (b i ) were determined by a linear correlation between the ship size and the average fuel consumption of the auxiliary machinery. The auxiliary machinery fuel consumption (B ij ) for a ship i operating on the route j (t) is calculated as follows: where a i is the ship i auxiliary engine daily average fuel consumption of the ship's auxiliary engine (t/d), b i is the ship i boiler daily average fuel consumption of the ship (t/d), and st ij is the round-trip sailing time (d). Next, the possible annual trip (n ij ) of a ship i serving the route j is defined as follows: where nx ij and ny ij are the possible annual trips considering the economic days and cargo movement, respectively. The yd is the number of days in a year, taken as 365 (d), md i is the assumed annual maintenance days (d), st ij is the round-trip sailing time (d), ps ij is the loading-unloading time (d), CD ij is the annual transported cargo movement of route j (t/y), and CM ij is the ship i cargo movement per trip serving the route j (t). Finally, we define the total GHG emissions, which covered the GHG emissions of the main engine and auxiliary machinery. We calculate the annual GHG emissions (GH ij ) for a ship i operating in the route j (t) as follows: where A ij is the main engine fuel consumption (t), B ij is the auxiliary machinery fuel consumption (t), CF and c f are the HFO and marine diesel oil (MDO) carbon factors (t-CO2/t-fuel), respectively, and n ij is the possible annual trip number. Following the IMO GHG study [28], CF and c f were taken as 3.114 and 3.206 t-CO2/t-fuel, respectively.

COST Calculation
In this study, the ship's total COST considers the fuel, operation, depreciation, and stockpile costs [1]. The fuel cost represents the cost of the fuel needed to transport specific cargo from the origin port to the destination port. The operation cost describes the cost of operating the ship. It comprises the crew, consumables, repairs, and insurance costs. The depreciation cost defines the reduction in the ship resale price after a year of operation. The stockpile cost depicts the cost of stockpiling transported cargo at the destination port until the cargo is consumed based on the port's daily cargo movement demand.

Fuel Cost Calculation
The fuel cost is calculated from the fuel consumption. The fuel cost (FC ij ) for a ship i serving the route j (USD) is calculated as the GHG emissions equation, by replacing the GHG emissions carbon factor (CF, c f ) with fuel price constants (FP, f p), as follows: where A ij is the main engine fuel consumption (t), B ij is the auxiliary machinery fuel consumption (t), and FP and f p are the HFO and MDO prices (USD/t-fuel), respectively. After observing the fluctuations in fuel prices [44], FP and f p were assumed to be 300 and 600 USD, respectively.

Operation Cost and Depreciation Cost Calculations
We composed response surfaces using the DWT and age of the dry bulk carrier ships to calculate the operation and depreciation costs. Considering Stopford [1], we expanded the correlation between the Capesize dry bulk carrier operation cost to other bulk carrier size classes [29,45,46]. Similarly, the depreciation cost was acquired by referring to the bulk carrier newbuilding and secondhand prices [44]. Figure 4a shows the daily operation cost for the given ship DWT and ship age which is defined by referring to the following attributes:  Correspondingly, the depreciation cost is determined by referring to the bulk carrier newbuilding and secondhand ship prices in 2019 [44]. The ship prices for the given ship DWT and ship age are shown in Figure 4b. The annual depreciation cost is obtained from the difference between the ship price in the current year and the previous year. Then, the  Correspondingly, the depreciation cost is determined by referring to the bulk carrier newbuilding and secondhand ship prices in 2019 [44]. The ship prices for the given ship DWT and ship age are shown in Figure 4b. The annual depreciation cost is obtained from the difference between the ship price in the current year and the previous year. Then, the daily depreciation cost was taken as annual depreciation cost divided by the number of days in a year, taken as 365 days. Finally, we represent these attributes into the bulk carrier operation cost and depreciation cost response surfaces, as shown in Figure 5. Correspondingly, the depreciation cost is determined by referring to the bulk carrier newbuilding and secondhand ship prices in 2019 [44]. The ship prices for the given ship DWT and ship age are shown in Figure 4b. The annual depreciation cost is obtained from the difference between the ship price in the current year and the previous year. Then, the daily depreciation cost was taken as annual depreciation cost divided by the number of days in a year, taken as 365 days. Finally, we represent these attributes into the bulk carrier operation cost and depreciation cost response surfaces, as shown in Figure 5. Thus, the operation and depreciation costs (OC ij , DC ij ) for a ship i operating in the route j (USD) are calculated using the following equation: where oc i is the ship i daily operation cost (USD/d), st ij is the round-trip sailing time (d), dc i is the ship i daily depreciation cost (USD/d), and yd is the number of days in a year (d), taken as 365. oc i and dc i for each ship i are shown in Figure 5a,b, respectively, based on each ship's DWT and age.

Stockpile Cost Calculation
Additionally, this study establishes the stockpile cost such that it resembles the cargo stockpile cost at the destination port. Origin and destination ports of route j are defined as k and l, respectively. The stockpile cost is determined by referring to the daily destination port cargo movement demand (CD l , t/d), which is calculated as the sum of cargo movements (∑ k (CA kl )) arriving at the destination port l divided by 365 days (yd). Then, the accumulated period (sd ij , d) is calculated by rounding off the division of cargo movement (CM ij ) carried by a ship i to the destination port l by the daily destination port cargo movement demand (CD l ) as follows: Hence, the stockpile cost (SC ij ) for a ship i serving the route j (USD) is calculated as follows: where sc l is the stockpile cost unit (USD/t/d), which depends on the port and cargo type. In this study, sc l was defined as a constant by referring to the average stockpile cost of the Japan Port Association [47]. Hence, sc l for iron ore and coal were taken as 0.012 and 0.042 USD/t/day, respectively.

Total COST Calculation
Finally, we define the total COST, which includes the fuel, operation, depreciation, and stockpile costs. The annual total COST (TC ij ) for a ship i operating in the route j (USD) is obtained as follows: where FC ij is the fuel cost (USD), OC ij is the operation cost (USD), SC ij is the stockpile cost (USD), DC ij is the depreciation cost (USD), and n ij is the possible annual trip number. A sample of the COST composition is shown in Figure 6.

Evaluation of GHG and COST Calculation
In this section, we will evaluate the GHG and COST calculation schemes. An evaluation of the GHG calculation is shown in Table 4. We calculated the main engine GHG emissions using the 2018 worldwide average ship age and draught rate of the Capesize dry bulk carriers [39] and the average variables depicted in the IMO GHG study [28]. For the following size categories, an error margin of less than 4% indicated the applicability of our GHG calculation scheme.  Figure 6. A sample of the COST composition.

Evaluation of GHG and COST Calculation
In this section, we will evaluate the GHG and COST calculation schemes. An evaluation of the GHG calculation is shown in Table 4. We calculated the main engine GHG emissions using the 2018 worldwide average ship age and draught rate of the Capesize dry bulk carriers [39] and the average variables depicted in the IMO GHG study [28]. For the following size categories, an error margin of less than 4% indicated the applicability of our GHG calculation scheme.  The error margin is between 0.7% and 4.1% for the COST item parameters, which justifies the practicability of our COST calculation scheme.

Algorithm 1: Ship Replacement while Preserving the Existing Ship Allocation
Our algorithm proposes an optimization by offering a direct clone of the existing ship, which is a new ship with the exact specifications and operation conditions that serve the same annual allocations. The existing ship is replaced if the new ship offers practically lower COST or GHG emissions. Locally, a new ship justifies the replacement of an existing ship with the capability to reduce COST or GHG emissions. Globally, a collection of existing ships is replaced by new ship formations.
For each instance of an existing ship i, we offered a new ship i with the same specifications, operation conditions, and annual allocations as the ship i but with the admiralty coefficient of the year 2021. Additionally, we proposed two optimization schemes: COSTand GHG-optimized. In this section, we have provided an example of a COST-optimized case. Along with the calculation of the annual COST (TC ij ) per ton cargo (CM ij ) of a ship i, we also calculated the annual COST (TC i j ) per ton cargo (CM i j ) of an offered ship i . Thus, the KPI (w ij , w i j ) for each ship i and i serving the route j is defined as follows: Finally, we acquired the global merit (W ij ), for an offered new ship i replacing an existing ship i by fulfilling the following conditions: where W ij represents a COST-effective ship. A ship i is determined to replace the ship i only if it offers a lower KPI. These steps resulted in a new ship allocation after conducting them for all the ships. Additionally, the seized COST reduction can be examined for an existing ship i replaced by a new ship i .

Algorithm 2: An Optimization for Reconstructing Ship Allocation
Based on the greedy algorithm approach, this algorithm reconstructs the ship allocation to transport the worldwide cargo movement demand using the available ships in the existing allocation or with the new ships instance. Using the total COST and GHG emissions for each ship and route, the allocation algorithm transports the same cargo amount as the existing ship allocation in a time-charter contract manner. This algorithm allocates ships conforming to the limitations of port dimensions for certain ships.
A ship with the highest global merit is sequentially designated to transport cargo on a certain route on a one-year time-charter contract by replicating the ship bidding process to obtain the most advantageous ship. The allocation algorithm was adopted separately for both cases to achieve a global reduction in the total COST and GHG emissions. Next, we have provided an example of the ship allocation in the COST-optimized case only, since the allocation scheme for both the COST-and GHG-optimized cases are the same.
The annual COST (TC ij ) per ton cargo (CM ij ) was set as an individual ship KPI (y ij ). Additionally, the average KPI for the route j (y j ) is known; therefore, the merit for each ship i operating in the route j (m ij ) is defined as follows: Finally, the global merit (M) is obtained as follows: A set of ships and routes with the global merit M is the most COST-effective allocation. Thereafter, the annual cargo movement demand of the route j (CA j ) is subtracted by the annual cargo carried by that ship (CM ij ). The allocation algorithm is then repeated until all the cargo movement demands are assigned to a particular ship, while the allocated ship (i) is removed from the available ship list (I). All these processes result in the operation-level ship allocation. The annual CARGO, total COST, and GHG emissions were elaborated for each ship i allotted on a specific route j. Therefore, we were able to explain the reduction in the total COST and GHG emissions and visualize the results geographically.

Case Studies and Discussions
We conducted several simulations using algorithms 1 and 2: ship replacement while preserving existing ship allocation (case study 1), optimization of ship allocation using existing ships (case study 2), and optimization of ship allocation using new ships instance (case study 3). First, we analyzed the actual ship allocation to understand the current ship allocation characteristics. Case study 1 analyzed new ships to replace existing ships without changing their allocation. Case study 2 reconstructed the ship allocation using only the existing ships. Case study 3 reconstructed the ship allocation using the existing ships and new ships instance (Table 6). Finally, we summarized the results of the average DWT and sailing speed graphically in a great circle format to give a better overview [48][49][50]. The following were the objectives of the simulations:

Actual Ship Allocation
Before conducting simulations, it is necessary to analyze the typical features of the actual ship allocation. Figure 7 proportionally illustrates the ships' average DWT and sailing speed, such that the line thickness indicates the annual cargo movement on each route, and the sorted colors represent the variance of the DWT and sailing speed.

Case Studies and Discussions
We conducted several simulations using algorithms 1 and 2: ship replacement while preserving existing ship allocation (case study 1), optimization of ship allocation using existing ships (case study 2), and optimization of ship allocation using new ships instance (case study 3). First, we analyzed the actual ship allocation to understand the current ship allocation characteristics. Case study 1 analyzed new ships to replace existing ships without changing their allocation. Case study 2 reconstructed the ship allocation using only the existing ships. Case study 3 reconstructed the ship allocation using the existing ships and new ships instance (Table 6). Finally, we summarized the results of the average DWT and sailing speed graphically in a great circle format to give a better overview [48][49][50]. The following were the objectives of the simulations:

Actual Ship Allocation
Before conducting simulations, it is necessary to analyze the typical features of the actual ship allocation. Figure 7 proportionally illustrates the ships' average DWT and sailing speed, such that the line thickness indicates the annual cargo movement on each route, and the sorted colors represent the variance of the DWT and sailing speed. In the West Australia-China routes, which have the highest cargo movement, the average DWT is 180,000-200,000, and the average sailing speed is 8.5-9.5 knots, which indicates the use of smaller ships at typically slower speeds. Additionally, the long-distance routes, such as Brazil-China, show the ships operating with an average DWT of 200,000-260,000, and an average sailing speed exceeding 9.5 knots, which indicates that the routes are served by larger ships at higher speeds.

Definition of Case Study 1
This case study attempted to directly replace each existing ship with a new ship without changing the actual ship allocation, and then to observe the changes in the total COST and GHG emissions, including COST-optimized ship replacement while preserving the existing ship allocation (case study 1a) and GHG-optimized ship replacement while preserving the existing ship allocation (case study 1b). The following were the parameters of case study 1: In the West Australia-China routes, which have the highest cargo movement, the average DWT is 180,000-200,000, and the average sailing speed is 8.5-9.5 knots, which indicates the use of smaller ships at typically slower speeds. Additionally, the longdistance routes, such as Brazil-China, show the ships operating with an average DWT of 200,000-260,000, and an average sailing speed exceeding 9.5 knots, which indicates that the routes are served by larger ships at higher speeds. This case study attempted to directly replace each existing ship with a new ship without changing the actual ship allocation, and then to observe the changes in the total COST and GHG emissions, including COST-optimized ship replacement while preserving the existing ship allocation (case study 1a) and GHG-optimized ship replacement while preserving the existing ship allocation (case study 1b). The following were the parameters of case study 1: • Method: Allocation algorithm 1; • Offered ships: A new ship for each existing ship in the actual ship allocation.
We offered new ships with the same principal particulars but with the admiralty coefficient of 2021.

Results of Case Study 1
The allocated new ship rates for cases 1a and 1b are shown in Figure 8. Both the cases show a typically high allocation rate for new ships. Figure 8a shows the major routes that the allotted new ships take for over 40%-60% of the operation numbers. Meanwhile, Figure 8b shows that the new ships accounted for more than 60% of the worldwide operation numbers.

Results of Case Study 1
The allocated new ship rates for cases 1a and 1b are shown in Figure 8. Both the cases show a typically high allocation rate for new ships. Figure 8a shows the major routes that the allotted new ships take for over 40%-60% of the operation numbers. Meanwhile, Figure 8b shows that the new ships accounted for more than 60% of the worldwide operation numbers. As shown in Table 7, minor reductions in the total COST and GHG emissions can be seen in case study 1. In case study 1a, new ships replaced 58.8% of the total number of existing ship operations, which reduced the total COST by 3.5%, and the GHG emissions were reduced by 5.6% compared to the actual ship allocation. Moreover, case study 1b allotted 87.5% of its operations to new ships, which allowed a 9.7% reduction in the GHG emissions, but the total COST increased by 0.4%.

Definition of Case Study 2
In this case study, we discussed optimizing the ship allocation using only the existing ships. The two case studies covered are: COST-optimized ship allocation using existing ships (case study 2a) and GHG-optimized ship allocation using existing ships (case study 2b). The following were the parameters for case study 2: • Method: Allocation algorithm 2; • Offered ships: The existing ships operated in the current ship allocation.

Results of Case Study 2
The results obtained for case study 2 are illustrated in Figure 9. From Figure 9a, we can ascertain that the Brazil-China routes allotted larger ships with a DWT of more than 260,000 and an average sailing speed of more than 9.5 knots. The average DWT of ships designated to the Brazil-China routes increased, compared to the actual ship allocation, due to the use of larger ships. As shown in Table 7, minor reductions in the total COST and GHG emissions can be seen in case study 1. In case study 1a, new ships replaced 58.8% of the total number of existing ship operations, which reduced the total COST by 3.5%, and the GHG emissions were reduced by 5.6% compared to the actual ship allocation. Moreover, case study 1b allotted 87.5% of its operations to new ships, which allowed a 9.7% reduction in the GHG emissions, but the total COST increased by 0.4%. Table 7. Results of case studies. In this case study, we discussed optimizing the ship allocation using only the existing ships. The two case studies covered are: COST-optimized ship allocation using existing ships (case study 2a) and GHG-optimized ship allocation using existing ships (case study 2b). The following were the parameters for case study 2:

Total Operation Numbers (Reduction in %)
• Method: Allocation algorithm 2; • Offered ships: The existing ships operated in the current ship allocation.

Results of Case Study 2
The results obtained for case study 2 are illustrated in Figure 9. From Figure 9a, we can ascertain that the Brazil-China routes allotted larger ships with a DWT of more than 260,000 and an average sailing speed of more than 9.5 knots. The average DWT of ships designated to the Brazil-China routes increased, compared to the actual ship allocation, due to the use of larger ships. In Figure 9b, the major routes are allocated to ships with an average sailing speed of 8.5-9.5 knots. Subsequently, Figure 9c depicts a similar average DWT as the actual ship allocation, and the average sailing speed shown in Figure 9d is generally slower. This analysis indicated the importance of speed reduction in optimizing the fuel consumption and GHG emissions of the ships.
As summarized in Table 7, case study 2a reduced the total COST by 7.3% compared to the actual ship allocation, whereas the total GHG emissions were reduced by 14.8% in case study 2b. Additionally, the ship operation numbers of case studies 2a and 2b were reduced by 1.1% and 2.2%, respectively. Such changes in the ship allocation were made possible by the operation-level alterations, despite simply employing the existing ships. Nevertheless, a significant reduction in the total COST and GHG emissions could not be obtained using only the existing ships.

Definition of Case Study 3
In addition to the analysis of case study 2 (Section 7.3), this case study presents a set of new ships. We examined the competitive new ships that could potentially replace the existing ships by conducting simulations to optimize the ship allocation with these new ships. The following were the parameters of case study 3: • Method: Allocation algorithm 2; • Offered ships: The existing ships of case study 2 and new ships.
The new ship specifications were sampled from the existing ships of the 2018 ship allocation. The set of sampled ships is defined as . We defined these ships as more efficient since they had the admiralty coefficient of the year 2021: 697.29 ( ), as shown in Figure 10. In Figure 9b, the major routes are allocated to ships with an average sailing speed of 8.5-9.5 knots. Subsequently, Figure 9c depicts a similar average DWT as the actual ship allocation, and the average sailing speed shown in Figure 9d is generally slower. This analysis indicated the importance of speed reduction in optimizing the fuel consumption and GHG emissions of the ships.
As summarized in Table 7, case study 2a reduced the total COST by 7.3% compared to the actual ship allocation, whereas the total GHG emissions were reduced by 14.8% in case study 2b. Additionally, the ship operation numbers of case studies 2a and 2b were reduced by 1.1% and 2.2%, respectively. Such changes in the ship allocation were made possible by the operation-level alterations, despite simply employing the existing ships. Nevertheless, a significant reduction in the total COST and GHG emissions could not be obtained using only the existing ships. In addition to the analysis of case study 2 (Section 7.3), this case study presents a set of new ships. We examined the competitive new ships that could potentially replace the existing ships by conducting simulations to optimize the ship allocation with these new ships. The following were the parameters of case study 3: • Method: Allocation algorithm 2; • Offered ships: The existing ships of case study 2 and new ships.
The new ship specifications were sampled from the existing ships of the 2018 ship allocation. The set of sampled ships is defined as S. We defined these ships as more efficient since they had the admiralty coefficient of the year 2021: 697.29 (ADM 2021 ), as shown in Figure 10.
Hence, the main engine power (P s ) was reduced by retaining a constant ship s design speed as follows [51]: where P s is the main engine power of ship s (kW) and ADM s is the admiralty coefficient of ship s. Table 6 lists the specifications of the new ships. These ships had the same principal particulars as the sampled ships but lower main engine power. As new ships could be allotted indefinitely, we discussed the following case studies: COST-optimized ship allocation with new ships instance (case study 3a) and GHG-optimized ship allocation with new ships instance (case study 3b).  [41,51,52].
Hence, the main engine power ( ) was reduced by retaining a constant ship design speed as follows [51]: where is the main engine power of ship (kW) and is the admiralty coefficient of ship . Table 6 lists the specifications of the new ships. These ships had the same principal particulars as the sampled ships but lower main engine power. As new ships could be allotted indefinitely, we discussed the following case studies: COST-optimized ship allocation with new ships instance (case study 3a) and GHG-optimized ship allocation with new ships instance (case study 3b). This case study optimized the ship allocation by allowing the existing ships to be replaced by new ships. This allowed us to ascertain the average DWT and sailing speed patterns simultaneously on a global scale. Figure 11 shows the ship allocation results of case study 3. From Figure 11a,c, we can observe that both the case studies allocated ships with a similar average DWT. The major routes displayed in Figure 11b have a uniform average sailing speed of more than 8.5 knots, similar to case study 2. The allocated ships' lower average sailing speed compared to the actual ship allocation can be seen in Figure  11d. The short-and long-distance routes allocated ships with an average sailing speed of less than 8.5 knots and 8.5-9.5 knots, respectively. This indicated that the importance of sailing speed varied depending on the sailing distance.

Results of Case Study 3
This case study optimized the ship allocation by allowing the existing ships to be replaced by new ships. This allowed us to ascertain the average DWT and sailing speed patterns simultaneously on a global scale. Figure 11 shows the ship allocation results of case study 3. From Figure 11a,c, we can observe that both the case studies allocated ships with a similar average DWT. The major routes displayed in Figure 11b have a uniform average sailing speed of more than 8.5 knots, similar to case study 2. The allocated ships' lower average sailing speed compared to the actual ship allocation can be seen in Figure 11d. The short-and long-distance routes allocated ships with an average sailing speed of less than 8.5 knots and 8.5-9.5 knots, respectively. This indicated that the importance of sailing speed varied depending on the sailing distance.
In Table 7, we observe that case study 3a reduces the total COST by 9.5% compared to the actual ship allocation. These reductions were accounted for by new ships, replacing 40.9% of the existing operations. In case study 3b, the GHG emissions were reduced by 22.8%, but the total COST increased by 4.6%. To conclude, we identified a significant reduction in the total COST and GHG emissions by adding new ships.

Discussion 1: Significance of the Total COST and GHG Emissions Reductions
The results of each case study are compiled in Table 7. The case studies 1a and 1b employed new ships covering 58.8% and 87.5%, respectively, of the worldwide operation numbers using allocation algorithm 1. However, only minor changes were typically observed compared to the actual ship allocation, with a 3.5% total COST reduction in case study 1a and a 9.7% GHG emissions reduction in case study 1b.
Case study 2 optimized the ship allocation using only the existing ships operating in the current ship allocation by deploying allocation algorithm 2. Case study 2a realized a 7.3% total COST reduction, two times that of case study 1a, compared to the actual ship allocation. For the GHG emissions aspect, case study 2b achieved a 14.8% reduction, which was 1.5 times that of case study 1b. A more significant reduction in the total COST and GHG emissions was achieved by case study 2 than case study 1.
Case study 3 proposed new ships in the application of allocation algorithm 2 under the same constraints as that of case study 2. The new ships were defined as being able to be allocated indefinitely, which allowed such ships with higher merits to replace the existing ships as needed. In case study 3a, 40.9% of the operation numbers were served by the new ships, allowing a 9.5% total COST reduction compared to the actual ship allocation. For GHG emissions, 48.2% of the operation numbers in case study 3b were served by the new ships. This enabled a considerable reduction of 22.8% in the GHG emissions. Additionally, case study 3 allocated new ships at a considerably lower rate than case study 1. Nevertheless, the most significant reduction in the total COST and GHG emissions were actualized by presenting the new ships for the ship allocation optimization. In Table 7, we observe that case study 3a reduces the total COST by 9.5% compared to the actual ship allocation. These reductions were accounted for by new ships, replacing 40.9% of the existing operations. In case study 3b, the GHG emissions were reduced by 22.8%, but the total COST increased by 4.6%. To conclude, we identified a significant reduction in the total COST and GHG emissions by adding new ships.

Discussion 1: Significance of the Total COST and GHG Emissions Reductions
The results of each case study are compiled in Table 7. The case studies 1a and 1b employed new ships covering 58.8% and 87.5%, respectively, of the worldwide operation numbers using allocation algorithm 1. However, only minor changes were typically observed compared to the actual ship allocation, with a 3.5% total COST reduction in case study 1a and a 9.7% GHG emissions reduction in case study 1b.

Discussion 2: New Ship Demand
In the previous section, case study 3 specified that notable results could be obtained by optimizing the ship allocation with the new ships instance. The allocated new ship rates of case studies 3a and 3b are shown in Figure 12. The allocated new ship rate was observed to vary depending on the sailing distance. Similarly, in the COST-and GHG-optimized cases, more than 60% of the new ships were allocated in the West Australia-China routes, the major iron ore routes with 646 million tons of cargo movement demands. Case study 2 optimized the ship allocation using only the existing ships operating in the current ship allocation by deploying allocation algorithm 2. Case study 2a realized a 7.3% total COST reduction, two times that of case study 1a, compared to the actual ship allocation. For the GHG emissions aspect, case study 2b achieved a 14.8% reduction, which was 1.5 times that of case study 1b. A more significant reduction in the total COST and GHG emissions was achieved by case study 2 than case study 1.
Case study 3 proposed new ships in the application of allocation algorithm 2 under the same constraints as that of case study 2. The new ships were defined as being able to be allocated indefinitely, which allowed such ships with higher merits to replace the existing ships as needed. In case study 3a, 40.9% of the operation numbers were served by the new ships, allowing a 9.5% total COST reduction compared to the actual ship allocation. For GHG emissions, 48.2% of the operation numbers in case study 3b were served by the new ships. This enabled a considerable reduction of 22.8% in the GHG emissions. Additionally, case study 3 allocated new ships at a considerably lower rate than case study 1. Nevertheless, the most significant reduction in the total COST and GHG emissions were actualized by presenting the new ships for the ship allocation optimization.

Discussion 2: New Ship Demand
In the previous section, case study 3 specified that notable results could be obtained by optimizing the ship allocation with the new ships instance. The allocated new ship rates of case studies 3a and 3b are shown in Figure 12. The allocated new ship rate was observed to vary depending on the sailing distance. Similarly, in the COST-and GHGoptimized cases, more than 60% of the new ships were allocated in the West Australia-China routes, the major iron ore routes with 646 million tons of cargo movement demands. Next, we examined the distinct specifications of the demanded new ships in the West Australia-China routes. Figure 13 shows the composition of operations in the COST-and GHG-optimized ship allocations with the new ships instance. The new ship E (DWT 290,000) was the most allocated, accounting for 47.9% and 54.7% of the operations of the COST-and GHG-optimized ship allocation, respectively.   Next, we examined the distinct specifications of the demanded new ships in the West Australia-China routes. Figure 13 shows the composition of operations in the COST-and GHG-optimized ship allocations with the new ships instance. The new ship E (DWT 290,000) was the most allocated, accounting for 47.9% and 54.7% of the operations of the COST-and GHG-optimized ship allocation, respectively.
This proved that the new ship E specification improvement, through the admiralty coefficient of 2021, presented a larger benefit than the existing ships in both the optimization cases. Moreover, looking at the significant main port restrictions of West Australia-China (Table 8), the largest new ship that entered was the new ship E, thus confirming that the demanded ship was a large iron ore bulk carrier. study 3a (COST-optimized) and (b) case study 3b (GHG-optimized).
Next, we examined the distinct specifications of the demanded new ships in the West Australia-China routes. Figure 13 shows the composition of operations in the COST-and GHG-optimized ship allocations with the new ships instance. The new ship E (DWT 290,000) was the most allocated, accounting for 47.9% and 54.7% of the operations of the COST-and GHG-optimized ship allocation, respectively.

Concluding Remarks
In this study, we have proposed an enhancement to our basic ship-planning support system. We expanded the scope of our simulations to a worldwide scale, which allowed a broad comprehension of the demand ship specifications. In addition, we developed a ship allocation algorithm to consider the detailed COST and GHG emissions attributes. The system itself consisted of three distinct models and an algorithm to resemble the ship allocation bidding process: a global network model, cargo movement model, ship model, and ship allocation algorithm. The global network model defined the main ports and port clusters, which resulted in the following: a cluster of ports and their attributes, port and route limitations, and sailing distances. The cargo movement model specified the cargo movement and cargo types between the ports. The ship model represented each ship and its specifications. This model predicted the operating conditions: draught rate, average sailing speed, and port staying time for each ship serving a particular route.
Ultimately, as the core of our system, the ship allocation algorithm reconstructed the operation-level ship allocation. It calculated COST and GHG emissions for each ship and a greedy algorithm mechanism was used to resemble the ship allocation. This study proposed two algorithms: ship replacement while preserving the existing ship allocation and an optimization for reconstructing the ship allocation. Allocation algorithm 1 offered a new ship with enhanced performance for each existing ship without changing the current ship allocation. Additionally, we reallocated the ship by applying allocation algorithm 2, which optimized the total COST and GHG emissions. We conducted two discrete simulations: optimization using only the existing ships operating in the actual ship allocation and with new ships instance. The ship allocation optimization indicated that significant reductions in the total COST and GHG emissions were not achievable using only the existing ships. Using the developed system, we were able to recreate an operation-level ship allocation considering various scenarios. Finally, we confirmed the demanded new ship specifications by presenting the new ships instance. However, this study is only feasible for the current constraints of the Capesize dry bulk carrier operating in the time-charter contract, despite the presented results. The application of the proposed system for other dry bulk carrier size categories such as Panamax and Handymax, and the voyage-charter contracts are considered future tasks. Likewise, considering the importance of weather routing, the weather correction factor in the main engine actual power calculation, which was previously assumed constant, is possibly varied [53]. Hence, further studies are crucial to assess the applicability of our system beyond these limitations.