# Bottom-Up Approach Ship Emission Inventory in Port of Incheon Based on VTS Data

Authors: Hyangsook Lee; Hoang T. Pham; Maowei Chen; Sangho Choo
Journal: Journal of Advanced Transportation
 (2021)

---

## Abstract

As a result of the rapid growth of international trade, atmospheric pollution from transportation has been more topical than ever, especially in dense hub port-cities. The shipping industry should pay more attention corresponding to its contribution to local atmospheric pollution. This paper supports the application of data collected from the vessel tracking service system with a bottom-up approach to generate a comprehensive 2019 local ship emission inventory at Port of Incheon. The calculated emission inventory presented the dominance of CO2 emission and the considerable contribution of NOx and SOx emissions, the significant contribution of auxiliary engines during the hotelling at berth during the year of 2019. Then, based on calculated emission inventory, this study suggested and simulated applicable green policies in the practice: (1) local emission control area realization, (2) vessel speed reduction program, (3) application of cold ironing, and (4) establishment of a national integrated emission platform. The combination of the three first policies could help reduce the significant volume of emitted CO (29%), NOx (30%), SOx (93%), PM10 and PM2.5 (64%), VOC (28%), NH3 (30%), and CO2 (30%).


---

## Body

## 1. Introduction

Transportation has been considered a remarkable contributor to atmospheric pollution [1]. The recent decades have detected a growing consideration of a leading global anthropogenic emissions contribution from in-port traffic [2]. With 90% of the global transport volume being carried by ship, maritime air pollution has exacerbated local society’s health because of the high volume of emissions exhausted as well as the geographical condition. Approximately 70% of ship emissions are discharged within a radius of 400 km from the coastline [2]. As the result of the hub-port trend, ships often concentrate on well-located ports, where they often are situated close to busy industrial zones and densely populated cities. In-port ship emissions could supply 55–77% of total local emissions [3, 4]. Therefore, these port-cities and their hinterlands are grappling with a great burden from transport emissions and infrastructure [5]. It is estimated that approximately 230 million people, living in the top 100 world ports, are directly affected by exhaust shipping emissions [6]. Thus, the port industry has kept the momentum in green reform by reducing ship emissions and performs more of their social responsibility under higher pressure from public authorities [7, 8].More much than it seems, it is stated that there are 450 different air pollutants which are emitted by ship engine combustion [9]. Greenhouse gases (GHGs), carbon monoxide (CO), nitrogen oxides (NOx), sulphur oxides (SOx), and particulate matter (PM) are stated as key ship-source air pollutants [10, 11]. Shipping SOx and NOx contribute markedly roles compared to other air pollutants in total national emissions [12]. Besides a negative affection to global climate, numerous research articles have proved a considerable negative tie between closed-to-land ship emissions and the local community’s health (lead to asthma, respiratory and cardiovascular diseases, lung cancer, and premature death) [13]. Exhausted PM emissions from shipping activities are considered the main reason for the annual 60,000 cardiopulmonary and lung cancer deaths in Europe, East Asia, and South Asia littoral regions [14].To reduce air pollution from ships, emission control legislations were contemplated by international regulation-making agencies. Likewise, the International Maritime Organization (IMO) is reducing ship emission through its International Convention on the Prevention of Pollution from Ships (MARPOL) Annex VI. The sulphur content in the fuel used onboard is reduced from 3.5% m/m to 0.5% m/m from 1 January 2020 [15]. The sulphur directive in the European Union (EU), the fourth federal standard, also called Tier 4, in the United States (U.S.), and emission control areas (ECAs) in China are promulgated to follow this regulation [16–18]. It is reported that the implementation of these regulations helped reduce 77% SOx, 23% NOx, and 27% PM10 emissions in the EU by 2017, compared with 2000 [19]. Despite the expected expansion in shipping, stricter emissions regulations could reduce the volume of SOx, PM2.5, and NOx emissions by 87%, 92%, and 56%, respectively, between 2015 and 2050 [20].In Korea, the public has been enhancing national air quality through its 10-year comprehensive plans since the 1990s. Air quality management basic plan during 2015–2024 is established to manage an integrated air management system. Besides, in 2017, an other comprehensive plan of particular matter management was announced with an epicenter of 30% PM emission reduction and reinforcement control of the volatile organic compounds (VOCs) until 2022. [21] From December 2019, a vessel speed reduction (VSR) program has been conducted to reduce the PM level during winter in five key ports of Korea, including Port of Incheon (POI). A VSR area in each port spans 20 nm in radius [22]. However, almost all national policies and green measures are decided with the support of annual national emission inventories provided in the clean air policy support system (CAPSS), using national fuel supply statistics. However, these inventories revealed inconsistency and uncertainty with other local academic research articles, because of no actual traffic data [23]. Quantifying emitted emissions from shipping activities is a stride to explore the cause and scale of marine pollution [24]. It also provides valuable statistics to anticipate the future trend of pollution and then establishes counter-measures and policies [25]. Therefore, generating a more reliable in-port ship emission inventory is a pressing issue for producing appropriate regional policies and measures related to air pollutant management at ports.In this study, a comprehensive annual ship emission inventory of target pollutants of CAPSS at POI in 2019 was estimated, according to geographical areas, operational phases, and ship types. This study applied another ship tracking system named vessel tracking service (VTS), operated by port authorities or coast guards to keep track of ship traffic and ensure navigational safety in the port area. This emission inventory applied a bottom-up approach with both EO and FC methods to see the difference in the estimating results of the two methods and then achieved a comprehensive multidimensional understanding of air pollutants from in-port ship operation. Then, applicable green policies for the POI situation in the practice will be suggested and simulated to evaluate their effects based on the calculated emission inventory. The remainder of this study is grouped as follows: Section2 summarizes the previous studies related to ship-related emissions and then points out the research demand of this paper. Materials preparation is introduced in Section 3. Section 4 describes two applied methodologies for ship emission prediction. Section 5 presents and discusses calculation results. Several policies to reduce the volume of emissions emitted from ships operated at POI are suggested and simulated in Section 6 and Section 7 is the conclusion of this study.


## 2. Literature Review

