User:Milton Beychok/Sandbox: Difference between revisions

From Citizendium
Jump to navigation Jump to search
imported>Milton Beychok
No edit summary
No edit summary
 
(489 intermediate revisions by one other user not shown)
Line 1: Line 1:
In the past, the Owens Lake dry bed area was characterized as being the largest single source of  PM<sub>10</sub> (particulate matter less than 10 microns in size) in the [[United States]] by the [[Great Basin Unified Air Pollution Control District]] (GBUAPCD)<ref name=PM-HealthEffects>[http://www.gbuapcd.org/Information/OwensLakeParticulateMatterHealthEffects.htm Survey of Reported Health Effects of Owens Lake Particulate Matter] January, 2000. Retrieved from the website of the Great Basin Unified Air Pollution Control District (GBUAPCD) on 2011-01-31.</ref><ref>'''Note:''' The primary air pollution regulatory agency in the state of California is the [[California Air Resources Board]] (CARB). The GBUAPCD is one of the 35 air pollution control districts, which report to CARB and which address air pollution problems at the local level.</ref>. That characterization was based on an approximate, unscientific estimate made in about 1986 that the annual PM<sub>10</sub> emissions were 300,000 tons.<ref>Personal correspondence in an email from Mr. Ted Schade of the GBUAPCD on February 3, 2011.</ref> Subsequently, the GBUAPCD monitored and scientifically tested the PM<sub>10</sub> emissions over a long period, including daily measurements from the period of July 2000 to 2001, and determined the annual PM<sub>10</sub> emissions were only 79,000 tons for that year-long test period.<ref name=EPA-Conference>[http://www.epa.gov/ttn/chief/conference/ei12/fugdust/present/ono.pdf Quantifying Particulate Matter Emissions from Wind Blown Dust Using Real-time Sand Flux Measurements] Duane Ono and Scott Weaver of the Great Basin Unified Air Pollution Control District, and Ken Richmond of MFG, Inc., U.S. EPA Emission Inventory Conference, San Diego, California, April 2003. From website of the U.S. EPA. Retrieved on 2011-01-31.</ref> That is equivalent to an average of 217 tons per day. During that test period, a single  peak daily value of 7,200 tons  was encountered.  
{{AccountNotLive}}
__NOTOC__
[[File:Crude oil-fired power plant.jpg|thumb|right|225px|Industrial air pollution source]]
Atmospheric dispersion modeling is the mathematical simulation of how air pollutants disperse in the ambient atmosphere. It is performed with computer programs that solve the mathematical equations and algorithms which simulate the pollutant dispersion. The dispersion models are used to estimate or to predict the downwind concentration of air pollutants emitted from sources such as industrial plants, vehicular traffic or accidental chemical releases.  


During the 1980s and the early 1990s, the [[U.S. Environmental Protection Agency]] (U.S. EPA) had designated the Owens Lake area as having not attained the [[National Ambient Air Quality Standards]] (NAAQS) for PM<sub>10</sub>. The wind-blown sand and  PM<sub>10</sub> air pollution affected approximately 40,000 residents living near Owens Lake. Regionally the blowing sand and dust also affected [[U.S. Forest Service]] employees and the [[U.S. Navy]]'s air weapons station located at nearby [[ China Lake]].<ref name=EPA-Region9>[http://epa.gov/region9/air/owens/history.html Owens Valley, CA Particulate Matter Plan: History] From website of the U.S. EPA. Retrieved on 2011-01-31.</ref> During those same years, nearby residents in [[Keeler, California]] reportedly were exposed to 25 unhealthful days a year while residents in the town of [[Ridgecrest, California]] experienced 10 unhealthful days per year due to the PM<sub>10</sub> air pollution. An emergency physician at the Ridgecrest Community Hospital was reported as saying "When we see the white cloud headed down through the pass, the ER and doctors' offices fill up with people who suddenly got worse. It's a pretty straightforward cause and effect.”&thinsp;<ref name=USGS>[http://geochange.er.usgs.gov/sw/impacts/geology/owens/ Owens (Dry) Lake, California: A Human-Induced Dust Problem] Abstracted from: Reheis, M.C., in press, ''Dust deposition downwind of Owens (Dry) Lake, 1991-1994: Preliminary findings''. Retrieved from website of the [[U.S. Geological Survey]] on 2011-01-31.</ref>
Such models are important to governmental agencies tasked with protecting and managing the ambient air quality. The models are typically employed to determine whether existing or proposed new industrial facilities are or will be in compliance with the National Ambient Air Quality Standards (NAAQS) in the United States or similar regulations in other nations. The models also serve to assist in the design of effective control strategies to reduce emissions of harmful air pollutants. During the late 1960's, the Air Pollution Control Office of the U.S. Environmental Protection Agency (U.S. EPA) initiated research projects to develop models for use by urban and transportation planners.<ref>J.C. Fensterstock et al, "Reduction of air pollution potential through environmental planning", ''JAPCA'', Vol. 21, No. 7, 1971.</ref>
In the late 1990s, after many years of litigation and negotiations that included the city of Los Angeles through its Department of Water and Power (LADWP), the Great Basin Unified Air Pollution Control District (GBUAPCD), the U.S. EPA and other interested parties, the [[California Air Resources Board]] (CARB) submitted a [[State Implementation Plan]] (SIP)<ref>'''Note: '''An SIP is a detailed plan of how to achieve a healthy ambient air quality as required  by the United States'  Clean Air Act.</ref> to the U.S. EPA as required by the United States'  [[Clean Air Act Amendments of 1990|Clean Air Act]] (CAAA).<ref name=EPA-Region9/> In accordance with that SIP, the LADWP was made responsible for installing and operating  dust control measures on Owens Lake so as to attain the ambient air quality standards of the NAAQS. As of December 17, 2010, the air quality of the Owens Valley had not yet attained the NAAQS standards for (PM<sub>10</sub>)<ref>'''Note:''' The current NAAQS for PM<sub>10</sub> is a concentration in the ambient air of no more than 150 μg/m<sup>3</sup></ref> and its non-attainment status is still designated by the EPA as being "Serious".<ref>[http://www.epa.gov/airquality/greenbk/ancl2.html Criteria Pollutant Area Summary Report] From website of the U.S. EPA Green Book. Retrieved on 2011-02-02</ref>


{{Image|Owens Valley Mitigation Veg.jpg|right|250px|Re-vegetated section of the Owens Lake area as of early 2009.}}
Air dispersion models are also used by emergency management personnel to develop emergency plans for accidental chemical releases. The results of dispersion modeling, using worst case accidental releases and meteorological conditions, can provide estimated locations of impacted areas and be used to determine appropriate protective actions. At industrial facilities in the United States, this type of consequence assessment or emergency planning is required under the Clean Air Act (CAA) codified in Part 68 of Title 40 of the Code of Federal Regulations.
{{Image|Owens Valley IBA 2010.jpg|right|250px|Owens Lake, 2010, designated as an "Important Bird Area" by the National Audubon Society}}


However, the LADWP had made much progress:<ref>[http://www.ovcweb.org/docs/Rainshadow1208.pdf  President's Message] Scroll down to the section entitled "Recent News From Owens Lake" on page 5. Retrieved from the website of the Owens Valley Committee on 2011-02-02.</ref><ref>[http://iba.audubon.org/iba/viewSiteProfile.do?siteId=213&navSite=state Owens Lake] Retrieved from an Audubon Society website on 2011-02-02.</ref><ref name=OVC-OwensLake>[http://www.ovcweb.org/OwensValley/OwensLake.html Owens Lake] Retrieved from the  website of the Owens Valley Committee on 2011-02-02.</ref><ref>[http://www.waterplan.water.ca.gov/docs/cwpu2009/0310final/v3_southlahontan_cwp2009.pdf California Water Plan: Update 2009] Bulletin 160-09, Volume 3, Regional Reports, pdf page 32.</ref>
The dispersion models vary depending on the mathematics used to develop the model, but all require the input of data that may include:


*The Owens Lake dust control project currently stretches across approximately 30 of the lake’s 100 square miles. About 27 square miles are covered with ponded water or are shallow-flooded and about  3.5 square miles have been re-vegetated and covered with native [[Saltgrass|saltgrass]] (see the adjacent photo). The saltgrass is irrigated by a water drip system.
* Meteorological conditions such as wind speed and direction, the amount of atmospheric turbulence (as characterized by what is called the "stability class"), the ambient air temperature, the height to the bottom of any inversion aloft that may be present, cloud cover and solar radiation.
*About 64 acres (0.1 square miles) has been covered with gravel as a dust control measure.
* The emission parameters such the type of source (i.e., point, line or area), the mass flow rate, the source location and height, the source exit velocity, and the source exit temperature.
*An additional 9 square miles of ponding and shallow-flooding were planned to be added by the LADWP some time in 2010.
* Terrain elevations at the source location and at receptor locations, such as nearby homes, schools, businesses and hospitals.
*The ponding and shallow-flooding have resulted in a substantial increase in wetlands habitat.
* The location, height and width of any obstructions (such as buildings or other structures) in the path of the emitted gaseous plume as well as the terrain surface roughness (which may be characterized by the more generic parameters "rural" or "city" terrain).
*Owens Lake has once again become a bird sanctuary and a major migratory stop-over site for shorebirds and waterfowl. Each spring and fall, brine flies on the lake support thousands of shorebirds, mainly sandpipers, and thousands of ducks utilize the wetlands habitat. An estimated 63,000 American Avocets stop at the lake during fall and a group of 300 to 400 Snow Geese now winters at Owens Lake.
* The [[National Audubon Society]] has now designated the Owens Lake bird sanctuary (see the adjacent photo) as  a "Important Bird Area" (IBA).


Audubon-California, Nature Conservancy, and Eastern Sierra Audubon are coordinating an effort to develop a comprehensive lake-wide management plan for Owens Lake with the cooperation and participation of LADWP, the California Department of Fish and Game, and the Great Basin Air Pollution Control District. The plan will examine the dust control project as well as springs and wetlands around the shoreline of the lake and is intended to help managers and conservation groups preserve the lake's rich wildlife resources.
Many of the modern, advanced dispersion modeling programs include a pre-processor module for the input of meteorological and other data, and many also include a post-processor module for graphing the output data and/or plotting the area impacted by the air pollutants on maps. The plots of areas impacted usually include isopleths showing areas of pollutant concentrations that define areas of the highest health risk. The isopleths plots are useful in determining protective actions for the public and first responders.


The LADWP is also scheduled to complete a plan for the long-term habitat management in the entire dust control project area. To that end, the LADWP is cooperating with the GBUAPCD, the [[California Department of Fish and  Game]], the Audubon Society and the[[ Nature Conservancy]] in the development of a long-term habitat management plan for the entire Lake Owens project. The plan will examine and consider the dust control project as well as the springs and wetlands around the shoreline of the lake and is intended to help managers and conservation groups preserve the lake's rich wildlife resources.<ref name=OVC-OwensLake/>
The atmospheric dispersion models are also known as atmospheric diffusion models, air dispersion models, air quality models, and air pollution dispersion models.


==Atmospheric layers==


Discussion of the layers in the Earth's atmosphere is needed to understand where airborne pollutants disperse in the atmosphere. The layer closest to the Earth's surface is known as the ''troposphere''. It extends from sea-level up to a height of about 18 km and contains about 80 percent of the mass of the overall atmosphere. The ''stratosphere'' is the next layer and extends from 18 km up to about 50 km. The third layer is the ''mesosphere'' which extends from 50 km up to about 80 km. There are other layers above 80 km, but they are insignificant with respect to atmospheric dispersion modeling.
The lowest part of the troposphere is called the ''atmospheric boundary layer (ABL)'' or the ''planetary boundary layer (PBL)'' and extends from the Earth's surface up to about 1.5 to 2.0 km in height. The air temperature of the atmospheric boundary layer decreases with increasing altitude until it reaches what is called the ''inversion layer'' (where the temperature increases with increasing altitude) that caps the atmospheric boundary layer. The upper part of the troposphere (i.e., above the inversion layer) is called the ''free troposphere'' and it extends up to the 18 km height of the troposphere.
The ABL is the most important layer with respect to the emission, transport and dispersion of airborne pollutants. The part of the ABL between the Earth's surface and the bottom of the inversion layer is known as the ''mixing layer''. Almost all of the airborne pollutants emitted into the ambient atmosphere are transported and dispersed within the mixing layer. Some of the emissions penetrate the inversion layer and enter the free troposphere above the ABL.
In summary, the layers of the Earth's atmosphere from the surface of the ground upwards are: the ABL made up of the mixing layer capped by the inversion layer; the free troposphere; the stratosphere; the mesosphere and others. Many atmospheric dispersion models are referred to as ''boundary layer models'' because they mainly model air pollutant dispersion within the ABL. To avoid confusion, models referred to as ''mesoscale models'' have dispersion modeling capabilities that can extend horizontally as much as  a few hundred kilometres. It does not mean that they model dispersion in the mesosphere.
==Gaussian air pollutant dispersion equation==
The technical literature on air pollution dispersion is quite extensive and dates back to the 1930s and earlier. One of the early air pollutant plume dispersion equations was derived by Bosanquet and Pearson.<ref>C.H. Bosanquet and J.L. Pearson, "The spread of smoke and gases from chimneys", ''Trans. Faraday Soc.'', 32:1249, 1936.</ref> Their equation did not assume Gaussian distribution nor did it include the effect of ground reflection of the pollutant plume.
Sir Graham Sutton derived an air pollutant plume dispersion equation in 1947<ref>O.G. Sutton, "The problem of diffusion in the lower atmosphere", ''QJRMS'', 73:257, 1947.</ref><ref>O.G. Sutton, "The theoretical distribution of airborne pollution from factory chimneys", ''QJRMS'', 73:426, 1947.</ref> which did include the assumption of Gaussian distribution for the vertical and crosswind dispersion of the plume and also included the effect of ground reflection of the plume.
Under the stimulus provided by the advent of stringent environmental control regulations, there was an immense growth in the use of air pollutant plume dispersion calculations between the late 1960s and today. A great many computer programs for calculating the dispersion of air pollutant emissions were developed during that period of time and they were commonly called "air dispersion models". The basis for most of those models was the '''Complete Equation For Gaussian Dispersion Modeling Of Continuous, Buoyant Air Pollution Plumes''' shown below:<ref name=Beychok>{{cite book|author=M.R. Beychok|title=Fundamentals Of Stack Gas Dispersion|edition=4th Edition| publisher=author-published|year=2005|isbn=0-9644588-0-2}}.</ref><ref>{{cite book|author=D. B. Turner| title=Workbook of atmospheric dispersion estimates: an introduction to dispersion modeling| edition=2nd Edition |publisher=CRC Press|year=1994|isbn=1-56670-023-X}}.</ref>
<math>C = \frac{\;Q}{u}\cdot\frac{\;f}{\sigma_y\sqrt{2\pi}}\;\cdot\frac{\;g_1 + g_2 + g_3}{\sigma_z\sqrt{2\pi}}</math>
{| border="0" cellpadding="2"
|-
|align=right|where:
|&nbsp;
|-
!align=right|<math>f</math> 
|align=left|= crosswind dispersion parameter
|-
!align=right|&nbsp;
|align=left|= <math>\exp\;[-\,y^2/\,(2\;\sigma_y^2\;)\;]</math>
|-
!align=right|<math>g</math>
|align=left|= vertical dispersion parameter = <math>\,g_1 + g_2 + g_3</math>
|-
!align=right|<math>g_1</math>
|align=left|= vertical dispersion with no reflections
|-
!align=right|&nbsp;
|align=left|= <math>\; \exp\;[-\,(z - H)^2/\,(2\;\sigma_z^2\;)\;]</math>
|-
!align=right|<math>g_2</math>
|align=left|= vertical dispersion for reflection from the ground
|-
!align=right|&nbsp;
|align=left|= <math>\;\exp\;[-\,(z + H)^2/\,(2\;\sigma_z^2\;)\;]</math>
|-
!align=right|<math>g_3</math>
|align=left|= vertical dispersion for reflection from an inversion aloft
|-
!align=right|&nbsp;
|align=left|= <math>\sum_{m=1}^\infty\;\big\{\exp\;[-\,(z - H - 2mL)^2/\,(2\;\sigma_z^2\;)\;]</math>
|-
!align=right|&nbsp;
|align=left|&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; <math>+\, \exp\;[-\,(z + H + 2mL)^2/\,(2\;\sigma_z^2\;)\;]</math>
|-
!align=right|&nbsp;
|align=left|&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; <math>+\, \exp\;[-\,(z + H - 2mL)^2/\,(2\;\sigma_z^2\;)\;]</math>
|-
!align=right|&nbsp;
|align=left|&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; <math>+\, \exp\;[-\,(z - H + 2mL)^2/\,(2\;\sigma_z^2\;)\;]\big\}</math>
|-
!align=right|<math>C</math>
|align=left|= concentration of emissions, in g/m³, at any receptor located:
|-
!align=right|&nbsp;
|align=left|&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; x meters downwind from the emission source point
|-
!align=right|&nbsp;
|align=left|&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; y meters crosswind from the emission plume centerline
|-
!align=right|&nbsp;
|align=left|&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; z meters above ground level
|-
!align=right|<math>Q</math>
|align=left|= source pollutant emission rate, in g/s
|-
!align=right|<math>u</math>
|align=left|= horizontal wind velocity along the plume centerline, m/s
|-
!align=right|<math>H</math>
|align=left|= height of emission plume centerline above ground level, in m
|-
!align=right|<math>\sigma_z</math>
|align=left|= vertical standard deviation of the emission distribution, in m
|-
!align=right|<math>\sigma_y</math>
|align=left|= horizontal standard deviation of the emission distribution, in m
|-
!align=right|<math>L</math>
|align=left|= height from ground level to bottom of the inversion aloft, in m
|-
!align=right|<math>\exp</math>
|align=left|= the exponential function
|}
The above equation not only includes upward reflection from the ground, it also includes downward reflection from the bottom of any inversion lid present in the atmosphere.
The sum of the four exponential terms in <math>g_3</math> converges to a final value quite rapidly. For most cases, the summation of the series with '''''m''''' = 1, '''''m''''' = 2 and '''''m''''' = 3 will provide an adequate solution.
<math>\sigma_z</math> and <math>\sigma_y</math> are functions of the atmospheric stability class (i.e., a measure of the turbulence in the ambient atmosphere) and of the downwind distance to the receptor. The two most important variables affecting the degree of pollutant emission dispersion obtained are the height of the emission source point and the degree of atmospheric turbulence. The more turbulence, the better the degree of dispersion.
Whereas older models rely on stability classes for the determination of <math>\sigma_y</math> and <math>\sigma_z</math>, more recent models increasingly rely on Monin-Obukhov similarity theory to derive these parameters.
==Briggs plume rise equations==
The Gaussian air pollutant dispersion equation (discussed above) requires the input of ''H'' which is the pollutant plume's centerline height above ground level. ''H'' is the sum of ''H''<sub>s</sub> (the actual physical height of the pollutant plume's emission source point) plus Δ''H'' (the plume rise due the plume's buoyancy).
[[File:Gaussian Plume.png|thumb|right|333px|Visualization of a buoyant Gaussian air pollutant dispersion plume]]
To determine Δ''H'', many if not most of the air dispersion models developed between the late 1960s and the early 2000s used what are known as "the Briggs equations." G.A. Briggs first published his plume rise observations and comparisons in 1965.<ref>G.A. Briggs, "A plume rise model compared with observations", ''JAPCA'', 15:433–438, 1965.</ref> In 1968, at a symposium sponsored by CONCAWE (a Dutch organization), he compared many of the plume rise models then available in the literature.<ref>G.A. Briggs, "CONCAWE meeting: discussion of the comparative consequences of different plume rise formulas", ''Atmos. Envir.'', 2:228–232, 1968.</ref> In that same year, Briggs also wrote the section of the publication edited by Slade<ref>D.H. Slade (editor), "Meteorology and atomic energy 1968", Air Resources Laboratory, U.S. Dept. of Commerce, 1968.</ref> dealing with the comparative analyses of plume rise models.  That was followed in 1969 by his classical critical review of the entire plume rise literature,<ref>G.A. Briggs, "Plume Rise", ''USAEC Critical Review Series'', 1969.</ref> in which he proposed a set of plume rise equations which have become widely known as "the Briggs equations".  Subsequently, Briggs modified his 1969 plume rise equations in 1971 and in 1972.<ref>G.A. Briggs, "Some recent analyses of plume rise observation", ''Proc. Second Internat'l. Clean Air Congress'', Academic Press, New York, 1971.</ref><ref>G.A. Briggs, "Discussion: chimney plumes in neutral and stable surroundings", ''Atmos. Envir.'', 6:507–510, 1972.</ref>
Briggs divided air pollution plumes into these four general categories:
* Cold jet plumes in calm ambient air conditions
* Cold jet plumes in windy ambient air conditions
* Hot, buoyant plumes in calm ambient air conditions
* Hot, buoyant plumes in windy ambient air conditions
Briggs considered the trajectory of cold jet plumes to be dominated by their initial velocity momentum, and the trajectory of hot, buoyant plumes to be dominated by their buoyant momentum to the extent that their initial velocity momentum was relatively unimportant.  Although Briggs proposed plume rise equations for each of the above plume categories, '''''it is important to emphasize that "the Briggs equations" which become widely used are those that he proposed for bent-over, hot buoyant plumes'''''.
In general, Briggs's equations for bent-over, hot buoyant plumes are based on observations and data involving plumes from typical combustion sources such as the flue gas stacks from steam-generating boilers burning fossil fuels in large power plants.  Therefore the stack exit velocities were probably in the range of 20 to 100 ft/s (6 to 30 m/s) with exit temperatures ranging from 250 to 500 °F (120 to 260 °C).
A logic diagram for using the Briggs equations<ref name=Beychok/> to obtain the plume rise trajectory of bent-over buoyant plumes is presented below:
[[Image:BriggsLogic.png|none]]
:{| border="0" cellpadding="2"
|-
|align=right|where:
|&nbsp;
|-
!align=right| Δh
|align=left|= plume rise, in m
|-
!align=right| F<sup>&nbsp;</sup> <!-- The HTML is needed to line up characters. Do not remove.-->
|align=left|= buoyancy factor, in m<sup>4</sup>s<sup>−3</sup>
|-
!align=right| x
|align=left|= downwind distance from plume source, in m
|-
!align=right| x<sub>f</sub>
|align=left|= downwind distance from plume source to point of maximum plume rise, in m
|-
!align=right| u
|align=left|= windspeed at actual stack height, in m/s
|-
!align=right| s<sup>&nbsp;</sup> <!-- The HTML is needed to line up characters. Do not remove.-->
|align=left|= stability parameter, in s<sup>−2</sup>
|}
The above parameters used in the Briggs' equations are discussed in Beychok's book.<ref name=Beychok/>
==References==
{{reflist}}
{{reflist}}
== Further reading==
*{{cite book | author=M.R. Beychok| title=Fundamentals Of Stack Gas Dispersion | edition=4th Edition | publisher=author-published | year=2005 | isbn=0-9644588-0-2}}
*{{cite book | author=K.B. Schnelle and P.R. Dey| title=Atmospheric Dispersion Modeling Compliance Guide  | edition=1st Edition| publisher=McGraw-Hill Professional | year=1999 | isbn=0-07-058059-6}}
*{{cite book | author=D.B. Turner| title=Workbook of Atmospheric Dispersion Estimates: An Introduction to Dispersion Modeling | edition=2nd Edition | publisher=CRC Press | year=1994 | isbn=1-56670-023-X}}
*{{cite book | author= S.P. Arya| title=Air Pollution Meteorology and Dispersion | edition=1st Edition | publisher=Oxford University Press | year=1998 | isbn=0-19-507398-3}}
*{{cite book | author=R. Barrat| title=Atmospheric Dispersion Modelling | edition=1st Edition | publisher=Earthscan Publications | year=2001 | isbn=1-85383-642-7}}
*{{cite book | author=S.R. Hanna and R.E. Britter| title=Wind Flow and Vapor Cloud Dispersion at Industrial and Urban Sites  | edition=1st Edition | publisher=Wiley-American Institute of Chemical Engineers | year=2002 | isbn=0-8169-0863-X}}
*{{cite book | author=P. Zannetti| title=Air pollution modeling : theories, computational methods, and available software | edition= | publisher= Van Nostrand Reinhold | year=1990 | isbn=0-442-30805-1 }}

Latest revision as of 04:25, 22 November 2023


The account of this former contributor was not re-activated after the server upgrade of March 2022.


Industrial air pollution source

Atmospheric dispersion modeling is the mathematical simulation of how air pollutants disperse in the ambient atmosphere. It is performed with computer programs that solve the mathematical equations and algorithms which simulate the pollutant dispersion. The dispersion models are used to estimate or to predict the downwind concentration of air pollutants emitted from sources such as industrial plants, vehicular traffic or accidental chemical releases.

Such models are important to governmental agencies tasked with protecting and managing the ambient air quality. The models are typically employed to determine whether existing or proposed new industrial facilities are or will be in compliance with the National Ambient Air Quality Standards (NAAQS) in the United States or similar regulations in other nations. The models also serve to assist in the design of effective control strategies to reduce emissions of harmful air pollutants. During the late 1960's, the Air Pollution Control Office of the U.S. Environmental Protection Agency (U.S. EPA) initiated research projects to develop models for use by urban and transportation planners.[1]

Air dispersion models are also used by emergency management personnel to develop emergency plans for accidental chemical releases. The results of dispersion modeling, using worst case accidental releases and meteorological conditions, can provide estimated locations of impacted areas and be used to determine appropriate protective actions. At industrial facilities in the United States, this type of consequence assessment or emergency planning is required under the Clean Air Act (CAA) codified in Part 68 of Title 40 of the Code of Federal Regulations.

The dispersion models vary depending on the mathematics used to develop the model, but all require the input of data that may include:

  • Meteorological conditions such as wind speed and direction, the amount of atmospheric turbulence (as characterized by what is called the "stability class"), the ambient air temperature, the height to the bottom of any inversion aloft that may be present, cloud cover and solar radiation.
  • The emission parameters such the type of source (i.e., point, line or area), the mass flow rate, the source location and height, the source exit velocity, and the source exit temperature.
  • Terrain elevations at the source location and at receptor locations, such as nearby homes, schools, businesses and hospitals.
  • The location, height and width of any obstructions (such as buildings or other structures) in the path of the emitted gaseous plume as well as the terrain surface roughness (which may be characterized by the more generic parameters "rural" or "city" terrain).

Many of the modern, advanced dispersion modeling programs include a pre-processor module for the input of meteorological and other data, and many also include a post-processor module for graphing the output data and/or plotting the area impacted by the air pollutants on maps. The plots of areas impacted usually include isopleths showing areas of pollutant concentrations that define areas of the highest health risk. The isopleths plots are useful in determining protective actions for the public and first responders.

The atmospheric dispersion models are also known as atmospheric diffusion models, air dispersion models, air quality models, and air pollution dispersion models.

Atmospheric layers

Discussion of the layers in the Earth's atmosphere is needed to understand where airborne pollutants disperse in the atmosphere. The layer closest to the Earth's surface is known as the troposphere. It extends from sea-level up to a height of about 18 km and contains about 80 percent of the mass of the overall atmosphere. The stratosphere is the next layer and extends from 18 km up to about 50 km. The third layer is the mesosphere which extends from 50 km up to about 80 km. There are other layers above 80 km, but they are insignificant with respect to atmospheric dispersion modeling.

The lowest part of the troposphere is called the atmospheric boundary layer (ABL) or the planetary boundary layer (PBL) and extends from the Earth's surface up to about 1.5 to 2.0 km in height. The air temperature of the atmospheric boundary layer decreases with increasing altitude until it reaches what is called the inversion layer (where the temperature increases with increasing altitude) that caps the atmospheric boundary layer. The upper part of the troposphere (i.e., above the inversion layer) is called the free troposphere and it extends up to the 18 km height of the troposphere.

The ABL is the most important layer with respect to the emission, transport and dispersion of airborne pollutants. The part of the ABL between the Earth's surface and the bottom of the inversion layer is known as the mixing layer. Almost all of the airborne pollutants emitted into the ambient atmosphere are transported and dispersed within the mixing layer. Some of the emissions penetrate the inversion layer and enter the free troposphere above the ABL.

In summary, the layers of the Earth's atmosphere from the surface of the ground upwards are: the ABL made up of the mixing layer capped by the inversion layer; the free troposphere; the stratosphere; the mesosphere and others. Many atmospheric dispersion models are referred to as boundary layer models because they mainly model air pollutant dispersion within the ABL. To avoid confusion, models referred to as mesoscale models have dispersion modeling capabilities that can extend horizontally as much as a few hundred kilometres. It does not mean that they model dispersion in the mesosphere.

Gaussian air pollutant dispersion equation

The technical literature on air pollution dispersion is quite extensive and dates back to the 1930s and earlier. One of the early air pollutant plume dispersion equations was derived by Bosanquet and Pearson.[2] Their equation did not assume Gaussian distribution nor did it include the effect of ground reflection of the pollutant plume.

Sir Graham Sutton derived an air pollutant plume dispersion equation in 1947[3][4] which did include the assumption of Gaussian distribution for the vertical and crosswind dispersion of the plume and also included the effect of ground reflection of the plume.

Under the stimulus provided by the advent of stringent environmental control regulations, there was an immense growth in the use of air pollutant plume dispersion calculations between the late 1960s and today. A great many computer programs for calculating the dispersion of air pollutant emissions were developed during that period of time and they were commonly called "air dispersion models". The basis for most of those models was the Complete Equation For Gaussian Dispersion Modeling Of Continuous, Buoyant Air Pollution Plumes shown below:[5][6]


where:  
= crosswind dispersion parameter
  =
= vertical dispersion parameter =
= vertical dispersion with no reflections
  =
= vertical dispersion for reflection from the ground
  =
= vertical dispersion for reflection from an inversion aloft
  =
           
           
           
= concentration of emissions, in g/m³, at any receptor located:
            x meters downwind from the emission source point
            y meters crosswind from the emission plume centerline
            z meters above ground level
= source pollutant emission rate, in g/s
= horizontal wind velocity along the plume centerline, m/s
= height of emission plume centerline above ground level, in m
= vertical standard deviation of the emission distribution, in m
= horizontal standard deviation of the emission distribution, in m
= height from ground level to bottom of the inversion aloft, in m
= the exponential function

The above equation not only includes upward reflection from the ground, it also includes downward reflection from the bottom of any inversion lid present in the atmosphere.

The sum of the four exponential terms in converges to a final value quite rapidly. For most cases, the summation of the series with m = 1, m = 2 and m = 3 will provide an adequate solution.

and are functions of the atmospheric stability class (i.e., a measure of the turbulence in the ambient atmosphere) and of the downwind distance to the receptor. The two most important variables affecting the degree of pollutant emission dispersion obtained are the height of the emission source point and the degree of atmospheric turbulence. The more turbulence, the better the degree of dispersion.

Whereas older models rely on stability classes for the determination of and , more recent models increasingly rely on Monin-Obukhov similarity theory to derive these parameters.

Briggs plume rise equations

The Gaussian air pollutant dispersion equation (discussed above) requires the input of H which is the pollutant plume's centerline height above ground level. H is the sum of Hs (the actual physical height of the pollutant plume's emission source point) plus ΔH (the plume rise due the plume's buoyancy).

Visualization of a buoyant Gaussian air pollutant dispersion plume

To determine ΔH, many if not most of the air dispersion models developed between the late 1960s and the early 2000s used what are known as "the Briggs equations." G.A. Briggs first published his plume rise observations and comparisons in 1965.[7] In 1968, at a symposium sponsored by CONCAWE (a Dutch organization), he compared many of the plume rise models then available in the literature.[8] In that same year, Briggs also wrote the section of the publication edited by Slade[9] dealing with the comparative analyses of plume rise models. That was followed in 1969 by his classical critical review of the entire plume rise literature,[10] in which he proposed a set of plume rise equations which have become widely known as "the Briggs equations". Subsequently, Briggs modified his 1969 plume rise equations in 1971 and in 1972.[11][12]

Briggs divided air pollution plumes into these four general categories:

  • Cold jet plumes in calm ambient air conditions
  • Cold jet plumes in windy ambient air conditions
  • Hot, buoyant plumes in calm ambient air conditions
  • Hot, buoyant plumes in windy ambient air conditions

Briggs considered the trajectory of cold jet plumes to be dominated by their initial velocity momentum, and the trajectory of hot, buoyant plumes to be dominated by their buoyant momentum to the extent that their initial velocity momentum was relatively unimportant. Although Briggs proposed plume rise equations for each of the above plume categories, it is important to emphasize that "the Briggs equations" which become widely used are those that he proposed for bent-over, hot buoyant plumes.

In general, Briggs's equations for bent-over, hot buoyant plumes are based on observations and data involving plumes from typical combustion sources such as the flue gas stacks from steam-generating boilers burning fossil fuels in large power plants. Therefore the stack exit velocities were probably in the range of 20 to 100 ft/s (6 to 30 m/s) with exit temperatures ranging from 250 to 500 °F (120 to 260 °C).

A logic diagram for using the Briggs equations[5] to obtain the plume rise trajectory of bent-over buoyant plumes is presented below:

BriggsLogic.png
where:  
Δh = plume rise, in m
F  = buoyancy factor, in m4s−3
x = downwind distance from plume source, in m
xf = downwind distance from plume source to point of maximum plume rise, in m
u = windspeed at actual stack height, in m/s
s  = stability parameter, in s−2

The above parameters used in the Briggs' equations are discussed in Beychok's book.[5]

References

  1. J.C. Fensterstock et al, "Reduction of air pollution potential through environmental planning", JAPCA, Vol. 21, No. 7, 1971.
  2. C.H. Bosanquet and J.L. Pearson, "The spread of smoke and gases from chimneys", Trans. Faraday Soc., 32:1249, 1936.
  3. O.G. Sutton, "The problem of diffusion in the lower atmosphere", QJRMS, 73:257, 1947.
  4. O.G. Sutton, "The theoretical distribution of airborne pollution from factory chimneys", QJRMS, 73:426, 1947.
  5. 5.0 5.1 5.2 M.R. Beychok (2005). Fundamentals Of Stack Gas Dispersion, 4th Edition. author-published. ISBN 0-9644588-0-2. .
  6. D. B. Turner (1994). Workbook of atmospheric dispersion estimates: an introduction to dispersion modeling, 2nd Edition. CRC Press. ISBN 1-56670-023-X. .
  7. G.A. Briggs, "A plume rise model compared with observations", JAPCA, 15:433–438, 1965.
  8. G.A. Briggs, "CONCAWE meeting: discussion of the comparative consequences of different plume rise formulas", Atmos. Envir., 2:228–232, 1968.
  9. D.H. Slade (editor), "Meteorology and atomic energy 1968", Air Resources Laboratory, U.S. Dept. of Commerce, 1968.
  10. G.A. Briggs, "Plume Rise", USAEC Critical Review Series, 1969.
  11. G.A. Briggs, "Some recent analyses of plume rise observation", Proc. Second Internat'l. Clean Air Congress, Academic Press, New York, 1971.
  12. G.A. Briggs, "Discussion: chimney plumes in neutral and stable surroundings", Atmos. Envir., 6:507–510, 1972.

Further reading

  • M.R. Beychok (2005). Fundamentals Of Stack Gas Dispersion, 4th Edition. author-published. ISBN 0-9644588-0-2. 
  • K.B. Schnelle and P.R. Dey (1999). Atmospheric Dispersion Modeling Compliance Guide, 1st Edition. McGraw-Hill Professional. ISBN 0-07-058059-6. 
  • D.B. Turner (1994). Workbook of Atmospheric Dispersion Estimates: An Introduction to Dispersion Modeling, 2nd Edition. CRC Press. ISBN 1-56670-023-X. 
  • S.P. Arya (1998). Air Pollution Meteorology and Dispersion, 1st Edition. Oxford University Press. ISBN 0-19-507398-3. 
  • R. Barrat (2001). Atmospheric Dispersion Modelling, 1st Edition. Earthscan Publications. ISBN 1-85383-642-7. 
  • S.R. Hanna and R.E. Britter (2002). Wind Flow and Vapor Cloud Dispersion at Industrial and Urban Sites, 1st Edition. Wiley-American Institute of Chemical Engineers. ISBN 0-8169-0863-X. 
  • P. Zannetti (1990). Air pollution modeling : theories, computational methods, and available software. Van Nostrand Reinhold. ISBN 0-442-30805-1.