Evaluating the Impact of Assimilating Aerosol Optical Depth Observations on Dust Forecasts Over North Africa and the East Atlantic Using Different Data Assimilation Methods

Performance comparisons of the three data assimilation methods for improved predictability of PM2·5: Ensemble Kalman filter, ensemble square root filter, and three-dimensional variational methods

To improve the predictability of concentrations of atmospheric particulate matter, a data assimilation (DA) system using ensemble square root filter (EnSRF) has been developed for the Community Multiscale Air Quality (CMAQ) model. The EnSRF DA method is a deterministic variant of the ensemble Kalman filter (EnKF) method, which means that unlike the EnKF method, it does not add random noise to the observations. To compare the performances of the EnSRF with those of other DA methods, such as EnKF and 3DVAR (three-dimensional variational), these three methods were applied to the same CMAQ model simulations with identical experimental settings. This is the first attempt in the field of chemical DA to compare the EnKF and EnSRF methods. An identical set of surface fine particulate matter (PM_2.5) were assimilated every 6 h by all the DA methods over a CMAQ domain of East Asia, during the period from 01 May to 11 June 2016. In parallel with 'reanalysis experiments', we also carried out '48 h prediction experiments' using the optimized initial conditions produced by the three DA methods. Detailed analyses among the three DA methods were then carried out by comparing both the reanalysis and the prediction outputs with the observed surface PM_2.5 over four regions (i.e., South Korea, the Beijing-Tianjin-Hebei (BTH) region, Shandong province, and Liaoning province). The comparison results revealed that the EnSRF produced the best reanalysis and prediction fields in terms of several statistical metrics. For example, when the 3DVAR, EnKF, and EnSRF methods were used, averaged normalized mean biases (NMBs) decreased by (57.6, 85.6, and 91.8) % in reanalyses and (39.7, 87.6, and 91.5) % in first-day predictions, compared to the CMAQ control experiment (i.e., without DA) over South Korea, respectively. Also, over the three Chinese regions, the EnSRF method outperformed the EnKF and 3DVAR methods.

Vertical aerosol data assimilation technology and application based on satellite and ground lidar: A review and outlook

Observations and numerical models are mainly used to investigate the spatiotemporal distribution and vertical structure characteristics of aerosols to understand aerosol pollution and its effects. However, the limitations of observations and the uncertainties of numerical models bias aerosol calculations and predictions. Data assimilation combines observations and numerical models to improve the accuracy of the initial, analytical fields of models and promote the development of atmospheric aerosol pollution research. Numerous studies have been conducted to integrate multi-source data, such as aerosol optical depth and aerosol extinction coefficient profile, into various chemical transport models using various data assimilation algorithms and have achieved good assimilation results. The definition of data assimilation and the main algorithms will be briefly presented, and the progress of aerosol assimilation according to two types of aerosol data, namely, aerosol optical depth and extinction coefficient, will be presented. The application of vertical aerosol data assimilation, as well as the future trends and challenges of aerosol data assimilation, will be further analysed.

Seasonal Dependence of Aerosol Data Assimilation and Forecasting Using Satellite and Ground-Based Observations

This study examines the performance of a data assimilation and forecasting system that simultaneously assimilates satellite aerosol optical depth (AOD) and ground-based PM10 and PM2.5 observations into the Weather Research and Forecasting model coupled with Chemistry (WRF-Chem). The data assimilation case for the surface PM10 and PM2.5 concentrations exhibits a higher consistency with the observed data by showing more correlation coefficients than the no-assimilation case. The data assimilation also shows beneficial impacts on the PM10 and PM2.5 forecasts for South Korea for up to 24 h from the updated initial condition. This study also finds deficiencies in data assimilation and forecasts, as the model shows a pronounced seasonal dependence of forecasting accuracy, on which the seasonal changes in regional atmospheric circulation patterns have a significant impact. In spring, the forecast accuracy decreases due to large uncertainties in natural dust transport from the continent by north-westerlies, while the model performs reasonably well in terms of anthropogenic emission and transport in winter. When the south-westerlies prevail in summer, the forecast accuracy increases with the overall reduction in ambient concentration. The forecasts also show significant accuracy degradation as the lead time increases because of systematic model biases. A simple statistical correction that adjusts the mean and variance of the forecast outputs to resemble those in the observed distribution can maintain the forecast skill at a practically useful level for lead times of more than a day. For a categorical forecast, the skill score of the data assimilation run increased by up to 37% compared to that of the case with no assimilation, and the skill score was further improved by 10% through bias correction.