Numerous previous studies have suggested estimation methods to generate ship emissions. The mainstream methods to prepare ship emission inventory can be broadly categorized into two main approaches: top-down approach (fuel-based) and bottom-up approach (activity-based). In the top-down approach, exhausted ship emissions are estimated from analyzed statistically marine fuel sales and fuel-related emission factors [26]. The marine fuel sale reports are mainly published by the Energy Information Administration, the International Energy Agency, and the United Nations Framework Convention on Climate Change [27]. This approach is recommended for a low level of traffic data availability situation. Due to the positive correlation between fuel consumption and emissions discharged from engine combustion, this approach would be the most accurate method if researchers could be confident about the marine fuel sales data collected. However, it proved that there is a significant discrepancy between banker fuel sales statistics and the actual fuel used by the global fleets [28]. The main reasons behind this difference could be the practice of offshore tinkering or the inaccurate fuel statistics in several countries [27].Due to difficulties in collecting reliable inputs for the top-down approach, the bottom-up approach is recommended as another choice if accurate sailing statistics (e.g., actual travel distance with speed and port calling records with real-time operations) are available. Also, this method requires a higher level of input parameters such as detailed ship technical data (e.g., ship types, engine characteristics, and design information). An amount of emitted emissions is calculated for each specific ship activity and then scaled up over activities and trips to figure out the total volume of emissions [29]. The generic emission estimation equation of bottom-up can be expressed in the following equation:(1)E=Energy.EF.CF,where E is the emission amount of certain pollutants from the ship’s engines; Energy is the energy demand; EF is the emission factor; and CF is a controlling factor used when the ship is equipped with reduction technologies or fuel correction factors. In general, the estimation of energy demand can be dichotomized into the following:(a)
The total energy output (EO) of operated engines during the operating time, by applying the following equation:(2)Energy=P.LF.T.(b)
The total engines’ fuel consumption (FC) during the operating time, using the following equation:(3)Energy=SFOC.LF.T,where P denotes engine power; LF denotes load factor; T denotes activity time; and SFOC denotes specific fuel consumption.Bottom-up approach ship-related emission studies since 2009 are summarized in Table1. As shown in Table 1, the EO method dominated the bottom-up approach for estimating ship-related emissions, especially from 2014. The most considered pollutants in previous studies are CO2, NOx, CO, SOx, and PM.Table 1
Summary of bottom-up approach studies.
AuthorMethodOperation modesConsidered pollutantsGeneral estimating equationKwon et al. [1]EOCru, Man, HotNOx, CO, SOx, PM10E = P. LF. T. EFSong [2]EOCru, Man, HotNOx, CO, SOx, CO2, CH4, N2O, HC, PM10, PM2.5E = P. LF. T. EF. CFSanabra et al. [5]EOCru, Man, HotSO2, PM2.5, NOx, VOCsE = P. LF. T. EFMerk [6]EOCru, Man, HotCH4, CO, CO2, NOx, PM10, PM2.5, SOxE = P. LF. T. EFEEA [12]EOCru, Man, HotNOx, CO, SOx, NMVOC, TSP, PM10, PM2.5,…E = P. LF. T. EFSong and Shon [23]EOCru, Man, HotSO2, NOx, CO2, VOC, PME = P. LF. T. EFICF International [30]EOCru, Dec, Man, HotNOx, CO, SOx, CO2, HC, CH4, PM10, PM2.5, N2OE = P. LF. T. EFCorbett et al. [31]FCN/aCO2E = SFOC. LF. T. EFJoseph et al. [32]EOCru, Man, HotNOx, SO2, TSP, PM10E = P. LF. T. EFDeniz and Kilic [33]FCCru, Man, HotNOx, SO2, CO, CO2, VOC, PME = SFOC. T. EFDeniz et al. [34]EOCru, Man, HotNOx, SO2, CO2, HC, PME = P. LF. T. EFHowitt et al. [35]FCN/aCO2E = P. LF. SFOC. T. EFKiliç and Deniz [36]EOCru, Man, HotNOx, SO2, CO2, HC, PME = P. T. EFLonati et al. [37]FCMan, HotNOx, SOx, CO, VOC, PM10E = SFOC. T. EFTzannatos [38]EOMan, HotNOx, SO2, PME = P. LF. T. EFShin and Cheong [39]FCMan, HotCO2, N2O, CH4Man: E = D. FE-1. EFHot: E = (FC. 0.2). 0.79. EFVillalba and Gemechu [40]EOMan, HotCO2E = P. LF. T. EFChang and Wang [41]FCDecNOx, SO2, CO2, HC, PME = SFOC. LF. T. EFBerechman and Tseng [42]EOHotNOx, CO, SO2, CO2, HC, VOC, PM10, PM2.5E = P. LF. T. EFYau et al. [43]EOCru, Dec, Man, HotNOx, SO2, PM10E = P. LF. T. EFChang et al. [44]FCCru, Man, HotCO2E = SFC. LF. T. EFMcArthur and Osland [45]EOCru, Dec, Man, HotNOx, SO2, PM10E = P. LF. T. EFNg et al. [46]EOCru, Dec, Man, HotSO2, NOx, CO, VOC, PM10E = P. LF. T. EFSaraçoğlu et al. [47]EOCru, Man, HotNOx, SO2, CO2, HC, PME = P. LF. T. EFTai and Lin [48]FCCru, Man, HotNOx, SO2, CO2, HC, PME = T. FE. EFChang et al. [49]FCCru, Man, HotNOx, SO2, PME = SFC. LF. T. EFLiu et al. [50]FC, EOMan, HotSO2E = FC. EFE = P. LF. T. EF. CFGoldsworthy and Goldsworthy [51]EOCru, Man, HotNOx, CO, SO2, CO2e, PAH, VOC, PM10, PM2.5E = P. LF. T. EFCoello et al. [52]EON/aCO2, NOx, CO, NMVOC, SOx, PME = P. LF. T. EFTichavska and Tovar [53]EOCru, Man, HotCO2, NOx, CO, SOx, PM2.5E = P. LF. T. EFMaragkogianni et al. [54]EOMan, HotNOx, SOx, PME = P. LF. T. EFCullinane et al. [55]EOMan, HotNOx, CO, SO2, CO2, HC, PM10, PM2.5E = P. LF. T. EFFan et al. [56]EON/aNOx, CO, SO2, NMVOC, PM10, PM2.5, OC, EC, V, NiE = P. LF. T. EF. CFPapaefthimiou et al. [57]EOMan, HotNOx, SO2, PM2.5E = P. LF. T. EFChen et al. [58]EOCru, Man, HotNOx, CO, SO2, HC, PM10, PM2.5E = P. LF. T. EF. CFStyhre et al. [59]EOCru, Dec, Man, HotCO2eE = P. LF. T. EFAlver et al. [7]EOCru, Man, HotNOx, SO2, HC, PM10E = P. LF. T. EFKnežević et al. [60]EOMan, HotNOx, SOx, VOC, PME = P. LF. T. EFKhan et al. [61]FCCru, Man, HotCO2E = SFC. LF. T. EFSun et al. [62]EOCru, Dec, Man, HotNOx, CO, HC, CO2E = P. LF. T. EFCao et al. [63]EOCru, Man, HotNOx, CO, SO2, CO2, HC, PM10, PM2.5E = P. LF. T. EF. CFZhang et al. [64]EOCru, Dec + Man, HotNOx, CO, SOx, CO2, HC, PM10, PM2.5E = P. LF. T. EF. CFIvce et al. [65]EOCru, Man, HotCO2E = P. LF. T. EFWan et al. [66]EOCru, Dec, Man, HotSOx, NOx, PM10E = P. LF. T. EF. CFStazić et al. [67]EOCru, Man, HotSOx, CO2, VOC, PME = P. LF. T. EFWan et al. [68]EOCru, Dec, Man, HotNOx, CO, SOx, CO2, HC, CH4, NMVOC, PM10, PM2.5E = P. LF. T. EF. CFEkmekçioğlu et al. [69]EOCru, Man, HotCO, CO2, NOx, SO2, PM, VOCE = P. LF. T. EFLee et al. [70]EOCru, Man, HotCO, NOx, SOx, PM, VOC, NH3E = P. LF. T. EFIt is widely agreed that the bottom-up approach is generally more accurate than the former [52]. Therefore, over the last two decades, the latter has been adopted frequently by researchers [26]. Inaccuracy estimation from the top-down approach at the national level could lead to a completely different understanding and measures at the regional air quality management level. Therefore, at the regional level, the bottom-up approach is recommended to be a reliable source to establish effective green policies for better air condition [40, 53, 71, 72]. However, for larger scales, applications of the bottom-up approach have been finite because of data gaps and anomalies [52, 54]. Also, the diversity and details in input information could help improve the reliability of results; however, its corresponding complexity also increases the level of uncertainty, which is associated with the different applied models and assumptions, and how to collect the set of inputs [27, 29]. Especially, on a global scale, the use of average input parameters such as different load factors, operational activity time, and emission factors for a different size, age, and type of ships leads to considerable uncertainty in estimation [54, 73].Recently, a vessel tracking system named automatic identification system (AIS), proposed in the IMO International Convention of Safety of Life at Sea (SOLAS) for better identification of ships, has been widely applied in the stream to enhance estimations [44, 46, 49, 51, 52, 56, 74] and considered as “the best method to report the activities and movements of ships” [26]. This system is equipped compulsively on passenger ships of all sizes and commercial ships with 300 gross tonnage (GT) or more operated at sea. Therefore, high-resolution ship movement AIS data could be a fount of reliable relative ship operational profile as vessel travel time and average speed between waypoints on the sea at short-time intervals and identify ship routes. Also, AIS collects and transmits ship characteristics as IMO identification number, size, weight, etc. [52]. Hence, ship’s activities are better geographically characterized and analyzed, and consequently, it improves the reliability of emission inventory investigation. Despite the abovementioned benefits, the AIS system does not cover ships with less than 300 GT. Besides, with ships greater than or equals to 300 GT, in several cases, AIS data could not be fully accessible. Time gaps also occur in several cases if temporary signal disruption happens, resulting in erroneous position reports.During the last decades, several investigations about Korean port-related ship emissions have been conducted domestically. At the national level, under CAPSS, national air pollutants emission inventories, including in-port ship emissions, had been updated annually, and the latest update is for the year 2016 [75]. However, CAPSS only considers hotelling at anchorage and berth and maneuvering process in port emission. The equations with values of parameters, provided by NIER [76], follow a top-down approach based on fuel consumption and are shown in the following equations:(4)E=∑phFCph×EF,(5)FChotelling=N×DF×0.79×0.2,(6)FCmanoeuvring=∑2×N×DFE,where FC represents fuel consumption; N represents annual total ship call; DF is average daily fuel consumption by ship type (ton/day); D represents average travel distance in port (km); FE represents average fuel economy by ship tonnage (km/kL); 0.79 is the average hotel time of a ship call per day; and 0.2 means the fuel consumption of hotelling at anchorage and berth will be assumed as 20% of the fuel consumption of full operation.At the regional level, Shin and Cheong [39] generated the GHGs ship emissions at Port of Busan (POB) with a top-down approach. Chang et al. [44] firstly used the bottom-up approach to assess the GHGs emissions from ship operations at the Port of Incheon (POI) from January to October 2012 based on the FC method. Later, Chang et al. [49] continued to apply the FC method to update NOx, SOx, and PM emissions released from ships at POI during the same period. With the same method, Khan et al. [61] collected ship traffic data at POI from the AIS system to figure out ship GHG emissions, however, during only October 2014. Kwon et al. [1] developed a system for emission estimation and then considered the amount of CO, NOx, SOx, and PM10 emitted from ships around POI only in December 2017. Lee et al. [70] estimated a non-CO2 emission inventory in POI during 2017 using VTS data and the EO method.This paper supports the idea of the application of VTS data in the bottom-up approach studies as a promising confident alternative suggested by Lee et al. [70]. However, the contribution of this paper is that(1)
On one face, this paper considers an FC method, which was used to estimate emissions in POI several times before by Chang et al. [44], Chang et al. [49], and Khan et al. [61]. However, all of these three papers did not provide an up-to-date and full-year estimation. On the other hand, this paper also considers the EO method applied by Lee et al. [70]; however, this paper investigates additionally CO2 emission, which is the most important pollutant but was not considered by Lee et al. [70]. Then, the paper compares the results from the two methods to figure out the gap between them. Hence, it provides a comprehensive annual ship emission inventory by different calculating methods, which covers all ships operated at the port during all months in one calendar year as well as focuses on all key pollutants of CAPSS [77]: CO2, CO, NOx, SOx, PM (including PM10 and PM2.5), VOCs, and ammonia (NH3);(2)
This paper suggests several available policy implications to reduce the local ship-related emissions and then simulated the effects of them if they will be conducted in the future.


