Home About HyMeX
Motivations
Science questions
Observation strategy
Modelling strategy
Target areas
Key documents
Organisation
International coordination
Working groups
Task teams
National contributions
Endorsements
Resources
Database
Data policy
Publications
Education and summer schools
Drifting balloons (BAMED)
SOP web page
Google maps data visualisation
Workshops Projects
ASICS-MED
MOBICLIMEX
MUSIC
IODA-MED
REMEMBER
FLOODSCALE
EXAEDRE
Offers Links Contacts
Science & Task teams
Science teams
Task teams
Implementation plan
Coordination
International Scientific Steering Committee (ISSC)
Executive Committee for Implementation and Science Coordination (EC-ISC)
Executive Committee - France (EC-Fr)
HyMeX France
HyMeX Italy
HyMeX Spain
Archive
by Legrand, R., Michel, Y. and Montmerle, T.
Abstract:
In numerical weather prediction, the problem of estimating initial conditions with a variational approach is usually based on a Bayesian framework associated with a Gaussianity assumption of the probability density functions of both observations and background errors. In practice, Gaussianity of errors is tied to linearity, in the sense that a nonlinear model will yield non-Gaussian probability density functions. In this context, standard methods relying on Gaussian assumption may perform poorly. This study aims to describe some aspects of non-Gaussianity of forecast and analysis errors in a convective-scale model using a Monte Carlo approach based on an ensemble of data assimilations. For this purpose, an ensemble of 90 members of cycled perturbed assimilations has been run over a highly precipitating case of interest. Non-Gaussianity is measured using the K2 statistics from the D'Agostino test, which is related to the sum of the squares of univariate skewness and kurtosis. Results confirm that specific humidity is the least Gaussian variable according to that measure and also that non-Gaussianity is generally more pronounced in the boundary layer and in cloudy areas. The dynamical control variables used in our data assimilation, namely vorticity and divergence, also show distinct non-Gaussian behaviour. It is shown that while non-Gaussianity increases with forecast lead time, it is efficiently reduced by the data assimilation step especially in areas well covered by observations. Our findings may have implication for the choice of the control variables.
Reference:
Legrand, R., Michel, Y. and Montmerle, T., 2016: Diagnosing non-Gaussianity of forecast and analysis errors in a convective scale modelNonlinear Processes in Geophysics, 23, 1-12.
Bibtex Entry:
@Article{Legrand2016,
  Title                    = {Diagnosing non-Gaussianity of forecast and analysis errors in a convective scale model},
  Author                   = {Legrand, R. and Michel, Y. and Montmerle, T.},
  Journal                  = {Nonlinear Processes in Geophysics},
  Year                     = {2016},

  Month                    = {January},
  Number                   = {1},
  Pages                    = {1-12},
  Volume                   = {23},

  Abstract                 = {In numerical weather prediction, the problem of estimating initial conditions with a variational approach is usually based on a Bayesian framework associated with a Gaussianity assumption of the probability density functions of both observations and background errors. In practice, Gaussianity of errors is tied to linearity, in the sense that a nonlinear model will yield non-Gaussian probability density functions. In this context, standard methods relying on Gaussian assumption may perform poorly.
This study aims to describe some aspects of non-Gaussianity of forecast and analysis errors in a convective-scale model using a Monte Carlo approach based on an ensemble of data assimilations. For this purpose, an ensemble of 90 members of cycled perturbed assimilations has been run over a highly precipitating case of interest. Non-Gaussianity is measured using the K2 statistics from the D'Agostino test, which is related to the sum of the squares of univariate skewness and kurtosis.
Results confirm that specific humidity is the least Gaussian variable according to that measure and also that non-Gaussianity is generally more pronounced in the boundary layer and in cloudy areas. The dynamical control variables used in our data assimilation, namely vorticity and divergence, also show distinct non-Gaussian behaviour. It is shown that while non-Gaussianity increases with forecast lead time, it is efficiently reduced by the data assimilation step especially in areas well covered by observations. Our findings may have implication for the choice of the control variables.},
  Copublication            = {3: 3 Fr},
  Doi                      = {10.5194/npg-23-1-2016},
  Owner                    = {hymexw},
  Timestamp                = {2016.03.17},
  Url                      = {http://www.nonlin-processes-geophys.net/23/1/2016/npg-23-1-2016.html}
}