A Line Ship Emissions while Manoeuvring and Hotelling—A Case Study of Port Split

: Strategically, the Republic of Croatia, with its economy focused on tourism, is directly connected to the sea and coastal area, and integrated management of this area contributes to the sustainable development strategy. Worldwide, the problem of atmospheric pollution from maritime tra ﬃ c is a poorly researched area, especially when this type of tra ﬃ c is continuously growing. On the example of Port Split, the paper aims to present the following emission, carbon dioxide (CO 2 ), nitrogen oxides (NO X ), carbon monoxide (CO), volatile organic compounds (VOC), sulphur dioxide (SO 2 ), particulate matter (PM) and black carbon, of line vessels during manoeuvring and hotelling phase for 2017, 2018 and 2019. Furthermore, the statistical analysis and appropriate conclusions have been performed on CO 2 since all other emissions are linearly dependent. From the analyses in the hotelling and manoeuvring phase of line ships, it can be concluded that during 2019 there was a slight increase in emissions, but overall there was no signiﬁcant increase in the number of line vessels and increased tra ﬃ c. The obtained results of case study of port Split provide recommendations leading to further reduction of harmful gas emission, monitoring them, and integrating it into management of urban ports.


Introduction
Maritime transport is considered to be the most energy-efficient mode of transport since it can carry the largest amount of cargo with the least energy consumed. Unfortunately, it has bad effects on the environment and human health [1]. Pollution from ships is not limited only to maritime accidents but also to the regular navigation and ship operations. The main ship impacts on the environment are sea discharges, gas emissions, and noise [2].
Sea pollution from ships has a visible impact on the surrounding area and countries since is transferred transboundary in the atmosphere and globally affects the air quality [3]. Air emissions include pollutants and greenhouse gases. The main focus of this paper are major pollutants (nitrogen oxides (NO X ), carbon monoxide (CO), volatile organic compounds (VOC), sulphur dioxide (SO 2 ), particulate matter (PM) and greenhouse gas (carbon dioxide (CO 2 ) [4]. Furthermore, there are a lot of sources of pollution that can come from the ship, such as oil spills, ballast waters, grey waters, black waters, anti-fouling paint, noise, solid waste, as shown in Figure 1 [5]. According to the World Health Organization (WHO), 4.2 million people died in 2016 due to air pollution [6]. Comparing the emission values of five commonly used inventories (EMEP; TNO-MACC_III; E-PRTR; EDGAR and STEAM), the contribution of shipping to overall emissions is 16% for NOX 11% for SOX and 5% for PM10 [7]. Shipping emissions contribute 1-7% of ambient air PM10 levels, 1-14% of PM2.5, and at least 11% of PM1 in European coastal areas. Shipping emissions contribute with 7-24% of ambient air NO2. The highest values of NO2 have been recorded in the Netherlands and Denmark [8]. In the Mediterranean area locally released NOX is mainly responsible for the production of ozone. Excluding NOX emissions NOX emissions from ships in model would reduce the surface ozone concentration by 15% [8]. Ship emissions are affecting air quality and human health in the coastal communities. NO2 and CO-emissions in ports are connected to bronchitic symptoms, and exposure to SO2 emissions is connected with respiratory issues and premature births [9]. If the legislation and regulations on land result in a reduction of emissions from sources on land, the impact of maritime transport at the global level will increase, especially taking into account its continued growth. Therefore, by placing limits on harmful gases from maritime transport, this impact will be minimized.
Although the International Convention for the Prevention of Pollution from Ships (MARPOL) has a major role in regulating shipping pollution, several directives are adopted to supplement or clarify specific fields of interest. Table 1 shows the organizations, conventions, and laws related to ship emissions.  According to the World Health Organization (WHO), 4.2 million people died in 2016 due to air pollution [6]. Comparing the emission values of five commonly used inventories (EMEP; TNO-MACC_III; E-PRTR; EDGAR and STEAM), the contribution of shipping to overall emissions is 16% for NO X 11% for SO X and 5% for PM 10 [7]. Shipping emissions contribute 1-7% of ambient air PM 10 levels, 1-14% of PM 2.5 , and at least 11% of PM 1 in European coastal areas. Shipping emissions contribute with 7-24% of ambient air NO 2 . The highest values of NO 2 have been recorded in the Netherlands and Denmark [8]. In the Mediterranean area locally released NO X is mainly responsible for the production of ozone. Excluding NO X emissions NO X emissions from ships in model would reduce the surface ozone concentration by 15% [8]. Ship emissions are affecting air quality and human health in the coastal communities. NO 2 and CO-emissions in ports are connected to bronchitic symptoms, and exposure to SO 2 emissions is connected with respiratory issues and premature births [9]. If the legislation and regulations on land result in a reduction of emissions from sources on land, the impact of maritime transport at the global level will increase, especially taking into account its continued growth. Therefore, by placing limits on harmful gases from maritime transport, this impact will be minimized.
Although the International Convention for the Prevention of Pollution from Ships (MARPOL) has a major role in regulating shipping pollution, several directives are adopted to supplement or clarify specific fields of interest. Table 1 shows the organizations, conventions, and laws related to ship emissions. The first column shows the organizations, conventions, and laws. The second column contains abbreviations. Further, the third column contains years when they are founded/entered into force, and the fourth column briefly indicates their role. For example, IMO is an organization with responsibility for the safety and security of shipping and the prevention of marine and atmospheric pollution by ships.
Carbon dioxide in shipping occurs during the combustion of fossil fuels. According to the fourth IMO study (2020) on greenhouse gases emissions, maritime contribution to greenhouse gas emissions is 2.89% [10]. The same study shows that transporting goods by ships is responsible for approximately 1,056 million tons of CO 2 annually. The projections show that shipping emissions could increase between 90% and 130% of 2008 emissions by 2050 [10].
The combustion of fossil fuels emits various sulphur oxides (SO X ). As reported by several authors, shipping contributes to SO X 's overall anthropogenic emission from 5 to 8% [11,12]. The percentage of sulphur in fossil fuels can vary, depending on the fuel type. Sulphur oxides have a negative impact both on human health and on the environment [13,14]. The IMO is reducing the percentage of sulphur used in marine fossil fuels from year to year. The sulphur content of any fuel used on board shall not exceed the following limits: 4.50% by weight before 1 January 2012; 3.50% by weight from 1 January 2012 onwards; 0.50% by weight from 1 January 2020 onwards. According to the EU Directive 2005/33/EC, while at berths in ports, all ships must use fuel with sulphur content less than 0.1 by weight [15]. The same strict limit of 0.10% m/m. has already been applied in the emission control areas (ECAS), set by the International Maritime Organization [16]. Another way of air pollutants limitation is by installing exhaust gas cleaning systems ("scrubbers"). Ships with installed scrubbers can continue use heavy fuel oil of 3.5% sulphur content [17]. Nitrogen oxide emissions (NO X ) from ships are forming when fuel burns at high temperatures in the ship's internal combustion engine. The overall ship sector contributes to the anthropogenic emission of NO X by 15% [11]. NO X emissions affect the environment by causing acid rain. When combined with VOC, ground-level ozone is formed, and impacts human health [13,14]. NO X emissions are regulated by the MARPOL Annex VI. The Different levels (Tiers) of control apply based on the ship construction date. The Tier I regulation refers to ships built after 1 January 2000, the Tier II regulation refers to ships built after 2011, and Tier III refers to ships built after 2016. The NO X emission limits vary for the slow-speed engines (<130 rpm) the high-speed engines (>2000 rpm) and the intermediate speed engines (130 < n < 2000 rpm) [18].
The particulate matter (PM) is the aerosols, consisting of mixtures of solid particles and liquid droplets found in the air, and they are defined by size. PM 10 are defined by size. PM 10 are inhalable particles with a diameter larger than 2.5 micrometres and smaller than 10 micrometres and PM 2.5 are fine particles that are 2.5 micrometres and smaller. The PMs are affecting human health by the lungs, causing inflammation, and restricting the passage of oxygen to the blood [13,14]. According to the European Environment Agency (EEA) PM 2.5 concentrations in 2016 were responsible for more than 412,000 premature deaths due to long-term exposure in Europe [19].
The non-methane volatile organic compounds (NMVOCs) are a collection of organic compounds that differ widely in chemical composition when emitted into the atmosphere from a large number of sources, including combustion. NMVOCs have a negative impact on the environment and human health [13,19].
Air pollution from ships can be estimated on a global and/or local scale. The global scale impact implies emissions during the ship navigation, while the local implies emissions in the ports or nearby ports [13]. Approximately 70% of the ship emissions are estimated to occur within 400 km of land and can significantly influence the air quality of a coastal area [5].
Several studies on port emissions are related to shipping, but it is difficult to compare their results since they use different methodologies [9]. Several methodologies are used for estimating emissions, which can be summed into a bottom-up approach and a fully top-down approach [20]. Reviewing the literature [21][22][23][24][25][26][27][28][29][30][31][32][33][34], comparing the methodologies used in bottom-up and top-down methods, and taking into account that our research problem was based on ship activity data, the bottom-up method has been used. The bottom-up method uses more data from the Automatic Identification System (AIS) such as ship characteristics, ship phase, loading factors, and the time spent in each stage [13]. In contrast, the top-down method uses sold fuel and the fuel emission factor [35]. The uncertainties of top down methodology are based on question are bunker fuel sale statistics representative [5]. According to the 4th IMO GHG study sources of uncertainties can be fuels reported under different categories or placed in both categories (national and international navigation) [10].
The goals of this paper are: to examine correlations between emissions over the observed period, in manoeuvring and hotelling phase, to identify seasonal oscillations, and to give recommendations on how to reduce the gas emissions to improve the quality of living in city port area.

Materials and Methods
The main variables used are the time spent in hotelling/manoeuvring and emissions. Concerning the previously mentioned research problem of emission, the following hypotheses can be defined: (1) there has been no significant increase in the emission of harmful gases in the observed period; (2) there is no change in the trend of CO 2 emission in a period of three months for one year.
The city of Split is the economic and cultural centre of the Croatian region Dalmatia and the second-largest city in Croatia. It is also the greatest passenger port in Croatia, with 2,800,502 passenger arrivals in 2019 [38]. In 2016 the port of Split was the leading port in the ferry, hydrofoil, and fast catamaran traffic on the Adriatic Sea, compared to other Adriatic ports from other countries. The Split-Supetar route is the main passenger transport route in Croatia [39]. The passenger port is located on the south side of the Split peninsula, while on the north side of the Split peninsula is the base of the Croatian Navy and Split Cargo Port.
As a consequence of its position, Split has been a transit city for decades. Recent tourism growth profiled Split as a top destination and tourist record holder by numerous indicators. According to the Split Tourist Board Statistics page, in 2019, Split city had a total of 932,722 tourist arrivals [40]. In comparison to 2018, an 8.15% increase in tourist arrivals has been noted. With its rich historical heritage and favourable climate, it has become a trending location for tourists, whether they are arriving to Split, or just passing by on their way to the Central Dalmatian.
With the increase on tourists, the need for a better connection between the mainland and islands also increases. According to the Port Authority data, the number of ship arrivals is recording growth from 2010 onwards. This growth is graphically presented in Figure 2. In this paper, the port gas emission estimation is calculated for the line ships since these ships make the majority in the number of arrivals and, therefore, impact air pollution.

Emission Estimation for City Port of Split
The approach used in estimating the emissions is consistent with the methodology for quantifying ship emissions in the EMEP/EEA air pollutant emission inventory guidebook. As mentioned in the introduction, the bottom-up approach uses more detailed information than topdown in estimating the emissions.
Although there are several methods for emission calculation, the Tier 3 method was used since it provides insight into the emission of different ships activity like manoeuvring, hotelling and cruising [41].
By comparing several references for estimating emission (ENTEC 2002, ENTEC 2007, and EMEP/EEA 2019 Shipping Tier 3-Ship movement calculation), older emission factors are used due to the fact that in newer literature are only given emission factors for NOX, NMVOC and PM. By comparing newer and older literature, NMVOC factors were the same, and PM and NOX emission factors had slight differences. NOX factors are the same for the ships older than year 2000, making the majority of observed vessels in this paper.
When estimating ship emissions, it is necessary to determine the ship activity. The ship activities are divided into three phases: at the sea, the manoeuvring and the hotelling. The hotelling is the phase when the ships are berthed, while they await their next voyage or cargo load/discharge. The total ship emissions are the sum of the emissions in the abovementioned activities. In this paper, due to the field of interest (the City port of Split), and the data availability, only the manoeuvring and the hotelling phase emissions are calculated. The so-called "At sea" phase is not taken into account in this paper since the area of interest is the City port of Split and emissions in the harbour. Split Port Authority provided the number of ship arrivals and the time spent hotelling as daily based. It must  In this paper, the port gas emission estimation is calculated for the line ships since these ships make the majority in the number of arrivals and, therefore, impact air pollution.

Emission Estimation for City Port of Split
The approach used in estimating the emissions is consistent with the methodology for quantifying ship emissions in the EMEP/EEA air pollutant emission inventory guidebook. As mentioned in the introduction, the bottom-up approach uses more detailed information than top-down in estimating the emissions.
Although there are several methods for emission calculation, the Tier 3 method was used since it provides insight into the emission of different ships activity like manoeuvring, hotelling and cruising [41].
By comparing several references for estimating emission (ENTEC 2002, ENTEC 2007, and EMEP/ EEA 2019 Shipping Tier 3-Ship movement calculation), older emission factors are used due to the fact that in newer literature are only given emission factors for NO X , NMVOC and PM. By comparing newer and older literature, NMVOC factors were the same, and PM and NO X emission factors had slight differences. NO X factors are the same for the ships older than year 2000, making the majority of observed vessels in this paper.
When estimating ship emissions, it is necessary to determine the ship activity. The ship activities are divided into three phases: at the sea, the manoeuvring and the hotelling. The hotelling is the phase when the ships are berthed, while they await their next voyage or cargo load/discharge. The total ship emissions are the sum of the emissions in the abovementioned activities. In this paper, due to the field of interest (the City port of Split), and the data availability, only the manoeuvring and the hotelling phase emissions are calculated. The so-called "At sea" phase is not taken into account in this paper since the area of interest is the City port of Split and emissions in the harbour. Split Port Authority provided the number of ship arrivals and the time spent hotelling as daily based. It must be pointed out that data have been converted into hours for the estimation formula's purpose.
Although ENTEC gives a value of 0.8 h for passenger ships spent in the manoeuvring phase, passenger ships in the port of Split on average of 20 min (0.33 h) spent in manoeuvring phase were provided by Jadrolinija officeholder [42].
The emissions are related to the engine and the fuel type. For each ship installed main or auxiliary engine power data is provided form CRS (Croatian Registry of Shipping) or E-vessel portal, which gives access to electronic services of the Ministry of the Sea, Transport and Infrastructure [43,44]. Hence, in 2017, 35 ships were observed, 2018, 32 ships were observed, and in 2019, 36 ships were observed. The auxiliary engine power is unknown for ten ships, which affects the amount of emissions.
It has to be noted that the main and the auxiliary engines installed in the passenger ships are assumed to be using MDO (marine diesel oil) to comply with the sulphur limits of the Sulphur Content of Marine fuels (SCMF) directive for fuels used by the passenger vessels on the regular services between EU ports.
The engine load factor is defined as the engine's actual power output relative to its Maximum Continuous Rating (MCR).
The Emission Factors are taken from the ENTEC study for estimated pollutants. LF ME is the main engine load factor, LF AE is auxiliary engine load factor, and TO ME is the main engine time of operation during the phase of manoeuvring and hoteling. Their values are shown in Table 2 [44]. The emission factors depend on several factors, such as the main engine type, the auxiliary engine type, and fuel type. Furthermore, ships are divided by the engine speed (slow speed diesel (SSD), medium-speed diesel (MSD), high-speed diesel (HSD), gas turbine, and steam turbine) and the fuel types (RO "Residual Oil" (heavy fuel oil), MDO "Marine Diesel Oil" and MGO "Marine Gas Oil"). There are different NO X emission factors for the main engine depending if the ships are built before or after 2000. The newer engines, which comply with the NO X Technical Code requirements, have roughly 17% lower NO X emissions than the pre-2000 engines [42]. The requirements of the NO X Technical Code have roughly 17% lower NO X emissions than the pre-2000 engines [42]. The emission factors used in this paper are shown in Table 3 [42].  In order to use the bottom up methodology, detailed data is required, including the engine type, the installed power, the hours spent in different phases and the fuel type [42]. The main engine load factors, the auxiliary engine load factors and the main engine time of operation are taken from ENTEC 2002 [45].

Mathematical Backgrounds
The emissions are calculated by multiplying manoeuvring and hotelling time with the sum of the installed main and auxiliary engine power, the load factors for the main and auxiliary engine, the load factors for the main engine, and the operation's main engine time and emission factors. The time spent hotelling was provided in days and converted into hours for the estimation formula's purposes. The same formula is used for estimating air pollution in Ancona harbor [41].
The Emissions for the phase are calculated as follows: where ME is the main engine power (kW); LF ME is the main engine load factor (%); EF ME is the main engine emission factor (g/kWh); TO ME is the main engine time of operation (%); AE is the auxiliary engine (kW); LF AE is the auxiliary engine load factor (%); EF AE is the auxiliary engine emission factor (%); T is the time spent in port (h) or manoeuvring (h); E is emissions (g).
Presented statistical measures would give an insight into relationships between CO 2 emission variables of the Split port. The correlation between the variables is shown with the matrix C A correlation. To get correlations between months a trend or moving average operation needs to be performed on raw data.
The mathematical foundations that are further used in the paper are based on the application from the literature, with the following variables X 2017 , X 2018 , X 2019 being added to the vectors: where X 2017 , X 2018 , X 2019 represents random emission variables of ships date samples obtained from 2017 to 2019.
Standard statistical metrics such as expectation or the average value, standard deviation, and correlation coefficient are used to study random variables. The average value of the random variable x [46].
where E[x] represents the expectation of a random variable x, and N is the number of measurement samples. The following equation can represent the standard deviation of the random variable x: where σ x represents the standard deviation of the random variable x.
The statistical metric used to quantify the similarity and/or dependence among variables, x p , and x v , is the correlation coefficient between random variables. The correlation coefficient can be calculated from the following equation [28]: where r represents the correlation coefficient, N-is the number of measurements while σ xp and σ xv represent the standard deviations of the random variables x p and x v . Furthermore, a model matrix A can be created using the following equation: where X 2017 , X 2018 , X 2019 are vectors of random variables. The correlation matrix C A is defined by the following equation: Using (2) and (6) the correlation matrix C A is defined as follows: In order to perform the smoothing. different moving average algorithms can be used: simple moving average (SMA). a weighted moving average (WMA). an exponential moving average (EMA), and an exponential weighted moving average (EWMA). To perform a simple moving average (SMA) filtering on data following equitation is used [47]: where SMA(n) denotes the moving-average filtering of a vector x. A moving-average filter slides a window of length (WL) along the data and computes averages of the data contained in the (WL) window size.

Case Study
The port of Split is the largest passenger port in the Republic of Croatia, where the arrivals of ship depend on seasonality. The assumption that the larger number of arrivals will produce more emissions, will be tested using the bottom-up method. Because the number of all ship types arrivals is increasing, the paper's main objective is to establish the relationships between the number of arrivals and gas emissions.
In this paper, 34 line ships are observed during 2017, the total number of hours spent in the port is 70,699.97, and the number of calls is 12,330.
In 2018, the number of line ships was 33, the total number of hours spent in the port was 84,519.816, and the number of calls is 13,639.
The total number of line ships was 36 in 2019, the total number of hours spent in the port was 65,908.61, and the number of calls is 14,522.
The engine powers of line ships are different and range from 220 kW-15,015 kW. The calculation for the ship Biokovo is performed by Equation (1), as it shown in Table 4. Biokovo ferry has the main engine power of 1968 kW and the auxiliary engine power of 532 kW. The emission factors for the main engine are taken from Table 3. The main engine emission factor for NO X is 8.8 g/kWh due to the fact that the engine is post-2000. The Auxiliary engine emission factors are calculated from Table 4. The main and the auxiliary engine load factors and the operation's main engine are taken from Table 2. The Port Authority provided the arrival time and the departure time. as well as the days spent in the port.

Results
In this study, the ship emissions are calculated using the activity-based emission estimation for the City port of Split, which is the most significant passenger port in Croatia. The ship emissions are estimated for the ship manoeuvring and hotelling phase. The hotelling phase is responsible for the largest emissions in the port: NO X 90.1%; PM 2.5 78.0% and SO X 88.5% [28]. These percentages can vary depending mainly on time spent hotelling and duration of manoeuvring phase [13].
Total emissions in tons of emitted parameters are shown in Table 5. for each month of the observed years. As seen in Table 5. the highest amount of emissions in 2018 was calculated during ship manoeuvring and hotelling operations for all investigated pollutants.  From recent literature, the traffic data for passenger ships other than a cruise, together with annual emissions (tons per year) of NO X and PM 10 for the ports Barcelona, Hong Kong, Copenhagen, Venice, Elsinore, St Petersburg, Las Palmas, Genoa, and Marseille are given [13]. These ports have no similarities in size or number of ship arrivals/departures, but the aim was to take into account emissions in urban cities that have ports within the city, such as Split. The minimum NO X emissions are 20 tons per year in the Hong Kong and 1300 tons per year calculated in Marseille port [13]. For PM 10 emissions in ports, the range is between 1 ton per year in the Hong Kong port and 80 tons per year in Marseille port [13]. Exhaust gas emissions from ships are calculated for Izmir Port using the ship activity-based methodology. Total emissions from ships in the port is estimated as 1923 tons per year for NO X , 1405 tons per year for SO 2 , 82,753 tons per year for CO 2 , 1 ton per year for HC, and 165 tons per year for PM in the year 2007 [30].
The Emissions for each pollutant, divided into manoeuvring and hotelling phase, are shown in Figures 3-7. The emissions are divided into years and months and are represented in tons (t). Blue plots represent the hotelling phase (H) and the orange plots represent the manoeuvring phase (M). If side by side comparison is done, the emissions were higher for all pollutants during the hotelling phase than during manoeuvring phase, due to longer time spent in port. Figure 8 the hotelling phase contributes to most SO 2 emissions with several peaks during the year. On the contrary, the emission trough the manoeuvring phase shows seasonality each year, and it is always relatively smaller than emission during the hotelling phase. The total annual SO 2 emission varies from 118 to 127 tons per year, which is way lower than 1405 tons per year in Izmir port [30].

As shown in
Furthermore, NO X emission varies between 215 tons per year and 235 tons per year (Figure 3), which is almost equal to Barcelona NO X emissions. Ports Hong Kong. Copenhagen, Venice, and Elsinore have lower NO X emissions in ports. Ports St Petersburg, Las Palmas, Genoa and Marseille, have higher NO X emissions. Croatian port of Zadar has total annual cruise ships emissions 310.23 tons per year for NOx, and 9.62 tons per year for PM [35]. Port of Split has lower NO X emissions and higher PM emissions compared to the port of Zadar. tons per year for NOx, and 9.62 tons per year for PM [35]. Port of Split has lower NOX emissions and higher PM emissions compared to the port of Zadar. The study of PM10 emission in Ancona's port reveals that most emissions (70%) happen during the hotelling phase [41]. These emissions correlate with the PM10 emissions from the port of Split (Figure 4). During the hoteling PM emission represents 69% in 2017, 68% in 2018 and 65% in 2019 of total PM emissions in Venice port, Ports Hong Kong, Elsinore, St Petersburg, and Las Palmas have lower PM emissions. Contrary, ports with higher PM10 emissions are Barcelona, Genoa, Copenhagen, and Marseille [13]. From Figure 5, the total annual NMVOC emission varies from 12 to 13 tons per year. In manoeuvring phase emissions are slightly higher during the tourist season period than hotelling phase emissions. This can be explained that in the tourist season period, the number of ship lines has been increased, and therefore, ships spend less time in the ports. As a consequence, the hotelling time is reduced and, therefore, emissions. Additionally, more port arrivals/departures increase The study of PM 10 emission in Ancona's port reveals that most emissions (70%) happen during the hotelling phase [41]. These emissions correlate with the PM 10 emissions from the port of Split (Figure 4) tons per year for NOx, and 9.62 tons per year for PM [35]. Port of Split has lower NOX emissions and higher PM emissions compared to the port of Zadar. The study of PM10 emission in Ancona's port reveals that most emissions (70%) happen during the hotelling phase [41]. These emissions correlate with the PM10 emissions from the port of Split (Figure 4). During the hoteling PM emission represents 69% in 2017, 68% in 2018 and 65% in 2019 of total PM emissions in Venice port, Ports Hong Kong, Elsinore, St Petersburg, and Las Palmas have lower PM emissions. Contrary, ports with higher PM10 emissions are Barcelona, Genoa, Copenhagen, and Marseille [13]. From Figure 5, the total annual NMVOC emission varies from 12 to 13 tons per year. In manoeuvring phase emissions are slightly higher during the tourist season period than hotelling phase emissions. This can be explained that in the tourist season period, the number of ship lines has been increased, and therefore, ships spend less time in the ports. As a consequence, the hotelling time is reduced and, therefore, emissions. Additionally, more port arrivals/departures increase From Figure 5, the total annual NMVOC emission varies from 12 to 13 tons per year. In manoeuvring phase emissions are slightly higher during the tourist season period than hotelling phase emissions. This can be explained that in the tourist season period, the number of ship lines has been increased, and therefore, ships spend less time in the ports. As a consequence, the hotelling time is reduced and, therefore, emissions. Additionally, more port arrivals/departures increase manoeuvring time and thus manoeuvring phase emissions. CO2 emissions contribute majorly to total emissions. From Figure 6, the total annual CO2 emission varies from 12.5 to 13.5 kt. From study [48] and monthly results, there is a difference between the in-season and off-season periods. All emissions slightly increase starting from May and then decrease from October towards the end of the year. The highest manoeuvring emissions are in July and August for all observed years. Such emission trend corresponds to the Croatian shipping company "Jadrolinija" which has the largest number of ship arrivals in the Split city port. The number of ferry lines increases since the end of May, with the peak arrivals in July, August, and September. From the beginning of October, the number of lines decreases. CO 2 emissions contribute majorly to total emissions. From Figure 6, the total annual CO 2 emission varies from 12.5 to 13.5 kt. From study [48] and monthly results, there is a difference between the in-season and off-season periods. All emissions slightly increase starting from May and then decrease from October towards the end of the year. The highest manoeuvring emissions are in July and August for all observed years. Such emission trend corresponds to the Croatian shipping company "Jadrolinija" which has the largest number of ship arrivals in the Split city port. The number of ferry lines increases since the end of May, with the peak arrivals in July, August, and September. From the beginning of October, the number of lines decreases. CO2 emissions contribute majorly to total emissions. From Figure 6, the total annual CO2 emission varies from 12.5 to 13.5 kt. From study [48] and monthly results, there is a difference between the in-season and off-season periods. All emissions slightly increase starting from May and then decrease from October towards the end of the year. The highest manoeuvring emissions are in July and August for all observed years. Such emission trend corresponds to the Croatian shipping company "Jadrolinija" which has the largest number of ship arrivals in the Split city port. The number of ferry lines increases since the end of May, with the peak arrivals in July, August, and September. From the beginning of October, the number of lines decreases. Unlike manoeuvring emissions, hotelling emissions are not showing seasonality. When hotelling emissions are estimated monthly, there are no significant differences between the in-season and off-season periods. All emissions increase between May and October and emission increase correlates with more time spent in port. The highest hotelling phase emissions in 2017 were in February. It is well known that February is a month when the ship spends almost all the time in the port. Consequently, when hotelling phase emission data are observed, it is evident that the months with a reduction in emissions are the months with a time reduction in port retention.

Analysis of Data
In the next section, a trend and statistical analysis will be performed on CO 2    To check the relationship between independent (month) variable and dependent variables (CO2 emissions for 2017, 2018, and 2019), an F-test is performed, as shown in Table 6.    To check the relationship between independent (month) variable and dependent variables (CO2 emissions for 2017, 2018, and 2019), an F-test is performed, as shown in Table 6.  To check the relationship between independent (month) variable and dependent variables (CO 2 emissions for 2017, 2018, and 2019), an F-test is performed, as shown in Table 6. From Table 6 it can be seen that all p values are significant, which implies that there are no hidden correlations between independent variable, i.e., month, and all dependent variables, i.e., CO 2 emissions for 2017, 2018, and 2019.
Next, trend analysis on the dependent variable, i.e., CO 2 emissions, has been performed using regression. The trend analysis is performed to check hypothesis two. Figure 9 shows a graphical presentation of trend analysis of CO 2 emissions through observed years, using Equation (9) and regressions. From Table 6 it can be seen that all p values are significant, which implies that there are no hidden correlations between independent variable, i.e., month, and all dependent variables, i.e., CO2 emissions for 2017, 2018, and 2019.
Next, trend analysis on the dependent variable, i.e., CO2 emissions, has been performed using regression. The trend analysis is performed to check hypothesis two. Figure 9 shows a graphical presentation of trend analysis of CO2 emissions through observed years, using Equation (9) and regressions. From Figure 9, it can be observed that blue circles show the values after performing moving average filtering (Equation (9)) of CO2 emissions for observed years. A period of three months is chosen as windows length (WL) for the moving average algorithm. The choice for WL = 3 is derived from the fact that in Republic Croatia, we observe three (3) months, which corresponds with in-season and off-season periods and to eliminate extreme values. Furthermore, the red line shows the trend (regression) of CO2 emissions for the observed years. From the trend of CO2 emissions in 2017, it can be observed that the emissions are constant, although the slope constant has a small negative value (−0.0085). Additionally, the same conclusion could be drawn for 2018. Contrary, for 2019, it can be seen that the slope constant has a small positive value (+0.0068), but since the constant is negligible, From Figure 9, it can be observed that blue circles show the values after performing moving average filtering (Equation (9)) of CO 2 emissions for observed years. A period of three months is chosen as windows length (WL) for the moving average algorithm. The choice for WL = 3 is derived from the fact that in Republic Croatia, we observe three (3) months, which corresponds with in-season and off-season periods and to eliminate extreme values. Furthermore, the red line shows the trend (regression) of CO 2 emissions for the observed years. From the trend of CO 2 emissions in 2017, it can be observed that the emissions are constant, although the slope constant has a small negative value (−0.0085). Additionally, the same conclusion could be drawn for 2018. Contrary, for 2019, it can be seen that the slope constant has a small positive value (+0.0068), but since the constant is negligible, it can be concluded that CO 2 emission is constant over the year. Overall, the average emission of CO 2 stays constant, although the number of vessels during 2019 is increased.
To get further insight into CO 2 emission, for observed years, a correlation analysis is performed (Equations (6)-(9)), as shown in Table 7.  Table 7, it can be observed there are strong correlations between all variables. For example, the correlation coefficient between the CO 2 emission from 2017 and CO 2 emission from 2018 is 0.906. Additionally, the "weakest" correlation is between the CO 2 emission from 2017 and CO 2 emission from 2019 (0.670). Further, ANOVA analysis is performed to check the hidden correlations within variables and see. The ANOVA analysis (F = 0.521, Fcrit = 3.285, p = 0.598) shows there are no variations between variables, which implies that all variables come from the same population. That can be explained by the fact that the Split city port operates within full capacity.
It is obvious that some organizational and technological solutions have to be applied to reduce harmful gas emissions, such as better voyage planning (more vessels in less time, especially during the season), usage of more environmentally friendly fuel, etc.
Besides these measures, it is important to find new improvements and solutions such as better cooperation with agents, cargo suppliers. i.e., the creation of integrated management. In terms of technological options, one should be open to renewable sources (fuel cells, wind energy, solar, . . . ) if possible, use environmentally friendly materials of more environmentally efficient technology and find more efficient solutions when building a ship.

Conclusions
This paper investigates the ship emissions of pollutants and gases developed by line ships during manoeuvring and hotelling phases in cities that have ports in their centres. The analysis was conducted on the example of the Port of Split.
In the hotelling phase, when the ship is on the berth, the auxiliary engines have higher contribution to emissions than during the "at sea" phase. During manoeuvring phase propulsion engines are operating at low loads and auxiliary engine loads are at their highest load for additional onboard equipment such as thrusters. At berth propulsion engines are usually off and auxiliary engine loads are high during discharging cargo, cars, etc. it is presumed that vessels have their auxiliary engines on during the whole of the time spent in the port.
Through the manoeuvring, it can be observed the emissions show seasonality. That seasonality corresponds to increasing the ferry lines number from the end of May, with the largest arrivals in July, August, and September. Hotelling phase emissions are not shown seasonal patterns, but they depend on ship's number and port retention.
Comparing all emissions, carbon dioxide (CO 2 ), nitrogen oxides (NO X ), volatile organic compounds (VOC), sulphur dioxide (SO 2 ), and particulate matter (PM) it can be seen they are linearly correlated. From the data, it is evident that only CO 2 emission had a major impact on pollution.
Trend analysis for CO 2 emissions in 2017, 2018, and 2019 shows that the average emissions are constant, despite the number of vessels and staying time in the port.
Correlation analysis shows there are strong correlations between the CO 2 emission variables. Additionally, the ANOVA test confirms the findings. That implies all variables come from the same population, which points out that the Split city port operates within full capacity.
The conducted analysis shows that the variables observed during the manoeuvring and hoteling phases in the city port must be taken into account when developing city strategies to make maritime transport (line ships) more efficient, which could contribute in the future to reduce emissions.