scholarly journals Evaluation of a Spatial/Spectral Covariance Localization Approach for Atmospheric Data Assimilation

2012 ◽  
Vol 140 (2) ◽  
pp. 617-636 ◽  
Author(s):  
Mark Buehner

In this study, several approaches for estimating background-error covariances from an ensemble of error realizations are examined, including a new spatial/spectral localization approach. The new approach shares aspects of both the spatial localization and wavelet-diagonal approaches. This approach also enables the use of different spatial localization functions for the covariances associated with each of a set of overlapping horizontal wavenumber bands. The use of such scale-dependent spatial localization (more severe localization for small horizontal scales) is shown to reduce the error in spatial correlation estimates. A comparison of spatial localization, spatial/spectral localization, and wavelet-diagonal approaches shows that the approach resulting in the lowest estimation error depends on the ensemble size. For a relatively large ensemble (48 members), the spatial/spectral localization approach produces the lowest error. When using a much smaller ensemble (12 members), the wavelet-diagonal approach results in the lowest error. Qualitatively, the horizontal correlation functions resulting from spatial/spectral localization appear smoother and less noisy than those from spatial localization, but preserve more of the heterogeneous and anisotropic nature of the raw sample correlations than the wavelet-diagonal approach. The new spatial/spectral localization approach is compared with spatial localization in a set of 1-month three-dimensional variational data assimilation (3D-Var) experiments using a full set of real atmospheric observations. Preliminary results show that spatial/spectral localization provides a nearly similar forecast quality, and in some regions improved forecast quality, as spatial localization while using an ensemble of half the size (48 vs 96 members).

2014 ◽  
Vol 142 (10) ◽  
pp. 3756-3780 ◽  
Author(s):  
Yujie Pan ◽  
Kefeng Zhu ◽  
Ming Xue ◽  
Xuguang Wang ◽  
Ming Hu ◽  
...  

Abstract A coupled ensemble square root filter–three-dimensional ensemble-variational hybrid (EnSRF–En3DVar) data assimilation (DA) system is developed for the operational Rapid Refresh (RAP) forecasting system. The En3DVar hybrid system employs the extended control variable method, and is built on the NCEP operational gridpoint statistical interpolation (GSI) three-dimensional variational data assimilation (3DVar) framework. It is coupled with an EnSRF system for RAP, which provides ensemble perturbations. Recursive filters (RF) are used to localize ensemble covariance in both horizontal and vertical within the En3DVar. The coupled En3DVar hybrid system is evaluated with 3-h cycles over a 9-day period with active convection. All conventional observations used by operational RAP are included. The En3DVar hybrid system is run at ⅓ of the operational RAP horizontal resolution or about 40-km grid spacing, and its performance is compared to parallel GSI 3DVar and EnSRF runs using the same datasets and resolution. Short-term forecasts initialized from the 3-hourly analyses are verified against sounding and surface observations. When using equally weighted static and ensemble background error covariances and 40 ensemble members, the En3DVar hybrid system outperforms the corresponding GSI 3DVar and EnSRF. When the recursive filter coefficients are tuned to achieve a similar height-dependent localization as in the EnSRF, the En3DVar results using pure ensemble covariance are close to EnSRF. Two-way coupling between EnSRF and En3DVar did not produce noticeable improvement over one-way coupling. Downscaled precipitation forecast skill on the 13-km RAP grid from the En3DVar hybrid is better than those from GSI 3DVar analyses.


2019 ◽  
Vol 12 (9) ◽  
pp. 4031-4051 ◽  
Author(s):  
Shizhang Wang ◽  
Zhiquan Liu

Abstract. A reflectivity forward operator and its associated tangent linear and adjoint operators (together named RadarVar) were developed for variational data assimilation (DA). RadarVar can analyze both rainwater and ice-phase species (snow and graupel) by directly assimilating radar reflectivity observations. The results of three-dimensional variational (3D-Var) DA experiments with a 3 km grid mesh setting of the Weather Research and Forecasting (WRF) model showed that RadarVar was effective at producing an analysis of reflectivity pattern and intensity similar to the observed data. Two to three outer loops with 50–100 iterations in each loop were needed to obtain a converged 3-D analysis of reflectivity, rainwater, snow, and graupel, including the melting layers with mixed-phase hydrometeors. It is shown that the deficiencies in the analysis using this operator, caused by the poor quality of the background fields and the use of the static background error covariance, can be partially resolved by using radar-retrieved hydrometeors in a preprocessing step and tuning the spatial correlation length scales of the background errors. The direct radar reflectivity assimilation using RadarVar also improved the short-term (2–5 h) precipitation forecasts compared to those of the experiment without DA.


2010 ◽  
Vol 138 (10) ◽  
pp. 3946-3966 ◽  
Author(s):  
Jean-François Caron ◽  
Luc Fillion

