Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Asian Dust Storm Elevates Children’s Respiratory Health Risks: A Spatiotemporal Analysis of Children’s Clinic Visits across Taipei (Taiwan)

  • Hwa-Lung Yu,

    Affiliation Department of Bioenvironmental Systems Engineering, National Taiwan University, Taipei, Taiwan

  • Lung-Chang Chien,

    Affiliation Department of Internal Medicine, Division of Health Behavior Research, Washington University School of Medicine, St. Louis, Missouri, United States of America

  • Chiang-Hsing Yang

    yangch@ntunhs.edu.tw

    Affiliation Department of Health Care Management, National Taipei University of Nursing and Health Sciences, Taipei, Taiwan

Abstract

Concerns have been raised about the adverse impact of Asian dust storms (ADS) on human health; however, few studies have examined the effect of these events on children’s health. Using databases from the Taiwan National Health Insurance and Taiwan Environmental Protection Agency, this study investigates the documented daily visits of children to respiratory clinics during and after ADS that occurred from 1997 to 2007 among 12 districts across Taipei City by applying a Bayesian structural additive regressive model controlled for spatial and temporal patterns. This study finds that the significantly impact of elevated children’s respiratory clinic visits happened after ADS. Five of the seven lagged days had increasing percentages of relative rate, which was consecutively elevated from a 2-day to a 5-day lag by 0.63%∼2.19% for preschool children (i.e., 0∼6 years of age) and 0.72%∼3.17% for school children (i.e., 7∼14 years of age). The spatial pattern of clinic visits indicated that geographical heterogeneity was possibly associated with the clinic’s location and accessibility. Moreover, day-of-week effects were elevated on Monday, Friday, and Saturday. We concluded that ADS may significantly increase the risks of respiratory diseases consecutively in the week after exposure, especially in school children.

Introduction

Every winter and spring season, Asian dust storms (ADS) are frequently a major health problem in many cities throughout the world. Annual dust emissions total approximately 43 million tons, with springtime dust emissions accounting for almost half of the annual totals [1]. ADS usually originate in the deserts of Mongolia and China, lift soil particles into the atmosphere, and carry them east and southeast towards Korea, Japan, Taiwan, and even the Philippines [2], [3]. In addition to highly-elevated ambient particulate matter (PM) concentration, the PM composition during ADS differs from that during non-ADS periods. Reports from central Taiwan show that the concentration of the crustal elements and sea salt species in PM2.5-10 (i.e., particulate matters with granulometric diameter between 2.5 and 10 µm) during ADS exceeds the mean concentration of these same entities during non-ADS periods by at least 2 units of factors [4]. Moreover, during ADS in China, higher PM concentrations were measured in Korea and Japan, and the ADS particles were also observed to be rich in some elements that include aluminum, iron, and calcium [5][7].

Several environmental epidemiologic studies have demonstrated evidence of an adverse relationship between airborne particles and human health. Airborne particles uniquely impact the mortality and hospitalization rates of several diseases, most notably respiratory diseases. Most research has concluded that the effects of PM, especially composed of fine particles, are more severe than that of other air pollutants because they can be inhaled deeply into the lungs [8]. A European study revealed that each 10 µg/m3 increase of PM10 (i.e., particulate matters with granulometric diameter less than 10 µm) raises the relative rate of clinic visit by 1.20% for asthma and 1.00% for chronic obstructive pulmonary disease [9]. In addition, the effect of fine particulate air pollution on mortality may contribute to higher instances of stroke and respiratory deaths [10]. Besides PM, the adverse health impact of CO, O3, NO2 and SO2 has also been well-documented historically [11][13].

The health impact of ADS has been discussed by comparing mortality rates or emergency room admissions between ADS and non-ADS events [14], [15]. In addition, researchers have investigated health related lag effects associated with ADS, and they have discovered deferred influence on hospital admissions [16][18]. Yang et al. [19] found a statistically significant association between ADS events and primary intracerebral hemorrhagic stroke admissions 3 days after the dust storms. ADS also increased the risk of both asthma and cerebrovascular admissions from a 1-day lag to a 3-day lag [20], [21]. Among the research population, children are particularly sensitive to airborne exposure [22]. However, few studies have focused on children’s health from the perspective of ADS events [23], [24]. Most of the previous studies regarding the health impact analysis of ADS utilize relatively few health observations, and therefore, inferences from these studies may be limited and conservative [25]. Chien et al. [24] showed the elevated rate of children’s respiratory clinic usage during one week following ADS; however, no temporal lag effect structure of health impact was discussed.

