scholarly journals GSI Three-Dimensional Ensemble–Variational Hybrid Data Assimilation Using a Global Ensemble for the Regional Rapid Refresh Model

2017 ◽  
Vol 145 (10) ◽  
pp. 4205-4225 ◽  
Author(s):  
Ming Hu ◽  
Stanley G. Benjamin ◽  
Therese T. Ladwig ◽  
David C. Dowell ◽  
Stephen S. Weygandt ◽  
...  

The Rapid Refresh (RAP) is an hourly updated regional meteorological data assimilation/short-range model forecast system running operationally at NOAA/National Centers for Environmental Prediction (NCEP) using the community Gridpoint Statistical Interpolation analysis system (GSI). This paper documents the application of the GSI three-dimensional hybrid ensemble–variational assimilation option to the RAP high-resolution, hourly cycling system and shows the skill improvements of 1–12-h forecasts of upper-air wind, moisture, and temperature over the purely three-dimensional variational analysis system. Use of perturbation data from an independent global ensemble, the Global Data Assimilation System (GDAS), is demonstrated to be very effective for the regional RAP hybrid assimilation. In this paper, application of the GSI-hybrid assimilation for the RAP is explained. Results from sensitivity experiments are shown to define configurations for the operational RAP version 2, the ratio of static and ensemble background error covariance, and vertical and horizontal localization scales for the operational RAP version 3. Finally, a 1-week RAP experiment from a summer period was performed using a global ensemble from a winter period, suggesting that a significant component of its multivariate covariance structure from the ensemble is independent of time matching between analysis time and ensemble valid time.

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.


2015 ◽  
Vol 2015 ◽  
pp. 1-17 ◽  
Author(s):  
Govindan Kutty ◽  
Xuguang Wang

The impact of observations can be dependent on many factors in a data assimilation (DA) system including data quality control, preprocessing, skill of the model, and the DA algorithm. The present study focuses on comparing the impacts of observations assimilated by two different DA algorithms. A three-dimensional ensemble-variational (3DEnsVar) hybrid data assimilation system was recently developed based on the Gridpoint Statistical Interpolation (GSI) data assimilation system and was implemented operationally for the National Center for Environmental Prediction (NCEP) Global Forecast System (GFS). One question to address is, how the impacts of observations on GFS forecasts differ when assimilated by the traditional GSI-three dimensional variational (3DVar) and the new 3DEnsVar. Experiments were conducted over a 6-week period during Northern Hemisphere winter season at a reduced resolution. For both the control and data denial experiments, the forecasts produced by 3DEnsVar were more accurate than GSI3DVar experiments. The results suggested that the observations were better and more effectively exploited to increment the background forecast in 3DEnsVar. On the other hand, in GSI3DVar, where the observation will be making mostly local, isotropic increments without proper flow dependent extrapolation is more sensitive to the number and types observations assimilated.


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.


2009 ◽  
Vol 137 (3) ◽  
pp. 1046-1060 ◽  
Author(s):  
Daryl T. Kleist ◽  
David F. Parrish ◽  
John C. Derber ◽  
Russ Treadon ◽  
Ronald M. Errico ◽  
...  

Abstract The gridpoint statistical interpolation (GSI) analysis system is a unified global/regional three-dimensional variational data assimilation (3DVAR) analysis code that has been under development for several years at the National Centers for Environmental Prediction (NCEP)/Environmental Modeling Center. It has recently been implemented into operations at NCEP in both the global and North American data assimilation systems (GDAS and NDAS, respectively). An important aspect of this development has been improving the balance of the analysis produced by GSI. The improved balance between variables has been achieved through the inclusion of a tangent-linear normal-mode constraint (TLNMC). The TLNMC method has proven to be very robust and effective. The TLNMC as part of the global GSI system has resulted in substantial improvement in data assimilation at NCEP.


2019 ◽  
Vol 147 (1) ◽  
pp. 135-151 ◽  
Author(s):  
Jean-François Caron ◽  
Yann Michel ◽  
Thibaut Montmerle ◽  
Étienne Arbogast

Following the recent development of a three-dimensional ensemble–variational (3DEnVar) data assimilation algorithm for the AROME-France NWP system, this paper examines different approaches to reduce the sampling noise in the ensemble-derived background error covariances in this new scheme without modifying the background ensemble generation strategy. We first examine two variants of scale-dependent localization: one method consists of applying different amounts of localization to different ranges of background error covariance spatial scales, while simultaneously assimilating all of the available observations. Another separate approach uses time-lagged forecasts in order to increase the effective ensemble size, up to a factor of 3 here. This approach of time-lagged forecasts is considered both on its own and together with scale-dependent localization. When the background error covariances are derived from the most recent 25-member ensemble forecasts, the results from data assimilation cycles over a 33-day winter period show that avoiding cross covariances between scales in the scale-dependent localization formulation first proposed by Buehner performs better than the more recent formulation of Buehner and Shlyaeva. However, when increasing the effective ensemble size to 75 members with time-lagged forecasts, the two scale-dependent formulations provide similar forecast improvements overall. It is also found that the lagged-members approach outperforms scale-dependent localization on its own. The largest forecast improvements are obtained when combining the two approaches.


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.


