WRF-Chem v3.9 simulations of the East Asian dust storm in May 2017: modeling sensitivities to dust emission and dry deposition schemes
- 1Institute for Climate and Global Change Research and School of Atmospheric Sciences, Nanjing University, 210023 Nanjing, China
- 2Joint International Research Laboratory of Atmospheric and Earth System Sciences & Institute for Climate and Global Change Research, Nanjing University, Nanjing, China
- 3School of Earth and Space Sciences, University of Science and Technology of China, Hefei, 230026, China
- 4Key Laboratory for Semi-Arid Climate Change of the Ministry of Education, College of Atmospheric Sciences, Lanzhou University, Lanzhou, 730000, China
- 5Key Laboratory of Marine Environment and Ecology, Ministry of Education/Institute for Advanced Ocean Study, Ocean University of China, Qingdao, 266100, China
Correspondence: Minghuai Wang (firstname.lastname@example.org)
Dust aerosol plays an important role in the radiative budget and hydrological cycle, but large uncertainties remain for simulating dust emission and dry deposition processes in models. In this study, we investigated dust simulation sensitivity to two dust emission schemes and three dry deposition schemes for a severe dust storm during May 2017 over East Asia using the Weather Research and Forecasting model coupled with chemistry (WRF-Chem). Results showed that simulated dust loading is very sensitive to different dry deposition schemes, with the relative difference in dust loading using different dry deposition schemes ranging from 20 %–116 %. Two dust emission schemes are found to produce significantly different spatial distributions of dust loading. The difference in dry deposition velocity in different dry deposition schemes comes from the parameterization of collection efficiency from impaction and rebound effect. An optimal combination of dry deposition scheme and dust emission scheme has been identified to best simulate the dust storm in comparison with observation. The optimal dry deposition scheme accounts for the rebound effect and its collection efficiency from impaction changes with the land use categories and therefore has a better physical treatment of dry deposition velocity. Our results highlight the importance of dry deposition schemes for dust simulation.
Dust aerosol is an important component in the atmosphere, and it can impact many processes of the Earth system. Through absorbing and scattering shortwave and longwave radiative fluxes, dust can alter the radiative budgets, which is called the direct effect (Chen et al., 2013; Kok et al., 2017; Zhao et al., 2010, 2011, 2012). Acting as cloud condensation nuclei (CCN) and ice nuclei (IN), dust can change cloud properties and precipitation, which are called the indirect effects (Creamean et al., 2013; Demott et al., 2010). Besides, dust aerosol can absorb solar radiation and change the atmospheric stability and therefore cloud formation, which is known as the semidirect effect (Hansen et al., 1997). Furthermore, natural dust is important for air quality assessments and has significant impacts on human health (Abuduwaili et al., 2010; Chen et al., 2019; Hofer et al., 2017; Jiménez-Guerrero et al., 2008; Ozer et al., 2007). Although great progress has been made in dust models and dust simulations in recent decades, large uncertainties remain in dust simulations (Huneeus et al., 2011; Prospero et al., 2010; Todd et al., 2008; Uno et al., 2006; Zender et al., 2004; Zhao et al., 2013).
A complete description of dust events includes dust emission, deposition, and transport processes. The differences in dust simulation mainly result from the uncertainties in dust emission, deposition, and transport processes in models. One uncertainty is from dry deposition processes. Dry deposition refers to the transport of particles from the atmosphere to the Earth's surface in the absence of precipitation (Seinfeld and Pandis, 2006). In most aerosol modeling, dry deposition velocity Vd is used to calculate the dry deposition flux, and Vd is usually modeled using the resistance-based approach (Pryor et al., 2008). In the resistance-based approach, Vd is determined by gravitational settling, aerodynamic resistance, and surface resistance. Surface resistance is determined by collection efficiency from Brownian diffusion, impaction, and interception and is corrected for particle rebound. Slinn (1982) proposed a semi-analytical description of particle collection efficiencies based on the wind tunnel studies, and many dry deposition schemes since then are variants of this model (Binkowski and Shankar, 1995; Giorgi, 1986; Peters and Eiden, 1992; Zhang et al., 2001). As the formulations for collection efficiencies from different dry deposition schemes are derived from measurements that have been obtained under different meteorological conditions and land surface types, there remains a large discrepancy of these formulations between different dry deposition models (Petroff et al., 2008).
At present, the comparisons of different dry deposition schemes with reliable field measurements are mainly focused on 1D dry deposition models (Hicks et al., 2016; Khan and Perlinger, 2017; Petroff et al., 2008; Ruijrok et al., 1995). For example, Hicks et al. (2016) compared five deposition models with observations and found that Vd predicted for particles less than 0.2 µm is consistent with the measurements, but predicted Vd can vary greatly in the size range of 0.3 to about 5 µm. However, few studies have been conducted to study how different dry deposition schemes affect aerosol concentrations and their spatial distribution in the 3D numerical models. Wu et al. (2018) compared the effects of different dry deposition schemes on black carbon simulation in a global climate model (CESM-CAM5) but did not examine how different dry deposition schemes affect aerosol concentrations for large-size aerosol particles (e.g., diameters>2.5 µm), such as dust.
Another uncertainty in dust simulation is the treatment of the dust emission process in models. Natural dust is typically emitted from dry, erodible surfaces when the wind speed is high. Dust emission process is closely related to soil texture, soil moisture content, surface conditions, atmospheric stability, and the wind velocity (Marticorena and Bergametti, 1995). Dust emission schemes are used to predict the dust emission flux and to describe the dust size distribution. Many studies have compared and evaluated the performance of different dust emission schemes (Kang et al., 2011; LeGrand et al., 2019; Su and Fung, 2015; Wu and Lin, 2013, 2014; Yuan et al., 2019; Zhao et al., 2010, 2006). These studies show a large diversity of simulated dust emission flux among different dust emission schemes. Zhao et al. (2006) implemented two dust emission schemes in the NARCM (Northern Aerosol Regional Climate Model) regional model and found that both schemes captured the dust mobilization episodes and produced the similar spatial distributions of dust loading over East Asia, but significant differences exist in the dust emission fluxes and surface concentrations. Kang et al. (2011) compared three dust emission schemes in WRF-Chem and found that the difference between the vertical dust fluxes derived from the three emission schemes can reach to several orders of magnitude. Yuan et al. (2019) found that one scheme strongly underestimated the dust emission, while another two schemes can better show the spatial and temporal variation in dust aerosol optical depth (AOD) based on WRF-Chem simulation of a storm outbreak in Central Asia. In another WRF-Chem study, Chen et al. (2017a) concluded that the dust emission differences mainly come from the dust emission flux parameterizations and differences in soil and surface input parameters in different dust emission schemes.
While dust emission schemes have been studied quite extensively, few studies have examined dust emission and dry deposition schemes simultaneously. As both dust emission schemes and dry deposition schemes contribute significantly to the uncertainties in dust simulations, evaluating dust emission schemes based on a single dry deposition scheme may be problematic, especially if the dry deposition schemes employed have deficiencies. For example, as a widely used regional model that has been coupled with a variety of dust emission schemes, the WRF-Chem model has been used in many studies to evaluate the performance of the dust emission schemes (LeGrand et al., 2019; Su and Fung, 2015; Wu and Lin, 2013, 2014; Yuan et al., 2019). But most of these studies use the GOCART aerosol scheme, and only one dry deposition scheme (Wesely et al., 1985) is coupled within the GOCART aerosol scheme. Zhang et al. (2019) compared the modeled dust deposition using the GOCART aerosol scheme in WRF-Chem with observed dust deposition and found that modeled dust deposition is highly underestimated by more than 1 order of magnitude compared to the observed deposition. This indicates that the dry deposition scheme (Wesely et al., 1985) in the GOCART aerosol scheme may not be suitable for dust simulation and needs to be further improved.
In this study, we adopted the MOSAIC aerosol scheme coupled within the WRF-Chem model to study how dry deposition schemes and dust emission schemes affect dust simulations by evaluating model results against observations. As the MOSAIC aerosol scheme includes several different dry deposition schemes, this allows us to choose more advanced dry deposition schemes. As the default MOSAIC aerosol scheme only includes the GOCART dust emission scheme, we further implemented the dust emission scheme Shao2011 (Shao et al., 2011) in the MOSAIC aerosol scheme, which allows us to compare these two widely used dust schemes along with multiple dry deposition schemes. The goals of this study are (1) to study dust simulation sensitivity to different dust emission schemes and dry deposition schemes and (2) to explore which combination of dust emission scheme and dry deposition scheme can better simulate dust storms in East Asia. The paper is organized as follows. Section 2 introduces the WRF-Chem model, dust emission schemes and dry deposition schemes used, experiments design and measurements. Section 3 analyzes the dust simulation sensitivity to dust emission schemes and dry deposition schemes and the comparisons with observations. Section 4 is the summary and discussion.
2.1 Model description
In this study, WRF-Chem version 3.9 is used. WRF-Chem is built based on the regional mesoscale model WRF and fully coupled with gas and aerosol chemistry module (Grell et al., 2005). A summary of the settings used to configure the model is listed in Table 1. The Noah land surface model (Chen and Dudhia, 2001) and the Yonsei University (YSU) planetary boundary scheme (Hong et al., 2006) are used in this study. The global soil categorization data set from the United States Geological Survey (USGS) with 24 land categories are used. The Rapid Radiative Transfer Model for General Circulation (RRTMG) radiation scheme (Iacono et al., 2008) is used to calculate the longwave and shortwave radiation. The Grell–Freitas convective scheme (Grell and Freitas, 2014) and the Morrison two-moment microphysics scheme (Morrison et al., 2008) are used. The gas-phase chemistry module used is the Carbon Bond Mechanism version Z (CBMZ; Zaveri and Peters, 1999). The aerosol module used here is the Model for Simulating Aerosol Interactions and Chemistry with four bins (MOSAIC 4-bin) (Zaveri et al., 2008). The MOSAIC 4-bin aerosol scheme divides airborne particles into four size bins by their effective diameter (0.039–0.156, 0.156–0.625, 0.625–2.5, and 2.5–10.0 µm) to represent aerosol size distribution. The first three bins represent the Aitken mode and accumulation mode of aerosol. The last bin represents the coarse mode of aerosol. The MOSAIC aerosol scheme includes sulfate, methane sulfonate, nitrate, chloride, carbonate, ammonium, sodium, calcium, black carbon (BC), primary organic mass (OC), liquid water, and other inorganic mass (OIN). The OIN species include silica, other inert minerals, and trace metals. The emitted dust is assigned to the OIN class of MOSAIC to simulate the major aerosol processes. To study the sensitivity of dust simulation to different dust emission schemes and dry deposition schemes, we test two different dust emission schemes (see Sect. 2.2) and three dry deposition schemes (see Sect. 2.3) within MOSAIC.
2.2 Dust emission schemes
Dust emission schemes include empirical schemes and schemes based on dust physical processes. Because of differences in input parameters and formulas to calculate dust flux, dust emission varies among different dust emission schemes. The Goddard Chemistry Aerosol Radiation and Transport (GOCART) dust emission scheme (Ginoux et al., 2001) is an empirical scheme and was implemented in MOSAIC by Zhao et al. (2010). The GOCART dust emission scheme within the MOSAIC aerosol scheme is called by setting dust_opt=13. The University of Cologne (UoC) dust emission schemes (Shao, 2001, 2004; Shao et al., 2011) (Shao schemes) are size-resolved dust emission scheme based on the wind erosion physical theory. The UoC dust emission scheme within the GOCART aerosol scheme is called by setting dust_opt=4. When the UoC dust emission scheme is selected, the user should also choose one of the UoC sub-options by setting dust_scheme=1 for Shao2001, dust_schme=2 for Shao2004, or dust_schme=3 for Shao2011. The Shao dust emission schemes are widely used for dust simulations in East Asia and have been found to perform well in simulating dust emission fluxes (Shao et al., 2011; Su and Fung, 2015; Wu and Lin, 2014), and Shao2011 (Shao et al., 2011) is a simplified version of Shao2004 (Shao, 2004). To test the sensitivity of dust simulation to different dust emission schemes, we implemented the Shao2011 dust emission scheme in the MOSAIC aerosol scheme. Each dust emission scheme is described in detail below.
The formula of vertical dust flux in GOCART is approximated as
where C is an empirical proportionality constant, and S is the source function that is determined by the erodibility factor (see Fig. 1). sp is the fraction of each size class of the emitted dust. u10 is the horizontal wind speed at 10 m. ut is the threshold velocity below which the dust emission does not occur. ut is calculated as
where ut0 is the threshold velocity for dry soil, and w is the soil surface wetness. The formula of ut0 is not from the original GOCART paper (Ginoux et al., 2001) but rather from Marticorena and Bergametti (1995).
where ρp is the density of particles, ρa is the density of air, dp is particle diameter, a equals 1331, x equals 1.56, and b equals 0.38. The original GOCART dust emission scheme in the GOCART aerosol scheme (dust_opt=1) calculates the dust emission flux from 0.2 to 20 µm. For the GOCART dust scheme in the MOSAIC aerosol scheme (dust_opt=13), the total dust emissions from 0.2 to 20 µm are redistributed to the size bins of MOSAIC (0.039–0.156, 0.156–0.625, 0.625–2.5, and 2.5–10.0 µm) with mass fractions of 0 %, 0.38 %, 8.8 %, and 68.0 % (Kok, 2011; Zhao et al., 2013). We note that in addition to the size distribution, the values of empirical proportionality constant C are also different for the two GOCART dust emission scheme options. For dust_opt=13, C value is set to kg s2 m−5, which is consistent with the original GOCART dust emission scheme paper (Ginoux et al., 2001). For dust_opt=1, C value is set to kg s2 m−5.
The Shao2011 dust emission scheme is a size-resolved dust emission scheme based on the wind erosion physical theory. The dust flux is determined by
where F(di) is the dust emission rate of particle size di; cy is the dimensionless coefficient; ηmi is the mass fraction of free dust for a unit soil mass; σm is bombardment efficiency; Q is the saltation flux averaged over the range of sand particle sizes. In Shao2011, the erodibility factor is only used to constrain the potential emission regions. Dust emission is permitted in Shao2011 where the erodibility factor is greater than zero. As the Shao2011 scheme is a size-resolved dust emission scheme, it first calculates the emitted dust from 0.98 um to 20 µm with 40 size bins. Dust emissions from these 40 size bins are then grouped into the four size bins of the MOSAIC aerosol scheme (0.039–0.156, 0.156–0.625, 0.625–2.5, and 2.5–10 µm). The details of the Shao2011 dust emission scheme are described in Appendix A. There is a bug in calculating dust emission flux in Shao2011 scheme reported after WRF-Chem v3.9, and we have already corrected it in our simulation (see Appendix A). We should mention that the Shao2011 dust emission scheme used in this study is based on WRF-Chem v3.9 with some modifications from WRF-Chem v3.7.1. The differences of Shao2011 among different WRF-Chem versions are documented in Appendix B.
2.3 Dry deposition schemes
For dry deposition schemes, dry deposition velocity (Vd) is used to calculate dry deposition flux. Vd is determined by gravitational settling velocity (Vg), aerodynamic resistance (Ra), and surface resistance (Rs). There are three dry deposition schemes available in WRF-Chem coupled with the MOSAIC module and used in this study, and they are referred to as BS95 (Binkowski and Shankar, 1995), PE92 (Peters and Eiden, 1992), and Z01 (Zhang et al., 2001). Each dry deposition scheme will be described in detail below.
In the BS95 scheme (Binkowski and Shankar, 1995), Vd is expressed as
where Ra and Rs are aerodynamic and surface resistance; Vg is the gravitational settling velocity and is given as
where Cc is the Cunningham correction factor as a function of dp and mean free path of air (λ), and μ is the viscosity dynamic of air. The surface resistance is calculated as
where EB is collection efficiency from Brownian diffusion. EB is calculated as follows:
where Sc is the Schmidt number, given by . ν is the kinematic viscosity of air, and D is the particle Brownian diffusivity. EIM is the collection efficiency due to impaction of the particle with the collecting surface (Gallagher, 2002). Impaction occurs when there are changes in the direction of airflow, and particles that cannot follow the flow will collide with the obstacle and stay on the surface due to the inertia (Giardina and Buffa, 2018). EIM is given by
where St is the Stokes number, given by
St is the ratio of the particle stop distance to the characteristic length of the flow and describes the ability of particles to adopt the fluid velocity (Pryor et al., 2008; Seinfeld and Pandis, 2006).
In PE92 scheme (Peters and Eiden, 1992), the dry deposition velocity (Vd) is expressed as
The formula of Vg and Ra is the same as in BS95, but the way to calculate Rs is quite different. In PE 92, Rs is parameterized as
where EIN is collection efficiency from interception, and R is the factor for particle rebound. EIM, EIN, and R are expressed as
z0 is the roughness length, and dp is particle diameter. The Stokes number is given by
u is the horizontal wind velocity and dc is the diameter of the obstacle.
In the Z01 scheme (Zhang et al., 2001), the formula of Vd is the same as in the BS95 scheme (Eq. 5). Surface resistance Rs is calculated as
where γ depends on land use categories (LUCs) and lies between 0.50 and 0.58.
EIM is expressed as
where β equals to 2. α depends on LUC and lies between 0.6 and 100.0. The Stokes number is given by
over vegetated surfaces (Slinn, 1982) and
over smooth surfaces or surfaces with bluff roughness elements (Giorgi, 1988). EIN is the collection efficiency based on the relative dimensions of the particle to the collector diameter (Gallagher, 2002). Interception occurs when particles moving with the mean flow and the distance between an obstacle and particle center is less than half of the diameter. Then the particles will collide with and be collected by the obstacle. EIN is expressed as
over vegetated surfaces and EIN=0 for nonvegetated surfaces, where A is the characteristic radius of collectors. A depends on LUC and lies between 2.0 and 10.0 mm. R is expressed as
The main differences in formulas used to calculate dry deposition velocity for three different dry deposition schemes are listed in Table 2. For Rs, PE92 and Z01 include the collection efficiency from interception (EIN) and the rebound effect (R), while these two are neglected in BS95. For the EIM parameterization, all three schemes use St to parameterize EIM, but the formulas are quite different. BS95 has a different formula from PE92 and Z01, while the PE92 and Z01 have the same formula but with different coefficients. For PE92, the coefficient for EIM is constant for all of the surface types. For Z01, the coefficients α and β for EIM change with different surface types. For the EIN parameterization, BS95 ignores this effect; PE92 and Z01 use different formulas and variables to calculate EIN. When large particles (usually >5 µm) hit the nonsticky surface, they are liable to rebound from the surface if they have sufficient kinetic energy. The rebound factor R represents the fraction of particles that stick to the surface (Seinfeld and Pandis, 2006). For rebound effect, BS95 does not consider it; PE92 and Z01 use the same e-exponential form to calculate the rebound effect with a different coefficient b. For PE92, b is 2.0; for Z01, b is 1.0. In addition, the parameterization of St is quite different for different dry deposition schemes. For BS95, the formulation of St tends to emphasize the nature of the flow field (Binkowski and Shankar, 1995; Pryor et al., 2008). For Z01, the formulation of St is from Slinn (1982) over vegetated surfaces and from Binkowski and Shankar (1995) over smooth surfaces. The formulation of St from Slinn (1982) and Peters and Eiden (1992) are focus on the individual obstacles (Pryor et al., 2008).
2.4 Experiments design
We use WRF-Chem v3.9 with 20 km×20 km horizontal resolution and 35 vertical levels with model top pressure at 50 hPa. The domain covers most of East Asia (14–60∘ N, 74–130∘ E) as shown in Fig. 1. The simulation period is from 26 April to 7 May 2017 with time step of 60 s and frequency of output every hour. The time step between radiation physics calls is 20 min. During this period, a severe dust storm event originated from northwestern China and Outer Mongolia, and air quality deteriorated dramatically in a very short time in downwind areas (Guo et al., 2019; Zhang et al., 2018). Meteorological conditions are initialized and forced at the lateral boundaries using the 6-hourly National Center for Environmental Prediction Final (NCEP/FNL) Operational Global Analysis data at a resolution of . For meteorological conditions (such as wind speed and temperature), we reinitialized every 24 h using NCEP/FNL reanalysis data. For chemistry, the output of the aerosol field (such as the concentration of different aerosol species) from the previous 1 d run was used as the initial chemical conditions for the next 1 d run. Our simulation period is from 26 April to 7 May 2017, and one experiment consists of 12 runs, each of 1 d. In this way, the chemical fields are continuous, and we can also get more reliable meteorological conditions. The MOSAIC aerosol scheme was used for all of the simulations. Simulation results prior to 28 April are treated as model spinup for chemical initial condition and are not included in results presented in Sect. 3. The model results from 1 to 7 May are used for the dust loading and concentration analysis. And the model results from 28 April to 7 May are used for the dust emission analysis as the dust emissions before 1 May also have an influence on the dust concentration during 1–7 May. To study the dust simulation sensitivity to dust emission and dry deposition schemes, we run six experiments with two different dust emission schemes and three dry deposition schemes (See Table 3). The corresponding model configuration for dry deposition processes of the six experiments also listed in Table 3. We note here that the USGS LUCs should be selected for the Z01 dry deposition scheme.
Hourly surface observed PM10 is used to compare with the simulated PM10 from WRF-Chem. In China, hourly surface PM10 concentrations were collected from more than 1000 environmental monitoring stations (locations shown in results section) maintained by the Ministry of Environmental Protection (MEP). The hourly PM10 data from 1 to 7 May 2017 were downloaded from https://quotsoft.net/air/ (last access: 3 April 2020). We colocated the PM10 data with WRF-Chem simulation grids to evaluate model performance with different configurations.
2.5.2 MODIS AOD
Daily aerosol optical depth (AOD) from the Moderate Resolution Imaging Spectroradiometer (MODIS) is used to compare with our simulated AOD from WRF-Chem. The MODIS on board the Aqua satellite was launched by the NASA in 2002, and Aqua is a part of the A-Train satellite constellation. To compare modeled AOD with observations, we use AOD retrievals at 550 nm from MODIS AOD products on Aqua with daily gridded data at a resolution of (MYD08_D3, Collection 6, combined dark target and deep blue AOD). The MODIS Aqua collection daily MYD08_D3 files were obtained from https://ladsweb.nascom.nasa.gov (last access: 23 April 2020). As Aqua passes through every region of Earth at around 13:30 local time, we extract the model simulation results at 13:00 to compare with the daily MODIS AOD. For the model results, first we divided the domain into different time zones according to the longitude. Then the model results at corresponding UTC when the local time is 13:00 are extracted. The colocated model AOD results for each day are used in comparison with daily MODIS AOD.
2.5.3 CALIPSO data
The vertical profile of aerosol extinction coefficient at a wavelength of 532 nm from the Cloud-Aerosol Lidar and Infrared Pathfinder Satellite observation (CALIPSO) satellite is used to evaluate model results. The CALIPSO launched on 28 April 2006 equipped with CALIOP (Cloud-Aerosol Lidar with Orthogonal Polarization). The CALIOP lidar provides an along-track observation of aerosol and cloud vertical profile. The vertical and horizontal resolutions for CALIOP from the surface to 8.2 km are 30 and 333 m, respectively. Above 8.2 km, the vertical and horizontal resolutions are 60 m and 1 km, respectively. We use the CALIPSO level 2 APro product (V4.20) to obtain the aerosol extinction coefficient (CAL_LID_L2_05kmAPro-Standard-V4-20). The CALIPSO data are available at https://www-calipso.larc.nasa.gov/ (last access: 23 April 2020).
3.1 Meteorological conditions
Dust emission and transport processes are closely related to the meteorological conditions. So we first evaluated the model performance in simulating the synoptic conditions. Figure 2 shows the surface meteorological conditions during the dust event. Figure 2a, d, g, and j show the daily mean wind field at 10 m and daily mean temperature at 2 m from NCEP/FNL reanalysis data. The meteorological conditions at 700 hPa are shown in the Supplement (Fig. S1). This dust storm was triggered by the development of a Mongolian cyclone (Figs. S1c and 2d). With the strong northwest and southwest wind near the dust source region, emitted dust was transported to the southeast and northeast of China (Figs. S1c, e and 2d, g). Figure 2b, e, h, and k show the WRF-Chem-simulated daily mean wind field at 10 m and daily mean temperature field at 2 m. Figure. 2c, f, i, and l show the difference in daily mean wind speed at 10 m between WRF-Chem simulation and NCEP/FNL reanalysis data. The WRF-Chem model was able to simulate the wind speed well over the dust source regions (the Taklimakan Desert and the Gobi Desert) and eastern and southern China, where the differences were mostly in the range of −2.0–2.0 m s−1. The wind speed is slightly underestimated near the center of the cyclone (Fig. 2c, f, i, l) and as this is away from dust source regions, we do not expect this underestimation causes large bias in dust emissions. The correlation coefficient (R) and root mean square error (RMSE) between WRF-Chem simulation and FNL reanalysis data for temperature at 2 m, U component of wind, V component of wind, and wind speed at 10 m during simulation period are shown in Table 4. The R for time-averaged temperature at 2 m, U component of wind, V component of wind, and wind speed at 10 m from 1 to 7 May are 1.0, 0.90, 0.86, and 0.82, respectively. The RMSE for time-averaged temperature at 2 m, U component of wind, V component of wind, and wind speed at 10 m from 1 to 7 May are 1.03, 1.08, 0.98, and 1.11, respectively. The R for temperature, U component of wind, V component of wind, and wind speed at 700 hPa from 1 to 7 May are 1.0, 0.94, 0.91, and 0.95, respectively (Table S1 in the Supplement). The RMSE for temperature, U component of wind, V component of wind, and wind speed at 700 hPa from 1 to 7 May are 0.67, 2.34, 2.70, and 1.76, respectively (Table S1). Overall, the correlation coefficients are generally large, and the RMSEs are generally small. This indicates that WRF-Chem performed well in simulating the meteorological conditions. We also compared the difference in the meteorological conditions in our six experiments and found that the difference is negligible (Fig. S2 and Table S2).
3.2 Dust simulation sensitivity to dust emission schemes
In this section, we examine the changes in the simulated dust loading using different dust emission schemes. Figure 3 shows simulated mean dust loading for six experiments over the 7 d simulation period 1–7 May 2017. When using the same dry deposition scheme (BS95, PE92, or Z01), different dust emission schemes give very different dust spatial distributions. Compared with the Shao2011 scheme, GOCART has higher dust loading over the Taklimakan desert (TD) but has relatively lower dust loading over the Gobi Desert (GD), the south of Outer Mongolia, and most parts of northern China. The difference in the spatial distribution of dust loading is mainly caused by the different spatial distributions of dust emission flux from dust emission schemes, as shown in Fig. 4. As the dust emissions before 1 May also have an influence on the dust loading during 1–7 May, the total dust emissions from 28 April to 7 May are analyzed. The total dust emissions from 00:00 UTC on 28 April to 23:00 UTC on 7 May over the GD from GOCART and Shao2011 are 4.90 and 13.88 Tg, respectively. The total dust emissions from 00:00 UTC on 28 April to 23:00 UTC on 7 May over the TD from GOCART and Shao2011 are 7.16 and 2.75 Tg, respectively. Over the GD, the Shao2011 scheme has higher dust emission than GOCART, while over the TD, the GOCART scheme has higher dust emission than Shao2011 (Fig. 4c).
Figure 5a, c, e, and g show the spatial distribution of friction velocity, threshold friction velocity, the difference between friction velocity and threshold friction velocity, and the dust emission flux from Shao2011 at 06:00 UTC on 3 May. The areas where the friction velocity is greater than the threshold friction velocity is mainly located in the west of Inner Mongolia and the south of Outer Mongolia (Fig. 5e). This is consistent with Fig. 5g. When the friction velocity is larger than threshold friction velocity, dust can be emitted from the surface. Figure 5b, d, f, and h show the spatial distribution of wind speed at 10 m, threshold velocity, the difference between wind speed at 10 m and threshold velocity, and the dust emission flux from the GOCART dust emission scheme. Different from Shao2011, the dust emission regions from GOCART are not only determined by wind speed but also constrained by the erodibility factor (Eq. 1). From Fig. 5f, the threshold velocity is much smaller than the wind speed at 10 m in most areas. In these areas, GOCART uses Eq. (1) to calculate the dust emission flux, and the source function S depends on the erodibility factor. The dust emission flux in GOCART is directly scaled by erodibility factor. Figure 1 shows the erodibility factor, which describes the fraction of erodible surface in each grid cell. As shown in Fig. 5h, dust emission occurs where the wind speed is high and the erodibility factor is larger than 0.
Over the TD, Shao2011 produces lower dust emission flux than GOCART. One reason may be the formula used to calculate the threshold velocity (Eq. 3). The formula used to calculate threshold velocity is from Marticorena and Bergametti (1995), which was originally designed to calculate threshold friction velocity (see LeGrand et al., 2019 for details). This inconsistency leads to a very small threshold velocity in GOCART, which may result in dust emission at low wind speed. Another reason may be the incorrect soil particle size distribution over the TD (Wu and Lin, 2014). The incorrect soil particle size distribution can lead to the unreasonable dust emission flux in Shao2011 over the TD. Over the GD, the GOCART scheme has lower dust emission than the Shao2011 scheme. As mentioned by Su and Fung (2015), the erodibility factor over the GD is highly underestimated and needs to be improved for the GOCART dust emission scheme.
As we mentioned in Sect. 2.2.2, Shao2011 used in this study is based on WRF-Chem v3.9 with some modifications from WRF-Chem v3.7.1. The modified Shao2011 simulates better dust loading than the original Shao2011 scheme in WRF-Chem v3.9 (not shown). Simulated dust emission fluxes are quite different when using two versions of the Shao2011 scheme in WRF-Chem v3.9 and WRF-Chem v3.7.1, which is mainly caused by different soil particle size distributions in two versions. The differences in Shao2011 among different WRF-Chem versions are documented in Appendix B.
3.3 Dust simulation sensitivity to dry deposition schemes
In this section, we analyze dust simulation sensitivity to different dry deposition schemes using the six experiments. For simulated dust loading using the GOCART dust emission scheme (Fig. 3a–c), compared to the BS95 dry deposition scheme, PE92 and Z01 produce higher dust loading over the dust source regions and remote regions. The relative differences in mean dust loading from PE92 and Z01 relative to BS95 are 20 % and 59 %, respectively. As for the simulated dust loading using the Shao2011 dust emission scheme (Fig. 3d–f), PE92 and Z01 schemes also produce higher dust loading than the BS95 scheme, and the differences relative to BS95 are 72 % and 116 %, respectively. This indicates that dust simulation is very sensitive to dry deposition schemes.
Figure 6a shows the modeled dry deposition velocity over desert surface. As desert dust mass is mainly concentrated in the large particle size range, our dry deposition analysis focuses on the coarse mode (2.5–10 µm) (Kok, 2011; Zhao et al., 2013). The reference diameter of the coarse mode is defined at 5 µm (Fig. 6). BS95 produces larger Vd than PE92 and Z01 in the coarse aerosol mode. Larger Vd leads to larger dry deposition and thus lower dust loading, consistent with the lower simulated dust loading from the BS95 scheme discussed above (Fig. 3). In Eq. (5), the dry deposition velocity is comprised of gravitational velocity, aerodynamic resistance, and surface resistance. The diversity of different dry deposition schemes mainly comes from the way to parameterize surface resistance, and differences from gravitational settling and aerodynamics resistance are small (not shown), consistent with previous studies (e.g., Bergametti et al., 2018). Figure 6b shows the surface resistance from different schemes as a function of particle diameter (dp). In the coarse aerosol mode, Z01 produces the largest surface resistance, followed by PE92 and BS95. Larger surface resistance causes smaller dry deposition velocity in Z01, thus larger dust concentration as shown in Fig. 3.
The surface collection efficiency is comprised of Brownian diffusion, impaction, and interception and is corrected for particle rebound (see Eq. 12). Collection from Brownian diffusion is most important for the smaller particles, while collection from impaction and interception play a more important role for large particles in surface collection processes. Figure 6c shows the surface collection efficiency from impaction (EIM) from different schemes as a function of particle diameter. BS95 gives the largest EIM, and Z01 gives the smallest. Based on field observation data, Slinn (1982) used a semiempirical fit for smooth surface (Eq. 9), and Binkowski and Shankar (1995) adopted this formula for EIM and used it for all land surface types. Peters and Eiden (1992) uses Eq. (19) to describe EIM, with α equals 0.8, and β equals 2 to get the best fit for the data collected over a spruce forest (Eq. 13). In Zhang et al. (2001) scheme, α varies with LUC, and β is chosen as 2 (Eq. 19). For BS95 and PE92, the formula of EIM is derived from a specific land surface type, but they have been applied to all land surface types in WRF-Chem. This may lead to large uncertainties for dry deposition over the whole domain with different surface types. As the EIM of Z01 varies with LUC, Z01 may have a better physical treatment of EIM than the other two dry deposition schemes.
Figure 6d shows the surface collection efficiency from interception (EIN). EIN depends on the particle diameter and the characteristic radius of the collectors (Seinfeld and Pandis, 2006). EIN is important for large particles on hairs at the leaf surface and is negligible over nonvegetated surfaces such as the desert surface we analyzed here (Chamberlain, 1967; Slinn, 1982; Zhang et al., 2001). In BS95, the effect of interception is not considered. In the original PE92 scheme as described in Peters and Eiden (1992), EIN is also not considered. But in the PE92 scheme used in WRF-Chem, EIN increases with particle diameter as in Eq. (14). In Z01, the effect of interception is considered as Eq. (22) over vegetated surface and is not considered for nonvegetated surface (as shown in Fig. 6d over desert surface type). The parameterization of EIN partially results in the difference in surface resistance between PE92 and the other two dry deposition schemes.
Figure 6e shows the rebound factor from different dry deposition schemes. Rebound and resuspension have long been recognized as a mechanism by which the surface can act as sources of particles (Pryor et al., 2008). Due to limited knowledge of particle rebound and resuspension processes, most dry deposition models adopted the form of the rebound effect as suggested by Slinn (1982) (Zhang and Shao, 2014; Zhang et al., 2001), while some dry deposition schemes do not include the rebound effect with R=1.0 (Binkowski and Shankar, 1995; Petroff and Zhang, 2010; Zhang and He, 2014). BS95 does not consider the rebound effect. b is equal to 2.0 for the PE92 scheme and 1.0 for the Z01 scheme. Another difference between PE92 and Z01 is the threshold particle diameter for including the rebound effect. Rebound effect is included for PE92 when particles are larger than 0.625 µm and for Z01 when particles are larger than 2.5 µm. In summary, the smaller EIM and rebound factor lead to larger Rs in Z01, while the larger EIM leads to smaller Rs in BS95, and the moderate EIM and rebound effect give a moderate Rs for PE92.
Figure 6f shows the Stokes number from different dry deposition schemes. Over smooth surfaces, the formula of St for BS95 and Z01 is the same, as shown in Eq. (10). In PE92, St is calculated using Eq. (16), which is similar to the formula used in Slinn (1982). BS95 and Z01 schemes give a larger St than PE92. The Stokes number is used to calculate both R and EIM. The difference in Stokes numbers and the different formulas of R and EIM lead to the different R and EIM among different dry deposition schemes (Fig. 6c and e).
Our discussion indicates that Z01 has a better physical treatment of dry deposition velocity, as Z01 considers the rebound effect, and EIM changes with LUC. The Z01 scheme has also been documented to agree better with measured dry deposition fluxes and dry deposition velocity (e.g., Zhang et al., 2012; Connan et al., 2018). Zhang et al. (2012) compared the dry deposition fluxes measured at five sites in Taiwan with the modeled dry deposition fluxes and found that the measured dry deposition fluxes can be reproduced reasonably well using the Z01 scheme. Connan et al. (2018) conducted experimental campaigns on-site to determine dry deposition velocity of aerosols and found that the Z01 scheme is most suitable for operational use in the size range 0.2–10 µm. All of these indicate that the Z01 dry deposition scheme is more physically meaningful than other two dry deposition schemes.
3.4 Comparisons with observations
To better evaluate the performance of different experiments, we compared the model results with observations. Figure 7 shows hourly observed PM10 concentrations over observational sites at 02:00 UTC on 4 May 2017 (10:00 Beijing Time, BJT, on 4 May 2017). Very high PM10 values (>1000 µg m−3) are observed in northern China. Figure 8 compares simulated PM10 in six experiments with observed PM10. During the comparison, the observational sites closest to the model grids are paired up. The correlation coefficients (R), root mean square errors (RMSE) between model and observations, and the mean simulated and observed PM10 for all of the sites over the five regions during the 7 d period 1–7 May are marked in Fig. 8. The simulated PM10 of all of the six experiments have obviously underestimated the observations. Among all of these experiments, GOBS95 has the lowest average PM10 concentration, with a value of 26.45 µg m−3, and S11Z01 has the largest one, with a value of 105.17 µg m−3, the closest one to the observed mean value of 172.70 µg m−3. S11Z01 gives a large R of 0.77 and the smallest RMSE of 96.14 compared to other experiments. Table 5 shows the R and RMSE between the model and observations for PM10 for six experiments over five subregions and over the whole of China. Over the TD, GOBS95 gives the largest R and smallest RMSE. Over the GD, GOZ01 and S11Z01 give a better performance compared with other experiments. For other regions (North China Plain, NCP, Northeast Plain, NEP, and the middle and lower reaches of the Yangtze River plain, YR), S11Z01 gives a relatively larger R and smallest RMSE. For all of the stations in total, S11Z01 gives a larger R of 0.83 and the smallest RMSE of 82.98. Overall, the S11Z01 experiment has the best performance for simulating this dust storm.
Figure 9 shows the MODIS-observed daily mean AOD and WRF-Chem-simulated AOD over the simulation period 1–5 May. For strong dust storms like the one we examined here, dust particles contribute the most to AOD, and AOD therefore can represent the dust loading in the atmosphere. To match the MODIS AOD observation time, simulated AOD at 13:00 local time is used for comparison (see Sect. 2.5 for details). For each grid with observed AOD from MODIS, the average value of simulated AOD from WRF-Chem in this grid is calculated. Grid points without valid MODIS AOD retrieval are masked for both observational and model results in Fig. 9. A major dust emission event occurred over the GD on 3 May (Fig. 9c). Shao2011 simulated well the dust emission event over the GD on 3 May (Fig. 9m), while GOCART obviously underestimated dust emission over the GD (Fig. 9h). On 4 May, emitted dust from the GD was transported to northeast China, and the highest AOD values for this case study were observed in northern China (Fig. 9d). As the GD is the main dust source region of this dust storm, Shao2011 correctly captured the emission phase of this dust event. Simulated AOD values from the S11Z01 configuration produced the closest match to the observed daily MODIS AOD with respect to the magnitude and spatial pattern (Figs. 9n and S3). For a more quantitative comparison, Table 6 shows the correlation coefficient (R) and root mean square error (RMSE) between the model and observed AOD for six experiments during 1–7 May. Overall, the S11Z01 experiment gives a larger correlation coefficient, and the RMSE is almost the same among different experiments; the correlation coefficient is still lower than 0.5. The low correlation may partly come from the spatial and temporal limitation of satellites and the difficulties to retrieve aerosol in the vicinity of clouds for satellites.
To evaluate the model performance in simulating the vertical profile of dust aerosol, we compared the extinction coefficient from the model and from CALIPSO (Fig. 10). Figure 10 shows the simulated and observed aerosol extinction profiles at 532 nm at 18:00 UTC on 4 May. The trajectory of CALIPSO passes East Asia (Fig. 10d). All of the six experiments show the similar dust location in the atmosphere, which is consistent with the CALIPSO observation. However, the magnitude of dust concentration differs substantially. The simulated extinction coefficients using GOCART dust emission schemes are significantly underestimated compared to the CALIPSO observation (Fig. 10a, b, and c), while the modeled extinction coefficients using Shao2011 dust emission scheme agrees better with observation though they are still underestimated (Fig. 10e, f and g). Among all of the six experiments, results from S11Z01 agree the best with observation.
In summary, both ground and satellite observations indicate that the S11Z01 experiment yields the best performance in simulating this dust storm. As we discussed in Sect. 3.2, the Z01 dry deposition scheme indeed has a better physical treatment and performs better than some other dry deposition schemes.
In this study, we analyzed the dust simulation sensitivity to different dust emission schemes and dry deposition schemes. In order to compare different dust emission schemes, the Shao2011 dust emission scheme has been implemented into the MOSAIC aerosol scheme in WRF-Chem v3.9. Six model experiments were conducted to simulate the dust storm in May 2017 over East Asia, with two dust emission schemes (GOCART and Shao2011) and three dry deposition schemes (BS95, PE92, and Z01). The simulation results of different experiments were evaluated against surface and satellite observations.
Our results show that dust loading is very sensitive to different dry deposition schemes. The relative difference in dust loading in different experiments range from 20 % to 116 % when using different dry deposition schemes. The difference in dry deposition velocity in different dry deposition schemes comes from the parameterization of surface resistance, and difference in surface resistance mainly comes from the parameterization of collection efficiency from impaction and rebound effect. In addition, different dust emission schemes result in different spatial distributions of dust loading, as dust emission fluxes in dust source regions differ substantially among different dust emission schemes, which is mainly attributed to differences in the threshold conditions for dust emission and in formulas and parameters for calculating dust emission flux. We noted that the Shao2011 dust emission scheme is different among different WRF-Chem versions, and a significant difference exists in the simulated dust emission fluxes between WRF-Chem v3.9 and WRF-Chem v3.7.1, which is mainly caused by differences in soil particle size distributions used in two versions (see Appendix B).
Compared with both surface PM10 station observations and MODIS AOD, the Shao2011 dust emission scheme coupled with the Z01 dry deposition scheme produces the best simulation for the dust storm in East Asia. For PM10, the S11Z01 experiment gives a larger R of 0.83 and the smallest RMSE of 82.98 of all of the stations (Table 5). The spatial distribution of AOD during the simulation period obtained by S11Z01 agrees the best with MODIS AOD (Fig. 9), with the largest R and a relatively small RMSE (Table 6). Our analysis indicates Z01 accounts for the rebound effect, and EIM changes with LUC and therefore has a better physical treatment of dry deposition velocity than the two other dry deposition schemes. Previous studies have also shown that the Z01 scheme agrees better with measured dry deposition fluxes and dry deposition velocity (e.g., Zhang et al., 2012; Connan et al., 2018). The Shao2011 dust emission scheme has larger dust emission fluxes than the GOCART dust emission scheme over the Gobi Desert, and the transport of dust emitted from the Gobi Desert is the most important source of dust weather in northern China (Chen et al., 2017b). Compared with daily MODIS AOD (Fig. 9), our results indicate that dust emission from Shao2011 is better for this dust event, in terms of dust spatial and temporal distributions. We note that our results are obtained from simulations of a dust storm over a short period, and longer simulations are desirable in the future to test whether the optimal scheme here still produces best simulations.
This study highlights the importance of dry deposition process in dust simulation. Future studies on dust simulation should pay attention to improve dry deposition schemes as well as the dust emission schemes. Additional field measurements of dry deposition process and comparisons with model results are required to reduce the uncertainties in dust simulation.
Here we describe the Shao2011 dust emission scheme in more detail as a supplement to the Sect. 2.2.2 of this article. The total saltation flux Q in Eq. (4) is calculated as
where d1 and d2 define the upper and lower limits of saltation particle size. pm(d) is the minimally disturbed particle-size distribution. The saltation flux Q for each particle size d is calculated as
where cf is the fraction of vegetation cover, (1−cf) means the fraction of erodible surface area, and c0 is a coefficient. u∗t is the threshold friction velocity, and u∗ is the friction velocity. When u∗ is larger than u∗t, it calculates the dust emission flux. Before WRF-Chem v4.0, there is a bug in calculating the saltation flux Q(ds) in Shao2011. The last term is miscalculated as 1 in WRF-Chem codes (LeGrand et al., 2019). In WRF-Chem v4.0 and later versions, this bug is fixed, and we also fixed this bug in our simulations.
The threshold friction velocity u∗t is calculated as
where u∗t0 is the ideal threshold friction velocity when soil is dry, bare, and free of crust and salt, fλ is the correction functions for surface roughness, and fθ is the correction functions for soil moisture. The ideal threshold friction u∗t0 is calculated as
where a1 and a2 are constant. ρp and ρa are particle and air density. d is the particle diameter.
The correction functions for surface roughness fλ are calculated as
where m is a constant, σ is the ratio of roughness-element basal area to frontal area, λ is the frontal area index, and β0 is the ratio of the drag coefficient of an isolated roughness element on the surface to the drag coefficient of the substrate surface itself.
The mass fraction of free dust ηmi is calculated as
where pm(d) is the minimally disturbed particle-size distribution, which is regarded as a composite of several lognormal distributions, and pm(d) is expressed as
Soil samples collected from experiment sites are used to determine the particle-size distribution.
As we noted in Sect. 3.2, the Shao2011 scheme in different versions of WRF-Chem can produce significantly different dust emission fluxes. Here we document differences in Shao2011 among different WRF-Chem versions.
The first difference is c0, where c0 is a coefficient used to calculate the saltation flux as in Eq. (A2). In versions before WRF-Chem v3.8, c0 is equal to 0.5; in WRF-Chem v3.8 and later versions, c0 is equal to 2.3 (Table B1).
The second difference is β0. β0 is a coefficient used to calculate the correction function for surface roughness fλ in Eq. (A5). In versions before WRF-Chem v3.8, β0 is 90; in WRF-Chem v3.8 and later versions, β0 is 200 (Table B1).
The third difference is caused by the minimally disturbed particle-size distribution pm(d) (see Eq. A6). pm(d) is used to calculate the free dust fraction ηmi (see Eq. A6). Free dust fraction is the fraction of dust that has low enough binding energy so that it can be easily lifted from the surface by either aerodynamic forces or mechanical abrasion (Shao, 2001). The ηmi is used to calculate the dust emission rate in Eq. (4). Twelve soil types are included in all WRF-Chem versions. In WRF-Chem v3.8 and later versions, each soil type has a corresponding pm(d) as listed in Table B1 from Shao et al. (2010); in versions before WRF-Chem v3.8, there are only four pm(d) as listed in Table 1 from Shao (2004) for 12 soil types (Fig. B1). For example, (f) sand and (g) loamy sand soil types use the same free dust fraction distribution in versions before WRF-Chem v3.8. As shown in Table B2, the loam and clay loam are the two soil types with the largest percentage, while the other soil types account for a very small percentage. From Fig. B1c and e, for loam and clay loam soil types, the free dust fraction is very small in the particle size range 0–10 um in WRF-Chem v3.8 and later versions, almost all close to 0; however, in the versions before WRF-Chem v3.8, the free dust fraction is relatively high. In different WRF-Chem versions, the total saltation flux Q is the same, but dust emission flux F(di) is different due to different free dust fractions (see Eq. 4). With a smaller free dust fraction, the dust emission flux is smaller in WRF-Chem v3.8 and later versions.
To examine the importance of these changes, we run four experiments to quantify the contribution of each factor (Table B3). For the control run, c0 is 2.3, β0 is 200, and pm(d) has 12 distributions based on WRF-Chem v3.8 or later versions. For the case 1 experiment, β0 is changed to 90, the one used in WRF-Chem v3.7.1 and all other parameters are kept the same as in the control run. The dust emission in case 1 is 1.35 times higher than the control run. For the case 2 experiment, c0 is changed to 0.5, the one used in WRF-Chem v3.7.1, and all other parameters remain the same as in the control run. The dust emission in case 2 is 21 % of the dust emission of the control run. For the case 3 experiment, pm(d) is adopted from WRF-Chem v3.7.1 and has four distributions, and all other parameters remain the same as in the control run. The dust emission in case 3 is 13 times higher than the control run. This indicates that the difference in dust emission between different versions of Shao2011 scheme is mainly caused by the change in pm(d). As pm(d) is determined by soil particle size distribution, this also highlights the need to improve the accuracy of soil texture.
We should mention that the Shao2011 dust emission scheme we used in this study is based on WRF-Chem v3.9 with the soil particle size distribution from WRF-Chem v3.7.1, which simulates better dust loading compared with observations. Compared with the original Shao2011 scheme in WRF-Chem v3.9, the total dust emission simulated in our experiments during 1–7 May is 13 times higher.
The source code of WRF-Chem is available at http://www2.mmm.ucar.edu/wrf/users/download/get_sources.html (last access: 23 April 2020; WRF Users page, 2020). The modified WRF-Chem v3.9 with the Shao2011 dust emission scheme implemented in the MOSAIC aerosol scheme is available upon request to the corresponding author.
The 6-hourly National Center for Environmental Prediction Final (NCEP/FNL) Operational Global Analysis data at a resolution of can be obtained from https://doi.org/10.5065/D6M043C6 (NCEP/NWS/NOAA/U.S. Department of Commerce, 2000). The observed PM10 data are collected from the national air quality real time release platform at http://188.8.131.52:20035/ (last access: 23 April 2020; China National Environmental Monitoring Centre, 2020a). The historical data of air quality used in this study can be downloaded from https://quotsoft.net/air/ (last access: 23 April 2020; China National Environmental Monitoring Centre, 2020b). Daily MYD08_D3 files from the MODIS on board the Aqua satellite can be obtained from https://doi.org/10.5067/MODIS/MOD08_D3.006 (Platnick et al., 2015).
The supplement related to this article is available online at: https://doi.org/10.5194/gmd-13-2125-2020-supplement.
MW and YZ conceived the idea and designed the model experiments. YZ performed the simulations, conducted the analysis, and wrote the paper. CZ and SC provided guidelines for the dust simulations in WRF-Chem and helped in the interpretation of the results. XH and ZL provided guidelines for the WRF-Chem simulations and contributed to the model experiments design. YG helped in the PM10 data processing and usage. Everyone edited the paper.
The authors declare that they have no conflict of interest.
This research is supported by the Collaborative Innovation Center of Climate Change, Jiangsu Province. The numerical calculations in this paper have been done on the computing facilities in the High Performance Computing Center (HPCC) of Nanjing University.
This research has been supported by the Minister of Science and Technology of China (grant nos. 2016YFC0200503 and 2017YFA0604002), the National Natural Science Foundation of China (grant nos. 41575073, 41621005, 41775146, and 91744208), and the Fundamental Research Funds for the Central Universities.
This paper was edited by Samuel Remy and reviewed by Sandra LeGrand and Sagar Parajuli.
Abuduwaili, J., Liu, D., and Wu, G.: Saline dust storms and their ecological impacts in arid regions, J. Arid Land, 2, 144–150, https://doi.org/10.3724/sp.j.1227.2010.00144, 2010.
Bergametti, G., Marticorena, B., Rajot, J. L., Foret, G., Alfaro, S. C., and Laurent, B.: Size-Resolved Dry Deposition Velocities of Dust Particles: In Situ Measurements and Parameterizations Testing, J. Geophys. Res.-Atmos., 123, 11080–11099, https://doi.org/10.1029/2018JD028964, 2018.
Binkowski, F. S. and Shankar, U.: The Regional Particulate Matter Model: 1. Model description and preliminary results, J. Geophys. Res., 100, 26191, https://doi.org/10.1029/95JD02093, 1995.
Chamberlain, A. C.: Transport of Lycopodium spores and other small particles to rough surfaces, P. Roy. Soc. Lond. Ser. A., 296, 45–70, https://doi.org/10.1098/rspa.1967.0005, 1967.
Chen, F. and Dudhia, J.: Coupling an Advanced Land Surface–Hydrology Model with the Penn State–NCAR MM5 Modeling System. Part I: Model Implementation and Sensitivity, Mon. Weather Rev., 129, 569–585, https://doi.org/10.1175/1520-0493(2001)129<0569:CAALSH>2.0.CO;2, 2001.
Chen, S., Huang, J., Zhao, C., Qian, Y., Leung, L. R., and Yang, B.: Modeling the transport and radiative forcing of taklimakan dust over the tibetan plateau: A case study in the summer of 2006, J. Geophys. Res.-Atmos., 118, 797–812, https://doi.org/10.1002/jgrd.50122, 2013.
Chen, S., Huang, J., Qian, Y., Zhao, C., Kang, L., Yang, B., Wang, Y., Liu, Y., Yuan, T., Wang, T., Ma, X., and Zhang, G.: An overview of mineral dust modeling over East Asia, J. Meteorol. Res., 31, 633–653, https://doi.org/10.1007/s13351-017-6142-2, 2017a.
Chen, S., Huang, J., Kang, L., Wang, H., Ma, X., He, Y., Yuan, T., Yang, B., Huang, Z., and Zhang, G.: Emission, transport, and radiative effects of mineral dust from the Taklimakan and Gobi deserts: comparison of measurements and model results, Atmos. Chem. Phys., 17, 2401–2421, https://doi.org/10.5194/acp-17-2401-2017, 2017b.
Chen, S., Zhang, X., Lin, J., Huang, J., Zhao, D., Yuan, T., Huang, K., Luo, Y., Jia, Z., Zang, Z., Qiu, Y., and Xie, L.: Fugitive Road Dust PM2.5 Emissions and Their Potential Health Impacts, Environ. Sci. Technol., 53, 8455–8465, https://doi.org/10.1021/acs.est.9b00666, 2019.
China National Environmental Monitoring Centre: Site-observed surface hourly PM10 concentration across China, available at: http://184.108.40.206:20035/, last access: 23 April 2020a.
China National Environmental Monitoring Centre: Site-observed surface hourly PM10 concentration across China, available at: https://quotsoft.net/air/, last access: 23 April 2020b.
Connan, O., Pellerin, G., Maro, D., Damay, P., Hébert, D., Roupsard, P., Rozet, M., and Laguionie, P.: Dry deposition velocities of particles on grass?: Field experimental data and comparison with models, J. Aerosol Sci., 126, 58–67, https://doi.org/10.1016/j.jaerosci.2018.08.004, 2018.
Creamean, J. M., Suski, K. J., Rosenfeld, D., Cazorla, A., DeMott, P. J., Sullivan, R. C., White, A. B., Ralph, F. M., Minnis, P., Comstock, J. M., Tomlinson, J. M. and Prather, K. A.: Dust and Biological Aerosols from the Sahara and Asia Influence Precipitation in the Western U.S, Science, 339, 1572–1578, https://doi.org/10.1126/science.1227279, 2013.
DeMott, P. J., Prenni, A. J., Liu, X., Kreidenweis, S. M., Petters, M. D., Twohy, C. H., Richardson, M. S., Eidhammer, T. and Rogers, D. C.: Predicting global atmospheric ice nuclei distributions and their impacts on climate, P. Natl. Acad. Sci. USA, 107, 11217–11222, https://doi.org/10.1073/pnas.0910818107, 2010.
Gallagher, M. W.: Measurements and parameterizations of small aerosol deposition velocities to grassland, arable crops, and forest: Influence of surface roughness length on deposition, J. Geophys. Res., 107, 4154, https://doi.org/10.1029/2001JD000817, 2002.
Giardina, M. and Buffa, P.: A new approach for modeling dry deposition velocity of particles, Atmos. Environ., 180, 11–22, https://doi.org/10.1016/j.atmosenv.2018.02.038, 2018.
Ginoux, P., Chin, M., Tegen, I., Prospero, J. M., Holben, B., Dubovik, O., and Lin, S.-J.: Sources and distributions of dust aerosols simulated with the GOCART model, J. Geophys. Res.-Atmos., 106, 20255–20273, https://doi.org/10.1029/2000JD000053, 2001.
Giorgi, F.: A particle dry-deposition parameterization scheme for use in tracer transport models, J. Geophys. Res., 91, 9794, https://doi.org/10.1029/JD091iD09p09794, 1986.
Giorgi, F.: Dry deposition velocities of atmospheric aerosols as inferred by applying a particle dry deposition parameterization to a general circulation model, Tellus B, 40B, 23–41, https://doi.org/10.1111/j.1600-0889.1988.tb00210.x, 1988.
Grell, G. A. and Freitas, S. R.: A scale and aerosol aware stochastic convective parameterization for weather and air quality modeling, Atmos. Chem. Phys., 14, 5233–5250, https://doi.org/10.5194/acp-14-5233-2014, 2014.
Grell, G. A., Peckham, S. E., Schmitz, R., McKeen, S. A., Frost, G., Skamarock, W. C., and Eder, B.: Fully coupled “online” chemistry within the WRF model, Atmos. Environ., 39, 6957–6975, https://doi.org/10.1016/j.atmosenv.2005.04.027, 2005.
Guo, P., Yu, S., Wang, L., Li, P., Li, Z., Mehmood, K., Chen, X., Liu, W., Zhu, Y., Yu, X., Alapaty, K., Lichtfouse, E., Rosenfeld, D. and Seinfeld, J. H.: High-altitude and long-range transport of aerosols causing regional severe haze during extreme dust storms explains why afforestation does not prevent storms, Environ. Chem. Lett., 17, 1333–1340, https://doi.org/10.1007/s10311-019-00858-0, 2019.
Hansen, J., Sato, M., and Ruedy, R.: Radiative forcing and climate response, J. Geophys. Res.-Atmos., 102, 6831–6864, https://doi.org/10.1029/96JD03436, 1997.
Hicks, B. B., Saylor, R. D., and Baker, B. D.: Dry deposition of particles to canopies-A look back and the road forward, J. Geophys. Res.-Atmos., 121, 14691–14707, https://doi.org/10.1002/2015JD024742, 2016.
Hofer, J., Althausen, D., Abdullaev, S. F., Makhmudov, A. N., Nazarov, B. I., Schettler, G., Engelmann, R., Baars, H., Fomba, K. W., Müller, K., Heinold, B., Kandler, K., and Ansmann, A.: Long-term profiling of mineral dust and pollution aerosol with multiwavelength polarization Raman lidar at the Central Asian site of Dushanbe, Tajikistan: case studies, Atmos. Chem. Phys., 17, 14559–14577, https://doi.org/10.5194/acp-17-14559-2017, 2017.
Hong, S.-Y., Noh, Y., and Dudhia, J.: A New Vertical Diffusion Package with an Explicit Treatment of Entrainment Processes, Mon. Weather Rev., 134, 2318–2341, https://doi.org/10.1175/MWR3199.1, 2006.
Huneeus, N., Schulz, M., Balkanski, Y., Griesfeller, J., Prospero, J., Kinne, S., Bauer, S., Boucher, O., Chin, M., Dentener, F., Diehl, T., Easter, R., Fillmore, D., Ghan, S., Ginoux, P., Grini, A., Horowitz, L., Koch, D., Krol, M. C., Landing, W., Liu, X., Mahowald, N., Miller, R., Morcrette, J.-J., Myhre, G., Penner, J., Perlwitz, J., Stier, P., Takemura, T., and Zender, C. S.: Global dust model intercomparison in AeroCom phase I, Atmos. Chem. Phys., 11, 7781–7816, https://doi.org/10.5194/acp-11-7781-2011, 2011.
Iacono, M. J., Delamere, J. S., Mlawer, E. J., Shephard, M. W., Clough, S. A., and Collins, W. D.: Radiative forcing by long-lived greenhouse gases: Calculations with the AER radiative transfer models, J. Geophys. Res.-Atmos., 113, 2–9, https://doi.org/10.1029/2008JD009944, 2008.
Jiménez-Guerrero, P., Pérez, C., Jorba, O., and Baldasano, J. M.: Contribution of Saharan dust in an integrated air quality system and its on-line assessment, Geophys. Res. Lett., 35, 1–6, https://doi.org/10.1029/2007GL031580, 2008.
Kang, J.-Y., Yoon, S.-C., Shao, Y. and Kim, S.-W.: Comparison of vertical dust flux by implementing three dust emission schemes in WRF/Chem, J. Geophys. Res., 116, 1–18, https://doi.org/10.1029/2010JD014649, 2011.
Khan, T. R. and Perlinger, J. A.: Evaluation of five dry particle deposition parameterizations for incorporation into atmospheric transport models, Geosci. Model Dev., 10, 3861–3888, https://doi.org/10.5194/gmd-10-3861-2017, 2017.
Kok, J. F.: A scaling theory for the size distribution of emitted dust aerosols suggests climate models underestimate the size of the global dust cycle, P. Natl. Acad. Sci. USA, 108, 1016–1021, https://doi.org/10.1073/pnas.1014798108, 2011.
Kok, J. F., Ridley, D. A., Zhou, Q., Miller, R. L., Zhao, C., Heald, C. L., Ward, D. S., Albani, S., and Haustein, K.: Smaller desert dust cooling effect estimated from analysis of dust size and abundance, Nat. Geosci., 10, 274–278, https://doi.org/10.1038/ngeo2912, 2017.
LeGrand, S. L., Polashenski, C., Letcher, T. W., Creighton, G. A., Peckham, S. E., and Cetola, J. D.: The AFWA dust emission scheme for the GOCART aerosol model in WRF-Chem v3.8.1, Geosci. Model Dev., 12, 131–166, https://doi.org/10.5194/gmd-12-131-2019, 2019.
Marticorena, B. and Bergametti, G.: Modeling the atmospheric dust cycle: 1. Design of a soil-derived dust emission scheme, J. Geophys. Res., 100, 16415, https://doi.org/10.1029/95JD00690, 1995.
Morrison, H., Thompson, G., and Tatarskii, V.: Impact of Cloud Microphysics on the Development of Trailing Stratiform Precipitation in a Simulated Squall Line: Comparison of One- and Two-Moment Schemes, Mon. Weather Rev., 137, 991–1007, https://doi.org/10.1175/2008mwr2556.1, 2008.
NCEP/NWS/NOAA/U.S. Department of Commerce: NCEP FNL Operational Model Global Tropospheric Analyses, continuing from July 1999, Research Data Archive at the National Center for Atmospheric Research, Computational and Information Systems Laboratory, https://doi.org/10.5065/D6M043C6 (last access: 23 April 2020), 2000.
Ozer, P., Laghdaf, M. B. O. M., Lemine, S. O. M., and Gassani, J.: Estimation of air quality degradation due to Saharan dust at Nouakchott, Mauritania, from horizontal visibility data, Water. Air. Soil Pollut., 178, 79–87, https://doi.org/10.1007/s11270-006-9152-8, 2007.
Peters, K. and Eiden, R.: Modelling the dry deposition velocity of aerosol particles to a spruce forest, Atmos. Environ. Part A, 26, 2555–2564, https://doi.org/10.1016/0960-1686(92)90108-W, 1992.
Petroff, A. and Zhang, L.: Development and validation of a size-resolved particle dry deposition scheme for application in aerosol transport models, Geosci. Model Dev., 3, 753–769, https://doi.org/10.5194/gmd-3-753-2010, 2010.
Petroff, A., Mailliat, A., Amielh, M., and Anselmet, F.: Aerosol dry deposition on vegetative canopies. Part I: Review of present knowledge, Atmos. Environ., 42, 3625–3653, https://doi.org/10.1016/j.atmosenv.2007.09.043, 2008.
Platnick, S., et al.: MODIS Atmosphere L3 Daily Product, NASA MODIS Adaptive Processing System, Goddard Space Flight Center, USA, https://doi.org/10.5067/MODIS/MOD08_D3.006, 2015.
Prospero, J. M., Landing, W. M., and Schulz, M.: African dust deposition to Florida: Temporal and spatial variability and comparisons to models, J. Geophys. Res., 115, D13304, https://doi.org/10.1029/2009JD012773, 2010.
Pryor, S. C., Gallagher, M., Sievering, H., Larsen, S. E., Barthelmie, R. J., Birsan, F., Nemitz, E., Rinne, J., Kulmala, M., Grönholm, T., Taipale, R., and Vesala, T.: A review of measurement and modelling results of particle atmosphere–surface exchange, Tellus B, 60, 42–75, https://doi.org/10.1111/j.1600-0889.2007.00298.x, 2008.
Ruijrok, W., Davidson, C. I., and Nicholson, K. W.: Dry deposition of particles, Tellus B Chem. Phys. Meteorol., 47, 587–601, https://doi.org/10.3402/tellusb.v47i5.16074, 1995.
Seinfeld, J. H. and Pandis, S. N.: Atmospheric chemistry and physics: from air pollution to climate change, John Wiley & Sons, 2006.
Shao, Y.: A model for mineral dust emission, J. Geophys. Res.-Atmos., 106, 20239–20254, https://doi.org/10.1029/2001JD900171, 2001.
Shao, Y.: Simplification of a dust emission scheme and comparison with data, J. Geophys. Res.-Atmos., 109, 1–6, https://doi.org/10.1029/2003JD004372, 2004.
Shao, Y., Fink, A. H., and Klose, M.: Numerical simulation of a continental-scale Saharan dust event, J. Geophys. Res., 115, D13205, https://doi.org/10.1029/2009JD012678, 2010.
Shao, Y., Ishizuka, M., Mikami, M., and Leys, J. F.: Parameterization of size-resolved dust emission and validation with measurements, J. Geophys. Res., 116, D08203, https://doi.org/10.1029/2010JD014527, 2011.
Slinn, W. G.: Predictions for particle deposition to vegetative canopies, Atmos. Environ., 16, 1785–1794, https://doi.org/10.1016/0004-6981(82)90271-2, 1982.
Su, L. and Fung, J. C. H.: Sensitivities of WRF-Chem to dust emission schemes and land surface properties in simulating dust cycles during springtime over East Asia, J. Geophys. Res.-Atmos., 120, 11215–11230, https://doi.org/10.1002/2015JD023446, 2015.
Todd, M. C., Bou Karam, D., Cavazos, C., Bouet, C., Heinold, B., Baldasano, J. M., Cautenet, G., Koren, I., Perez, C., Solmon, F., Tegen, I., Tulet, P., Washington, R., and Zakey, A.: Quantifying uncertainty in estimates of mineral dust flux: An intercomparison of model performance over the Bodélé Depression, northern Chad, J. Geophys. Res., 113, D24107, https://doi.org/10.1029/2008JD010476, 2008.
Uno, I., Wang, Z., Chiba, M., Chun, Y. S., Gong, S. L., Hara, Y., Jung, E., Lee, S.-S., Liu, M., Mikami, M., Music, S., Nickovic, S., Satake, S., Shao, Y., Song, Z., Sugimoto, N., Tanaka, T., and Westphal, D. L.: Dust model intercomparison (DMIP) study over Asia: Overview, J. Geophys. Res., 111, D12213, https://doi.org/10.1029/2005JD006575, 2006.
Wesely, M. L., Cook, D. R., and Hart, R. L.: Measurements and parameterization of particulate sulfur dry deposition over grass, J. Geophys. Res., 90, 2131–2143, 1985.
WRF Users page: WRF Source Codes and Graphics Software Downloads, available at: http://www2.mmm.ucar.edu/wrf/users/download/get_source.html, last access: 23 April 2020.
Wu, C.-L. and Lin, Z.-H.: Uncertainty in Dust Budget over East Asia Simulated by WRF/Chem with Six Different Dust Emission Schemes, Atmos. Ocean. Sci. Lett., 6, 428–433, https://doi.org/10.1080/16742834.2013.11447120, 2013.
Wu, C.-L. and Lin, Z.-H.: Impact of two different dust emission schemes on the simulation of a severe dust storm in East Asia using the WRF/Chem model, Clim. Environ. Res., 19, 419–436, https://doi.org/10.3878/j.issn.1006-9585.2013.13041, 2014 (in Chinese).
Wu, M., Liu, X., Zhang, L., Wu, C., Lu, Z., Ma, P. L., Wang, H., Tilmes, S., Mahowald, N., Matsui, H., and Easter, R. C.: Impacts of Aerosol Dry Deposition on Black Carbon Spatial Distributions and Radiative Effects in the Community Atmosphere Model CAM5, J. Adv. Model. Earth Syst., 10, 1150–1171, https://doi.org/10.1029/2017MS001219, 2018.
Yuan, T., Chen, S., Huang, J., Zhang, X., Luo, Y., Ma, X., and Zhang, G.: Sensitivity of simulating a dust storm over Central Asia to different dust schemes using the WRF-Chem model, Atmos. Environ., 207, 16–29, https://doi.org/10.1016/j.atmosenv.2019.03.014, 2019.
Zaveri, R. A. and Peters, L. K.: A new lumped structure photochemical mechanism for large-scale applications, J. Geophys. Res.-Atmos., 104, 30387–30415, https://doi.org/10.1029/1999JD900876, 1999.
Zaveri, R. A., Easter, R. C., Fast, J. D., and Peters, L. K.: Model for Simulating Aerosol Interactions and Chemistry (MOSAIC), J. Geophys. Res.-Atmos., 113, 1–29, https://doi.org/10.1029/2007JD008782, 2008.
Zender, C. S., Miller, R. L. R. L., and Tegen, I.: Quantifying mineral dust mass budgets:Terminology, constraints, and current estimates, Eos, Trans. Am. Geophys. Union, 85, 509–512, https://doi.org/10.1029/2004EO480002, 2004.
Zhang, J. and Shao, Y.: A new parameterization of particle dry deposition over rough surfaces, Atmos. Chem. Phys., 14, 12429–12440, https://doi.org/10.5194/acp-14-12429-2014, 2014.
Zhang, L. and He, Z.: Technical Note: An empirical algorithm estimating dry deposition velocity of fine, coarse and giant particles, Atmos. Chem. Phys., 14, 3729–3737, https://doi.org/10.5194/acp-14-3729-2014, 2014.
Zhang, L., Gong, S., Padro, J., and Barrie, L.: A size-segregated particle dry deposition scheme for an atmospheric aerosol module, Atmos. Environ., 35, 549–560, https://doi.org/10.1016/S1352-2310(00)00326-5, 2001.
Zhang, L., Fang, G. C., Liu, C. K., Huang, Y. L., Huang, J. H., and Huang, C. S.: Dry deposition fluxes and deposition velocities of seven trace metal species at five sites in central Taiwan – a summary of surrogate surface measurements and a comparison with model estimations, Atmos. Chem. Phys., 12, 3405–3417, https://doi.org/10.5194/acp-12-3405-2012, 2012.
Zhang, X.-X., Sharratt, B., Liu, L.-Y., Wang, Z.-F., Pan, X.-L., Lei, J.-Q., Wu, S.-X., Huang, S.-Y., Guo, Y.-H., Li, J., Tang, X., Yang, T., Tian, Y., Chen, X.-S., Hao, J.-Q., Zheng, H.-T., Yang, Y.-Y., and Lyu, Y.-L.: East Asian dust storm in May 2017: observations, modelling, and its influence on the Asia-Pacific region, Atmos. Chem. Phys., 18, 8353–8371, https://doi.org/10.5194/acp-18-8353-2018, 2018.
Zhang, X. X., Sharratt, B., Lei, J. Q., Wu, C. L., Zhang, J., Zhao, C., Wang, Z. F., Wu, S. X., Li, S. Y., Liu, L. Y., Huang, S. Y., Guo, Y. H., Mao, R., Li, J., Tang, X., and Hao, J. Q.: Parameterization schemes on dust deposition in northwest China: Model validation and implications for the global dust cycle, Atmos. Environ., 209, 1–13, https://doi.org/10.1016/j.atmosenv.2019.04.017, 2019.
Zhao, C., Liu, X., Leung, L. R., Johnson, B., McFarlane, S. A., Gustafson Jr., W. I., Fast, J. D., and Easter, R.: The spatial distribution of mineral dust and its shortwave radiative forcing over North Africa: modeling sensitivities to dust emissions and aerosol size treatments, Atmos. Chem. Phys., 10, 8821–8838, https://doi.org/10.5194/acp-10-8821-2010, 2010.
Zhao, C., Liu, X., Ruby Leung, L., and Hagos, S.: Radiative impact of mineral dust on monsoon precipitation variability over West Africa, Atmos. Chem. Phys., 11, 1879–1893, https://doi.org/10.5194/acp-11-1879-2011, 2011.
Zhao, C., Liu, X., and Leung, L. R.: Impact of the Desert dust on the summer monsoon system over Southwestern North America, Atmos. Chem. Phys., 12, 3717–3731, https://doi.org/10.5194/acp-12-3717-2012, 2012.
Zhao, C., Chen, S., Leung, L. R., Qian, Y., Kok, J. F., Zaveri, R. A., and Huang, J.: Uncertainty in modeling dust mass balance and radiative forcing from size parameterization, Atmos. Chem. Phys., 13, 10733–10753, https://doi.org/10.5194/acp-13-10733-2013, 2013.
Zhao, T. L., Gong, S. L., Zhang, X. Y., Abdel-Mawgoud, A., and Shao, Y. P.: An assessment of dust emission schemes in modeling east Asian dust storms, J. Geophys. Res.-Atmos., 111, 1–11, https://doi.org/10.1029/2004JD005746, 2006.