Development of a Lightning NOx Algorithm for WRF-Chem

1. Development of a Lightning NOx Algorithm for WRF-Chem Amanda Hopkins Hansen Department of Meteorology Florida State University ahopkins@met.fsu.edu Advisor: Henry E. Fuelberg Florida State University

2. OUTLINE • 1. What is Lightning NOx (LNOx)? • 2. Rationale - why is LNOx important ? • Lightning Parameterization • Previous Research • FSU LNOx Research Plan-necessary ingredients • 4. Validation

3. O2 N2 N2 N2 N2 N2 O2 N2 O2 O2 N2 N2 • High temps and pressure break up N2 and O2 • As the hot channel cools from the outside in, • NO is produced • the rate of cooling determines concentration • NO2 produced during the return stroke N2 Air is ~78%N2 and ~21% O2 N2 + O2 NO + NO Zel’dovich mechanism (Zel-dovich and Raizer 1966) NOx(NO+NO2)

4. Rationale - why is LNOx important ? June, July, and August 1. Lightning is an important source of NOx in the relatively clean upper troposphere Dec, Jan, and Feb Pressure (hPa) Zonal Mean Lightning NOx production (10-2kg/s) From the NASA GISS GCM Hopkins, 2003

5. Rationale - Why is LNOx important ? • 2. LNOx indirectly affects local air quality and global climate • Has a strong influence on Ozone (O3) and hydroxyl • radical (OH) concentration NOx [NO+NO2]: • Is a primary pollutant found in photochemical smog • Is a precursor for tropospheric ozone formation TROPOSPHERIC OZONE: • Is the third most important • greenhouse gas • Impacts the Earth’s radiation budget • and can cause changes in • atmospheric circulation patterns • Is toxic to humans, plants and • animals Photo courtesy of University of California at Berkley

6. Global NOx Budget 3. LNOx is difficult to model realistically, but is Important in the global budget of NOx Figures from Global Emissions Inventory Activity data sets Note: lightning figures are calculated by model internally

7. An LNOx algorithm for WRF-Chem with Parameterized Convection What are the necessary ingredients? • Flash rate parameterization • 2) Average NO production per flash • 3) Method of specifying the vertical distribution • of LNOx emissions

8. Schemes Often Used for Flash Rate Production 1. Cloud Top Height [Price and Rind, 1992; Michalon et al., 1999] 2. Convective Precipitation Rate [Meijer et al., 2001; Allen and Pickering, 2002] Separate parameterizations are needed to simulate land and oceanic lightning 3. Upward Cloud Mass Flux/Updraft Velocity [Grewe et al., 2001; Allen and Pickering, 2002] • This scheme cannot be developed based on observations. Model output must be used. Therefore, the relationship between these variables and flashrate is model-dependent. None of these schemes yield global flash distributions that compare well with satellite flash observations. What is needed is a more microphysically-based scheme. Possibilities include the schemes by Deierling et al. (2005) and Futyan and Del Genio (2007).

9. Flux Hypothesis Deierling et al., 2005 They propose that lightning frequency, f, is proportional to the product of the precipitation rate, p, and the mass upflux of ice crystals, Fi. The charging mechanism used in their study involves rebounding collisions between heavy graupel pellets and the lighter ice crystals . This scheme would work very well on the cloud-scale We need a Regional scale relationship!!! After Saunders et al. 1993

10. Convective Radar Storm Height Futyan and Del Genio, 2007 Cloud top height may be several kilometers higher than the height of significant radar signal Fifth order power law for radar top height (above surface) Second order for radar top height above 0 degree isotherm.

11. LNOx algorithm for WRF-Chem What are the necessary ingredients? Lightning Flash Rate Parameterization • Flash rate (F) can be calculated from the following relation: • F=AHn • Relationship between Radar storm height above the freezing level • and LDAR data is needed to determine A and n above. • Radar Reflectivity calculated within WRFChem using hydrometeors • from microphysical scheme (WSM6) coupled with the Kain-Fritsch • cumulus parameterization. • Convective storm height (H) is found using the 20 dBz contour and t • the freezing level is obtained from the WRF temperature field.

12. WSR-88 Doppler Radar and the Lightning Detection and Ranging (LDAR) data Futyan and Del Genio (2007) used TRMM data for their research: We will use Doppler radar and LDAR data to form a relationship between flash rate and radar storm height. 3-D Lightning Mapping Locations Kennedy Space Center, Florida** - LDAR/Vaisala Huntsville, Alabama – NM Tech Dallas, Texas – LDAR/Vaisala Washington D.C. – NM Tech **We will start with KSC