## 3. Materials Preparation

### 3.1. Geographical Scope

Incheon metropolitan city is South Korea’s third most populous city after Seoul and Busan and also is the gateway of the northwestern area and Seoul capital area, the world’s fourth-largest metropolitan area by population, with the Yellow Sea. After becoming an international port in 1883, POI has been the country’s logistics hub that handles and distributes both general cargoes and containers. It is ranked as the global 27th busiest port in terms of cargo volume and the 50th biggest container port [78]. Currently, POI has expanded into 125 berths, with a total of 26,031 meters of berth length, distributed to five main component ports as North Port, Inner Port, Coastal Port, South Port, and New Port. The North Port is an industrial port specialized in managing raw materials as scrap iron, feed-by products, hardwood, etc. With a lock-gate, a calm water level in the Inner Port is available for handling semiconductor equipment, automobiles, and precision machine parts. It also is a multipurpose port that serves other cargoes as grain, fruit, and general cargo. The South Port is available for handling small and medium containers and general cargoes. The New Port is under-constructed and specialized for handling medium and large containers. The Coastal Port, with international passenger and ferry terminals, focuses on serving passenger cruise ships and car ferries. Also, around these component ports, POI operates three smaller specialized ports named Geocheom-do Port, Song-do Port, and Yeongheung-do Port. The Song-do Port is specialized in fossil energy products, while the Geocheom-do Port is used for handling sand. The Yeongheung-do Port is constructed to support Yeongheung Thermal Power Plants. The capacity of POI is described in Table 2.Table 2
Summary of POI’s capacity. Source: [79].
No.PortHandling CapacityMain productsMax. Ship DWTBerthsBerth length (m)Bulk (1000RT)Container (1000TEU)1Inner Port50,000469,83838,161General cargo, iron, grains2South Port100,000253,64217,600762Chemical, cement, sand3Coastal Port50,00091,429Passenger, oil, LPG4North Port100,000266,42113,900Oil, general cargo, wood product5New Port3,00061,6002,100Container6Song-do75,00041,300LPG, oil7Yeongheung-do200,00051,12614,690Bituminous coal, limestone8Geocheom-do5,00046758,320SandSummary12526,03192,6712,86228,735.5The POI geographical segments are presented in Figure1. The study not only covered ship operation inside the port boundary but also the “affected zone” (within 5 km from the port boundary) to consider the effect of close-to-port emitted ship emissions, following EPA [80].Figure 1
Geographical segments of Port of Incheon. Note: images were downloaded from the Internet and Incheon Port Authority (IPA) website.


### 3.2. Activity Phase

A typical ship call often contains consecutive activity phases, which consist of a series of continuous activities that have similar features, categorized as “cruising, “maneuvering,” “hotelling at anchorage,” and “hotelling at berth.” Due to the VSR program was applied for a few last days in December 2019, the “decelerating phase” is assumed as not considered in this study. As defined in ICF [30] and by Song [2], at the “cruising” phase, the ship moves inside the port boundary, and all engines keep running. “Maneuvering” is the time that the ship transits between the breakwater (intersection of open sea and inland waterway) and berths at a slower speed. Even with tug assist, the propulsion engines are still in operation. “Hotelling at anchorage” is a period of waiting for berth call while “hotelling at berth” is the time that ships are docking at berth. In both “hotelling at anchorage” and “hotelling at berth” phases, only auxiliary engines (including boilers) still work with a peak load for providing on-board power. In the “hotelling at berth” phase, the ship can use shore power instead of turning on auxiliary engines, called cold ironing; however, cold ironing has been still not implemented in POI yet. The average speed and travel distance in each phase were provided by a local pilot company. Table 3 summarizes detailed activity phases information in POI.Table 3
Ship activity segments and features in POI.
Phase categoryPropulsion engineAuxiliary engineBoilerAvg. speed (knots)Travel distance (nm)1Hotelling at anchorageOffOnOn002CruisingOnOnOff12Varies by berth3ManeuveringOnOnOnAround 3.514Hotelling at berthOffOnOn00


### 3.3. Ship Classification

Ship characteristics (e.g., speed, engine size, and usage) vary considerably depending on the ship type. Previous studies classified ships in different ways; however, the common standard applied is a type of cargoes carried on the ship. This study followed ICF’s classification with 11 types of ships [30]. However, in POI, there are eight types available, involving bulk carrier, container ship, cruise ship, general cargo, reefer, roll-on/roll-off (RORO), tanker, and miscellaneous.


### 3.4. Data Collection and Analysis

#### 3.4.1. Data Source

Due to the limitations of the AIS system mentioned above, in this study, in-port ship traffic data were collected from the VTS-based Korean Port Management Information System (Port-MIS) to consider ships with GT less than 300. A total of 16,677 ship calls were recorded from Port-MIS, as shown in Table4. Tanker fleet dominated ship calls at POI in 2019 with 42.6% of the total ship calls, following by the general cargo fleet and container ship fleet with 20.2% and 17.3%, respectively. Reefer fleet only contributed negligibly to total ship calls with 0.02% of total ship calls.Table 4
Ship call statistics (unit: ship call). Source: Port-MIS.
Ship typesNorth PortInner PortCoastal PortSouth PortNew PortOthersTotalBulk carrier253304—104——661Container ship271379681,871—2,886Cruise ship6367—123General cargo1,0566961491,40647203,374RoRo2777738330——1,847Reefer2—11——4Tanker4,1444736398361658477,104Miscellaneous113617720429294778Total5,6032,3271,6173,8562,1121,16216,677


#### 3.4.2. Data Cleaning and Make-Up

A total of 47,626 statuses describing the entire marine traffic at POI during 2019 were collected. These statuses were noted in order of occurrence inside POI, which is different from the ship call’s timeline format. Two consecutive events in the same ship call usually are several hours apart; then, a ship call can last from a few hours to a few days. Also, a high level of traffic exaggerated the complexity level of the original data and made it an arduous task to handle nagging details. Therefore, a code, combined from call sign and time points, was applied to resort original data corresponding to the timeline of each ship call.Waypoints are used to present ship activities in data. A typical ship call contains waypoints as (1) entrance, clearance from the port area; (2) anchoring, collecting anchor; (3) berthing, unberthing; and (4) shift within the port area between berths (if need). A combination of two consecutive waypoints identifies a corresponding activity phase of the ship call. The time information of each waypoint is provided in the YYYY/MM/DD HH: MM format. The time of an activity phase is identified as the subtractions of 2 corresponding waypoints’ time. The average moving time and distance of the cruising and maneuvering process for each target port were calculated. In several cases, when only arrival and departure times were reported in the data, the total time of ship spent on the port is the difference between arrival time and departure time. Later, the corresponding average moving time and distance for the target port are applied, and then the hotelling time at the berth of that ship is the subtraction of average cruising and maneuvering time from the total time.


#### 3.4.3. Ship Basic Data