2012 ◽  
Vol 140 (2) ◽  
pp. 587-600 ◽  
Author(s):  
Meng Zhang ◽  
Fuqing Zhang

A hybrid data assimilation approach that couples the ensemble Kalman filter (EnKF) and four-dimensional variational (4DVar) methods is implemented for the first time in a limited-area weather prediction model. In this coupled system, denoted E4DVar, the EnKF and 4DVar systems run in parallel while feeding into each other. The multivariate, flow-dependent background error covariance estimated from the EnKF ensemble is used in the 4DVar minimization and the ensemble mean in the EnKF analysis is replaced by the 4DVar analysis, while updating the analysis perturbations for the next cycle of ensemble forecasts with the EnKF. Therefore, the E4DVar can obtain flow-dependent information from both the explicit covariance matrix derived from ensemble forecasts, as well as implicitly from the 4DVar trajectory. The performance of an E4DVar system is compared with the uncoupled 4DVar and EnKF for a limited-area model by assimilating various conventional observations over the contiguous United States for June 2003. After verifying the forecasts from each analysis against standard sounding observations, it is found that the E4DVar substantially outperforms both the EnKF and 4DVar during this active summer month, which featured several episodes of severe convective weather. On average, the forecasts produced from E4DVar analyses have considerably smaller errors than both of the stand-alone EnKF and 4DVar systems for forecast lead times up to 60 h.


2011 ◽  
Vol 139 (10) ◽  
pp. 3243-3247 ◽  
Author(s):  
Thomas M. Hamill ◽  
Jeffrey S. Whitaker ◽  
Daryl T. Kleist ◽  
Michael Fiorino ◽  
Stanley G. Benjamin

Abstract Experimental ensemble predictions of tropical cyclone (TC) tracks from the ensemble Kalman filter (EnKF) using the Global Forecast System (GFS) model were recently validated for the 2009 Northern Hemisphere hurricane season by Hamill et al. A similar suite of tests is described here for the 2010 season. Two major changes were made this season: 1) a reduction in the resolution of the GFS model, from 2009’s T384L64 (~31 km at 25°N) to 2010’s T254L64 (~47 km at 25°N), and some changes in model physics; and 2) the addition of a limited test of deterministic forecasts initialized from a hybrid three-dimensional variational data assimilation (3D-Var)/EnKF method. The GFS/EnKF ensembles continued to produce reduced track errors relative to operational ensemble forecasts created by the National Centers for Environmental Prediction (NCEP), the Met Office (UKMO), and the Canadian Meteorological Centre (CMC). The GFS/EnKF was not uniformly as skillful as the European Centre for Medium-Range Weather Forecasts (ECMWF) ensemble prediction system. GFS/EnKF track forecasts had slightly higher error than ECMWF at longer leads, especially in the western North Pacific, and exhibited poorer calibration between spread and error than in 2009, perhaps in part because of lower model resolution. Deterministic forecasts from the hybrid were competitive with deterministic EnKF ensemble-mean forecasts and superior in track error to those initialized from the operational variational algorithm, the Gridpoint Statistical Interpolation (GSI). Pending further successful testing, the National Oceanic and Atmospheric Administration (NOAA) intends to implement the global hybrid system operationally for data assimilation.


2010 ◽  
Vol 138 (2) ◽  
pp. 563-578 ◽  
Author(s):  
Jean-François Caron ◽  
Luc Fillion

Abstract The differences in the balance characteristics between dry and precipitation areas in estimated short-term forecast error fields are investigated. The motivation is to see if dry and precipitation areas need to be treated differently in atmospheric data assimilation systems. Using an ensemble of lagged forecast differences, it is shown that perturbations are, on average, farther away from geostrophic balance over precipitation areas than over dry areas and that the deviation from geostrophic balance is proportional to the intensity of precipitation. Following these results, the authors investigate whether some improvements in the coupling between mass and rotational wind increments over precipitation areas can be achieved by using only the precipitation points within an ensemble of estimated forecast errors to construct a so-called diabatic balance operator by linear regression. Comparisons with a traditional approach to construct balance operators by linear regression show that the new approach leads to a gradually significant improvement (related to the intensity of the diabatic processes) of the accuracy of the coupling over precipitation areas as judged from an ensemble of lagged forecast differences. Results from a series of simplified data assimilation experiments show that the new balance operators can produce analysis increments that are substantially different from those associated with the traditional balance operator, particularly for observations located in the lower atmosphere. Issues concerning the implementation of this new approach in a full-fledged analysis system are briefly discussed but their investigations are left for a following study.


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.


Sign in / Sign up

Export Citation Format

Share Document