Practical Ensemble-Based Approaches to Estimate Atmospheric Background Error Covariances for Limited-Area Deterministic Data Assimilation

2018 ◽  
Vol 146 (11) ◽  
pp. 3717-3733 ◽  
Author(s):  
Joël Bédard ◽  
Mark Buehner ◽  
Jean-François Caron ◽  
Seung-Jong Baek ◽  
Luc Fillion

Abstract High-resolution flow-dependent background error covariances can allow for a better usage of dense observation networks in applications of data assimilation for numerical weather prediction. The generation of high-resolution ensembles, however, can be computationally cost prohibitive. In this study, practical and low-cost ensemble generation methods are presented and compared against both global and regional ensemble Kalman filters (G-EnKF and R-EnKF, respectively). The goal is to provide limited-area deterministic assimilation schemes with higher-resolution flow-dependent background error covariances that perform at least as well as those from the G-EnKF when assimilating the same observations. The low-cost methods are based on short-range regional ensemble forecasts initialized from 1) deterministic analysis plus balanced perturbations (filter free approach) and 2) a simplified ensemble square root filter (S-EnSRF), centered on deterministic analyses. The resulting ensembles from the different approaches are used within a 4D ensemble–variational (4D-EnVar) assimilation system covering most of Canada and the northern United States. Diagnostic results show that the mean is an important component of the ensembles. Results also show that the persistence of the homogeneous characteristics of the perturbations in the filter free approach makes this method unsuited for short assimilation time windows since some error structures take longer to develop. The S-EnSRF approach overcomes this limitation by recycling part of the prior perturbations. Results from 1-month assimilation experiments show that the S-EnSRF and R-EnKF experiments provide forecasts of similar quality to those from G-EnKF. Furthermore, results from precipitation verification indicate that the R-EnKF experiment provides the best precipitation accumulation predictions over 24-h periods.

2011 ◽  
Vol 139 (7) ◽  
pp. 2025-2045 ◽  
Author(s):  
Zhiyong Meng ◽  
Fuqing Zhang

Abstract Ensemble-based data assimilation is a state estimation technique that uses short-term ensemble forecasts to estimate flow-dependent background error covariance and is best known by varying forms of ensemble Kalman filters (EnKFs). The EnKF has recently emerged as one of the primary alternatives to the variational data assimilation methods widely used in both global and limited-area numerical weather prediction models. In addition to comparing the EnKF with variational methods, this article reviews recent advances and challenges in the development and applications of the EnKF, including its hybrid with variational methods, in limited-area models that resolve weather systems from convective to meso- and regional scales.


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.


2018 ◽  
Vol 11 (1) ◽  
pp. 257-281 ◽  
Author(s):  
Piet Termonia ◽  
Claude Fischer ◽  
Eric Bazile ◽  
François Bouyssel ◽  
Radmila Brožková ◽  
...  

Abstract. The ALADIN System is a numerical weather prediction (NWP) system developed by the international ALADIN consortium for operational weather forecasting and research purposes. It is based on a code that is shared with the global model IFS of the ECMWF and the ARPEGE model of Météo-France. Today, this system can be used to provide a multitude of high-resolution limited-area model (LAM) configurations. A few configurations are thoroughly validated and prepared to be used for the operational weather forecasting in the 16 partner institutes of this consortium. These configurations are called the ALADIN canonical model configurations (CMCs). There are currently three CMCs: the ALADIN baseline CMC, the AROME CMC and the ALARO CMC. Other configurations are possible for research, such as process studies and climate simulations. The purpose of this paper is (i) to define the ALADIN System in relation to the global counterparts IFS and ARPEGE, (ii) to explain the notion of the CMCs, (iii) to document their most recent versions, and (iv) to illustrate the process of the validation and the porting of these configurations to the operational forecast suites of the partner institutes of the ALADIN consortium. This paper is restricted to the forecast model only; data assimilation techniques and postprocessing techniques are part of the ALADIN System but they are not discussed here.


2010 ◽  
Vol 10 (6) ◽  
pp. 1129-1149 ◽  
Author(s):  
M. Milelli ◽  
M. Turco ◽  
E. Oberto

Abstract. The forecast in areas of very complex topography, as for instance the Alpine region, is still a challenge even for the new generation of numerical weather prediction models which aim at reaching the km-scale. The problem is enhanced by a general lack of standard observations, which is even more evident over the southern side of the Alps. For this reason, it would be useful to increase the performance of the mathematical models by locally assimilating non-conventional data. Since in ARPA Piemonte there is the availability of a great number of non-GTS stations, it has been decided to assimilate the 2 m temperature, coming from this dataset, in the very-high resolution version of the COSMO model, which has a horizontal resolution of about 3 km, more similar to the average resolution of the thermometers. Four different weather situations have been considered, ranging from spring to winter, from cloudy to clear sky. The aim of the work is to investigate the effects of the assimilation of non-GTS data in order to create an operational very high-resolution analysis, but also to test the option of running in the future a very short-range forecast starting from these analyses (RUC or Rapid Update Cycle). The results, in terms of Root Mean Square Error, Mean Error and diurnal cycle of some surface variables such as 2 m temperature, 2 m relative humidity and 10 m wind intensity show a positive impact during the assimilation cycle which tends to dissipate a few hours after the end of it. Moreover, the 2 m temperature assimilation has a slightly positive or neutral impact on the vertical profiles of temperature, eventhough some calibration is needed for the precipitation field which is too much perturbed during the assimilation cycle, while it is unaffected in the forecast period. So the stability of the planetary boundary layer, on the one hand, has not been particularly improved by the new-data assimilation, but, on the other hand, it has not been destroyed. It has to be pointed out that a correct description of the planetary boundary layer, even only the lowest part of it, could be helpful to the forecasters and, in general, to the users, in order to deal with meteorological hazards such as snow (in particular snow/rain limit definition), or fog (description of temperature inversions).


