Super dust storms re-occurred over East Asia in 2021 spring and casted great health damages and property losses. It is essential to achieve an accurate dust forecast to reduce the damage for early warning. The forecasting system fundamentally relies on a numerical model which can forecast the full evolution of dust storms. However, large uncertainties exist in model forecasts. Meanwhile, various near-real-time observations are available that contain valuable dust information. A dust storm forecasting system is here developed through coupling a chemical transport model, LOTOS-EUROS, and Localized EnKF (LEnKF) assimilation approach. The assimilations are carried out via an interface of our self-designed assimilation toolbox, PyFilter v1.0. Ground-based PM₁₀ measurements from air quality monitoring network are assimilated. Sequential assimilation tests are carried out over the 2021 spring super dust storms. The results show that the assimilation-based forecasting system produces a promising dust forecast than model-only forecast, and the improvements is also validated through comparing to the independent MODIS aerosol optical depth (AOD). Superior performance is obtained when LEnKF is implemented, as the localization helps EnKF in resolving the PM₁₀ measurements that have a large spatial variability with limited ensemble members. In addition, sensitivity experiments are conducted to exploit the distance-dependent localization for the LEnKF. Considering both cases, the optimal choice of the distance is tested to be around 500 km: the larger distance is less effective in removing the spurious correction, while the smaller one easily falls into the local optimum and the model would become divergent rapidly.

@en

Atmospheric ammonia has been hazardous to the environment and human health for decades. Current inventories are usually constructed in a bottom-up manner and subject to uncertainties and incapable of reproducing the spatiotemporal characteristics of ammonia emission. Satellite measurements, for example, Infrared Atmospheric Sounder Interferometer (IASI) and Cross-Track Infrared Sounder, which provide global coverage of ammonia distribution, have gained popularity in ammonia emission estimation through data assimilation methods. However, satellite-based emission inversion studies on China are limited. In this study, we propose a four-dimensional ensemble variational-based ammonia emission inversion system to optimize ammonia emissions in China. It was developed by assimilating the IASI ammonia retrievals onboard Meteorological Operational satellite A and B into a chemical transport model Goddard Earth Observing System Chemical model (GEOS-Chem). Monthly inversion experiments were conducted in April, July, and October 2016 to test the performance. The inversion result indicated that the prior inventory from the MEIC model captured ammonia spreads in general; however, it heterogeneously underrated the emission intensity. The increments obtained in the assimilation were as high as 50% in North, East, and Northwest China. The posterior emission inventory presented a regional emission flux consistent with relevant studies. Driven by the optimized source estimate, GEOS-Chem provides superior results than using the prior in the evaluation of the assimilated IASI retrievals and the surface ammonia concentration measured by the ground-based Ammonia Monitoring Network in China.

@en

Last spring, super dust storms reappeared in East Asia after being absent for one and a half decades. The event caused enormous losses in both Mongolia and China. Accurate simulation of such super sandstorms is valuable for the quantification of health damage, aviation risks, and profound impacts on the Earth system, but also to reveal the climatic driving force and the process of desertification. However, accurate simulation of dust life cycles is challenging, mainly due to imperfect knowledge of emissions. In this study, the emissions that lead to the 2021 spring dust storms are estimated through assimilation of MODIS AOD and ground-based PM10 concentration data simultaneously. With this, the dust concentrations during these super storms could be reproduced and validated with concentration observations. The multi-observation assimilation is also compared against emission inversion that assimilates AOD or PM10 concentration measurements alone, and the added values are analyzed. The emission inversion results reveal that wind-blown dust emissions originated from both China and Mongolia during spring 2021. Specifically, 19.9×106 and 37.5×106ĝ€¯t of particles were released in the Chinese and Mongolian Gobi, respectively, during these severe dust events. By source apportionment it was revealed that the Mongolian Gobi poses more severe threats to the densely populated regions of the Fenwei Plain (FWP) and the North China Plain (NCP) located in northern China than does the Chinese Gobi. It was estimated that 63ĝ€¯% of the dust deposited in FWP was due to transnational transport from Mongolia. For NCP, the long-distance transport dust from Mongolia contributes about 69ĝ€¯% to the dust deposition.