Basic information about the ship, involving name, ship type, engine type, propulsion engine power, weight tonnage, and design speed, was obtained from the Korea Ship Safety Technology Authority and Korean Register. For missing propulsion engine power data, a simple linear approach to modeling the relationship between ship weight tonnage and propulsion engine power is applied by ship types. Linear models are shown in Figure2. The coefficients of determination (R2) are relatively high (>0.85), showing the high reliability of prediction.Figure 2
Linear regression of missing main engine data by ship type.Besides, installed auxiliary engine power also was not collected completely because there is a lack of information from manufacturers. However, it is impossible to apply a linear approach to estimate auxiliary engine power because there is no connection between installed auxiliary engine power and speed [51]. Also, it is difficult to recognize the using level of auxiliary engines in each activity phase during a ship call in the practice. Therefore, in previous studies, installed auxiliary engine power is commonly estimated by applying the default ratios by ship types with the total propulsion engine power of the ship [51]. This study applied the ratios revised from ICF [30], which is shown in Table 5.Table 5
The ratio for estimating auxiliary engine power. Source: ICF [30].
Ship typeAuxiliary to propulsion ratioBulk carrier0.222Container ship0.220Cruise ship0.278General cargo0.191RORO0.259Reefer0.406Tanker0.211Others0.100Lastly, the real values of revolutions per minute (RPM) of engines also were obtained from the Korea Ship Safety Technology Authority and Korean Register. RPM value of the engine helps to classify those engine speed designations: high-speed diesel (HSD), medium-speed diesel (MSD), and slow-speed diesel (SSD). Then, based on engine speed designations, proper EF values are matched in the calculation. In the case of missing values, they were made up of average values, as shown in Table6.Table 6
Average RPM values.
Ship typeAverage RPM values of propulsion enginesAverage RPM values of auxiliary enginesBulk carrier103843Container ship131847Cruise ship484895General cargo2711,225RoRo105787Reefer3111,043Tanker2861,147Miscellaneous4901,037


#### 3.4.4. Activity Time

The ship transit time was calculated with the actual speed of the ship and the representative travel distance of each port (from border to berth). The average transit time by the port is shown in Table7. Among ports in POI, Inner Port shows the largest value because of a lock-gate.Table 7
Average transit time in POI 2019 (unit: hour).
PortAverage transit timeNorth Port2.2Inner Port3.0Coastal Port1.6South Port1.7New Port1.4Song-do Port1.4Geocheom-do Port2.5Yeongheung-do Port0.75The actual times in hotelling at anchorage and berth were obtained from the difference in time between anchoring, collecting anchor, and berthing, unberthing. To avoid anomalies, the ceilings are applied as 7 days (168 hours) for hotelling at anchorage and 14 days (336 hours) for hotelling at berth. The average docking time by ship type after applying ceilings is shown in Table8. Almost average values by ship type (except container ship’s average time) are considerably greater than the national average time applied in CAPSS (0.79 days∼18.96 hours). These differences may lead to considerable variances between estimated results and national estimations, especially tankers, general cargo ships, and container ships dominated the total number of ship calls at POI with the north of 80%.Table 8
Average time in hotelling at anchorage and berth in POI 2019 (unit: hour).
Types of shipAverage time in hotelling (hour)Compare to the national avg. time in CAPSS (%)Bulk carrier83.3439.3Container ship13.671.7Cruise Ship50.3265.3General cargo39.9210.4RORO25.9136.6Reefer89.0469.4Tanker22.0116.0Miscellaneous51.9273.7


## 3.1. Geographical Scope

Incheon metropolitan city is South Korea’s third most populous city after Seoul and Busan and also is the gateway of the northwestern area and Seoul capital area, the world’s fourth-largest metropolitan area by population, with the Yellow Sea. After becoming an international port in 1883, POI has been the country’s logistics hub that handles and distributes both general cargoes and containers. It is ranked as the global 27th busiest port in terms of cargo volume and the 50th biggest container port [78]. Currently, POI has expanded into 125 berths, with a total of 26,031 meters of berth length, distributed to five main component ports as North Port, Inner Port, Coastal Port, South Port, and New Port. The North Port is an industrial port specialized in managing raw materials as scrap iron, feed-by products, hardwood, etc. With a lock-gate, a calm water level in the Inner Port is available for handling semiconductor equipment, automobiles, and precision machine parts. It also is a multipurpose port that serves other cargoes as grain, fruit, and general cargo. The South Port is available for handling small and medium containers and general cargoes. The New Port is under-constructed and specialized for handling medium and large containers. The Coastal Port, with international passenger and ferry terminals, focuses on serving passenger cruise ships and car ferries. Also, around these component ports, POI operates three smaller specialized ports named Geocheom-do Port, Song-do Port, and Yeongheung-do Port. The Song-do Port is specialized in fossil energy products, while the Geocheom-do Port is used for handling sand. The Yeongheung-do Port is constructed to support Yeongheung Thermal Power Plants. The capacity of POI is described in Table 2.Table 2
Summary of POI’s capacity. Source: [79].
No.PortHandling CapacityMain productsMax. Ship DWTBerthsBerth length (m)Bulk (1000RT)Container (1000TEU)1Inner Port50,000469,83838,161General cargo, iron, grains2South Port100,000253,64217,600762Chemical, cement, sand3Coastal Port50,00091,429Passenger, oil, LPG4North Port100,000266,42113,900Oil, general cargo, wood product5New Port3,00061,6002,100Container6Song-do75,00041,300LPG, oil7Yeongheung-do200,00051,12614,690Bituminous coal, limestone8Geocheom-do5,00046758,320SandSummary12526,03192,6712,86228,735.5The POI geographical segments are presented in Figure1. The study not only covered ship operation inside the port boundary but also the “affected zone” (within 5 km from the port boundary) to consider the effect of close-to-port emitted ship emissions, following EPA [80].Figure 1
Geographical segments of Port of Incheon. Note: images were downloaded from the Internet and Incheon Port Authority (IPA) website.


## 3.2. Activity Phase

A typical ship call often contains consecutive activity phases, which consist of a series of continuous activities that have similar features, categorized as “cruising, “maneuvering,” “hotelling at anchorage,” and “hotelling at berth.” Due to the VSR program was applied for a few last days in December 2019, the “decelerating phase” is assumed as not considered in this study. As defined in ICF [30] and by Song [2], at the “cruising” phase, the ship moves inside the port boundary, and all engines keep running. “Maneuvering” is the time that the ship transits between the breakwater (intersection of open sea and inland waterway) and berths at a slower speed. Even with tug assist, the propulsion engines are still in operation. “Hotelling at anchorage” is a period of waiting for berth call while “hotelling at berth” is the time that ships are docking at berth. In both “hotelling at anchorage” and “hotelling at berth” phases, only auxiliary engines (including boilers) still work with a peak load for providing on-board power. In the “hotelling at berth” phase, the ship can use shore power instead of turning on auxiliary engines, called cold ironing; however, cold ironing has been still not implemented in POI yet. The average speed and travel distance in each phase were provided by a local pilot company. Table 3 summarizes detailed activity phases information in POI.Table 3
Ship activity segments and features in POI.
Phase categoryPropulsion engineAuxiliary engineBoilerAvg. speed (knots)Travel distance (nm)1Hotelling at anchorageOffOnOn002CruisingOnOnOff12Varies by berth3ManeuveringOnOnOnAround 3.514Hotelling at berthOffOnOn00


## 3.3. Ship Classification

Ship characteristics (e.g., speed, engine size, and usage) vary considerably depending on the ship type. Previous studies classified ships in different ways; however, the common standard applied is a type of cargoes carried on the ship. This study followed ICF’s classification with 11 types of ships [30]. However, in POI, there are eight types available, involving bulk carrier, container ship, cruise ship, general cargo, reefer, roll-on/roll-off (RORO), tanker, and miscellaneous.


## 3.4. Data Collection and Analysis

### 3.4.1. Data Source

Due to the limitations of the AIS system mentioned above, in this study, in-port ship traffic data were collected from the VTS-based Korean Port Management Information System (Port-MIS) to consider ships with GT less than 300. A total of 16,677 ship calls were recorded from Port-MIS, as shown in Table4. Tanker fleet dominated ship calls at POI in 2019 with 42.6% of the total ship calls, following by the general cargo fleet and container ship fleet with 20.2% and 17.3%, respectively. Reefer fleet only contributed negligibly to total ship calls with 0.02% of total ship calls.Table 4
Ship call statistics (unit: ship call). Source: Port-MIS.
Ship typesNorth PortInner PortCoastal PortSouth PortNew PortOthersTotalBulk carrier253304—104——661Container ship271379681,871—2,886Cruise ship6367—123General cargo1,0566961491,40647203,374RoRo2777738330——1,847Reefer2—11——4Tanker4,1444736398361658477,104Miscellaneous113617720429294778Total5,6032,3271,6173,8562,1121,16216,677


### 3.4.2. Data Cleaning and Make-Up

