An overview of mineral dust modeling over East Asia

East Asian dust (EAD) exerts considerable impacts on the energy balance and climate/climate change of the earth system through its influence on solar and terrestrial radiation, cloud properties, and precipitation efficiency. Providing an accurate description of the life cycle and climate effects of EAD is therefore critical to better understanding of climate change and socioeconomic development in East Asia and even worldwide. Dust modeling has undergone substantial development since the late 1990s, associated with improved understanding of the role of EAD in the earth system. Here, we review the achievements and progress made in recent decades in terms of dust modeling research, including dust emissions, long-range transport, radiative forcing (RF), and climate effects of dust particles over East Asia. Numerous efforts in dust/EAD modeling have been directed towards furnishing more sophisticated physical and chemical processes into the models on higher spatial resolutions. Meanwhile, more systematic observations and more advanced retrieval methods for instruments that address EAD related science issues have made it possible to evaluate model results and quantify the role of EAD in the earth system, and to further reduce the uncertainties in EAD simulations. Though much progress has been made, large discrepancies and knowledge gaps still exist among EAD simulations. The deficiencies and limitations that pertain to the performance of the EAD simulations referred to in the present study are also discussed.


Introduction
Dust, one of the major aerosol species contributing to global aerosol burden and optical depth, is a highly active component of the physical, chemical, and biogeochemical cycles of the earth system (Qian et al., 1999). As the second largest contributor of dust aerosols in the world, approximately 2000 Mt of desert dust from East Asia is injected into the atmosphere annually. Of this amount, approximately 30% is re-deposited onto the region's deserts, 20% is transported at the regional scale primarily within continental China, and the remaining 50% is transported farther eastwards to Korea, Japan, the Pacific islands, and across the Pacific Ocean to the United States, Canada, and even Greenland (e.g., Zhang et al., 1997;Liao and Seinfeld, 1998;Qian et al., 1999;Gong et al., 2003bGong et al., , 2006Huang et al., 2008aHuang et al., , 2012Huang et al., , 2016Uno et al., 2008;Eguchi et al., 2009;Fu et al., 2009;Liao et al., 2009;Shao et al., 2011;Li et al., 2011Li et al., , 2012Wang S. H. et al., 2012;. Such a large amount of East Asian dust (EAD) has a considerable impact on the regional climate by altering the radiative balance between incoming solar and outgo-ing terrestrial radiation in the atmosphere (the "direct effect") (Han et al., 2004(Han et al., , 2012Su et al., 2008;Han, 2010;Bi et al., 2011;Wang J. et al., 2012;Wang Z. L. et al., 2013). Moreover, EAD can modify the microphysical properties of clouds and the precipitation efficiency in the earth system (the "indirect effect" or "semi-direct effect") (Yin et al., 2002;Wang H. et al., 2004;Huang et al., 2006aHuang et al., , b, 2010Huang et al., , 2014Yin and Chen, 2007;Su et al., 2008;Qian et al., 2009;Chen et al., 2010;Jia et al., 2010;Guo and Yin, 2015). Dust particles also provide key mineral supplement to influence oceanic CO 2 uptake and nutrients in the ocean and rain forest via long-distance transport, and then modulate the ecosystem and carbon cycle at both regional and global scales (Shao et al., 2011). In addition, EAD induces highly adverse effects on social and economic development. For example, a severe dust storm in South Korea caused enormous economic loss and human suffering, costing approximately US $5600 million (0.8% of GDP). The total economic loss due to dust storms in Beijing ranged from 2268.5 million RMB (US $273.3 million) to 5796 million RMB (US $698.3 million), according to Ai and Polenske (2008). Increased concentrations of particulate matter in the atmosphere during dust storms can also lead to acute damage to respiratory systems and long-term damage such as desert pneumoconiosis Kameda et al., 2016).
This paper provides a review of the modeling framework for EAD-related physical processes and the associated radiative forcing (RF) and climate effects. A stateof-the-art model, the Weather Research and Forecasting model with Chemistry (WRF-Chem), is used to show the regional and seasonal variations of dust mass loading and its radiative effect over East Asia. Section 2 describes recent achievements in estimating dust emissions over East Asia. Section 3 summarizes the modeling of the longrange transport of EAD. Section 4 discusses the climate effects of EAD. Finally, several issues related to the achievements, uncertainties, and challenges involved in dust modeling are discussed in Section 5.

