scholarly journals Source backtracking for dust storm emission inversion using adjoint method: case study of northeast China

2020 ◽  
Author(s):  
Jianbing Jin ◽  
Arjo Segers ◽  
Hong Liao ◽  
Arnold Heemink ◽  
Richard Kranenburg ◽  
...  
Keyword(s):  
Data Assimilation ◽  
Dust Storm ◽  
Northeast China ◽  
Source Region ◽  
Background Error ◽  
Dust Transport ◽  
Error Covariance ◽  
Trace Back ◽  
Using Data ◽  
Sensitivity Method

Abstract. Emission inversion using data assimilation fundamentally relies on having the correct assumptions on the emission background error covariance. A perfect covariance accounts for the uncertainty based on prior knowledge, and is able to explain differences between model simulations and observations. In practice, emission uncertainties are constructed empirically, hence a partially unrepresentative covariance is unavoidable. Concerning its complex parameterization, dust emissions are a typical example where the uncertainty could be induced from many underlying inputs, e.g., information on soil composition and moisture, landcover and erosive wind velocity, and these can hardly be taken into account together. This paper describes how an adjoint model can be used to detect errors in the emission uncertainty assumptions. This adjoint based sensitivity method could serve as a supplement of a data assimilation inverse modeling system to trace back the error sources, in case that large observation-minus-simulation residues remain after assimilation based on empirical background covariance. The method follows on application of a data assimilation emission inversion for an extreme severe dust storm over East Asia (Jin et al., 2019b). The assimilation system successfully resolved observation-minus-simulation errors using satellite AOD observations in most of the dust-affected regions. However, a large underestimation of dust in northeast China remained despite the fact the assimilated measurements indicated severe dust plumes there. An adjoint implementation of our dust simulation model is then used to detect the most likely source region for these unresolved dust loads. The backward modeling points to the Horqin desert as source region, which was indicated as a non-source region by the existing emission scheme. The reference emission and uncertainty are then reconstructed over the Horqin desert by assuming higher surface erodibility. After the emission reconstruction, the emission inversion is performed again and the posterior dust simulations and reality are now in much closer harmony. Based on our results, it is advised that emission sources in dust transport models include Horqin desert as a more active source region.

scholarly journals Source backtracking for dust storm emission inversion using an adjoint method: case study of Northeast China

2020 ◽  
Vol 20 (23) ◽  
pp. 15207-15225
Author(s):  
Jianbing Jin ◽  
Arjo Segers ◽  
Hong Liao ◽  
Arnold Heemink ◽  
Richard Kranenburg ◽  
...  
Keyword(s):  
Data Assimilation ◽  
Dust Storm ◽  
Northeast China ◽  
Source Region ◽  
Background Error ◽  
Dust Transport ◽  
Error Covariance ◽  
Trace Back ◽  
Using Data ◽  
Sensitivity Method

Abstract. Emission inversion using data assimilation fundamentally relies on having the correct assumptions about the emission background error covariance. A perfect covariance accounts for the uncertainty based on prior knowledge and is able to explain differences between model simulations and observations. In practice, emission uncertainties are constructed empirically; hence, a partially unrepresentative covariance is unavoidable. Concerning its complex parameterization, dust emissions are a typical example where the uncertainty could be induced from many underlying inputs, e.g., information on soil composition and moisture, land cover and erosive wind velocity, and these can hardly be taken into account together. This paper describes how an adjoint model can be used to detect errors in the emission uncertainty assumptions. This adjoint-based sensitivity method could serve as a supplement of a data assimilation inverse modeling system to trace back the error sources in case large observation-minus-simulation residues remain after assimilation based on empirical background covariance. The method follows an application of a data assimilation emission inversion for an extreme severe dust storm over East Asia (Jin et al., 2019b). The assimilation system successfully resolved observation-minus-simulation errors using satellite AOD observations in most of the dust-affected regions. However, a large underestimation of dust in Northeast China remained despite the fact that the assimilated measurements indicated severe dust plumes there. An adjoint implementation of our dust simulation model is then used to detect the most likely source region for these unresolved dust loads. The backward modeling points to the Horqin desert as the source region, which was indicated as a non-source region by the existing emission scheme. The reference emission and uncertainty are then reconstructed over the Horqin desert by assuming higher surface erodibility. After the emission reconstruction, the emission inversion is performed again, and the posterior dust simulations and reality are now in much closer harmony. Based on our results, it is advised that emission sources in dust transport models include the Horqin desert as a more active source region.

Aerosol Data Assimilation Using Data from Fengyun-3A and MODIS: Application to a Dust Storm over East Asia in 2011

2018 ◽  
Vol 36 (1) ◽  
pp. 1-14 ◽  
Author(s):  
Xiaoli Xia ◽  
Jinzhong Min ◽  
Feifei Shen ◽  
Yuanbing Wang ◽  
Chun Yang
Keyword(s):  
Data Assimilation ◽  
East Asia ◽  
Dust Storm ◽  
Using Data

scholarly journals Background Error in WRF Model

2020 ◽  
Vol 16 ◽  
Keyword(s):  
Data Assimilation ◽  
Physical Properties ◽  
Covariance Matrix ◽  
Numerical Weather Prediction ◽  
Weather Prediction ◽  
Wrf Model ◽  
Background Error ◽  
Numerical Weather ◽  
Error Covariance ◽  
First Time

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.

Reducing TC Position Uncertainty in an Ensemble Data Assimilation and Prediction System: A Case Study of Typhoon Fanapi (2010)

2018 ◽  
Vol 33 (2) ◽  
pp. 561-582 ◽  
Author(s):  
Kuan-Jen Lin ◽  
Shu-Chih Yang ◽  
Shuyi S. Chen
Keyword(s):  
Data Assimilation ◽  
Negative Impact ◽  
Wrf Model ◽  
Track Forecast ◽  
Background Error ◽  
Position Uncertainty ◽  
Background Error Covariance ◽  
Error Covariance ◽  
The Impact

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.

