Air Force Weather Agency Integrity Service Excellence Probabilistic

Air Force Weather Agency Integrity - Service - Excellence Probabilistic Lightning Forecasts Using Deterministic Data Evan Kuchera and Scott Rentschler 16 Aug 2007

Motivation n Air Force operators require skillful and objective probabilistic weather information to maximize efficiency and minimize loss n Typically this is accomplished with ensembles for grid scale phenomena n However, sub grid scale processes are probabilistic in nature even with deterministic data n We believe that ensemble forecast skill will be higher if a probabilistic approach is taken with each ensemble member for sub grid scale phenomena n Addresses both sub-grid scale and flow uncertainties Integrity - Service - Excellence 2

Motivation n Example—lightning forecast with SPC SREF method n 10 ensemble members n CAPE values of 130, 125, 120, 115, 110, 105, 103, 102, 101 n With a forecast threshold of 100 J/kg, this gives a 100% chance of lightning n However, with values so close to the threshold, the true probability is likely much closer to 50% than 100% n This can be accounted for somewhat with real-time calibration after the ensemble is created (as SPC does with success), but this is not necessarily an option for the Air Force (resource constraints, lack of calibration data) Integrity - Service - Excellence 3

Background n Lightning background: Need graupel and ice particle collisions to transfer negative charge to the larger particles n Thunderstorm updrafts need to grow large graupel particles with enough fall speed to cause a separation of charge in the vertical n The theoretical value of CAPE required to do this is only 25 J/kg n Integrity - Service - Excellence 4

Background n CAPE background: Accepted parcel theory assumption is that as the parcel rises, all condensate is immediately removed, and that there is no latent heat of freezing n However, lightning is caused by frozen condensates in an updraft! n We decided to test CAPE both ways—the traditional way, and with condensates/latent heat of freezing n Integrity - Service - Excellence 5

Image from NASA-GHCC Worldwide lightning climatology Integrity - Service - Excellence 6

Traditional Lifted Index Integrity - Service - Excellence 7

TEST Lifted Index Integrity - Service - Excellence 8

Methodology n Goal: create a probabilistic lightning algorithm using a large set of CONUS observations and physical assumptions relevant worldwide n 2006 3 -hourly 20 km RUC analyses NLDN lightning in the RUC grid box (0 -3 hr after analysis) n 3 hour precipitation from METARS n n Find which forecast parameters are the best, then curve fit the probability of lightning given a binned value of that parameter Integrity - Service - Excellence 9

2006 Results Integrity - Service - Excellence 10

NLDN 3 -hourly lightning climatology for a 16 km grid box (2003 -2006) Integrity - Service - Excellence 11

Results GL CAPE is calculated from the LFC to -20 C Set to zero if equilibrium level is warmer than -20 C TEST is condensate and latent heat of freezing included Integrity - Service - Excellence 12

CAPE > 0, Precipitation > 0. 01 Integrity - Service - Excellence 13

CAPE=0, Precipitation > 0. 01 Integrity - Service - Excellence 14

CAPE > 0, Precipitation=0 Integrity - Service - Excellence 15

Results Climatology=0. 155 Integrity - Service - Excellence 16

Results Perfect Reliability Integrity - Service - Excellence 17

Results No Skill Forecast Integrity - Service - Excellence 18

Results SPC method: forecast 100% chance of lightning if GL CAPE is greater than 100 J/kg and precipitation is greater than 0. 01 inches. Forecast 0% otherwise. NULL method: Always forecast 0% chance of lightning. TEST method: Algorithm presented here. BSS: Brier skill score, compares mean squared error of forecast to mean squared error of climatology. 1 is perfect, 0 is no skill, negative is worse than climatology. ROC area: Total integrated area underneath ROC curve. 1 is perfect, 0. 5 is no skill. Integrity - Service - Excellence 19

Summary n Algorithm has been developed to forecast lightning probability given observed instability (RUC analysis) and precipitation (METARS) n Algorithm is somewhat sharp, reliable at all forecast probabilities, and has good resolution of events and non-events n Buoyancy calculations probably need to account for condensate and latent heat of freezing—but our data are not conclusive on this point Integrity - Service - Excellence 20

Other/Future Work n Equations have been developed (not shown here) to forecast strikes per unit area for application to any model resolution n After knowing strikes per unit area, can forecast probabilities for smaller areas (i. e. Air Force base warning criteria area) based on downscaling climatology—equation has been developed for this purpose as well n Just beginning to look at algorithm with model data and in ensembles—issues with model precipitation forecasts n Acknowledgments: ARM data archive, Dr. Tony Eckel, Stephen Augustyn, Bill Roeder, Dr. David Bright, Jeff Cunningham Integrity - Service - Excellence 21

Questions? GFS 66 hour grid point lightning probability forecast valid this afternoon Integrity - Service - Excellence 22

Backup Slides Integrity - Service - Excellence 23

Backup Slides n Adjustments for changes in model resolution or area of interest First, re-calculate total number of strikes for the new model grid box area n If model grid is finer than RUC, re-calculate probabilities using inverse of strikes equation n If model grid is coarser than RUC, increase probabilities using special upscaling equation n If area of interest is smaller than area of model grid, recalculate strikes and use downscaling equation to get probabilities n Integrity - Service - Excellence 24

Backup Slides n Downscaling equation details n Inputs: n Strikes (S) n n n horizontal resolution of coarse area in km (C) n horizontal resolution of fine area in km (F) Equation: 1 -[1 -(F^2/C^2)]^(S^A) n Where A is a “fudge factor” depending on F n A=1 -0. 17*LN(F-1) A equals unity when F is 2 km, and slowly decreases toward zero as F approaches ~350 km In nature, lightning tends to be randomly distributed at 2 km (storm scale) but more clustered at higher resolutions. “A” attempts to account for this Best to use this equation from 2 to 128 km grid sizes If strikes is less than one, calculate equation using 1 strike, then multiply result times number of strikes Integrity - Service - Excellence 25

Backup Slides n Upscaling n Probability added to: n [1 -probability]*[1 -(F^2/C^2)]*downscaled probability This ensures high probabilities will only occur when the original probability was high, or the area has increased substantially with moderately high initial probabilities n No testing as to whether this is calibrated n Integrity - Service - Excellence 26

NWS Topeka forecast taken from the web on 15 Aug: Friday, August 17 at 7 pm Temperature: 89°F Thunder: <10% Backup Slides Integrity - Service - Excellence 27

Backup Slides Integrity - Service - Excellence 28

Backup Slides Integrity - Service - Excellence 29