AJ Designer

Atmospheric Dispersion Model Calculator

Effective stack height equals physical stack height plus plume rise

Solution

—
Share:

Gaussian Plume Model

The Gaussian plume model estimates pollutant concentration at any downwind point by assuming the plume spreads in a normal distribution both horizontally and vertically. It requires emission rate, wind speed, dispersion coefficients, and effective stack height.

C(x,y,z) from Q, u, σy, σz, H

Effective Stack Height

Effective stack height combines the physical chimney height with the extra rise the hot plume gains from buoyancy and momentum.

H = hp + Δh

Wind Speed at Elevation

The power law wind profile extrapolates a wind speed measured at one height to any other elevation. Widely used in air quality modeling and wind energy assessments.

u = u₀(h/h₀)ⁿ

How It Works

Atmospheric dispersion models predict how pollutants spread downwind from a source like a smokestack. The Gaussian plume model, developed by Pasquill, calculates the concentration of a contaminant at any point in space by combining the emission rate, wind speed, and the horizontal and vertical spread of the plume (σy and σz). The effective stack height combines the physical chimney height with the extra rise the hot plume gains from buoyancy and momentum. Plume rise formulas differ by atmospheric stability: superadiabatic (unstable), neutral, and subadiabatic (inversion) conditions each have distinct coefficients. The power-law wind profile adjusts a measured wind speed at one elevation to any other height.

Example Problem

A power plant has a 60 m physical stack and the plume rise is estimated at 25 m. What is the effective stack height?

  1. Identify the knowns. Physical stack height hp = 60 m, plume rise Δh = 25 m (estimated from buoyancy and momentum at the stack exit).
  2. Identify what we're solving for. We want the effective stack height H — the elevation used as the release point in the Gaussian dispersion model.
  3. Write the effective stack height equation: H = hp + Δh. The effective height captures how high the hot, buoyant plume actually rises above the chimney before traveling downwind.
  4. Substitute the values: H = 60 m + 25 m.
  5. Simplify the arithmetic: 60 + 25 = 85.
  6. **The effective stack height is H = 85 m** — this is the release elevation fed into the Gaussian plume equation, not the 60 m physical chimney height.

When to Use Each Variable

  • Gaussian Plume: Point in Space — when you need the pollutant concentration at a specific (x, y, z) location downwind of a stack, e.g., assessing exposure at a nearby school.
  • Gaussian Plume: Ground Level — when you need concentration at ground level (z = 0) at a lateral offset y, e.g., predicting worst-case exposure along a property boundary.
  • Gaussian Plume: Plume Centerline — when you need the maximum concentration directly downwind at plume height, e.g., determining the peak concentration footprint.
  • Gaussian Plume: Ground Source — when the emission comes from ground level with no stack, e.g., modeling fugitive dust from a construction site.
  • Solve for Effective Stack Height — when you know the physical chimney height and plume rise, e.g., combining stack measurements with meteorological data for a permit application.
  • Solve for Wind Speed at Elevation — when you need to extrapolate a surface wind measurement to stack-top height, e.g., adjusting weather station data for dispersion modeling.
  • Solve for Plume Rise — when you know the stack exit conditions and need to estimate how high the plume rises above the stack, e.g., calculating effective release height for an EIS.

Key Concepts

The Gaussian plume model assumes pollutant concentrations follow a bell-curve distribution in both the horizontal and vertical directions as the plume travels downwind. Effective stack height — the sum of physical stack height and buoyancy-driven plume rise — determines the initial release elevation. Atmospheric stability (classified from A through F by Pasquill) controls how quickly the plume disperses: unstable conditions spread pollutants rapidly, while stable inversions trap them near the surface.

Applications

  • Air quality permitting: demonstrating compliance with National Ambient Air Quality Standards (NAAQS) for new industrial sources
  • Environmental impact assessment: predicting downwind pollutant concentrations from power plants, refineries, and incinerators
  • Emergency response: estimating chemical plume extent after an accidental release at a chemical plant
  • Urban planning: siting schools, hospitals, and residential developments at safe distances from major emission sources
  • Wind energy: using the power-law wind profile to estimate wind resources at turbine hub height from ground-level weather data

