Assimilating lightning data into numerical forecast models use of the ensemble kalman filter
1 / 44

Assimilating Lightning Data Into Numerical Forecast Models: Use of the Ensemble Kalman Filter - PowerPoint PPT Presentation

  • Uploaded on

Assimilating Lightning Data Into Numerical Forecast Models: Use of the Ensemble Kalman Filter. Greg Hakim, Cliff Mass, Phil Regulski, Ryan Torn Department of Atmospheric Sciences University of Washington. Vaisala ILMC Meeting Tucson, April 24-25, 2008.

I am the owner, or an agent authorized to act on behalf of the owner, of the copyrighted work described.
Download Presentation

PowerPoint Slideshow about 'Assimilating Lightning Data Into Numerical Forecast Models: Use of the Ensemble Kalman Filter' - karah

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.While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server.

- - - - - - - - - - - - - - - - - - - - - - - - - - E N D - - - - - - - - - - - - - - - - - - - - - - - - - -
Presentation Transcript
Assimilating lightning data into numerical forecast models use of the ensemble kalman filter
Assimilating Lightning Data Into Numerical Forecast Models: Use of the Ensemble Kalman Filter

Greg Hakim, Cliff Mass, Phil Regulski, Ryan Torn

Department of Atmospheric Sciences

University of Washington

Vaisala ILMC Meeting

Tucson, April 24-25, 2008

Use of lightning data in numerical weather prediction nwp previous studies
Use of Lightning Data in Numerical Weather Prediction (NWP): Previous Studies

  • Earlier studies have generally used fairly primitive assimilation approaches or were completed during earlier periods without the massive amounts of observations that are now available from satellite and aircraft.

  • Several of these studies have noted substantial forecast improvements using lightning data.

Poorly forecast 1993 superstorm
Poorly Forecast 1993 Superstorm Previous Studies

Lots of lightning during its early developmental stages over the Gulf

Alexander study 1993 superstorm
Alexander Study: 1993 Superstorm Previous Studies

  • A relationship between lightning flash rate and convective precipitation was used to alter the latent heating rate in the MM5 during a spin-up period.

  • Precipitation based on satellite microwave information was also used.

  • The model was then run in forecast mode, improving predictions when satellite and lightning data were used.

Assimilating lightning data into numerical forecast models use of the ensemble kalman filter

Assimilation of Pacific Lightning Data into a Mesoscale NWP Model

Antti Pessi, Steven Businger, and Tiziana Cherubini

University of Hawaii

K. Cummins, N. Demetriades, and T. Turner

Vaisala Thunderstorm Group Inc. Tucson, AZ

Assimilating lightning data into numerical forecast models use of the ensemble kalman filter

Conversion of Lightning Rate to Moisture Profile Model

  • Determined the relationship between convective rainfall and lightning rate.

  • Determined the relationship of rainfall with the moisture profile using MM5 data.

  • Thus, Lightning rate => rainfall rate => moisture profile

  • Nudged moisture in MM5 model towards the moisture profile

Assimilating lightning data into numerical forecast models use of the ensemble kalman filter

L Model983


Reducing Forecast Error over the Eastern Pacific


Assimilation of lightning data results in a significantly improved forecast of storm central pressure (December 18-19, 2002).

The big questions
The Big Questions Model

What is the potential impact of lightning, particularly over the oceans, now that there are massive amounts of satellite information from cloud and moisture track winds, as well as increasing number of satellite vertical soundings and scatterometer winds. Plus, increasing aircraft observations.

The big questions1
The Big Questions Model

  • Does lightning data provide information content and potential forecast improvements that are not available from conventional and satellite assets?

  • What is the impact of new data assimilation approaches that allow better use of conventional and non-conventional data? Will this allow a more effective use of lightning data? Or will it make lightning data redundant with other data sources?

The university of washington lightning assimilation project
The University of Washington Lightning Assimilation Project Model

  • The UW has has been working for several years, both in research and operational modes, with a new type of data assimilation that has a number of potential advantages over more traditional types of data assimilation, such as nudging and 3D-VAR.

  • Known as the Ensemble Kalman Filter (EnKF), this approach is essentially probabilistic and makes use of the modeling system as a central component of the data assimilation process.

Enkf primer
EnKF Primer Model

  • Modern data assimilation systems combine the background (or model first guess) fields and observations to produce an optimal analysis.

  • A key element of such data assimilation systems is the background error covariance matrix, which spreads errors in the background fields both spatially and among other parameters.

  • Current data assimilation approaches, such a 3D-Var spread the errors using simplified structures and functions that are not necessarily realistic.

Covariance structures in 3dvar
Covariance structures in 3dvar Model



Data assimilation
Data Assimilation Model

  • Data assimilation should be probabilistic, providing uncertainty information regarding the analyses and the forecasts derived from them.

  • Data assimilation should also spread information among parameters, say using a precipitation (or lightning) observation to update other parameters such as wind or temperature.

  • Ensemble-based data assimilation and particularly the Ensemble Kalman filter offers a way to do this.

  • Makes use of an ensemble of forecasts to produce state-dependent error covariance structures, uncertainty information for analyses and forecasts, and allows the spread of information among parameters.

