Next Article in Journal
Large-Scale Atmospheric Circulation Variability and Its Climate Impacts
Previous Article in Journal
Atmospheric Forcing of Coastal Upwelling in the Southern Baltic Sea Basin
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Feasibility of the Inverse-Dispersion Model for Quantifying Drydock Emissions

Department of Civil & Environmental Engineering, University of New Orleans, New Orleans, LA 70148, USA
*
Author to whom correspondence should be addressed.
Atmosphere 2019, 10(6), 328; https://doi.org/10.3390/atmos10060328
Submission received: 13 May 2019 / Revised: 6 June 2019 / Accepted: 12 June 2019 / Published: 17 June 2019
(This article belongs to the Section Air Quality)

Abstract

:
Important processes within the shipbuilding and ship repair industry include metal cutting, welding, surface preparation, and painting. When dealing with ship repair, ships are brought into drydocks to carry out necessary repairs. Typical repairs include but are not limited to dry or wet abrasive blasting for removing the old paint and rust followed by repainting of the external hull. Also, the painting of superstructure is carried out as necessary. Additionally, many metal cutting and welding operations are carried out. Air pollutant emissions generated from repair operations carried out within drydock are challenging to quantify, particularly if some of these repair activities do not have reliable emission factors. This paper investigates the feasibility of the inverse dispersion model for quantifying drydock emissions in a shipyard environment. The authors use a well-established Gaussian dispersion model that is used as a regulatory model in the United States and many other countries in a two-step process using a code developed in MATLAB: (1) Source-to-Receptor modeling to compute ambient concentrations using assumed emissions from various sources and meteorological conditions, and (2) The utilization of the computed ambient concentrations at various receptors to compute emissions at those sources (assumed in the first step) using the inverse Gaussian code developed.

1. Introduction

According to the 2015 United States Maritime Administration (MARAD), there are 124 shipyards in the United States that actively engage in shipbuilding and more than 200 shipyards that only engage in ship repairs are capable of building ships (but not actively engaged in shipbuilding). In 2013, the U.S. private shipbuilding and repairing industry created 399,420 jobs, $25.1 billion of labor income, and $37.3 billion in the gross domestic product (GDP) [1]. The shipbuilding industry can be divided into two categories: shipbuilding facilities and ship repair facilities [2]. The shipbuilding facilities are where shipbuilders concentrate on constructing new vessels. The ship repair facilities concentrate on repair and maintenance activities and can be further subdivided into topside repair yards and full-service repair yards. Topside repair yards are restricted to performing repairs above the waterline of the vessel, whereas dry-docking facilities or full-service repair yards, allow repairs and maintenance even below a ship’s water line.

1.1. Processes at a Shipbuilding and Repair Industry

The shipbuilding and repair industry consists of several processes that include but are not limited to surface preparation, metal cutting and surface finishing/treatment, painting and coating, machine shop operations, welding, and fiberglass operations [3,4]. Processes in the shipbuilding and repair industry are complex and of a large scale, and these activities produce significant amounts of multimedia emissions (solid, liquid, and air) [5]. Various production processes in the shipbuilding and repair industry consume a large amount of raw materials and generate significant amounts of waste materials and emissions. The raw materials used in the shipbuilding industry are steel and other metals, abrasives, paints, solvents, solvent degreasers, acid and alkaline cleaners, etc. [6]. All the processes in the industry generate waste such as leftover raw material and containers with residual raw materials. Other wastes generated are organic paint sludges, heavy metals, stripping chemicals, volatile organic compound (VOC) emissions, paints, thinners, degreasers, solvents, resins, and others. Among the air pollutants, metal-ion bearing mists (including hexavalent chromium, cadmium, arsenic, lead, nickel, and zinc), acid mists, acetone, xylene, and toluene are the most prominent pollutants at the shipbuilding and repair industry [7]. Many researchers have focused on air pollutant emissions from various processes at shipyards. Kura et al. (1999) investigated air emissions (heavy metals) from the oxyfuel cutting of steel plate ship structures [8]. Kura et al. (2000) studied hazardous air pollutant emissions from shipyard cutting operations [9]. Celebi et al. (2017) investigated air pollutant emissions from the manual metal arc welding of shipbuilding steel [10]. Mert and Ekinci (2017) studied the fume formation rate analysis of shipbuilding steel with shielded metal arc welding [11]. Celebi and Vardar (2008) investigated VOC emissions from indoor and outdoor painting processes in shipyards [6]. Kura et al. (2006) investigated atmospheric particulate emissions from dry abrasive blasting [12].
Out of all the air pollutant emissions from the industry, 98 percent are released into the air from fugitive or point sources. VOCs accounted for about 86 percent of the shipbuilding and repair industry’s reported toxic release inventory (TRI) releases [13]. In 2005, the criteria air pollutants (CAPs) from the shipbuilding and repair industry in the United States amounted to 5900 tons, and all TRI hazardous air pollutant (HAP) emissions amounted to 0.8 million pounds (0.36 million kgs) [14]. The workers are exposed to various VOCs and airborne metal fumes (e.g., nickel, zinc, and more) during processes such as painting, welding, and blasting and are prone to acute and chronic health problems [2].

