Stochastic Air Quality Dispersion Model for Deﬁning Queuing Ships Seaport Location

: This work develops a stochastic air quality dispersion to predict the pollution concentration originating from ships queuing in a seaport. The Gaussian dispersion model for ﬁve ships operating in the Black Sea queuing in the front of the port of Varna as sources of gas emission of NO x , SO x and PM 10 is used to deﬁne the air pollution concentration at receptors (crowded areas of the port and other reference points) and consequently the distance to the seaport queuing location. Uncertainties, which are inherent in the input data and mathematical model, are accounted for to estimate the propagating uncertainties of the emission concentration at the receptors accounting for the wind speed, horizontal and vertical dispersion parameters as a function of the geographical location of the emission sources (ships), effective emission height and weather conditions. The estimated uncertainties of the air quality prediction are of signiﬁcant importance for the decision-making on the regulatory purposes, and the probability of exceeding the threshold limits needs to be quantiﬁed. The most expected value and the probability of exceeding the acceptable limits of pollution concentration are deﬁned by employing the ﬁrst-order reliability method. The target reliability level is deﬁned as the failure cause and mode used for identifying the safety calibration factors that may be employed for deﬁning the most suitable location of the ship queuing seaport. Several conclusions about the applicability of the developed stochastic model and its use for regulatory purposes are also provided.


Introduction
Ports are one of the most important parts of the hinterland links and short sea shipping since international trade is mostly transported using ships [1]. Sea cargo transportation is constantly expanding resulting in constant growth in the emissions of nitrogen oxides NO x , sulfur oxides SO x , particulate matter PM and others originating from the diesel combustion in ocean-going vessels [2,3], and their control becomes dominant in keeping the green environmental conditions [4]. The EU Thematic Strategy on Air Pollution was identified in [5] and confirmed in [6] the importance for the health and environment the reduction of sulfur dioxide SO 2 , NO x and PM from ships.
The predominant fuel oil in shipping [7] is high-sulfur fuel oil (HSFO) and it is still used in many low-to-medium speed engines [8]. HSFO was the dominant fuel in international shipping about 79% of the total fuel consumption by energy content in 2018 [9]. The share of bunker sales reached 90% in the first quarterof 2019 in Singapore, which is the world's largest bunkering hub [10]. After the quick adjustment to the global sulfur cap (from 1 January 2020), the share of HSFO declined to about 17% in the first quarter of 2020 in Singapore. The decline later stopped and the steady rise in scrubber-fitted ships has supported the demand for HSFO (23% in the first quarter of 2021 in Singapore) and this seems to remain until new solutions and future fuels are widely introduced.
The artificial neural networks methodologies have been also applied recently in analyzing the impact of environmental pollution [26], associative memories [27], and factor analysis [28] among others.
The study here uses the Gaussian dispersion model in defining the pollution concentration at receptor locations (crowded port areas and other reference points). Uncertainties of the emission concentration at the receptors account for the wind speed, horizontal and vertical dispersion parameters as a function of the geographical location of the emission sources (ships), effective emission height, and weather conditions. A target level is defined related to failure cause and mode is used for identifying safety calibration factors in defining the most suitable location and distance from the receptors (terminal port) and the source of pollution represented by the ships queuing at the seaport.