A total of 47,626 statuses describing the entire marine traffic at POI during 2019 were collected. These statuses were noted in order of occurrence inside POI, which is different from the ship call’s timeline format. Two consecutive events in the same ship call usually are several hours apart; then, a ship call can last from a few hours to a few days. Also, a high level of traffic exaggerated the complexity level of the original data and made it an arduous task to handle nagging details. Therefore, a code, combined from call sign and time points, was applied to resort original data corresponding to the timeline of each ship call.Waypoints are used to present ship activities in data. A typical ship call contains waypoints as (1) entrance, clearance from the port area; (2) anchoring, collecting anchor; (3) berthing, unberthing; and (4) shift within the port area between berths (if need). A combination of two consecutive waypoints identifies a corresponding activity phase of the ship call. The time information of each waypoint is provided in the YYYY/MM/DD HH: MM format. The time of an activity phase is identified as the subtractions of 2 corresponding waypoints’ time. The average moving time and distance of the cruising and maneuvering process for each target port were calculated. In several cases, when only arrival and departure times were reported in the data, the total time of ship spent on the port is the difference between arrival time and departure time. Later, the corresponding average moving time and distance for the target port are applied, and then the hotelling time at the berth of that ship is the subtraction of average cruising and maneuvering time from the total time.


### 3.4.3. Ship Basic Data

Basic information about the ship, involving name, ship type, engine type, propulsion engine power, weight tonnage, and design speed, was obtained from the Korea Ship Safety Technology Authority and Korean Register. For missing propulsion engine power data, a simple linear approach to modeling the relationship between ship weight tonnage and propulsion engine power is applied by ship types. Linear models are shown in Figure2. The coefficients of determination (R2) are relatively high (>0.85), showing the high reliability of prediction.Figure 2
Linear regression of missing main engine data by ship type.Besides, installed auxiliary engine power also was not collected completely because there is a lack of information from manufacturers. However, it is impossible to apply a linear approach to estimate auxiliary engine power because there is no connection between installed auxiliary engine power and speed [51]. Also, it is difficult to recognize the using level of auxiliary engines in each activity phase during a ship call in the practice. Therefore, in previous studies, installed auxiliary engine power is commonly estimated by applying the default ratios by ship types with the total propulsion engine power of the ship [51]. This study applied the ratios revised from ICF [30], which is shown in Table 5.Table 5
The ratio for estimating auxiliary engine power. Source: ICF [30].
Ship typeAuxiliary to propulsion ratioBulk carrier0.222Container ship0.220Cruise ship0.278General cargo0.191RORO0.259Reefer0.406Tanker0.211Others0.100Lastly, the real values of revolutions per minute (RPM) of engines also were obtained from the Korea Ship Safety Technology Authority and Korean Register. RPM value of the engine helps to classify those engine speed designations: high-speed diesel (HSD), medium-speed diesel (MSD), and slow-speed diesel (SSD). Then, based on engine speed designations, proper EF values are matched in the calculation. In the case of missing values, they were made up of average values, as shown in Table6.Table 6
Average RPM values.
Ship typeAverage RPM values of propulsion enginesAverage RPM values of auxiliary enginesBulk carrier103843Container ship131847Cruise ship484895General cargo2711,225RoRo105787Reefer3111,043Tanker2861,147Miscellaneous4901,037


### 3.4.4. Activity Time

The ship transit time was calculated with the actual speed of the ship and the representative travel distance of each port (from border to berth). The average transit time by the port is shown in Table7. Among ports in POI, Inner Port shows the largest value because of a lock-gate.Table 7
Average transit time in POI 2019 (unit: hour).
PortAverage transit timeNorth Port2.2Inner Port3.0Coastal Port1.6South Port1.7New Port1.4Song-do Port1.4Geocheom-do Port2.5Yeongheung-do Port0.75The actual times in hotelling at anchorage and berth were obtained from the difference in time between anchoring, collecting anchor, and berthing, unberthing. To avoid anomalies, the ceilings are applied as 7 days (168 hours) for hotelling at anchorage and 14 days (336 hours) for hotelling at berth. The average docking time by ship type after applying ceilings is shown in Table8. Almost average values by ship type (except container ship’s average time) are considerably greater than the national average time applied in CAPSS (0.79 days∼18.96 hours). These differences may lead to considerable variances between estimated results and national estimations, especially tankers, general cargo ships, and container ships dominated the total number of ship calls at POI with the north of 80%.Table 8
Average time in hotelling at anchorage and berth in POI 2019 (unit: hour).
Types of shipAverage time in hotelling (hour)Compare to the national avg. time in CAPSS (%)Bulk carrier83.3439.3Container ship13.671.7Cruise Ship50.3265.3General cargo39.9210.4RORO25.9136.6Reefer89.0469.4Tanker22.0116.0Miscellaneous51.9273.7


## 3.4.1. Data Source

Due to the limitations of the AIS system mentioned above, in this study, in-port ship traffic data were collected from the VTS-based Korean Port Management Information System (Port-MIS) to consider ships with GT less than 300. A total of 16,677 ship calls were recorded from Port-MIS, as shown in Table4. Tanker fleet dominated ship calls at POI in 2019 with 42.6% of the total ship calls, following by the general cargo fleet and container ship fleet with 20.2% and 17.3%, respectively. Reefer fleet only contributed negligibly to total ship calls with 0.02% of total ship calls.Table 4
Ship call statistics (unit: ship call). Source: Port-MIS.
Ship typesNorth PortInner PortCoastal PortSouth PortNew PortOthersTotalBulk carrier253304—104——661Container ship271379681,871—2,886Cruise ship6367—123General cargo1,0566961491,40647203,374RoRo2777738330——1,847Reefer2—11——4Tanker4,1444736398361658477,104Miscellaneous113617720429294778Total5,6032,3271,6173,8562,1121,16216,677


## 3.4.2. Data Cleaning and Make-Up

A total of 47,626 statuses describing the entire marine traffic at POI during 2019 were collected. These statuses were noted in order of occurrence inside POI, which is different from the ship call’s timeline format. Two consecutive events in the same ship call usually are several hours apart; then, a ship call can last from a few hours to a few days. Also, a high level of traffic exaggerated the complexity level of the original data and made it an arduous task to handle nagging details. Therefore, a code, combined from call sign and time points, was applied to resort original data corresponding to the timeline of each ship call.Waypoints are used to present ship activities in data. A typical ship call contains waypoints as (1) entrance, clearance from the port area; (2) anchoring, collecting anchor; (3) berthing, unberthing; and (4) shift within the port area between berths (if need). A combination of two consecutive waypoints identifies a corresponding activity phase of the ship call. The time information of each waypoint is provided in the YYYY/MM/DD HH: MM format. The time of an activity phase is identified as the subtractions of 2 corresponding waypoints’ time. The average moving time and distance of the cruising and maneuvering process for each target port were calculated. In several cases, when only arrival and departure times were reported in the data, the total time of ship spent on the port is the difference between arrival time and departure time. Later, the corresponding average moving time and distance for the target port are applied, and then the hotelling time at the berth of that ship is the subtraction of average cruising and maneuvering time from the total time.


## 3.4.3. Ship Basic Data

Basic information about the ship, involving name, ship type, engine type, propulsion engine power, weight tonnage, and design speed, was obtained from the Korea Ship Safety Technology Authority and Korean Register. For missing propulsion engine power data, a simple linear approach to modeling the relationship between ship weight tonnage and propulsion engine power is applied by ship types. Linear models are shown in Figure2. The coefficients of determination (R2) are relatively high (>0.85), showing the high reliability of prediction.Figure 2
Linear regression of missing main engine data by ship type.Besides, installed auxiliary engine power also was not collected completely because there is a lack of information from manufacturers. However, it is impossible to apply a linear approach to estimate auxiliary engine power because there is no connection between installed auxiliary engine power and speed [51]. Also, it is difficult to recognize the using level of auxiliary engines in each activity phase during a ship call in the practice. Therefore, in previous studies, installed auxiliary engine power is commonly estimated by applying the default ratios by ship types with the total propulsion engine power of the ship [51]. This study applied the ratios revised from ICF [30], which is shown in Table 5.Table 5
The ratio for estimating auxiliary engine power. Source: ICF [30].
Ship typeAuxiliary to propulsion ratioBulk carrier0.222Container ship0.220Cruise ship0.278General cargo0.191RORO0.259Reefer0.406Tanker0.211Others0.100Lastly, the real values of revolutions per minute (RPM) of engines also were obtained from the Korea Ship Safety Technology Authority and Korean Register. RPM value of the engine helps to classify those engine speed designations: high-speed diesel (HSD), medium-speed diesel (MSD), and slow-speed diesel (SSD). Then, based on engine speed designations, proper EF values are matched in the calculation. In the case of missing values, they were made up of average values, as shown in Table6.Table 6
Average RPM values.
Ship typeAverage RPM values of propulsion enginesAverage RPM values of auxiliary enginesBulk carrier103843Container ship131847Cruise ship484895General cargo2711,225RoRo105787Reefer3111,043Tanker2861,147Miscellaneous4901,037


## 3.4.4. Activity Time