1.2. Regulations for Air Pollutants in the Shipbuilding and Repair Industry

An emissions factor is a representative value that attempts to relate the quantity of a pollutant released to the atmosphere with an activity associated with the release of that pollutant. These factors are usually expressed as the weight of the pollutant divided by a unit weight, volume, distance, or duration of the activity emitting the pollutant (e.g., kilograms of particulate emitted per megagram of coal burned) [15].
Under Title III of the 1990 Clean Air Act Amendments (CAAA), the United States Environmental Protection Agency (USEPA) is required to develop national emission standards for 189 hazardous air pollutants (NESHAP). The USEPA is developing maximum achievable control technology (MACT) standards for all new and existing sources. The National Emission Standards for Shipbuilding and Repair Operations (Surface Coating) were finalized in 1995 and apply to major source shipbuilding and ship repairing facilities that carry out surface coating operations [13].

1.3. Importance of Research

Although USEPA does define limits to the amount of air pollutant emissions from all the point and fugitive sources at the ship repair industries, to accurately quantify emissions from a drydock is a big challenge for all. The reason drydock emissions are difficult to quantify because a drydock is large in the area which is used to carry out various repair/maintenance activities such as welding, cutting, painting and solvent degreasing. Drydocks allow shipyards to carry out ship repair and construction activities. These activities release criteria and hazardous air pollutants. However, these air pollutants are not easy to monitor and quantify due to the complexity of repair and construction operations.
The purpose of this paper is to check the feasibility of the inverse dispersion model [16,17,18,19,20,21,22,23] for quantifying the emissions from Drydocks. The following sections discuss the model in more detail. This model will not only help the industry identify the actual emissions from drydocks but also help them comply with regulations. Another concept of Equivalent Source is introduced in this article. According to this experimental study, any area source can be represented by a single “equivalent source” instead of having several sources to represent the area source such as a drydock in this case.

2. Methodology

2.1. Gaussian Dispersion Model

The Gaussian dispersion model has been the basis of many air quality models used for air quality permitting, compliance evaluation, and health impact assessments. The USEPA regulatory models such as ISCST3 [24], AERMOD [25] principally use the basic concept of the Gaussian dispersion equation with many advanced features/modifications to compute the ambient concentrations resulting from an emitting source. The Gaussian dispersion model used to compute the ambient concentration of a pollutant at a receptor downwind from a point source is illustrated in Figure 1.
The ambient concentration (µg/m3) of a pollutant at a receptor downwind from a point source is calculated using the following equation [26]. The equation is well known as the “Gaussian Dispersion Equation.”
C = Q 2 π u σ y σ z exp 1 2 y 2 σ y 2 exp ( z H e ) 2 2 σ z 2 + exp ( z + H e ) 2 2 σ z 2 ,
where
C = Concentration at a receptor (x, y, z) distance from a point source (μg/m3),
Q = Emission rate (μg/s),
σy and σz = Horizontal and vertical dispersion parameters (m) that are functions of the downwind distance (x) and atmospheric stability,
u = Average wind speed at the physical stack height (m/s),
y = Crosswind distance from ground level (m),
z = Vertical distance above the ground (m), and
He = Effective stack height (physical stack height plus plume rise) (m).
For most modeling scenarios, the emission rate is determined either through source monitoring or using reliable emission factors. Emission rate, meteorological parameters, and other input parameters are used to compute the ambient concentrations at various receptors from an emitting source using the Gaussian dispersion equation. The determination of the emission rate is challenging for a drydock as it does not have a defined source location. The drydocks are more like an area source, where source monitoring will not provide accurate results since the sources are not defined. This paper attempts to develop a model to compute the emission rate(s) from drydocks using the monitored ambient concentrations using the inverse dispersion model.
Usually, in the case of shipyards, air pollutant emission rates (Qs) from drydocks are unknown due to a lack of reliable emission factors for many repair activities carried out in drydocks. Shipyards, however, do have several shops such as blasting houses, welding shops, metal cutting shops, and paint booths which have stacks attached to which the air pollutant emissions can be quantified. The quantification of air pollutant emissions (Qs) for these shops/operations is either carried out by (a) source monitoring (stack testing) and/or (b) application of reliable emission factors.

2.2. Inverse Gaussian Dispersion Model

