Skip to main content

Mortality associated with wildfire smoke exposure in Washington state, 2006–2017: a case-crossover study

Abstract

Background

Wildfire events are increasing in prevalence in the western United States. Research has found mixed results on the degree to which exposure to wildfire smoke is associated with an increased risk of mortality.

Methods

We tested for an association between exposure to wildfire smoke and non-traumatic mortality in Washington State, USA. We characterized wildfire smoke days as binary for grid cells based on daily average PM2.5 concentrations, from June 1 through September 30, 2006–2017. Wildfire smoke days were defined as all days with assigned monitor concentration above a PM2.5 value of 20.4 μg/m3, with an additional set of criteria applied to days between 9 and 20.4 μg/m3. We employed a case-crossover study design using conditional logistic regression and time-stratified referent sampling, controlling for humidex.

Results

The odds of all-ages non-traumatic mortality with same-day exposure was 1.0% (95% CI: − 1.0 - 4.0%) greater on wildfire smoke days compared to non-wildfire smoke days, and the previous day’s exposure was associated with a 2.0% (95% CI: 0.0–5.0%) increase. When stratified by cause of mortality, odds of same-day respiratory mortality increased by 9.0% (95% CI: 0.0–18.0%), while the odds of same-day COPD mortality increased by 14.0% (95% CI: 2.0–26.0%). In subgroup analyses, we observed a 35.0% (95% CI: 9.0–67.0%) increase in the odds of same-day respiratory mortality for adults ages 45–64.

Conclusions

This study suggests increased odds of mortality in the first few days following wildfire smoke exposure. It is the first to examine this relationship in Washington State and will help inform local and state risk communication efforts and decision-making during future wildfire smoke events.

Peer Review reports

Background

Wildfires are increasing in the western United States during the summer and fall months, emphasizing the importance of understanding the health impacts of wildfire smoke exposure [1, 2]. It is estimated that the total forest fire area burned in the western U.S. nearly doubled during 1984–2015 compared to the area projected to have burned without climate change [3]. This trend is expected to worsen, with climate projections indicating wildfires in the western U.S. will increase in frequency and intensity [4, 5]. The Intergovernmental Panel on Climate Change (IPCC) estimates that climate change will increase the length of wildfire season in North America by 10–30% [6], which is expected to result in worsening air quality during wildfire season in the coming decades [7].

Wildfire smoke contains a wide range of compounds known to be harmful to human health, including fine particulate matter (PM2.5), acrolein, benzene, carbon monoxide, and polycyclic aromatic hydrocarbons [8, 9]. Exposure to these toxic compounds is of concern near the source, and extending several hundred to thousands of kilometers away [10,11,12]. While it has been shown that the toxic compounds from wildfire smoke travel long distances from the source, potentially exposing thousands of individuals, the health effects associated with wildfire smoke exposure are just beginning to be understood [13].

Exposure assessment is challenging, as there is no standard approach for defining what constitutes a wildfire smoke-affected day or period in the health effects literature. Common methods utilize area monitoring particulate matter (PM) measurements, satellite data, chemical transport models, or a combination of these approaches [14,15,16,17]. Differences in effect estimates across studies in part may be due to differences in their exposure assessment approaches, limiting useful comparison within the growing published literature of wildfire smoke epidemiology. More research to define wildfire smoke-affected time periods is needed to validate methods currently in use.

Using the above exposure assessment methods, evidence for an association between exposure to wildfire smoke and all-age mortality is mixed. Four of five studies reported small increases in wildfire smoke associated mortality among all ages, however, all confidence intervals included the null [18,19,20,21]. The fifth study reported positive odds ratios (ORs) on lag days 0, 2, 3, and 4, but not on lag day 1 [22]. A study in Finland reported a 0.8% (95% CI: −3.5 – 5.3%) increase in all-cause mortality per 10 μg/m3 increase in same-day PM2.5 concentration [19], while another study in Australia estimated that wildfire events were associated with a 2.0% (95% CI: −2.0 – 5.0%) increase in the odds of same-day non-accidental mortality [21]. Several studies also report estimates for respiratory and cardiovascular mortality, finding limited evidence for an association with respiratory mortality [18, 21, 23] and small increases in the risk of cardiovascular mortality [18, 21].

A subset of these studies report larger effect estimates in groups 65 years and older [17, 24,25,26]. Analitis et al. reported that the effect of respiratory mortality in Greece was greater in adults age 75 and over during large fires [24]. Further, Haikerwal et al. observed an increase in risk of cardiac arrests, especially in older adults in Australia, although not all resulted in death [17], and Nunes et al. reported that older adults in their study in Brazil had the strongest association between exposure to biomass burning and circulatory disease mortality [26]. However, there are few U.S.-based mortality studies, providing little evidence for U.S.-specific mortality associated with wildfire smoke exposure. Additionally, there are no published mortality studies in Washington State, necessitating further research on the association between wildfire smoke exposure and non-traumatic mortality in the region.