Abstract This study examines the modification to the balance properties of the analysis increments in a global three-dimensional variational data assimilation scheme when using flow-dependent background-error covariances derived from an operational ensemble Kalman filter instead of static homogenous and isotropic background-error covariances based on lagged forecast differences. It is shown that the degree of balance in the analysis increments is degraded when the former method is used. This change can be attributed in part to the reduced degree of rotational balance found in short-term ensemble Kalman filter perturbations as compared to lagged forecast differences based on longer-range forecasts. However, the use of a horizontal and vertical localization technique to increase the rank of the ensemble-based covariances are found to have a significant deleterious effect on the rotational balance with the largest detrimental impact coming from the vertical localization and affecting particularly the upper levels. The examination of the vertical motion part of the analysis increments revealed that the spatial covariance localization technique also produces unrealistic vertical structure of vertical motion increments with abnormally large increments near the surface. A comparison between the analysis increments from the ensemble Kalman filter and from the ensemble-based three-dimensional variational data assimilation (3D-Var) scheme showed that the balance characteristics of the analysis increments resulting from the two systems are very similar.


2015 ◽  
Vol 143 (8) ◽  
pp. 3087-3108 ◽  
Author(s):  
Aaron Johnson ◽  
Xuguang Wang ◽  
Jacob R. Carley ◽  
Louis J. Wicker ◽  
Christopher Karstens

Abstract A GSI-based data assimilation (DA) system, including three-dimensional variational assimilation (3DVar) and ensemble Kalman filter (EnKF), is extended to the multiscale assimilation of both meso- and synoptic-scale observation networks and convective-scale radar reflectivity and velocity observations. EnKF and 3DVar are systematically compared in this multiscale context to better understand the impacts of differences between the DA techniques on the analyses at multiple scales and the subsequent convective-scale precipitation forecasts. Averaged over 10 diverse cases, 8-h precipitation forecasts initialized using GSI-based EnKF are more skillful than those using GSI-based 3DVar, both with and without storm-scale radar DA. The advantage from radar DA persists for ~5 h using EnKF, but only ~1 h using 3DVar. A case study of an upscale growing MCS is also examined. The better EnKF-initialized forecast is attributed to more accurate analyses of both the mesoscale environment and the storm-scale features. The mesoscale location and structure of a warm front is more accurately analyzed using EnKF than 3DVar. Furthermore, storms in the EnKF multiscale analysis are maintained during the subsequent forecast period. However, storms in the 3DVar multiscale analysis are not maintained and generate excessive cold pools. Therefore, while the EnKF forecast with radar DA remains better than the forecast without radar DA throughout the forecast period, the 3DVar forecast quality is degraded by radar DA after the first hour. Diagnostics revealed that the inferior analysis at mesoscales and storm scales for the 3DVar is primarily attributed to the lack of flow dependence and cross-variable correlation, respectively, in the 3DVar static background error covariance.


2020 ◽  
Author(s):  
Sojin Lee ◽  
Chul Han Song ◽  
Kyung Man Han ◽  
Daven K. Henze ◽  
Kyunghwa Lee ◽  
...  

Abstract. For the purpose of improving PM prediction skills in East Asia, we estimated a new background error covariance matrix (BEC) for aerosol data assimilation using surface PM2.5 observations that accounts for the uncertainties in anthropogenic emissions. In contrast to the conventional method to estimate the BEC that uses perturbations in meteorological data, this method additionally considered the perturbations using two different emission inventories. The impacts of the new BEC were then tested for the prediction of surface PM2.5 over East Asia using Community Multi-scale Air Quality (CMAQ) initialized by three-dimensional variational method (3D-VAR). The surface PM2.5 data measured at 154 sites in South Korea and 1,535 sites in China were assimilated every six hours during the Korea-United States Air Quality Study (KORUS-AQ) campaign period (1 May–14 June 2016). Data assimilation with our new BEC showed better agreement with the surface PM2.5 observations than that with the conventional method. Our method also showed closer agreement with the observations in 24-hour PM2.5 predictions with ~ 44 % fewer negative biases than the conventional method. We conclude that increased standard deviations, together with horizontal and vertical length scales in the new BEC, tend to improve the data assimilation and short-term predictions for the surface PM2.5. This paper also suggests further research efforts devoted to estimating the BEC to improve PM2.5 predictions.


2011 ◽  
Vol 2011 ◽  
pp. 1-10 ◽  
Author(s):  
Chien-Ben Chou ◽  
Huei-Ping Huang

This work assesses the effects of assimilating atmospheric infrared sounder (AIRS) observations on typhoon prediction using the three-dimensional variational data assimilation (3DVAR) and forecasting system of the weather research and forecasting (WRF) model. Two major parameters in the data assimilation scheme, the spatial decorrelation scale and the magnitude of the covariance matrix of the background error, are varied in forecast experiments for the track of typhoon Sinlaku over the Western Pacific. The results show that within a wide parameter range, the inclusion of the AIRS observation improves the prediction. Outside this range, notably when the decorrelation scale of the background error is set to a large value, forcing the assimilation of AIRS data leads to degradation of the forecast. This illustrates how the impact of satellite data on the forecast depends on the adjustable parameters for data assimilation. The parameter-sweeping framework is potentially useful for improving operational typhoon prediction.