Figueroa et al. (2008) successfully demonstrated the Inverse Gaussian Dispersion Model that uses the concentrations calculated or measured in the downwind direction of a source to calculate the emission rates from unknown sources in a landfill case study [17]. The Gaussian equation discussed previously can be described as a single function linking C and Q. This relationship between C and Q is a function of x (source-to-receptor downwind distance), and y (source-to-receptor crosswind distance) as all the variables in the function are dependent on x and y.
F = 1 2 π u σ y σ z exp ( 1 2 y 2 σ y 2 ) exp ( z H e ) 2 2 σ z 2 + exp ( z + H e ) 2 2 σ z 2 ,
Assuming ground level receptors (z = 0), Equation (2) reduces to Equation (3) as shown below
F = 1 π u σ y σ z exp ( y 2 2 σ y 2 ) exp H e 2 2 σ z 2 ,
Since F is a function of x and y, it can be written as f(x,y)r,s. Hence, the concentration can be re-written as
Cr,s = f(x,y)r,s × Qs,
where
Cr,s = Concentration on the receptor “r” due to source “s”,
f(x,y)r,s = Function f(x,y) calculated for source “s” affecting receptor “r”, and
Qs = Emissions from source “s”.
In order to demonstrate the inverse dispersion model, let us assume C is a 3 × 1, F is a 3 × 2, and Q is a 2 × 1 matrix as follows
      C =   C r 1 C r 2 C r 3       F =   a 11 a 12 a 21 a 22 a 31 a 32       Q =   Q s 1 Q s 2
where
Cr1 = Concentration measured/modeled at receptor 1
ars = Effect of source “s” on the receptor “r,” calculated by the function “F.”
Qs1 = Emissions due to source 1
As discussed earlier, generally the concentrations are calculated using the known emission rates Q and the function F. In the inverse dispersion model, in order to compute the emission rates Q with known concentrations C, the following relationship can be used by inverting the function f(x,y):
Q s = F * C measured ,
where F* is the inverse of the function f(x,y). Since F is not a square matrix, we perform the following to invert it.
F * =   F T F 1 F T ,
where
F T =   a 11 a 12 a 13 a 21 a 22 a 23
F T F 1 =   b 11 b 12 b 21 b 22
F T F 1 F T =   b 11 b 12 b 21 b 22 a 11 a 12 a 13 a 21 a 22 a 23
F * =   F 11 F 12 F 13 F 21 F 22 F 23
Finally, the emission rate Q is computed using F* and the measured concentrations C as shown in Equation (7).
Q s =   Q s 1 =   C r 1 F 11 +   C r 2 F 12 +   C r 3 F 13 Q s 2 =   C r 1 F 21 +   C r 2 F 22 +   C r 3 F 23 ,

2.3. Best Shipyard Case

The Best Shipyard is a hypothetical concept created by the Maritime Environmental Resources and Information Center (MERIC) at the University of New Orleans (UNO) for educational purposes to demonstrate various shipyard processes and multimedia emissions to the students and environmental engineers/professionals, and to develop software applications. This was helpful in communicating among environmental engineering and software developers. For demonstrating the inverse dispersion model in this paper, the same hypothetical shipyard, “the Best Shipyard” was used with a layout as shown in Figure 2. This Best Shipyard has three drydocks (unknown emission rates, Qs) and seven other stack sources (known emission rates, Qs).

2.4. Inverse Dispersion Model for Drydocks

To evaluate the feasibility of the inverse dispersion model for drydocks, the authors in this research consider two approaches as follows (a) assuming one equivalent point source in each drydock and (b) assuming nine point sources in each drydock, considering it as an area source. The model involved in demonstrating the inverse dispersion model to drydock scenario can be described in six (6) distinct steps mentioned in the following section. For this exercise, the emission rates from all the sources in the shipyard (sources include the drydocks and other stack sources) are assumed as depicted in Table 1. For the purpose of this paper and the intended demonstration, the “actual” emissions from the ten point sources were assumed as shown in Table 1. Sources S1–S3 are the drydocks and S4–S10 are the stack sources.
A grid was laid over the Best Shipyard layout as depicted in Figure 3 with the triangular markers showing the source locations emitting air pollutants (three drydock sources and seven stack sources) and the circular markers represent the receptor locations where the concentrations are monitored. The receptor points are the points where the air pollutant concentrations can be either measured using an air quality monitor or modeled using air quality models such as atmospheric dispersion models. In this paper, modeled values at the receptors are being used but in the actual case, one must use measured/monitored air pollutant concentrations which can be achieved by either grab sampling method or real-time monitoring.
To use the Gaussian dispersion equation, the required meteorological conditions used/ assumed are as follows:
Wind Speed: 5 m/sec
Wind Direction (degree): 270 (wind blowing from West to East)
Pasquill (Atmospheric) Stability Class: C.
The six steps in discussing the model are as follows:

2.4.1. Step 1: Ten (10) Point Sources with Known Emission Rates (Qs) to Calculate Concentrations (Cs) at Seventy (70) Receptors

The ambient concentrations resulting from known emissions from the ten sources were calculated at all 70 receptor locations. The concentrations were calculated using the Gaussian dispersion model as depicted by Equation (1). The coordinates (x and y) were computed using the source and receptor grids and are shown in Table 1 (source co-ordinates) and Table 2 (receptor co-ordinates). The drydock emission height was considered to be “0 meters” as it is believed that the most emissions from drydocks are emitted ground level. The emission height for all seven other point sources (stack sources) was considered to be “10 meters” as all the sources are connected to a stack to spew the emissions. Since there are 70 receptors and ten sources, there will be 700 (70 × 10) calculations. The MATLAB program was used to compute the pollutant concentration in ambient air at seventy (70) receptors as shown in Table 2.
The source locations (three drydocks and seven stack sources) along with receptor (70) locations as depicted using MATLAB is shown in Figure 4.