Only two efforts have examined the health effects of wildfire smoke in Washington State, the 13th most populous state in the U.S., with an estimated 7.5 million people in 2018 [27]. Both of these efforts focused on the wildfire season of 2012, which was concentrated in Central Washington, resulting in limited population exposure [15, 28], and both studies examined hospitalizations, emergency department (ED) visits, and outpatient visits. No studies have been conducted on the risk of mortality associated with wildfire smoke exposure in Washington. Our study examines the association between wildfire smoke exposure and the odds of non-traumatic mortality in Washington State over 12 years. We hypothesized a priori that we would find a same-day effect and an effect at a lag of one or more days, within respiratory and cardiovascular mortality, and among individuals age 65 and older.

Methods

Mortality data

We conducted this study using Washington State geo-coded mortality and exposure data for wildfire season, defined as June through September, for our study period, 2006 to 2017. Historically, wildfire smoke has been documented in Washington during these months, and is described in an interagency operating plan as “peak fire season” in the Pacific Northwest [29]. Washington geo-coded mortality data includes latitude and longitude of residence, underlying cause of death, date of death, and decedent age, sex, and race. We examined non-traumatic causes of death (ICD-10 codes: A01-V99), including cardiovascular (ICD-10 codes: I05-I52); respiratory (ICD-10 codes: J01-J99); and cerebrovascular causes (ICD-10 codes: I60–67). About 2.7% of total cases were excluded due to missing latitude and longitude data (see Additional file 1: Figure S1).

We obtained median household income data at the census tract level from the U.S. Census Bureau for 2010–2017 and extrapolated to 2006 based on the five-year average percent change in growth within each census tract [30]. For the median household income analysis only, we omitted 12.9% of cases due to missing or incomplete household income data (see Additional file 1: Figure S1).

Exposure data

Common approaches to assessment of outdoor PM2.5 exposure include assigning subjects’ data from their nearest stationary monitor and a variety of modeling and interpolation techniques [31]. The “nearest monitor” approach is not adequate in Washington due to the complexity of Pacific Northwest meteorology and terrain [32]. Since accurate daily particulate matter concentrations were required, modeled or interpolated products alone could not provide sufficient daily accuracy for such an analysis, particularly during wildfire smoke conditions [33].

For the exposure assessment, we used 4 × 4 km grid cells from the Air Indicator Report for Public Awareness and Community Tracking (AIRPACT-4) model domain [34]. We assigned each grid cell in the domain to one of three regulatory air quality monitors closest to it, of the 75 regulatory air quality monitors in Washington (see Additional file 1: Figure S2), using an interpolated surface of summer mean PM2.5 concentrations in smoke-free conditions. This method matched grid cells to nearby monitors based on typical agreement between the interpolated PM2.5 at the grid cell and the monitor. Grid cells were assigned to a secondary or tertiary monitor on days when data at the primary monitor were not available, which decreased data loss due to missing exposure data (see Additional file 1: Text S1, for additional detail). We then assigned each grid cell the daily PM2.5 concentration from its assigned monitor on that day and meteorological variables from the monitor’s nearest National Weather Service meteorological site. With this approach, we could leverage the accuracy of daily measurement data while improving upon the typical approach of simply assigning the nearest monitor. The result is a dataset with the following for each day and each grid cell: 24-h average PM2.5 concentration and humidex, a measure of apparent temperature calculated from air temperature and dew point [35].

About 1.5% of person-days across the study time period were not considered to be represented by any monitoring site, and were excluded from the analysis (see Additional file 1: Figure S3, for a map showing areas excluded). This resulted in the exclusion of 1.2% of all non-traumatic deaths in Washington during the study period (see Additional file 1: Figure S1). We joined this dataset with the above described mortality data using a spatial join in ArcGIS (version 10.5.1; Esri, Redlands, CA), assigning the latitude and longitude of the residence of each decedent to the nearest grid cell and corresponding PM2.5 concentration and humidex value.

Wildfire smoke day classification

In order to identify wildfire smoke-affected days, we considered a number of approaches. First, a statewide PM2.5 concentration was tested as the threshold between wildfire smoke and non-wildfire smoke days, set at 15 μg/m3, corresponding to the 99th percentile of measured PM2,5 concentrations across two relatively smoke free years. The majority of days that exceeded 15 μg/m3 were recorded in urban areas without wildfire smoke contributions. Thus, a more nuanced and area-specific approach was needed to minimize false positives in urban areas with higher background particulate matter and to minimize false negatives in rural areas with lower background particulate matter [36, 37]. We defined wildfire smoke-affected days as grid cell days with a 24-h average PM2.5 concentration greater than 20.4 μg/m3, with an additional set of criteria for days between 9 and 20.4 μg/m3. A concentration of 20.4 μg/m3 corresponds to the threshold between Moderate and Unhealthy for Sensitive Groups Air Quality by the Washington Air Quality Advisory [38], and we found background anthropogenic particulate matter across the study period to be below this level. For days with a 24-h average PM2.5 concentration between 9 and 20.4 μg/m3, we applied the following criteria:

  1. 1)

    The day must be part of an event in which at least 2 of 3 consecutive days are greater than 9 μg/m3;

  2. 2)

    One of the days in the 3-day event window must be greater than 15 μg/m3;

  3. 3)

    For urban areas (Seattle, Tacoma, Spokane), at least 50% of the air monitors in those areas must be greater than 9 μg/m3 (see Additional file 1: Text S2, for a definition of urban areas).

