Injection height for biomass burning emissions from boreal forest fires Fok-Yan Leung

Injection height for biomass burning emissions from boreal forest fires

• Fok-Yan Leung

• April 12, 2007.

• Harvard University

• Special thanks to:

• Jennifer Logan, Rokjin Park, and Dominic Spracklen (Harvard)

• Edward Hyer and Eric Kasischke (UMD)

• Leonid Yurganov

• David Diner, Dominic Mazzoni, David Nelson, and Ralph Kahn (NASA/JPL)

• Funding from the NSF and EPA

We began by looking at emissions estimates for 1998 boreal fires, which vary significantly.

Comparison of surface and column data from 1998 with results of GEOS-Chem simulations

During intense boreal fires, intense heat can result in lofting of emissions well above the boundary layer

Putting large fraction of emissions in free troposphere reconciled model results with both surface and column data

• Putting large fraction of emissions in free troposphere reconciled model results with both surface and column data

• For both surface and column data (anomaly data): KAS05 seems to perform better in capturing the CO behavior using parameterization

Conclusions from study of 1998 study: Injection of biomass burning emissions in the free troposphere are necessary to reconcile ground and column data Injection of biomass burning emissions in free troposphere results in higher tropospheric ozone throughout the northern hemisphere due to longer sequestration of NOx by PAN formation. Preliminary studies suggest that model results are not particularly sensitive to the exact fractional split of emissions (Turquety et al., [2007], and unpublished work) We were motivated to move beyond the “sensitivity analysis” level, and to our ongoing study of plume heights using the MISR instrument

Using the MISR Instrument

We look at discrete plumes from fires. Algorithm of Mazzoni and Nelson:

• We look at discrete plumes from fires. Algorithm of Mazzoni and Nelson:

• Detects plumes by trained plume shape recognition algorithm
• Uses MODIS hotspots to narrow down number of plumes
• Determines the maximum plume height

• Using the algorithm, 66 discrete plumes were found in Alaska and Northern Canada during summer of 2002

• Example of algorithm at work…

Kahn et al, [2006] observed clear relationship between atmospheric stability and observed plume heights. We compare the stability profiles calculated using the coarser GEOS4 data with the finer resolution BRAMS data, at the 66 sites

• GEOS4

• 1°x1.25° horizontal resolution

• 30 vertical levels in troposphere

• Reanalysis product.

• Pressure is instantaneous pressure

• Temperature is 6 hour average

Stability profiles: BRAMS (courtesy Marcos Longo) and GEOS4: “Neutral” profiles and profiles with regions of high stability

• In general, same vertical structural characteristics in BRAMS and GEOS4 stability profiles
• Vertical “Offset” between GEOS4 and BRAMS profiles

Example of plume in trapped in a layer of high stability

Directions: Moving forward

• Preliminary results suggests that stability profiles may provide a way to parameterize injection heights

• First pass analysis of coarser grid, 2°x2.5° data show very similar stability profiles to those calculated using data from 1°x1.25° grid.

• Ultimately interested in relationship between plume heights and height of diffuse smoke

• We are moving towards a parameterization for injection heights of emissions from boreal forest fires in GEOS-Chem

S1: Modeling fire plumes is actually a quite well defined problem

• Essentially plume rise is governed by the characteristics of the fire itself (rate fuel consumption determines buoyant energy) and on local meteorology (wind direction, convection, stability of the atmosphere)

• However, challenge is in parameterizing plume injection height on a coarse grid

• In a coarse grid

• Meteorology is averaged over a large geographical area.
• Other factors are highly uncertain at best (e.g. fuel loading)

• Need a statistical method, preferably one that can be done online during simulations

S2: Effect of PAN on ozone chemistry

Implications for ozone chemistry – the effect of PAN carried aloft.

S3: Comparing 66 plume histograms to stability profiles derived from GEOS4 data:

S4: Comparing 66 plume histograms to stability profiles derived from BRAMS data:

S5:Comparing 66 plume histograms to stability profiles derived from BRAMS and GEOS4 data: BRAMS/GEOS4

S6: GEOS4 vs. BRAMS stability profiles

• Generally, similar vertical structures

• More levels of high stability at lower latitudes in GEOS4

• Tropopause tends to be higher in BRAMs data

• Both stability profiles calculated by simple forward method – however, BRAMS has higher vertical resolution (150m) and horizontal resolution (40km)

• Assumption of standard US atmosphere in calculation of GEOS4 data

• Difference in terrain levels not sufficient to account for “vertical shift”