13. WSR-88 Doppler Radar and the Lightning Detection and Ranging (LDAR ) data from Kennedy Space Center (KSC) Doppler Radar Melbourne, Florida KSC LDAR. Sensor 0 is the central LDAR receiver. (After Poehler and Lennon 1979 and Vollmer 2002)

14. WSR-88 Doppler Radar and the Lightning Detection and Ranging (LDAR ) Data from Kennedy Space Center (KSC) The location of the Melbourne, Florida National Weather Service Office WSR-88D radar in relation to the Kennedy Space Center. Range rings are provided at 10 km intervals. (After Nelson 2002).

15. LDAR Flash Algorithm Created by Nelson (2002), and used recently by Stano (2007) Turns “sparks” into flashes • Warm season months of May through September 1997-2005 For a spark to be included in a flash it must meet Temporal and Spatial Constraints: Occurs within 3 s of the first observed spark Each successive spark of a flash had to be within 0.5 s of the previous spark. Occurs within the ellipse if included in the flash. LDAR is 97% accurate out to 100 km And 99% accurate to 25 km. (Applied Meteorology Unit 1996) Nelson 2002

16. WDSS: Warning Decision Support System This will provide a way to visualize lightning flashes with radar reflectivity The image depicts individual Lightning Detection and Ranging (LDAR) sparks as measured by the Kennedy Space Center LDAR network, cloud-to ground strikes measured by the National Lightning Detection Network (NLDN), and quality controlled reflectivity data from the Melbourne NWS radar. From the LDAR sparks we can calculate flash rate.

17. WDSS: Warning Decision Support System Cross section view of cell 49 in previous slide

18. Composite Radar Reflectivity • from KMLB for 3 Aug 2005 • Black contours are LDAR flash • counts in the corresponding • 6 minute radar scan • Red rings are distance (km) from • KMLB radar located: • Lat=28.109 • Lon=-80.650 • LDAR central receiver location: • Lat=28.5386 • Lon=-80.6431

19. Cross section through approx • 28.4N, Radar Reflectivity • from KMLB for 3 Aug 2005 • Black contours are LDAR flash • counts in the corresponding • 6 minute radar scan

20. LNOx algorithm for WRF-Chem Flash rate (F) may be calculated based on the following relationship: F=AHn The plot below shows a relationship between WSR88D radar storm height And LDAR flash rate on 1km2 bins NOTE: FD07 plot was normalized over TRMM scan time and 300km2 This is a much larger area than mine and a much longer time scale

21. LNOx algorithm for WRF-Chem What are the necessary ingredients? Lightning Flash Rate Parameterization • Flash rate (F) can be calculated from the following relation: • F=AHn • Relation between Radar storm height above the freezing level • and LDAR data is needed to determine A and n above. • Radar Reflectivity calculated within WRF-Chem using hydrometeors • from microphysical scheme (WSM6) coupled with the Kain-Fritsch • cumulus parameterization. • Convective storm height (H) is found using the 20 dBz contour, and • the freezing level is obtained from the WRF temperature field.

22. LNOx algorithm for WRF-Chem what are the necessary ingredients? LNOx parameterization • Production rate of NO from both IC and CG lightning: • -500 moles per flash (Ott et al.,in progress) • Vertical Distribution of NO: • Pickering et al., 1998: Wind fields from Goddard Cumulus • Ensemble (GCE) model were used to redistribute LNOx throughout • the duration of the storm. Profiles were constructed for mid-latitude • continental, tropical continental, and tropical marine regimes • based on profiles computed for individual storms in each regime

23. WRF-Chem “Online” Chemical Model • Consistent: All transport done by meteorology • model • * Same vertical and horizontal coordinates (no • horizontal and vertical interpolation) • * Same physics parameterization for subgrid scale • * No interpolation in time • Easy Data Handling • Makes computation time faster

24. Directly Involved in Major WRF-Chem Development NOAA/ESRL Georg Grell, Steven Peckham, Stuart McKeen PNNL Jerome Fast, Bill Gustafson, Rahul A. Zaveri, James C. Barnard NCAR Bill Skamarock Rainer Schmitz (University of Chile – Santiago, Chile) Marc Salzmann (Max Planck Institute for Chemistry – Mainz, Germany) And Many more national and international collaborators About 250 registered users