Length Scale Analyses of Background Error Covariances for EnKF and EnSRF Data Assimilation

Data assimilation (DA) combines incomplete background values obtained via chemical transport model predictions with observational information. Several 3-Dimensional variational (3DVAR) and sequential methods (e.g., ensemble Kalman filter (EnKF)) are used to define model errors and build a background error covariance (BEC) and are important factors affecting the prediction performance of DA. The BEC determines the spatial range, where observation concentration is reflected in the model when DA is applied to an air pollution transport model. However, studies investigating the characteristics of BEC using air quality models remain lacking. In this study, horizontal length scale (HLS) and vertical length scale (VLS) analyses of a BEC were applied to EnKF and ensemble square root filter (EnSRF), respectively, and two ensemble-based DA methods were performed; the characteristics were compared with those of a BEC applied to 3DVAR. The results of 6 h PM2.5 predictions performed for 42 days were evaluated for a control run without DA (CTR), 3DVAR, EnKF, and EnSRF. HLS and VLS respectively exhibited a high correlation with the ground wind speed and with the planetary boundary layer height for diurnal and daily variations; EnKF and EnSRF exhibited superior performances among all the methods. The root mean square errors were 11.9 μg m−3 and 11.7 μg m−3 for EnKF and EnSRF, respectively, while those for 3DVAR and CTR were 12.6 μg m−3 and 18.3 μg m−3, respectively. Thus, we proposed a simple method to find a Gaussian function that best described the error correlation of the BEC based on the physical distance.

The impacts of background error covariance on particulate matter assimilation and forecast: An ideal case study with a modal aerosol model over China

Accurate estimation of background error covariance (BEC) is the key to successful data assimilation (DA). In aerosol three-dimensional variational (3DVAR) DA, the National Meteorological Center (NMC) method is typically applied to estimate BEC, which uses the difference between forecasts of dissimilar lengths valid at the common time. The difference will be considerably small when the underestimation of the aerosol is caused by the lack of emissions or the missing of chemical progress, which makes the aerosol concentration field too difficult to constrain. In this study, a modified module for adjusting the BEC of individual aerosol species was developed in the Gridpoint Statistical Interpolation (GSI) 3DVAR system. This module was mainly utilized to expand the standard deviations and the horizontal length scales of BEC for the specified aerosol components by multiplying an adjustment factor. The results of the impacts of BEC on PM₁₀ 24-hour forecast indicated that the horizontal length scales take a relatively more important role than the standard deviations. The horizontal length scales affect the influence sphere more significantly, which might be crucial for the longer length forecast. Moreover, the larger and the wider differences of the aerosol initial conditions produced by DA, the longer duration of DA benefits. Using the original BEC, the 24-hour forecast of PM₁₀ reduced fractional error by 13%, while using the modified BEC in DA can decline fractional error by 29%. More work needs to be conducted to investigate how to modify the aerosol BEC in 3DVAR, or how to generate a suitable BEC, which is crucial for aerosol forecast and analysis, especially during the aerosol-polluted period.

Data Assimilation of AOD and Estimation of Surface Particulate Matters over the Arctic

In this study, more accurate information on the levels of aerosol optical depth (AOD) was calculated from the assimilation of the modeled AOD based on the optimal interpolation method. Additionally, more realistic levels of surface particulate matters over the Arctic were estimated using the assimilated AOD based on the linear relationship between the particulate matters and AODs. In comparison to the MODIS observation, the assimilated AOD was much improved compared with the modeled AOD (e.g., increase in correlation coefficients from −0.15–0.26 to 0.17–0.76 over the Arctic). The newly inferred monthly averages of PM10 and PM2.5 for April–September 2008 were 2.18–3.70 μg m−3 and 0.85–1.68 μg m−3 over the Arctic, respectively. These corresponded to an increase of 140–180%, compared with the modeled PMs. In comparison to in-situ observation, the inferred PMs showed better performances than those from the simulations, particularly at Hyytiala station. Therefore, combining the model simulation and data assimilation provided more accurate concentrations of AOD, PM10, and PM2.5 than those only calculated from the model simulations.