We selected these criteria due to the nature of wildfire smoke events in Washington State. The first two criteria were informed by the historical observation that nearly all smoke events span multiple days, and the third criterion was informed by the observation that smoke events tend to affect nearby monitors in a region.

Statistical analysis

We employed a time-stratified case-crossover design to examine the association between wildfire smoke exposure and non-traumatic mortality, using conditional logistic regression. This study design compares wildfire smoke exposure, defined as the binary wildfire smoke day classification described above, on the day of death, the day prior to death, and on the 4 days prior to death, to wildfire smoke exposure on referent (non-event, or control) days for the same decedent. We selected referent days using time-stratified sampling, where we defined the strata as the same day of the week, month, and year of death, yielding 3.39 referent days per decedent, on average. By design, this technique controls for time-invariant confounders, including sex, age, race, pre-existing health conditions, and other individual characteristics and risk factors, as each person serves as their own control [39]. This design also controls for some time-dependent variables based on the referent selection method, including day of the week, and seasonal trends in air pollution [40]. We reported results as a percent change in the odds of non-traumatic mortality for wildfire smoke-affected days versus non-wildfire smoke-affected days, after controlling for humidex. We adjusted for humidex by adding a term into the conditional logistic regression estimating equation. We used the clogit function in the survival R package to conduct the regression analysis [41].

To examine the effect of wildfire smoke exposure by characteristics of interest, we conducted subgroup analyses, stratifying by sex, age group, race category, cause of death, location (urban and rural), and census tract median household income, stratified by income groups shown in Table 1, for both same-day and previous-day exposures. We also conducted a lag analysis using an unconstrained distributed lag model from days 0–4, with day 0 modeled as the day of death, lag day 1 as the previous day, and so on.

Table 1 Non-traumatic mortality characteristics

We conducted secondary subgroup analyses for same-day exposures. We reported estimates by age group and race within respiratory causes of death, and by age group for Chronic Obstructive Pulmonary Disease (COPD) causes of death (Table 4). We also conducted a sensitivity analysis, setting a PM2.5 concentration of 20.4 μg/m3 as the wildfire smoke-affected day threshold, without additional criteria for days between 9 and 20.4 μg/m3, to assess whether our results were sensitive to the exposure definition (see Additional file 1: Table S1).

All analyses were performed using R 3.4.3 [41].

Results

Table 1 summarizes characteristics of the 170,985 non-traumatic deaths included in our study in Washington from June through September for 2006–2017. Most non-traumatic deaths occurred in those 65 years and older (76.6%), most were white race (90.2%), and most lived in census tracts with a median household income of less than $75,000 (70.5%). About a quarter of deaths were due to cardiovascular causes (26.0%) and less than 10% were due to respiratory causes (7.8%). Table 1 also reports the number and percent of deaths that contribute to the inferential analysis, defined as belonging to a stratum with exposure variation, i.e. containing both wildfire smoke and non-wildfire smoke days.

In Table 2, we display exposure characteristics across the study period, including mean PM2.5 and humidex on wildfire smoke days and non-wildfire smoke days, and average PM2.5 on event and referent days.

Table 2 Daily PM2.5 characteristics for mortality days and referent days

In Fig. 1, we display the results of our unconstrained distributed lag analysis, examining the effect of exposure to wildfire smoke in the 4 days prior to death. The pattern suggests that previous day exposure conveys the highest risk and that it diminishes rapidly such that there is no evidence of increased risk after 2 days. We estimate a 1.3% (95% CI: 0.2–2.4%) increase in the odds of non-traumatic mortality given wildfire smoke exposure on the previous day, while holding constant the humidex and exposure on lag days 0 and 2–4 (Fig. 1). The results indicate some evidence for an effect of exposure at 2 days prior to death, with inconclusive evidence for an effect of exposure in the preceding days on death.

Fig. 1
figure 1

ORs of non-traumatic mortality on wildfire smoke days compared to non-wildfire smoke days from an unconstrained distributed lag model by lag day, adjusted by humidex

In Table 3, we present the results of the inferential analysis, indicating a 1.0% (95% CI: − 1.0 – 4.0%) increase in the odds of all same-day non-traumatic mortality on wildfire smoke days compared to non-wildfire smoke days, controlling for humidex. We further observed a 1.0% (95% CI: − 6.0 – 4.0%) decrease in the odds of same-day cardiovascular mortality, a 9.0% (95% CI: 0.0–18.0%) increase in the odds of same-day respiratory mortality, and a 14.0% (95% CI: 2.0–26.0%) increase in the odds of same-day COPD mortality. Among ages 65–84, we observed a 2.0% (95% CI: − 2.0 – 6.0%) increase in the odds of all same-day non-traumatic mortality. With previous day exposure, we observed a 2.0% (95% CI, 0.0–5.0%) increase in the odds of all non-traumatic mortality, and a 5.0% (95% CI: − 3.0 – 15.0%) increase in the odds of respiratory mortality.