Common Mistakes

  • Using ground-level wind speed directly in the Gaussian equation instead of adjusting it to stack height with the power-law profile — this overestimates concentrations because the actual wind speed at release height is usually higher
  • Confusing physical stack height with effective stack height — ignoring plume rise can overpredict ground-level concentrations by a large margin, especially for hot buoyant plumes
  • Applying the wrong stability class — using neutral stability when conditions are actually stable can underpredict peak concentrations by a factor of 5 to 10
  • Forgetting that the Gaussian model assumes steady-state, uniform wind — it is not valid for calm winds, complex terrain, or rapidly shifting meteorological conditions

Frequently Asked Questions

What is effective stack height in air pollution modeling?

Effective stack height H is the physical chimney height plus the plume rise produced by buoyancy and exit momentum. It is the elevation Gaussian models use as the release point — a 60 m stack with 25 m of plume rise behaves like a 85 m release, dramatically reducing predicted ground-level concentrations.

How does atmospheric stability change dispersion behavior?

Unstable (superadiabatic) air mixes pollutants vertically and dilutes them quickly. Stable (subadiabatic or inversion) layers suppress vertical motion and can trap a plume near the surface, producing fumigation events. The stability parameter n in the wind-profile law and the Pasquill class (A–F) capture this effect.

How does the Gaussian plume model estimate concentrations?

It assumes the pollutant cloud spreads as a bivariate normal distribution around the plume centerline as it moves downwind. Concentration at any (x, y, z) point depends on emission rate Q, wind speed u, lateral and vertical spread σy and σz, and the effective release height H. A ground-reflection term accounts for pollutants that bounce off the surface.

What is the difference between superadiabatic, neutral, and subadiabatic plume rise?

Superadiabatic (unstable) air gives the largest Δh because convection adds to buoyant rise — coefficients are 3.47 and 5.15. Neutral conditions use smaller coefficients (0.35 and 2.64), and subadiabatic (inversion) conditions suppress rise further with 1.04 and 2.24 and a sign change on the velocity term. Each regime is fit separately to plume-rise data.

Why is the power-law wind profile used in dispersion modeling?

Weather stations measure wind at ~10 m, but the dispersion equation needs wind at stack-top height. The power law u = u₀ × (h/h₀)ⁿ extrapolates upward, with the exponent n encoding stability (roughly 0.10 for unstable, 0.15 neutral, 0.30+ for stable). It is much simpler than a full logarithmic boundary-layer profile.

What units does the Gaussian model use for emission rate?

Q is the mass emission rate in grams per second (g/s). When wind speed is in m/s and σy, σz, H are in meters, the resulting concentration C comes out in g/m³ — convert to µg/m³ for comparison with most air quality standards (the NAAQS for fine particles, for example, is 35 µg/m³ over 24 hours).

What are the main limitations of the Gaussian plume approach?

It assumes steady-state, uniform wind, level terrain, no chemical reactions, and continuous emissions. It fails for calm winds (u below ~1 m/s), complex terrain (mountains, urban canyons), short-duration releases, and reactive pollutants. For those cases, regulators use models like AERMOD, CALPUFF, or computational fluid dynamics.

Worked Examples

Air-Quality Permitting

What is the effective stack height for an 80 m physical stack with 20 m of plume rise?

A combustion source has a physical chimney height of 80 m above grade and the buoyant plume rises an additional 20 m above the stack top before leveling off. What effective stack height H should be used as the release elevation in the Gaussian plume model?

  • Knowns: hp = 80 m, Δh = 20 m
  • H = hp + Δh
  • H = 80 + 20

H = 100 m

EPA modeling guidance (AERMOD, ISC) uses effective stack height — not the physical chimney height — because the buoyancy-driven rise lifts the centerline well above the stack tip before the plume bends over and travels horizontally.

Met-Tower Wind Extrapolation

What is the wind speed at 60 m if 4 m/s is measured at the 10 m anemometer?