Geographic heterogeneity has been a salient factor in ambient pollutant distributions [26], [27] as well as its associations with health outcomes [12], [28]. Nevertheless, few studies have assessed the spatial variation of an ADS’s impact on human health. In order to address this issue, this study applies a unidirectional approach [29] under a spatiotemporal model framework to diagnose the space-time disparity of children’s respiratory clinic visits. The influence of ADS on children’s health was examined by considering temporal lag effects starting from the end of each ADS event as well as the spatial variation over study areas. This study specifically investigates the daily clinic visits of children with respiratory diseases in 12 districts in Taipei City from 1997 to 2007.

Materials and Methods

Children’s Clinic Data

Initiated in March 1995, Taiwan’s National Health Insurance (NHI) program contacted more than 97% of hospitals and clinics nationwide within its first year of inception, enrolling more than 96% of Taiwanese residents. The Taiwan National Health Research Institute maintains the NHI program database, and has established a standard procedure that assures the quality and accuracy of claims data [30]. The NHI database includes ambulatory care expenditures by visit as well as the registries of contracted medical facilities nationwide. The procedure and diagnostic codes are used to retrieve cause-specific data according to diagnosis-related groups or International Classification of Diseases, Ninth Revision, Clinical Modification (ICD-9-CM) classification codes by the Bureau of National Health Insurance. Regarding personal privacy and confidentiality, all individually identifiable health information has been encrypted prior to release, e.g., personal identification or hospital identification numbers.

For this study, the following respiratory diseases were highlighted: acute respiratory infections (ICD-9∶460–466, allergic rhinitis (ICD-9∶477), other diseases of upper respiratory tract (ICD-9∶478), pneumonia and influenza (ICD-9∶480–488), asthma (ICD-9∶493), bronchiectasis (ICD-9∶494), and extrinsic allergic alveolitis (ICD-9∶495). This study obtained a population-based database containing space-time data for clinic and hospital visits (i.e., hospital location and appointment times) for all-cause respiratory diseases of children under 14 years old in Taipei City from 1997–2007, including both ambulatory and emergency visits. We split these clinic visits data into preschool children (0∼6 years of age) and school children (7∼14 years of age) in this analysis.

Dust Storm Data

ADS often occur in northern and northwestern China, impacting Taiwan only under certain atmospheric circumstances. Before 2000, the Department of Atmospheric Science at Chinese Culture University (CCU) was tasked with the responsibility of characterizing, defining, and monitoring ADS events in Taiwan. They came up with the following criteria for identifying an ADS event: 1) dust storm events with PM10 concentrations >100 µg/m3 observed by any air quality monitoring stations located in Wanli, Guanyin, Danshui, and Yilan, and 2) dust storm events with visibility less than 1 km for 24 hours in any of three neighboring First Global GARP Experiment-type ground stations [31]. After 2000, the Taiwan Environmental Protection Agency (TWEPA) became the official organization to define, monitor, and predict ADS. They utilize three distinct steps in order to categorize storm events in Taiwan as ADS and to predict their arrival time. First, the Weather Integration and Nowcasting System is consulted to confirm the occurrence of ADS in Mongolia and China. Second, the Moderate Resolution Imaging Spectroradiometer remote data and several models of ADS are used to track the transport of ADS and to predict the probability of the arrival of ADS in Taiwan. Third, if the data confirm that ADS may blow to Taiwan, the TWEPA will issue an early warning with the estimated arrival date [32]. This study highlights and analyzes 76 storm events that were thus categorized as ADS during the period of 1997–2007 (see Table 1). These ADS events span a total period of 172 dust storm days. Proportionally, these dust storm days account for 4.28% of the entire study period that lasted for a total time span of 4017 days. In addition to the ADS data, ambient pollutants concentrations and temperature have been regularly monitored at the TWEPA stations across Taiwan since 1994. The temperature measurements used in this analysis are the daily observations at the Jhongshan air quality monitoring station located in the most populated area of Taipei City (see Figure 1).

thumbnail
Figure 1. Spatial distribution of the districts, the medical centers, the MRT lines, and the Jhongshan air quality monitoring station in Taipei City.

https://doi.org/10.1371/journal.pone.0041317.g001

Study Area

Taipei City, located in northern Taiwan, is the country’s capital and largest metropolitan area with a population of more than six million inhabitants. Topographically, it is the second largest basin in Taiwan, bounded by the Yangming Mountains to the north, the Linkou mesa to the west, and Snow Mountains to the southeast. The basin region along the rivers is more populated than regions near the mountains, creating serious air pollution problems due to its heavy traffic. In the past decade, Taipei’s metropolitan rapid transit system (MRT) has effectively reduced cross-district traveling times within Taipei, and therefore, it has increased the accessibility to major medical facilities. Figure 1 displays the Taipei City map and essential urban information relevant to this study, including geographic topography, districts, major medical facilities, and the MRT. The clinic visit data are aggregated with respect to the 12 districts across Taipei City.

