Estimated Global Mortality Attributable to Smoke from Landscape Fires
Background: Forest, grass, and peat fires release approximately 2 petagrams of carbon into the atmosphere each year, influencing weather, climate, and air quality.
Objective: We estimated the annual global mortality attributable to landscape fire smoke (LFS).
Methods: Daily and annual exposure to particulate matter ≤ 2.5 μm in aerodynamic diameter (PM2.5) from fire emissions was estimated globally for 1997 through 2006 by combining outputs from a chemical transport model with satellite-based observations of aerosol optical depth. In World Health Organization (WHO) subregions classified as sporadically affected, the daily burden of mortality was estimated using previously published concentration–response coefficients for the association between short-term elevations in PM2.5 from LFS (contrasted with 0 μg/m3 from LFS) and all-cause mortality. In subregions classified as chronically affected, the annual burden of mortality was estimated using the American Cancer Society study coefficient for the association between long-term PM2.5 exposure and all-cause mortality. The annual average PM2.5 estimates were contrasted with theoretical minimum (counterfactual) concentrations in each chronically affected subregion. Sensitivity of mortality estimates to different exposure assessments, counterfactual estimates, and concentration–response functions was evaluated. Strong La Niña and El Niño years were compared to assess the influence of interannual climatic variability.
Results: Our principal estimate for the average mortality attributable to LFS exposure was 339,000 deaths annually. In sensitivity analyses the interquartile range of all tested estimates was 260,000–600,000. The regions most affected were sub-Saharan Africa (157,000) and Southeast Asia (110,000). Estimated annual mortality during La Niña was 262,000, compared with 532,000 during El Niño.
Conclusions: Fire emissions are an important contributor to global mortality. Adverse health outcomes associated with LFS could be substantially reduced by curtailing burning of tropical rainforests, which rarely burn naturally. The large estimated influence of El Niño suggests a relationship between climate and the burden of mortality attributable to LFS.
Landscape fires (encompassing wild and prescribed forest fires, tropical deforestation fires, peat fires, agricultural burning, and grass fires) release approximately 2 petagrams (2 × 1012 kg) of carbon into the atmosphere annually (van der Werf et al. 2010). These emissions affect planetary processes such as radiative forcing (which influences average global temperatures) and hydrological cycles (which influence regional cloud formation and rainfall) (Bowman et al. 2009; Cochrane and Laurance 2008; Fargione et al. 2008; Langmann et al. 2009; Tosca et al. 2010; Yokelson et al. 2007). Most emissions originate from fires set in tropical rainforests and savannas, where they cause recurrent episodes of severe pollution that affect some of the poorest regions of the world (van der Werf et al. 2010). Despite extensive literature describing the harmful effects of air pollution, the health impacts of landscape fire smoke (LFS) are rarely highlighted in discussions about fires and their role in the earth system (Lohman et al. 2007).
Smoke from the combustion of biomass is composed of hundreds of chemicals, many of which are known to be harmful to human health (Naeher et al. 2007). The most important risk-related measure of smoke is particulate matter (PM) with an aerodynamic diameter ≤ 2.5 μm (PM2.5). This PM primarily consists of organic carbon and black carbon components, along with smaller contributions from inorganic species (Naeher et al. 2007; Reid et al. 2005). PM is also produced by the combustion of fossil fuels, and most health evidence for PM2.5 comes from studies in urban environments (Pope and Dockery 2006). Urban PM has been associated with a wide range of adverse health outcomes including all-cause, neonatal and cardiorespiratory mortality, exacerbations of respiratory and cardiovascular conditions, and pathophysiological changes such as inflammation, oxidative stress, and procoagulation (Pope and Dockery 2006). The effects of PM derived from burning biomass have been less extensively investigated, and much of the evidence comes from studies of air pollution from household solid fuel use (Naeher et al. 2007). A handful of toxicological studies suggest that biomass smoke particles elicit pathophysiological effects similar to those of urban PM (Barregard et al. 2006; Danielsen et al. 2009; Kocbach et al. 2008). Although there are relatively few epidemiological studies on smoke-related PM, they also report outcomes consistent with those elicited by urban PM, including increased all-cause mortality and exacerbations of respiratory conditions (Delfino et al. 2009; Hänninen et al. 2009; Johnston et al. 2007, 2011; Morgan et al. 2010; Sastry 2002). However, evidence concerning cardiovascular outcomes of smoke-related PM remains scarce and inconclusive (Naeher et al. 2007; Sanhueza et al. 2009). Results from several studies of the extensive rainforest and peat fires in Southeast Asia in 1997 through 1998 suggest substantial health and economic impacts of LFS (Jayachandran 2009; Mott et al. 2005; Sastry 2002; Schweithelm et al. 2006). Further, fires are becoming more widespread and frequent in some regions (Turetsky et al. 2011; Westerling et al. 2006), and this source of air pollution is likely to continue to grow in magnitude and consequent health impacts (Confalonieri et al. 2007; Denman et al. 2007; Langmann et al. 2009). Because fire emissions contribute to radiative forcing, there is potential for the development of a positive feedback between a warming climate and increasingly severe fire events in several biomes (Bowman et al. 2009). In this context, a global assessment of the mortality impacts of LFS is required.
Materials and Methods
Studying the magnitude of health impacts from LFS presents several technical challenges, including estimation of the exposure to smoke-specific PM for each spatial unit of analysis, selection of the most appropriate concentration–response functions, and consideration of what theoretic minimum (counterfactual) exposure values to apply. Moderate to high levels of uncertainty are associated with many of these steps, so our objectives were to provide a reasonable principal estimate given the available data and then to evaluate the sensitivity of the principal estimate to the assumptions used in the principal analysis. The World Health Organization (WHO) Global Burden of Disease (GBD) Comparative Risk Assessment framework provides a standard set of methods for this and has previously been used to evaluate the annual mortality attributable to urban air pollution and to indoor air pollution from household solid fuel use (Ezzati et al. 2002; Lopez et al. 2006a). Methods for estimating the global mortality associated with particulate air pollution are being revised in the light of new epidemiological evidence and exposure assessment methods, and new cause-specific results are expected in 2012 (Institute for Health Metrics and Evaluation 2010). However, the epidemiological evidence concerning LFS remains limited, and evidence concerning LFS and cause-specific mortality is not currently available. For this reason, our analyses evaluate all-cause mortality.
Input data. Exposure estimates. We combined information from satellite-derived observations of global fire activity, geographic area burned, and type of vegetation burned in a global atmospheric three-dimensional (3-D) chemical transport model. We then combined output from that model with satellite-based measurements of aerosol optical depth (AOD) to estimate annual PM2.5 emissions from landscape fires. For a detailed description of the exposure estimates, see Supplemental Material, pp. 3–8, Table 1, and Figures 1 and 2 ( http://dx.doi.org/10.1289/ehp.1104422). A summary is presented below.
|Scenario||Global||Sub-Saharan Africa>a||Southeast Asiab||South Americac|
|Annual average (1997–2006)||339,000||157,000||110,000||10,000|
|EL Niño year (September 1997–August 1998)||532,000||137,000||296,000||19,000|
|La Niña year (September 1999–August 2000)||262,000||157,000||43,000||11,000|
|Results are shown for the three most severely smoke-affected regions. These estimates are based on the assumptions used in the principal analysis (see Table 2). aWHO subregions 18–21. bWHO subregion 5 only. cWHO subregions 11–14.|
Monthly resolved emissions estimates were obtained from the Global Fire Emission Database (Global Fire Data 2012; van der Werf et al. 2006), which combines satellite observations of burned area (in square kilometers) with estimates of fuel loads obtained from a biogeochemical model (Giglio et al. 2006). These emissions estimates were used in the GEOS-Chem global 3-D chemical transport model (Bey et al. 2001), which simulates the transport, transformation, and deposition of organic carbon and black carbon aerosols. The model had a 2° (latitude) × 2.5° (longitude) horizontal resolution ~ 222 × 278 km at the equator and 30 vertical layers (Bey et al. 2001). We performed two sets of GEOS-Chem simulations spanning a 10-year period (1997 through 2006). The first included all aerosol emission sources (fossil fuel, biofuel, landscape fires, natural dust, and sea salt), whereas the second excluded landscape fire emissions to separate the contribution from this source.
Finally, we scaled the modeled PM2.5 estimates using two sets of AOD observations from the Moderate Resolution Imaging Spectroradiometer (MODIS) and the Multiangle Imaging Spectroradiometer (MISR) aboard the U.S. National Aeronautics and Space Administration (NASA) Terra satellite (Martonchik et al. 2009; Remer et al. 2005). We maintained the same seasonal, regional, and vertical aerosol distributions as predicted by the GEOS-Chem simulations. Our best estimate of surface PM2.5 (1997 through 2006 average shown in Figure 1) combined information from the model estimates along with the two satellite AOD-scaled estimates:
LFS PM2.5 = [(2 × MODEL) + MODIS + MISR]/4, 
where MODEL is the estimate of PM2.5 from LFS derived from GEOS-Chem and MODIS and MISR are the two satellite AOD-scaled estimates. We multiplied the model contribution by 2 so that our best estimate gave equal weight to the a priori atmospheric model distribution and the sum of the two satellite-scaled estimates. The total aerosol emissions from fires used in the model simulations was 23.5 teragrams (Tg; 1 Tg = 109 kg) per year averaged over 1997 through 2006. Comparable estimates for the MISR and MODIS AOD-based optimizations were 55.0 and 45.5 Tg/year, respectively [see Supplemental Material, Table 1 ( http://dx.doi.org/10.1289/ehp.1104422)] and were within the range of previously published estimates (see Supplemental Material, Table 2). Our best estimate, defined according to Equation 1, was 36.9 Tg/year.
|Source of uncertainty/principal analysis assumption and variations||Annual mortality proportion|
|Estimated PM2.5 concentrations|
|Principal analysis: LFS PM2.5 concentrations estimated from the combination of a global chemical transport model GEOS-Chem and satellite-derived aerosol data from MODIS and MISR||1.00|
|MODEL: PM2.5 concentrations estimated from the GEOS-Chem global chemical transport model||0.68|
|MODIS: MODEL estimate optimized using satellite-derived aerosol data from MODIS||1.47|
|MISR: MODEL estimates optimized using satellite-derived aerosol data from MISR||1.20|
|Pattern of exposure|
|Principal analysis: mortality in sporadically affected subregions estimated using daily average exposure estimates and response functions; mortality in chronically affected WHO subregions estimated using yearly mean exposure estimates and response functions||1.00|
|Sporadic only: mortality in all subregions estimated using daily average exposure estimates and response functions||0.41|
|Chronic only: mortality in all subregions estimated using yearly average exposure estimates and response functions||1.54|
|Shape of concentration–response function|
|Principal analysis: mortality response calculated as a linear function of the PM2.5 concentration||1.00|
|Log-linear: mortality response calculated as a function of the logarithm of the PM2.5 concentration||2.31|
|Counterfactual exposure estimates for chronically affected regions|
|Principal analysis: the counterfactual estimated for each WHO subregion as the mean of the minimum 12-month running-average smoke-specific PM2.5 concentration for each exposure cell within the subregion||1.00|
|Zero: a global value of 0 μg/m3||1.44|
|La Niña: cell-by-cell average for a La Niña year, September 1999–August 2000 inclusive||0.45|
|La Niña regional average: regional average of the values from La Niña||0.81|
|Cell-by-cell minimum: minimum of the 12-month running averages of each cell||0.78|
|Cell-by-cell categorization: global categorization of the values above at the 90th, 97th, and 99th percentiles, applying the average of the category to all cells in the category||0.82|
|Maximum yearly average exposure used for estimating chronic mortality impacts|
|Principal analysis: maximum exposure of 50 μg/m3 was used for estimating the mortality associated with chronic exposure||1.00|
|Maximum exposure of 30 μg/m3 was used for estimating the mortality associated with chronic exposure||0.99|
|Range of minimum and maximum daily exposures used for estimating sporadic exposure impacts|
|Principal analysis: range of exposure assessed was 5–200 μg/m3||1.00|
|Most restrictive range tested: 10–100 μg/m3||0.98|
|Least restrictive range tested: 1–300 μg/m3||1.01|
Evaluation of exposure estimates. Surface measurements of PM2.5 are not available for most regions with high fire emissions. To evaluate the quality of the global exposure estimates, we used ground-based AOD from National Aeronautics and Space Administration’s (NASA) Aerosol Robotic Network (AERONET; NASA 2012) (Holben et al. 1998), PM2.5 measurements from the U.S. Environmental Protection Agency IMPROVE (Interagency Monitoring of Protected Visual Environments) program (Chow and Watson 2002), and visibility data in tropical regions from the National Climatic Data Center Global Summary of the Day (National Oceanic and Atmospheric Administration 2009). Our exposure estimates correlated well with these other measures in regions with high fire activity [see Supplemental Material, Figures 3–6 ( http://dx.doi.org/10.1289/ehp.1104422)]. Correlations (Pearson’s r) of estimated AOD with monthly mean AODs from AERONET were 0.81 in southern Africa (n = 119), 0.90 in northern Africa (n = 74), and 0.76 in Southeast Asia (n = 148; see Supplemental Material, Figure 4). Median correlations between PM2.5 and visibility were 0.57 for sub-Saharan Africa (n = 58), 0.60 for South America (n = 47), and 0.68 for Southeast Asia (n = 13; see Supplemental Material, Figure 6).
Gridded mortality estimates. Country-specific estimates for all-cause all-age mortality in the year 2002 were obtained from the WHO Global Health Observatory (2011). Estimates from the Gridded Population of the World (GPW; version 3) project were used to map country-specific mortality onto the 2° × 2.5° exposure cells (Sociodemographic Data and Applications Centre 2011). The spatial resolution of the GPW data is 2.5 arc-min (~ 4.6 × 4.6 km at the equator), meaning that each exposure cell encompassed 2,880 population cells. To distribute mortality between the population cells, we assigned each cell to the underlying country that contained most of it, summed the GPW population for each country and calculated the percentage of the total population in each cell, and then assigned that percentage of the national mortality to the cell. In the < 1% of cases where population cells were assigned to countries that do not belong to the WHO, we followed the same steps for the 21 WHO subregions and assigned those values instead. The mortality in each exposure cell was estimated by summing the mortality in the 2,880 underlying population cells.
Global burden calculations. Pattern of exposure: subregions of sporadic and chronic impact. Fire activity varied widely across the globe during the 1997 through 2006 period. Some areas were affected sporadically, with a limited number of smoky days in any given year; some areas were affected chronically, with whole seasons being smoke-affected in multiple years. Our principal analysis treats these areas as fundamentally different because acute and chronic PM exposures have independent health effects (Pope and Dockery 2006; Schwartz 2000). We began by classifying each of the 21 WHO subregions as sporadically affected or chronically affected.
The complete set of smoke-specific PM2.5 estimates (12 months × 10 years × 4,208 exposure cells = 504,960) was log-normally distributed with a 90th percentile value of 3 μg/m3. When concentration estimates were rounded to integers, most exposure cells had a value of zero in most months (331,035 of 504,960), indicating low smoke-specific PM2.5. An exposure cell with a 1-month smoke-specific PM2.5 estimate > 3 μg/m3 was classified as being smoke affected during that month. Exposure cells with ≥ 3 smoke-affected months in ≥ 5 of the years were classified as chronically affected (732 of 4,208). Exposure cells that were not chronically affected were classified as sporadically affected (3,476 of 4,208). A WHO subregion was classified as chronically affected if > 50% of its population and/or > 50% of its land area was covered by chronically affected exposure cells (7 of 21; Figure 2). All other WHO regions were classified as sporadically affected (14 of 21; Figure 2).
Burden for sporadically affected subregions. For sporadically affected subregions, we estimated effects of short-term (daily) fluctuations in smoke-specific PM2.5 concentrations on mortality. Daily output from GEOS-Chem was used to estimate the number of days per year that PM2.5 concentrations exceeded a set of threshold values (300, 200, 100, 50, 40, 30, 20, 10, 5, 4, 3, 2, and 1 μg/m3). These threshold values were chosen to provide a range of possible concentrations for sensitivity analyses, and because they reflect clinically relevant increments (10 μg/m3, 100 μg/m3) reported in the literature.
The annual mortality attributable to LFS in each sporadically affected 2° × 2.5° exposure cell was calculated as
Sporadically affected attributable mortality =
where PM is one smoke-specific PM2.5 threshold concentration out of n possible threshold values (see above), DPM is the number of days between PM and the next highest concentration, M is the annual number of deaths in the exposure cell, and RRSI is a relative rate estimate for all-cause mortality due to short-term PM exposure. Although annual mortality is not evenly distributed among the 365 days of the year, there are insufficient data to estimate seasonal mortality on a global scale.
For the principal analysis, a linear RRSI estimate of 0.11% [95% confidence interval (CI): 0, 0.26%] per increase of 1 μg/m3 was used with minimum and maximum concentrations of 5 and 200 μg/m3. This means that cells with daily exposure estimates of < 5 μg/m3 were not included, and cells with exposure estimates > 200 μg/m3 were fixed at a value of 200 μg/m3. The RRSI was calculated using the average (weighted by the inverse of the standard errors) of values from studies reporting associations between all-cause mortality and short-term elevations of ambient PM10 during fire events (Morgan et al. 2010; Sastry 2002) and PM2.5 (Hänninen et al. 2009). Associations with ambient PM10 were converted to associations with PM2.5 by assuming that 75% of all particles < 10 μm were also < 2.5 μm. This is halfway between the 90% ratio measured during fire events (Ward and Hardy 1991) and the 60% ratio used by Cohen et al. (2004) in the initial GBD estimate for urban air pollution.
Burden for chronically affected subregions. No studies have yet reported on the mortality impacts of long-term exposure to LFS. As such, we estimated all-cause mortality in chronically affected exposure cells by assuming the effects of smoke-related PM to be the same as those of urban PM. Specifically, for the principal analysis we assumed a linear 0.64% (95% CI: 0.35%, 0.94%) increase in annual all-cause mortality for each 1-μg/m3 increase in the long-term smoke-specific PM2.5 average, as reported in the American Cancer Society study on urban air pollution (Pope et al. 1995). This is one of the most conservative concentration–response estimates that has been reported in multiple studies of urban PM (Pope and Dockery 2006). The maximum concentration of effect was assumed to be 50 μg/m3. This means that cells with annual exposure estimates > 50 μg/m3 were fixed at a value of 50 μg/m3. The annual mortality attributable to LFS in each chronically affected exposure cell was calculated as
Chronically affected attributable mortality = M × ([RRCI(PM – CF)] – 1), 
where PM is the estimated average annual smoke-specific PM2.5 concentration in the exposure cell based on estimates for 1997 through 2006, CF is the counterfactual concentration for the WHO subregion in which the exposure cell was located, M is the annual number of deaths in the exposure cell, and RRCI is the relative rate of all-cause mortality for long-term PM exposure (i.e., 0.64% for the principal analysis).
The counterfactual concentration is the theoretical minimum annual smoke-specific PM2.5 concentration under ideal conditions. For example, if landscape fires were completely eliminated worldwide, the global counterfactual value would be zero. Given that fire is a natural part of the earth system, we used a more data-driven approach to set counterfactual values for chronically affected WHO subregions. We used a subregion-wide approach because emissions from similar landscapes in neighboring countries can vary widely because of different land management practices, so the theoretical minimum exposure estimated for a single exposure cell might not truly reflect the minimum exposure possible for that particular landscape (Bowman et al. 2011). Specifically, we determined the smallest 12-month running average smoke-specific PM2.5 concentration for each exposure cell within a WHO subregion, and averaged the minimum annual concentrations across all exposure cells to determine the counterfactual value for that WHO subregion.
Sensitivity analyses. There are several sources of uncertainty in our inputs, and we addressed these through multiple sensitivity analyses. First, we assumed both linear and log-linear forms for the concentration–response functions (i.e., RRSI in Equation 2 and RRCI in Equation 3). Although there is increasing evidence of a log-linear association for cardiovascular mortality related to urban air pollution (Pope et al. 2011), we used the linear assumption for the principal analysis because studies on the cardiovascular effects of LFS have been inconclusive. We also tested a range of different exposure limits. For the sporadic assumption, the minimum concentration was varied between 1 and 10 μg/m3 and the maximum was varied between 50 and 300 μg/m3. For the chronic assumption, five alternative counterfactual definitions [a global value of 0 μg/m3; cell-by-cell average for a La Niña year, September 1999–August 2000 inclusive; regional average of the values from La Niña; minimum of the 12-month running averages of each cell; and global categorization of the values above at the 90th, 97th, and 99th percentiles, applying the average of the category to all cells in the category] were tested with maximum yearly average concentrations at 30 and 50 μg/m3. We repeated analyses using the GEOS-Chem and satellite AOD-scaled exposure estimates separately. To assess the effect of our assumptions concerning the combination of sporadic and chronic exposures, all analyses were repeated with all subregions classified as being sporadically affected and with all subregions being classified as chronically affected. There is large interannual variation in emissions from landscape fires mostly driven by changes in climatic conditions (van der Werf et al. 2008). To assess the influence of interannual climatic variability, analyses were repeated with concentration estimates for a strong El Niño year that occurred between September 1997 and August 1998 (inclusive) and a strong La Niña year that occurred between September 1999 and August 2000 (inclusive) (van der Werf et al. 2004).
Exposure. Estimated annual average concentrations ranged from 0 to 45 μg/m3 annually (mean = 1.8 μg/m3; Figure 1). The population-weighted annual average was 2.1 μg/m3, ranging from 0.2 μg/m3 in the Caribbean subregion to 12.2 μg/m3 in sub-Saharan Africa. The population-weighted average number of annual days > 5 μg/m3 was 28, ranging from 6 in the Caribbean subregion to 141 in sub-Saharan Africa.
Burden of mortality. Our principal estimate for the average annual mortality associated with exposure to LFS was 339,000 worldwide, including 157,000 in sub-Saharan Africa and 110,000 in Southeast Asia (Figure 3). The estimates for mortality due to LFS exposure compared with no LFS exposure at all (i.e., a zero exposure counterfactual) were 286,000 in sub-Saharan Africa and 119,000 in Southeast Asia, reflecting much higher background fire activity in sub-Saharan Africa than in Southeast Asia. During the El Niño year, the estimated mortality was higher, particularly in Southeast Asia, where El Niño is associated with dry conditions and more fires (Table 1).
Outputs from all tested models (n = 2,192) had a median of 379,000 and interquartile range of 260,000–600,000 [see Supplemental Material, Figure 7 ( http://dx.doi.org/10.1289/ehp.1104422)]. Results of the sensitivity analyses are shown in Table 2. If a log-linear, rather than linear, concentration–response function was assumed, the mortality estimates more than doubled. The results were also sensitive to the exposure estimates, the assumed pattern of exposure (sporadic vs. chronic), and the choice of the counterfactual exposure estimation, all of which caused the estimated mortality to vary between 0.41 and 1.54 times the principal estimate (Table 2). Results were minimally influenced by the maximum and minimum exposures of effect, which caused the estimates to vary just 0.98 to 1.01 times the principal estimate (Table 2).
Our estimate of 339,000 annual deaths attributable to exposure to LFS is lower than estimates for urban air pollution (800,000) and much lower than estimates for household solid fuel use (1,600,000) (Lopez et al. 2006b). Similar to other environmental risk factors such as unsafe water and indoor and urban air pollution, the mortality burden attributable to LFS falls disproportionately on low-income regions of the world (Figure 4) (Ezzati et al. 2002).
The major strengths of these analyses lie in the use of existing global data sets for terrestrial fire emissions, meteorology, population density, and mortality. Using the WHO geographic subregions and mortality estimates helped make our findings comparable with previously reported estimates for other environmental risk factors. However, there are many limitations inherent in compiling and modeling data at a global scale. A major source of uncertainty comes from the emission factors for fire-derived aerosols that were used to model the exposure estimates. We used emission factors at the lower end of the range in the literature [see Supplemental Material, Table 2 ( http://dx.doi.org/10.1289/ehp.1104422)] even though larger emission factors have been shown to improve model estimates of PM2.5 compared with satellite and surface network observations (Chin et al. 2009; Reid et al. 2009). In addition, the sum of the black carbon and organic carbon emissions factors was often lower than the observed PM2.5 emissions factors, likely resulting in GEOS-Chem underestimates of smoke specific PM2.5. We also chose to be conservative in applying a linear concentration–response function because other studies have suggested higher slopes at lower PM2.5 concentrations (Pope et al. 2009).
In the absence of empirical PM data for many regions most severely affected by LFS, we evaluated our results against global data sets of visibility and ground-based AOD, both of which are proxies for particulate air pollution. Although there was considerable regional variation in the degree of correlation with these independent measures, the estimated PM2.5 performed comparatively well in sub-Saharan Africa and Southeast Asia (the two global regions with highest mortality contributions). Further reductions in uncertainty of the daily exposures could be achieved with the use of higher temporal resolution fire emission inventories. For example, Mu et al. (2011) used active fire observations from Aqua, Terra, and GOES satellites to develop a daily and 3-hourly fire emissions product for the 2002–2010 period.
The WHO subregions with the highest mortality were those we identified as being chronically affected by LFS (Figure 2). The principal estimate of 339,000 annual deaths is composed of 81% mortality due to chronic exposure and 19% due to sporadic exposure. When the analysis was run under the sporadic-only and chronic-only assumptions (Table 2), WHO subregions identified as chronically affected contributed 53% of the total estimates (138,000 and 520,000, respectively) in both cases.
Previous estimates of the global mortality associated with urban air pollution (Cohen et al. 2005) and smoke from household solid fuel use (Lopez et al. 2006b) assumed purely chronic exposure to PM. Our distinction between chronic and sporadic impacts is a departure from this approach, reflecting the current state of epidemiological evidence and the nature of LFS exposure. On the one hand, only a few studies have reported on the mortality effects of LFS (Hänninen et al. 2009; Morgan et al. 2010; Sastry 2002), and all have estimated associations with short-term fluctuations in PM concentrations. On the other hand, urban air pollution studies have clearly demonstrated that chronic exposure to PM is associated with greater increases in mortality than are short-term fluctuations (Pope and Dockery 2006). LFS is episodic in many parts of the world, and annual average exposures are not appropriate for estimating smoke-related mortality in those regions. Similarly, fire smoke exposure is more chronic (because of high seasonal averages) in some regions, and mortality estimates based on short-term fluctuations might be overly conservative. To date, the short-term mortality impacts for PM from landscape fires have been consistent with those of urban PM. Thus, we considered it reasonable to estimate the chronic effects of PM from LFS using conservative values for the chronic effects of PM from urban sources until more specific evidence becomes available. We were also unable to account for different population responses to air pollution. Although our coefficient for acute exposure was driven by a study in Southeast Asia, no studies conducted in sub-Saharan Africa were available.
Estimates of counterfactual exposures are highly uncertain. Human influence on landscape fire activity varies considerably between ecoclimatic regions. We set the theoretical minimum for PM2.5 from LFS as the lowest estimated for each chronically affected WHO subregion over the decade-long study period. However, variation in fire activity during the last decade will not necessarily capture the reduction in fire activity that could be achieved in each environment. For example, tropical rainforests and peat swamps, the primary source of fire emissions in Southeast Asia, rarely burn without human instigation. If such deforestation fires were to be halted, fire activity in this region (and the associated mortality) would be minimal. However, the role of human fire management in savannas, the primary source of emissions in Africa, is less well understood because fire is an integral part of these landscapes (van der Werf et al. 2008). The large estimated influence of El Niño on mortality related to LFS implies that the burden may change in the future if climate change modifies the El Niño Southern Oscillation or drier conditions occur in places with adequate fuels and ignition sources.
Landscape fire activity has been recognized as a global-scale environmental challenge because plumes transgress international boundaries and component gases and particles contribute to climate change (Bowman et al. 2009; Pope and Dockery 2006; van der Werf et al. 2008). This first attempt to quantify the global burden of mortality attributable to LFS has demonstrated important impacts at regional and global scales. We anticipate that subsequent estimates will be improved by better exposure assessment (particularly as empirical PM data become more globally available), further epidemiological studies on mortality and morbidity associated with LFS (particularly in regions with high exposure), and improved understanding of how fire regimes can be modified to reduce smoke emissions. Reducing population level exposure to air pollution from landscape fires is a worthwhile endeavor that is likely to have immediate and measurable health benefits. Such interventions could also potentially provide benefits for the mitigation of climate change and slowing the loss of biodiversity.
This project was supported in part by grants from the Australian Research Council (LP LP0882048), the University of Tasmania (J0018237), and the U.S. National Aeronautics and Space Administration (NNX08AF64G). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
The authors declare they have no actual or potential competing financial interests.
Barregard L, Sällsten G, Gustafson P, Andersson L, Johansson L, Basu S. 2006. Experimental exposure to wood-smoke particles in healthy humans: effects on markers of inflammation, coagulation, and lipid peroxidation.Inhal Toxicol 18(11):845-863 16864402. Crossref, Medline, Google Scholar Bey I, Jacob DJ, Yantosca RM, Logan JA, Field BD, Fiore AM. 2001. Global modeling of tropospheric chemistry with assimilated meteorology: model description and evaluation.J Geophys Res Atmos 106(D19):23073-23095. Crossref, Google Scholar Bowman DMJS, Balch JK, Artaxo P, Bond WJ, Carlson JM, Cochrane MA. 2009. Fire in the Earth system.Science 324(5926):481-484 19390038. Crossref, Medline, Google Scholar Bowman DMJS, Balch J, Artaxo P, Bond WJ, Cochran MA, D’Antonio CM. 2011. The human dimension of fire regimes on earth.J Biogeogr 38(12):2223-2236 22279247. Crossref, Medline, Google Scholar Chin M, Diehl T, Dubovik O, Eck TF, Holben BN, Sinyuk A. 2009. Light absorption by pollution, dust and biomass burning aerosols: a global model study and evaluation with AERONET measurements.Ann Geophys 27:3439-3464. Crossref, Google Scholar Chow JC, Watson JG. 2002. PM2.5 carbonate concentrations at regionally representative Interagency Monitoring of Protected Visual Environment sites.J Geophys Res Atmos 107(D21):8344-8353. Crossref, Google Scholar Cochrane MA, Laurance WF. 2008. Synergisms among fire, land use, and climate change in the Amazon.Ambio 37(7–8):522-527 19205173. Crossref, Medline, Google Scholar Cohen AJ, Anderson HR, Ostro B, Pandey KD, Krzyzanowski M, Künzli N. 2004. Urban air pollution. In Comparative Quantification of Health Risks: Global and Regional Burden of Disease Due to Selected Major Risk Factors, Vol 2 (Ezzati M, Lopez AD, Rodgers A, Murray CJL, eds). Geneva:World Health Organization, 1353–1435. Google Scholar Cohen AJ, Anderson HR, Ostro B, Pandey KD, Krzyzanowski M, Künzli N. 2005. The global burden of disease due to outdoor air pollution.J Toxicol Environ Health A 68(13–14):1301-1307 16024504. Crossref, Medline, Google Scholar Confalonieri U, Menne B, Akhtar R, Ebi KL, Hauengue M, Kovats RS. 2007. Human health. Climate change 2007: impacts, adaptation and vulnerability. In: Contribution of Working Group II to the Fourth Assessment Report of the Intergovernmental Panel on Climate Change (Parry ML, Canziani OF, Palutikof JP, van der Linden PJ, Hanson CE, eds). Cambridge, UKCambridge University Press, 391–431. Google Scholar Danielsen PH, Loft S, Kocbach A, Schwarze PE, Møller P. 2009. Oxidative damage to DNA and repair induced by Norwegian wood smoke particles in human A549 and THP-1 cell lines.Mut Res 674(1–2):116-122 19041418. Crossref, Medline, Google Scholar Delfino RJ, Brummel S, Wu J, Stern H, Ostro B, Lipsett M. 2009. The relationship of respiratory and cardiovascular hospital admissions to the Southern California wildfires of 2003.Occup Environ Med 66(3):189-197 19017694. Crossref, Medline, Google Scholar Denman KL, Brasseur G, Chidthaisong A, Ciais P, Cox PM, Dickinson RE. 2007. Couplings between changes in the climate system and biogeochemistry. In Climate Change 2007: The Physical Science Basis. Contribution of Working Group I to the Fourth Assessment Report of the Intergovernmental Panel on Climate Change (Solomon S, Qin D, Manning M, Chen Z, Marquis M, Averyt KB, et al., eds). Cambridge:Cambridge University Press, 499–587. Google Scholar Ezzati M, Lopez AD, Rodgers A, Vander Hoorn S, Murray CJ. 2002. Selected major risk factors and global and regional burden of disease.Lancet 360(9343):1347-1360 12423980. Crossref, Medline, Google Scholar Fargione J, Hill J, Tilman D, Polasky S, Hawthorne P.. 2008. Land clearing and the biofuel carbon debt.Science 319(5867):1235-1238 18258862. Crossref, Medline, Google Scholar Giglio L, van der Werf GR, Randerson JT, Collatz GJ, Kasibhatla P. 2006. Global estimation of burned area using MODIS active fire observations.Atmos Chem Phys 6:957-974. Crossref, Google Scholar
- Global Fire Data. 2012. Global Fire Data (GFED). Version 2.Available: http://www.globalfiredata.org/index.html [accessed 30 March 2012]. Google Scholar
Hänninen OO, Salonen RO, Koistinen K, Lanki T, Barregard L, Jantunen M. 2009. Population exposure to fine particles and estimated excess mortality in Finland from an East European wildfire episode.J Expo Sci Environ Epidemiol 19(4):414-422 18523459. Crossref, Medline, Google Scholar Holben BN, Eck TF, Slutsker I, Tanre D, Buis JP, Setzer A. 1998. AERONET—a federated instrument network and data archive for aerosol characterization.Remote Sens Environ 66:1-16. Crossref, Google Scholar
- Institute for Health Metrics and Evaluation. 2010. Global Burden of Diseases Study.Available: www.globalburden.org/ [accessed 2 August 2011]. Google Scholar
Jayachandran S.. 2009. Air quality and early-life mortality during Indonesia’s massive wildfires in 1997.J Human Resour 44(4):916-954. Google Scholar Johnston FH, Bailie R, Pilotto L, Hanigan I. 2007. Ambient biomass smoke and cardio-respiratory hospital admissions.BMC Public Health 7(1):240; doi: 10.1186/1471-2458-7-240[Online 13 September 2007] 17854481. Crossref, Medline, Google Scholar Johnston FH, Hanigan IC, Henderson SB, Morgan GG, Bowman DMJS. 2011. Extreme air pollution events from bushfires and dust storms and their association with mortality in Sydney, Australia 1994–2007.Environ Res 111(6):811-816 21601845. Crossref, Medline, Google Scholar Kocbach A, Namork E, Schwarze PE. 2008. Pro-inflammatory potential of wood smoke and traffic-derived particles in a monocytic cell line.Toxicology 247(2–3):123-132 18406506. Crossref, Medline, Google Scholar Langmann B, Duncan B, Textor C, Trentmann J, van der Werf GR. 2009. Vegetation fire emissions and their impact on air pollution and climate.Atmos Environ 43(1):107-116. Crossref, Google Scholar Lohman DJ, Bickford D, Sodhi NS. 2007. The burning issue.Science 316(5823):376 17446376. Crossref, Medline, Google Scholar Lopez AD, Mathers CD, Ezzati M, Jamison DT, Murray CJL. 2006a. Global and regional burden of disease and risk factors, 2001: systematic analysis of population health data.Lancet 367(9524):1747-1757 16731270. Crossref, Medline, Google Scholar
- Lopez AD, Mathers CD, Ezzati M, Jamison DT, Murray CJL, Abramson M, eds. 2006b. Global Burden of Disease and Risk Factors. New York:Oxford University Press and World Bank Publications. Google Scholar
Martonchik JV, Kahn RA, Diner DJ. 2009. Retrieval of aerosol properties over land using MISR observations. In Satellite Aerosol Remote Sensing over Land (Kokhanovsky AA, de Leeuw G, eds). Berlin:Springer, 267–293. Google Scholar Morgan G, Sheppeard V, Khalaj B, Ayyar A, Lincoln D, Jalaludin B. 2010. The effects of bushfire smoke on daily mortality and hospital admissions in Sydney, Australia, 1994 to 2002.Epidemiology 21(1):47-55 19907335. Crossref, Medline, Google Scholar Mott JA, Mannino DM, Alverson CJ, Kiyu A, Hashim J, Lee T. 2005. Cardiorespiratory hospitalizations associated with smoke exposure during the 1997, Southeast Asian forest fires.Int J Hyg Environ Health 208(1–2):75-85 15881981. Crossref, Medline, Google Scholar Mu M, Randerson JT, van der Werf GR, Giglio L, Kasibhatla P, Morton DC. 2011. Daily and hourly variability in global fire emissions and consequences for atmospheric model predictions of carbon monoxide.J Geophys Res Atmos 116(D24303):1-19. Google Scholar Naeher LP, Brauer M, Lipsett M, Zelikoff JT, Simpson CD, Koenig JQ. 2007. Woodsmoke health effects: a review.Inhal Toxicol 19(1):67-106 17127644. Crossref, Medline, Google Scholar
- NASA (National Aeronautics and Space Administration)Global Change Master Directory. The AErosol RObotic NETwork (AERONET) from NASA/GSFC. Available: http://gcmd.nasa.gov/records/GCMD_AERONET_NASA.html [accessed 30 March 2012]. Google Scholar
- National Oceanic and Atmospheric Administration. 2009. National Climatic Data Center.Available: http://www.ncdc.noaa.gov/oa/ncdc.html [accessed 17 May 2011]. Google Scholar
Pope CA, Burnett RT, Krewski D, Jerrett M, Shi Y, Calle EE. 2009. Cardiovascular mortality and exposure to airborne fine particulate matter and cigarette smoke: shape of the exposure–response relationship.Circulation 120(11):941-948 19720932. Crossref, Medline, Google Scholar Pope CA, Burnett RT, Turner MC, Cohen AJ, Krewski D, Jerrett M. 2011. Lung cancer and cardiovascular disease mortality associated with ambient air pollution and cigarette smoke: shape of the exposure–response relationships.Environ Health Perspect 119:1616-1621 21768054. Link, Google Scholar Pope CA, Dockery DW. 2006. Health effects of fine particulate air pollution: lines that connect.J Air Waste Manag Assoc 56(6):709-742 16805397. Crossref, Medline, Google Scholar Pope CA, Thun MJ, Namboodiri MM, Dockery DW, Evans JS, Speizer FE. 1995. Particulate air pollution as a predictor of mortality in a prospective study of U.S. adults.Am J Respir Crit Care Med 151(3 Pt 1):669-674 7881654. Crossref, Medline, Google Scholar Reid JS, Hyer EJ, Prins EM, Westphal DL, Zhang J, Wang J. 2009. Global monitoring and forecasting of biomass-burning smoke: Description of and lessons from the Fire Locating and Modeling of Burning Emissions (FLAMBE) program.IEEE J Sel Top Appl Earth Obs Remote Sens 2(3):144-162. Crossref, Google Scholar Reid JS, Koppmann R, Eck TF, Eleuterio DP. 2005. A review of biomass burning emissions part II: intensive physical properties of biomass burning particles.Atmos Chem Phys 5:799-825. Crossref, Google Scholar Remer LA, Kaufman YJ, Tanré D, Mattoo S, Chu DA, Martins JV. 2005. The MODIS aerosol algorithm, products, and validation.J Atmos Sci 62:947-973. Crossref, Google Scholar Sanhueza PA, Torreblanca MA, Diaz-Robles LA, Schiappacasse LN, Silva MP, Astete TD. 2009. Particulate air pollution and health effects for cardiovascular and respiratory causes in Temuco, Chile: a wood-smoke-polluted urban area.J Air Waste Manag Assoc 59(12):1481-1488 20066914. Crossref, Medline, Google Scholar Sastry N.. 2002. Forest fires, air pollution, and mortality in southeast Asia.Demography 39(1):1-23 11852832. Crossref, Medline, Google Scholar Schwartz J.. 2000. Harvesting and long term exposure effects in the relation between air pollution and mortality.Am J Epidemiol 151(5):440-448 10707911. Crossref, Medline, Google Scholar Schweithelm J, Glover D, Jessup T. 2006. Conclusions and policy recommendations. In Indonesia’s Fire and Haze: The Cost of Catastrophe (Glover D, Jessup T, eds). Ottawa, Ontario, Canada, and Singapore:International Development Research Centre and the Institute of Southeast Asian Studies, 130–143. Google Scholar
- Sociodemographic Data and Applications Centre. 2011. Gridded Population of the World and the Global Rural-Urban Mapping Project. Colombia University.Available: http://sedac.ciesin.columbia.edu/gpw/ [accessed 17 May 2011]. Google Scholar
Tosca MG, Randerson JT, Zender CS, Flanner MG, Rasch PJ. 2010. Do biomass burning aerosols intensify drought in equatorial Asia during El Niño?Atmos Chem Phys 10:3515-3528. Crossref, Google Scholar Turetsky MR, Kane ES, Harden JW, Ottmar RD, Manies KL, Hoy E. 2011. Recent acceleration of biomass burning and carbon losses in Alaskan forests and peatlands.Nat Geosci 4:27-31. Crossref, Google Scholar van der Werf GR, Randerson JT, Collatz GJ, Giglio L, Kasibhatla PS, Arellano AF. 2004. Continental-scale partitioning of fire emissions during the 1997 to 2001 El Niño/La Niña period.Science 303(5654):73-76 14704424. Crossref, Medline, Google Scholar van der Werf GR, Randerson JT, Giglio L, Collatz GJ, Kasibhatla PS, Arellano AF. 2006. Interannual variability in global biomass burning emissions from 1997 to 2004.Atmos Chem Phys 6:3423-3441. Crossref, Google Scholar van der Werf GR, Randerson JT, Giglio L, Collatz GJ, Mu M, Kasibhatla PS. 2010. Global fire emissions and the contribution of deforestation, savanna, forest, agricultural, and peat fires (1997–2009).Atmos Chem Phys 10:11707-11735. Crossref, Google Scholar van der Werf G, Randerson J, Giglio L, Gobron N, Dolman A.. 2008. Climate controls on the variability of fires in the tropics and subtropics.Global Biogeochem Cycles 22(3):GB3028; doi: 10.1029/2007GB003122[Online 5 September 2008]. Crossref, Google Scholar Ward DE, Hardy CC. 1991. Smoke emissions from wildland fires.Environ Int 17(2–3):117-134. Crossref, Google Scholar Westerling AL, Hidalgo HG, Cayan DR, Swetnam TW. 2006. Warming and earlier spring increase western U.S. forest wildfire activity.Science 313(5789):940-943 16825536. Crossref, Medline, Google Scholar
- WHO (World Health Organization). Global Health Observatory Data Respository. 2011. Global Health Observatory Data Repository.Available: http://apps.who.int/ghodata/ [accessed 30 March 2012]. Google Scholar
Yokelson RJ, Karl T, Artaxo P, Blake DR, Christian TJ, Griffith DWT. 2007. The Tropical Forest and Fire Emissions Experiment: overview and airborne fire emission factor measurements.Atmos Chem Phys 7:5175-5196. Crossref, Google Scholar