Table 3 Odds ratios (ORs) and 95% confidence intervals (CIs) for non-traumatic mortality for lag day 0 and lag day 1

We conducted additional same-day and previous-day analyses, stratifying by age group, location, census tract median household income, race, and sex (Table 3). The majority of these analyses indicate little evidence for an effect, with odds ratios near 1.0 and confidence intervals consistent with null effects. However, when stratifying by race category, we observed a 19.0% (95% CI: −33.0 – −1.0%) decrease in the odds of same-day non-traumatic mortality among Hispanics, on wildfire smoke days compared to non-wildfire smoke days. Stratifying by census tract median household income, we observed a 7.0% (95% CI: 2.0–13.0%) increase, and a 7.0% (95% CI: 0.0–15.0%) increase in the odds of non-traumatic mortality with previous day exposure among those living in a census tract with median household income of $35,000 - $50,000 and $75,000–$100,000, respectively. We also observed a 3.0% (95% CI: 0.0–6.0) increase and a 21.0% (95% CI: −36.0 – −2.0%) decrease in non-traumatic mortality with previous day exposure among Whites and Native Americans, respectively. It is important to note, however, that many of the stratified groups have small sample sizes, limiting the power to detect an effect. Thus, the effect estimates should be interpreted with an appropriate degree of caution.

We conducted secondary analyses stratified by age group and race category within respiratory causes of death with same-day exposure (Table 4). For ages 45–64, we observed a 35.0% (95% CI: 9.0–67.0%) increase in the odds of same-day respiratory mortality, and a 33.0% (95% CI: 0.0–78.0%) increase in the odds of same-day COPD mortality. Within all respiratory mortality, we observed a 12.0% (95% CI: 2.0–22.0%) increase in the odds of same-day respiratory mortality among Whites.

Table 4 Odds ratios (ORs) and 95% confidence intervals (CIs) for stratification of same-day respiratory causes of death

Discussion

Overall, our findings suggest evidence of a small increase in odds of all-ages non-traumatic mortality with same-day and with previous-day wildfire smoke exposure. Other mortality studies also provide evidence of an increase in the risk of non-accidental or natural mortality with same-day and previous day wildfire smoke exposure [18, 21, 23], including Johnston et al., that report a 1.0% (95% CI: −2.0 – 5.0%) and a 5.0% (95% CI: 0–10.0%) increase in the odds of non-accidental mortality with same-day and previous-day wildfire smoke exposure, respectively [21].

From our distributed lag model, we observed a 1.3% (95% CI: 0.0–2.0%) increase in the odds of next day all-ages, non-traumatic mortality for previous day wildfire smoke exposure versus non-wildfire smoke exposure, controlling for both humidex and exposure for the 4 days prior to death. Other studies report evidence for a lagged effect of non-traumatic mortality from wildfire smoke exposure from 1 to 5 days prior to death [18, 21], although few studies employ a distributed lag model [18]. Future wildfire smoke studies should employ distributed lag models, as is commonplace in air pollution epidemiology, to better capture the multiday lingering effects of exposure on health outcomes [42].

In examining specific causes of death, we observed estimates for all-ages cardiovascular mortality that are consistent with either an increase or decrease in risk, with same-day or previous-day wildfire smoke exposure. Other studies find evidence of an increased risk of all-ages cardiovascular mortality with both same-day and previous-day wildfire smoke exposure [21, 23]. However, Morgan et al. employed different exposure metrics, using PM10 rather than PM2.5 and reporting risk for each 10 μg/m3 increase in PM10, rather than comparing smoke days to non-smoke days [23]. Thus, we cannot easily compare risk estimates across studies.

Many of our subgroup analyses are limited by their more exploratory nature and wide confidence intervals, and thus must be interpreted with caution. However, we discuss two subgroup analyses that merit additional research to better understand the risk to each population.

In one subgroup analysis, we found evidence of an increase in the odds of all-ages COPD mortality with same-day and with previous-day wildfire smoke exposure. To our knowledge, no studies examine the association between wildfire smoke exposure and COPD mortality. However, populations with underlying health conditions, and in particular, asthma and COPD, have been found to be more susceptible to wildfire smoke compared to healthy populations in several studies examining hospital admissions and ED visits [23, 43,44,45]. Further research into the association between wildfire smoke exposure and COPD mortality is needed to confirm this observation.