Spatiotemporal Modeling

In this study, assuming is the number of daily children’s respiratory clinic visits at calendar time in district , this outcome variable would follow a Poisson distribution by with the expected value and variance , where is an over-dispersion parameter representing the variation of clinic visits unable to be calculated by statistical models. A vector of dust storm lag index (DSLI) contains eight dummy variables to represent 0-lag to 7-lag days of ADS. Moreover, a vector of day-of-week (DOW) dummy variables from Monday to Saturday is used to control short-term temporal autoregressive correlations. The B-spline with a second order random walk prior is used in a time smoother for calendar time to consider long-term autoregressive correlations and in a temperature smoother to control nonlinear weather confounding effect. Therefore, a model framework can be constructed by applying a Bayesian structural additive regression (STAR) modeling approach:

(1)

where is an intercept for interpreting the overall association for all districts. The parameters and are a 1×6 vector with six coefficients of day-of-week variables and a 1×8 vector with eight coefficients of DSLI variables. The offset is the logarithm of the district-level population based on the 2000 Census.

In particular, a spatial function was appended to consider potential spatial autocorrelations among the 12 districts. It is actually a Markov random field [33] achieved by a conditional autoregressive prior to following a normal distribution , where is the number of adjacent districts s’ connected to district s, and represents that the district s’ belongs to the set of the neighboring districts of district s. Note that the geographic information system in this study is boundary data, and that the definition of two connected districts means that they share parts of boundaries. In the spatial function, the unknown variance, , and smoothing parameter are assumed to follow an inverse Gamma distribution IG(0.001, 0.001). The spatial effect can be explained by the relative rate (RR) in district s compared with the mean value for the whole population controlling for spatial autocorrelations [34]. Maps of the spatial function at the district level were presented to visualize the geographic distribution of RR. Spatial effects were also classified into three groups according to their posterior probabilities with respect to the number 1 in the following manner: (i) 80% of the posterior distribution below 1 representing a significantly lower RR than the mean level for the Taipei City area; (ii) 80% of the posterior distribution above 1 representing a significantly higher RR than the mean level for the Taipei City area; (iii) the other districts representing a non-significant difference of RR compared to the mean level for the Taipei City area.

All unknown parameters in this STAR model were estimated by an empirical Bayesian approach using a restricted maximum likelihood (REML) estimation technique [35]. The REML method is a substitution of the maximum likelihood method that avoids the loss of degrees of freedom while estimating fixed effects, and prevents estimators toward zero with biases. Demographic analyses were generated using SAS v9.2 software [36], and the model-based estimation procedure was accomplished in the software package BayesX v2.01 [37]. A p-value <0.05 was considered statistically significant in estimated coefficients.

Results

Table 2 depicts the regional distribution of children’s average daily respiratory clinic visits in Taipei City from 1997–2007. During ADS events (lag 0), the daily visit average reached its highest value with 1484.30 clinic visits (SD = 477.40) in the Shihlin District, while the Datong District only had 399.67 (SD = 161.72) clinic visits per day. The daily average of clinic visits increased during the 1-day lag at the end of ADS events, with the largest increment occurring in the Daan District (153.62 daily clinic visits). The Jhongjheng District only featured a 27.58 increase. After a 1-day lag, clinic visits decreased in most areas, however this falling trend was inconsistent. Some areas had a higher average during the 5-day lag or 7-day lag. Moreover, according to the daily records at the Jhongshan air quality monitoring station, the average concentration of PM10 was 90.64 µg/m3 during dust storm days, which was significantly higher than the average concentration during non-dust storm days (p-value <0.0001), as shown in Table 3. The average concentration of O3 during dust storm days was 5.61 ppb higher than that during non-dust storm days significantly (p-value <0.0001). This finding confirms that during ADS events as verified by the CCU and TWEPA, higher concentrations of two important air pollutants (i.e., PM10 and O3) were measured.

thumbnail
Table 2. Daily average children’s respiratory clinic visits within 12 Taipei City districts, 1997–2007.

https://doi.org/10.1371/journal.pone.0041317.t002

thumbnail
Table 3. Mean levels of air pollutants on dust storm days and non-dust storm days in Taipei City, 1997–2007.

https://doi.org/10.1371/journal.pone.0041317.t003