2.4.2. Step 2: Step 2: Known Cs at 70 Receptors to Calculate Three Unknown Drydock Qs

In a real-world scenario, the concentrations are known using monitoring instruments of a different kind, depending on the type of pollutant at hand. Once these measured/monitored concentrations at various locations (receptors) are known, one has to identify the possible sources and then use the inverse dispersion model [17] to back-calculate the emission rates for these identified sources. In this paper, the calculated ambient concentrations in Step 1 will be treated as measured concentrations from this point. In Step 2, the demonstration of reverse calculation is implemented by treating the three drydock sources as unknown sources, while emissions from the stack sources are known. The concentrations from Table 2 were used to implement the inverse dispersion model to compute the emission rates from the three drydocks. Again, MATLAB was used to compute the unknown emission rates, which are shown in Table 3. These “computed” emission rates were compared with “actual” emission rates from Table 1 as can be seen in Table 3.

2.4.3. Step 3: Using only the Three Drydocks as the Sources with Known Emission Rates (Qs) to Calculate Cr, DS at 70 Receptors

In Step 3, three drydocks with emission rates (Qs) calculated in Step 2 were used to compute the ambient concentrations at the seventy receptors using the direct Gaussian dispersion model. The purpose was to see whether the concentrations at these 70 receptors will change with a change in the number and location of sources within the facility. These computed concentrations are depicted in Table 4. The term Cr, DS is used to explain that the concentration is calculated for the receptor ‘r’ which is caused due to the drydock sources ‘DS.’ The results of this exercise are shown in Table 4.

2.4.4. Step 4: Seventy Receptors with Known Cr, DS to Compute Qs from Nine Sources in Each of Three Drydocks (Total of 27 Sources in 3 Drydocks)

The second aspect of the research is considering each drydock as an area source. As mentioned earlier, there is no one particular emission source for a drydock. Often, repair work is carried out at many locations in the drydock, and the actual emissions take place from multiple locations within a drydock. In order to simulate this situation, a total of twenty-seven (27) point sources were considered with nine (9) sources in each of the three (3) drydocks. Figure 5 shows the distribution of sources in each drydock. The pollutant concentrations in ambient air available from Step 3 were used in combination with the inverse dispersion model to compute the unknown Qs at twenty-seven (27) sources. Emission rates obtained from this exercise are depicted in Table 5. The reason for this exercise is to see the difference in total emissions from each drydock between the two alternatives considered, (a) one source in each drydock and (b) nine sources in each drydock.

2.4.5. Step 5: Known Qs at 27 Sources (3 Drydocks) to Compute Cs at 30 New Receptors

In Step 5, all 27 sources within three drydocks were considered, and ambient concentrations were computed using the Gaussian equation (MATLAB code) at 30 receptors as shown in Figure 6 (new receptor grid; different from the previous 70 receptor grid). All other parameters such as wind speed and stability class were kept the same. The purpose is to see how the concentrations at these 30 receptors will change from (a) total of 3 sources, one source in each of three drydocks to (b) total of 27 sources, nine sources in each of the three drydocks. Source-receptor geometry is illustrated in Figure 6. The computed concentrations are shown in Table 6. This demonstration will help study the health impacts of the facility at various receptor locations that are of interest.

2.4.6. Step 6: Known Qs at 3 Sources (3 Drydocks) to Compute Cs at 30 New Receptors

Now knowing the concentrations at the 30 receptors, the concept of the equivalent source is implemented using the same set of 30 receptors by combining the nine sources at each drydock to consider the drydock as a single source. Step 6 is the same as Step 5 except that only three sources (one source in each drydock) were considered for computing the concentration in ambient air at 30 newly selected receptors. The results are shown in Table 7, and the source-receptor geometry is illustrated in Figure 7. This step shows that the concept of the equivalent source can be implemented since the drydock can be considered as one source. As seen in Table 6, the error margin between the two scenarios (considering drydocks as a point source and area source) is of the magnitude 10−4.

3. Results and Discussion

