Source Allocation of Long-Range Asian Dusts Transportation across the Taiwan Strait by Innovative Chemical-Assisted Identification Methods

Abstract

This study used the backward trajectory calculation to obtain the transportation routes of Asian dusts and further combined the chemical composition with the enrichment factor (EF) and the grey relational analysis (GR) to identify the potential sources of eighteen Asian dust storm (ADS) events. The results showed that the chemical compositions of atmospheric particles sampled at the Pescadores Islands were very similar to source soils fugitively emitted from Inner Mongolia, which could assist in identifying the source regions of Asian dusts. This study further compared the source allocation of Asian dusts obtained from EF, GR, and backward trajectory, which showed that the source regions of Asian dusts obtained from these three methods were quite similar. The similarity of backward trajectory and GR reached as high as 83.3%. Moreover, the similarity of backward trajectory calculation and EF or GR was up to 77.8% while that of the GR and EF was up to 83.3%. Overall, these three methods can successfully allocate the source regions of Asian dusts by 66.7%. Moreover, these innovative chemical-assisted methods can be successfully applied to identify the source regions of Asian dusts for 18 ADS events.

(ProQuest: ... denotes non-US-ASCII text omitted.)

Yi-Hsiu Jen 1 and Yi-Chi Liu 1 and Iau-Ren Ie 1 and Chung-Shin Yuan 1 and Chung-Hsuang Hung 2

Academic Editor:Hesham El-Askary

1, Institute of Environmental Engineering, National Sun Yat-Sen University, Kaohsiung 80424, Taiwan
2, Department of Safety, Health and Environmental Engineering, National Kaohsiung First University, Kaohsiung 824, Taiwan

Received 13 June 2014; Accepted 8 August 2014; 27 August 2014

This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

1. Introduction