Pollutant Emissions Abatement Control
The abatement cost for pollutant emissions is established based on the Gaussian dispersion model and it is associated with the applied technology for removing the air pollutants. For each pollutant, the total cost, conditional on the abatement targets is estimated using the semi-infinite programming model based on the unabated emissions and the unit cost of the cleaning technology.
The relationship between the concentrations emission rate accounting for meteorological conditions can be defined by the vertical and the crosswind Gaussian dispersion model [29][30][31], estimating the air pollution concentration for one specific period: where Q m is the emission rate, σ m,y , σ m,z are the dispersion coefficients in crosswind and vertical directions, x and y are the longitudinal and transverse distance from the source, u m is the wind speed, and H is the effective stack height, showing that the concentration in each point of the space is proportional to the emission rate Q m and inversely proportional to the wind speed u m (see Figure 1). The dispersion coefficients σ m,y , σ m,z are related to a specific receptor position and weather conditions.
solutions and in the case of air pollutants in a continuous value between 0 and 1 demonstrating a satisfactory solution, where for assessing air quality environmental level index has been introduced [23][24][25]. Most of the studies do not analyze statistically the impact of sampling deviations. The artificial neural networks methodologies have been also applied recently in analyzing the impact of environmental pollution [26], associative memories [27], and factor analysis [28] among others.
The study here uses the Gaussian dispersion model in defining the pollution concentration at receptor locations (crowded port areas and other reference points). Uncertainties of the emission concentration at the receptors account for the wind speed, horizontal and vertical dispersion parameters as a function of the geographical location of the emission sources (ships), effective emission height, and weather conditions. A target level is defined related to failure cause and mode is used for identifying safety calibration factors in defining the most suitable location and distance from the receptors (terminal port) and the source of pollution represented by the ships queuing at the seaport.