scholarly journals Assimilation of Stratospheric Temperature and Ozone with an Ensemble Kalman Filter in a Chemistry–Climate Model

2011 ◽  
Vol 139 (11) ◽  
pp. 3389-3404 ◽  
Author(s):  
Thomas Milewski ◽  
Michel S. Bourqui
Keyword(s):  
Kalman Filter ◽  
Data Assimilation ◽  
Covariance Matrix ◽  
Ensemble Kalman Filter ◽  
Climate Model ◽  
Observation Error ◽  
Background Error ◽  
Error Covariance Matrix ◽  
Error Covariance ◽  
Analysis Error

Abstract A new stratospheric chemical–dynamical data assimilation system was developed, based upon an ensemble Kalman filter coupled with a Chemistry–Climate Model [i.e., the intermediate-complexity general circulation model Fast Stratospheric Ozone Chemistry (IGCM-FASTOC)], with the aim to explore the potential of chemical–dynamical coupling in stratospheric data assimilation. The system is introduced here in a context of a perfect-model, Observing System Simulation Experiment. The system is found to be sensitive to localization parameters, and in the case of temperature (ozone), assimilation yields its best performance with horizontal and vertical decorrelation lengths of 14 000 km (5600 km) and 70 km (14 km). With these localization parameters, the observation space background-error covariance matrix is underinflated by only 5.9% (overinflated by 2.1%) and the observation-error covariance matrix by only 1.6% (0.5%), which makes artificial inflation unnecessary. Using optimal localization parameters, the skills of the system in constraining the ensemble-average analysis error with respect to the true state is tested when assimilating synthetic Michelson Interferometer for Passive Atmospheric Sounding (MIPAS) retrievals of temperature alone and ozone alone. It is found that in most cases background-error covariances produced from ensemble statistics are able to usefully propagate information from the observed variable to other ones. Chemical–dynamical covariances, and in particular ozone–wind covariances, are essential in constraining the dynamical fields when assimilating ozone only, as the radiation in the stratosphere is too slow to transfer ozone analysis increments to the temperature field over the 24-h forecast window. Conversely, when assimilating temperature, the chemical–dynamical covariances are also found to help constrain the ozone field, though to a much lower extent. The uncertainty in forecast/analysis, as defined by the variability in the ensemble, is large compared to the analysis error, which likely indicates some amount of noise in the covariance terms, while also reducing the risk of filter divergence.

An improved variational Data Assimilation method for ocean models with limited number of observations.

2020 ◽  
Author(s):  
Lewis Sampson ◽  
Jose M. Gonzalez-Ondina ◽  
Georgy Shapiro
Keyword(s):  
Data Assimilation ◽  
Covariance Matrix ◽  
Observational Data ◽  
Inner Product ◽  
Product Analysis ◽  
Background Error ◽  
Error Covariance Matrix ◽  
Error Covariance ◽  
Ocean Prediction ◽  
The Stability

<p>Data assimilation (DA) is a critical component for most state-of-the-art ocean prediction systems, which optimally combines model data and observational measurements to obtain an improved estimate of the modelled variables, by minimizing a cost function. The calculation requires the knowledge of the background error covariance matrix (BECM) as a weight for the quality of the model results, and an observational error covariance matrix (OECM) which weights the observational data.</p><p>Computing the BECM would require knowing the true values of the physical variables, which is not feasible. Instead, the BECM is estimated from model results and observations by using methods like National Meteorological Centre (NMC) or the Hollingsworth and Lönnberg (1984) (H-L). These methods have some shortcomings which make them unfit in some situations, which includes being fundamentally one-dimensional and making a suboptimal use of observations.</p><p>We have produced a novel method for error estimation, using an analysis of observations minus background data (innovations), which attempts to improve on some of these shortcomings. In particular, our method better infers information from observations, requiring less data to produce statistically robust results. We do this by minimizing a linear combination of functions to fit the data using a specifically tailored inner product, referred to as an inner product analysis (IPA).</p><p>We are able to produce quality BECM estimations even in data sparse domains, with notably better results in conditions of scarce observational data. By using a sample of observations, with decreasing sample size, we show that the stability and efficiency of our method, when compared to that of the H-L approach, does not deteriorate nearly as much as the number of data points decrease. We have found that we are able to continually produce error estimates with a reduced set of data, whereas the H-L method will begin to produce spurious values for smaller samples.</p><p>Our method works very well in combination with standard tools like NEMOVar by providing the required standard deviations and length-scales ratios. We have successfully ran this in the Arabian Sea for multiple seasons and compared the results with the H-L (in optimal conditions, when plenty of data is available), spatially the methods perform equally well. When we look at the root mean square error (RMSE) we see very similar performances, with each method giving better results for some seasons and worse for others.</p>

scholarly journals A Comparison of Multiscale GSI-Based EnKF and 3DVar Data Assimilation Using Radar and Conventional Observations for Midlatitude Convective-Scale Precipitation Forecasts

2015 ◽  
Vol 143 (8) ◽  
pp. 3087-3108 ◽  
Author(s):  
Aaron Johnson ◽  
Xuguang Wang ◽  
Jacob R. Carley ◽  
Louis J. Wicker ◽  
Christopher Karstens
Keyword(s):  
Data Assimilation ◽  
Multiscale Analysis ◽  
Multiple Scales ◽  
Three Dimensional ◽  
Forecast Period ◽  
Background Error ◽  
Cold Pools ◽  
Error Covariance ◽  
Convective Scale

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.

Sign in / Sign up

Export Citation Format

Share Document