System of Multigrid Nonlinear Least-squares Four-dimensional Variational Data Assimilation for Numerical Weather Prediction (SNAP): System Formulation and Preliminary Evaluation

A new forecasting system—the System of Multigrid Nonlinear Least-squares Four-dimensional Variational (NLS-4DVar) Data Assimilation for Numerical Weather Prediction (SNAP)—was established by building upon the multigrid NLS-4DVar data assimilation scheme, the operational Gridpoint Statistical Interpolation (GSI)-based data-processing and observation operators, and the widely used Weather Research and Forecasting numerical model. Drawing upon lessons learned from the superiority of the operational GSI analysis system, for its various observation operators and the ability to assimilate multiple-source observations, SNAP adopts GSI-based data-processing and observation operator modules to compute the observation innovations. The multigrid NLS-4DVar assimilation framework is used for the analysis, which can adequately correct errors from large to small scales and accelerate iteration solutions. The analysis variables are model state variables, rather than the control variables adopted in the conventional 4DVar system. Currently, we have achieved the assimilation of conventional observations, and we will continue to improve the assimilation of radar and satellite observations in the future. SNAP was evaluated by case evaluation experiments and one-week cycling assimilation experiments. In the case evaluation experiments, two six-hour time windows were established for assimilation experiments and precipitation forecasts were verified against hourly precipitation observations from more than 2400 national observation sites. This showed that SNAP can absorb observations and improve the initial field, thereby improving the precipitation forecast. In the one-week cycling assimilation experiments, six-hourly assimilation cycles were run in one week. SNAP produced slightly lower forecast RMSEs than the GSI 4DEnVar (Four-dimensional Ensemble Variational) as a whole and the threat scores of precipitation forecasts initialized from the analysis of SNAP were higher than those obtained from the analysis of GSI 4DEnVar.

[1]  M. Buehner,et al.  Atmospheric Data Assimilation with an Ensemble Kalman Filter: Results with Real Observations , 2005 .

[2]  Ming Hu,et al.  A GSI-Based Coupled EnSRF-En3DVar Hybrid Data Assimilation System for the Operational Rapid Refresh Model: Tests at a Reduced Resolution , 2014 .

[3]  Xiaobing Feng,et al.  A non-linear least squares enhanced POD-4DVar algorithm for data assimilation , 2015 .

[4]  G. Evensen.  Data Assimilation: The Ensemble Kalman Filter , 2006 .

[5]  Xiangjun Tian,et al.  An Efficient Local Correlation Matrix Decomposition Approach for the Localization Implementation of Ensemble‐Based Assimilation Methods , 2018 .

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

[7]  J. M. Lewis,et al.  The use of adjoint equations to solve a variational adjustment problem with advective constraints , 1985 .

[8]  Wei Li,et al.  Application of the Multigrid Method to the Two-Dimensional Doppler Radar Radial Velocity Data Assimilation , 2010 .

[9]  A. Simmons,et al.  The ECMWF operational implementation of four‐dimensional variational assimilation. I: Experimental results with simplified physics , 2007 .

[10]  Brian A. Colle,et al.  Numerical Simulations of the Extratropical Transition of Floyd (1999): Structural Evolution and Responsible Mechanisms for the Heavy Rainfall over the Northeast United States , 2003 .

[11]  M. Buehner,et al.  Intercomparison of Variational Data Assimilation and the Ensemble Kalman Filter for Global Deterministic NWP. Part I: Description and Single-Observation Experiments , 2010 .

[12]  C. Bishop,et al.  Cloud-Resolving Hurricane Initialization and Prediction through Assimilation of Doppler Radar Observations with an Ensemble Kalman Filter , 2009 .

[13]  J. McWilliams,et al.  Impacts of distinct observations during the 2009 Prince William Sound field experiment: A data assimilation study , 2012 .

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

[15]  Li Wei,et al.  Fitting model fields to observations by using singular value decomposition: An ensemble‐based 4DVar approach , 2007 .

[16]  Liang Xu,et al.  Development of NAVDAS-AR: non-linear formulation and outer loop tests , 2006 .

[17]  Yuanfu Xie,et al.  Nonlinear Least Squares En4DVar to 4DEnVar Methods for Data Assimilation: Formulation, Analysis and Preliminary Evaluation , 2018 .

[18]  G. Grell,et al.  A North American Hourly Assimilation and Model Forecast Cycle: The Rapid Refresh , 2016 .

[19]  William C. Skamarock,et al.  A time-split nonhydrostatic atmospheric model for weather research and forecasting applications , 2008, J. Comput. Phys..

[20]  Andrew C. Lorenc,et al.  The potential of the ensemble Kalman filter for NWP—a comparison with 4D‐Var , 2003 .

[21]  Yuanfu Xie.  A sequential variational analysis approach for mesoscale data assimilation , 2005 .

[22]  J. Dudhia.  Numerical Study of Convection Observed during the Winter Monsoon Experiment Using a Mesoscale Two-Dimensional Model , 1989 .

[23]  P. L. Houtekamer,et al.  Ensemble Kalman filtering , 2005 .