Pollutant Emissions Abatement Control
The abatement cost for pollutant emissions is established based on the Gaussian dispersion model and it is associated with the applied technology for removing the air pollutants. For each pollutant, the total cost, conditional on the abatement targets is estimated using the semi-infinite programming model based on the unabated emissions and the unit cost of the cleaning technology.
The relationship between the concentrations emission rate accounting for meteorological conditions can be defined by the vertical and the crosswind Gaussian dispersion model [29][30][31], estimating the air pollution concentration for one specific period: , , where is the emission rate, , , , are the dispersion coefficients in crosswind and vertical directions, and are the longitudinal and transverse distance from the source, is the wind speed, and is the effective stack height, showing that the concentration in each point of the space is proportional to the emission rate and inversely proportional to the wind speed (see Figure 1). The dispersion coefficients , , , are related to a specific receptor position and weather conditions.  The average concentration over a long period can be estimated by assuming that the period is composed of a sum of a series of stationary conditions, not accounting for the possible correlation between the sequence conditions, where assuming that the concentration C m is estimated for the meteorological condition m in one specific period, which may occur with a probability p m of occurrence, where the average concentration is defined as: where C b is the background concentration.
If the emission rate of the source Q m is a constant for a relatively long period, and the probability of occurrence of any possible weather condition is known p m , the average pollution at the receptor can be estimated as: where α m is constant depending only on the weather condition for a distance from the receptor to the source.
In general, the air pollution is dispersed in a three-dimensional space Ω generated from n sources (ships), where at each of the m point receptor locations x i , y i , z i the concentration on the contaminant is: where C j x ij , y ij , z ij is the concentration provided by the jth source and ith receptor and based on the meteorological and geographic conditions and the geographical distance between the resources and receptors x ij , y ij , z ij are calculated as: where θ represents the angle of the wind direction, x S,i , y S,i are the coordinates of the point source and x R,j , y R,j are the coordinates of the receptor. If w j is the reduction factor of the air contamination from the jth source, when j = 1 . . . n, and w = (w 1 , . . . , w n ) ∈ [0,1] n and assuming that the total cost G(w) is a linear function of w j defined as: the problem of minimizing the air pollution control costs into the space Ω can be defined by satisfying the following conditions: which is constrained to where C MAC,i is the maximum allowable concentration at the ith receptor location. The result of minimizing the air pollution cost will be the maximum emission rate that the sources (ships) may have without crossing the threshold of the intolerable average air pollution concentration at the location of receptors (ports).
In general, it can consider the existence of finite sources (ships), and an infinite number of constraints, as many as the number of receptors leading to a semi-infinite programming problem. This problem allows the compliance of the constraints observed at each receptor point. The optimal values of the w j indicate the contamination reduction of each source at the least cost [32]. Some discussions about the cost associated with the air pollutants at the source can be found in [33]. In the case of cyclones, momentum separators, filters, scrubbers, electrostatic precipitators, the cost is about 0.71 Euro per kg removed particulates. For scrubber, catalytic conversion the cost is about 0.55 Euro per kg removed SO 2 . For scrubber, catalytic conversion the cost is about 0.55 Euro per kg removed NO x .
The variation of the pollution concentration at any distance from the source depends on the variation of the input parameters and the sensitivity of the output estimate on the sensitivities of the input.
A study about the requirements and cargo transportation demand in the Black Sea as part of a multimodal transportation frame was performed, estimating the potential need for a shipping fleet of multipurpose ships [34,35]. Five of these ships, transporting cargoes and queuing very closely in the front of the port of Varna are considered as a one-point source of pollution, n = 1.
The emission rate generated by the source of pollution is 4347 NO x kg/day, 2344 SO x kg/day and 276 PM 10 kg/day located at x = 0.00 m, y = 0.00 m, y = 0.00 m, m = 1, C b = 0, as indicated in Equation (1). The ocean sea is assumed as a flat surface, and the effective stack height is 12 m. The wind speed is 2 m/s, oriented from the source to the receptors [36]. The environmental stability parameters are based on night-cloudiness used for calculating the horizontal and vertical dispersion parameters, σ y and σ z .
The original pollution concentration generated by the queuing ships cross the acceptable limit and the receptor locations and cleaning devices were employed to reduce it, as shown in Figure 2. The Maximum Allowable Concentrations, C MAC,i , considered here are: C MAC,NO x = 30 µg/m 3 , C MAC,SO x = 20 µg/m 3 , C MAC,PM 10 = 20 µg/m 3 more information can be seen in [37].
In general, it can consider the existence of finite sources (ships), and an infinite number of constraints, as many as the number of receptors leading to a semi-infinite programming problem. This problem allows the compliance of the constraints observed at each receptor point. The optimal values of the indicate the contamination reduction of each source at the least cost [32].
Some discussions about the cost associated with the air pollutants at the source can be found in [33]. In the case of cyclones, momentum separators, filters, scrubbers, electrostatic precipitators, the cost is about 0.71 Euro per kg removed particulates. For scrubber, catalytic conversion the cost is about 0.55 Euro per kg removed . For scrubber, catalytic conversion the cost is about 0.55 Euro per kg removed . The variation of the pollution concentration at any distance from the source depends on the variation of the input parameters and the sensitivity of the output estimate on the sensitivities of the input.
A study about the requirements and cargo transportation demand in the Black Sea as part of a multimodal transportation frame was performed, estimating the potential need for a shipping fleet of multipurpose ships [34,35]. Five of these ships, transporting cargoes and queuing very closely in the front of the port of Varna are considered as a onepoint source of pollution, 1. The emission rate generated by the source of pollution is 4347 kg/day, 2344 kg/day and 276 kg/day located at 0.00 m, 0.00 m, 0.00 m, 1, 0, as indicated in Equation (1). The ocean sea is assumed as a flat surface, and the effective stack height is 12 m. The wind speed is 2 m/s, oriented from the source to the receptors [36]. The environmental stability parameters are based on night-cloudiness used for calculating the horizontal and vertical dispersion parameters, and . The original pollution concentration generated by the queuing ships cross the acceptable limit and the receptor locations and cleaning devices were employed to reduce it, as shown in Figure   The problem analyzed here can be extended to define what is the maximum emission rate that may guarantee a tolerable air pollution concentration with a given probability in a determined period, what is the minimum abatement cost in the presence of several air pollution sources (ships) that may guarantee a tolerable air pollution concentration at The problem analyzed here can be extended to define what is the maximum emission rate that may guarantee a tolerable air pollution concentration with a given probability in a determined period, what is the minimum abatement cost in the presence of several air pollution sources (ships) that may guarantee a tolerable air pollution concentration at receptors located at the port. It can also be solved an inverse engineering problem related to how large are the pollution emissions at given locations that determine the existing concentration at a specific source location sharing and to find the optimal abatement rates for several ships to guarantee a tolerable air pollution concentration at a series of receptors. Finally, is to find the best position of a seaport for queuing ships with emission rates at a distance from a port for cargo operation, which is analyzed in the next sessions.