2015 ◽  
Vol 144 (1) ◽  
pp. 149-169 ◽  
Author(s):  
Juanzhen Sun ◽  
Hongli Wang ◽  
Wenxue Tong ◽  
Ying Zhang ◽  
Chung-Yi Lin ◽  
...  

Abstract The momentum variables of streamfunction and velocity potential are used as control variables in a number of operational variational data assimilation systems. However, in this study it is shown that, for limited-area high-resolution data assimilation, the momentum control variables ψ and χ (ψχ) pose potential difficulties in background error modeling and, hence, may result in degraded analysis and forecast when compared with the direct use of x and y components of wind (UV). In this study, the characteristics of the modeled background error statistics, derived from an ensemble generated from Weather Research and Forecasting (WRF) Model real-time forecasts of two summer months, are first compared between the two control variable options. Assimilation and forecast experiments are then conducted with both options for seven convective events in a domain that encompasses the Rocky Mountain Front Range using the three-dimensional variational data assimilation (3DVar) system of the WRF Model. The impacts of the two control variable options are compared in terms of their skills in short-term qualitative precipitation forecasts. Further analysis is performed for one case to examine the impacts when radar observations are included in the 3DVar assimilation. The main findings are as follows: 1) the background error modeling used in WRF 3DVar with the control variables ψχ increases the length scale and decreases the variance for u and υ, which causes negative impact on the analysis of the velocity field and on precipitation prediction; 2) the UV-based 3DVar allows closer fits to radar wind observations; and 3) the use of UV control variables improves the 0–12-h precipitation prediction.


2020 ◽  
Vol 35 (2) ◽  
pp. 309-324
Author(s):  
Susan Rennie ◽  
Lawrence Rikus ◽  
Nathan Eizenberg ◽  
Peter Steinle ◽  
Monika Krysta

Abstract The impact of Doppler radar wind observations on forecasts from a developmental, high-resolution numerical weather prediction (NWP) system is assessed. The new 1.5-km limited-area model will be Australia’s first such operational NWP system to include data assimilation. During development, the assimilation of radar wind observations was trialed over a 2-month period to approve the initial inclusion of these observations. Three trials were run: the first with no radar data, the second with radial wind observations from precipitation echoes, and the third with radial winds from both precipitation and insect echoes. The forecasts were verified against surface observations from automatic weather stations, against rainfall accumulations using fractions skill scores, and against satellite cloud observations. These methods encompassed verification across a range of vertical levels. Additionally, a case study was examined more closely. Overall results showed little statistical difference in skill between the trials, and the net impact was neutral. While the new observations clearly affected the forecast, the objective and subjective analyses showed a neutral impact on the forecast overall. As a first step, this result is satisfactory for the operational implementation. In future, upgrades to the radar network will start to reduce the observation error, and further improvements to the data assimilation are planned, which may be expected to improve the impact.


Author(s):  
Magnus Lindskog ◽  
Adam Dybbroe ◽  
Roger Randriamampianina

AbstractMetCoOp is a Nordic collaboration on operational Numerical Weather Prediction based on a common limited-area km-scale ensemble system. The initial states are produced using a 3-dimensional variational data assimilation scheme utilizing a large amount of observations from conventional in-situ measurements, weather radars, global navigation satellite system, advanced scatterometer data and satellite radiances from various satellite platforms. A version of the forecasting system which is aimed for future operations has been prepared for an enhanced assimilation of microwave radiances. This enhanced data assimilation system will use radiances from the Microwave Humidity Sounder, the Advanced Microwave Sounding Unit-A and the Micro-Wave Humidity Sounder-2 instruments on-board the Metop-C and Fengyun-3 C/D polar orbiting satellites. The implementation process includes channel selection, set-up of an adaptive bias correction procedure, and careful monitoring of data usage and quality control of observations. The benefit of the additional microwave observations in terms of data coverage and impact on analyses, as derived using the degree of freedom of signal approach, is demonstrated. A positive impact on forecast quality is shown, and the effect on the precipitation for a case study is examined. Finally, the role of enhanced data assimilation techniques and adaptions towards nowcasting are discussed.


2008 ◽  
Vol 136 (3) ◽  
pp. 945-963 ◽  
Author(s):  
Jidong Gao ◽  
Ming Xue

Abstract A new efficient dual-resolution (DR) data assimilation algorithm is developed based on the ensemble Kalman filter (EnKF) method and tested using simulated radar radial velocity data for a supercell storm. Radar observations are assimilated on both high-resolution and lower-resolution grids using the EnKF algorithm with flow-dependent background error covariances estimated from the lower-resolution ensemble. It is shown that the flow-dependent and dynamically evolved background error covariances thus estimated are effective in producing quality analyses on the high-resolution grid. The DR method has the advantage of being able to significantly reduce the computational cost of the EnKF analysis. In the system, the lower-resolution ensemble provides the flow-dependent background error covariance, while the single-high-resolution forecast and analysis provides the benefit of higher resolution, which is important for resolving the internal structures of thunderstorms. The relative smoothness of the covariance obtained from the lower 4-km-resolution ensemble does not appear to significantly degrade the quality of analysis. This is because the cross covariance among different variables is of first-order importance for “retrieving” unobserved variables from the radar radial velocity data. For the DR analysis, an ensemble size of 40 appears to be a reasonable choice with the use of a 4-km horizontal resolution in the ensemble and a 1-km resolution in the high-resolution analysis. Several sensitivity tests show that the DR EnKF system is quite robust to different observation errors. A 4-km thinned data resolution is a compromise that is acceptable under the constraint of real-time applications. A data density of 8 km leads to a significant degradation in the analysis.


2018 ◽  
Vol 99 (7) ◽  
pp. 1415-1432 ◽  
Author(s):  
Yong Wang ◽  
Martin Belluš ◽  
Andrea Ehrlich ◽  
Máté Mile ◽  
Neva Pristov ◽  
...  

AbstractThis paper describes 27 years of scientific and operational achievement of Regional Cooperation for Limited Area Modelling in Central Europe (RC LACE), which is supported by the national (hydro-) meteorological services of Austria, Croatia, the Czech Republic, Hungary, Romania, Slovakia, and Slovenia. The principal objectives of RC LACE are to 1) develop and operate the state-of-the-art limited-area model and data assimilation system in the member states and 2) conduct joint scientific and technical research to improve the quality of the forecasts.In the last 27 years, RC LACE has contributed to the limited-area Aire Limitée Adaptation Dynamique Développement International (ALADIN) system in the areas of preprocessing of observations, data assimilation, model dynamics, physical parameterizations, mesoscale and convection-permitting ensemble forecasting, and verification. It has developed strong collaborations with numerical weather prediction (NWP) consortia ALADIN, the High Resolution Limited Area Model (HIRLAM) group, and the European Centre for Medium-Range Weather Forecasts (ECMWF). RC LACE member states exchange their national observations in real time and operate a common system that provides member states with the preprocessed observations for data assimilation and verification. RC LACE runs operationally a common mesoscale ensemble system, ALADIN–Limited Area Ensemble Forecasting (ALADIN-LAEF), over all of Europe for early warning of severe weather.RC LACE has established an extensive regional scientific and technical collaboration in the field of operational NWP for weather research, forecasting, and applications. Its 27 years of experience have demonstrated the value of regional cooperation among small- and medium-sized countries for success in the development of a modern forecasting system, knowledge transfer, and capacity building.


2011 ◽  
Vol 139 (6) ◽  
pp. 2008-2024 ◽  
Author(s):  
Brian C. Ancell ◽  
Clifford F. Mass ◽  
Gregory J. Hakim

Abstract Previous research suggests that an ensemble Kalman filter (EnKF) data assimilation and modeling system can produce accurate atmospheric analyses and forecasts at 30–50-km grid spacing. This study examines the ability of a mesoscale EnKF system using multiscale (36/12 km) Weather Research and Forecasting (WRF) model simulations to produce high-resolution, accurate, regional surface analyses, and 6-h forecasts. This study takes place over the complex terrain of the Pacific Northwest, where the small-scale features of the near-surface flow field make the region particularly attractive for testing an EnKF and its flow-dependent background error covariances. A variety of EnKF experiments are performed over a 5-week period to test the impact of decreasing the grid spacing from 36 to 12 km and to evaluate new approaches for dealing with representativeness error, lack of surface background variance, and low-level bias. All verification in this study is performed with independent, unassimilated observations. Significant surface analysis and 6-h forecast improvements are found when EnKF grid spacing is reduced from 36 to 12 km. Forecast improvements appear to be a consequence of increased resolution during model integration, whereas analysis improvements also benefit from high-resolution ensemble covariances during data assimilation. On the 12-km domain, additional analysis improvements are found by reducing observation error variance in order to address representativeness error. Removing model surface biases prior to assimilation significantly enhances the analysis. Inflating surface wind and temperature background error variance has large impacts on analyses, but only produces small improvements in analysis RMS errors. Both surface and upper-air 6-h forecasts are nearly unchanged in the 12-km experiments. Last, 12-km WRF EnKF surface analyses and 6-h forecasts are shown to generally outperform those of the Global Forecast System (GFS), North American Model (NAM), and the Rapid Update Cycle (RUC) by about 10%–30%, although these improvements do not extend above the surface. Based on these results, future improvements in multiscale EnKF are suggested.


Sign in / Sign up

Export Citation Format

Share Document