In an additional subgroup analysis, we reported evidence of an increase in the odds of respiratory mortality among individuals ages 45–64 with same-day exposure. Several papers find higher effects of wildfire smoke exposure among adults under 65 compared to adults over 65, albeit for different health endpoints, none examining mortality as an endpoint for any age group [45,46,47,48]. Henderson et al. report the largest ORs for respiratory physician visits among adults ages 30–40 year [46], while Mott et al. find asthma and COPD hospital admissions greatest among adults ages 40–64 [47]. Rappold et al. find higher increases in respiratory ED visits (including asthma, COPD, pneumonia, and acute bronchitis) among individuals under 65 compared to those over 65 [48], and Reid et al. report higher COPD ED visits associated with PM2.5 in those ages 20–64 compared to over 65 [45]. However, most of the effect estimates we and others report are from secondary subgroup analyses with wide confidence intervals, meriting cautious interpretation. Nonetheless, people 45–64 with COPD are less likely to be using oxygen than people 65 and over with COPD, and thus more likely to be mobile and exposed to wildfire smoke [49]. Our findings suggest that underlying respiratory health conditions may contribute to the increased risk of respiratory and COPD-related mortality in the 45–64 age group. Thus, we recommend risk messaging target those of all ages with underlying health conditions, and specifically respiratory and COPD health conditions. Additional research is needed to further examine the risk to this population.

This analysis is limited by the challenges of separating anthropogenic PM2.5 and wildfire smoke PM2.5. We put considerable effort into determining a viable threshold between wildfire smoke affected PM2.5 and non-wildfire smoke affected PM2.5, but it is likely some misclassification exists [50]. We conducted a sensitivity analysis using 20.4 μg/m3 as the wildfire smoke-affected day threshold (see Additional file 1: Table S1), which showed the effects from wildfire smoke are sensitive to the exposure definition, and that the effect seen in this study may be capturing both the risk due to anthropogenic PM2.5 and wildfire smoke PM2.5. Further, existing wildfire smoke and mortality studies employ a wide variety of exposure assessment methods, which impedes direct comparison of effect estimates across studies. Additional research is needed to develop an ideal method to identify wildfire smoke-affected time periods.

Additional study limitations include the assumption that PM2.5 concentrations at each monitor represent the exposure for the population attributed to each monitor. Further, our method of assigning exposure allows individuals to be assigned to a different monitor on each day if a monitor is not working or is malfunctioning. This method, while not optimal, recoups data that would otherwise be dropped. Additionally, the air quality monitors do not represent the true exposure experienced by each person in the monitor grid cell area, likely resulting in some degree of exposure misclassification, a common limitation in air pollution epidemiology studies [51]. Further, exposures assigned from ambient PM2.5 concentrations do not reflect the reality that most people spend about 90% of their lives indoors [52], and additionally do not account for the steps people take to reduce exposures, where the main public health guidance during smoke episodes is to go indoors and keep indoor air clean [53].

Another limitation is that Washington State geo-coded mortality data may not correspond to the location of exposure for some decedents. We linked the geo-coded location of each decedent’s home residence to the corresponding grid cell. If the location of residence was not the location of exposure, this may result in some misclassification of exposure.

A final limitation of our study is regarding median household income, where an ecological indicator was employed, assigning the census tract median household income to each decedent. This method misclassifies individual income, but serves as a proxy for neighborhood-level socioeconomic status. Due to inadequate median household income data for 2006–2009, estimates for that time period are prone to higher rates of misclassification. Further analyses should develop more accurate proxies for individual or neighborhood-level socioeconomic status.

Conclusion

This study is the first to estimate mortality risk associated with wildfire smoke exposure in Washington State. This study uses a tiered approach to exposure assessment, minimizing false allocation of anthropogenic PM-dominated days in urban areas as wildfire smoke days, as well as false allocation of smoke-dominated days in rural areas as anthropogenic PM-dominated days. This work will support local and state risk communication efforts and decision-making during future wildfire smoke events, especially for susceptible subpopulations identified in this study: those of all ages with COPD and other underlying respiratory conditions. Additional research is needed in Washington State to characterize the association between wildfire smoke exposure and less severe health endpoints of interest, including hospitalizations and ED visits, the health effects among vulnerable populations, as well as the health effects of prolonged smoke exposure.

Availability of data and materials

The mortality data that support these findings are publicly available for a fee from the Washington State Department of Health.

The air quality data that support these findings are publicly available from the Air Quality Program at the Washington State Department of Ecology upon request.

Abbreviations

AIRPACT:

Air Indicator Report for Public Awareness and Community Tracking model

CI:

Confidence interval

COPD:

Chronic Obstructive Pulmonary Disease

ED:

Emergency Department

ICD:

International Classification of Diseases

OR:

Odds ratio

PM2.5 :

Particulate matter less than 2.5 μm in diameter

SD:

Standard deviation