Sensitivity and Uncertainty Analysis
Two types of sensitivity analysis are usually distinguished: local and global sensitivity. The local sensitivity analysis focuses on the local impact of the factors on the model outputs, and it is carried out by computing partial derivatives of the output concerning the input factors. With this kind of approach, the factors are allowed to vary within small intervals around the nominal values. These intervals are not related to the uncertainty in the factor values.
The global sensitivity analysis can be used to identify the importance of any individual factor concerning the uncertainty in the output. The global sensitivity analysis involves the distributions of the input factors and computes the output variation. The computation involves the sensitivity index for each input variable. These indices are estimated to vary the factors over their whole uncertainty ranges. The global sensitivity analysis allows identifying the variables that need an accurate measure or estimation.
First, a local sensitivity analysis is performed to evaluate the Gaussian dispersion model output sensitivities caused by the variation of the input parameters identifying the most significant ones. The specified important parameters will require the greatest accuracy and precision in their quantification in reducing the uncertainty in model predictions.
The local sensitivity of the governing variables is calculated by: where x= {Q, u, σ y , σ z and H}. The local sensitivities estimates are shown in Figure 3, conditional to y = z = 0. In defining the air pollution concentration, the most sensitive parameters vary with the distance between the source of pollution (ships) and the receptor (port), where the most sensitive factors are, u, σ z , H, σ y , Q and as can be seen from Figure 3.
concentration at a specific source location sharing and to find the optimal abatement rates for several ships to guarantee a tolerable air pollution concentration at a series of receptors. Finally, is to find the best position of a seaport for queuing ships with emission rates at a distance from a port for cargo operation, which is analyzed in the next sessions.