At the beginning (February-May) of each year, Asian dusts are frequently blown from hyperarid, arid, and semiarid lands located at Northwest and North China, South Mongolia, and South Siberia [1-3]. For the past half-century, the frequency of Asian dust storms (ADS) occurring in China increased significantly from 0.5 (1950') to 14 (2000) times per year, especially the strong and/or super ADS. This phenomenon is correlated mainly to climate change both locally and globally [4, 5]. It is estimated that approximately 800 million metric tons of sands were emitted to the atmosphere by ADS, which accounts for almost half of natural fugitive sands globally [6, 7].

ADS could transport in the upper troposphere easterly to Korea, Japan, and even North America and may also transport southerly to Okinawa, Taiwan, Hong Kong, and Philippine due to anticyclone of atmospheric circulation [8-11]. Although many sands might be settled down to the ocean and provided nutrients (e.g., Fe²⁺ ) to the oceanic biosphere, which actually benefits the circulation of global geochemistry, ADS could deteriorate ambient air quality, reduce atmospheric visibility, cause adverse effects on human health, and result in poor scenery and public welfare [12-14]. Additionally, great losses of national economy due to ADS are constantly reported in the West Pacific regions, including crop damage, soil productivity losses, livestock losses, mass migration, human health impacts, and impacts on climate [15-18].

In general, ADS invades Taiwan 2-3 times per year and lasts for about 1-2 days for each invasion. The concentration of suspended particles (PM₁₀ ) in the events of ADS increases approximately 2-3-fold when compared to normal days [7, 19, 20]. During the ADS periods, significant increases of atmospheric aerosols, particularly PM₁₀ , have been frequently observed in Taiwan [7, 8, 19, 21]. Particularly, yellow rain episodes have been reported both in Japan and Taiwan at the end of March 2000 [22]. The morphology and apparent color of dust particles filtered from the yellow rain were similar to yellow sands at their source [22]. Furthermore, the results obtained from the investigation on source apportioning atmospheric aerosols at an Asian dust episode in metro Kaohsiung indicated that the percentage of fugitive dusts increased significantly from 8.5% (regular periods) to 26.2% (Asian dust periods) [23]. Moreover, during the transportation of ADS across continents, the physicochemical characteristics of Asian dusts might be changed by mixing with anthropogenic particles emitted from typical stationary and mobile sources [24, 25]. Yuan et al. [19] reported that the size distribution of atmospheric aerosols sampled at the Pescadores Islands located in the middle of the Taiwan Strait during the period of ADS switched from bimode to single mode with a peak in coarse particles. Besides, the sulfate ( ^{_{SO 4} 2 -} ) and ammonia ( ^{_{NH 4} +} ) of coarse particles (PM_2.5-10 ) increase significantly during the period of ADS when compared to non-ADS periods [7, 26, 27].

In order to investigate the feasibility of chemical-assisted ADS source identification methods which can be used to trace back the soil source regions of ADS, this study compared the transportation routes of ADS event from the backward trajectory calculation to the chemical composition of atmospheric particles and soil sources. Top soils collected in three soil source regions of Inner Mongolia were resuspended and sampled to analyze their chemical composition for further conducting the enrichment factor analysis (EF) and grey relational analysis (GR). Ultimately, this study revealed and discussed the similarity percentages from the backward trajectory calculation, EF, and GR and further obtained the source regions of Asian dusts for 18 ADS events.

2. Methodologies

2.1. Sampling Sites of Atmospheric Particles and Soil Dusts

Located at the middle of Taiwan Strait between southeastern China and Taiwan Island, the Pescadores Islands have an area of 127.97 km² , which has the subtropical weather being mainly influenced by East Asian monsoons. During the 18 ADS events, the atmospheric particles were sampled at Hsiao-men site (23°38[variant prime]471[variant prime][variant prime] N; 119°30[variant prime]316[variant prime][variant prime] E), which is located at the northwestern coastline of the Pescadores Islands and is approximately 12 m above the ground and 500 m and 50 m far from the coastline and the major roads, respectively. At the sampling sites, this study used to simultaneously collect PM₁₀ (including PM_2.5-10 and PM_2.5 ) with dichotomous samplers (Anderson, Model Series 241). The sampling flow rate of the sampler was 16.7 L/min for a total of 24 hours, and the filters used in this study were 37 mm quartz filters.

The soil samples were collected from three ADS source regions in the western (38°13[variant prime]19 to 39°06[variant prime]45 N; 115°21[variant prime]47 to 116°53[variant prime]38 E), central (40°06[variant prime]49 to 41°39[variant prime]0 N; 110°28[variant prime]35 to 111°32[variant prime]20 E), and eastern (43°03[variant prime]27 to 43°28[variant prime]37 N; 116°41[variant prime]41 to 119°02[variant prime]09 E) Inner Mongolia, China [7]. Among them, strong winds and dry climate result in serious soil erosion and thus cause heavy ADS frequently (Inner Mongolia Web Site, 2003). A total of thirty soil samples were taken at the top 10-15 cm surface of open sand lands and were then separately stored in labeled sampling bags. The soil samples were then sent back to the Air Pollution Laboratory at the Institute of Environmental Engineering in National Sun Yat-Sen University, Taiwan, for further chemical analyses. The locations of all soil sampling sites are illustrated in Figure 1.

Figure 1: Location of atmospheric particles (*) at the Pescadores Islands and soil source regions ([white circle]) in the Inner Mongolia.

[figure omitted; refer to PDF]

2.2. Backward Trajectory Calculation

In order to trace the transportation of air parcel, backward trajectories from a specific receptor site were commonly used in this study, which could be used to identify the air mass transportation routes from specific source regions. A hybrid single-particle Lagrangian integrated trajectory (HYSPLIT) has been a widely used model that plots the trajectory of a single air parcel from a specific location and height above ground over a period of time [28, 29]. Backward trajectory calculation used the interpolated measured or modeled meteorological fields to estimate the most likely route over the geographical areas that provided air parcels from its source to a receptor at a given time. The method essentially follows a parcel of air backward in hourly steps for a specified length of time [10, 30]. Moreover, HYSPLIT uses archived 3-dimensional meteorological fields generated from observations and short-term forecasts. HYSPLIT can be run to generate either forward or backward trajectories using several available meteorological data archives, which is developed by the National Oceanic and Atmospheric Administration's (NOAA) Air Resources Laboratory (ARL). All backward trajectories were drawn at an arrival height of 50 m above the sea level in this study. By using the NOAA-HYSPLIT model, the dates with high particle concentration at the Pescadores Islands were determined and the transportation routes of air parcels toward the Pescadores Islands during the sampling periods were then simulated. This information was further applied to examine the possible long-range transportation routes conveying ADS from several long-range soil source regions to the Pescadores Islands.

2.3. Enrichment Factor Analysis

Soil samples collected from the eastern, central, and western parts of the Inner Mongolia were sent back to the Air Pollution Laboratory at National Sun Yat-Sen University for further soil resuspension and chemical analyses. The soil samples were resuspended with a dry powder atomizer at the top of a resuspension chamber (2 m × 1 m² ), which is illustrated in Figure 2. The resuspended dusts were then sampled as PM_2.5-10 and PM_2.5 for a total of 0.5 hours by two dichotomous samplers with a total volumetric flow rate of 16.7 L/min, which was located at the bottom of the resuspension chamber. Particulate matter with diameters larger than 10 μ m tended to settle down at the bottom of the resuspension chamber. Besides, the resuspended soil dusts were further characterized chemically to apply for the identification of Asian dusts at the downwind sites by using an enrichment factor (EF). The enrichment factor (EF) defined by (1) was used to correlate the downwind atmospheric particles sampled at the Pescadores Islands, Taiwan, to ADS originated from Inner Mongolia, China: [figure omitted; refer to PDF] where EF is the enrichment factor; Tr is the trace element; Ref is the reference element; and subscripts PM and crust represents atmospheric particles and crustal soils. Among them, the reference element has to be initially determined. Previous studies reported that the mass ratio of Mg, K, and Fe (or Al) to Al (or Fe) showed a relatively narrow range for coarse particles, but not for fine particles, which can be used as the fingerprints of Asian dusts originated from the Inner Mongolia [7]. Generally, the enrichment factor close to unity (EF [approximate] 1) indicates that the elements consider very high possibility from the specific sources. It is also defined that the elements might come from the specific source when the EF ranged from 0.1 to 10.

Figure 2: Schematic diagram of the resuspension chamber.

[figure omitted; refer to PDF]

2.4. Grey Relational Analysis

The technique of grey relational analysis is suitable for unascertained problems with poor information [31]. It provides an efficient solution to the uncertainty, multi-input, and discrete data problem. This kind of interaction is mainly through the connection among parameters and some conditions that are already known. Also, it will indicate the relational degree between two sequences with the help of the grey relational analysis. Moreover, the grey relational grade could utilize the discrete measurement method to judge the distance.

When the range of the sequence is too large or the standard value is too enormous, it will cause the influence of some factors to be neglected. Also, in the sequence, if the factors' goals and directions are different, the grey relational analysis might also produce unexpected incorrect results. Therefore, preprocessing of all the data is highly required. This process is called the grey relational generating. An ideal sequence is _{x 0} ( k ) ( k = 1,2 , 3 ) for three responses. The definition of the grey relational grade in the grey relational analysis is to show the relational degree between n sequences ( _{x 0} ( k ) and _{x i} ( k ) , i = 1,2 , ... , n ; k = 1,2 , 3 , ... , n ). The grey relational coefficient _{ξ i} ( k ) could be calculated by the following equation: [figure omitted; refer to PDF] where _{Δ 0 i} ( k ) = | _{x 0} ( k ) - _{x i} ( k ) | is the difference of absolute value between _{x 0} ( k ) and _{x i} ( k ) ; ξ is the distinguishing coefficient ( 0,1 ) ; _{Δ min ...} is the smallest value of _{Δ 0 i} ( k ) ; and _{Δ max ...} is the largest value of _{Δ 0 i} ( k ) . After averaging the grey relational coefficients, the grey relational grade _{Y i} can be obtained as [figure omitted; refer to PDF] where m is the number of process responses and _{β k} is the weight value.

3. Results and Discussion

3.1. Correlation of Chemical Composition with Backward Trajectory Calculation

A total of 18 ADS events were observed at the Pescadores Islands from November 2002 to April 2006. Results obtained from the backward trajectory calculation indicated that three main types of routes dominated the Asian dusts transported from the Inner Mongolia to the Pescadores Islands, which included the eastward transportation and retraced route (ETR), southeastern transportation and circumrotated route (STC), and straight southeastern transportation route (SST), as illustrated in Figure 3.

Figure 3: Backward trajectories for three types of ADS transportation routes from the Inner Mongolia to the Pescadores Islands.

[figure omitted; refer to PDF]

Table 1 summarizes that the ERT-type ADS events accounted for approximately 33.3% of the total of 18 ADS events, which were mostly accompanied with strong northeastern monsoons occurring from late February to April. When the ADS transported toward the east, the transportation routes would detour toward the southwest due to atmospheric circulation, which usually passed through the Pescadores Islands. The STC-type ADS events were mainly affected by the strength of the northeastern monsoons mostly observed from November to next April, accounting for about 22.2% of the total of 18 ADS events. Weakening the northeastern monsoons would result in a smooth curvature of the backward trajectories. The SST-type ADS events occurred mostly from December to next April, in which air parcels transported almost directly from northern China to the Pescadores Islands, Taiwan.

Table 1: The transportation types and soil sources of 18 ADS events in this study.

Transportation route types	Events	Routes probability	Predilection periods	Asian dust source region
ETR	n = 18	33.3%	Late Feb.~April	Central (83%)
STC	22.2%	Nov.~next April	West (75%)
SST	44.5%	Dec.~next April	East (88%)

(1) ETR is the eastward transportation and retraced route.

(2) STC is the southeastern transportation and circumrotated route.

(3) SST is the straight southeastern transportation route.

According to the probability of transportation routes, it showed that the SST-type ADS was dominant in all ADS transportation routes, accounting for approximately 44.5% of the 18 ADS events. The ERT- and STC-type ADS were ordered the second and third routes, respectively. Moreover, this study further aggregated the routes of 18 ADS events and tried to trace them back to their major source regions. All of the ADS routes were transported along the coastline or passed through the major mega cities of North and/or East China, which could alter the chemical composition of atmospheric dusts while transporting from the Inner Mongolia to the Pescadores Islands. During the transportation processes, Asian dusts would mix with particulate matter emitted from local anthropogenic, industrial, and mobile sources. This study further compared the chemical composition of resuspended source soils and atmospheric dusts for similar transportation route type, as illustrated in Figure 4.

Figure 4: Chemical fingerprints of resuspended source soils and atmospheric dusts during the ADS events.

[figure omitted; refer to PDF]

The results indicated that the most abundant chemical composition of source soils were EC, ^{_{SO 4} 2 -} , Cl^- , ^{_{NO 3} -} , OC, and Ca²⁺ in the eastern source region of Inner Mongolia; EC, Fe, Al, ^{_{SO 4} 2 -} , Ca, and Ca²⁺ were in the central source region; EC, ^{_{SO 4} 2 -} , OC, ^{_{NO 3} -} , Cl^- , and Ca were in the western source region, respectively. On the other hand, the most abundant chemical composition for different transportation routes were ^{_{SO 4} 2 -} , Ca²⁺ , ^{_{NO 3} -} , Cl^- , Na⁺ , and EC for the ERT-type ADS; ^{_{SO 4} 2 -} , Cl^- , Na⁺ , ^{_{NO 3} -} , Al, and Fe for the STC-type ADS; and ^{_{SO 4} 2 -} , EC, ^{_{NO 3} -} , Na⁺ , Cl^- , and Ca²⁺ for the SST-type ADS, respectively. The chemical composition between the soils collected in the eastern and western soil sources of Inner Mongolia and Asian dusts sampled at the Pescadores Islands were similar, except the water-soluble Ca²⁺ and metal Ca. The chemical composition of soils collected at each soil source region was similar but not identical, which could cause the chemical composition variation of atmospheric dusts. Thus, the chemical composition from each soil source was similar to the chemical composition of atmospheric dust on their corresponding transportation route. Among them, the air masses traveled through several major industrial cities of China and transported along the coastal regions by the backward trajectories. Previous studies [29, 30, 32, 33] showed that the chemical proportion on the atmospheric dust was changed by the anthropogenic pollutions or nature sources emission, such as ^{_{SO 4} 2 -} and ^{_{NO 3} -} originated from direct emissions of major industries (such as petrochemical, power plant, and incinerator) or mobile sources (including onshore and offshore). Calcium and aluminum were dominated by stone processing, construction, and the cement industry. EC was originated primarily from direct emissions of fuel and coal combustion.

This study further compared the chemical composition of the source soils and atmospheric dusts sampled for different transportation routes. The results showed that the chemical fingerprints of the most soil sources and atmospheric dusts were coincided, except some chemical species were different. Actually, the chemical composition of atmospheric dusts could be influenced by the anthropogenic emission sources for different ADS transportation routes, thus resulting in the increase of the portion of ^{_{SO 4} 2 -} , ^{_{NO 3} -} , and EC. Moreover, Asian dusts could be long-range transported across the ocean to the Pescadores Islands, which could increase the proportion of Na⁺ and Cl^- due to oceanic spray. As a whole, the source soils were not exactly the same in their chemical composition between the eastern, central, or western Inner Mongolia, but the chemical species of the source soils in the Inner Mongolia and atmospheric dusts in the Pescadores Islands were similar in the same transportation routes. It could be explained that both backward trajectory calculation and chemical composition can be used to decide the ADS transportation routes concurrently.

3.2. Correlation of Chemical Composition with Enrichment Factors and Grey Relational Analysis

The enrichment factor was further applied to identify the ADS source region in the Inner Mongolia toward the Pescadores Islands in this study. During the transportation processes, water-soluble ions and carbonaceous components of Asian dusts could be interfered by the anthropogenic emissions. However, the characteristics of metallic elements are less likely varied in the ambient air. Consequently, this study selected Al or Fe as trace metals, which are fairly representative in the earth crust, and can be used to calculate the EF of 15 metals of Asian dusts originated from the Inner Mongolia during the 18 ADS events, as shown in Table 2. The EF values of 15 metal compositions in the range of 0.1 to 10 were shown to be highly relevant, as illustrated in Figure 5. Source identification results showed that 72% of ADS out of 18 events were originated from the central soil source region of the Inner Mongolia. However, three ADS events (numbers 3, 7, and 8) were not consistently identified for their source regions by using Al and Fe as reference elements. The source regions identified by EF values were either the western or the central Inner Mongolia. It might be attributed that Asian dusts could be mixed with other anthropogenic particles during the ADS transportation processes.

Table 2: Source allocation of 18 ADS events in the Inner Mongolia by enrichment factor analysis.

ADS events	Dates of ADS events	Reference elements	ADS source region	Identified sources
East	Central	West
1	Mar. 8-9, 2002	Al	14	9	12	East
Fe	13	9	12
2	Mar. 18-20, 2002	Al	14	11	13	East
Fe	14	13	13
3*	Mar. 31, 2002	Al	11	12	13	West, central
Fe	11	12	11
4	Apr. 13, 2002	Al	7	9	8	Central
Fe	7	9	7
5	Apr. 17-19, 2002	Al	11	13	12	East
Fe	11	8	9
6	Feb. 25, 2003	Al	12	15	12	Central
Fe	13	14	12
7*	Mar. 7-8, 2003	Al	12	15	15	West, central
Fe	13	13	13
8*	Mar. 26-27, 2003	Al	13	15	15	West, central
Fe	13	13	13
9	Apr. 26-27, 2003	Al	12	14	14	Central
Fe	12	13	12
10	Feb. 6, 2004	Al	10	12	10	Central
Fe	10	12	9
11	Feb. 15-16, 2004	Al	10	13	10	Central
Fe	11	11	11
12	Apr. 2, 2004	Al	11	13	10	Central
Fe	10	13	11
13	Apr. 13-15, 2004	Al	11	14	14	Central
Fe	12	14	12
14	Nov. 29, 2005	Al	12	15	14	Central
Fe	12	15	13
15	Dec. 21-26, 2005	Al	13	14	15	West
Fe	13	11	13
16	Mar. 14-17, 2006	Al	9	13	9	Central
Fe	10	12	12
17	Mar. 19, 2006	Al	8	13	8	Central
Fe	8	13	12
18	Apr. 24-25, 2006	Al	10	15	10	East
Fe	12	15	12

^* Different source regions determined by EF values using Al and Fe as reference elements.

Figure 5: EF values of 15 metals with Al or Fe as reference elements at the Pescadores Islands. For example, (a1)~(a3) are the EF of eastern, central, and western source regions calculated by using aluminum as reference element; (b1)~(b3) are the EF of eastern, central, and western source regions calculated by using iron as reference element.

[figure omitted; refer to PDF]

This study further applied the grey relational analysis to allocate the source regions of 18 ADS events, as summarized in Table 3. The results indicated that approximately 27, 23, and 50% of Asian dusts were originated from the eastern, western, and central Inner Mongolia for the 18 ADS events, respectively. Moreover, if the results of grey relational analysis were close to 1, the soil source of ADS events could be identified with relatively high correlation. Unlike the results of EF analysis, the source regions would be determined decisively by the grey relational analysis. However, if the grey relational grades of three soil source regions were close to each other, it could be difficult to allocate their sources in this specific case. For example, number 13 ADS event was identified from the central region of the Inner Mongolia by backward trajectory calculation. Its GR grades were shown as 0.78, 0.91, and 0.89 for the eastern, central, and western soil source regions, respectively. Similar results were obtained from EF analysis, indicating that the chemical composition of Asian dusts played a crucial role in identifying the soil source regions of Asian dusts. Overall, the chemical-based enrichment factor and the grey relational analysis both showed that the soil source region of number 13 ADS event was highly likely emitted from the central Inner Mongolia.

Table 3: Source allocation of 18 ADS events in the Inner Mongolia by grey relational calculation.

ADS events	Dates of ADS events	ADS source regions	Identified sources
East	Central	West
1	Mar. 8-9, 2002	0.90	0.87	0.82	East
2	Mar. 18-20, 2002	0.81	0.76	0.74	East
3	Mar. 31, 2002	0.79	0.84	0.87	West
4	Apr. 13, 2002	0.81	0.74	0.79	East
5	Apr. 17-19, 2002	0.79	0.78	0.71	East
6	Feb. 25, 2003	0.85	0.92	0.84	Central
7	Mar. 7-8, 2003	0.87	0.91	0.89	Central
8	Mar. 26-27, 2003	0.78	0.91	0.83	Central
9	Apr. 26-27, 2003	0.65	0.87	0.83	Central
10	Feb. 6, 2004	0.77	0.83	0.85	West
11	Feb. 15-16, 2004	0.70	0.84	0.78	Central
12	Apr. 2, 2004	0.83	0.86	0.93	West
13	Apr. 13-15, 2004	0.78	0.91	0.89	Central
14	Nov. 29, 2005	0.76	0.87	0.83	Central
15	Dec. 21-26, 2005	0.72	0.78	0.89	West
16	Mar. 14-17, 2006	0.63	0.87	0.84	Central
17	Mar. 19, 2006	0.69	0.79	0.80	Central
18	Apr. 24-25, 2006	0.92	0.89	0.83	East

3.3. Comparison of Three Source Identification Methods

Figure 6 illustrates and compares three methods for identifying the source regions of 18 ADS events. Through the backward trajectory calculation, enrichment factor and grey relational analysis could be applied to identify the ADS events numbers 1, 2, 6, 8, 9, 11, and 13-18, which accounted for 66.7% of the 18 ADS events. Comparing backward trajectory calculation with enrichment factor, the similarity percentage was 77.8% for 18 ADS events, including ADS events numbers 1, 2, 5-9, 11, 12, and 14-18 in this study. In terms of enrichment factor and grey relational analysis, it was shown that they were coincident for 83.3% of 18 ADS events, including ADS events numbers 1, 2, 4-6, 8-11, 12, and 14-18. In terms of grey relational analysis and backward trajectory calculation, it was shown that they were coincident for 77.8% of 18 ADS events, including ADS events numbers 1, 2, 4-6, 8-11, and 14-18. Generally speaking, the source allocation methods of ADS transported from Inner Mongolia to the Pescadores Islands were relatively accurate, especially for simultaneous application of enrichment factor and grey relational analysis. More importantly, accurate chemical fingerprints of source soils and atmospheric particles are highly required for chemically allocating the source regions of the ADS events.

Figure 6: Comparison of the correlation percentages among backward trajectory calculation, enrichment factor, and grey relational analysis.

[figure omitted; refer to PDF]

In order to further confirm the above three optimal methods feasible for allocating the mostly possible sources of ADS events, this study used the past general method to compare with the above three methods, which calculated the mass ratios of Al and Sr (Al/Sr), which had the highest differences in the 15 metals for 18 ADS events. Table 4 lists the mass ratio of Al/Sr for resuspended soils collected in the Inner Mongolia and atmospheric particles sampled at the Pescadores Islands. The results showed that the resuspended soil samples in the central Inner Mongolia had relatively higher Al/Sr ratio ( 325.0 ± 118.0 ) and observed simultaneously with the higher Al/Sr ratio ( 161.1 ± 65.1 ) of atmospheric particles at the Pescadores Islands. Furthermore, the resuspended soils and Asian dusts in the eastern Inner Mongolia both had high relative standard deviation (RSD) of 77 and 283%, respectively. The Al/Sr ratios of the resuspended soils and Asian dusts in the eastern and central Inner Mongolia had their significant characteristics, which could be used to identify the soil source regions of ADS events, while the Al/Sr ratios in the western Inner Mongolia were relatively broad. Overall, although the Al/Sr ratio could assist in identifying the sources of ADS, its confidence level was relatively lower than the backward trajectory calculation, enrichment factor, and grey relational analysis. Finally, this study has selected the mostly possible analysis methods to identify the sources of ADS events, which could be used to trace the source areas of ADS event in the future.

Table 4: Comparison of the Al/Sr ratios for source soils and atmospheric particles sampled in the Inner Mongolia and at the Pescadores Islands.

Source regions	Soil samplesAl/Sr	Dust samplesAl/Sr	Number of ADS events
East	55.6 ± 4.3	33.0 ± 93.3	1, 2, 4, 5, 18
Central	325.0 ± 118.0	161.1 ± 65.1	6-9, 11-14, 16, 17
West	2.5 ± 0.8	132.9 ± 42.6	7, 10, 12, 15

4. Conclusions

This study successfully combined backward trajectory calculation, enrichment factor, and grey relational analysis to allocate the potential source(s) of 18 ADS events, and identified the transportation routes of Asian dusts transported from the Inner Mongolia to the Pescadores Islands. Three major transportation routes observed for 18 ADS events included ETR, STC, and SST. The variations of chemical components of the atmospheric dusts or source soils were correlated to their transportation routes, in which the chemical characteristics of Asian dusts could be influenced by anthropogenic particles emitted from local stationary and mobile sources during the ADS transportation processes. However, most of chemical components between source soils and atmospheric dusts for the same transportation routes were coincident. Moreover, source identification results showed that 50-72% out of 18 ADS events were originated from the central region of Inner Mongolia. Comparison of the three methods showed that the chemical-assisted methods demonstrated the accuracy of 66.7% for identifying the source regions of ADS events. Particularly, simultaneous application of enrichment factor and grey relational analysis raised the accuracy up to 83.3%. Overall, these three chemical-assisted tracing methods were feasible for identifying the soil source regions of ADC events.

Acknowledgments

The authors gratefully acknowledge the financial support from the National Science Council and kind assistance from the Institute of Environmental Engineering at National Sun Yat-Sen University. The authors would like to express their sincere appreciation for its financial and technological support to accomplish this study.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

References

[1] H. Bian, X. Tie, J. Cao, Z. Ying, S. Han, Y. Xue, "Analysis of a severe dust storm event over China: application of the WRF-dust model," Aerosol and Air Quality Research , vol. 11, no. 4, pp. 419-428, 2011.

[2] J. J. Cao, J. C. Chow, J. G. Watson, F. Wu, Y. M. Han, Z. D. Jin, Z. X. Shen, Z. S. An, "Size-differentiated source profiles for fugitive dust in the Chinese Loess Plateau," Atmospheric Environment , vol. 42, no. 10, pp. 2261-2275, 2008.

[3] M. T. Cheng, W. C. Chou, C. P. Chio, S. C. Hsu, Y. R. Su, P. H. Kuo, B. J. Tsuang, S. H. Lin, C. C. K. Chou, "Compositions and source apportionments of atmospheric aerosol during Asian dust storm and local pollution in central Taiwan," Journal of Atmospheric Chemistry , vol. 61, no. 2, pp. 155-173, 2008.

[4] R. A. Duce, C. K. Unni, B. J. Ray, J. M. Prospero, J. T. Merrill, "Long-range atmospheric transport of soil dust from Asia to the tropical North Pacific: temporal variability," Science , vol. 209, no. 4464, pp. 1522-1524, 1980.

[5] X. Y. Zhang, S. L. Gong, Z. X. Shen, F. M. Mei, X. X. Xi, L. C. Liu, Z. J. Zhou, D. Wang, Y. Q. Wang, Y. Cheng, "Characterization of soil dust aerosol in China and its transport and distribution during 2001 ACE-Asia: 1. Network observations," Journal of Geophysical Research D: Atmospheres , vol. 108, no. 9, pp. 4261-4274, 2003.

[6] X. Y. Zhang, R. Arimoto, Z. S. An, "Dust emission from Chinese desert sources linked to variations in atmospheric circulation," Journal of Geophysical Research D , vol. 102, no. 23, pp. 28041-28047, 1997.

[7] C. S. Yuan, C. X. Hai, M. Zhao, "Source profiles and fingerprints of fine and coarse sands resuspended from the soils sampled in the central inner mongolia," Particuology , vol. 4, pp. 304-311, 2006.

[8] Y. I. Tsai, C. L. Chen, "Characterization of Asian dust storm and non-Asian dust storm PM_2.5 aerosol in southern Taiwan," Atmospheric Environment , vol. 40, no. 25, pp. 4734-4750, 2006.

[9] I. Mori, M. Nishikawa, T. Tanimura, H. Quan, "Change in size distribution and chemical composition of kosa (Asian dust) aerosol during long-range transport," Atmospheric Environment , vol. 37, no. 30, pp. 4253-4263, 2003.

[10] S. Tan, G. Shi, H. Wang, "Long-range transport of spring dust storms in Inner Mongolia and impact on the China seas," Atmospheric Environment , vol. 46, pp. 299-308, 2012.

[11] F. Tsai, G. T. J. Chen, T. H. Liu, W. D. Lin, J. Y. Tu, "Characterizing the transport pathways of Asian dust," Journal of Geophysical Research , vol. 113, 2008.

[12] R. D. Brook, "Is air pollution a cause of cardiovascular disease? Updated review and controversies," Reviews on Environmental Health , vol. 22, no. 2, pp. 115-137, 2007.

[13] C. Kang, W. Kim, H. Ko, S. Hong, "Asian Dust effects on Total Suspended Particulate (TSP) compositions at Gosan in Jeju Island, Korea," Atmospheric Research , vol. 94, no. 2, pp. 345-355, 2009.

[14] S. Wang, S. Tsay, N. Lin, N. C. Hsu, S. W. Bell, C. Li, Q. Ji, M. Jeong, R. A. Hansell, E. J. Welton, B. N. Holben, G. Sheu, Y. Chu, S. Chang, J. Liu, W. Chiang, "First detailed observations of long-range transported dust over the northern South China Sea," Atmospheric Environment , vol. 45, no. 27, pp. 4804-4808, 2011.

[15] M. V. K. Sivakumar, "Impacts of sand storms/dust storms on agriculture," Natural Disasters and Extreme Events in Agriculture , pp. 159-177, 2005.

[16] D. Y. Jeong, "Socio-economic costs from yellow dust damages in South Korea," Transborder Environmental and Natural Resource Management , pp. 185-198, 2008.

[17] H. Mu, S. Otani, M. Shinoda, Y. Yokoyama, K. Onishi, T. Hosoda, M. Okamoto, Y. Kurozawa, "Long-term effects of livestock loss caused by dust storm on mongolian inhabitants: a survey 1 year after the dust storm," Yonago Acta Medica , vol. 56, no. 1, pp. 39-42, 2013.

[18] H. Lee, Y. Honda, Y. H. Lim, Y. L. Guo, M. Hashizume, H. Kim, "Effect of Asian dust storms on mortality in three Asian cities," Atmospheric Environment , vol. 89, pp. 309-317, 2014.

[19] C. S. Yuan, C. C. Sau, M. C. Chen, M. H. Huang, S. W. Chang, Y. C. Lin, C. G. Lee, "Mass concentration and size-resolved chemical composition of atmospheric aerosols sampled at the Pescadores Islands during Asian dust storm periods in the years of 2001 and 2002," Terrestrial, Atmospheric and Oceanic Sciences , vol. 15, no. 5, pp. 857-879, 2004.

[20] B. Lee, H. K. Lee, N. Jun, "Analysis of regional and temporal characteristics of PM₁₀ during an Asian dust episode in Korea," Chemosphere , vol. 63, no. 7, pp. 1106-1115, 2006.

[21] S. Chen, L. Hsieh, M. Kao, W. Lin, K. Huang, C. Lin, "Characteristics of particles sampled in southern Taiwan during the Asian dust storm periods in 2000 and 2001," Atmospheric Environment , vol. 38, no. 35, pp. 5925-5934, 2004.

[22] C. J. Ma, M. Kasahara, S. Tohno, R. Holler, T. Kamiya, "Chemical composition of single raindrops fallen in yellow rainfall episode in Japan," in Proceedings of the 2nd Asian Aerosol Conference, pp. 185-186, Pusan, Republic of Korea, 2001.

[23] S. H. Liu, C. S. Yuan (Chinese) Source apportionment and optical properties of atmospheric aerosols in metro Kaohsiung [M.S. thesis] , Institute of Environmental Engineering, National Sun Yat-sen University, Kaohsiung, Taiwan, 2000.

[24] G. McTainsh, Y. Chan, H. McGowan, J. Leys, K. Tews, "The 23rd October 2002 dust storm in eastern Australia: characteristics and meteorological conditions," Atmospheric Environment , vol. 39, no. 7, pp. 1227-1236, 2005.

[25] Z. X. Shen, J. J. Cao, R. Arimoto, Z. W. Han, R. J. Zhang, Y. M. Han, S. X. Liu, T. Okuda, S. Nakao, S. Tanaka, "Ionic composition of TSP and PM 2.5 during dust storms and air pollution episodes at Xi'an, China," Atmospheric Environment , vol. 43, pp. 2911-2918, 2009.

[26] T. Nakamura, K. Matsumoto, M. Uematsu, "Chemical characteristics of aerosols transported from Asia to the East China Sea: An evaluation of anthropogenic combined nitrogen deposition in autumn," Atmospheric Environment , vol. 39, no. 9, pp. 1749-1758, 2005.

[27] R. Zhang, Z. Shen, T. Cheng, M. Zhang, Y. Liu, "The elemental composition of atmospheric particles at Beijing during Asian dust events in spring 2004," Aerosol and Air Quality Research , vol. 10, no. 1, pp. 67-75, 2010.

[28] S. F. Kong, B. Han, Z. P. Bai, L. Chen, J. W. Shi, Z. Xu, "Receptor modeling of PM_2.5 , PM₁₀ and TSP in different seasons and long-range transport analysis at a coastal site of Tianjin, China," Science of the Total Environment , vol. 408, no. 20, pp. 4681-4694, 2010.

[29] L. Zhu, X. Huang, H. Shi, X. Cai, Y. Song, "Transport pathways and potential sources of PM₁₀ in Beijing," Atmospheric Environment , vol. 45, no. 3, pp. 594-604, 2011.

[30] H. H. Tsai, C. S. Yuan, C. H. Hung, Y. C. Lin, "Physicochemical properties of PM_2.5 and PM_2.5-10 at inland and offshore sites over Southeastern coastal region of Taiwan strait," Aerosol and Air Quality Research , vol. 11, no. 6, pp. 664-678, 2011.

[31] J. L. Deng A Course on Grey System Theory , HUST Press, Wuhan, China, 1990.

[32] G. Dongarrà, E. Manno, D. Varrica, M. Lombardo, M. Vultaggio, "Study on ambient concentrations of PM₁₀ , PM_10-2.5 , PM_2.5 and gaseous pollutants: trace elements and chemical speciation of atmospheric particulates," Atmospheric Environment , vol. 44, no. 39, pp. 5244-5257, 2010.

[33] T. C. Li, W. H. Chen, C. S. Yuan, S. P. Wu, X. H. Wang, "Physicochemical characteristics and source apportionment of atmospheric aerosol particles in Kinmen-Xiamen Airshed," Aerosol and Air Quality Research , vol. 13, no. 1, pp. 308-323, 2013.