Step 1 provided the results of the ambient pollutant concentration at 70 receptors, which are shown in Table 2. These ambient concentrations are computed using a simple Gaussian dispersion method described by Equation (1).
Step 2 utilizes the inverse Gaussian model proposed by Figueroa et al. (2008). This step used the ambient concentrations computed (in this case result of Step 1; field measured ambient concentrations can be used in actual applications) to compute emission rates. The results are obtained using the MATLAB code utilizing Equation (3)–(5). Here, 7 out of the ten sources considered have known emission rates since they are stack (point) sources and can be measured easily. The three drydocks are considered as unknown sources, and the emissions are computed using this inverse dispersion model described in this paper. The results of modeled Qs for the three unknown drydock sources are shown in Table 3. From Step 2, it can be seen that the modeled results are the same as the actual emission rates (in this case, the assumed emission rates which were used in Step 1). Thus, this model proves to be very promising. Here, the actual ambient concentrations were unknown; the ambient concentrations at various receptors are calculated in Step 1. In the field, however, the concentrations will be measured using a simple hand-held instrument and carry on with Step 2 to compute the emission rate which is unknown.
Step 3 provided the results of the modeled concentrations (Cr, DS) at receptor “r” due to source “s” at 70 receptors from only three sources (one in each drydock). The results of this step are shown in Table 4. The concentration at each receptor is the sum of effect due to all the three sources at that receptor.
In Step 4, the authors provided the results of the emission rates (Qs) at 27 sources (nine sources in each drydock) that were calculated using the ambient concentrations at 70 receptors. These results are presented in Table 5.
Steps 5 and 6 resulted in modeled ambient concentrations at 30 new receptors as a result of emissions from 27 sources located in three drydocks (nine sources in each drydock) and emissions from three sources located in three drydocks (one source in each drydock) which are given in Table 6.
The concept of “equivalent source” is proved in Step 4–7. The comparison shown in Table 7 proves that treating the drydock as a single point source gives comparable results with the case when the drydock is considered to have nine source locations throughout the area of the drydock as shown in Figure 5. However, caution should be used in extrapolating this approach and should be used for situations that are tested through additional validation.
The inverse dispersion model was evaluated for quantifying drydock emissions. This model is scalable and can be applied to quantify emissions from animal farms, landfills, wastewater treatment plants, forest fires and so on.
The paper mainly provided a novel method of using the Gaussian dispersion model in a reverse way to solve complex emission inventory problems that exist in the maritime industry. The authors believe that this model has significant benefits as stated below:
  • The proposed model uses the proven scientific method of the Gaussian Dispersion Model, which is used by the USEPA as a regulatory model. Many countries around the World use the same model for regulatory compliance and approving air permits. The model’s limitations are well recognized, and efforts are constantly being made to improve the assumptions and algorithm. For now, the model is being used “as-is” for regulatory and litigation purposes.
  • The quantification of drydock emissions is complex due to the type of operations and their variations from shipyard to shipyard. Thus, the use of emission factors available from the USEPA’s AP-42 is insufficient. There is no better method to estimate emissions from drydocks unless this area is enclosed, and source testing is conducted for each situation.
  • The inverse dispersion model proposed is relatively easier to use/implement as the Gaussian Dispersion Model is already being used for regulatory purposes and is also less resource demanding.
  • The lack of a rational emission quantification method will continue to result in erroneous emission estimates, thus erroneous risk assessments.
  • The proposed model can be easily applied to other complex emission sources such as landfills, superfund sites, raw material yards, large mine sites, forest fires, wastewater treatment plants and many more.

4. Conclusions

This paper demonstrates the feasibility of the inverse dispersion model for quantifying emissions from drydocks that are not easily quantifiable. The results prove that this model is reliable in estimating the emission rates from sources that are not quantifiable through source monitoring methods and those that do not have reliable emission factors. The same model can be applied to quantify emission rates from landfills, animal farms, and wastewater treatment plants where the reliable emission factors are absent. Emission quantification from these sources can help various sources comply with regulations, perform impact analysis, and conduct health risk assessments in the nearby community. In most countries, it is a big challenge for any governmental agency to quantify the actual emissions from forest fires. By implementing the inverse dispersion model, it will be easier to quantify the emissions from forest fires, and this can help agencies in reporting the emissions to develop emission control strategies and health risk assessments. However, the authors caution the use of the proposed inverse dispersion model as this method relies on the applicability of the Gaussian dispersion model. Whenever and wherever the Gaussian dispersion model cannot be utilized, the proposed inverse model will have the same limitations.

Author Contributions

B.K. conceived the idea, designed the study outline, and guided A.J.; A.J. performed the literature review, developed the MATLAB code, and drafted the paper.

Funding

This research received no external funding.

Acknowledgments

