Numerical Studies on a Severe Dust Storm in East Asia Using WRF-Chem ()
1. Introduction
East Asia has been a region with frequent dust storm which is considered to be a seasonal meteorological phenomenon. Usually, the dust originates in the deserts in Mongolia, China and Kazakhstan with high-speed surface winds. There are three key source regions of deserting China, which are Taklamakan in Xinjiang, Hunshandak in Inner Mongolia and the Loess Plateau. The intense dust storms can kick up dense clouds of fine, dry soil particles which can be carried eastward by prevailing winds and pass over China. As one of the severe disaster weather, dust storm always occurs in North China during the springtime. It can do harm to industrial and agricultural activities, transportation, climate change, atmospheric environment and human health [1] [2] [3] [4] .
Many studies have been focused on dust storm in recent years. Especially, numerical models are widely applied to investigate the physical mechanism and long range transport of dust aerosols. Dust emission and transport studies have been undertaken since 1980s [5] [6] [7] . Gillette and Hanson (1989) combined atmospheric data with land surface data for wind erosion assessment to investigate the spatial and temporal variations of dust production in the United States [8] . Marticorena and Bergametti (1995), Shao et al. (1996) and Marticorena et al. (1997) have developed wind erosion schemes accounting for the impact of surface properties on sand drift and dust emission [9] [10] [11] . Shao et al. (1997) and Lu et al. (2001) have developed and implemented an almost fully integrated wind erosion modeling and prediction system [12] [13] . Alfaro and Gomes (2001) performed a modeling work on mineral aerosol production by wind erosion, focusing on emission intensities and aerosol size distributions in source areas [14] . Song et al. (2004) conducted a dust simulation and prediction in China during the spring of 2002 using numerical model [15] . Gong et al. (2008) reported an integrated sand and dust storm forecasting system CUACE-Dust (Chinese Unified Atmospheric Chemistry Environment for Dust) [16] . Based on the 1991-2010 simulation of dust aerosol emissions with the global climate model CAM5.1, Liu et al. (2015) analyzed the globally spatial and temporal variations of dust aerosol emissions over major deserts [17] . To ensure the accuracy of simulation, the parameterization scheme should be considered all of the uncertain factor: wind speed, surface properties. Gherboudj et al. (2015) showed that soil moisture and soil texture have significant roles in the dust emissions over the Middle East and North Africa [18] . However, few study has been conducted on dust storm simulation in China using the new generation weather research, forecast and chemistry model WRF-Chem. Wu et al. (2013) had clarified the uncertainty in dust budget over East Asia simulated by WRF-Chem with six different dust emission schemes [19] . By using WRF-Chem, Liu et al. (2016) gave the radiative forcing of dust aerosol and the radiative feedbacks on the planetary boundary layer in North China during a typical Asian dust storm [20] .
In Mar., 2010, starting from Mongolia, a dust storm occurred with widely impacting area and heavy air pollution. This dust storm event has close relations with the air temperature increases and strong wind prevailing during the development of Mongolia cyclone. Since Spring of 2010, the air temperature increases gradually, the soil is unfrozen, the surface heating becomes fast, and resulting in strong instability in the whole atmosphere and strong wind with the intersection of warm and cold air,which provide dynamical conditions for the breakout of dust storm. In addition, the Mongolia cyclone is strongly developed and located near source regions of dust, leading to strong upward movement of the air which is beneficial for the occurrence of dust storm.
In this paper, the characteristics of the severe dust storm were investigated by both observational and numerical studies. The models and data are descripted in Section 2. Surface observation of PM10 concentration and aerosol optical parameters from AERONET (Aerosol Robotic Network) are analyzed. Satellite data from MODIS (Moderate Resolution Imaging Spectroradiometer) and CALIPSO (Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observations) are analyzed. Modeling studies were performed using WRF-Chem (Whether Research and Forecasting model with Chemistry) in consideration of two different dust schemes as proposed by Shaw (2008) and Chin (2002). The model results were compared with those from CUACE-Dust and RegCCMS, which are models developed by Chinese scientists and had been applied for dust simulation. These results are presented and discussed in Section 3. The final conclusions and future perspective are given in Section 4.
2. Models and Data
2.1. Models
2.1.1. WRF-Chem
WRF-Chem is a new generation model where the air quality component (Chem) and meteorological component (WRF) are fully coupled “online”. The meteorological model WRF is a next-generation mesoscale numerical weather prediction system designed to serve both operational forecasting and atmospheric research needs. It is suitable for use in a broad spectrum of applications across scales ranging from meters to thousands of kilometers. The Chem model contains relatively comprehensive physical and chemical processes, includes convective, adjective and turbulent chemical transport, biogenic and anthropogenic emissions, deposition (dry and wet), photolysis, gas and aqueous chemical transformation, and aerosol chemistry and dynamics (including inorganic and organic aerosols) as well as the processes of aerosol direct and indirect effect through interaction with atmospheric radiation, photolysis, and microphysics.
In WRF-Chem version 3.2, two dust production schemes are implemented in the model so that it can be applied on dust simulations. Shaw (2008) used a dust generating module and the land cover data sets to describe how particles will be emitted with winds, which are included in the Dust Regional Atmospheric Model (DREAM) [21] . DREAM can be configured for multiple particle size bins. Chin (2002) adopted a new approach which has been used in the Georgia Institute of Technology-Goddard Global Ozone Chemistry Aerosol Radiation and Transport (GOCART) model to identify the dust source [22] . Dust particles ranging from 0.1 to 6 mm in radius are considered in the model. The emission flux is expressed as functions of a dimensional factor, probability-source function, fraction of size group within the soil, 10-m wind speed, and the threshold velocity of wind erosion determined by particle size and surface wetness. For more detail, previous works from Shaw et al. (2008) and Chin et al. (2002) can be referred to.
2.1.2. CUACE-Dust
CUACE-Dust (Chinese Unified Atmospheric Chemistry Environment for Dust) is an integrated sand and dust storm forecasting system which consists of a comprehensive dust aerosol module with emission, dry/wet depositions and other atmospheric dynamic processes, and a data assimilation system using observational data from the CMA (China Meteorological Administration) ground dust monitoring network and retrieved dust information from a Chinese geostationary satellite-FY-2C (Gong et al., 2008; Niu et al., 2007) [16] [23] . CUACE-Dust has been an operational model for dust forecast. Here, the model is applied for comparison studies with WRF-Chem.
2.1.3. RegCCMS
RegCCMS is a Regional Climate and Chemistry Modeling System developed by Nanjing University. It was designed for studies on the climatic effect of different aerosols, including sulfate, nitrate, black carbon, organic carbon, sea salt as well as dust (Wang et al., 2009, 2010; Zhuang et al., 2010; Li et al., 2009) [24] [25] [26] [27] . The dust scheme applied in RegCCMS is referred to Marticorena and Bergametii (1995)’s work [9] . The coarse, fine and iken modes were treated in three bins. Heterogeneous reactions on dust particles were considered as described by Deng et al. (2011) [28] . The indirect effect of dust aerosols on regional climate over East Asia was investigated by Shen et al. (2011) on the base of RegCCMS [29] .
2.2. Observations/Measurements
2.2.1. Surface Data
The Air Pollution Index (API) data as well as PM10 surface concentration data from 86 stations were collected to analyze the spatial distribution of air pollution and long range transport of dust aerosols during Mar. 17-24,2010. API is a number to describe daily air quality, which is calculated by the observed surface concentrations of three air pollutants (SO2, NO2, PM10) based on Chinese Air Quality Standard. There are five grades for API. They are <50, 50 - 100, 100 - 200, 200 - 300, and >300. API <50 and 50 - 100 represents for the good and better air quality, while the API >100 means that the air is polluted in different levels. The API and daily averaged PM10 data is reported by Chinese Environment Monitoring Station (http://www.cnemc.cn).
AERONET data at three sites of Xianghe (Beijing), NUIST (Nanjing), and HokTsui (Hong Kong) were collected (in Figure 1), including aerosol optical depth (AOD) at 340, 380, 440, 500, 675, 870, 1020, 1640 nm. We use fine mode AOD at 500 nm, coarse mode AOD at 500 nm, fine mode fraction at 500 nm and angstrom parameter. The optical characteristics of dust aerosols can be revealed from the surface-based observations. All these data can be download at http://aeronet.gsfc.nasa.gov/new_web/data.html. It’s worth noting that AERONET data represent the state of the overall atmospheric column.
Figure 1. Three AERONET stations in China.
2.2.2. Satellite Data
MODIS AOD (550 nm) data were selected for investigation on geographical/spatial patterns of dust aerosol. Global daily Level 3 data are used from the Terra platform with the spatial resolution of 1 × 1 degree. Although the data only covers a few part of China in space and time, it does provide some useful information on how aerosols are transported during the different stage of dust storm. The data set can be download at http://disc.sci.gsfc.nasa.gov/giovanni/ and analyzed and visualized with the Giovanni online data system, developed and maintained by the NASA GES DISC(Acker et al., 2007) [30] .
CALIPSO (Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observations) Lidar Level 2 Version 3.01 was collected for this study. It includes backscatter coefficients, depolarization ratio profile and extinction coefficients. Extinction coefficient was applied for investigation vertical optical characteristics of dust aerosol. CALIPSO data are provided by NASA Langley Research Center Atmospheric Science Data Center. Shun et al. (2016) give the statistical characteristics of aerosol extinction coefficient profile in East Asia from CALIPSO [31] .
3. Results and Discussions
3.1. Observations
3.1.1. API and PM10
As affected by dust storm, the air was heavily contaminated with more particles. The API is an effective index to describe how the air is polluted. In China, API is divided into five levels (<50: excellent; 50 - 100: good; 100 - 150: slightly polluted; 150 - 200: low grade polluted; 200 - 250: medium grade polluted; 250 - 300: medium to heavily polluted; >300: heavily polluted). During Mar. 19-23, 2010, the severe dust storm occurred in China. Figure 2 shows the evolution of the API at 86 stations during the event. Dust started from west and south Mongolia in the afternoon of Mar. 18 as influenced by
Mongolia cyclone and cold air. Dust storm hit China from Mar. 19. About 16 provinces or municipals, including Xinjiang, Inner Mongolia, Qinghai, Gansu, Ningxia, Shanxi, Shanxi, Hebei, Beijing, Tianjin, Shandong, Henan, Jiangsu, Anhui, Hubei, Sichuan, were suffered from this dust event on Mar. 20. API increase significantly and abruptly to 300 and exceeded 500 in many northern cities. The API maximum shifted to the southeast area, suggesting the transport of dust aerosols eastward and southward. The dust is transported from source regions to coastal areas, including Jiangsu, Zhejiang, Hubei, Anhui and Shanghai, on the day of Mar. 21. API remains high at North and East China with maximum of 500. In Mar. 22 and 23, dust aerosol continued to be transported to southeast China, far reaching Taiwan area. API of more than 500 in Zhejiang, Fujian and Guangdong provinces were observed. API higher than 400 in Hong Kong was recorded, which was never reported before. On Mar. 24, API at all stations decrease below 200, showing the weakness of dust transport in continental China. The spatial distribution of API during Mar. 18 - 23 show the spread of dust to southeast with heavy pollution induced. This dust event has confirmed to be the most widely-spread and intense one since 2009.
During the episode of dust storm, the primary pollutant in the air is inhalable particulate named as PM10. Figure 3 illustrates the daily average surface concentration of PM10 at 86 stations operated by China Environment Monitoring Center. Figure 3 shows that the maximum of surface hourly PM10 concentration for each site ranges from 500 to 2000 μg/m3. PM10 of 944 μg/m3 and 1990 μg/m3 in Shanghai and Nanjing were observed. Staring from the noon of Mar. 21, hourly PM10 concentration increased rapidly in Hong Kong with maximum of 783 μg/m3. The maximum hourly PM10 level of 1724 μg/m3 in Taiwan was also recorded. PM2.5/PM10 mass concentration ratio accounted for 14% and 3% in Shanghai and Nanjing, showing the domination of coarse particles affected by dust aerosols.
3.1.2. Aerosol Optical Parameters
The optical parameters observed at Xianghe, Beijing on Mar. 21 were illustrated in Figure 4. The dust storm hit Beijing from Mar. 20 and intensified on Mar. 21. High AOD (2.5) was measured with low fine mode fraction (0.2) and small angstrom parameter (~0.5). According to Tanre et al. (2001), angstrom parameter ranges between −1 to 0.5 for dust aerosol [32] . In general condition, the urban air is polluted due to emissions from industry and transportation, which is characterized with high fine mode fraction and large angstrom parameter due to the domination of fine particles which induce much strong scattering. Angstrom parameter is between 1.1 to 2.4 for industrial aerosol (Dubovik et al., 2001) [33] . The observed low angstrom parameter suggested that coarse particles are dominated during the dust event.
Figure 5 depicts the optical parameters observed at HokTsui, Hong Kong. It shows that AOD remains high on Mar. 22. Coarse particle is dominated in total AOD, and fine mode fraction only account for 20% - 40%. AOD (500 nm) varied at a mean level of 0.8 and the peak can be more than 1. Most of the Spectral Deconvolution Algorithm (SDA) is below 1 Angstrom parameter (440 - 870 nm) is below 0.7 and as low as 0.2, suggesting the weak scattering due to domination of coarse particles. The observed optical parameters confirmed that Hong Kong was hit by the dust storm, which has
(a)(b)
Figure 3. (a) PM10 surface concentration and (b) MODIS AOD at 86 cities.
never occurred before in the Hong Kong history.
3.1.3. MODIS AOD Data
The spatial/geographical distributions of AOD (550 nm) from MODIS during Mar. 18- 23 were presented in Figure 6. Although MODIS does not give retrieval AOD value in most area in northwest China, variation of AOD in east China can provide evidence of long range transport of dust form northwest to southeast. On Mar. 21, much high AOD was observed in central and southeast China, extending to East Ocean. Coverage area with high AOD decreases in the following two days. Considering mass concentration, the coarse particles contribute much in the total aerosols. However, both coarse and fine particles have large numbers during the dust days, resulting in strong aerosol extinction coefficient and large AOD with maximum of more than 2.5.
Figure 7 is the time series of AOD averaged over the region covering 70E - 140E, 10N - 60N. As expected, AOD shows gradual increase before Mar. 21, reaching 0.6 on the day of Mar. 21, and then decrease after Mar. 21. The aerosol small mode optical depth remains small and only accounts for 41% of total AOD, showing the dominance of coarse particles. The mass concentration ratio of PM2.5/PM10 in urban area is about 60% - 70% in non dust days. Angstrom parameter also exhibits lower value of 1, suggesting that the contribution of coarse particle is significant. On Mar. 22, although
Figure 6. MODIS AOD (550 nm) distribution during Mar. 18-23.
AOD decrease, small mode fraction as well as angstrom parameter increase, indicating the importance of fine particles when the impact of dust became weaker.
Figure 8 presents time-longitude cross section of aerosol optical parameters. From Figure 8(a), high AOD center move from mid-latitude (40N - 45N) to low latitude (20N - 25N), and then back to mid-latitude (35N - 40N), showing the transport of dust
Figure 7. Area averaged time series of optical parameters during Mar. 17-24 (a) AOD, (b) Aerosol fraction mode fraction, (c) Angstrom parameter.
from north to south. On the same days, aerosol small mode optical depth decrease with the increase of latitude (Figure 8(b)), showing more fine particles in the south and more coarse particles in the north. It reaches below 0.2 - 0.4 at 45N and 0.6 - 0.8 at 10N, suggesting the coarse particles in the north and fine particles in the south. Angstrom parameter shows consistent pattern with aerosol small mode optical depth (Figure 8(c)).
3.1.4. Extinction Coefficient Vertical Profile
The CALIPSO is the last in a sequence of several level 1 and level 2 algorithms that retrieve profiles of aerosol backscatter and extinction coefficients [34] . The profiles of mean extinction coefficient at 532 nm during Mar. 20 and 21 are illustrated in Figure 9. It shows that the dust can be transported up to 5 Km in vertical direction. The extinction coefficient maximum of 0.35 - 0.4 was observed at 1 - 3 Km, suggesting that strong
Figure 8. Time-longitude cross section of aerosol optical parameters (a) AOD, (b) Aerosol fraction mode fraction, (c) Angstrom parameter.
convection can lead to transport of dust aerosol to the higher level over the top of boundary layer.
3.1.5. Synopticweather
Depicts the weather chart (http://218.94.36.199:5050/dmsg/map.htm) on Mar. 19-20. The Mongolia cyclone appeared on Mar. 17. A warm center of 21 in the northwest was found at 850 hPa on Mar. 17 and 18. On the other hand, the cold high pressure generated near lake Balkhash on surface weather map. On Mar. 19, strong cold center appeared at the Caspian sea to lake Balkhash where north branch jet occurred. At 500 hPa, the north branch trough in Mongolia and Hetao area split, cold high pressure at the surface is located at north Xinjiang and west Mongolia with central pressure greater than 1040 hPa. At the mean time, warm low pressure was developing in the east of Mongolia. The extremely unstable stratification condition was formed, which is bene- ficial to the exchange of cold and warm air in upper and lower levels. The strong
Figure 9. Profile of CALIPSO mean extinction at 532 nm and lidar moving path on Mar. 20 ((a), (b)) and Mar. 21((c), (d)).
intersection of cold and warm front results in the great gradient of temperature and pressure, surface wind speed becomes very strong, resulting in breakout of dust storm in Mongolia and central/west area of Inner Mongolia. Accompanying by the eastward of the synoptic system, dust can be transported to the south and east. This cold front is one of the important synoptic weather systems which are favorable for formation of dust storm [35] .
3.1.6. Backward Trajectory
Three day backward trajectories arriving at Beijing, Nanjing and Hong Kong were calculated to explore the movement of the air mass. HYSPLIT4 was used to calculate backward trajectory with meteorological data from GDAS (Global Data Assimilation System) from NCEP (National Centers for Environmental Prediction). As shown in Figure 10, the air mass was originated from northwest, passing by Xinjiang, Inner
Figure 10. 72 h backward trajectories to Beijing, Nanjing and Hong Kong on (a) Mar. 20 and (b) Mar. 21.
Mongolia, Gansu, Hebei etc., and then reaching Beijing and Nanjing which were hit by the dust storm in Mar. 20 and Mar. 21, respectively. The dust storm affected Hong Kong in the Mar. 22. Three source regions are responsible of this strong dust event. Grini et al. (2005) had given the effects of soil erodibility and wind speed variability to the dust transport for the global atmosphere [36] . They are sandland in west and south of Mongolia, desert and sandland in west Inner Mongolia and its surrounding area, Hunshandak sandland in central Inner Mongolia. The calculated backward trajectories clearly exhibited the originations and transport path of dust aerosols.
3.2. Numerical Simulations
3.2.1. WRF-Chem Simulation
For numerical simulation of the dust event, two nesting domains were defined using Lambert projection. The large domain has grid spacing of 81 km with horizontal grids of 88 × 75, and the small domain has grid spacing of 27 km with horizontal grids of 85 × 70. The model top pressure is set as 50 hPa, with 27 sigma layers in vertical. We used NCEP global reanalysis data with 1˚ × 1˚ resolution as initial conditions and boundary conditions, the later were forced every six-hour. The ideal concentration profiles were used as initial and boundary conditions of chemical species. The RADM2 gas chemistry was chosen. MADE/SORGAM and GOCART aerosol module were selected for Shaw and Chin scheme, respectively. Other choices for physical processes are similar to Jiang (2008)’s work [37] .
WRF-Chem (V3.2) model was integrated for 7 days staring from Mar. 17. Figure 11 illustrated the surface PM10 concentration simulated using two different dust schemes. Both Shaw and Chin schemes show similar patterns of surface PM10 concentration and dust moving path. However, Shaw scheme predicted larger coverage of dust aerosol and stronger dust emission in Mongolia compared to Chin scheme.
3.2.2. Comparison with Observational Data
The simulated daily average surface concentrations of PM10 at 86 stations were compared with the observations. These results are shown in Figure 12. Both schemes over predicted the PM10 concentration to certain extent. The PM10 coorelation coefficient between the model and the observation is 0.41 and 0.56 for Shaw scheme and Chin Scheme, respectively, with relative error of 60% and 57%, respectively. It shows that Chin (2002) scheme is slightly better than Shaw (2008) scheme if taking PM10 concentration into account. There may be several reasons to explain the difference. First, the wind erosion scheme is very important in dust simulation, which determines how much dust was emitted into the air. Second, meteorology simulation play crucial role in dust transport and deposition, which determine where the dust can be transported. The input parameters, such as landuse, landcover and soil erosion data, become the third factors on the estimation of dust production [38] . The three aspects need to be improved to get more satisfied results.
3.2.3. Comparison with CUACE-Dust and RegCCMS
Figure 13 presents comparison of the simulated surface concentration of PM10 between CUACE-Dust and RegCCMS. Both of the model show similar patterns of surface PM10
Figure 11. PM10 surface daily mean concentration Shaw ((a), (b), (c)) and Chin ((d), (e), (f)) in WRF-Chem during Mar. 20-22 (μg/m3).
Figure 12. Daily PM10 surface concentration from observation and WRF-Chem for Shaw scheme and Chin scheme (μg/m3).
Figure 13. PM10 surface concentration by CUACE-Dust and RegCCMS during Mar. 20-22 (μg/m3).
concentration with WRF-Chem over dust source regions as well as downstream area, especially in Xinjiang and Mongolia. This dust storm is primarily originated from Mongolia, however there is also sand dust coming from Xinjiang. In general, CUACE- Dust predicted high PM10 concentration with narrow coverage, while PM10 level is underestimated and dust is widely spread and transported to Southeast China for RegCCMS. Further, dust emission in Mongolia is better estimated for RegCCMS than CUACE-Dust in consideration of the spatial distributions. The total dust of 120 Mt was roughly estimated by RegCCM, in which 64% was deposited and 36% was stayed in the air for this dust storm episode.
4. Conclusions
Dust storm has become one of the important sources of severe air pollution episode in China, especially during spring season. A serious dust storm formed in Mongolia on Mar. 19, 2010. Then it was transported to the northern China, and inhalable particles increased in Gansu, Inner Mongolia, Beijing and Shandong provinces. The dust then spread to the south of Yangtze River region two days later. The area with large API shifted from Northwest to Southeast. Dust aerosols were transported to Hubei, Anhui, Jiangsu and Zhejiang provinces in Central and East China on Mar. 21 and then Guangdong, Fujian provinces and Hong Kong in South China on Mar. 22. Air pollution induced by dust storm was alleviated in coastal regions and other areas of China since Mar. 23 due to the end of long range transport of dust aerosol. This serious dust event has close relations to the increasing temperature, strong wind and the development of Mongolia cyclone.
The Version 3.2 of WRF-Chem was applied to simulate this dust storm case. The surface observations of PM10 concentration and AOD from MODIS were used to validate the model performance. It shows that WRF-Chem has the ability to simulate the production, transportation and deposition of the dust storm. Considering the two dust schemes in WRF-Chem, Chin (2002) scheme presents much better results in terms of PM10 surface concentration simulation. However, Shaw (2008) scheme gives more reasonable spatial distribution of dust aerosols. WRF-Chem as well as CUACE-Dust and RegCCMS show performance in simulation of surface PM10 concentration and transport path of dust aerosols. Although large amount of dust was emitted into the air, most of which was deposited during the long range transport to downstream regions.
The studies in this paper suggested that estimation of sand emission is crucial in numerical simulation of dust intensity, which can result in different PM10 and AOD pattern in terms of airborne particulate. WRF-Chem is assessed to be a useful tool in simulation of dust event. However, in order to get more realistic results, assimilation of different aerosol data into the model are highly recommended to improve the simulation on concentration and optical parameters. This work will be addressed in another paper.
Acknowledgements
We acknowledge the MODIS mission scientists and associated NASA personnel for the production of the data used in this research effort. Analyses and visualizations on MODIS data used in this study was produced with the Giovanni online data system, developed and maintained by the NASA GES DISC. The CALIPSO data were obtained from the NASA Langley Research Center Atmospheric Science Data Center. AERONET Principal Investigators are thanked for making their data available. This work was supported by the State Key Research and Development Program (2016YFC0208504, 2016YFC0203303) and National Natural Science Foundation of China (91544230, 41575145 and 41621005).