Table 4 shows the increased percentage of RR of children’s respiratory clinic visits in certain DOW variables. Monday had the significantly highest increased percentage 37.55% (p-value <0.0001; 95% CI = 37.44, 37.65) compared to Sunday in preschool children and 37.88% (p-value <0.0001; 95% CI = 37.73, 38.08) in school children. Saturday was the second leading DOW variable, and the percentage increase of RR in school children (10.76%, 95% CI = 10.63, 10.88) was almost twice the percentage in preschool children (5.66%, 95% CI = 5.56, 5.75). Wednesday had the lowest percentage increase of RR of only 0.56% (p-value <0.0001; 95% CI = 0.47, 0.65) in preschool children, while the percentage largely inflated to 7.64% (p-value <0.0001; 95% CI = 7.51, 7.77) in school children. For all children, the greatest percentage increase of RR was 37.64%, occurred on Monday (p-value <0.0001; 95% CI = 37.55, 37.72).

thumbnail
Table 4. Percentage change in rates of daily children’s respiratory clinic visits in Taipei City, 1997–2007 [% (95% CI)].

https://doi.org/10.1371/journal.pone.0041317.t004

The association between ADS and children’s respiratory clinic visits was not positive until the second day after ADS, suggesting the percentage increase of RR for preschool children was −2.53% (p-value <0.0001; 95% CI = −2.69, −2.36), while it was much lower in school children by −6.28% (p-value <0.0001; 95% CI = −6.50, −6.06). The negative association lasted through the 1-day lag, and became positive from the 2-day lag to the 7-day lag, except for the 6-day lag. Among lag days with positive associations, preschool children had the highest 2.19% (p-value <0.0001; 95% CI = 1.95, 2.43) at the 3-day lag; meanwhile for school children, RR at the 7-day lag reached its highest percentage increase by 3.20% (p-value <0.0001, 95% CI = 2.81, 3.60). Regardless of age stratification, the strongest association happened at the 3-day lag with a 2.40% (p-value <0.0001; 95% CI = 2.20, 2.59) increase in RR for all children.

Figure 1 depicts the distribution of spatial effects attributed to children’s respiratory clinic visits in Taipei City. In most districts with preschool children or school children, positive risk of increased clinic visits was prevalent. The range of spatial effect in preschool children was (−0.37, 0.27), and it was wider than that for school children (−0.22, 0.13). Out of all the 12 districts studied, the Jhongshan District displayed the strongest spatial effect contributed to respiratory clinic visits for both preschool and school children. The positive spatial effect in preschool children’s clinic visits was uniformly distributed in the following districts: the Beitu, Shihlin, Neihu, Nangang, and Wanshan District; however, no specific pattern described the spatial heterogeneity in school children’s clinic visits. The maps of 80% posterior probability show that 6 of 12 districts demonstrated a significantly positive spatial effect in both preschool children and school children. Combing two groups, 7 of 12 districts had a significantly positive spatial effect, which locations are identical to the finding in preschool children.

Discussion

Due to high PM concentrations and unusual PM compositions during ADS, ADS and their occurrences have been considered to pose a high risk to human health. Recent studies have demonstrated potential health risks associated with ADS in terms of higher mortality rates and hospital admissions [21], [38]. Some studies have evaluated the biological plausibility of the ADS to induce adverse health effects. These studies have shown that the particles in ADS can exert toxicological effects on the respiratory system, such as causing pulmonary inflammation and inducing cytotoxicity in rat alveolar cells [39][44]. However, other epidemiological studies have noted that the relationship between ADS and adverse health effects, particularly respiratory diseases, is at best uncertain or statistically insignificant [3], [21], [45][48]. One plausible explanation of these inconsistencies may be that the health assessment measures used to evaluate the health impact of the ADS did not adequately capture the health effects. For instance, some of the health measures utilized only accounted for severe cases that required inpatient care, and consequentially, negative environmental events such as ADS may not necessarily induce such severe medical conditions. Furthermore, these previous analyses were based upon observations from limited hospitals [21], [25], [45], [48], and reflected severe health conditions that were exhibited by the most vulnerable individuals within a general population. In contrast, this study provides complete ambulatory and emergency service utilization information from the NHI, and hopefully, better captures the potential health impact of ADS events on the general population.

In recent decades, the importance of spatiotemporal analysis has been emphasized in environmental epidemiological research, especially in quantifying uncertainties in space-time health and exposure data as well as capturing the resulting impact on the estimates of these associations [28], [49], [50]. Although temporal health impacts caused by ADS have been extensively investigated, previous studies have seldom considered the spatial heterogeneity of such health data in which the geographic disparity may be derived from geographical variations of medical resources and exposure levels.