The ship transit time was calculated with the actual speed of the ship and the representative travel distance of each port (from border to berth). The average transit time by the port is shown in Table7. Among ports in POI, Inner Port shows the largest value because of a lock-gate.Table 7
Average transit time in POI 2019 (unit: hour).
PortAverage transit timeNorth Port2.2Inner Port3.0Coastal Port1.6South Port1.7New Port1.4Song-do Port1.4Geocheom-do Port2.5Yeongheung-do Port0.75The actual times in hotelling at anchorage and berth were obtained from the difference in time between anchoring, collecting anchor, and berthing, unberthing. To avoid anomalies, the ceilings are applied as 7 days (168 hours) for hotelling at anchorage and 14 days (336 hours) for hotelling at berth. The average docking time by ship type after applying ceilings is shown in Table8. Almost average values by ship type (except container ship’s average time) are considerably greater than the national average time applied in CAPSS (0.79 days∼18.96 hours). These differences may lead to considerable variances between estimated results and national estimations, especially tankers, general cargo ships, and container ships dominated the total number of ship calls at POI with the north of 80%.Table 8
Average time in hotelling at anchorage and berth in POI 2019 (unit: hour).
Types of shipAverage time in hotelling (hour)Compare to the national avg. time in CAPSS (%)Bulk carrier83.3439.3Container ship13.671.7Cruise Ship50.3265.3General cargo39.9210.4RORO25.9136.6Reefer89.0469.4Tanker22.0116.0Miscellaneous51.9273.7


## 4. Methodology

To provide a multidimensional overview of in-port ship pollution, this study applied both EO and FC methods to figure out port emissions emitted in POI during 2019 with the activity time data collected from the VTS system. The methodology of the bottom-up approach was developed by improving the general estimation equations from the previous works, which are reviewed in the literature review. The total volume of emissions discharged from a ship call is the aggregation of emissions from the combustion of all engines operated in all activity phases. Then, emission results from calls were scaled up to reach the entire POI emission in 2019. Also, CF is used for adjusting emission factors, and in this study, CF will be assumed to be 1, meaning there are no emission reduction technologies installed on ships.


### 4.1. EO Method

EO method adheres to the generic equation mentioned in the literature review part. The energy demand here considers the energy output of the engine over the operating time. The equation for the EO method applied in this study considers all activity phases, air pollutants, and engines in a ship call as follows:(7)Es,i=∑ph,jTs,ph∑jPs,j×LFs,j×EFi,j,ph,where E represents the total volume of emission emitted over a complete ship call (g); T represents activity time (hour); P represents engine power (kW); LF represents load factor (%); EF represents emission factor (g/kWh); s represents ship call; i represents pollutant; ph represents the activity phase of a ship call (hotelling at anchorage, cruising, maneuvering, and hotelling at berth); and j represents engine type (propulsion engine and auxiliary engine).The propeller law is applied to calculate the propulsion engine load factor (LFm), as shown in the following equation:(8)LFm=ASMS3,where AS = ship actual speed (knots) and MS = ship maximum speed (knots), which is defined by the manufacturer. In the case of the auxiliary engine load factor (LFa), because it varies by ship type and activity phases, and the data sources are limited, therefore, ICF [30] also suggested default LFa assumptions in his study. This study referred to these values for calculation.The non-CO2 EFs used here were referred from Lee et al. [70]. For CO2 emission, the carbon fraction in the fuel used is defined at 86.8% by weight, and a factor of molecular weight of C and CO2 equals 44/12. Using SFCs referred from ICF [30], the EFs for CO2 emission are calculated as(9)CO2EF=0.868.4412.SFC.As an international hub port, most ships visited POI are defined as international ocean-going ships; therefore, international EFs estimation methods suggested in them are deemed adaptable for this study. Applied EFs in the EO method are summarized in Table9.Table 9
Applied EFs for EO method (unit: g/kWh). Source: Lee et al. [70].
Engine typePhaseCOCO2NOxSOxPM10PM2.5VOCsNH3BSFCPropulsion-HSDCruising1.1646.1412.03.970.470.430.210.0014203Propulsion-MSDCruising1.1646.1413.23.970.470.430.630.0014203Propulsion-SSDCruising0.5588.8517.03.620.450.420.530.0013185Propulsion-HSDManeuvering2.2709.89.64.360.500.460.630.0016223Propulsion-MSDManeuvering2.2709.810.64.360.500.461.580.0016223Propulsion-SSDManeuvering1.0649.3213.63.990.470.441.900.0014204AuxiliaryAll1.1690.7113.94.240.490.450.420.0015217


### 4.2. FC Method

Similar to the EO method, the FC method’s equation adheres to the generic equation and considers fuel consumption of the engine over the operating time to figure out energy demand. The detailed equation for fuel consumption applied in this study is referred from Corbett et al. [31]; then, the final equation for the FC method is formed and shown in the following equation:(10)Es,i,ph=∑ph,jFCph,j×EFj,i,ph=∑phTph×MFph×ASMS3×EFm,i,ph+AFph×EFa,i,ph.


## 4.1. EO Method

EO method adheres to the generic equation mentioned in the literature review part. The energy demand here considers the energy output of the engine over the operating time. The equation for the EO method applied in this study considers all activity phases, air pollutants, and engines in a ship call as follows:(7)Es,i=∑ph,jTs,ph∑jPs,j×LFs,j×EFi,j,ph,where E represents the total volume of emission emitted over a complete ship call (g); T represents activity time (hour); P represents engine power (kW); LF represents load factor (%); EF represents emission factor (g/kWh); s represents ship call; i represents pollutant; ph represents the activity phase of a ship call (hotelling at anchorage, cruising, maneuvering, and hotelling at berth); and j represents engine type (propulsion engine and auxiliary engine).The propeller law is applied to calculate the propulsion engine load factor (LFm), as shown in the following equation:(8)LFm=ASMS3,where AS = ship actual speed (knots) and MS = ship maximum speed (knots), which is defined by the manufacturer. In the case of the auxiliary engine load factor (LFa), because it varies by ship type and activity phases, and the data sources are limited, therefore, ICF [30] also suggested default LFa assumptions in his study. This study referred to these values for calculation.The non-CO2 EFs used here were referred from Lee et al. [70]. For CO2 emission, the carbon fraction in the fuel used is defined at 86.8% by weight, and a factor of molecular weight of C and CO2 equals 44/12. Using SFCs referred from ICF [30], the EFs for CO2 emission are calculated as(9)CO2EF=0.868.4412.SFC.As an international hub port, most ships visited POI are defined as international ocean-going ships; therefore, international EFs estimation methods suggested in them are deemed adaptable for this study. Applied EFs in the EO method are summarized in Table9.Table 9
Applied EFs for EO method (unit: g/kWh). Source: Lee et al. [70].
Engine typePhaseCOCO2NOxSOxPM10PM2.5VOCsNH3BSFCPropulsion-HSDCruising1.1646.1412.03.970.470.430.210.0014203Propulsion-MSDCruising1.1646.1413.23.970.470.430.630.0014203Propulsion-SSDCruising0.5588.8517.03.620.450.420.530.0013185Propulsion-HSDManeuvering2.2709.89.64.360.500.460.630.0016223Propulsion-MSDManeuvering2.2709.810.64.360.500.461.580.0016223Propulsion-SSDManeuvering1.0649.3213.63.990.470.441.900.0014204AuxiliaryAll1.1690.7113.94.240.490.450.420.0015217


## 4.2. FC Method

Similar to the EO method, the FC method’s equation adheres to the generic equation and considers fuel consumption of the engine over the operating time to figure out energy demand. The detailed equation for fuel consumption applied in this study is referred from Corbett et al. [31]; then, the final equation for the FC method is formed and shown in the following equation:(10)Es,i,ph=∑ph,jFCph,j×EFj,i,ph=∑phTph×MFph×ASMS3×EFm,i,ph+AFph×EFa,i,ph.


## 5. Results and Discussion