References

  1. Dennison PE, Brewer SC, Arnold JD, Moritz MA. Geophysical research letters. Geophys Prospect. 2014;41:2928–33.

    Google Scholar 

  2. Westerling AL. Increasing western US forest wildfire activity: sensitivity to changes in the timing of spring. Philos Trans R Soc B. 2016;371:20150178.

    Article  Google Scholar 

  3. Abatzoglou JT, Williams AP. Impact of anthropogenic climate change on wildfire across western US forests. PNAS. 2016;113(42):11770–175.

    Article  CAS  Google Scholar 

  4. Joyce L, Running S, Breshears D, Dale V, Malmsheimer R, Sampson R, et al. Ch. 7: Forests. Climate change impacts in the United States: the third national climate assessment. 2014.

    Book  Google Scholar 

  5. Littell J, McKenzie D, Peterson D, Westerling A. Climate and wildfire area burned in western U.S. Ecol Appl. 2009;19(4):1003–21.

    Article  Google Scholar 

  6. Parry ML. Technical summary. Climate change 2007: impacts, adaptation, and vulnerability. Cambridge: Cambridge University Press; 2007.

    Google Scholar 

  7. Interagency Working Group on Climate Change and Health. A human health perspective on climate change: a report outlining the research needs on the human health effects of climate change. 2010.

    Google Scholar 

  8. Naeher LP, Brauer M, Lipsett M, Zelikoff JT, Simpson CD, Koenig JQ, et al. Woodsmoke health effects: a review. Inhal Toxicol. 2010;19(1):67–106.

    Article  CAS  Google Scholar 

  9. Urbanski S, Hao W, Baker S. Chapter 4: chemical composition of wildland fire emissions. Dev Environ Sci. 2008;8:79–107.

    Google Scholar 

  10. Moeltner K, Kim M, Zhu E, Yang W. Wildfire smoke and health impacts: a closer look at fire attributes and their marginal effects. J Environ Econ Manag. 2013;66(3):476–96.

    Article  Google Scholar 

  11. Munoz-Alpizar R, Pavlovic R, Moran MD, Chen J, Gravel S, Henderson SB, et al. Exposure over North America as determined from operational air quality forecasts. Atmosphere (Basel). 2017;8:179.

    Article  CAS  Google Scholar 

  12. Kollanus V, Tiittanen P, Niemi JV, Lanki T. Effects of long-range transported air pollution from vegetation fires on daily mortality and hospital admissions in the Helsinki metropolitan area, Finland. Environ Res. 2016;151:351–8.

    Article  CAS  Google Scholar 

  13. Black C, Tesfaigzi Y, Bassein JA, Miller LA. Wildfire smoke exposure and human health: significant gaps in research for a growing public health issue. Environ Toxicol Pharmacol. 2017;55(August):186–95.

    Article  CAS  Google Scholar 

  14. Alman BL, Pfister G, Hao H, Stowell J, Hu X, Liu Y, et al. The association of wildfire smoke with respiratory and cardiovascular emergency department visits in Colorado in 2012: a case crossover study. Environ Heal A Glob Access Sci Source. 2016;15(1):1–9.

    Google Scholar 

  15. Gan RW, Ford B, Lassman W, Pfister G, Vaidyanathan A, Fischer E, et al. Associations with cardiopulmonary-related hospital admissions. Geohealth. 2017;1(3):122–36.

    Article  Google Scholar 

  16. Yao J, Eyamie J, Henderson SB. Evaluation of a spatially resolved forest fire smoke model for population-based epidemiologic exposure assessment. J Expo Sci Environ Epidemiol. 2016;26(3):233–40.

    Article  CAS  Google Scholar 

  17. Haikerwal A, Akram M, Del MA, Smith K, Sim MR, Meyer M, et al. Impact of fine particulate matter (PM2.5) exposure during wildfires on cardiovascular health outcomes. J Am Heart Assoc. 2015;4(7):1–10.

    Article  Google Scholar 

  18. Faustini A, Alessandrini ER, Pey J, Perez N, Samoli E, Querol X, et al. Short-term effects of particulate matter on mortality during forest fires in southern Europe: results of the MED-PARTICLES project. Occup Environ Med. 2015;72(5):323–9.

    Article  Google Scholar 

  19. Hänninen OO, Salonen RO, Koistinen K, Lanki T, Barregard L, Jantunen M. Population exposure to fine particles and estimated excess mortality in Finland from an east European wildfire episode. J Expo Sci Environ Epidemiol. 2009;19(4):414–22.

    Article  CAS  Google Scholar 

  20. Hong KY, King GH, Saraswat A, Henderson SB. Seasonal ambient particulate matter and population health outcomes among communities impacted by road dust in British Columbia, Canada. J Air Waste Manage Assoc. 2017;67(9):986–99.

    Article  CAS  Google Scholar 

  21. Johnston F, Hanigan I, Henderson S, Morgan G, Bowman D. Extreme air pollution events from bushfires and dust storms and their association with mortality in Sydney, Australia 1994-2007. Environ Res. 2011;111(6):811–6.

    Article  CAS  Google Scholar 

  22. Kollanus V, Prank M, Gens A, Soares J, Vira J, Kukkonen J, et al. Mortality due to vegetation fire-originated PM2.5 exposure in Europe - assessment for the years 2005 and 2008. Enviromental Heal Perpectives. 2017;125(1):30–7.

    Article  CAS  Google Scholar 

  23. Morgan G, Sheppeard V, Khalaj B, Ayyar A, Lincoln D, Jalaludin B, et al. Effects of bushfire smoke on daily mortality and hospital admissions in Sydney, Australia. Epidemiology. 2010;21(1):47–55.

    Article  Google Scholar 

  24. Analitis A, Georgiadis I, Katsouyanni K. Forest fires are associated with elevated mortality in a dense urban setting. Occup Environ Med. 2012;69(3):158–62.

    Article  Google Scholar 

  25. Shaposhnikov D, Revich B, Bellander T, Bedada GB, Bottai M, Kharkova T, et al. Mortality related to air pollution with the Moscow heat wave and wildfire of 2010. Epidemiology. 2014;25(3):359–64.

    Article  Google Scholar 

  26. Nunes KVR, Ignotti E, Hacon SDS. Circulatory disease mortality rates in the elderly and exposure to PM2.5 generated by biomass burning in the Brazilian Amazon in 2005. Cad Saude Publica. 2013;29(3):589–98.

    Article  Google Scholar 

  27. U.S. Census Burea. Annual estimates of the resident population: April 1, 2010 to July 1, 2018. 2019. Available from: https://factfinder.census.gov/faces/tableservices/jsf/pages/productview.xhtml?src=bkmk

    Google Scholar 

  28. Washington State Department of Health. Surveillance investigation of the cardiopulmonary health effects of the 2012 wildfires in North Central Washington State 2012.

    Google Scholar 

  29. 2019 Northwest area fire weather annual operating plan. 2019. Available from: https://gacc.nifc.gov/nwcc/content/pdfs/weather/2019%20Region%206%20AOP_Final_p.pdf.

  30. U.S. Census Bureau. Selected economic characteristics: American comunity survey 5-year estimates. Available from: https://factfinder.census.gov/faces/tableservices/jsf/pages/productview.xhtml?pid=ACS_17_5YR_DP03&src=pt. Accessed 10 Mar 2019.

  31. Hoek G. Methods for assessing long-term exposures to outdoor air pollutants. Curr Envir Heal Rpt. 2017;4:450–62.

    Article  CAS  Google Scholar 

  32. Lassman W, Ford B, Gan RW, Fischer EV, Pierce JR. Spatial and temporal estimates of population exposure to wildfire smoke during the Washington state 2012 wildfire season using blended model, satellite, and in situ data. GeoHealth. 2017;1:106–21.

    Article  Google Scholar 

  33. Herron-Thorpe F, Mount G, Emmons L, Lamb B, Jaffe D, Wigder N, et al. Air quality simulations of wildfires in the Pacific northwest evaluated with surface and satellite observations during the summers of 2007 and. Atmos Chem Phys. 2008;2014:12533–51.

    Google Scholar 

  34. Washington State University. AIRPACT. 2019. Available from: http://lar.wsu.edu/airpact/

    Google Scholar 

  35. Masterton J, Richardson F. Humidex: a method of quantifying human discomfort due to excessive heat and humidity. Downsview: Environment Canada, Atmospheric Environment; 1979.

    Google Scholar 

  36. Dockery D, Pope C, Xu X, Spengler J, Ware J, Fay M, et al. An association between air pollution and mortality in six U.S. cities. N Engl J Med. 1993;329:1753–9.

    Article  CAS  Google Scholar 

  37. U. S. Environmental Protection Agency. Quantitative health risk assessment for particulate matter. U S Environ Prot Agency; 2010. p. 1–596. Available from: http://www.epa.gov/ttn/naaqs/standards/pm/data/PM_RA_FINAL_June_2010.pdf%5Cnpapers2://publication/uuid/36AC0D6B-29A0-4A0F-A8B5-699617A33A15

    Google Scholar 

  38. WA DOH; WA ECY. WAQA: Washington air quality advisory. Available from: https://fortress.wa.gov/ecy/publications/documents/1802024.pdf. Accessed 4 Apr 2019.

  39. Lumley T, Levy D. Bias in the case-crossover design: implications for studies of air pollution. Environmetrics. 2000;11:689–704.

    Article  CAS  Google Scholar 

  40. Janes H, Sheppard L, Lumley T. Case-crossover analyses of air pollution exposure data: referent selection strategies and their implications for bias. Epidemiology. 2005;16(6):717–26.

    Article  Google Scholar 

  41. R Development Core Team. R: A language and environment for statistical computing. Vienna: R Foundation for Statistical Computing; 2018. Available from: http://www.r-project.org

    Google Scholar 

  42. Schwartz J. The distributed lag between air pollution and daily deaths. Epidemiology. 2000;11(3):320–6.

    Article  CAS  Google Scholar 

  43. Delfino RJ, Brummel S, Wu J, Stern H, Ostro B, Lipsett M, et al. The relationship of respiratory and cardiovascular hospital admissions to the southern California wildfires of 2003. Occup Environ Med. 2009;66(3):1–24.

    Article  Google Scholar 

  44. Johnston FH, Purdie S, Jalaludin B, Martin KL, Henderson SB, Morgan GG. Air pollution events from forest fires and emergency department attendances in Sydney, Australia 1996-2007: a case-crossover analysis. Environ Heal A Glob Access Sci Source. 2014;13(1):1–9.

    Google Scholar 

  45. Reid CE, Jerrett M, Tager IB, Petersen ML, Mann JK, Balmes JR. Differential respiratory health effects from the 2008 northern California wildfires: a spatiotemporal approach. Environ Res. 2016;150:227–35.

    Article  CAS  Google Scholar 

  46. Henderson SB, Brauer M, Kennedy S, MacNab Y. Three measures of forest fire smoke exposure and their association with respiratory and cardiovascular physician visits and hospital admissions. Epidemiology. 2009;20(9):S82.

    Article  Google Scholar 

  47. Mott JA, Mannino DM, Alverson CJ, Kiyu A, Hashim J, Lee T, et al. Cardiorespiratory hospitalizations associated with smoke exposure during the 1997 southeast Asian forest fires. Int J Hyg Environ Health. 2005;208(1–2):75–85.

    Article  Google Scholar 

  48. Rappold AG, Stone SL, Cascio WE, Neas LM, Kilaru VJ, Carraway MS, et al. Peat bog wildfire smoke exposure in rural North Carolina is associated with cardiopulmonary emergency department visits assessed through syndromic surveillance. Enviromental Heal Perpectives. 2011;119(10):1415–20.

    Article  Google Scholar 

  49. Maselli DJ, Bhatt SP, Anzueto A, Bowler RP, DeMeo DL, Diaz AA, Dransfield MT, Fawzy A, Foreman MG, Hanania NA, Hersh CP, Kim V, Kinney GL, Putcha N, Wan ES, Wells JM, Westney GE, Young KA, Silverman EKBJ. Clinical epidemiology of COPD – Insights from 10 years of the COPDGene study. Chest. 2019;156(2):228–38. https://doi.org/10.1016/j.chest.2019.04.135.

    Article  Google Scholar 

  50. Doubleday A. Mortality associated with wildfire smoke exposure in Washington state, 2006–2017. University of Washington; 2019.

    Google Scholar 

  51. Nieuwenhuijsen M. Exposure assessment in environmental epidemiology. New York: Oxford University Press; 2015.

    Book  Google Scholar 

  52. U.S. EPA. Report to congress on indoor air quality: volume 2. EPA/400/1–89/001C. Washington, DC. 1989.

    Google Scholar 

  53. U.S. EPA, U.S. Forest Service, U.S. Centers for disease control and prevention, California air resources board. Wildfire smoke: a guide for public health officials, revised May 2016; 2016. p. 1–27.

    Google Scholar 