Sensitivity and Uncertainty Analysis
Two types of sensitivity analysis are usually distinguished: local and global sensitivity. The local sensitivity analysis focuses on the local impact of the factors on the model outputs, and it is carried out by computing partial derivatives of the output concerning the input factors. With this kind of approach, the factors are allowed to vary within small intervals around the nominal values. These intervals are not related to the uncertainty in the factor values.
The global sensitivity analysis can be used to identify the importance of any individual factor concerning the uncertainty in the output. The global sensitivity analysis involves the distributions of the input factors and computes the output variation. The computation involves the sensitivity index for each input variable. These indices are estimated to vary the factors over their whole uncertainty ranges. The global sensitivity analysis allows identifying the variables that need an accurate measure or estimation.
First, a local sensitivity analysis is performed to evaluate the Gaussian dispersion model output sensitivities caused by the variation of the input parameters identifying the most significant ones. The specified important parameters will require the greatest accuracy and precision in their quantification in reducing the uncertainty in model predictions.
The local sensitivity of the governing variables is calculated by: where x= { , , , and }. The local sensitivities estimates are shown in Figure 3, conditional to 0. In defining the air pollution concentration, the most sensitive parameters vary with the distance between the source of pollution (ships) and the receptor (port), where the most sensitive factors are, , , , , and as can be seen from   The global Gaussian dispersion model uncertainties are evaluated using the standard deviation around the mean value of the output air pollution concentration employing the Taylor series expansion. Several assumptions are made where the estimations of the parameters included in Equation (1) are uncertain including the wind speed, u, horizontal dispersion parameter, σ y , vertical dispersion parameter, σ z , emission rate, Q, and effective emission height, H. The error related to u, σ y , σ z , Q and H are random, uncorrelated, and are normally distributed with mean values of E u , E σ y , E σ z , E H , E Q and standard deviation equal to S u , S σ y , S σ z , S H , S Q , respectively.
The second-order approximation of the variance based on the Taylor series expansion is used to describe the uncertainty in the modeled dispersion. A non-linear relationship between the input variables and the output estimate exists. It has to be pointed out that the dispersion parameters, σ y (x), σ z (x) are calculated as functions of downwind distance x and y = z = 0 for discrete stability classes as defined by Pasquill (1961).
The uncertainty related to C(x, 0, 0) when the input parameters are in a non-linear relationship to the output, the higher-order terms of the Taylor series equation are included, and the second-order approximation of variance is given [38] as: where ∂C ∂x i and ∂ 2 C ∂x i 2 are the first and second-order sensitivities of the input parameters. Here for the concentration rate, an uncertainty of 10 percent is assumed. The wind speed and direction are measured at a 10 m height with the wind vanes and anemometers within a tolerance of 0.25 m/s, 15 percent uncertainty is assumed.
The original data on which the Pasquill [30] dispersion curves are based, it is difficult to estimate the measurement uncertainties for σ y , σ z . For the base case uncertainties, a conservative value of 20 percent is used to represent the measurement portion of the total uncertainty.
The effective emission height, H (height of plume above ground), depends on the release height and plume rise, expressed as a function of wind speed, ambient and stack gas temperatures, and exit velocity. For simplicity, H is treated as a measured variable and assigned a base case uncertainty to 15 percent of its value. The stochastic model of the basic variables considered in this study is presented in Table 1. A Normal distributed random variable B xi is introduced to take into account the model uncertainty on Q, u, σ y , σ z , H. The uncertainty calculations can be performed using the total uncertainty on any individual variable estimate and model uncertainty represented by the random variable with a mean value of B x = 1 and coefficient of variation COV Bx = 0.05 for B Q , B u , B σ y , B σ z and B H are determined as: The emission rate E Q = 4347 NO x kg/day, E SO x = 2344 SO 2 kg/day and E PM 10 = 276 PM 10 kg/day, E u = 2 m/s, E H = 12, m, and E σ y E σ z as a function of the distance and weather conditions as defined by Pasquill [30]. The propagating second-order global uncertainties are shown in  The emission rate = 4347 kg/day, = 2344 kg/day and kg/day, = 2 m/sec, = 12, m, and as a function of the distan weather conditions as defined by Pasquill [30]. The propagating second-order glo certainties are shown in Figure 4.