2013 ◽  
Vol 2013 ◽  
pp. 1-6 ◽  
Author(s):  
Hongze Leng ◽  
Junqiang Song ◽  
Fengshun Lu ◽  
Xiaoqun Cao

This study considers a new hybrid three-dimensional variational (3D-Var) and ensemble Kalman filter (EnKF) data assimilation (DA) method in a non-perfect-model framework, named space-expanded ensemble localization Kalman filter (SELKF). In this method, the localization operation is directly applied to the ensemble anomalies with a Schur Product, rather than to the full error covariance of the state in the EnKF. Meanwhile, the correction space of analysis increment is expanded to a space with larger dimension, and the rank of the forecast error covariance is significantly increased. This scheme can reduce the spurious correlations in the covariance and approximate the full-rank background error covariance well. Furthermore, a deterministic scheme is used to generate the analysis anomalies. The results show that the SELKF outperforms the perturbed EnKF given a relatively small ensemble size, especially when the length scale is relatively long or the observation error covariance is relatively small.


2014 ◽  
Vol 142 (5) ◽  
pp. 1852-1873 ◽  
Author(s):  
Eric Wattrelot ◽  
Olivier Caumont ◽  
Jean-François Mahfouf

AbstractThis paper presents results from radar reflectivity data assimilation experiments with the nonhydrostatic limited-area model Application of Research to Operations at Mesoscale (AROME) in an operational context. A one-dimensional (1D) Bayesian retrieval of relative humidity profiles followed by a three-dimensional variational data assimilation (3D-Var) technique is adopted. Several preprocessing procedures of raw reflectivity data are presented and the use of the nonrainy signal in the assimilation is widely discussed and illustrated. This two-step methodology allows the authors to build up a screening procedure that takes into account the evaluation of the results from the 1D Bayesian retrieval. In particular, the 1D retrieval is checked by comparing a pseudoanalyzed reflectivity to the observed reflectivity. Additionally, a physical consistency between the reflectivity innovations and the 1D relative humidity increments is imposed before assimilating relative humidity pseudo-observations with other observations. This allows the authors to counteract the difficulty of the current 3D-Var system to correct strong differences between model and observed clouds from the crude specification of background-error covariances. Assimilation experiments of radar reflectivity data in a preoperational configuration are first performed over a 1-month period. Positive impacts on short-term precipitation forecast scores are systematically found. The evaluation shows improvements on the analysis and also on objective conventional forecast scores, in particular for the model wind field up to 12 h. A case study for a specific precipitating system demonstrates the capacity of the method for improving significantly short-term forecasts of organized convection.


2013 ◽  
Vol 141 (8) ◽  
pp. 2721-2739 ◽  
Author(s):  
Chengsi Liu ◽  
Qingnong Xiao

Abstract A four-dimensional ensemble-based variational data assimilation (4DEnVar) algorithm proposed in Part I of the 4DEnVar series (denoted En4DVar in Part I, but here we refer to it as 4DEnVar according to WMO conference recommendation to differentiate it from En4DVar algorithm using adjoint model) uses a flow-dependent background error covariance calculated from ensemble forecasts and performs 4DVar optimization based on an incremental approach and a preconditioning algorithm. In Part II, the authors evaluated 4DEnVar with observing system simulation experiments (OSSEs) using the Advanced Research Weather Research and Forecasting Model (ARW-WRF, hereafter WRF). The current study extends the 4DEnVar to assimilate real observations for a cyclone in the Antarctic and the Southern Ocean in October 2007. The authors performed an intercomparison of four different WRF variational approaches for the case, including three-dimensional variational data assimilation (3DVar), first guess at the appropriate time (FGAT), and ensemble-based three-dimensional (En3DVar) and four-dimensional (4DEnVar) variational data assimilations. It is found that all data assimilation approaches produce positive impacts in this case. Applying the flow-dependent background error covariance in En3DVar and 4DEnVar yields forecast skills superior to those with the homogeneous and isotropic background error covariance in 3DVar and FGAT. In addition, the authors carried out FGAT and 4DEnVar 3-day cycling and 72-h forecasts. The results show that 4DEnVar produces a better performance in the cyclone prediction. The inflation factor on 4DEnVar can effectively improve the 4DEnVar analysis. The authors also conducted a short period (10-day lifetime of the cyclone in the domain) of analysis/forecast intercomparison experiments using 4DEnVar, FGAT, and 3DVar. The 4DEnVar scheme demonstrates overall superior and robust performance.


Sign in / Sign up

Export Citation Format

Share Document