In response, this study implemented the STAR modeling approach for the spatiotemporal analysis of children’s clinic visits related to ADS. The model identifies temporal patterns of a time process by accounting for linear and nonlinear explanatory variables similar to many time series models, such as the generalized additive model [51]. Moreover, the STAR model reveals the spatial heterogeneity independent of temporal variations by using Markov random fields. In addition, the Bayesian framework of the STAR model allows feasibility to account for the parameter’s uncertainty. This novel approach provides a more comprehensive perspective on the impact of ADS on children’s clinic visits for respiratory illnesses.

Previous studies have demonstrated different results in the temporal lag effects of adverse human health related ADS. For instance, hospital admissions were prominent 2 days after ADS in one asthma study [45]. Also, in another study, a positive influence was noted between ADS and ischemic stroke hospital admissions on the third day following a dust storm event [19]. At the 1-day lag, the relative risk of the association between ADS and cardiovascular disease hospital admissions is also increased [14]. However, these associations were statistically insignificant. In contrast, hospital admissions records in Taipei City noted a significant increase in ischaemic heart disease admissions at the 2-day lag and asthma admissions at the 3-day lag [21]. Meanwhile, total respiratory diseases at the 3-day lag and upper respiratory tract infection in males at the 4-day lag were significant in Minqin City, China [48]. As noted, these findings were mostly based upon the analysis of more severe health measures. This study conducted a population-based study and found that children’s respiratory health can be affected by ADS. This impact significantly occurred during most days within a week after a dust storm event. The elevated rates for children’s respiratory clinic visits after a dust storm began at the 2-day lag and attained its highest impact at the 3-day lag for preschool children and the 7-day lag for school children.

Several of the following considerations may provide plausible explanations for such an increased rate of children’s respiratory clinic visits: First, the impact of the ADS may not necessarily incite immediate respiratory illness or illness severe enough to necessitate the patient to seek medical services. A latency period may exist between the adverse environmental influence and the onset of respiratory symptoms requiring the need for medical services. Second, the adverse weather conditions that exist during ADS, such as strong winds and low visibility [52], often prevent citizens from going out. Third, the increasing popularity of ADS forecasting by the media and governmental agency may increase the population awareness of ADS and their potential health effects. Fourth, since over-the-counter pharmaceuticals are easily accessible and inexpensive in Taiwan, many Taiwanese residents may prefer initiating treatment of their symptoms and their children’s symptoms with these products before seeking medical treatment at a clinic, especially if they perceive that their condition is not serious. However, if unsuccessful, treatment with these over-the-counter medications could also account for the lag time noticed in the children’s respiratory clinic visits following ADS. Table 4 notes that the consecutive elevated risks may only apply to the children because of their vulnerability to ambient pollutants. Further studies should assess the ADS health impact on other age groups. Table 4 also shows that school children were affected by ADS much easier than preschool children. This may be explained by the fact that Taiwanese schools were not suspended during ADS, and school children may have experienced higher exposures and elevated concentrations of heavy metals and ambient PM compositions than their preschool counterparts who more likely stayed at home [4], [53], [54]. This type of exposure has been closely associated with the reduction of pulmonary functions in children [23]. Further studies should investigate the temporal fluctuation of daily relative risks that may result from the cross-infection of respiratory diseases among children or the influence payment policy of the Taiwanese NHI system.

Interestingly, the day-of-week has been considered as a meaningful confounding factor for clinic visits in Taiwan. It is important to know that local ambulatory service is commonly rendered on a “first-come, first-serve” basis, and that physician appointments are not necessary for the regular weekday schedule. Access to any level of healthcare facility or provider is therefore unconstrained during the week. However, weekend medical services must be justified by severe symptoms and result in higher co-payment requirements (out-of-pocket amount) from the NHI. Thus, the service schedule and payment system are important factors affecting the timing of medical-care-seeking behavior. Thus, the temporal pattern of clinic visits is closely associated with this weekend effect.

In Taiwan, the majority of the medical services in hospitals and clinics are closed from Saturday afternoon until Monday morning. Therefore, there is a strong incentive to visit clinics on Fridays and Saturday (before weekend effect), and on Monday (after weekend effect). The day-of-week clinic visit pattern observed in this study is quite consistent with medical care-seeking pattern that has resulted under the current national health care delivery system. Moreover, the highly elevated RR on Monday essentially comprised those patients seeking medical treatment from Saturday to Monday, especially in the case of children who are incapable of accessing clinic care independently and are dependent on a parent’s working schedule which allows for only nighttime availability [55].