The authors would like to acknowledge the Maritime Environmental Resources and Information Center (MERIC), University of New Orleans for the research support.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Maritime Administration (MARAD). The Economic Importance of the U.S. Shipbuilding and Repairing Industry. Available online: https://www.maritime.dot.gov/sites/marad.dot.gov/files/docs/resources/3641/maradeconstudyfinalreport2015.pdf (accessed on 9 May 2019).
  2. Dincer, I.; Midilli, A.; Hepbasli, A.; Karakoc, T.H. (Eds.) Global Warming Engineering Solutions Series: Green Energy and Technology; Springer: Berlin, Germany, 2010; pp. 579–590. [Google Scholar]
  3. Celebi, U.B.; Ekinci, S.; Alarcin, F.; Unsalan, D. The Risk of occupational safety and health in shipbuilding industry in Turkey. In Proceedings of the 3rd International Conference on Maritime and Naval Science and Engineering, Constantza, Romania, 3–5 September 2010. [Google Scholar]
  4. Papaioannou, D. Environmental implications related to the shipbuilding and ship repairing activity in Greece. J. Marit. Transp. Sci. 2003, 41, 241–252. [Google Scholar]
  5. Kura, B.; Lacoste, S.; Patibanda, P. Multimedia pollutant emissions from shipbuilding facilities. In Proceedings of the United States and Japan Natural Resources (UJNR) Conference, Washington, DC, USA, 25 October–4 November 1998. [Google Scholar]
  6. Celebi, U.B.; Vardar, N. Investigation of VOC Emissions from Indoor and Outdoor Painting Processes in Shipyards. Atmos. Environ. 2008, 42, 5685–5695. [Google Scholar] [CrossRef]
  7. Kura, B. Multi-media pollutants in the shipbuilding industry and the importance of environmental management. In Proceedings of the 1st international symposium on Naval Architecture and Maritime, Yildiz Technical University, Istanbul, Turkey, 24–25 October 2011. [Google Scholar]
  8. Kura, B.; Wisbith, A.S.; Stone, R.; Judy, T. Metal Cutting Operations: Emission Factors for Particulates, Metals and Metal Ions. In Proceedings of the Emission Inventory: Regional Strategies for the Future Proceedings, Raleigh, NC, USA, 26–28 October 1999; pp. 1–13. [Google Scholar]
  9. Kura, B.; Judy, T.; Wisbith, S.A.; Stone, R. Assessment of air emissions from shipyard cutting processes. In Proceedings of the United States Japan Natural Resources (UJNR) Conference, Tokyo, Japan, 17–18 May 2000; pp. 1–18. [Google Scholar]
  10. Celebi, U.B.; Mert, T.; Bilgili, L.; Senoz, K.M.; Ekinci, S. Characterization of sub-micrometer fume particles in MMA welding of shipbuilding steel with different types of electrodes. Fresenius Environ. Bull. 2017, 26, 140–145. [Google Scholar]
  11. Mert, T.; Ekinci, S. Fume formation rate analysis of shipbuilding steel with SMAW using Taguchi design and ANOVA. Acta Phys. Pol. A. 2017, 13, 495–499. [Google Scholar] [CrossRef]
  12. Kura, B.; Kambham, K.; Sangameswaran, S.; Potana, S. Atmospheric particulate emissions from dry abrasive blasting using coal slag. J. Air Waste Manag. Assoc. 2006, 56, 1205–1215. [Google Scholar] [CrossRef] [PubMed]
  13. U.S. Environmental Protection Agency. Office of Compliance Sector Notebook Project: Profile of the Shipbuilding and Repair Industry; EPA/310-R-97-008; U.S. Environmental Protection Agency: Washington, DC, USA, 1997.
  14. U.S. Environmental Protection Agency. Sector Performance Report. Shipbuilding & Ship Repair. Available online: https://archive.epa.gov/sectors/web/pdf/shipbuilding.pdf (accessed on 9 May 2019).
  15. U.S. Environmental Protection Agency. Technology Transfer Network Air Toxics Web Site. Emissions Factors & AP 42, Compilation of Air Pollutant Emission Factors. Available online: http://www.epa.gov/ttnchie1/ap42/ (accessed on 9 May 2019).
  16. Flesch, T.K.; Wilson, J.D.; Harper, L.A.; Crenna, B.P. Estimating gas emissions from a farm with an inverse-dispersion technique. Atmos. Environ. 2005, 39, 4863–4874. [Google Scholar] [CrossRef]
  17. Figueroa, V.K.; Cooper, C.D.; Mackie, K.R. Estimating landfill greenhouse gas emissions from measured ambient methane concentrations and dispersion modeling. In Proceedings of the 101st Air and Waste Management Association Annual Conference and Exhibition, Portland, OR, USA, 24–27 June 2008. [Google Scholar]
  18. Faulkner, W.B.; Goodrich, L.B.; Botlaguduru, V.S.V.; Capareda, S.C.; Parnell, C.B. Particulate matter emission factors for almond harvest as a function of harvester speed. J. Air Waste Manag. Assoc. 2009, 59, 943–949. [Google Scholar] [CrossRef] [PubMed]
  19. Figueroa, V.K.; Mackie, K.R.; Guarriello, N.; Cooper, C.D. A robust method for estimating landfill methane emissions. J. Air Waste Manag. Assoc. 2009, 59, 925–935. [Google Scholar] [CrossRef] [PubMed]
  20. Flesch, T.K.; Harper, L.A.; Desjardins, R.L.; Gao, Z.; Crenna, B.P. Multi-source emission determination using an inverse-dispersion technique. Bound. Layer Meteorol. 2009, 132, 11–30. [Google Scholar] [CrossRef]
  21. Lushi, E.; Stockie, J.M. An inverse Gaussian plume approach for estimating atmospheric pollutant emissions from multiple point sources. Atmos. Environ. 2010, 44, 1097–1107. [Google Scholar] [CrossRef] [Green Version]
  22. Schauberger, G.; Piringer, M.; Knauder, W.; Petz, E. Odour emissions from a waste treatment plant using an inverse dispersion technique. Atmos. Environ. 2011, 45, 1639–1647. [Google Scholar] [CrossRef]
  23. Bonifacio, H.F.; Maghirang, R.G.; Auvermann, B.W.; Razote, E.B.; Murphy, J.P.; Harner, J.P., III. Particulate matter emission rates from beef cattle feedlots in Kansas-reverse dispersion modeling. J. Air Waste Manag. Assoc. 2012, 62, 350–361. [Google Scholar] [CrossRef] [PubMed]
  24. U.S. Environmental Protection Agency. United States Environmental Protection Agency, “User’s Guide for the Industrial Source Complex (ISC3) Dispersion Models. Volume II—Description of Model Algorithms”; U.S. Environmental Protection Agency: Washington, DC, USA, 1995.
  25. Perry, S.G.; Cimorelli, A.J.; Lee, R.F.; Paine, R.J.; Venkatram, A.; Weil, J.C.; Wilson, R.B. AERMOD: A dispersion model for industrial source applications. In Proceedings of the 87th Annual Meeting Air and Waste Management Association, Pittsburgh, PA, USA, 19–24 June 1994; p. 16. [Google Scholar]
  26. Cooper, C.D.; Alley, F.C. Air Pollution Control, A Design Approach, 4th ed.; Waveland Press: Long Grove, IL, USA, 2011. [Google Scholar]