@en

When calibrating simulations of dust clouds, both the intensity and the position are important. Intensity errors arise mainly from uncertain emission and sedimentation strengths, while position errors are attributed either to imperfect emission timing or to uncertainties in the transport. Though many studies have been conducted on the calibration or correction of dust simulations, most of these focus on intensity solely and leave the position errors mainly unchanged. In this paper, a grid-distorted data assimilation, which consists of an image-morphing method and an ensemble-based variational assimilation, is designed for realigning a simulated dust plume to correct the position error. This newly developed grid-distorted data assimilation has been applied to a dust storm event in May 2017 over East Asia. Results have been compared for three configurations: a traditional assimilation configuration that focuses solely on intensity correction, a grid-distorted data assimilation that focuses on position correction only and the hybrid assimilation that combines these two. For the evaluated case, the position misfit in the simulations is shown to be dominant in the results. The traditional emission inversion only slightly improves the dust simulation, while the grid-distorted data assimilation effectively improves the dust simulation and forecasting. The hybrid assimilation that corrects both position and intensity of the dust load provides the best initial condition for forecasting of dust concentrations.@en

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 <span classCombining double low line"cit"idCombining double low line"xref_paren.1">(<a hrefCombining double low line"#bib1.bibx31">Jin et al.</a>, <a hrefCombining double low line"#bib1.bibx31">2019</a><a hrefCombining double low line"#bib1.bibx31">b</a>)</span>. 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.

@en

Numerical models of chemical transport have been used to simulate the complex processes involved in the formation and transport of air pollutants. Although these models can predict the spatiotemporal variability of a variety of chemical species, the accuracy of these models is often limited. Therefore, in the past two decades, data assimilation methods have been applied to use the available measurements for improving the forecast. Nowadays, machine learning techniques provide new opportunities for improving the air quality forecast. A case study on PM ₁₀ concentrations during a dust storm is performed. It is known that the PM ₁₀ concentrations are caused by multiple emission sources, e.g., dust from desert and anthropogenic emissions. An accurate modeling of the PM ₁₀ concentration levels owing to the local anthropogenic emissions is essential for an adequate evaluation of the dust level. However, real-time measurement of local emissions is not possible, so no direct data is available. Actually, the lack of in-time emission inventories is one of the main reasons that current numerical chemical transport models cannot produce accurate anthropogenic PM ₁₀ simulations. Using machine learning techniques to generate local emissions based on real-time observations is a promising approach. We report how it can be combined with data assimilation to improve the accuracy of air quality forecast considerably.

@en

Severe dust storms present great threats to the environment, property and human health over the areas in the downwind of arid regions. Several dynamical dust models have been developed to predict the dust concentrations in the atmosphere. Currently, the accuracy of these models is limited mainly due to the imperfect modeling of dust emissions. Along with the progress in the dust and aerosol modeling, the advances in sensor technologies have made large-scale aerosol measurements feasible. The rich measurements provide opportunities to estimate uncertain emission fields, and subsequently, to improve the forecast skill. Such process of emission optimization conditioned on measurements is usually referred as emission inversion. Here, the termof emission inversion specially represents the way of deriving estimates from observations through the use of an atmospheric chemical transport model and a data assimilationmethod. @en

