Adaptive covariance relaxation methods for ensemble data assimilation: experiments in the real atmosphere

Covariance inflation plays an important role in the ensemble Kalman filter because the ensemble-based error variance is usually underestimated due to various factors such as the limited ensemble size and model imperfections. Manual tuning of the inflation parameters by trial and error is computationally expensive; therefore, several studies have proposed approaches to adaptive estimation of the inflation parameters. Among others, this study focuses on the covariance relaxation method which realizes spatially dependent inflation with a spatially homogeneous relaxation parameter. This study performs a series of experiments with the non-hydrostatic icosahedral atmospheric model (NICAM) and the local ensemble transform Kalman filter (LETKF) assimilating the real-world conventional observations and satellite radiances. Two adaptive covariance relaxation methods are implemented: relaxation to prior spread based on Ying and Zhang (adaptive-RTPS), and relaxation to prior perturbation (adaptive-RTPP). Both adaptive-RTPS and adaptive-RTPP generally improve the analysis compared to a baseline control experiment with an adaptive multiplicative inflation method. However, the adaptive-RTPS and adaptive-RTPP methods lead to an over-dispersive (under-dispersive) ensemble in the sparsely (densely) observed regions compared with the adaptive multiplicative inflation method. We find that the adaptive-RTPS and adaptive-RTPP methods are robust to a sudden change in the observing networks and observation error settings.

[1]  Takemasa Miyoshi,et al.  Ensemble Kalman Filter and 4D-Var Intercomparison with the Japanese Operational Global Analysis and Prediction System , 2010 .

[2]  Paul Poli,et al.  Diagnosis of observation, background and analysis‐error statistics in observation space , 2005 .

[3]  Masaki Satoh,et al.  Nonhydrostatic icosahedral atmospheric model (NICAM) for global cloud resolving simulations , 2008, J. Comput. Phys..

[4]  Jeffrey L. Anderson.  Spatially and temporally varying adaptive covariance inflation for ensemble filters , 2009 .

[5]  Takemasa Miyoshi,et al.  Assimilation of TRMM Multisatellite Precipitation Analysis with a Low-Resolution NCEP Global Forecast System , 2016 .

[6]  G. Evensen.  Sequential data assimilation with a nonlinear quasi‐geostrophic model using Monte Carlo methods to forecast error statistics , 1994 .

[7]  Lars Nerger,et al.  On Serial Observation Processing in Localized Ensemble Kalman Filters , 2015 .

[8]  Takemasa Miyoshi,et al.  The Non-hydrostatic Icosahedral Atmospheric Model: description and development , 2014, Progress in Earth and Planetary Science.

[9]  Istvan Szunyogh,et al.  Efficient data assimilation for spatiotemporal chaos: A local ensemble transform Kalman filter , 2005, physics/0511236.

[10]  Koji Terasaki,et al.  Local Ensemble Transform Kalman Filter Experiments with the Nonhydrostatic Icosahedral Atmospheric Model NICAM , 2015 .

[11]  T. Miyoshi.  The Gaussian Approach to Adaptive Covariance Inflation and Its Implementation with the Local Ensemble Transform Kalman Filter , 2011 .

[12]  Jeffrey L. Anderson,et al.  An adaptive covariance inflation error correction algorithm for ensemble filters , 2007 .

[13]  D. Dee.  On-line Estimation of Error Covariance Parameters for Atmospheric Data Assimilation , 1995 .

[14]  Eugenia Kalnay,et al.  Ensemble Forecasting at NMC: The Generation of Perturbations , 1993 .

[15]  P. Houtekamer,et al.  Data Assimilation Using an Ensemble Kalman Filter Technique , 1998 .

[16]  Istvan Szunyogh,et al.  Use of the breeding technique to estimate the structure of the analysis "errors of the day" , 2003 .

[17]  J. Whitaker,et al.  Ensemble Data Assimilation without Perturbed Observations , 2002 .