This section reports a multidimensional bottom approach in-port ship emission inventory in POI during 2019, including two inventories estimated under the two methods mentioned above. These inventories figured out in-port ship emissions in seven geographical segments (at sea, five main ports, and other specialized ports), by eight types of air pollutants and by four ship activity phases. There was considerable variance between results from the two methods. EO method reported 201,612.26 tonnes of CO2, 323.85 tonnes of CO, 4,097.79 tonnes of NOx, 1,237.92 tonnes of SOx, 145.06 tonnes of PM10 (including PM2.5), 133.32 tonnes of PM2.5, 136.21 tonnes of VOCs, and 0.44 tonne of NH3, while FC method showed 193,981.71 tonnes, 452.88 tonnes, 3,573.27 tonnes, 1,224.01 tonnes, 88.21 tonnes, 82.33 tonnes, 116.09 tonnes, and 0.43 tonne, respectively. The total fuel consumption for ships at POI during 2019 (reported in FC method) was nearly 61,193 tonnes.In general, the EO method reported higher values than the FC method’s value with almost all pollutants, except for CO emissions. Also, there are big gaps (in percentage) between the calculated volumes of NOx, CO, and PM emissions of the two inventories. However, both inventories mutualized in the rank of emitted pollutants in POI’s emission inventory. CO2 was the dominating air pollutant in both inventories at POI during 2019, accounting for about 97% of the total amount of emitted emissions. Among the rest, NOx was the most serious air pollutant in both inventories at POI during 2019, which covers over 64% of the total amount of non-CO2 emissions, followed by SOx emissions with over 20%. The contribution of CO ranged from 5–8%. NH3 just accounted for nearly 0.01% of the total amount of non-CO2 emissions. Other pollutants’ shares in both inventories were also insignificant with less than 2.5% of the total non-CO2 volume. The comparison of EO with FC methods about non-CO2 emissions is shown in Figure 3.Figure 3
Non-CO2 emissions comparison between FC and EO methods (unit: tonnes).Figure4 presents the inventories of ship emission by geographical areas during 2019, considering the sea area (when ships moved on the area of the sea from affected zone to berths) and port areas (when ships docked at berths, involving five main component ports and three others). Both inventories agreed that almost all ship emissions were exhausted on sea area. Also, the big difference between the two inventories happened in sea area (7%), while they showed similar volumes at berths. The emission share on sea area is 53–55% of the total amount of emissions. Among port areas, Inner Port was the most polluted port, accounting for 14% of the ship emissions, followed by North Port and South Port with a ratio of 9% and 7%, respectively. Coastal Port shared the smallest proportion, compared to five main component ports, with almost 4% of the total amount of emission.Figure 4
POI 2019 ship emission inventories by geographical areas (Unit: tonnes).The great contribution of Inner Port, North Port, and South Port to the air pollution at POI has been anticipated because they are the most important and busiest ports at POI. With a lock-gate, the Inner Port is the ideal port for handling car ferries; therefore, the RORO fleet (mainly car ferries) supplies mainly for the first rank of Inner Port. The tanker fleet was distributed in all ports of POI; however, the concentration in a large number of tankers in North Port pushed it to become the second polluted port at POI. In the case of South Port, general cargo and container ship fleets are the main emission contributors there.Next, Figure5 illustrates the emission contribution by ship types at POI during 2019. Tanker fleet was agreed as the most polluted fleet at POI, with around 29% of total in-port emission volume, followed by general cargo, RORO, and container ship. However, two inventories show a dissimilar view in the ranks of these three groups. While EO inventory reports that second rank belongs to container ship fleet, which contributed up to 20.1% of emission volume, followed closely by RORO fleet and general cargo fleet with the ratio of 19.8% and 19.6%, respectively; in case of FC inventory, the RORO fleet raised to occupy the second position with 20.3%, pushed general cargo fleet down to the third rank with 19.8%, and container ship fleet only contributed 19.7% of the total amount of emissions, ranked fourth. The top-4 polluted fleets supplied 88.7% of the total emission volume. In both inventories, the bulk carrier fleet also was a considerable polluted source, which contributed over 8.8% of the total amount of emission, ranked fifth. Other fleets present minor contributions, fewer than 2%. Therefore, IPA should focus on emission reduction from tanker, general cargo, RORO, and container ship.Figure 5
POI 2019 ship emission inventory by ship type (unit: tonnes).Then, Figure6 shows ship emission inventories by activity phases. Although the rankings and contribution (in percentage) of cruising and hotelling at berth are interchanged in two inventories, however, it is undeniable that they overshadowed others in total emitted emissions. During hotelling at berth, the propulsion engine is in respite, and then auxiliary engines monopolize in emitting all emissions. Therefore, besides the implementation of ECA, to reduce the volume of in-port ship emissions, efforts to cut down emissions during cruising and hotelling at berth phases are indispensable.Figure 6
POI 2019 ship inventory by ship activity phase (unit: tonnes).


## 6. Ship-Related Emissions Reduction Policies at POI

As discussed above, several policies are designed to reduce the volume of emissions emitted from ships operated at POI as (1) local ECA realization, (2) VSR program motivation, (3) application of cold ironing, and (4) establishment of a national integrated emission platform. In this study, the above-calculated EO emission inventory is adopted as the baseline emission inventory and solutions (1)–(3) are simulated and evaluated based on EO method if assuming that these suggested policies are applied at POI. The potential benefit of solution (4) also is discussed.


### 6.1. Assessing the Environmental Benefit of Designating the Local ECA

The great contribution of NOx and SOx emissions among non-CO2 emissions to air pollution at POI suggests the necessity of designating a local NOx and SOx emission control area (ECA). In ECA, ships will be required to follow strictly Tier III in Regulation 13 of Annex VI of MARPOL 73/78 about reducing NOx emissions and the sulphur contents in the fuel used will be reduced to 0.1%. This idea was suggested firstly by Chang et al. [49]; however, it has not been implemented yet in Incheon or any other Korean ports. As of 2019, the required fuel used inside the POI area was MGO with the maximum sulphur content at 1.0% by mass. Then, with SOx ECA, ships have to switch from high sulphur content fuel (1.0%) to a very lower one (0.1%). Assuming that energy efficiency is unchanged, since the volumes of SOx and PM emissions released are dependent critically on the sulphur content in the fuel used according to estimating equations applied in Lee et al. [70], the SOx ECA could help decrease 10 times the amount of SOx emission and nearly 2.5 times the amount of PM emission. The cutdown volume of emissions will be 1114.12 tonnes of SOx and 71.07 tonnes of PM10 (including 65.25 tonnes of PM2.5) if the SOx ECA is established.Different from SOx emission, in the case of NOx emission, the abatement is implemented mostly by technology improvement in the engine’s combustion and selective catalytic reduction (SCR) system equipped onboard. Due to the limitation of the estimating method, it is difficult to simulate and estimate the expected reduction in NOx emission; therefore, the reduction of NOx is not considered in this study.


### 6.2. Assessing the Environmental Benefit of the VSR Program

In the last few days of December 2019, the Korean Ministry of Oceans and Fisheries started designating the “Vessel Speed Reduction Program” in five Korean main ports, including POI, to lessen PM levels during winter. Ships are motivated to transit slower than certain speeds (12 knots for container ships and car-carriers and 10 knots for others). To facilitate the analysis, it is assumed that there was no RSZ in POI during 2019 when calculating the baseline emission inventory and then simulating the effect of RSZ if it is designated. With lower cruising speed (10 knots), the ships will take more cruising time but lower load factor for propulsion engines (see equation (6)). Then, an RSZ can help reduce 7-8% amount of each pollutant: 22.77 tonnes of CO, 295.83 tonnes of NOx, 87.03 tonnes of SOx, 10.45 tonnes of PM10 (including 9.61 tonnes of PM2.5), 11.15 tonnes of VOC, 0.03 tonne of NH3, and 14,167.93 tonnes of CO2.


### 6.3. Assessing the Environmental Benefit of Applying Cold Ironing

In the case of hotelling at berth, emission could be cut down through reducing docking time or cut-downing emissions from docking time. In practice, reducing docking time is complicating and challenging because it is related to overhauling, optimizing, and scheduling the operation of a chain including other port facilities such as quay cranes, internal trucks, and yard cranes. In contrast, reducing emissions from docking time could be achieved easily by using an on-shore power supply. Assuming that 50% of hotelling time in each port is applied this technology, it will help reduce 22-23% amount of each pollutant: 72.26 tonnes of CO, 913.14 tonnes of NOx, 278.54 tonnes of SOx, 32.19 tonnes of PM10 (including 29.56 tonnes of PM2.5), 27.59 tonnes of VOC, 0.1 tonne of NH3, and 45,375.03 tonnes of CO2.If three of them are conducted together, this will help cut down 95.03 tonnes of CO, 1208.97 tonnes of NOx, 1150.68 tonnes of SOx, 92.51 tonnes of PM10 (including 84.97 tonnes of PM2.5), 38.74 tonnes of VOC, 0.13 tonne of NH3, and 59,542.96 tonnes of CO2 corresponding to 29%, 30%, 93%, 64% (64%), 28%. 30% and 30% of the total emitted amount of each pollutant in the 2019 baseline inventory.


### 6.4. Establishment of a National Integrated Emission Platform

Besides counter-measures, a good assessment and management system is indispensable for managing green reform at ports. With the widespread Internet, the combination of the platform industry and assessments of environmental situations around seaport would be an interesting and promising research topic shortly. With the platform, users can access and interact with the system actively and easily. This study suggests an idea about a national integrated platform that can standardize and systematize procedures in monitoring harmful ship-related emissions released from ports in Korea. This platform will provide three main functions: (1) data collection and analysis; (2) estimation, modification, and visualization; and (3) prediction. The abovementioned three functions are proposed to ensure and emphasize the systematicness to generate a better air pollutant assessment and management.The necessary input data for emission estimation is linked with and automatic-synchronized simultaneously from the national Port-MIS system. Then, based on collected data, the system analyses, estimates, and visualizes the volume of emissions emitted or environment indexes promptly to warn the port operator if the negative impacts on the local community from in-port emissions exceed a certain “safe” level. Then, input data adjustment in the modification function allows port operators to adjust or reschedule the port operation, through suppositions and corresponding input adjustment, to reduce geographic emissions to the “safe” level. Also, by machine-learning algorithms and statistical analyses, the prediction function helps port operators predict the trend of air pollution at ports and then plan port operations appropriately. Also, the researchers use modification and prediction function in annual or seasonal emission inventory with their recommended air management policies or counter-measures to assess the performances of them.