[24]  H. D. Orville,et al.  Bulk Parameterization of the Snow Field in a Cloud Model , 1983 .

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

[26]  M. Xue,et al.  A Prototype Regional GSI-based EnKF-Variational Hybrid Data Assimilation System for the Rapid Refresh Forecasting System: Dual-Resolution Implementation and Testing Results , 2018, Advances in Atmospheric Sciences.

[27]  William L. Briggs,et al.  A multigrid tutorial, Second Edition , 2000 .

[28]  Barry E. Schwartz,et al.  An Hourly Assimilation–Forecast Cycle: The RUC , 2004 .

[29]  Xiangjun Tian,et al.  Implementations of a square-root ensemble analysis and a hybrid localisation into the POD-based ensemble 4DVar , 2012 .

[30]  M. Buehner,et al.  Intercomparison of Variational Data Assimilation and the Ensemble Kalman Filter for Global Deterministic NWP. Part II: One-Month Experiments with Real Observations , 2010 .

[31]  J. Dudhia,et al.  Coupling an Advanced Land Surface–Hydrology Model with the Penn State–NCAR MM5 Modeling System. Part I: Model Implementation and Sensitivity , 2001 .

[32]  F. Karagiannis.  AN ECONOMICAL APPROACH , 2000 .

[33]  William L. Briggs,et al.  A multigrid tutorial , 1987 .

[34]  A. Lorenc,et al.  Operational implementation of a hybrid ensemble/4D‐Var global data assimilation system at the Met Office , 2013 .

[35]  J. Whitaker,et al.  Distance-dependent filtering of background error covariance estimates in an ensemble Kalman filter , 2001 .

[36]  Y. Kuo,et al.  An economical approach to four-dimensional variational data assimilation , 2010 .

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

[38]  Ming Hu,et al.  A Regional GSI-Based Ensemble Kalman Filter Data Assimilation System for the Rapid Refresh Configuration: Testing at Reduced Resolution , 2013 .

[39]  P. Courtier,et al.  A strategy for operational implementation of 4D‐Var, using an incremental approach , 1994 .

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

[41]  R. Purser,et al.  Three-Dimensional Variational Analysis with Spatially Inhomogeneous Covariances , 2002 .

[42]  E. Mlawer,et al.  Radiative transfer for inhomogeneous atmospheres: RRTM, a validated correlated-k model for the longwave , 1997 .

[43]  Monique Tanguay,et al.  Extension of 3DVAR to 4DVAR: Implementation of 4DVAR at the Meteorological Service of Canada , 2007 .

[44]  Andrew C. Lorenc,et al.  Modelling of error covariances by 4D‐Var data assimilation , 2003 .

[45]  Kayo Ide,et al.  A three‐dimensional variational data assimilation scheme for the Regional Ocean Modeling System: Implementation and basic experiments , 2008 .

[46]  Craig H. Bishop,et al.  Comparison of Hybrid Ensemble/4DVar and 4DVar within the NAVDAS-AR Data Assimilation Framework , 2013 .

[47]  J. Derber,et al.  Introduction of the GSI into the NCEP Global Data Assimilation System , 2009 .

[48]  Peter V. Hobbs,et al.  The Mesoscale and Microscale Structure and Organization of Clouds and Precipitation in Midlatitude Cyclones. XII: A Diagnostic Modeling Study of Precipitation Development in Narrow Cold-Frontal Rainbands , 1984 .

[49]  Shu‐Hua Chen,et al.  A One-dimensional Time Dependent Cloud Model , 2002 .

[50]  Xiangjun Tian,et al.  An ensemble-based explicit four-dimensional variational assimilation method , 2008 .

[51]  Hongqin Zhang,et al.  A Big Data‐Driven Nonlinear Least Squares Four‐Dimensional Variational Data Assimilation Method: Theoretical Formulation and Conceptual Evaluation , 2019, Earth and Space Science.

[52]  Fuqing Zhang,et al.  E4DVar: Coupling an Ensemble Kalman Filter with Four-Dimensional Variational Data Assimilation in a Limited-Area Weather Prediction Model , 2012 .

[53]  Qin Xu.  Generalized Adjoint for Physical Processes with Parameterized Discontinuities. Part I: Basic Issues and Heuristic Examples , 1996 .

[54]  T. Hamill,et al.  A Hybrid Ensemble Kalman Filter-3D Variational Analysis Scheme , 2000 .

[55]  Paul E. Bieringer,et al.  A Space-Time Multiscale Analysis System: A Sequential Variational Analysis Approach , 2011 .

[56]  John E. Dennis,et al.  Numerical methods for unconstrained optimization and nonlinear equations , 1983, Prentice Hall series in computational mathematics.

[57]  Hongqin Zhang,et al.  A Multigrid Nonlinear Least Squares Four‐Dimensional Variational Data Assimilation Scheme With the Advanced Research Weather Research and Forecasting Model , 2018 .

[58]  J. Dudhia,et al.  A New Vertical Diffusion Package with an Explicit Treatment of Entrainment Processes , 2006 .