Improved Initialization and Prediction of Clouds in Numerical - - PowerPoint PPT Presentation

improved initialization and prediction of clouds
SMART_READER_LITE
LIVE PREVIEW

Improved Initialization and Prediction of Clouds in Numerical - - PowerPoint PPT Presentation

Sixth Symposium on Data Assimilation. Washington, DC. Oct. 7-11, 2013 Improved Initialization and Prediction of Clouds in Numerical Weather Prediction Tom Aulign National Center for Atmospheric Research Acknowledgments: Gael Descombes,


slide-1
SLIDE 1

Sixth Symposium on Data Assimilation. Washington, DC. Oct. 7-11, 2013

Improved Initialization and Prediction of Clouds

in Numerical Weather Prediction

Tom Auligné National Center for Atmospheric Research

Acknowledgments: Gael Descombes, Francois Vandenberghe, Dongmei Xu, Thomas Nehrkorn, Brian Woods, Yann Michel, Greg Thompson.

slide-2
SLIDE 2

Initialization of clouds in NWP models

Model Observation Non-linear model & radiative transfer Complex balance Underdetermined problem Significant model errors …

Challenging NWP initialization

slide-3
SLIDE 3

Our approach to initializing clouds

  • WRF (Weather Research and Forecasting)

regional, non-hydrostatic model

  • All-sky satellite radiances
  • Expansion of analysis control variable
  • Total water + linearized physics
  • Microphysical parameters
  • Hybrid data assimilation (variational/ensemble)
slide-4
SLIDE 4

NEW ALGORITHM: update ensemble perturbations within variational analysis

3D/4D-Hybrid: ensemble covariance included via state augmentation

(Lorenc 2003, Wang et al. 2008, Fairbairn et al., 2012)

Ensemble/Variational Integrated Localized (EVIL)

dx = bcdxc + bedxe

dxe = (Pf Ca)1/2va =

dxc = B1/2v

with

Climatology Ensemble Localization

X a = X f I + zk qk

  • 1

2 -1

æ è ç ö ø ÷ zk

T k=1 K

å

æ è ç ö ø ÷

(Gratton et al., 2011)

J(v,va) = Jo + 1 2 vTv + 1 2 va

Tva[

where zk,qk

( ) are Ritz pairs from Lanczos algorithm

slide-5
SLIDE 5

Control Variable Transform

  • Multivariate covariances for qc, qr, qi, qsn
  • Binning using dynamical cloud mask
  • Vertical and Horizontal autocorrelations (Recursive Filters)
  • 3D Variance

Poster Descombes (A-p06)

slide-6
SLIDE 6

Displacement Pre-Processing

Forecast Calibration & Alignment (Grassotti et al. 1999) OSSE: Hurricane Katrina Synthetic observations

(Total Column Precipitable Water)

Balanced displacement

(Nehrkorn et al. 2013)

Poster Nehrkorn (H-p22)

slide-7
SLIDE 7
  • IR and MW radiance: AIRS, IASI, CrIS, MODIS, GOES,

AMSU-A/B, MHS, SSMI/S

  • VarBC: Variational Bias Correction
  • Revisited QC and thinning: to conserve cloudy information
  • Huber Norm: robust definition of observation error
  • Land Surface: T

skin, εs introduced as sink variable

  • Field of View: advanced interpolation scheme
  • CRTM Jacobians: rescaled base state

(floor and ceiling values for cloud parameters)

  • Middle Loop: Multiple re-linearizations of obs. operator

Processing All-Sky Satellite data

Normalized departures

slide-8
SLIDE 8

First Guess Second Guess Third Guess Observation

Update of qcloud, qice in WRF

Observations

AIRS Window Channel #787

Guess 1 Guess 2 Guess 3

slide-9
SLIDE 9

CONUS 15km, 20012/06/03 (12UTC) WRF-ARW model, Thompson microphysics First Guess = Mean of 50-member ensemble from EnKF experiment (courtesy Romine) No displacement pre-processing

CTRL = no DA

  • 3DVAR

Multivariate B matrix (5 middle-loops)

  • EVIL

3D-Hybrid-EnVar (5 middle-loops)

Experimental Demonstration

slide-10
SLIDE 10

qcloud

(level 10)

3DVAR EVIL qice

(level 20)

slide-11
SLIDE 11

3DVAR EVIL qcloud

(level 10)

qice

(level 20)

slide-12
SLIDE 12

0.12° 0.25° 0.5° 1° 2° 3°

Multi-scale verification

slide-13
SLIDE 13

Multi-scale verification Analysis Forecast

EVIL 3DVAR CTRL EVIL 3DVAR CTRL

GOES-Imager (channel 5)

slide-14
SLIDE 14
  • Expansion of analysis vector for clouds
  • Multivariate, flow-dependent background errors
  • Displacement pre-processing
  • Updated processing of all-sky satellite observations
  • Sustained impact in short-term forecast
  • More work required…

Conclusion