Aerosol optical depths (AODs) from the new Himawari-8 satellite instrument have been assimilated in a dust simulation model over East Asia. This advanced geostationary instrument is capable of monitoring the East Asian dust storms which usually have great spatial and temporal variability. The quality of the data has been verified through a comparison with AErosol RObotic NETwork AODs. This study focuses on extreme dust events only when dust aerosols are dominant; promising results are obtained in AOD assimilation experiments during a case in May 2017. The dust emission fields that drive the simulation model are strongly improved by the inverse modeling, and consequently, the simulated dust concentrations are in better agreements with the observed AOD as well as ground-based observations of PM ₁₀ . However, some satellite AODs show significant inconsistence with the simulations and the PM ₁₀ and AErosol RObotic NETwork observations, which might arise from retrieval errors over a partially clouded scene. The data assimilation procedure therefore includes a screening method to exclude these observations in order to avoid unrealistic results. A dust mask screening method is designed, which is based on selecting only those observations where the deterministic model produces a substantial amount of dust. This screen algorithm is tested to give more accurate result compared to the traditional method based on background covariance in the case study. Note that our screen method would exclude valuable information in case the model is not able to simulate the dust plume shape correctly; hence, applications in related studies require inspections of simulations and observations by user.

@en

Data assimilation algorithms rely on a basic assumption of an unbiased observation error. However, the presence of inconsistent measurements with nontrivial biases or inseparable baselines is unavoidable in practice. Assimilation analysis might diverge from reality since the data assimilation itself cannot distinguish whether the differences between model simulations and observations are due to the biased observations or model deficiencies. Unfortunately, modeling of observation biases or baselines which show strong spatiotemporal variability is a challenging task. In this study, we report how data-driven machine learning can be used to perform observation bias correction for data assimilation through a real application, which is the dust emission inversion using PM10 observations. PM10 observations are considered unbiased; however, a bias correction is necessary if they are used as a proxy for dust during dust storms since they actually represent a sum of dust particles and non-dust aerosols. Two observation bias correction methods have been designed in order to use PM10 measurements as proxy for the dust storm loads under severe dust conditions. The first one is the conventional chemistry transport model (CTM) that simulates life cycles of non-dust aerosols. The other one is the machine-learning model that describes the relations between the regular PM10 and other air quality measurements. The latter is trained by learning using 2 years of historical samples. The machine-learning-based non-dust model is shown to be in better agreement with observations compared to the CTM. The dust emission inversion tests have been performed, through assimilating either the raw measurements or the bias-corrected dust observations using either the CTM or machine-learning model. The emission field, surface dust concentration, and forecast skill are evaluated. The worst case is when we directly assimilate the original observations. The forecasts driven by the a posteriori emission in this case even result in larger errors than the reference prediction. This shows the necessities of bias correction in data assimilation. The best results are obtained when using the machine-learning model for bias correction, with the existing measurements used more precisely and the resulting forecasts close to reality.

@en

In previous studies, a number of model-based dust forecasts and early warning systems have been developed for the prevention of environmental impacts due to dusts. However, the accuracy of the model is limited by imperfect identification of dust emissions, in particular by the friction velocity threshold parameterization in the emission process. In this study, an integrated dust storm forecast system - LOTOS-EUROS/Dust coupled with reduced-tangent-linearization 4DVar data assimilation has been developed. In order to overcome the inflexibility and inaccuracy of the existing friction velocity threshold parameterization in large-scale models, a spatially varying multiplicative factor for the threshold is introduced. This parameter is estimated by assimilating measurements from a field station network developed by China Ministry of Environmental Protection. The data assimilation algorithm is adjoint-free, and its computational complexity increases with the number of uncertain parameters. Two model reducing techniques, sensitivity-based parameter filters and proper orthogonal decomposition, are sequentially implemented after each other, which lead to a reduction of parameter dimension from initially O(10⁴) to O(10²). Twin experiments are conducted to evaluate the impact of assimilation settings on the dust forecast accuracy. In addition, experiments with real observations are conducted. However, these observations also reflect the aerosol concentration from local emissions. To effectively use these observations as representative for dust concentrations, an observation bias correction and a variable representation error are designed. Improvements on the dust storm forecast with our system are demonstrated.

@en