Risk-Based Decision Making
A distance analysis is performed here in defining the best location of the seap queuing ships, identified as a source of air pollution. The analysis involves the p cargo operation (receptor), the seaport for queuing (source), the distance betwe cargo operation and queuing ports, and the metric defined either as a distance o between the ports for cargo operation and queuing.
The distances between the cargo operation and queuing ports can be calculated the shortest route between the ports. The distance analysis includes the pull prob minimize the distance between the ports, the capture problem that relates to the c sociated with the distance of the queuing ship (port services and cargo transpor push problem that considers pollution issues that maximize the distances betwe ports.
The final objective is to find a distance that minimizes the sum of the total ri associated cost. The push and pull distance strategies will be employed to estim risk and associated cost as a function of the distances between the cargo-operatio and queuing-seaport. The push strategy will minimize the emission exposure by m izing the distances between the ports such that the pollutant exposure to the cargo ation port is less than the established threshold for air pollution. The pull strategy f on finding a distance, where the total cost of queuing ships losses, and cargo deliv minimized. The pull strategy will tend to reduce the distance between the ports.
A decision about the location of the seaport for queuing ships can be made b lyzing the pollution concentration generated by the ships and resulting environ damage in minimizing the total expected risk. The risk associated with the ecologi lapse is estimated as the probability of collapse, defined as the pollution concen crosses acceptable limits, times the consequences in a monetary term defined as: where | , is the probability of failure, is the cost associated with the failure is tor of governing parameters, is the Beta reliability index, and is the period where th estimated.

Risk-Based Decision Making
A distance analysis is performed here in defining the best location of the seaport for queuing ships, identified as a source of air pollution. The analysis involves the port for cargo operation (receptor), the seaport for queuing (source), the distance between the cargo operation and queuing ports, and the metric defined either as a distance or time between the ports for cargo operation and queuing.
The distances between the cargo operation and queuing ports can be calculated using the shortest route between the ports. The distance analysis includes the pull problem to minimize the distance between the ports, the capture problem that relates to the cost associated with the distance of the queuing ship (port services and cargo transportation), push problem that considers pollution issues that maximize the distances between the ports.
The final objective is to find a distance that minimizes the sum of the total risk and associated cost. The push and pull distance strategies will be employed to estimate the risk and associated cost as a function of the distances between the cargo-operation port and queuing-seaport. The push strategy will minimize the emission exposure by maximizing the distances between the ports such that the pollutant exposure to the cargo-operation port is less than the established threshold for air pollution. The pull strategy focuses on finding a distance, where the total cost of queuing ships losses, and cargo delivery are minimized. The pull strategy will tend to reduce the distance between the ports.
A decision about the location of the seaport for queuing ships can be made by analyzing the pollution concentration generated by the ships and resulting environmental damage in minimizing the total expected risk. The risk associated with the ecological collapse is estimated as the probability of collapse, defined as the pollution concentration crosses acceptable limits, times the consequences in a monetary term defined as: where P f (t n |x, β ) is the probability of failure, G is the cost associated with the failure x is the vector of governing parameters, β is the Beta reliability index, and t n is the period where the risk is estimated. The reliability, R or the probability of failure, P f = 1 − R, analysis is incorporated into the decision-making procedure, where the statistical nature of the governing parameters and decision-making problem is defined in an objective function, and the probabilistic constraints specify the required reliability target level.
The reliability analysis performed here uses the First Order Reliability Method, FORM [39][40][41] that identifies a set of basic random variables that influence the failure mode or the limit-state under consideration. The failure probability, P f concerning a single failure mode can formally be defined as: All random variables are considered non-correlated Normal distributed. The limits state function, g(x) is minimized, and it is subjected to constraints, where x is the vector of the random variables.
The limit state function is defined as: where The Beta index, as a function of the distance between the source and receptor, is presented in Figure 5. The minimum Beta index is achieved at 1794 m and from the source and from that point ahead the Beta index group up since the pollution concentration reduces.
FORM [39][40][41] that identifies a set of basic random variables that influence the failure mode or the limit-state under consideration. The failure probability, concerning a single failure mode can formally be defined as: All random variables are considered non-correlated Normal distributed. The limits state function, is minimized, and it is subjected to constraints, where is the vector of the random variables.
The limit state function is defined as: The Beta index, as a function of the distance between the source and receptor, is presented in Figure 5. The minimum Beta index is achieved at 1794 m and from the source and from that point ahead the Beta index group up since the pollution concentration reduces.
The importance of the contribution of each random variable to the limit state function * at the design point is assessed by the sensitivity factors, which are determined as: 1 ∑ * * (18) and are shown in Figure 6.  The importance of the contribution of each random variable to the limit state function g * (x) at the design point is assessed by the sensitivity factors, which are determined as: and are shown in Figure 6. The uncertainties associated with the pollution and the potential ecological damage resulting from the air pollution gases, which is expressed as a probability of failure, P f is accounted for in the risk assessment.
The uncertainties associated with the pollution and the potential ecological damage resulting from the air pollution gases, which is expressed as a probability of failure, is accounted for in the risk assessment.
The optimal safety level is defined by employing a cost-benefit analysis. The goal is to establish an acceptable safety and reliability level by using the risk control option in determining the acceptable threshold of the air pollution concentration. The cost-benefit analysis is performed based on the total expected risk, , which is a product of the probability of failure and consequence cost [35,42]: where is the risk associated with the environmental failure and its consequence costs and is the cost of the implemented safety measures in keeping the environment safe.
A decision can be made analyzing the cost associated with the pollution abatement, ecological damage, and port service as a function of the distance where the ships are queuing in minimizing the total expected risk. The risk associated with the environmental collapse accounts for the probability of failure and the cost of consequences is defined as: where 0 is the probability of the air pollution concentration crossing the threshold limit and is the impact, which defines the consequence cost of the air concentration crossing the threshold limit.
The damage of the ecosystem estimated in a monetary unit, , due to the air pollution gases generated by the ships engine propulsion and auxiliary systems, once the abundant concentration is abated (using scrubbers or other cleaning systems) to satisfy the minimum allowable concentration, it is defined as: A recent study about the contribution of different pollutants to our health and economy reveals that , and contributes as 15%, 18%, and 50% respectively of the total annual damage per μg m ⁄ , assumed here as 2000 €/μg m ⁄ .
The cost related to the enhancement of environmental safety by taking measures is the target reliability level that depends on different factors and may vary from one The optimal safety level is defined by employing a cost-benefit analysis. The goal is to establish an acceptable safety and reliability level by using the risk control option in determining the acceptable threshold of the air pollution concentration. The cost-benefit analysis is performed based on the total expected risk, Risk t (x), which is a product of the probability of failure and consequence cost [35,42]: where Risk f (x) is the risk associated with the environmental failure and its consequence costs and Risk m (x) is the cost of the implemented safety measures in keeping the environment safe. A decision can be made analyzing the cost associated with the pollution abatement, ecological damage, and port service as a function of the distance where the ships are queuing in minimizing the total expected risk. The risk associated with the environmental collapse accounts for the probability of failure and the cost of consequences is defined as: (20) where P f {P[g(x) ≤ 0]} is the probability of the air pollution concentration crossing the threshold limit and G f (x) is the impact, which defines the consequence cost of the air concentration crossing the threshold limit. The damage of the ecosystem estimated in a monetary unit, G eco (x), due to the air pollution gases generated by the ships engine propulsion and auxiliary systems, once the abundant concentration is abated (using scrubbers or other cleaning systems) to satisfy the minimum allowable concentration, it is defined as: where G eco, k denotes the damage cost of the pollutant k, Euro per µg/m 3 and Q a,k (x|x, 0, 0 ) is the emission rate after the abatement.
A recent study about the contribution of different pollutants to our health and economy reveals that SO x , NO x and PM 10 contributes as 15%, 18%, and 50% respectively of the total annual damage per µg/m 3 , assumed here as 2000 €/µg/m 3 .
The cost related to the enhancement of environmental safety by taking measures is the target reliability level that depends on different factors and may vary from one industry to another, as a function of the seriousness of its consequence, or public and media sensitivity.
The target level defined here is related to failure cause and mode. Considering a reference period, t r of one year, the Beta target reliability level, β t (t r ), is defined as 2.71. A more detailed recommendation is provided in [43] where the Beta target reliability index is given for the service life and related not only to the consequences but also to the relative costs of safety measures [44].
The cost of the implemented safety measures in keeping the environment safe is defined as: Risk m (x) = G air (x)+G port (x|x ) (22) where G air (x) is the cost associated with the air pollutants at the source that needs to be cleaned (abated), and G port (x|x ) is the cost associated with the port services as a function of the distance between the port and queuing seaport where the ships are waiting to enter the port. The weight of pollution discharge that needs to be cleaned up (abated) is defined as C MAC,k and the cost, associated with that is: where G air,k (x)·is the cost associated with the air pollutants at the source that need to be cleaned [33].
The anchoring areas where the ships are queuing and generating pollution are established by examining environmental, navigational, and economic conditions converting each designated anchoring area to a polygon, conditional to local judgment for a specific port. Depending on the distance to the port, size and services need to the ship a specific cost is generated, which is defined here as G port (x|x ) and may vary from port to port.
As can be seen from industry to another, as a function of the seriousness of its consequence, or public and media sensitivity. The target level defined here is related to failure cause and mode. Considering a reference period, tr of one year, the Beta target reliability level, , is defined as 2.71. A more detailed recommendation is provided in [43] where the Beta target reliability index is given for the service life and related not only to the consequences but also to the relative costs of safety measures [44].
The cost of the implemented safety measures in keeping the environment safe is defined as: where is the cost associated with the air pollutants at the source that needs to be cleaned (abated), and | is the cost associated with the port services as a function of the distance between the port and queuing seaport where the ships are waiting to enter the port. The weight of pollution discharge that needs to be cleaned up (abated) is defined as , and the cost, associated with that is: where , •is the cost associated with the air pollutants at the source that need to be cleaned [33].
The anchoring areas where the ships are queuing and generating pollution are established by examining environmental, navigational, and economic conditions converting each designated anchoring area to a polygon, conditional to local judgment for a specific port. Depending on the distance to the port, size and services need to the ship a specific cost is generated, which is defined here as | and may vary from port to port.
As can be seen from Figures 7-9 the optimal distance where the ship as a source of pollution needs to queue to keep the safe environment is 1794 m, with a concentration of the pollution at the receptor (port) with a location (0, 0, 1794 m) 9.35 μg/m 3 respecting a Beta index of 2.71 Figure 10 shows the total cost as a function of distance including the impact of , and .      Partial safety factors are developed based on the characteristic values of the governing factors Q, u, σ y , σ z , C a and H calculated at a confidence level of the original probability density function. The probabilistic design values of all parameters involved in the limit state functions are Q * , u * , σ y * , σ z * , C a * , H * , B * Q , B * u , B * σ y , B * σ z , B * C a , B * H , are respecting the Beta target reliability index, and the partial safety factors are defined as: The resulting partial safety factors can be used in the decision-making for defining the queuing (anchoring) location as the best distance from the receptor (port) conditional to the target Beta index and minimum total risk (cost) by satisfying the following design criterion: where C uγ = γ Q Q c πγ σ y σ y c γ σ z σ z c γ u u c exp − 1 2 γ H H c γ σ z σ z c 2 (28) The estimated partial safety factors are: γ Q = 1.07, γ c a = 1.49, γ σ y = 0.67, γ σ z = 0.84, γ H = 0.95, γ u = 0.85.
As can be seen from Figure 11 the ideal distance solution, satisfying the minimum risk and target reliability index is associated with 1794 m. Any distance greater than that will be an acceptable but somehow more expensive solution and any distance lesser than that is not acceptable. The solution of the problem can be defined as: min x=(x 1 , ...,x n ) C aγ − C uγ (x), x ∈ R n (29) which is constrained to the Beta target reliability level, β t (t r ), defined here as 2.71.

Conclusions
This study developed a stochastic model for predicting the air pollution concentration of ships queuing in a seaport. The Gaussian dispersion model was used to define the pollution concentration at receptor locations (crowded port areas and other reference points). Uncertainties related to the input data and mathematical model were accounted

Conclusions
This study developed a stochastic model for predicting the air pollution concentration of ships queuing in a seaport. The Gaussian dispersion model was used to define the pollution concentration at receptor locations (crowded port areas and other reference points). Uncertainties related to the input data and mathematical model were accounted for to estimate the propagating uncertainties of the air pollution emission concentration at the receptors accounting for the wind speed, horizontal and vertical dispersion parameters as a function of the geographical location of the emission sources (ships), effective emission height and weather conditions. The most expected probability of exceeding the acceptable limits of pollution concentration was defined by employing the first-order reliability method. A target probability level was defined based on the failure cause and mode and was used to develop safety calibration factors. The safety partial factors together with the Gaussian dispersion model were developed to be used for defining the safety location of the queuing ships in a seaport avoiding air pollution at the receptor (port). The developed approach is flexible and demonstrated capacity to be used in defining the distance between the port and queuing seaport in preserving the environment from the air gas pollution generated by waiting ships. Different assumptions were made that are not essential for the approach but are needed for the example calculation.
Author Contributions: Y.G. and P.G. contributed to the design and implementation of the study, to the analysis of the results and to the writing of the manuscript. All authors have read and agreed to the published version of the manuscript.

Data Availability Statement:
The data presented in this study are available within the article.