Download references

Acknowledgements

The authors would like to thank the Washington State Department of Health for use of geo-coded mortality data.

Funding

This research did not receive any specific grant from funding agencies in the public, commercial, or not-for-profit sectors.

Author information

Authors and Affiliations

Authors

Contributions

AD: Conceptualization, Methodology, Software, Validation, Formal Analysis, Investigation, Data Curation, Writing – Original Draft, Writing – Review & Editing, Visualization. JS: Conceptualization, Methodology, Software, Validation, Data Curation, Writing – Review & Editing. LS: Methodology, Formal Analysis, Writing – Review & Editing. MK: Conceptualization, Methodology, Writing – Review & Editing. RD: Conceptualization, Methodology, Data Curation, Writing – Reviewing & Editing. JF: Conceptualization, Data Curation, Writing – Review & Editing. TBI: Conceptualization, Methodology, Investigation, Writing – Review & Editing, Supervision, Project Administration. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Annie Doubleday.

Ethics declarations

Ethics approval and consent to participate

The University of Washington Institutional Review Board reviewed the study design and determined the study to be exempt. The Washington State Institutional Review Board determined that we did not need to submit a form for official determination of exempt status for use of Washington State mortality data for this study.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary information

Additional file 1:

Figure S1. Number of cases of non-traumatic mortality at each stage of the study. Figure S2. Locations of the 75 regulatory air quality monitors in Washington State. Text S1. Exposure grid methods. Figure S3. Number of years with monitored exposure data for each 4 × 4 km grid cell. Text S2. Classification of monitors. Table S1. ORs and 95% CIs for all-ages, all same-day non-traumatic mortality associated with a > 20.4 μg/m3 threshold defining wildfire smoke days.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Doubleday, A., Schulte, J., Sheppard, L. et al. Mortality associated with wildfire smoke exposure in Washington state, 2006–2017: a case-crossover study. Environ Health 19, 4 (2020). https://doi.org/10.1186/s12940-020-0559-2

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12940-020-0559-2

Keywords