1 / 37

Probabilistic QPF

Probabilistic QPF. Steve Amburn, SOO WFO Tulsa, Oklahoma. Preview. Product review Graphic and text Forecaster involvement Examples of rainfall distributions Theory and Definitions Comparisons To previous work To real events. How do we distribute PQPF information? .

bryson
Download Presentation

Probabilistic QPF

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript


  1. Probabilistic QPF Steve Amburn, SOO WFO Tulsa, Oklahoma

  2. Preview • Product review • Graphic and text • Forecaster involvement • Examples of rainfall distributions • Theory and Definitions • Comparisons • To previous work • To real events

  3. How do we distribute PQPF information?

  4. Bar Chart – Osage County Average

  5. Distribute of text information Probability to exceed 0.10” = 67% Probability to exceed 0.50” = 46% Probability to exceed 1.00” = 30% Probability to exceed 2.00” = 12% 6-hr PoP = 70% 6-hr QPF = 0.82

  6. Premise • Rainfall events are typically characterized by a distribution or variety of rainfall amounts. • Every distribution has a mean. • Data indicates most rainfall events have exponential or gamma distributions. • Forecast the mean, and you forecast the distribution. • That forecast distribution lets us calculate exceedance probabilities (POEs), i.e., the probabilistic QPF.

  7. Characteristics of POE Method for Probabilistic QPF (exponential distributions) • As mean increases POE increases

  8. Forecaster Involvement • No additional workload • No extra grids to edit • No extra time to create the products • Forecasters simply issue their QPF, given the more specific definition. • New QPF definition ~ Old QPF definition

  9. Gamma Distribution • “There are a variety of continuous distributions that are bounded on the left by zero and positively skewed. One commonly used choice, used especially often for representing precipitation data, is the gamma distribution. ...” D.S. Wilks (1995), Statistical Methods in the Atmospheric Sciences, Academic Press, pp. 467.

  10. Gamma Distributions exponential distribution (special case of gamma) blue line => α =1red line => α=2black line => α=3

  11. Individual Spring Events

  12. Individual Summer Events

  13. Individual Autumn Events

  14. Individual Winter Events

  15. Gamma Distribution? YesUsually, the Exponential Distribution • Over time at a point (TUL, FSM...) • For individual events (nearly 700) • Virtually all were exponential distributions • A few were other forms of gamma distribution • Let’s look at the math...

  16. Gamma Function The gamma function is defined by: Г(α) = ∫ x (α-1) e-x dx ,for α > 0 , for 0 → ∞.

  17. Gamma Distributions exponential distribution (special case of gamma) blue line => α =1red line => α =2black line => α =3

  18. Gamma Function for α = 1 • Г(α) = ∫ x (α-1) e-x dx ,for α > 0 , for 0 → ∞. • So, for α = 1, • Г(1) =∫X(1-1)e-x dx • =∫e-x dx = - e-∞ - (-e-0) = 0 + 1 • Г(1) = 1. ∞ 0

  19. Gamma PDF  Exponential PDF Now, for α = 1, gamma distribution PDF simplifies. f(x) = { 1 / [ βα Γ(α) ] } • xα-1 e-x/β, (gamma density function) So, substituting α = 1 yields: f(x) ={ 1 / [β1 Γ(1) ] } • x1-1 e-x/β or, f(x) = (1/β) • e-x/β, (exponential density function)

  20. POE for theExponential Distribution So, to find the probability to exceed a value “x”, we integrate from x to infinity (∞). f(x) = POE(x) = ∫ (1/β) • e-x/β , from x → ∞. = -e-∞/β - (-e-x/β ) = 0 + e-x/β so, POE(x) = e-x/β,where β = mean

  21. Unconditional POE • Simply multiply conditional POE by PoP, • uPOE(x) = PoP x (e-x/μ) • , where μ is the conditional QPF. • Consider, the NWS PoP = uPOE (0.005)

  22. Individual event examples • Data from ABRFC • 4578 HRAP grid boxes in TSA CWFA. • Rainfall distributions created • Means calculated • POEs calculated from observed data • Actual means used to calculate POEs from PDFs (formula)

  23. Actual POE = computed from the observed data. Perfect POE = computed using the exponential equation and the observed mean

  24. Actual POE = computed from the observed data Perfect POE = computed using the exponential equation and the observed mean Climo POE = exp equation with climatological mean

  25. Actual POE = computed from the observed data Exponential Eq = computed using the exponential equation and the observed mean

  26. Actual POE computed from frequency data shown below. Using integrated exponential distribution (gamma distribution where alpha = 1) to compute POE. POE using integrated exponential distribution and climatological mean precipitation.

  27. Actual POE computed from frequency data shown below. POE using integrated gamma distribution (alpha=3) POE from exponential distribution (gamma distribution, alpha = 1). POE using integrated exponential distribution and climatological mean precipitation.

  28. Gamma Distribution Gamma function is defined by: Г(α) = x (α-1) e-x dx ,for α > 0 , for 0 → ∞. Gamma density function is given by: 1 ] _________ [ f(x) = • (xα-1 • e–x/β ) β • Γ(α)

  29. Gamma Distribution, A = 3 1 _________ [ ] • (xα-1 • e–x/β ) f(x) = β • Γ(α) Then, integrate, for α =3, from x to infinity to obtain the probability of exceedance for x. POE(x) =(0.5)•(e –x/β)•(x2/β + 2x/β + 2) ,where β = µ/α = (qpf / 3)

  30. Gamma, where a=3

  31. Gamma, where a=3

  32. Gamma, where a=3

  33. Gamma, where a=3

  34. Error in Exceedance Calculations (419 events) • Calculated at: 0.10, 0.25, 0.50, 1.0, 1.5” Areal Coverage

  35. Error in Exceedance Calculations (279 events) • Calculated at: 0.10, 0.25, 0.50, 1.0, 2.0, 3.0, 4.0 inches Areal Coverage

  36. Summary • Rainfall events have distributions • Typically exponential • Gamma for coverage > 90% • Each distribution has a mean. • Mean can be spatial for a single event • Mean can be at a point over a period of time • Forecast the mean (QPF) and you have effectively forecast the distribution. • So, QPF lets us calculate probabilities of exceedance, or PQPF.

  37. Questions? Steve Amburn, SOO WFO Tulsa, Oklahoma

More Related