Dust emission schemes
Dust particles can be dispersed by turbulence into the free atmosphere, transported to remote regions by winds, and eventually deposited back onto the earth's surface. Dust emissions are also a key component of the dust cycle, which greatly influences other dust-related processes such as the energy and carbon cycle at the global scale in the earth system (Shao et al., 2011). Han et al. (2004), Zender et al. (2003, Shao and Dong (2006), and Shao et al. (2011) classified dust-emission schemes into three categories. In this paper, we focus on two major dust emission schemes in EAD simulations, hereafter referred to as Scheme I and Scheme II.
(1) Scheme I: The dust emission flux is a function of wind speed or friction speed with an empirical size distribution upon the emitted dust. This scheme was used by Tegen and Fung (1994), Mahowald et al. (1999), Ginoux et al. (2001), Perlwitz et al. (2001), Zender et al. (2003), and Bian et al. (2011). Here, we use the Georgia Tech/Goddard Global Ozone Chemistry Aerosol Radiation and Transport (GOCART) dust scheme (Ginoux et al., 2001) to illustrate this approach. The entire bare ground surface is assumed as a potential dust source in this scheme. The dust emission flux is expressed as a function of the surface wind speed and soil wetness. The GOCART scheme calculates the dust emission flux G (μg m -2 s -1 ) as where C (µg s 2 m -5 ) is an empirical proportionality constant (the default value of C is set to 1 µg s 2 m -5 ); S is a semi-empirical dust source function that defines potential dust source regions based on underlying surface factors such as the vegetation and snow cover; s p is the fraction of emissions for each dust-size class; u 10m is the horizontal wind speed at 10-m height; and u t is the threshold wind velocity for triggering dust emissions, which is a function of the particle size, air density, and surface soil moisture. Notably, the value of the empirical proportionality constant C is highly tunable because it depends on regionally specific data. Generally, we tune C to keep the modeled aerosol optical depth (AOD) consistent with the observational data, such as AERONET measurements or Multi-angle Imaging Spectroradiometer (MISR) and Moderate Resolution Imaging Spectroradiometer (MODIS) retrievals (Zhao et al., 2010(Zhao et al., , 2011Chen et al., 2013Chen et al., , 2014a. (2) Scheme II: This scheme is physically explicit and more complex. It uses sophisticated mechanisms of dust production to predict the size-resolved dust emissions (Shao, 2001(Shao, , 2004(Shao, , 2008Kok, 2011). Shao (2008) divided this dust emission scheme into three parts: saltation bombardment, aggregate disintegration, and aerodynamic entrainment. In recent years, dust emission schemes have included one or both of the first two mechanisms (Lu and Shao, 1999;Shao, 2004) because these two mechanisms are efficient and the most important (Klose et al., 2014). The dust emission scheme of Shao (2004) is calculated as follows: F where is the dust emission flux (μg m -2 s -1 ) for the particle group of size d i , generated by the saltation of particles of size d s ; c y is a dimensionless coefficient; η fi is the total dust fraction; γ is a weighting function; Q is the streamwise saltation flux of d s ; g is acceleration due to gravity; u * is the friction velocity; α m is the ratio of the mass of the impacting particle to the mass ejected by bombardment; ρ b is the ratio of free dust to aggregated dust; and P is the soil plastic pressure.
The dust emission from bin i is given bỹ where d 1 and d 2 define the lower and upper size limit, respectively, for saltating particles. The total dust emission rate is estimated as follows: Aerodynamic dust entrainment has not been included in traditional dust emission schemes because the intensity of aerodynamic dust emissions is much weaker than that of strong dust storms. However, this process may occur frequently, which may be a major contributor to regional or global dust budgets. Klose and Shao (2012) proposed the first comprehensive CTDE (Convective Turbulent Dust Emission) model and further estimated the contribution of convective dust emissions to the regional/global dust budget. Klose and Shao (2013) developed a large-eddy dust model (WRF-LES/D) with this new CTDE scheme to investigate the CTDE processes for various conditions of boundary-layer stability and wind fields for a given soil type. Based on Shao (2012, 2013), Klose et al. (2014) developed an improved parameterization scheme for CTDE, which provides a new description of instantaneous momentum flux and a correction function for the cohesive force to consider the effect of soil moisture, and the lifting force to consider the effect of vegetation roughness. The dust emission flux is given as where D is the viscous sub-layer thickness, w t is the particle terminal velocity, m p is the particle mass, T p is the particle response time, f is a given lifting force, f i is the cohesive force, and d is the particle diameter. As shown in Klose and Shao (2012), α N = N p D, N p is the particle number concentration. Thus, α N has the dimensions of [m -2 ] and can be interpreted as the over D integrated particle number concentration per unit area. Accordingly, α N decreases with particle size as follows: is an empirical constant of dimensions [m -2 ] and d ref is the particle size for which the scheme is calibrated.

AUGUST 2017
CTDE is initiated only if f exceeds the retarding force f t = f i + m ρ g, where g is the gravitational acceleration. For given distribution functions for the lifting force p(f) and inter-particle cohesive force p j (f i ) corresponding to particle size d j , the dust emission flux can be obtained as follows: where F j represents the dust emissions from bin j, is the PDF of , and is the PDF of f. The total dust emission flux is where d max is the maximum dust particle diameter (20 µm in this study), and P A (d j ) is the area's particle size distribution (PSD). Note that P A is used instead of mass PSD, P M , because the surface shear stress acts on the particle area instead of the particle volume. The P M (d j ) is usually directly available from particle size analysis, but for a given P M , P A can be estimated as

Simulation of dust emissions
Dust simulation research has shown that modeled vertically integrated parameters (e.g., AOD) are reproduced within a factor of 2 or 3, whereas the dust emissions are reproduced within a factor of 10, on the global scale (Zender et al., 2003;Todd et al., 2008;Huneeus et al., 2011). Quantification of dust emissions over East Asia is hampered by the uncertainties of dust production mechanisms, complex local geomorphological and topographical features, and various parameters of surface properties (e.g., soil composition, soil texture, soil wetness, vegetation, and land use) (Chao and Alexander, 1984;Zhao et al., 2006;Yumimoto et al., 2008;Bian et al., 2011).
To better understand and demonstrate dust emission properties over East Asia, a number of large field campaigns [e.g., APEX (Asian Atmospheric Particulate Environmental Change Studies), ACE-Asia (Asian Pacific Regional Aerosol Characterization Experiment), and ADEC (Aeolian Dust Experiment on Climate Impact) (Takemura et al., 2003;Mikami et al., 2006;Uno et al., 2006;Huang et al., 2008bHuang et al., , 2011Ge et al., 2011;Chen et al., 2013Chen et al., , 2014bZhao et al., 2013b), together with ad-vanced retrieval methods for remote sensing (Wang and Min, 2008;Wang and Huang, 2009) that address EAD related problems, have been conducted over the last few decades, providing a benchmark test dataset for the development of EAD parameterization schemes in numerical models. Gong et al. (2003a) incorporated a detailed soil texture map and an up-to-date distribution of desertification in China in the dust emission scheme. The observed soil grain-size distribution, water moisture, and meteorology could be used to optimize the appropriate parameters and conditions for Chinese soil dust emissions and transport. The simulated dust concentrations in source regions and areas downwind of eastern China to North America compare well with ground-based and aircraft measurements, as well as satellite retrievals, in East Asia and North America (Gong et al., 2003a;Zhao et al., 2003). Chen et al. (2014a) used the GOCART dust emission scheme coupled with the WRF-Chem model to evaluate and compare with observed EAD in 2007-11. A comparison of MISR retrievals and WRF-Chem model simulations over the Taklimakan Desert (TD) and Gobi Desert (GD) in 2007-11 is shown in Fig. 1. The results indicate that the GOCART dust emission scheme adequately captures the monthly variations of dust AOD in the dust source regions.
However, the uncertainty in the simulated EAD emission flux is still severe in numerical models. Uno et al. (2006) compared eight dust models coupled with eight dust emission schemes over East Asia in an intercomparison study. They noted that numerical modeling correctly predicted the onset and ending of EAD events, but large differences with respect to dust emissions were apparent. The wide range of estimated dust emission fluxes over East Asia may be partly due to the differences in dust emission schemes among these simulations, but we cannot ignore the difference in other respects, such as meteorological conditions, dust size parameterization, and land-cover data used. Therefore, Kang et al. (2013) performed a parallel comparison among three dust emission schemes Bergametti, 1995, Lu andShao, 1999 (saltating bombardment only), and Shao, 2004 (both saltating bombardment and aggregate disintegration)-hereinafter referred to as MB, LS, and S04, respectively] with the same WRF-Chem model configuration. The large differences in vertical dust fluxes among different emission schemes ranged from a factor of approximately 100 for clay to 10 for sand.
Overall, the uncertainties in dust emissions over East Asia derive mainly from the following two causes: (1) Dust emission mechanism. Dust from soil grains produced by wind erosion involves complex nonlinear processes governed by meteorology. Scheme I is relatively simple and easier to apply in practice. The drawback is its empirical or semi-empirical nature, which lacks a physical mechanism for dust emissions. Scheme II considers sophisticated wind-erosion physics and dust emission mechanisms, which are more accurate and complex than Scheme I. This scheme improves in understanding of the dust emission mechanism and enhances the accuracy of dust emission prediction. However, a large number of input parameters that are difficult to estimate in Scheme II also add to the uncertainties in the simulated dust emission fluxes. For the LS and SO4 schemes, the dust emission flux is most sensitive to the plastic pressure. The appropriate tuning of parameters in Scheme II becomes important in dust modeling.
(2) Soil/surface information. The spatial distributions of land-surface properties over East Asia strongly affect the intensity of EAD emissions. Uno et al. (2006) argued that accurate soil/surface conditions and land surface information are more important than the complexity of the dust emission scheme. Therefore, the lack of detail with respect to land surface conditions (e.g., soil texture, surface roughness length, moisture content, soil particle size distribution, vegetation cover, land cover type, and snow cover), and the dearth of observations over the dust source regions of East Asia, have been a major obstacle to reducing the uncertainties in simulating dust emission fluxes (Sokolik and Golitsyn, 1993;Gong et al., 2003b;Zhao et al., 2006). Moreover, different size distributions of emitted dust may also add uncertainties in dust simulation (Zhao et al., 2013a).

Dust transport
Dust particles, along their long-range transport, can significantly change RF, atmospheric chemistry, and biogeochemical cycles of the atmosphere, land, and ocean, even into remote regions (Wang et al., 2000Liu et al., 2011Liu et al., , 2010. Despite its climatic and geochemical importance, dust transport, based on both ground-based and satellite observations, is difficult to be captured sufficiently because it ranges from short durations at the regional scale to long durations at the global scale. Numerical modeling provides an alternative and systematic approach for investigating horizontal and vertical distributions as well as the seasonal and interannual variability of EAD in longrange transport.

Eastward and southeastward transport
The eastward/southeastward transport of EAD exerts significant impacts on the environment, human health, and climate/climate change over East Asia, given EAD's high association to anthropogenic emissions (Huang et al., 2006c;Li et al., 2011;Han et al., 2012;Hsu et al., 2013). The long-distance transport of EAD supplies mineral elements to the northern South China Sea, and thus remarkably enhances phytoplankton growth by enriching the chlorophyll-a concentration (Wang S. H. et al., 2012). In addition, EAD is also adsorbed by nitric acid, suppressing the cycle of nitrogen oxides and decreasing the Asian pollution effect on surface ozone in America (Fairlie et al., 2010). Moreover, EAD also has important implications for atmospheric pollution, since the acidic surface of dust can accelerate AUGUST 2017 the nitration reactions of polycyclic aromatic hydrocarbons (PAHs) in the atmosphere, and thus the nitrated polycyclic aromatic hydrocarbons (NPAHs) increase dramatically during dust storms. NPAHs can enhance the toxicity of EAD in urban environments, which has adverse effects on people's health, e.g., through carcinogenicity (Kameda et al., 2016).
Dust storms over East Asia occur mainly in spring, due to the considerable instability of synoptic weather systems during this season. Results show that cold-air outbreaks induced by Mongolian cyclonic and frontal systems lead to about 78% of all dust storms in China (Sun et al., 2001;Huang et al., 2008b). Dust particles can be uplifted to high levels by upward motion and turbulent mixing, and are then transported farther to Korea, Japan, and across the Pacific Ocean to the United States by westerly winds (Gong et al., 2003a;Huang et al., 2012;Wang J. et al., 2012;Chen et al., 2016). The pathways of cold-air outbreaks can be divided into three broad categories: north (32%), northwest (41%), and west (27%). Evidently, the northwest route is the predominant one (Sun et al., 2001).
In particular, a strong dust storm originating from the GD moved eastwards, southeastwards, and southwards, successively, in March 2010, even reaching and attacking the Pearl River Delta (such as Hong Kong) and the South China Sea Wang S. H. et al., 2012), which was a unique case compared with previous ones that mostly moved eastwards. This southward-moving dust storm attracts considerable interest from researchers, yielding many studies on its environment and climatic effects. According to Li et al. (2011), it was the special synoptic conditions, i.e., the high pressure at the surface shifting to the Yellow Sea and then mixing with low pressure at the surface, inducing northerly and northeasterly winds, thus resulting in dust particles being transported back to the continent from the sea and influencing the Pearl River Delta, which was unique compared with the dust storms in March 2002 and 2006. This southeastward transport of EAD was extremely severest and the PM 10 concentration was above 400 mg m -3 over western China and the GD, and below 100 mg m -3 in remote regions such as Japan. Besides, dust aerosol was a main contributor to the PM 10 concentration. Furthermore, dust deposition during the dust storm of March 2010 was disproportionate. The relatively finer size of dust particles with a low dry deposition ratio led to the increase in dust concentration. The dry deposition velocities at peak time were about 0.2-0.6 cm s -1 , far lower than usual (1-2 cm s -1 ), which indicates that, in some cases, the synoptic conditions of the atmosphere can be more important than the strength of the dust storm, in terms of the southeastward transport of EAD to the South China Sea (Hsu et al., 2013). In addition, Chen et al. (2016) also noted that the GD was located in the warm zone before the approach of a cold front from the northwest, where enhanced convection increased the momentum transfer in the middle and lower troposphere. Moreover, the GD is located in a relatively flat, high-altitude region influenced by the confluence of the northern and southern westerly jets. Therefore, GD dust transmission was the ubiquitous contribution over East Asia. Compared with the quantity of the TD dust contribution to the dust sink in March 2010 (1.1 ton day -1 ), that of the GD dust was significantly greater (1.4 ton day -1 ).
The vertical distribution of dust aerosols during transport is crucial for estimating dust radiative heating rates (Sokolik and Golitsyn, 1993;Sokolik and Toon, 1996;Kim et al., 2006Kim et al., , 2011. However, the vertical distribution of dust is difficult to measure. The development of lidar and passive infrared remote sensors, e.g., IASI (the Infrared Atmospheric Sounding Interferometer) and AIRS (Atmospheric Infrared Sounder), have facilitated the detection of dust vertical structure. In particular, the Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observations (CALIPSO) satellite, which was launched in April 2006 and has two-wavelength depolarization lidar CALIOP devices on board, enable the accurate determination of the three-dimensional (3D) transport of aerosols at the global scale, especially for aerosol vertical distributions (Huang et al., 2008b;Liu et al., 2008).
Many studies have attempted to describe the spatial distribution and vertical structure of EAD by combining dust modeling with satellite/ground-based lidar observations, which provides comprehensive and reliable analyses of the 3D structure of EAD and sheds new insights into long-term EAD transport. These combinations can not only overcome the observational shortcomings such as the limited spatial and temporal coverage, but also help dig out the mechanisms involved in the long-range transport of EAD. Uno et al. (2009) used the SPRINTARS model and CALIOP retrievals to investigate the transport of dust clouds generated in the TD. The simulated dust extinction coefficients agreed well with the CALIOP retrievals in the transport pathway. Yumimoto et al. (2009Yumimoto et al. ( , 2010 and Eguchi et al. (2009) also reproduced the 3D structure of intercontinental dust transport from the TD to the Pacific Ocean, North America, and the Atlantic Ocean, using similar tools. Furthermore, studies have focused on EAD transport by using a 4DVar (four-dimensional variational) data assimilation system of regional dust modeling (RAMS/CFORS-4DVar) . Compared with 3DVar, 4DVar could simultaneously add to the temporal variations of lidar observations. The vertical distribution of simulated dust extinction coefficients was comparable to that from the CALIPSO results. An increase in EAD emissions (approximately 31%) in the data assimilation was observed, especially over the Mongolian region Yumimoto et al., 2008).

Southward transport
As discussed above, the eastward transport of EAD is common and has been elaborated on by many studies. Huang et al. (2007) firstly observed the southward transport of TD dust to the northern slopes of the Tibetan Plateau (TP) during summer based on CALIOP retrievals (Fig. 2). Figure 2 shows the dust layer that appeared along the northern slopes of the TP and extended from a height of 5-9 km above mean sea level. This study enhances our knowledge of EAD transport and furthers our understanding of the impact of the TP with transported dust aloft on regional and global climates. Xia et al. (2008) pointed out that, compared with springtime, the correlation of the seasonal variation of the AOD over the TP and those over the TD was closer in summertime, based on MISR retrievals. Jia et al. (2015) also found that meteorological conditions and topography benefit the transport of dust emitted adjacent to the TP towards the plateau in summer. Liu et al. (2015) further investigated the transport of summer dust and anthropogenic aerosols over the TP by implementing the SPRINTARS model with a nonhydrostatic regional model.
The mechanism of this observed southward transport of the TD dust in summer was investigated with the WRF-Chem model by Chen et al. (2013). Figure 3 shows a schematic representation of the streamlines and PM 2.5 dust dry mass concentrations during 26-29 July 2006, to indicate the flow trajectories along the mountainous terrain. The dust originated from the TD and propagated towards the TP under the steering of easterly winds. The flow trajectories were close to the ground on 26-27 July (Figs. 3a, b), but were no longer close to the ground after this date, with weakened dust emissions and increased southward transport of dust over the following days (Figs. 3c, d). A schematic representation of the emission and transport processes of TD dust near the TP is also shown. Specifically, the strong dust storm over the TD was generated by a mesoscale cold front system. Subsequently, large amounts of dust particles along the Kunlun Mountains were lifted above the northern slopes of the TP due to topographic lift and horizontal shear. The movement and development of the TP's low vortex led to strong vertical motion over the TP, resulting from the impact of the thermal effects of the TP (i.e., acting as a "heat pump"). Therefore, dust particles were lifted to the upper troposphere and formed a thick dust layer approximately 6-8 km above mean sea level, caused by strong thermal mixing. Figure 4 shows the contributions of dust emissions, transport, and dry/wet deposition to the dust mass balance in 7 regions over East Asia in 2007-11 as simulated by the WRF-Chem. The positive (negative) sign represents an increase (a decrease) in the dust concentration. The percentage is the contribution to the dust  AUGUST 2017 mass balance in the region. Among the four budget terms, the emission contribution is absolute positive. Overall, the contribution from dust transport (~80%) to the TP dust concentration was much larger than that from dust emissions (~20%), especially during summer and autumn (Fig. 4). Dry deposition (-51%) and wet deposition (-49%) were comparable sinks of dust in the atmosphere over the TP. TP dust could produce a net warming effect with a domain average of +0.4 W m -2 and a maximum of +2.2 W m -2 . Dust resulted in a slight cooling effect (-0.7 W m -2 ) at the surface. It could modify the atmospheric heating profile over the TP and further affect the atmospheric circulation and energy budget at both regional and global scales (Chen et al., 2013).

Contributions from TD and GD dust
Over the East Asian region, TD and GD are the major dust source regions, from which dust particles are transported eastwards to eastern and southern China, Korea, Japan, and parts of the Pacific Ocean Uno et al., 2011;Huang et al., 2013). However, relative importance of the TD and GD to EAD remains unclear, with an ongoing debate as to whether the TD is the foremost source of dust over East Asia Zhang L. et al., 2010). Using observations to distinguish the contributions from TD and GD to the dust concentration over East Asia is difficult. Numerical models can be used to effectively investigate the emission, long-term transport, and wet/dry deposition of TD and GD dust (Chen et al., , 2017. Dust emission (+100%) is the only factor contributing positively to the dust mass balance over the TD and GD, while dust transport, as well as dry and wet deposition, are sinks of dust in the atmosphere. Specifically, GD dust has a larger influence on the dust mass loading than TD dust at the regional scale. The contribution of dust transport over the GD (-36%) is larger than that over the TD (-26%). The contributions of dry and wet deposition over the GD (-42% and -22%, respectively) are less than those over the TD (-49% and -25%, respectively) (Fig. 4). The tendency of GD dust transport is approximately -24.8 g m -2 each season, which is 35.1% of GD dust emission during the spring .
The differences in the topography, elevation, and dynamic and thermodynamic environmental conditions of the TD and GD lead to differences in the dust emission, lifting height of dust, and long-distance transport of dust. The GD is located at a confluence of the southern and northern branches of the upper-level jet stream. The downward momentum transport produced by a high- altitude westerly jet with suitable vertical circulation is the main cause of the higher wind speeds in the middle and lower troposphere over the GD, which leads to GD dust being more susceptible to uplift and long-range transport. Therefore, the GD plays a more important role than the TD in terms of dust transport over continental China. Compared with GD dust, TD dust is not easily transported out of the basin because of its lower elevation, complex topography, and low wind speeds in the upper atmosphere, although the TD has the largest dust emission flux in East Asia. However, TD dust is still the predominant contributor to the dust transported over long distances when uplifted to high altitudes (> 5000 m) (Sun et al., 2001;Shi et al., 2012;Chen et al., 2016).

Climatic effect of EAD
East Asia experiences high dust concentrations, which may result in substantial effects on the energy budget and climate change via the "direct" and "indirect" effects of dust (Wang et al., 2005;Huang et al., 2006c;Liao et al., 2009). However, the overall climatic effects of dust can vary greatly among regions because of the complex nonlinear interplay between dust aerosols and climate systems. The climatic effect of EAD has rarely been explored. The direct and indirect effects of EAD remain one of the largest uncertainties in climate models. In this section, we focus on the climate change induced by EAD in the atmosphere and in snow, and further investigate its influences on a series of important physical quantities, such as atmospheric stability, cloud lifetimes, and precipitation efficiency.

RF of dust over East Asia
Dust direct radiative forcing (DDRF) over East Asia accounts for about 42% of the total aerosol forcing, which plays an important role in the radiation budget and regional climate. Studies of DDRF over East Asia have been conducted since the 1990s. Numerical models have been used in combination with broadband radiation observations, sun photometer readings, and lidar measurements associated with satellite retrievals, to predict the spatiotemporal variability of DDRF over East Asia Zhang H. et al., 2009aZhang H. et al., , b, 2010Zhang H. et al., , 2012Wang W. C. et al., 2010). Generally, the uncertainties in DDRF over East Asia derive mainly from the following two causes. First, biases from dust emission, chemical transformation, dust transport, and removal processes of dust, greatly affect the precision of DDRF assessments. Second, differences in the particle size distribution, vertical distribution of the dust layer, absorpt-ive characteristics, chemical composition, particle shape, and environmental conditions (e.g., surface albedo, solar zenith angle, and cloud height) can lead to large heterogeneities in the quantitative assessment of DDRF (Sokolik and Toon, 1996;Tegen et al., 1996;Takemura et al., 2003;Wang et al., 2004Wang et al., , 2007Wu et al., 2004). Therefore, more observational data of dust optical properties over East Asia were used to develop dust models. Takemura et al. (2003) coupled the SPRINTARS model with a 3-D aerosol transport-radiation model to simulate the DDRF over East Asia, based on observations from the APEC and ACE-Asia field campaigns (see Table 1). The imaginary part of the refractive index of soil dust was modified from 0.08 to 0.002 at the visible wavelength in the model, which agreed closely with the simulation and observation of single scattering albedo in the dust layer. Moreover, Wang et al. (2004) noted that EAD particles absorb less solar radiation than Saharan dust particles. Furthermore, the authors included the dust refractive index and dust size distribution from Cele (37°01′N, 80°44′E) in the TD in the radiative transfer model (Shi, 1984) and calculated the average DDRF over East Asia under clear-sky conditions (see Table 1). They further calculated the radiative heating rate of EAD from the deserts of China to the North Pacific Ocean using the same radiative transfer model (Wang et al., 2007). Huang et al. (2009) also used the Fu-Liou radiation model in combination with various surface observation data to investigate the EAD heating rate (Table 1).
In addition, separation of chemistry and meteorology increases the uncertainty in dust simulations (Grell et al., 2005). In recent years, increasing attention has been paid to online dust modeling, in which the feedbacks between dust and radiation were examined Wu et al., 2004;Zhao et al., 2008;Grell and Baklanov, 2011). More studies focusing on dust radiative and climate effects over East Asia have been conducted using two-way online coupled regional dust modeling Zhao et al., 2011;Chen et al., 2013Chen et al., , 2014a. Han et al. (2010) incorporated aerosol and gas chemistry modules into RIEMS (Regional Integrated Environmental Model System), based on MM5 (version 5 of the NCAR's Mesoscale Model), and developed a new online coupled regional climatechemistry-aerosol model (RIEMS-Chemaero). The dust emission flux was calculated for each size bin, while dust deflation was treated empirically based on several factors (e.g., threshold friction velocity, relative humidity, and the reduction factor of vegetation). This model was firstly applied to investigate the dust RF and transport of an extremely intense dust storm over East Asia (Han et al., 2012(Han et al., , 2013.  also developed a new radiative scheme within a two-way dust modeling procedure (integrating dust prediction and radiative effects) to investigate the DDRF over East Asia. The updated data on land desertification, land surface, soil texture, and dust optical properties over China's mainland were used. Chen et al. (2014b) applied the Rapid Radiative Transfer Model (RRTMG) to both shortwave (SW) and longwave (LW) radiation in WRF-Chem to calculate the DDRF over East Asia (Fig. 5). The authors believed that the online estimation of direct dust effects im-proved the prediction of meteorological variables (e.g., temperature, pressure, and wind) for the dust events. Moreover, these studies included LW radiation in the radiative scheme, which significantly decreased the biases in radiative budget estimation.

Climatic impact induced by direct and indirect/semidirect effect of EAD in the atmosphere
In addition to the direct effect of EAD on SW and LW radiation, the indirect effect of dust can modify the microphysical properties of clouds (e.g., cloud liquid water  2 February-March-April-May (FMAM); 3 Regional Climate Model Version 3 (RegCM3); 4 RIEMS-Chemaero represents an online-coupled regional climate-chemistry-aerosol model; 5 Santa Barbara DISORT Atmospheric Radiative Transfer model (SBDART model); 6 Regional Climate Model Version 4 (RegCM4).
content, cloud fraction, cloud particle size), as dust particles are efficient cloud condensation nuclei (CCN) and ice nuclei (IN), and further influence the precipitation efficiency. Yin et al. (2002) performed numerical simulations with a 2D nonhydrostatic cloud model with detailed microphysics to investigate the impacts of dust particles on development of cloud and precipitation. Their results indicated that dust particles can enhance the formation of rain droplets in continental clouds, while they have a small impact on precipitation in maritime clouds. The semi-direct effect involves dust particles absorbing solar radiation, which can be re-emitted as thermal radiation to heat the atmosphere, thus enhancing the evaporation of cloud droplets (Han et al., 2013). Using a dynamic cloud model with detailed microphysics of both warm-and ice-phase processes, Yin and Chen (2007) found that the dust layer at the cloud-base height and below 3 km could induce heating, suppress the formation of cloud droplets, and reduce cloud optical depth. In addition, dust particles within or above the clouds are likely to evaporate cloud droplets over the arid and semi-arid regions (Huang et al., 2006bWang W. C. et al., 2010;Zhao S. et al., 2015).
Combined with dust direct and indirect/semi-direct effects, many studies on the climate change induced by EAD have been conducted. Results show clear decreases in surface temperature accompanied by reduced sensible heat flux due to the scattering and absorption of solar radiation by dust particles (Wu et al., 2010;Wei and Zhang, 2011;Zhao S. et al., 2014). Moreover, the diurnal temperature range (DTR), as an important indicator of the climatic system, can reflect climatic anomalies well. Can mineral dust significantly influence the DTR over East Asia? How do dust aerosols contribute to the DTR variations? To address these questions, we employed the WRF-Chem model. Our results (Fig. 6) showed that dust aerosol could play a role in changing the DTR through modification of the surface energy balance, directly or through its impact on cloud properties, soil moisture, and snow cover over East Asia. Changes in DTR over East Asia induced by dust aerosol result from the diurnal asymmetry of surface temperature; that is, the maximum temperature (T max ) changes obviously while the minimum temperature (T min ) remains constant or is modified slightly (Fig. 6). The annual mean T max and T min decrease by 0.54 and 0.23°C, respectively, and the DTR decreases by 0.31°C, with the pattern similar to that of T max . The largest reduction in DTR, by up to 0.37°C, happened in springtime, which is the main period of dust production over East Asia.
In the daytime, the dust layers in the atmosphere can strongly absorb and scatter SW radiation, which reduces the amount of SW radiation reaching the ground. Dust  direct RF is negative because the SW cooling effect of dust outweighs the LW warming effect. The T max is dominated by surface SW solar radiation induced by dust. At nighttime, the SW RF of dust turns zero. The dust layers in the atmosphere can increase the LW radiation reaching the surface and induce a slight warming effect on surface temperature. However, the indirect effect of dust reduces the surface temperature through changing cloud fractions. Overall, dust reduces the surface temperature at nighttime, since the cooling effect caused by the dust indirect effect outweighs the LW warming effect of dust.
The cross-correlation coefficients of the changes in cloud water content, snow cover, soil moisture, precipitation, daily maximum temperature, daily minimum temperature, and DTR induced by EAD were calculated in this study. The correlation coefficient between the DTR and daily maximum temperature reached 0.68, while it was 0.15 between DDRF and daily minimum temperature. The variation in cloud water content and snow cover was closely related to the DTR, and the correlation coefficient reached -0.70 and -0.58, respectively, which indicated that the dust could indirectly affected the ground DTR by changing the cloud cover and snow cover over East Asia. EAD can accelerate the melting of snow on the TP and parts of North China, and reduce the snow cover in the area due to the dust absorption. The melting of snow weakened the surface albedo, and thus increased the DTR in the TP and parts of North China.
Precipitation responses induced by dust particles also vary in different areas and in different seasons. Guo and Yin (2015) investigated the impacts of EAD on the precipitation, summer monsoon, and sea surface temperature over East Asia using a regional coupled atmosphereocean-land model, i.e., the Regional Integrated Environment Model System (RCM RIEMS 2.0). The indication was that dust-induced atmospheric heating and surface cooling reduced precipitation by 0.03 mm day -1 in spring and increased precipitation by 0.28 mm day -1 in summer over East Asia. The simulation results also indicated that dust aerosol could weaken the East Asia summer monsoon by reducing the land-sea temperature contrast. Wu et al. (2010) pointed out that EAD helps to partly offset the "southern flood and northern drought" pattern in China, if only considering the direct effects of EAD. Sun et al. (2012) also found that the direct effect of EAD might lead to suppressed precipitation in the downwind areas (Northwest China and East China Sea) and enhance precipitation around its source. Moreover, as revealed in these studies, the "elevated heat pump" effect generated by the TP dust could alter moisture, cloud, and deep convection in northern India, which accelerated the snowmelt in the Himalayas, hence resulting in an earlier onset and intensification of the Indian summer monsoon. The increased rainfall in India, combined with the largescale pressure anomaly pattern at sea level, led to a northwestward shift of the rainbelt over East Asia and adjacent regions (Lau et al., 2006;2010;Huang et al., 2007;Qian et al., 2009Qian et al., , 2015. More recently, Gu et al. (2016) investigated the precipitation responses induced by dust in North Africa, South Asia, and East Asia using an atmospheric general circulation model. The direct and semi-direct effects of dust were considered in the simulation. Dust over the Saharan Desert enhances the upward motion and strengthens the cyclonic circulation to its south where precipitation occurs, leading to decreased precipitation over the Sahel regions and the tropical North Atlantic. The thermal and dynamic effects of East/South Asian dust correspond with North Africa. However, it differs in that dust particles are located in the upper troposphere over the major rainfall areas in South/East Asia, except southeastern China. The increasing upward motion and strengthened cyclonic circulation induced by dust aerosols lead to increased precipitation intensity and variability. As for southeastern China, where the monsoon regions are located, precipitation decreases during the premonsoon season because the dust layer is located to the north of it, which induces anomalously weakened upward motion of the atmosphere over southeastern China. The precipitation and cloud in southeastern China increase during the monsoon season because the stronger monsoon flow induced by dust aerosols bring moister and warmer air, although there is no obvious change in dust concentrations over these regions (Wu et al., 2013;Gu et al., 2016) (Fig. 7).

Climatic impact of EAD in seasonal snow
Snow is the most reflective natural surface on the earth. Small changes in snow albedo can greatly affect surface warming due to rapid feedbacks from changes in snow morphology, sublimation, and melting rates. Results have shown that the snow albedo, particularly within the visible bands, can be influenced by light-absorbing impurities deposited on or in the snowpack (Flanner et al., 2007(Flanner et al., , 2012Wang et al., 2014). In recent years, a number of modeling studies have indicated that the RF from black carbon (BC) in snow has been an important anthropogenic forcing for climate change during the past century (Lau et al., 2008(Lau et al., , 2010Ming et al., 2008;Qian et al., 2009Qian et al., , 2015Xu et al., 2009Xu et al., , 2012Cong et al., 2013). These studies illustrate that dust in snow causes a remarkable decline in its albedo, especially within the visible bands, which could have a significant effect on the atmospheric general circulation and hydrological cycle (Qian et al., 2009(Qian et al., , 2015Liu et al., 2013;Kuchiki et al., 2015). Kuchiki et al. (2015) further observed that the maximum dust mass concentration in surface snow in Japan was 260 ppmw, much larger than that of elemental carbon and organic carbon (OC) mass concentrations because (1) dust mass loading in the atmosphere is much larger than BC and OC, and (2) dust particles are more sensitive to dry deposition and melt amplification. Therefore, the radiative effect of dust in snow cannot be ignored in dust modeling.
The Snow, Ice, and Aerosol Radiation (SNICAR) model was developed by Flanner et al. (2007) to investigate the impact of light-absorbing aerosols on snow properties. This model applies the two-stream multilayer radiative approximation of Toon et al. (1989) for the radiative transfer solution, and theory from Warren and Wiscombe (1980) for estimating the albedo of pure or contaminated snow. Snow within the model is regarded as ice spheres following a log-normal distribution. The optical properties of aerosols and ice crystals are derived with Mie solutions. Dust particles are divided into four size bins to estimate the optical properties based on the Maxwell-Garnett approximation, with an assumption for ingredients of dust particles. This approach helps to better understand the dust RF in snow, one of the largest sources of uncertainty in the assessment of aerosol RF in climate change (Hansen et al., 1997;Aoki et al., 2006;Painter et al., 2007).  Dust in snow in different regions produces distinct radiative feedback. Therefore, observational studies of dust in snow over East Asia are necessary. Observational datasets of dust concentrations in midlatitude seasonal snow have been scarce. Many projects have been conducted to explore light-absorbing impurities in snow . Five field campaigns were conducted in 7 provinces in China from 2010 to 2014. These studies provided valuable observations of the albedo of seasonal snow in China. The results of these campaigns were reported by Huang et al. (2011) and . Based on these observational results, Zhao C. et al. (2014) coupled the SNICAR model (Flanner et al., 2007(Flanner et al., , 2012 with WRF-Chem, including sophisticated representation of snow metamorphism processes to simulate the RF of dust in seasonal snow over northern China. The study evaluated the simulated EAD in snow at a relatively high spatial resolution against field campaigns. The results showed that dust and BC mass concentrations in snow induced a similar magnitude of radiative warming (approximately 10 W m -2 ) in the snowpack, comparable to the magnitude of radiative cooling at the surface induced by absorbed aerosols in the atmosphere (Fig. 8).
Dust in snow over the TP has attracted considerable interest because its role in climate change can hardly be overstated. The averaged RF of dust in snow over the TP can increase from 1.1 to 8.6 W m -2 , exceeding the RF of BC in snow, during the melting seasons in the Zhadang Glacier. This can reduce the albedo of glaciers (also called the "dirtying" or "darkening" effect), exacerbate the mass loss of glaciers in the inner TP (Qu et al., 2014), and substantially change the surface albedo and radiative fluxes over East Asia (Qian et al., 2009). However, how do TP dust aerosols in snow impact the hydrological cycle on the TP, and even the monsoon climate over East Asia? Additional studies are needed to examine the impacts of TP dust, through direct and indirect effects, on snow, and hence its impacts on climate and the hydrological cycle.

Uncertainties in EAD simulations
East Asia is one of the largest dust source regions in the world. EAD modulates the energy balance and climate feedback through its influence on solar and terrestrial radiation, cloud properties, and precipitation on both regional and global scales. Detailed investigations of the life cycle and associated climate impacts of EAD are imperative because of its environmental, climatic, and geochemical importance. In this review, we summar- ized the achievements and progress in studies of the life cycle and associated DDRF and climate impacts of EAD over East Asia using numerical modeling. However, knowledge of the EAD's dust cycle and interactions with other biogeochemical cycles is far less complete than that for other dust source regions (e.g., the Sahara Desert). Many uncertainties and challenges in EAD simulations remain. Some crucial feedback mechanisms by which EAD modulates the cloud properties, cloud lifetimes, and rain formation are still under debate and highly uncertain.
Compared with other regions such as the Sahara Desert and North American deserts, accurate EAD simulations with a regional or global model still faces great challenges. Uncertainties in EAD modeling exist in various aspects as follow: (1) lack of systematic observation of dust aerosols from ground-based networks, aircraft observations, and satellite retrievals; (2) inaccurate description of dust emission, long-term transport, and dry/wet deposition in numerical models; and (3) poor understanding of the aggregation, coagulation, and heterogeneous chemical reaction of dust particles.
Moreover, large uncertainties also exist in physical parameterizations, including those for emission, cloud microphysics, planetary boundary layer, and convection (e.g., Yang et al., 2012Yang et al., , 2013Zhao et al., 2013a;Yan et al., 2014). Zou et al. (2014) and Yang et al. (2015a, b) showed that parametric uncertainties in cloud and precipitation processes in both regional and global climate models can evidently affect the simulations of cloud, precipitation, radiation, and circulation over the East Asian summer monsoon region. These could further influence the dust climatic effects because of the uncertainties associated with the simulated transport, radiative, and microphysical effects of dust. Therefore, the process by which dust acts as CCN or IN remains one of the most poorly understood issues for establishing a clear causal relationship between dust and precipitation. Substantial knowledge gaps in modeling the impacts of dust on clouds and precipitation over East Asia remain.

Future
EAD emission, transport, and radiative effects have received deserved attention in studies of radiative budget and climate change of the earth system. However, we also want to point out several directions in future research.
(1) The development of integrated systematic observations of EAD, especially the dust emission flux, dust dry/wet deposition, vertical structure of dust layers, aerosol size distribution, and dust mass concentrations in snow, is important for improving the dust modeling.
More specifically, powerful satellite retrievals (e.g., CALIOP), advanced instruments (e.g., ground-based lidar networks), and in situ measurements [e.g., the Aerosol Robotic Network (AERONET) and SACOL] are needed over East Asia to evaluate dust modeling results. Currently, barriers exist between observations and numerical models because of differences in research methods and technical backgrounds. Modeling studies have always passively accepted observed results but seldom considered observation errors and inverse problems of observations because of incomprehension regarding the principles and technical performances of instruments and observation methods. Observational studies have focused on revealing these phenomena but have often ignored some key parameters of dust modeling. Because of these inhibitors, the joint development of observational techniques and dust modeling is critical to EAD parameterizations (e.g., dust emission, dust deposition, radiation schemes, and land surface processes).
(2) Dust aerosols have various sources from both natural and anthropogenic processes (Tegen and Fung, 1994;Zender et al., 2003). All current studies of dust emission parameterization schemes have mainly focused on natural dust. Anthropogenic dust originates from human activities through modifying or disturbing soil, such as by agriculture, livestock herding, construction, and off-road vehicle use. Huang et al. (2015) found that anthropogenic dust accounts for more than 91.8%, 76.1%, and 73.9% of the total dust emission in East China, India, and North America, respectively (Fig. 9). Natural or anthropogenic dust plays an important role in the radiative budget and hydrological cycle. Therefore, developing an emission scheme for anthropogenic dust aerosols in numerical models is necessary for objective assessment of  (3) Huge uncertainties exist in the estimation of precipitation responses induced by EAD in model simulations. An increasing number of studies have noted that precipitation change is induced by dust from two aspects: surface cooling, which may reduce rainfall, and the heating of the atmosphere, which may enhance rainfall. Dust induces the redistribution of precipitation in different regions and different periods, rather than simply acts to increase or decrease precipitation. Most previous studies have investigated the climatic effect of dust by only considering dust direct effects. However, dust can also affect cloud properties and then climate through acting as IN in heterogeneous ice nucleation processes. The effect of EAD IN needs to be investigated as well.