25. WRF-Chem Structure http://ruc.fsl.noaa.gov/wrf/WG11/wrf_wksp_2007_minitutorial/02_totorial.wrf-chem.Peckham.pdf

26. WRF-Chem V2.2 Chemistry • Chemical mechanisms: • RADM2, Carbon Bond (CBMZ) • Photolysis (coupled with hydrometeors and aerosols): • Madronich, Fast-j (coupled to aerosols and microphysics) • Deposition: • Dry deposition (coupled with soil/veg scheme, “flux-resistance” analogy) • Simplified wet deposition by convective parameterization Biogenic emissions: • Guenther – online calculation based on USGS land use, T and radiation • BEISv3.11 (modify reference fields produced from complex land use data)

27. Comparison with 350 Ozone monitors during ICARTT/NEAQS experiment – 8hr peak Evaluation period: July 12 – July 31, 2004 Continuous model improvement over the last 2 years

28. Observed ozone vs model for the Thompson Farm surface station from 13-20 July 2002 1 hour max O3 for all EPA AIRNOW Surface stations from 21 July to 4 Aug 2004. Grell et al. 2005

29. Transects From STERAO 1996 Dye et al. 2000

30. Transects From STERAO 1996 Dye et al. 2000

31. WRF-Chem Model and Domain for FSULNOx Mother domain: 36km Nested domain: 12km • Within the WRF-Chem namelist two • User defined options will be added • lighting_opt 0,1 ; turn off, on lightning • Lightning_time_step Lightning Module will be placed in Phys/ directory and contains subroutine radar subroutine lightning Module will be called from dyn_em/solve_em.F right before Chemistry is called.

32. WRF-Chem Physics Microphysics-WSM6 Cumulus Parameterization-Kain Fritsch eta Long wave & Short wave radiation-RRTM/Dudhia Planetary Boundary Layer-Mellor-Yamada-Janjic (MYJ) Lightning Parameterization will likely be very sensitive to the Microphysics scheme chosen. We may need to perform a sensitivity Study to make sure WSM6 is the scheme to use.

33. Verification Simulations for Summer 2004 INTEX-NA period over eastern US and comparisons with aircraft data DC-8 flight track showing the path through the Huntsville, AL storms We have the same data for Kennedy Space Center Plot courtesy of Mike Porter

34. Verification SCIAMACHY satellite data SCanning Imaging Absorption SpectroMeter for Atmospheric CHartographY European Space agency Satellite in orbit in 2004 Primary mission to measure atmospheric trace gases Including tropospheric ozone and NO2 Measurements at relatively high resolution (0.2 nm to 0.5 nm) over a wide range: 240 nm to 1700 nm

35. Verification IONS-INTEX ozonesonde Network Study Since NOx is a precursor for tropospheric ozone formation; Model ozone will be a good indicator of NOx model performance Ozonesondes launched during the 2004 INTEX field campaign 10 Operational Sites in Eastern US Including Huntsville, AL

36. Summary • Flash rate parameterization being developed for WRF-Chem usingobserved radar and 3-D lightning mapping array data. • WRF-Chem first will be tested with the Futyan and Del Genio (2007) relationship and then with the newly developed scheme. • Existing NO production per flash and vertical distribution information will be used. • Testing of WRF-Chem with lightning will be conducted using aircraft NOx observations from the ICARTT (NASA and NOAA data) experiment from Summer 2004, ozonesondes, and possibly NO2 from satellite.

37. Additions in WRF/Chem V2.2 • Kinetic Pre-Processor (KPP), MPI Mainz • Improved convective (non-resolved) transport, coupling of convective parameterization with atmospheric and photolysis radiation (ESRL/GSD) • Non-resolved and resolved aqueous phase chemistry, wet deposition, (NOAA/ARL/EPA, PNNL,ESRL/GSD) • 2-way nesting (PNNL, ESRL/GSD) • Cloud-aerosol interaction (indirect effect) with Lin et al. 6-class microphysics scheme (PNNL) • Lateral boundary conditions from global models (U of Chile, ESRL/GSD) • Urban parameterizations (NCAR, coming soon Spain) • Positive definite advection (NCAR) • NMM and ARW dynamic cores (ESRL/GSD) • Offline version for the ARW core (to be released shortly (C-DAC, India and ESRL/GSD)