[18]  Jeffrey L. Anderson,et al.  A Monte Carlo Implementation of the Nonlinear Filtering Problem to Produce Ensemble Assimilations and Forecasts , 1999 .

[19]  J. Thepaut,et al.  The ERA‐Interim reanalysis: configuration and performance of the data assimilation system , 2011 .

[20]  K. Emanuel,et al.  Optimal Sites for Supplementary Weather Observations: Simulation with a Small Model , 1998 .

[21]  Jonathan Poterjoy,et al.  Comparison of Hybrid Four-Dimensional Data Assimilation Methods with and without the Tangent Linear and Adjoint Models for Predicting the Life Cycle of Hurricane Karl (2010) , 2016 .

[22]  Kumiko Takata,et al.  Development of the minimal advanced treatments of surface interaction and runoff , 2003 .

[23]  J. Whitaker,et al.  Ensemble Square Root Filters , 2003, Statistical Methods for Climate Scientists.

[24]  M. Kanamitsu,et al.  Observation system simulation experiments using water vapor isotope information , 2014 .

[25]  Jeffrey L. Anderson.  An Ensemble Adjustment Kalman Filter for Data Assimilation , 2001 .

[26]  J. Whitaker,et al.  Evaluating Methods to Account for System Errors in Ensemble Data Assimilation , 2012 .

[27]  Niels Bormann,et al.  Estimates of spatial and interchannel observation‐error characteristics for current sounder radiances for numerical weather prediction. I: Methods and application to ATOVS data , 2010 .

[28]  G. Evensen,et al.  Analysis Scheme in the Ensemble Kalman Filter , 1998 .

[29]  Matthias Zahn,et al.  A comparison of two identification and tracking methods for polar lows , 2012 .

[30]  Juanzhen Sun,et al.  Impacts of Initial Estimate and Observation Availability on Convective-Scale Data Assimilation with an Ensemble Kalman Filter , 2004 .

[31]  T. Hamill.  Interpretation of Rank Histograms for Verifying Ensemble Forecasts , 2001 .

[32]  P. Houtekamer,et al.  An Adaptive Ensemble Kalman Filter , 2000 .

[33]  T. Miyoshi,et al.  Nowcasting with Data Assimilation: A Case of Global Satellite Mapping of Precipitation , 2016 .

[34]  Geir Evensen,et al.  The Ensemble Kalman Filter: theoretical formulation and practical implementation , 2003 .

[35]  Takemasa Miyoshi,et al.  Using AIRS retrievals in the WRF-LETKF system to improve regional numerical weather prediction , 2012 .

[36]  Craig H. Bishop,et al.  Adaptive sampling with the ensemble transform Kalman filter , 2001 .

[37]  B. Hunt,et al.  A comparative study of 4D-VAR and a 4D Ensemble Kalman Filter: perfect model simulations with Lorenz-96 , 2007 .

[38]  F. Rabier,et al.  The potential of high‐density observations for numerical weather prediction: A study with simulated observations , 2003 .

[39]  E. Kalnay,et al.  Simultaneous estimation of covariance inflation and observation errors within an ensemble Kalman filter , 2009 .

[40]  Hirofumi Tomita,et al.  A new dynamical framework of nonhydrostatic global model using the icosahedral grid , 2004 .

[41]  Fuqing Zhang,et al.  An adaptive covariance relaxation method for ensemble data assimilation , 2015 .

[42]  Takemasa Miyoshi,et al.  Local Ensemble Transform Kalman Filtering with an AGCM at a T159/L48 Resolution , 2007 .

[43]  M. Rapp,et al.  Journal of Geophysical Research : Atmospheres Does Strong Tropospheric Forcing Cause Large-Amplitude Mesospheric Gravity Waves ? A DEEPWAVE Case Study , 2017 .

[44]  A. Arakawa,et al.  Interaction of a Cumulus Cloud Ensemble with the Large-Scale Environment, Part I , 1974 .