State dependent covariance matrices
State-dependent Covariance Matrices Model







Summary of ensemble kalman filter enkf algorithm
Summary of ModelEnsemble Kalman Filter (EnKF) Algorithm

  • Begin with a large ensemble of forecasts.

  • Ensemble forecast provides background error covariance statistics (B) for new analyses. (How to spread errors)

  • Ensemble analyses with new observations using these covariance structures. Many analyses and uncertainty information.

  • Make short forecasts of all ensemble members until the next observation time.

Mesoscale example cov v q rain
Mesoscale Example: cov(|V|, q Modelrain)

A nice example: the Puget Sound Convergence Zone

Experiment design

Eastern Pacific Ocean Model

Relatively low observation density; location of important storm tracks; errors propagate downstream to mainland United States

Other studies with similar domain

Pessi/Businger previously studied domain for lightning assimilation

Experiment Design

Experiment design1

Observations Model

Control case


Surface stations (ASOS, ship, buoy)


Cloud drift winds (no sat. radiances)

Experimental cases

Control observations


Experiment Design

Experiment design2

The WRF Model. Model

WRF 2.1.2 (Jan 27, 2006)

100 by 86 grid

45-km horizontal resolution

33 vertical levels

270 second timestep

Shortwave: Dudhia

Longwave: Rrtm

Surface: Noah land-sfc

PBL: MYJ TKE scheme

Cumulus: Kain-Fritsch (new Eta)

Experiment Design

Experiment design3

EnKF Setup Model

90 ensemble members

6-hr Analyses

24-hr Forecasts (starting every 12 hours)

8 assimilations period for “spin-up” before lightning assimilations

Square root filter (Whitaker and Hamill, 2002)

Horizontal localization – Gaspari and Cohn 5th order piecewise

Fixed covariance perturbations to lateral boundaries

Zhang covariance inflation method

Localization radius – 2000 km

Experiment Design

Experiment observations example acars observations spatial distribution
Experiment observations example ModelACARS observations spatial distribution

Experiment observations example cloud track wind observations spatial distribution
Experiment observations example ModelCloud track wind observations spatial distribution

Experiment observations example radiosonde surface station and buoy observations
Experiment observations example ModelRadiosonde, surface station and buoy observations

  • Radiosonde Obs

  • Surface Stations

  • Buoys

Experiment observations example lightning observations
Experiment observations example ModelLightning Observations

Test cases
Test Cases Model

  • Test Case #1

    • December 16-21, 2002 (already considered by Businger and Pessi)

  • Test Case #2

    • October 4-8, 2004

  • Test Case #3

    • November 8-12, 2006

Lightning assimilation techniques
Lightning Assimilation Techniques Model

Converted the density of lighting observations into convective rainfall using the Pessi/Businger Lightning rate/Convective rainfall rate relationship

Lightning assimilation
Lightning Assimilation Model

  • Then the convective rainfall was assimilated using the ensemble-based covariances to influence a wide variety of parameters.

  • We tried thinning and not thinning the lightning observations.

  • We tried assimilating the lightning over various periods.

  • We verified both the quality of the analyses and forecasts.

Lightning assimilation techniques1
Lightning Assimilation Techniques Model

  • Non-thinned Lightning Experiment

    • Lightning strike observations are converted into 30 minute lightning density rate from nearby LTNG observations.

    • Lightning rate converted into “observation” of convective rainfall rate using Pessi/Businger convective rain rate/lightning rate relationship

    • Convective rainfall (mm) is assimilated into WRF-EnKF

Lightning assimilation techniques2
Lightning Assimilation Techniques Model

  • Thinned Lightning Experiment

    • Same as the previous experiment except that any lightning strikes used in the density calculation are no longer allowed to be an assimilation point, resulting in a thinning out of the lightning “observations” (although strikes will be used to calculate nearby densities)

  • One hour and six hour lightning assimilation experiments. In all cases we calculate the lightning-based convective rainfall using lighting plus or minus one hour from the nominal observing time.

    • One hour-compared that rate to the one hour convective rainfall in model.

    • Six hour-scaled it to 6 hr and compared to six hour precipitation in model.

Results from the latest experiment
Results from the latest experiment Model

  • Thinned lightning

  • 1-hr precipitation assimilation (which should be more realistic)

  • Realistic error variance for lightning precipitation retrieval (5 mm)

  • Comparisons to GFS analysis

    • Although generally the best analysis provided by NCEP, the GFS analysis is certainly imperfect, especially for fine scale features.

Future work
Future Work compared to the no-lightning control?

  • Evaluation of other approaches to connecting lightning with meteorological variables:

    • One approach would be to connect lightning with graupel, or with some combination of strong vertical motion and cloud ice. Perhaps more general.

  • Improvements in the WRF EnKF, including experiments with varying EnKF settings (localization ratios, etc).

  • Increasing frequency to 3hr.

  • Weight lining with the lightning detection efficiencies.