1 / 30

3DVAR Storm-scale assimilation in real-time

This study focuses on the use of real-time 3DVAR storm-scale assimilation to improve tornado forecasting accuracy. By assimilating data from multiple sensors, including radar, the study aims to generate probabilistic tornado guidance and track potentially tornadic supercells. The system uses a high-resolution 3D analysis and ensemble forecasts to provide probabilities of hazardous weather. The results show improved forecast skill, with storm tracks and trends predicted up to an hour in advance.

sbrunelle
Download Presentation

3DVAR Storm-scale assimilation in real-time

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. 3DVAR Storm-scale assimilation in real-time Jidong Gao, Kristin Kuhlman, Travis Smith, David Stensrud

  2. Convective-scale Warn-on-Forecast Vision Radar and Initial Forecast at 2100 CST Radar at 2130 CST: Accurate Forecast Probabilistic tornado guidance: Forecast looks on track, storm circulation (hook echo) is tracking along centerline of highest tornadic probabilities An ensemble of storm-scale NWP models predict the path of a potentially tornadic supercell during the next 1 hour. The ensemble is used to create probabilistic tornado guidance. Most Likely Tornado Path Most Likely Tornado Path Developing thunderstorm 30% 30% 50% 50% 70% 70% T=2200 CST T=2200 CST T=2150 T=2150 T=2140 T=2140 T=2130 T=2130 T=2120 CST T=2120 CST Stensrud et al. 2009 (October BAMS)

  3. Probabilistic convective-scale analysis and forecast system • High-resolution 3-d analysis produced every 5 min • 1-h ensemble forecasts produced every 5 min • Provides probabilities of hazardous weather • Assimilation of all available observations

  4. Multiple sensors + Denser Observations! Radar #1 Radar #2 Radar #3 4

  5. Near-surface reflectivity Reflectivity @ -20 C Total Lightning Density Max Expected Size of Hail (~6.5 km AGL) Example: Multi-sensor data fields • Show physical relationships between data fields from multiple sensors • Storm tracks and trends can be generated at any spatial scale, for any data fields • Future state predicted through extrapolation shows skill out to about an hour 5

  6. A better alternative? Radar #1 Radar #1 Radar #2 Radar #2 Radar #3 Radar #3 Courtesy Lou Wicker

  7. Single analysis 3D Wind (U/V/W) Solve vorticity equation Not limited by edge of echo Pres/Temp/RH that can also be used to initialize forecast models for Warn-on-Forecast

  8. Radar Data Assimilation Tornado Tracks Dowell and Wicker (2009)

  9. Real-time assimilation / Hazardous Weather Testbed (Spring 2012)

  10. Getting started: R & D goals for Spring 2012 • To create realtime weather-adaptive 3DVAR analyses at high horizontal resolution and high time frequency with all operationally available radar data from the WSR-88D network. • To use the analysis product to help detect supercells and determine if these analyses can improve forecasters awareness of the hazardous weather event.

  11. 3D-Var: An Introduction • Given a set of observations and a model of some physical parameters, what does knowledge of observations tell us about the model state? • 3D-Var finds the optimal analysis that minimizes the distance to the observations weighted by the inverse of the error variance J(x) = cost function (to be minimized such that x = xa: ) Xa= optimal analysis of field of model variables Xb= background field Y0= observations H matrix transforms vectors in model space to observation space B = background error covariance R = observational error covariance

  12. 3D-Var: An Introduction • The minimum is found by numerical minimization, which is typically an iterative process • Schematic of the cost function minimization in two dimensions. • The minimum is found by performing several searches to move the control variable to areas where the cost-function is smaller, usually by looking at the local slope (the gradient) of the cost-function. Fig. from ECMWF (Bouttier & Courtier)

  13. 3D-Var: Advantages • The essence of the 3D-Var algorithm is to rewrite a least-squares problem as the minimization of a cost-function. The method was introduced in order to remove the local data selection in the Optimal Interpolation algorithm, thereby performing a global analysis of the 3-D meteorological fields, hence the name. • For 3D-Var: Observations and model variables can be non-linearly related • There is no data selection – all available data used simultaneously (no jumpiness in boundaries between regions) • 4D-Var is generalized as 3D-var for observations that are distributed in time. • 4D-Var is more complex: needs linearized perturbation forecast model and its adjoint to solve the cost function minimization problem efficiently

  14. Real-time system • 1pm to 9pm daily • 4 floating domains (user-controlled) • 1 km res • 31 vertical levels • 200x200 grid points • 5-min updates • Near-real-time (<5 min latency)

  15. Real-time system • NAM forecast is used for the background field • interpolated to the 3DVAR analysis time • Uses Reflectivity and Radial Velocity from all nearby WSR-88D radars • Cloud analysis • 3DVAR Wind analysis

  16. User-controlled domain w/ MRMS domains

  17. 3D Wind field Max updraft Simulated Reflectivity & 3D Wind Vector

  18. Vertical Vorticity

  19. Vorticity Tracks 3DVAR assimilation vorticity track 0-3 km MSL KTLX azimuthal shear track 0-3 km MSL

  20. May 10, 2010 Courtesy NWS-Norman

  21. May 10 – Stanley Draper EF4

  22. Vorticity Track / Max Updraft Track – Apr 30, 2010

  23. Hail Size (MR/MS) versus updraft intensity: May 16, 2010

  24. May 10 – Stanley Draper EF4

  25. April 14, 2011 – floating domains working together Floating Domain Vorticity Track (6-hour) Updraft Track (6-hour)

  26. May 16 Hail Storm / OKC Assimilation simulated reflectivity KTLX Reflectivity

  27. Notes • 20 prior cases (from 2010, all of which are supercells) showed: • max W range: 15 to 25 m/s • max Vorticity range: 0.012 to 0.032 s-1 • values seem to vary based on range from radar • More study is needed, here • Trends of vorticity match azimuthal shear • Updraft intensity frequently trends w/ MESH

  28. A typical session • You will use the 3DVAR analysis to help with warning decision assistance. • AWIPS / WarnGen; WDSS-II for some displays • Feedback: • real-time observation / discussion (recorded to EWP blog) • web-based surveys • daily post-mortem discussion

  29. Things to ponder this week • Look long-term (5-10-15 years into the future). • Think about how the 3DVAR storm structure and morphology compare to how you would analyze the data during typical forecast/warning operations in your office. • Consider how storm-scale ensemble prediction will affect warning decision-making and communication of warning information in 10-15 years.

  30. In Conclusion • You are getting the first look at these real-time data! • It’s the start of a > 10-year journey. • Web archive / real-time images: • http://www.nssl.noaa.gov/users/jgao/public_html/analysis/RealtimeAnalysis.htm • http://tiny.cc/3DVAR

More Related