Figure 1. Gaussian dispersion model [26].
Figure 1. Gaussian dispersion model [26].
Atmosphere 10 00328 g001
Figure 2. The best shipyard–A hypothetical layout for demonstration.
Figure 2. The best shipyard–A hypothetical layout for demonstration.
Atmosphere 10 00328 g002
Figure 3. The Best Shipyard case with the grid superimposed in MATLAB. The triangles and circles depict sources and receptors respectively.
Figure 3. The Best Shipyard case with the grid superimposed in MATLAB. The triangles and circles depict sources and receptors respectively.
Atmosphere 10 00328 g003
Figure 4. Grid showing ten Sources and 70 receptors as generated in MATLAB.
Figure 4. Grid showing ten Sources and 70 receptors as generated in MATLAB.
Atmosphere 10 00328 g004
Figure 5. (a) One source in each of the three drydocks vs. (b) nine sources in each of the three drydocks.
Figure 5. (a) One source in each of the three drydocks vs. (b) nine sources in each of the three drydocks.
Atmosphere 10 00328 g005
Figure 6. Source-receptor geometry for nine sources at each drydock and a total of 27 sources at the facility.
Figure 6. Source-receptor geometry for nine sources at each drydock and a total of 27 sources at the facility.
Atmosphere 10 00328 g006
Figure 7. New source-receptor grid with the three drydocks.
Figure 7. New source-receptor grid with the three drydocks.
Atmosphere 10 00328 g007
Table 1. Emissions assumed for each source.
Table 1. Emissions assumed for each source.
Sourcex-cord
(in m)
y-cord
(in m)
Height (H)
(in m)
Actual Emissions
(in μg/s)
S15010007000
S25015007520
S35030008540
S48040107480
S580340109400
S61601201012000
S7160280101359
S81803801016000
S92101001018800
S10320160101585
Table 2. Receptor locations and concentrations as calculated in Step 1.
Table 2. Receptor locations and concentrations as calculated in Step 1.
Receptorx-Coordy-coordConcentration (µg/m3)Receptorx-Coordy-CoordConcentration (µg/m3)
R110000R362502500.207
R2100500R372503001.798
R310010015.667R382503501.256
R410015016.830R392504000.098
R51002000R402504500
R61002500R4130000.253
R710030019.113R42300500.926
R81003500R433001003.685
R91004000R443001501.531
R101004500R453002000.240
R1115000R463002500.315
R12150500.179R473003001.457
R131501004.488R483003501.474
R141501504.820R493004000.987
R151502000.004R503004500
R161502500.005R5135000.288
R171503005.470R52350500.876
R181503500.225R533501004.431
R191504000R543501501.561
R201504500R553502000.276
R2120000.036R563502500.376
R22200500.851R573503001.211
R232001002.232R583503501.736
R242001502.386R593504001.433
R252002000.075R603504500.012
R262002500.085R6140000.293
R272003002.677R62400500.990
R282003501.067R634001003.877
R292004000.001R644001501.672
R302004500R654002000.315
R3125000.161R664002500.400
R32250501.001R674003001.026
R332501001.755R684003501.716
R342501501.592R694004001.419
R352502000.178R704004500.063
Table 3. Comparison of assumed and calculated emission rates in Step 2.
Table 3. Comparison of assumed and calculated emission rates in Step 2.
Sourcex-Coord
(in m)
y-Coord
(in m)
Height (H)
(in m)
Actual Emissions
(in μg/s)
Calculated Emissions
(in μg/s)
Difference
(in μg/s)
S150100070007000.011−0.011
S250150075207520.015−0.015
S350300085408539.9860.014
S480401074807480.286−0.286
S5803401094009400.021−0.021
S61601201012,00011,999.8910.109
S71602801013591358.9950.005
S81803801016,00016,000.189−0.189
S92101001018,80018,800.081−0.081
S103201601015851584.6050.396
Table 4. Concentrations at all 70 receptors due to the three drydocks alone, as calculated in Step 3.
Table 4. Concentrations at all 70 receptors due to the three drydocks alone, as calculated in Step 3.
ReceptorMeasured Concentration
(in µg/m3)
ReceptorMeasured Concentration
(in µg/m3)
ReceptorMeasured Concentration
(in µg/m3)
ReceptorMeasured Concentration
(in µg/m3)
R10R190R371.565R550.271
R20R200R380.201R560.291
R315.667R210R390R570.753
R416.830R220.070R400R580.278
R50R232.231R410.004R590.014
R60R242.386R420.220R600
R719.113R250.075R431.090R610.023
R80R260.085R441.138R620.244
R90R272.631R450.240R630.703
R100R280.085R460.267R640.722
R110R290R471.046R650.287
R120.004R300R480.264R660.293
R134.488R310R490.004R670.571
R144.820R320.165R500R680.268
R150.004R331.460R510.012R690.028
R160.005R341.543R520.241R700.001
R175.470R350.178R530.862
R180.005R360.201R540.891
Table 5. Emission rates calculated in Step 4 for drydocks as an area source with nine source locations at each drydock.
Table 5. Emission rates calculated in Step 4 for drydocks as an area source with nine source locations at each drydock.
Sourcex-Cord (in m)y-Cord (in m)Calculated Emission (µg/s)Sourcex-Cord (in m)y-Cord (in m)Calculated Emission (µg/s)Sourcex-Cord (in m)y-Cord (in m)Calculated Emission (µg/s)
S11301100S21301600S31303100
S12301000S22301500.158S32303000
S1330900S23301400S33302900
S14501100S24501600S34503100.098
S15501006997.936S25501507519.638S35503008538.738
S1650900.204S26501400.294S36502900.081
S17701100S27701600S37703100.128
S18701000.720S28701500.069S38703000.429
S1970900.860S29701400S39702900.353
S1TOTAL6999.729S2TOTAL7520.159S3TOTAL8539.827
Table 6. Comparison of calculated ambient concentrations.
Table 6. Comparison of calculated ambient concentrations.
ReceptorConcentration Using Single Equivalent Source (in µg/m3)Concentration Using Nine Sources in Each Drydock (in µg/m3)DifferenceReceptorConcentration Using Single Equivalent Source (in µg/m3)Concentration Using Nine Sources in Each Drydock (in µg/m3)Difference
R10.9150.9150R165.4705.4700
R22.3862.3860R17000
R30.0030.0030R180.0030.0030
R42.6312.6310R190.0440.0440
R50.7240.7240R200.1180.1180
R61.3521.3520R21000
R70.0550.0550R220.0360.0360
R81.3181.3180R23000
R90.5200.5200R240.1000.1000
R100.8170.8170R25000
R110.1770.1770R260.4200.4200
R120.6700.6700R270.5990.5990
R130.7610.7610R280.2330.2330
R144.8204.8200R290.4500.4500
R15000R300.0020.0020
Table 7. Comparison of total emissions at each drydock (single equivalent source vs. nine sources in each drydock).
Table 7. Comparison of total emissions at each drydock (single equivalent source vs. nine sources in each drydock).
Drydock NumberEmissions (µg/s)
(Single Equivalent Source in Each Drydock)
Emissions (µg/s)
(Nine Sources in Each Drydock)
Difference
17000.0116999.7190.292
27520.0157520.159−0.145
38539.9868539.8270.159

Share and Cite

MDPI and ACS Style

Kura, B.; Jilla, A. Feasibility of the Inverse-Dispersion Model for Quantifying Drydock Emissions. Atmosphere 2019, 10, 328. https://doi.org/10.3390/atmos10060328

AMA Style

Kura B, Jilla A. Feasibility of the Inverse-Dispersion Model for Quantifying Drydock Emissions. Atmosphere. 2019; 10(6):328. https://doi.org/10.3390/atmos10060328

Chicago/Turabian Style

Kura, Bhaskar, and Abhinay Jilla. 2019. "Feasibility of the Inverse-Dispersion Model for Quantifying Drydock Emissions" Atmosphere 10, no. 6: 328. https://doi.org/10.3390/atmos10060328

APA Style

Kura, B., & Jilla, A. (2019). Feasibility of the Inverse-Dispersion Model for Quantifying Drydock Emissions. Atmosphere, 10(6), 328. https://doi.org/10.3390/atmos10060328

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop