Inhomogeneous Background Error Modeling and Estimation over Antarctica

2010 ◽  
Vol 138 (6) ◽  
pp. 2229-2252 ◽  
Author(s):  
Yann Michel ◽  
Thomas Auligné

Abstract The structure of the analysis increments in a variational data assimilation scheme is strongly driven by the formulation of the background error covariance matrix, especially in data-sparse areas such as the Antarctic region. The gridpoint background error modeling in this study makes use of regression-based balance operators between variables, empirical orthogonal function decomposition to define the vertical correlations, gridpoint variances, and high-order efficient recursive filters to impose horizontal correlations. A particularity is that the regression operators and the recursive filters have been made spatially inhomogeneous. The computation of the background error statistics is performed with the Weather Research and Forecast (WRF) model from a set of forecast differences. The mesoscale limited-area domains of interest cover Antarctica. Inhomogeneities of background errors are shown to be related to the particular orography and physics of the area. Differences seem particularly pronounced between ocean and land boundary layers.

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.


WRF model have been tuned and tested over Georgia’s territory for years. First time in Georgia theprocess of data assimilation in Numerical weather prediction is developing. This work presents how forecasterror statistics appear in the data assimilation problem through the background error covariance matrix – B, wherethe variances and correlations associated with model forecasts are estimated. Results of modeling of backgrounderror covariance matrix for control variables using WRF model over Georgia with desired domain configurationare discussed and presented. The modeling was implemented in two different 3DVAR systems (WRFDA andGSI) and results were checked by pseudo observation benchmark cases using also default global and regional BEmatrixes. The mathematical and physical properties of the covariances are also reviewed.


2018 ◽  
Vol 33 (2) ◽  
pp. 561-582 ◽  
Author(s):  
Kuan-Jen Lin ◽  
Shu-Chih Yang ◽  
Shuyi S. Chen

Abstract Ensemble-based data assimilation (EDA) has been used for tropical cyclone (TC) analysis and prediction with some success. However, the TC position spread determines the structure of the TC-related background error covariance and affects the performance of EDA. With an idealized experiment and a real TC case study, it is demonstrated that observations in the core region cannot be optimally assimilated when the TC position spread is large. To minimize the negative impact from large position uncertainty, a TC-centered EDA approach is implemented in the Weather Research and Forecasting (WRF) Model–local ensemble transform Kalman filter (WRF-LETKF) assimilation system. The impact of TC-centered EDA on TC analysis and prediction of Typhoon Fanapi (2010) is evaluated. Using WRF Model nested grids with 4-km grid spacing in the innermost domain, the focus is on EDA using dropsonde data from the Impact of Typhoons on the Ocean in the Pacific field campaign. The results show that the TC structure in the background mean state is improved and that unrealistically large ensemble spread can be alleviated. The characteristic horizontal scale in the background error covariance is smaller and narrower compared to those derived from the conventional EDA approach. Storm-scale corrections are improved using dropsonde data, which is more favorable for TC development. The analysis using the TC-centered EDA is in better agreement with independent observations. The improved analysis ameliorates model shock and improves the track forecast during the first 12 h and landfall at 72 h. The impact on intensity prediction is mixed with a better minimum sea level pressure and overestimated peak winds.


2014 ◽  
Vol 21 (5) ◽  
pp. 971-985 ◽  
Author(s):  
C. Cardinali ◽  
N. Žagar ◽  
G. Radnoti ◽  
R. Buizza

Abstract. The paper investigates a method to represent model error in the ensemble data assimilation (EDA) system. The ECMWF operational EDA simulates the effect of both observations and model uncertainties. Observation errors are represented by perturbations with statistics characterized by the observation error covariance matrix whilst the model uncertainties are represented by stochastic perturbations added to the physical tendencies to simulate the effect of random errors in the physical parameterizations (ST-method). In this work an alternative method (XB-method) is proposed to simulate model uncertainties by adding perturbations to the model background field. In this way the error represented is not just restricted to model error in the usual sense but potentially extends to any form of background error. The perturbations have the same correlation as the background error covariance matrix and their magnitude is computed from comparing the high-resolution operational innovation variances with the ensemble variances when the ensemble is obtained by perturbing only the observations (OBS-method). The XB-method has been designed to represent the short range model error relevant for the data assimilation window. Spread diagnostic shows that the XB-method generates a larger spread than the ST-method that is operationally used at ECMWF, in particular in the extratropics. Three-dimensional normal-mode diagnostics indicate that XB-EDA spread projects more than the spread from the other EDAs onto the easterly inertia-gravity modes associated with equatorial Kelvin waves, tropical dynamics and, in general, model error sources. The background error statistics from the above described EDAs have been employed in the assimilation system. The assimilation system performance showed that the XB-method background error statistics increase the observation influence in the analysis process. The other EDA background error statistics, when inflated by a global factor, generate analyses with 30–50% smaller degree of freedom of signal. XB-EDA background error variances have not been inflated. The presented EDAs have been used to generate the initial perturbations of the ECMWF ensemble prediction system (EPS) of which the XB-EDA induces the largest EPS spread, also in the medium range, leading to a more reliable ensemble. Compared to ST-EDA, XB-EDA leads to a small improvement of the EPS ignorance skill score at day 3 and 7.


2006 ◽  
Vol 3 (6) ◽  
pp. 1977-1998 ◽  
Author(s):  
S. Dobricic ◽  
N. Pinardi ◽  
M. Adani ◽  
M. Tonani ◽  
C. Fratianni ◽  
...  

Abstract. This study presents the upgrade of the Optimal Interpolation scheme used in the basin scale assimilation scheme of the Mediterranean Forecasting System. The modifications include a daily analysis cycle, the assimilation of ARGO float profiles, the implementation of the geostrophic balance in the background error covariance matrix and the initialisation of the analyses. A series of numerical experiments showed that each modification had a positive impact on the accuracy of the analyses: The daily cycle improved the representation of the processes with a relatively high temporal variability, the assimilation of ARGO floats profiles significantly improved the salinity analyses quality, the geostrophically balanced background error covariances improved the accuracy of the surface elevation analyses, and the initialisation removed the barotropic adjustment in the forecast first time steps starting from the analysis.


2020 ◽  
Vol 20 (14) ◽  
pp. 8473-8500
Author(s):  
Yohanna Villalobos ◽  
Peter Rayner ◽  
Steven Thomas ◽  
Jeremy Silver

Abstract. This paper addresses the question of how much uncertainties in CO2 fluxes over Australia can be reduced by assimilation of total-column carbon dioxide retrievals from the Orbiting Carbon Observatory-2 (OCO-2) satellite instrument. We apply a four-dimensional variational data assimilation system, based around the Community Multiscale Air Quality (CMAQ) transport-dispersion model. We ran a series of observing system simulation experiments to estimate posterior error statistics of optimized monthly-mean CO2 fluxes in Australia. Our assimilations were run with a horizontal grid resolution of 81 km using OCO-2 data for 2015. Based on four representative months, we find that the integrated flux uncertainty for Australia is reduced from 0.52 to 0.13 Pg C yr−1. Uncertainty reductions of up to 90 % were found at grid-point resolution over productive ecosystems. Our sensitivity experiments show that the choice of the correlation structure in the prior error covariance plays a large role in distributing information from the observations. We also found that biases in the observations would significantly impact the inverted fluxes and could contaminate the final results of the inversion. Biases in prior fluxes are generally removed by the inversion system. Biases in the boundary conditions have a significant impact on retrieved fluxes, but this can be mitigated by including boundary conditions in our retrieved parameters. In general, results from our idealized experiments suggest that flux inversions at this unusually fine scale will yield useful information on the carbon cycle at continental and finer scales.


2016 ◽  
Author(s):  
Paolo Oddo ◽  
Andrea Storto ◽  
Srdjan Dobricic ◽  
Aniello Russo ◽  
Craig Lewis ◽  
...  

Abstract. A hybrid variational-ensemble data assimilation scheme to estimate the vertical and horizontal parts of the background-error covariance matrix for an ocean variational data assimilation system is presented and tested in a limited area ocean model implemented in the western Mediterranean Sea. An extensive dataset collected during the Recognized Environmental Picture Experiments conducted in June 2014 by the Centre for Maritime Research and Experimentation has been used for assimilation and validation. The hybrid scheme is used to both correct the systematic error introduced in the system from the external forcing (initial, lateral and surface open boundary conditions) and model parameterization and improve the representation of small scale errors in the Background Error Covariance matrix. An ensemble system is run off-line for further use in the hybrid scheme, generated through perturbation of assimilated observations. Results of four different experiments have been compared. The reference experiment uses the classical static formulation of the background error covariance matrix and has no systematic error correction. The other three experiments account, or not, for systematic error correction and hybrid daily estimates of the background error covariance matrix combining the static and the ensemble derived errors statistics. Results show that the hybrid scheme when used in conjunction with the systematic error correction reduce the mean absolute error of temperature and salinity misfit by 55 % and 42 % respectively versus statistics arising from standard climatological covariances without systematic error correction.


2014 ◽  
Vol 53 (10) ◽  
pp. 2287-2309 ◽  
Author(s):  
Hongli Wang ◽  
Xiang-Yu Huang ◽  
Juanzhen Sun ◽  
Dongmei Xu ◽  
Man Zhang ◽  
...  

AbstractBackground error modeling plays a key role in a variational data assimilation system. The National Meteorological Center (NMC) method has been widely used in variational data assimilation systems to generate a forecast error ensemble from which the climatological background error covariance can be modeled. In this paper, the characteristics of the background error modeling via the NMC method are investigated for the variational data assimilation system of the Weather Research and Forecasting (WRF-Var) Model. The background error statistics are extracted from short-term 3-km-resolution forecasts in June, July, and August 2012 over a limited-area domain. It is found 1) that background error variances vary from month to month and also have a feature of diurnal variations in the low-level atmosphere and 2) that u- and υ-wind variances are underestimated and their autocorrelation length scales are overestimated when the default control variable option in WRF-Var is used. A new approach of control variable transform (CVT) is proposed to model the background error statistics based on the NMC method. The new approach is capable of extracting inhomogeneous and anisotropic climatological information from the forecast error ensemble obtained via the NMC method. Single observation assimilation experiments show that the proposed method not only has the merit of incorporating geographically dependent covariance information, but also is able to produce a multivariate analysis. The results from the data assimilaton and forecast study of a real convective case show that the use of the new CVT improves synoptic weather system and precipitation forecasts for up to 12 h.


2007 ◽  
Vol 135 (4) ◽  
pp. 1506-1521 ◽  
Author(s):  
Haixia Liu ◽  
Ming Xue ◽  
R. James Purser ◽  
David F. Parrish

Abstract Anisotropic recursive filters are implemented within a three-dimensional variational data assimilation (3DVAR) framework to efficiently model the effect of flow-dependent background error covariance. The background error covariance is based on an estimated error field and on the idea of Riishøjgaard. In the anisotropic case, the background error pattern can be stretched or flattened in directions oblique to the alignment of the grid coordinates and is constructed by applying, at each point, six recursive filters along six directions corresponding, in general, to a special configuration of oblique lines of the grid. The recursive filters are much more efficient than corresponding explicit filters used in an earlier study and are therefore more suitable for real-time numerical weather prediction. A set of analysis experiments are conducted at a mesoscale resolution to examine the effectiveness of the 3DVAR system in analyzing simulated global positioning system (GPS) slant-path water vapor observations from ground-based GPS receivers and observations from collocated surface stations. It is shown that the analyses produced with recursive filters are at least as good as those with corresponding explicit filters. In some cases, the recursive filters actually perform better. The impact of flow-dependent background errors modeled using the anisotropic recursive filters is also examined. The use of anisotropic filters improves the analysis, especially in terms of finescale structures. The analysis system is found to be effective in the presence of typical observational errors. The sensitivity of isotropic and anisotropic recursive-filter analyses to the decorrelation scales is also examined systematically.


Ocean Science ◽  
2007 ◽  
Vol 3 (1) ◽  
pp. 149-157 ◽  
Author(s):  
S. Dobricic ◽  
N. Pinardi ◽  
M. Adani ◽  
M. Tonani ◽  
C. Fratianni ◽  
...  

Abstract. This study presents the upgrade of the Optimal Interpolation scheme used in the basin scale assimilation scheme of the Mediterranean Forecasting System. The modifications include a daily analysis cycle, the assimilation of ARGO float profiles, the implementation of the geostrophic balance in the background error covariance matrix and the initialisation of the analyses. A series of numerical experiments showed that each modification had a positive impact on the accuracy of the analyses: The daily cycle improved the representation of the processes with the temporal variability shorter than a week, the assimilation of ARGO floats profiles significantly improved the salinity analyses quality, the geostrophically balanced background error covariances improved the accuracy of the surface elevation analyses, and the initialisation removed the barotropic adjustment in the forecast first time steps starting from the analysis.


Sign in / Sign up

Export Citation Format

Share Document