The spatial heterogeneity of clinical visits also reflects children’s respiratory clinic visits (see Figure 2). Compared to Figure 1, the districts with elevated rates were closely linked to the areas with multiple urban medical centers, especially those along the MRT lines, implying that a high usage of ambulatory and emergency services for children’s health might logically occur in these districts. Because the NHI program is characterized by its low co-payments and open access to providers without choice restrictions, it encourages those insured under the current Taiwan NHI system to seek care in these medical centers with minimal personal financial impact. Consequently, each person in Taiwan averages 14.2 clinic visits per year. In addition, some people may seek treatment of common diseases in hospitals or even tertiary medical centers, rather than clinics [24]. Further study is required to investigate the relationship between the identified spatial heterogeneity and the locations of medical centers in the study area.

thumbnail
Figure 2. Spatial effects with 80% posterior probability for (a) preschool children, (b) school children, and (c) all children.

Districts shaded by white color showed significantly positive spatial effect, whereas districts shaded by black color depicted significantly negative spatial effect. Grey color represented non-significant spatial effect in the district.

https://doi.org/10.1371/journal.pone.0041317.g002

Conclusion

In summary, the spatiotemporal analysis presented in this study identifies the temporal pattern of the health risks during and after ADS and analyzes them day-by-day by considering the spatial confounding factor. The study results clearly show significant and increased rates for respiratory clinic visits in the studied population of children over time in 5 of 7 days after ADS. The findings of this population-based study can provide governmental agencies with an important reference source in order to plan and implement policies that can help to both protect children from the possible adverse health effects of ADS and to provide care related to such health effects.

Author Contributions

Conceived and designed the experiments: HLY LCC. Analyzed the data: HLY LCC. Contributed reagents/materials/analysis tools: HLY LCC CHY. Wrote the paper: HLY LCC CHY. Provided empirical data: CHY HLY.