A meteorological tower in rural open terrain records 4 m/s at the standard 10 m anemometer height. Using the power-law profile with stability exponent n = 0.143 (typical neutral, open country), what is the corresponding wind speed at the 60 m stack-top elevation?

  • Knowns: u₀ = 4 m/s, h₀ = 10 m, h = 60 m, n = 0.143
  • u = u₀ × (h / h₀)^n
  • u = 4 × (60 / 10)^0.143
  • u = 4 × 6^0.143 = 4 × 1.292

u ≈ 5.17 m/s at 60 m

The stability exponent n depends on roughness and atmospheric stability — common values are about 0.10 over water, 0.14 in rural open country (used here), and 0.22–0.35 in built-up urban terrain. EPA stability classes A–F map to specific n values for permitting work.

Ground-Level Receptor — Centerline Concentration

What is the ground-level centerline concentration from a Q = 10 g/s stack with H = 50 m?

A continuous point source emits 10 g/s of a non-reactive pollutant at an effective stack height of H = 50 m. Under the wind and stability conditions of interest, the lateral and vertical dispersion coefficients are σy = 20 m and σz = 15 m, with u = 5 m/s. What is the predicted ground-level concentration along the plume centerline?

  • Knowns: Q = 10 g/s, u = 5 m/s, σy = 20 m, σz = 15 m, H = 50 m
  • C = Q / (π × u × σy × σz) × exp(−H² / (2σz²))
  • C = 10 / (π × 5 × 20 × 15) × exp(−2500 / 450)
  • C = 0.002122 × exp(−5.556) = 0.002122 × 0.003866

C ≈ 8.2 × 10⁻⁶ g/m³ (≈ 8.2 μg/m³)

Real permitting models (AERMOD, CALPUFF) replace this single-formula evaluation with hour-by-hour calculations over a year of meteorology, but the underlying Gaussian point-source form is the conceptual basis. Conservative back-of-envelope checks like this are common screening tools.

Atmospheric Dispersion Formulas

Several related equations make up a Gaussian dispersion model — the plume itself, the effective stack height, the wind-profile extrapolation, and three plume-rise regimes by atmospheric stability:

C(x,y,z) = Q / (2π u σy σz) · exp(−y²/2σy²) · [exp(−(z−H)²/2σz²) + exp(−(z+H)²/2σz²)]Gaussian plume — concentration at point (x, y, z)
H = hp + ΔhEffective stack height
u = u0 × (h / h0)nPower-law wind profile
Δhsuper = 3.47 Vs d / u + 5.15 √Qh / uPlume rise — superadiabatic (unstable)
Δhneutral = 0.35 Vs d / u + 2.64 √Qh / uPlume rise — neutral stability
Δhsub = −1.04 Vs d / u + 2.24 √Qh / uPlume rise — subadiabatic (inversion)

Where:

  • C(x, y, z) — downwind pollutant concentration (g/m³)
  • Q — mass emission rate from the source (g/s)
  • u — wind speed at stack height (m/s)
  • σy, σz — horizontal and vertical dispersion coefficients (m), functions of downwind distance and Pasquill stability class
  • y, z — crosswind and vertical receptor coordinates (m), measured from the plume centerline and ground
  • H — effective stack height (m); release point in the model
  • hp — physical stack (chimney) height (m)
  • Δh — plume rise above the chimney from buoyancy and momentum (m)
  • u0, h0 — reference wind speed (m/s) and elevation (m), typically the weather-station 10 m measurement
  • n — stability parameter — ~0.10 unstable, ~0.15 neutral, ~0.30+ stable
  • Vs — stack exit velocity (m/s)
  • d — stack exit diameter (m)
  • Qh — sensible heat emission rate at the stack exit (kW)

These formulas assume a continuous, steady-state plume in a uniform wind field over flat terrain — the assumptions baked into screening models like SCREEN3. Refined regulatory models (AERMOD, CALPUFF) use the same Gaussian core but layer on boundary-layer physics, building downwash, and time-varying meteorology.

Related Calculators

Related Sites