## 6.1. Assessing the Environmental Benefit of Designating the Local ECA

The great contribution of NOx and SOx emissions among non-CO2 emissions to air pollution at POI suggests the necessity of designating a local NOx and SOx emission control area (ECA). In ECA, ships will be required to follow strictly Tier III in Regulation 13 of Annex VI of MARPOL 73/78 about reducing NOx emissions and the sulphur contents in the fuel used will be reduced to 0.1%. This idea was suggested firstly by Chang et al. [49]; however, it has not been implemented yet in Incheon or any other Korean ports. As of 2019, the required fuel used inside the POI area was MGO with the maximum sulphur content at 1.0% by mass. Then, with SOx ECA, ships have to switch from high sulphur content fuel (1.0%) to a very lower one (0.1%). Assuming that energy efficiency is unchanged, since the volumes of SOx and PM emissions released are dependent critically on the sulphur content in the fuel used according to estimating equations applied in Lee et al. [70], the SOx ECA could help decrease 10 times the amount of SOx emission and nearly 2.5 times the amount of PM emission. The cutdown volume of emissions will be 1114.12 tonnes of SOx and 71.07 tonnes of PM10 (including 65.25 tonnes of PM2.5) if the SOx ECA is established.Different from SOx emission, in the case of NOx emission, the abatement is implemented mostly by technology improvement in the engine’s combustion and selective catalytic reduction (SCR) system equipped onboard. Due to the limitation of the estimating method, it is difficult to simulate and estimate the expected reduction in NOx emission; therefore, the reduction of NOx is not considered in this study.


## 6.2. Assessing the Environmental Benefit of the VSR Program

In the last few days of December 2019, the Korean Ministry of Oceans and Fisheries started designating the “Vessel Speed Reduction Program” in five Korean main ports, including POI, to lessen PM levels during winter. Ships are motivated to transit slower than certain speeds (12 knots for container ships and car-carriers and 10 knots for others). To facilitate the analysis, it is assumed that there was no RSZ in POI during 2019 when calculating the baseline emission inventory and then simulating the effect of RSZ if it is designated. With lower cruising speed (10 knots), the ships will take more cruising time but lower load factor for propulsion engines (see equation (6)). Then, an RSZ can help reduce 7-8% amount of each pollutant: 22.77 tonnes of CO, 295.83 tonnes of NOx, 87.03 tonnes of SOx, 10.45 tonnes of PM10 (including 9.61 tonnes of PM2.5), 11.15 tonnes of VOC, 0.03 tonne of NH3, and 14,167.93 tonnes of CO2.


## 6.3. Assessing the Environmental Benefit of Applying Cold Ironing

In the case of hotelling at berth, emission could be cut down through reducing docking time or cut-downing emissions from docking time. In practice, reducing docking time is complicating and challenging because it is related to overhauling, optimizing, and scheduling the operation of a chain including other port facilities such as quay cranes, internal trucks, and yard cranes. In contrast, reducing emissions from docking time could be achieved easily by using an on-shore power supply. Assuming that 50% of hotelling time in each port is applied this technology, it will help reduce 22-23% amount of each pollutant: 72.26 tonnes of CO, 913.14 tonnes of NOx, 278.54 tonnes of SOx, 32.19 tonnes of PM10 (including 29.56 tonnes of PM2.5), 27.59 tonnes of VOC, 0.1 tonne of NH3, and 45,375.03 tonnes of CO2.If three of them are conducted together, this will help cut down 95.03 tonnes of CO, 1208.97 tonnes of NOx, 1150.68 tonnes of SOx, 92.51 tonnes of PM10 (including 84.97 tonnes of PM2.5), 38.74 tonnes of VOC, 0.13 tonne of NH3, and 59,542.96 tonnes of CO2 corresponding to 29%, 30%, 93%, 64% (64%), 28%. 30% and 30% of the total emitted amount of each pollutant in the 2019 baseline inventory.


## 6.4. Establishment of a National Integrated Emission Platform

Besides counter-measures, a good assessment and management system is indispensable for managing green reform at ports. With the widespread Internet, the combination of the platform industry and assessments of environmental situations around seaport would be an interesting and promising research topic shortly. With the platform, users can access and interact with the system actively and easily. This study suggests an idea about a national integrated platform that can standardize and systematize procedures in monitoring harmful ship-related emissions released from ports in Korea. This platform will provide three main functions: (1) data collection and analysis; (2) estimation, modification, and visualization; and (3) prediction. The abovementioned three functions are proposed to ensure and emphasize the systematicness to generate a better air pollutant assessment and management.The necessary input data for emission estimation is linked with and automatic-synchronized simultaneously from the national Port-MIS system. Then, based on collected data, the system analyses, estimates, and visualizes the volume of emissions emitted or environment indexes promptly to warn the port operator if the negative impacts on the local community from in-port emissions exceed a certain “safe” level. Then, input data adjustment in the modification function allows port operators to adjust or reschedule the port operation, through suppositions and corresponding input adjustment, to reduce geographic emissions to the “safe” level. Also, by machine-learning algorithms and statistical analyses, the prediction function helps port operators predict the trend of air pollution at ports and then plan port operations appropriately. Also, the researchers use modification and prediction function in annual or seasonal emission inventory with their recommended air management policies or counter-measures to assess the performances of them.


## 7. Conclusion

As a result of the quick spreading of global trade, atmospheric pollution from transportation has been more topical than ever, especially in dense hub port-cities. Given environmental legislation that is enforced currently, the shipping industry has not been paid adequate attention to its relatively considerable volume of emitted emissions. Besides published national reports, the awareness of atmospheric pollution from growing port traffic has invoked reliable and up-to-date regional comprehensive assessments of egregious environmental impacts from in-port ship traffic with timely counter-measures. This study contributes to the stream of in-port ship emission research both theoretically and practically. Besides the sole application of AIS for inputs in the bottom-up approach studies, VTS could be considered as a promising confident alternative. In Korea, with support from the government, the VTS-based Port-MIS system would be a reliable data source for local studies.This study generated a comprehensive annual ship emission inventory at POI in 2019, according to geographical areas, operational phases, and ship types by the bottom-up approach. Although there are small differences between the results of the two methods EO and FC, and it is impossible to demonstrate which one is better, however, two inventories agree about the dominance of CO2 emission and the considerable volume of NOx and SOx emissions in the total volume of port emissions, the significant contribution of auxiliary engines during hotelling at berth phase and the 4-most polluted ship types at POI during the year of 2019. Also, this study suggests and simulates the effects of four applicable green policies for the POI situation in the practice. The combination of the three abovementioned policies could help slash the significant volume of emitted CO (29%), NOx (30%), SOx (93%), PM10 and PM2.5 (64%), VOC (28%), NH3 (30%), and CO2 (30%).Given most ships visited POI could be considered as international waterborne navigation, using default inputs in the estimation process (e.g., EF, auxiliary characteristics data, and engine load factor), that referred from international studies, is deemed acceptable. Moreover, boilers also were not collected completely and considered because there is a lack of information from manufacturers and ship owners. Therefore, for better estimation, the authors also recommend the adoption of local values to reduce uncertainty. However, currently, Korean data are limited; therefore, it garners the attention of investigating local-specific inputs. Moreover, an idea of the automatic input collecting system be synchronized with our proposed integrated system could be another interesting topic for further research.Although both inventories also reported an insignificant volume of PM emissions, in terms of weight contribution, however, if compared their size and weight to other pollutants, they would be a potential threat to the port environment as well as the local community. Therefore, besides figuring out the weight of the total volume of PM emitted, the evaluation of spatial seasonal PM dispersion in the port area would be extremely necessary. This evaluation will contribute considerably to ship emission inventory and the deep understanding of the affection of in-port ship PM emissions.Finally, the port operation also contains other important activities with various types of cargo handling equipment, railcars, and drayage trucks, and they contribute significantly to port pollution. Therefore, a broader view of a comprehensive in-port emission inventory or a port-scale integrated system, which presents emissions not only from ships but also from other port-related land-based vehicles, should be considered soon to understand completely about port pollution. Of course, this complete inventory would be a great baseline for other in-port environment evaluating studies in Incheon as well as other similar regions.



---

*Source: 5568777-2021-04-10.xml*