References

  1. 1. Xuan J (1999) Dust emission factors for environment of Northern China. Atmos Environ 33: 1767–1776.
  2. 2. Chen SJ, Hsieh LT, Kao MJ, Lin WY, Huang KL, et al. (2004) Characteristics of particles sampled in southern Taiwan during the Asian dust storm periods in 2000 and 2001. Atmos Environ 38: 5925–5934.
  3. 3. Chen YS, Sheen PC, Chen ER, Liu YK, Wu TN, et al. (2004) Effects of Asian dust storm events on daily mortality in Taipei, Taiwan. Environ Res 95: 151–155.
  4. 4. Cheng MT, Lin YC, Chio CP, Wang CF, Kuo CY (2005) Characteristics of aerosols collected in central Taiwan during an Asian dust event in spring 2000. Chemosphere 61: 1439–1450.
  5. 5. Choi JC, Lee M, Chun Y, Kim J, Oh S (2001) Chemical composition and source signature of spring aerosol in Seoul, Korea. J Geophys Res 106: 18067–18074.
  6. 6. Ma C-J, Kasahara M, Höller R, Kamiya T (2001) Characteristics of single particles sampled in Japan during the Asian dust-storm period. Atmos Environ 35: 2707–2714.
  7. 7. Zhou M, Okada K, Qian F, Wu PM, Su L, et al. (1996) Characteristics of dust-storm particles and their long-range transport from China to Japan – case studies in April 1993. Atmos Res 40: 19–31.
  8. 8. Janssens HM, De Jongste JC, Hop WC, Tiddens HA (2003) Extra-fine particles improve lung delivery of inhaled steroids in infants: a study in an upper airway model. Chest 123: 2083–2088.
  9. 9. Atkinson RW, Anderson HR, Sunyer J, Ayres J, Baccini M, et al. (2001) Acute Effects of Particulate Air Pollution on Respiratory Admissions. Results from APHEA 2 Project. Am J Respir Crit Care Med 164: 1860–1866.
  10. 10. Zanobetti A, Schwartz J (2009) The effect of fine and coarse particulate air pollution on mortality: a national analysis. Environ Health Persp 117: 898–903.
  11. 11. Bell ML, McDermott A, Zeger SL, Samet JM, Dominici F (2004) Ozone and short-term mortality in 95 US urban communities, 1987–2000. JAMA-J Am Med Assoc 292: 2372–2378.
  12. 12. Dominici F, Daniels M, Zeger SL, Samet JM (2002) Air Pollution and Mortality: Estimating Regional and National Dose-Response Relationships. J Am Stat Assoc 97: 100–111.
  13. 13. Chan CC, Chuang KJ, Chien LC, Chen WJ, Chang WT (2006) Urban air pollution and emergency admissions for cerebrovascular diseases in Taipei, Taiwan. Eur Heart J 27: 1238–1244.
  14. 14. Chen YS, Yang CY (2005) Effects of Asian dust storm events on daily hospital admissions for cardiovascular disease in Taipei, Taiwan. J Toxicol Env Heal A 68: 1457–1464.
  15. 15. Chan CC, Chuang KJ, Chen WJ, Chang WT, Lee CT, et al. (2008) Increasing cardiopulmonary emergency visits by long-range transported Asian dust storms in Taiwan. Environ Res 106: 393–400.
  16. 16. Schwartz J (2000) The distributed lag between air pollution and daily deaths. Epidemiol 11: 320–326.
  17. 17. Zanobetti A, Wand MP, Schwartz J, Ryan LM (2000) Generalized additive distributed lag models: quantifying mortality displacement. Biostatistics 1: 279–292.
  18. 18. Zanobetti A, Schwartz J, Samoli E, Gryparis A, Touloumi G, et al. (2002) The temporal pattern of mortality responses to air pollution: a multicity assessment of mortality displacement. Epidemiol 13: 87–93.
  19. 19. Yang CY, Chen YS, Chiu HF, Goggins WB (2005) Effects of Asian dust storm events on daily stroke admissions in Taipei, Taiwan. Environ Res 99: 79–84.
  20. 20. Lee EC, Leem J, Hong YC, Kim H, Kim HC (2008) Effects of Asian Dust Storm Events on Daily Admissions for Asthma and Stroke in Seven Metropolitans of Korea. Epidemiol 19: S145.
  21. 21. Bell ML, Levy JK, Lin Z (2008) The effect of sandstorms and air pollution on cause-specific hospital admissions in Taipei, Taiwan. Occup Environ Med 65: 104–111.
  22. 22. Schwartz J (2004) Air pollution and children’s health. Pediatrics 113: 1037–1043.
  23. 23. Hong YC, Pan XC, Kim SY, Park K, Park EJ, et al. (2010) Asian Dust Storm and pulmonary function of school children in Seoul. Sci Total Environ 408: 754–759.
  24. 24. Chien LC, Yang CH, Yu HL (2012) Estimated Effects of Asian Dust Storms on Spatiotemporal Distributions of Clinic Visits for Respiratory Diseases in Taipei Children (Taiwan). Environ Health Persp. In press.
  25. 25. Middleton N, Yiallouros P, Kleanthous S, Kolokotroni O, Schwartz J, et al. (2008) A 10-year time-series analysis of respiratory and cardiovascular morbidity in Nicosia, Cyprus: the effect of short-term changes in air pollution and dust storms. Environ Health 7: 39.
  26. 26. Yu HL, Wang CH (2010) Retrospective prediction of intraurban spatiotemporal distribution of PM2.5 in Taipei. Atmos Environ 44: 3053–3065.
  27. 27. Bogaert P, Christakos G, Jerrett M, Yu HL (2009) Spatiotemporal modelling of ozone distribution in the State of California. Atmos Environ 43: 2471–2480.
  28. 28. Dominici F, McDermott A, Zeger SL, Samet JM (2002) National Maps of the Effects of Particulate Matter on Mortality: Exploring Geographical Variation. Environ Health Persp 111: 39–44.
  29. 29. Basu R, Dominici F, Samet JM (2005) Temperature and mortality among the elderly in the United States: a comparison of epidemiologic methods. Epidemiol 16: 58–66.
  30. 30. Tseng CH (2004) Mortality and causes of death in a national sample of diabetic patients in Taiwan. Diabetes Care 27: 1605–1609.
  31. 31. Nee JB, Chiang CW, Hu HL, Hu SX, Yu JY (2007) Lidar measurements of Asian dust storms and dust cloud interactions. J Geophys Res 112: D15202.
  32. 32. TWEPA (2011) Asian dust storm monitoring network. Available: http://dust.epa.gov.tw/dust/zh-tw/. Accessed: 2012 Jan 9.
  33. 33. Kindermann R, Snell JL (1980) Markov random fields and their applications. American Mathematical Society.
  34. 34. Fromont A, Binquet C, Sauleau EA, Fournel I, Bellisario A, et al. (2010) Geographic variations of multiple sclerosis in France. Brain 133: 1889–1899.
  35. 35. Fahrmeir L, Lang S (2001) Bayesian Inference for Generalized Additive Mixed Models Based on Markov Random Field Priors. J Roy Stat Soc C-App 50: 201–220.
  36. 36. SAS (2009) SAS OnlineDoc. 9.2 ed. Cary, NC: SAS Institute Inc. Available: http://support.sas.com/documentation/onlinedoc/base/index.html. Accessed 2012 Jan 9.
  37. 37. Brezger A, Kneib T, Lang S (2005) BayesX: Analysing Bayesian structured additive regression models. J Stat Softw 14(11): 1–22.
  38. 38. Meng ZQ, Lu B (2007) Dust events as a risk factor for daily hospitalization for respiratory and cardiovascular diseases in Minqin, China. Atmos Environ 41: 7048–7058.
  39. 39. He M, Ichinose T, Yoshida S, Yamamoto S, Inoue K, et al. (2012) Asian sand dust enhances murine lung inflammation caused by Klebsiella pneumoniae. Toxicol Appl Pharm 258: 237–247.
  40. 40. He M, Ichinose T, Yoshida S, Nishikawa M, Mori I, et al. (2010) Airborne Asian sand dust enhances murine lung eosinophilia. Inhal Toxicol 22: 1012–1025.
  41. 41. Ichinose T, Yoshida S, Hiyoshi K, Sadakane K, Takano H, et al. (2008) The effects of microbial materials adhered to Asian sand dust on allergic lung inflammation. Arch Environ Con Tox 55: 348–357.
  42. 42. Ichinose T, Yoshida S, Sadakane K, Takano H, Yanagisawa R, et al. (2008) Effects of Asian sand dust, Arizona sand dust, amorphous silica and aluminum oxide on allergic inflammation in the murine lung. Inhal Toxicol 20: 685–694.
  43. 43. Kishikawa RT, Imai T, Ichinose T, Okuno Y, Arashidani K, et al. (2012) Effects of Asian Dust and Spherical Particles Exposure on Human Health and Allergic Symptom, Fukuoka, Japan. J Allergy Clin Immun 129: Ab59–Ab59.
  44. 44. Yanagisawa R, Takano H, Ichinose T, Mizushima K, Nishikawa M, et al. (2007) Gene expression analysis of murine lungs following pulmonary exposure to Asian sand dust particles. Exp Biol Med 232: 1109–1118.
  45. 45. Yang CY, Tsai SS, Chang CC, Ho SC (2005) Effects of Asian dust storm events on daily admissions for asthma in Taipei, Taiwan. Inhal Toxicol 17: 817–821.
  46. 46. Yang CY, Cheng MH, Chen CC (2009) Effects of Asian dust storm events on hospital admissions for congestive heart failure in Taipei, Taiwan. Journal of toxicology and environmental health Part A 72: 324–328.
  47. 47. Chiu HF, Tiao MM, Ho SC, Kuo HW, Wu TN, et al. (2008) Effects of Asian dust storm events on hospital admissions for chronic obstructive pulmonary disease in Taipei, Taiwan. Inhal Toxicol 20: 777–781.
  48. 48. Meng Z, Lu B (2007) Dust events as a risk factor for daily hospitalization for respiratory and cardiovascular diseases in Minqin, China. Atmos Environ 41: 7048–7058.
  49. 49. Beyea J (1999) Geographic exposure modeling: a valuable extension of geographic information systems for use in environmental epidemiology. Environ Health Persp 107 Suppl 1181–190.
  50. 50. Yu HL, Chen JC, Christakos G, Jerrett M (2009) BME estimation of residential exposure to ambient PM10 and ozone at multiple time scales. Environ Health Persp 117: 537–544.
  51. 51. Hastie T, Tibshirani R (1990) Generalized additive models. London; New York: Chapman and Hall.
  52. 52. Chang SC, Lee CT (2007) Assessment of PM10 enhancement by yellow sand on the air quality of Taipei, Taiwan in 2001. Environ MonitAssess 132: 297–309.
  53. 53. Chou CK, Lin CY, Chen TK, Hsu SC, Lung SC, et al. (2004) Influence of long-range transport dust particles on local air quality: A case study on Asian dust episodes in Taipei during the spring of 2002. Terr Atmos Ocean Sci 15: 881–899.
  54. 54. Hsu SC, Liu SC, Lin CY, Hsu RT, Huang YT, et al. (2002) Metal Compositions of PM10 and PM2.5 Aerosols in Taipei during Spring, 2002. Terr Atmos Ocean Sci 15: 925–948.
  55. 55. Hsu RN (2003) Effect of Air Pollution on Daily Hospital Admissions for Respiratory Diseases in Kaohsiung City, 1997–2001. Thesis, Fooyin University. Kaohsiung, Taiwan.