The health impacts of Indonesian peatland fires
Environmental Health volume 21, Article number: 62 (2022)
Indonesian peatlands have been drained for agricultural development for several decades. This development has made a major contribution to economic development. At the same time, peatland drainage is causing significant air pollution resulting from peatland fires. Peatland fires occur every year, even though their extent is much larger in dry (El Niño) years. We examine the health effects of long-term exposure to fine particles (PM2.5) from all types of peatland fires (including the burning of above and below ground biomass) in Sumatra and Kalimantan, where most peatland fires in Indonesia take place.
We derive PM2.5 concentrations from satellite imagery calibrated and validated with Indonesian Government data on air pollution, and link increases in these concentrations to peatland fires, as observed in satellite imagery. Subsequently, we apply available epidemiological studies to relate PM2.5 exposure to a range of health outcomes. The model utilizes the age distribution and disease prevalence of the impacted population.
We find that PM2.5 air pollution from peatland fires, causes, on average, around 33,100 adults and 2900 infants to die prematurely each year from air pollution. In addition, peatland fires cause on average around 4390 additional hospitalizations related to respiratory diseases, 635,000 severe cases of asthma in children, and 8.9 million lost workdays. The majority of these impacts occur in Sumatra because of its much higher population density compared to Kalimantan. A main source of uncertainty is in the Concentration Response Functions (CRFs) that we use, with different CRFs leading to annual premature adult mortality ranging from 19,900 to 64,800 deaths. Currently, the population of both regions is relatively young. With aging of the population over time, vulnerabilities to air pollution and health effects from peatland fires will increase.
Peatland fire health impacts provide a further argument to combat fires in peatlands, and gradually transition to peatland management models that do not require drainage and are therefore not prone to fire risks.
Indonesia faces peat and forest fires every year, with the most severe fires occurring in El Niño years, including 2014, 2015 and 2019 . During these events, smoke extends over large parts of Indonesia as well as Singapore and parts of Malaysia . Much of the smoke is from peatland fires. Out of some 15 million hectares of peat in Indonesia, over half have been cleared and drained, in particular for plantation development (including oil palm and acacia for pulp and paper production) . Once drained, peatlands are susceptible to fires, which can have a significant impact on ambient air quality. Drained, but still wet peat soils burn incompletely, at relatively low temperatures, which results in relatively high emissions of a mix of pollutants including particulate matter, carbon monoxide and Polycyclic Aromatic Compounds (PACs) . The PM2.5 (particulate matter with a diameter less than 2.5 μm) concentrations reached in Indonesia during fire events are exceptionally high. For example, in 2015 in Central Kalimantan province, 24-h mean concentrations above 2000 μg/m3 were recorded [4, 5], far above short–term exposure levels considered hazardous for human health by the World Health Organization, i.e., 15 μg/m3 (24 hour mean) .
In the epimidiological literature, it has become common practice to use PM2.5 as an indicator of ambient air quality, and causal associations between PM2.5 exposure and health effects have been derived. The health effects of particulate air pollution include respiratory and cardiovascular morbidity (e.g., aggravation of asthma, respiratory symptoms and an increase in hospital admissions) and mortality from cardiovascular and respiratory diseases as well as contribute to lung cancer . Given that forest and peatland fires occur across Indonesia each year, and that smoke is retained in the ambient air for several months a year, forest and peatland fires not only lead to spikes in air pollution over the short-term, but also increase the population long-term exposure to PM2.5, which is known to have a greater impact on mortality than short-term episodes [8, 9]. Several studies have analyzed health effects of Indonesian peat fires. For example, the short-term health effects of the 1997 fires in Sumatra and Kalimantan were studied in  and the long-term health effects in Central Kalimantan province in . Air quality and health impacts of vegetation and peat fires in Equatorial Asia during 2004–2015 were studied in a comprehensive modelling study in . However none of these studies analysed average long-term health effects including both mortality and morbidity effects, based on a comprehensive spatial database of PM2.5 concentrations.
The objective of this paper is to analyse the long-term health impact of Indonesian peatland fires over the 5 year period extending between 2013 and 2017. We focus on the lowland provinces of Sumatra and Kalimantan which contain extensive peat areas and where people have the greatest exposures to the fires. We determine (i) the PM2.5 ambient air pollution resulting from fires occurring in peat; (ii) exposure of people to PM2.5; and (iii) the resulting heath impacts from long-term exposure including adult premature mortality, lost work days, cases of childhood asthma and chronic bronchitis, respiratory hospital admissions and infant mortality. Our paper focusses on Sumatra and Kalimantan where the majority of peatland fires take place and where the health effects of these fires are the highest. We include all types of fire in peatlands including fires affecting above ground vegetation and/or fires that occur in the form of burning and smouldering of below ground peat biomass. Our methodology is novel in the sense that it focusses on peatland fires, and integrates satellite data on air quality, calibrated with data from local monitoring stations, with state-of-the-art epidemiological modelling resulting in a high resolution (4.4 km) model for a large area (~ 1 million km2). We also assess the implications of using different concentration-response functions. We acknowledge that health effects of air pollution from peatland fires in Sumatra and Kalimantan occur in a broader area including Singapore and Malaysia (e.g. ), but we focus our study on Indonesia, where the fires originate.
Case study area
We focus on the health effects of peatland fires in Sumatra and Kalimantan, and while peat and forest fires also occur in Papua and West Papua provinces, their much lower population will result in smaller impacts from the fires. Sumatra covers 48 million hectares (ha). With around 21% of the Indonesian population, it is the most populous island after Java. Sumatra has in total around 6.4 million ha of peat according to government data , most of it located along the east coast of the island. Kalimantan is the Indonesian part of the island of Borneo and covers around 54 million ha. Based on government data of the Ministry of Agriculture , Kalimantan has approximately 4.9 million ha of peatland. Peat occurs in all provinces of Kalimantan: South, Central, West, East and North Kalimantan. Large scale deforestation and peat conversion started in the 1980s, initially in Sumatra. In recent years it has slowed down in response to government policies, but peat conversion is still ongoing in particular for the cultivation of oil palm. The peat extent, and land cover in the peatlands of Sumatra and Kalimantan are shown in Fig. 1. Peatland fires typically occur in the dry season – starting in June and extending to December in some parts of Indonesia – with most fires taking place in the period from August to October.
Peatland fire occurrence
The burned areas in Sumatra and Kalimantan in each year from 2013 to 2017 were derived from satellite imagery collected through MODIS (MCD64A1 Collection 6). The MODIS MCD64A1 product uses a burned area mapping algorithm to detect changes in vegetation using the optical MODIS sensors , providing daily data with a 500 m grid resolution. Daily MODIS MCD64A1 data were downloaded as shapefiles, aggregated by year, and overlaid with the Government of Indonesia Ministry of Agriculture peat map from 2011 , in order to identify which fires took place in areas that have peat soils (as opposed to mineral soils). The number of fires varied considerably between years, with over 100,000 individual fires recorded during the particularly dry year 2015.
In order to specify the type of land cover where the burned areas were identified, the results were subsequently overlaid with the land cover maps of the Ministry of Environment and Forestry from 2013 to 2017. These land cover maps show the spatial extent of 22 classes of land cover, which are aggregated here into five broad categories: 1) forests (intact and degraded peat swamp forest); 2) plantation (including oil palm, hevea rubber and coconut) and forest plantation; 3) degraded lands (bare ground, wet/dry shrubs, savannas, and grassy areas); 4) agricultural lands; and 5) other categories. For each of these classes, emission factors were calculated, in order to estimate which proportion of ambient PM2.5 is derived from fires in peatland (discussed below). All spatial analyses were conducted in the ArcMap 10.5 GIS package.
Ambient PM2.5 from peatland fires
Estimating ambient PM2.5 concentrations from satellites and simulation
Due to the limited availability of ground measurements of PM2.5 in Indonesia, we use estimates of surface PM2.5 concentrations derived from a combination of satellite remote sensing and chemical transport model simulation . Monthly estimates of surface PM2.5 concentrations are derived over Sumatra and Kalimantan for 2013 to 2017 at a spatial grid resolution of 4.4 km. PM2.5 concentrations are estimated using satellite-derived Aerosol Optical Depth (AOD). AOD is a measure of the extinction (through scattering and absorption) of light by aerosols in a column of air from the Earth’s surface up to space. The AOD is used to represent the abundance of total aerosol particles in an atmospheric column . The daily surface PM2.5 concentrations from each data source are obtained by applying the daily simulated AOD to PM2.5 ratios to the coincident daily calibrated AOD sources (see Supplementary Materials). Monthly means are calculated from the daily PM2.5 values. The monthly mean PM2.5 concentrations from each source are then combined based on their relative uncertainties with AERONET following . Details are presented in the Supplementary Materials. Modelled PM2.5 concentrations, from AOD, are calibrated with PM2.5 and PM10 data from Indonesia’s air quality monitoring network, that comprises 16 stations in Sumatra and 13 stations in Kalimantan, providing daily data from 2015 to 2017 (operated by the Meteorological, Climatological and Geophysical Agency and the Ministry of Environment and Forestry) [19, 20]. The majority (90%) of PM emitted from peatland fires is in the fine particulate fraction (PM2.5) .
Identifying the contribution of peatland fires to ambient PM2.5
The MODIS data show that fires almost exclusively take place in the period June to December. The Indonesian dry season varies somewhat between the southern part of Kalimantan and the northern part of Sumatra, but generally lasts from May to October. Hence, the peat burns in the dry season as well as in the beginning of the wet season, especially in El Nino years such as 2015 when major parts of Sumatra and Kalimantan were unusually dry from September to November . Consequently, in Sumatra and Kalimantan, the months with high PM2.5 concentration in the peat areas (above 50 μg/m3) are all during the period June to December, with the highest exposures occurring during the months from August to October (August being the peak of the dry season). This is contrary to Java, which is much more industrialized and more densely populated, and where concentrations exceeding 50 μg/m3 occur throughout the year. Note that, due to prevailing wind directions during the dry season, the peatland fire smoke generally does not reach Java . Hence, we assume that the increment in the PM2.5 concentration during the dry season, that we observe throughout Kalimantan and Sumatra, can be attributed to fires. To make sure all dry season peatland fires are included in the analysis, based on the rainfall distribution over the years of the analysis, the peat fire season is defined, for the purpose of this analysis, as the period between June and December. The period January to May is defined as the wet season, where we assume there are no peat fires. For each province, the mean concentration during the peat fire season and during the wet season is calculated. The difference between these two is assumed to represent the contribution of biomass fires to ambient PM2.5. This difference is then multiplied by 7/12 to calculate the annual average increase in PM2.5 concentration, which we use to quantify the health impact of peatland fires. This is because we relate health effects to the mean annual increase in PM2.5 as a consequence of peatland fires, in line with the epidemiological evidence that the health effects of long-term exposure are much more important than the impacts from peak exposure [8, 9]. The outcome of the analysis is not sensitive to the assumed length of the peat fire season. If it is assumed that the peat fire season is 6 months, for example, a higher monthly average increase during this season would be noticed, but this would be multiplied with 6/12, resulting in approximately the same annual increase in ambient PM2.5 due to peatland fires.
Next, we analyze which part of the incremental PM2.5 can be attributed to fires in peatlands, as opposed to PM2.5 originating from fires occurring in ecosystems on mineral soil (for instance from forest fires or the burning of biomass in crop fields). For each island, we assume that the increase in the concentration of PM2.5 from peatland fires is proportional to the share of PM2.5 emitted in peatland fires as opposed to PM2.5 emitted from fires on mineral soils. We estimate this proportion based on the number of hectares burned in peat versus mineral soil, taking into consideration differences between years, type of ecosystem and emission factors (EF). We acknowledge that the formation of secondary organic aerosols (SOA) may differ between peatland and forest fires, but we have not taken this difference into consideration in the analysis. We comment on the consequences of our simplification in the Discussion section. The burned areas are analyzed with MODIS images. The location of the fires is compared to the land cover map of MOEF , and the national peat map . The amount of biomass burned (above ground biomass and, in peat soil, soil biomass) per hectare was analysed for various ecosystem types (forests, plantations, shrubland, grassland and cropland). It was assumed that on average peat soils burn to a depth of 33 cm, and that per hectare of burned peat 505 t of underground biomass is lost [23,24,25]. The resulting EF for peat is 4.6 t PM2.5 per hectare burned . For above ground biomass, the EF ranges from 0.03 t PM2.5 per hectare for fires in agricultural land to 3.5 t PM2.5 per hectare for fires in forests [27, 28]). It is assumed that these numbers apply to both ecosystems on peat and on mineral soil. For peatland we add the PM2.5 emissions from peat burning (below ground) to the emissions from above ground biomass burning to calcuate the overall EF. Supplementary Materials (SM) Table A2 in the Annex provides detailed calculations.
Finally, for each province, month and year the population weighted increase in PM2.5 concentration due to peatland fires was calculated. This is important since the populations are not homogeneously distributed across provinces, and, consequently, the spatially average concentrations do not correlate perfectly with the true population exposure. We used population data from the NASA Socioeconomic Data and Applications Center (SEDAC), file: Gridded Population of the World 2015.Footnote 1 This dataset contains population data at 1 km grid cell resolution, world-wide. In our analysis, for every grid cell, the increase in PM2.5 concentration in the grid cell is multiplied with the population in the grid cell. Next, these numbers are added for all grid cells in a province, and the resulting number is divided by the number of people in the province. This results in the population-weighted average increase in PM2.5 concentration due to peatland fires by province.
Health effects from peatland fires
We follow a health impact assessment approach to calculate the health effects of peatland fires from long-term exposure (1-year or longer) to PM2.5 ambient air concentrations. We did not specifically calculate the health impacts attributable to episodic events of air pollution, although these short-term effects are implicitly included in the epidemiological assessment of the average long-term health consequences. When analyzing the impact of PM2.5 on human health, four main inputs are used to quantify effects on people: (1) current and counterfactual PM2.5 concentrations to determine the change in PM2.5; (2) size and age-composition of population groups exposed to current levels of air pollution; (3) baseline incidence of mortality and morbidity outcomes; and (4) risk functions relating a change in concentration to a change in the health outcome of interest. We analyze the annual average health effects of peatland fires, relating annual average health effects to annual average PM2.5 concentrations (in line with, e.g., [8, 9, 29, 30]. In reality, variation between years may occur since peak concentrations in PM2.5 that exacerbate health effects in the short term occur more frequently in dry years with more peatland fires. To capture the variability of the impact results, we also conduct a sensitivity analysis to the key input factors in the analysis, as described below. The health impact is calculated as follows:
Current and counterfactual PM2.5 concentrations
The health impacts of air pollution related to peatland fires are compared to the health effects of air pollution that would occur without peatland fires (in which case there would still be a certain level of air pollution, e.g. from industry or traffic). This background air pollution concentration is referred to as the “counterfactual”, the ambient air concentration that would be expected in the absence of peatland fires, and the dry season excess exposure above the counterfactual is used to evaluate the excess health impact attributed to peatland fires. In the assessment of the long-term adult mortality, the incremental PM2.5 from peatland fires is defined as the 5-year average of the population-weighted, province-specific concentrations from June to December (fire season) minus the concentrations from January to May (the counterfactual). On the other hand, infant mortality and morbidity impacts are calculated based on the mean PM2.5 incremental exposure above the counterfactual in a given year.
Size and age composition of population groups exposed to current levels of PM2.5
For each province in Sumatra and Kalimantan, annual population statistics, including size and age-group stratification, were calculated according to the Indonesian Statistics Agency and Bappenas population scenario . In 2013, the population of Sumatra and Kalimantan was 53.5 and 15.0 million persons, with an average annual population growth rate between 2013 and 2017 equal to 1.2 and 1.5% for each region, respectively. The average population composition over the 5-year period by broad age-group including infants, children (ages 1 to 14 years), young adults (15 to 30 years) and adults older than 30 years was 2.2, 27.7, 26.2, and 43.9%. Age structure was similar for both islands.
Baseline rates of mortality
Next, we estimate the underlying rate of the health effect (e.g., the baseline natural mortality rate, excluding accidental deaths, in the population in terms of deaths per thousand people). We found two sources for these data, so a sensitivity analysis was conducted. For the base case, the data used were taken from the Global Health Estimates (GHE) database by WHO , which provides all-cause and cause-specific mortality estimates based on official Indonesian national statistics as of year 2016 (the last year available for such data). For year 2016, the GHE indicates 12.56 natural deaths per 1000 adults aged 30 and older, and 9.53 deaths for non-communicable diseases and lower respiratory illness combined (NCD + LRI) per 1000 adults 25 and older in Indonesia. These baseline mortality rates have been applied, without further adjustment, to population estimates from the Bappenas population scenario for each of the 5 years between 2013 and 2017 to calculate the annual baseline mortality . For infant (< 1 year) mortality, modeled data were obtained from the United Nations Inter-agency Group for Child Mortality Estimation . The annual infant mortality rate in Indonesia over the period 2013 to 2017 was: 25.25, 24.34, 23.46, 22.62, and 21.86 deaths per 1000 live births. These data were used in the calculation of the annual infant mortality.
For the sensitivity calculations, we used modeled demographic data from the Global Health Data exchange (GHDx) database of the Institute for Health Metrics and Evaluation.Footnote 2 At the national-level, IHME estimates for year 2016 are 11.09 natural deaths per 1000 adults over 30 years, and 8.42 NCD + LRI deaths per 1000 adults aged 25 and over. Infant mortality figures are fairly close to the values by the UN IGME . The IHME database also provides modeled mortality rates by province (e.g., Aceh, West Kalimantan, etc.). For Sumatra and Kalimantan islands, respectively, the mortality rates are 5.50 and 6.41 natural deaths per 1000 population (compared to 6.71 per 1000 at the national-level from WHO’s GHE database), 10.05 and 11.40 natural deaths per 1000 adults over 30 years, and 7.56 and 8.49 NCD + LRI deaths per 1000 adults 25 and over.
The fourth input is based on statistical relationships from the epidemiological literature that relate ambient concentrations of PM2.5 to mortality and other health effects. Over the past decades, epidemiological studies on biomass-derived PM have investigated the mortality and morbidity effects of wildfires on respiratory and cardiovascular outcomes [34,35,36];;. While the toxicological evidence is suggestive that wildfire PM is more toxic than urban PM, we have adopted the conservative assumption that all particles are equitoxic, and have applied in our analysis the most robust CRFs recently derived on the basis of all PM sources. The accuracy of the risk functions depends in part on: (a) the data quality and methodology of the original studies; (b) the extent to which the concentration-response functions (CRF) include the range of concentrations for which they are being applied, and (c) the applicability of these functions to populations besides those from which they were originally estimated.
For the base-case health assessment of adult mortality from exposure to ambient PM2.5 resulting from peatland fires, a meta-analysis was used . This study included 53 cohort studies, 39 studies from North America, eight from Europe, and six from Asia. The mean concentration of PM2.5 across the studies was 15.7 μg/m3, with higher concentrations observed in the Asian studies with a mean of 30.5 μg/m3. The authors’ analysis indicated robust association between PM2.5 and premature mortality. At the mean concentration of 15.7 μg/m3, a 10 μg/m3 increase in PM2.5 was associated with a 10.3% (95% confidence interval (CI) = 9.7, 11.1%) increase in adult premature mortality. This risk function was used for the base case estimates since it included the greatest number of studies in total, with a few from Asia. The referenced health risk was applied to the mortality rate in the population over 30 years old .
A sensitivity analysis was performed by using three alternative concentration-response functions. The first is an estimate of adult mortality; a meta-analysis of 41 cohort studies from around the world that was conducted while allowing for flexibility in determining the shape of the CRF . The second study , utilized census data on 2.4 million Canadian adults to examine non-accidental and cause-specific mortality between 2001 and 2011. Of particular interest, the study compared impacts of using three different average exposure periods (1, 3, and 8 years), so it is relevant to the five-year average of concentrations due to peatland fires. The health risk from  was applied to the mortality rate in the population over 30 years old. The third study used in the sensitivity analysis is based on the recommendations from WHO’s HRAPIE revie w, which was based primarily on an earlier meta-analysis . The estimated risk for non-accidental mortality (for the population 30+) was based on eleven studies available at the time. The resulting risk was 6.2% (95% CI = 4, 8.2%) per 10 μg/m3 increment in PM2.5 concentration. As shown in Annex 2, the CRF that we use in our base case , lies in between those of the three other studies.
Infant mortality (age < 1 year) was assessed based on . This study considered almost 1 million births across sub-Saharan Africa using the Demographic and Health Surveys, a set of nationally representative samples from 30 sub-Saharan countries. PM2.5 was estimated at the residential location using remote sensing satellite data and assigned to each birth using exposures 9 months before and 12 months after birth. The study reported that a 10 μg/m3 change in PM2.5 was associated with 9.2% (95% CI = 4, 14%) increase in infant mortality. Previous studies in Mexico City [39, 40] generated risk estimates of 8.8 and 6.9%, respectively, for a 10 μg/m3 change in PM2.5. Finally, another large meta-analysis , included over a half million births from 69 nationally representative Demographic and Health Surveys conducted in 43 countries throughout Africa, the Middle East, Asia, Eastern Europe and the Caribbean. The study suggests an increase in infant mortality of 8.1% (95% CI = 2.5, 14.3%) for a 10 μg/m3 change in fire-related PM2.5, which is quite similar to the result of .
For morbidity, based on available data, the following were quantified: severe asthma attacks in children, asthma-related hospital admissions (all ages), all respiratory hospital admissions (all ages), and lost work days for the current labor force. Thus, outcomes such as cardiovascular morbidity (heart attacks, hospital admissions), respiratory symptoms leading to restricted activity and adverse birth outcomes (besides infant mortality) were not quantified because of uncertainties regarding baseline levels in Sumatra and Kalimantan. Most of the CRFs were based on WHO’s HRAPIE , except one study , which was used for asthma-related hospital admissions. Baseline prevalence rates were determined as follows: severe child asthma attacks were assumed to be 1.24 cases/year per child, assuming 62 asthma attacks a year, of which 2% are severe ([7, 43];, respiratory hospital admissions were assumed as 397 per 100,000 people ,and lost work days were 7 days/worker (UN ILO). For asthma-related hospital admissions, we were unable to find baseline evidence for Indonesia, instead, we used evidence from Europe, and assumed that 11% of all respiratory hospital admissions are asthma-related.
Peatland fire occurrence
Between 2013 and 2017, peatland fires occurred every year in Sumatra and Kalimantan (Table 1), with the largest burned areas detected in 2014 (0.72 million ha) and 2015 (0.42 million ha). The three provinces with the largest burned areas were: Riau, Central Kalimantan and South Sumatra. On average, over the period 2013–2017, 2.3% of all peatlands in Sumatra and 3.0% of all peatlands in Kalimantan burned each year. However, there were large variations between years. For instance, in 2014, 5.4% of all peatlands burned in Sumatra, and 7.8% in Kalimantan. Within the peatlands, most fires occurred in degraded lands (bare ground, wet and dry shrub, savanna and grassland), totaling 52 to 80% of total burned areas annually, followed by plantations (6 to 37%), agricultural lands (3 to 13%) and forests (including primary and degraded forest) (4 to 10%; see Annex 1, Table A1, for details). The largest fires in degraded lands occurred in 2014 and 2015 when 0.46 million ha and 0.34 million ha of degraded peatlands were burned, respectively.
Ambient PM2.5 from peatland fires
The monthly satellite-derived (SAT) surface PM2.5 concentration estimates for Indonesia are compared against monthly ground monitor data (Insitu) in Table 2. The comparison is shown for ‘All Sites’ (direct PM2.5 sites and the sites where PM2.5 is calculated from PM10), and for direct PM2.5 sites (Direct PM2.5 Sites). The ground measurements only cover the period of 2014 to 2017. The majority of the ground-based observations over Indonesia provided PM10 concentrations rather than PM2.5. The PM10 measurements were used to estimate PM2.5, based on sites where there were coincident PM2.5 and PM10 values. The column ‘All Sites’ in Table 2 contains the analysis for all sites (the sites measuring directly PM2.5 values as well as the sites where PM2.5 was estimated from PM10), while the column ‘Direct PM2.5 Sites’ is referring to the analysis for only sites which directly measured PM2.5.
There is a high degree of consistency (R2 = 0.77 for all sites; R2 = 0.90 for direct PM2.5 sites) between the PM2.5 estimates and the independent ground monitor data, with slopes of 0.78 for all sites and 1.07 for direct PM2.5 sites and no significant differences in performance between Kalimantan and Sumatra. The mean satellite-derived PM2.5 over the entire 2014–2017 period (33.81 μg/m3 for all sites; 33.53 μg/m3 for direct PM2.5 sites) demonstrate a good agreement with the mean ground monitor PM2.5 (28.39 μg/m3 for all sites; 24.39 μg/m3 for direct PM2.5 sites). This comparison with ground monitor data indicates that the satellite-derived PM2.5 estimates are accurately capturing the high PM2.5 concentrations associated with biomass burning events, and can be used with a high degree of confidence for estimating the impacts of PM2.5 due to biomass burning over Indonesia.
Figure 2 (A) showsthe 2013–2017 average PM2.5 concentrations for the core fire season (months of August to October) and (B) all other months in Sumatra and Kalimantan. As shown in Fig. 2 (A), hotspots are apparent predominantly in the fire season, in particular in the southeastern part of Sumatra (Riau, Jambi and South Sumatra provinces) and in the provinces of Central and West Kalimantan
Figure 3 shows how PM2.5 concentrations vary across the year, indicating both the 5-year average and the concentrations for a dry (2015) and wet (2016) year. During the fire season (June to December), average PM2.5 concentrations in Sumatra and Kalimantan reach almost 60 μg/m3 (with a maximum monthly average PM2.5 concentration of 143 μg/m3 in October 2015), while for the non-fire season months, average concentrations are around 20 μg/m3. Note that these maps show PM concentrations as determined by all emission sources including fires on peat and mineral land, industry, traffic, etc. The difference between the average PM2.5 concentration during the 5 months wet season (the ‘counterfactual’) and the 7 months dry season is assumed to be the contribution of biomass burning comprising both fires in areas with a mineral soil and in areas with a peat soil.
Supplementary Materials Table A5 presents the contribution of peatland fires to the overall PM2.5 emission from biomass burning. On average, over 5 years, this contribution is 76% in Kalimantan and 86% in Sumatra. Subsequently, the population-weighted exposure to PM2.5 from peatland fires by province and by year are calculated. This is shown in Table 3. The table shows the average PM2.5 exposure of people living in each province.
Health effects from peatland fires
Health impacts of fires in Indonesian peatlands, in Sumatra and Kalimantan, are shown in Fig. 4. Using the CRF curve of , each year around 33,100 adults (95% CI = 31,200, 35,700) and around 2900 infants die prematurely because of air pollution resulting from peatland fires. This mortality rate corresponds to around 10% of total adult mortality (over 30 years old) in Sumatra and Kalimantan and 7.7% of the annual infant mortality. For a different perspective, the loss of life expectancy due to peatland fires at birth per inhabitant in Sumatra and Kalimantan is, respectively, 0.9 and 1.2 years (life expectancy at birth is around 71 years). These estimates were calculated using a life table analysis , considering province-specific data on the population and mortality distribution by age group (sourced from IHME’s GHDx database). In addition to premature deaths, we calculate the morbidity incidences related to smoke inhalation. On average, for each year between 2013 and 2017, our analysis has identified an additional 4390 respiratory hospital admissions, 635 thousand severe asthma symptom days among asthmatic children, and 8.9 million workdays lost in the working population.
Our study is the most comprehensive analysis of the long-term health effects of peatland fires in Indonesia carried out to date, in terms of the spatial and temporal cover and the application of CRFs to analyse the adult long-term mortality, infant deaths, and the additional morbidity effects associated with smoke inhalation in the exposed population. This improved knowledge leads to a better understanding of the negative health effects linked to peat drainage. Nevertheless, the analysis is subject to several important sources of uncertainty, related to both data and the models that we use. In terms of data, the main uncertainties pertain to: (i) the areas covered by peatlands; (ii) the occurrence of fires; and (iii) the estimation of PM2.5. In the last 20 years, 5 studies have published estimates of peatland occurrence in Sumatra and Kalimantan (for an overview, see ). For Sumatra, these estimates range from 5.6 million ha , to 9.6 million ha . For Kalimantan, estimates range from 4.8 million ha  to 6.7 million ha . In this study, the numbers of , are used, in line with the numbers formally adopted by the Government of Indonesia. Based on these numbers, the contribution of peatland fires to overall PM2.5 emitted from fires is estimated at 86% in Sumatra (average over 2013–217, range: 75% in 2017 to 94% in 2013) and 76% in Kalimantan (range: 62% in 2016 to 88% in 2013). This result is comparable to Kiely et al. , who estimate that, across Indonesia, peatland fires contributed 71% of total PM2.5 emissions from fires in Indonesia during September–October 2015. Our study focusses on two islands with extensive peat areas, which can explain why we find somewhat higher values. However, the uncertainty in peat area means for our study that there is uncertainty in the attribution of PM2.5 to peatland fires (as opposed to fires in mineral soil). Since the numbers on peat extent that we use are relatively low , we may underestimate the contribution of peatland fires to overall increases in PM2.5, and thereby underestimate the associated health effects.
Next, there is uncertainty regarding the occurrence of fires as measured with MODIS. The MODIS MCD64A1 product presents an estimate of the area burned in a given time period, however it is likely to underreport the actual burned area. In particular, small, burned areas (less than 100 ha) and burned areas in croplands and degraded grasslands are not all registered . The estimates of burned areas presented in this report are subject to these same uncertainties. As for the measurement of PM2.5, we are combining ground observations of PM2.5 from the Government of Indonesia with models of ambient PM2.5 based on remote sensing observations. As explained in Section 3.2, our model has a good fit with the observed data, having an R2 of 0.90 with PM2.5 measurement stations and an R2 = 0.77 for all sites. However, we apply a simplification in the allocation of PM2.5 emissions to peatland fire versus forest fires on mineral land. In particular, we do not consider SOA formation in distinguishing the ratio PM2.5 from peatland fires versus forest fires. The incomplete burning of peat soils means that SOA formation may, on a per hectare basis, be higher in peatlands compared to forests on mineral land (e.g. ). Furtheromore, the EF we use for peat combustion (9.04 g/kg) is conservative. Other studies which have done field measurements have found up to 29.6 g/kg  and 22.3 g/kg . Both aspects mean that we may underestimate the amount of ambient PM2.5 resulting from peatland fires (and overestimate the PM2.5 resulting from forest fires on mineral land). This means we may therefore underestimate the health effects of peatland fires.
In terms of health impact models, a key source of uncertainty lies in the selection of the concentration response functions that we use to assess health effects, and in our underlying input data assumptions. A first important assumption that we make is that the exposure to PM2.5, as measured over the period 2013–2017, is representative for the long-term exposure in Sumatra and Kalimantan. The period we study has two dry years (2014, 2015) two wet years (2016, 2017) and a year with average rainfall (2013). This can be considered representative for long-term rainfall conditions. Another key assumption is that PM2.5 emissions from other sources, in particular transportation, industry, energy supply and households, are similar in magnitude during the wet and dry season, so that the total increase in the dry season can be attributed to fires. There is no reason to assume any of these sources emits significantly more PM2.5 in the dry season compared to the wet season, however we cannot exclude that this may be the case. Furthermore, it may be that different rainfall patterns in the dry and wet season influence the rate at which PM2.5 is removed from the atmosphere affecting our assumption that the difference between the dry and the wet season PM2.5 concentrations can be solely attributed to biomass burning. The PM2.5 data for the year 2016, the year with the fewest peatland fires, show that the associated error is likely to be small (see Fig. 3). For 2016, the average Jan to May = 16.75 μg/m3, average Jan to Jul = 16.36 μg/m3, and average Jan to Dec = 16.66 μg/m3 – in other words, the year with the least peatland fires experienced fairly constant average PM2.5 concentrations throughout the year. An additional, important assumption is that a peak in PM2.5 concentration during the dry season can be averaged out over the year in order to assess health impacts; in other words that the health effects of a seasonal peak in PM2.5 are similar to those of a (lower) annual average increase in PM2.5. Since we are analyzing the long-term impacts of PM2.5 exposure, we believe the resulting potential errror to be small, and we note that this assumption is also made in various other studies analyzing the health impacts of forest fires, e.g. [42, 52]. Finally, we assume that the human toxicity of PM2.5 derived from peatland fires is in line with the average toxicity of PM2.5 (as related to health effects in concentration-response curves world-wide). However, various studies indicate that PM2.5 from biomass burning is relatively toxic, compared to other sources of PM2.5 such as road dust or aerosols from sea salt spray [53, 54].
Furthermore, there is uncertainty related to applying concentration-response curves to Indonesia, in particular because many of the underlying data are from Western Europe and the US, even though the study that we used for our base case estimate, , also includes data from China. This is reflected in the spread of the values that we present in Table 4. The accuracy of the application depends on the similarity of factors such as: (a) the chemical composition of the PM2.5; (b) the activity patterns and time spent outdoors; (c) underlying population disease profile and health status, including availability and utilization of health care services; (d) socioeconomic status; (e) age distribution, and (f) exposure to other pollutants, both outdoors and indoors. Of particular relevance is that Indonesia’s population is much younger and hospital admissions are much lower compared to the US and Western Europe. In addition, Indonesia and many other low- and middle-income countries have a different proportion of communicable diseases (around 16% of natural deaths) than the U.S. and Western Europe (around 5–6%). To explore the implications of using different CRFs, Table 4 below presents a sensitivity analysis, using three different choices for the CRF, with further details provided in Table A1.5 in the SM.
Finally, we need to stress that a limitation of our analysis is that we only examine premature mortality and morbidity effects of peatlands fires in Sumatra and Kalimantan. In years with severe fires, such as 2014 and 2015, smoke from peatland fires is spread over considerably larger areas, including other islands of Indonesia, Malaysia, Singapore and even parts of Thailand. Even though the health effects are concentrated in the two islands that we examine, mortality and morbidity effects do occur in other areas, both as a consequence of peatland fires in Sumatra and Kalimantan and due to peatland fires in other areas (e.g. Papua and Sulawesi) [55,56,57]. Consequently, our numbers are an underestimate of the overall health effects of peatland fires in Indonesia.
Our results can be compared to two other studies that assessed the health effects of peatland fires in Sumatra and Kalimantan. Crippa et al.  find that short-term exposure to air pollution resulting from peat fires in 2015 may have caused 11,880 (6153–17,270) excess mortalities across SE Asia and state that “the estimated deaths represent only a fraction of the overall premature fatalities due to long-term exposure to unhealthy air quality conditions”. They estimate that ~ 75,600 excess premature mortality would occur each year if the population received long-term exposure to the pollutant concentrations experienced in fall 2015 . Kiely et al. (2020)  analysed excess mortality from health effects related to peat fires in six dry years with extensive peat fires in the period 2004 to 2015. For 2015, they found that exposure to PM resulted in 44,040 excess deaths in SE Asia, of which 61% (about 27,000 deaths) occurred in Sumatra and Kalimantan . Our estimate for 2015 is 45,300 premature deaths across Sumatra and Kalimantan due to peatland fires, a value that is halfway between the estimate of the two studies [12, 55]. A main factor driving the differences between the studies is the counterfactual, i.e. the amount of PM2.5 that is added to the background PM2.5 concentration due to peatland fires (for instance,  assumed this to be 25 μg/m3, whereas our counterfactual is 18 μg/m3), but also differences in input data such as share of total PM2.5 emissions from peatland compared to other fires, population exposure-level, and demographic data play a role.
The use and management of Indonesian peatlands is challenging. Peatlands cover large areas (between 15 and 20 million ha), often located in remote areas, in particular in Sumatra, Kalimantan and Papua. These peatlands are an important economic asset, given that they can be used for agriculture [47, 58, 59]. At present, only small areas of peatlands are being used for crops that do not require drainage, such as sago. The drainage of peatlands for other crops brings significant negative externalities. In addition to adverse health effects from peatland fires, these include CO2 emissions, loss of biodiversity, and soil subsidence that will, over time, increase flood risks across Indonesian lowlands [60,61,62]. Importantly, fire occurrence in 2016 and 2017 was markedly lower, which may be related to higher rainfall in 2016 and 2017, compared to 2014 and 2015. Fire occurrence increased again in 2019 .
We show that the health effects of draining and converting peatlands for agricultural development are substantial. Annual premature adult mortality due to peatland fires ranges from 19,900 to 64,800 cases, depending upon the concentration-response curve used (Table 4). Clearly, also the lower end of the model results indicates a substantial health effect attributable to peatland fires (6% of adult natural mortality). Furthermore, peatland fires put additional stress on the country’s health care system, and result in lost workdays due to air pollution related sicknesses. Both aspects are economic costs for Indonesia. There have also been reports that reduced lung functioning due to air pollution increases population vulnerability to COVID-19 . The relation between COVID-19 and air pollution from peat and forest fires remains to be further studied, but this may provide further policy incentive to address peatland fires.
The Indonesian government has already banned farmers from dry-season crop residue burning, and such efforts are very important in mitigating health risks. In addition, health risks should be considered in land use planning and decision making on peatland uses. It is relevant to consider in which land use types the fires predominantly take place. Our data show that the least fire-prone areas are forests, which under natural conditions burn very seldomly (0.1% of the forest land burns each year, an average over the period 2013–2017, in Sumatra and Kalimantan; part of this land is burned deliberately to clear space for agriculture). On the other end of the spectrum, each year some 4% of agricultural lands and some 6% of grasslands burn. Forest plantations and perennial (plantation) crops take an intermediate position; some 2% is burned each year on average across Sumatra and Kalimantan. Around 61% of PM2.5 emissions related to biomass are from grasslands and shrublands, and 23% from plantations (annual croplands cover a much smaller area and contribute only 7%).
To reduce fire risks and associated health effects, a first priority is the rehabilitation of degraded grasslands and shrublands in peat, involving rewetting, revegetating, and controlling fire, cf. . Second, in the near to medium term, there is a need to phase out oil palm, acacia, coconut and other crops on peat that require drainage, and replace these by other peatland uses that do not require drainage. It is noted that low water tables in plantations extend to adjacent areas ; drainage for plantations desiccates adjacent areas up to several km from the plantation boundary. Hence, drainage in plantations also increases fire risks in other land use types in peat including in forests and grasslands.
In the future, fire risks from plantation agriculture in peat involving drainage will only increase. Drainage leads to soil subsidence. As shown in , most drained plantations in peat have an expected lifecycle of 1 to 3 rotations (25 to 75 years) before they are subsided to the level that seasonal flooding kills the trees (both oil palm and acacia are sensitive to high water tables). Subsequently, the plantation will be abandoned. Abandoned plantations have a very high fire risk: their drained status and ample dead biomass mean fires will propagate quickly in the dry season. When the area covered by abandoned plantations increases in the coming decades, fire risks will further increase. This again points to the need to develop new plantation models in peat, using paludiculture (no drainage) crops such as sago . This should be done as soon as possible, before ongoing peat subsidence and increasing flood and fire risks constrain the possibilities for paludiculture, since most paludiculture crops such as sago and jelutong are also sensitive to flooding and fire.
Indonesian peatlands have been drained for agricultural development for several decades. Although this has contributed to economic development through large-scale and, to a lesser degree, smallholder plantations, peat drainage is increasingly causing environmental externalities with adverse economic impacts. One of these is air pollution resulting from peatland fires. Peatland fires occur every year, even though their extent is much larger in dry (El Niño) years. Peatland fires, therefore, increase the long-term exposure of people to air pollution. We find that PM2.5 air pollution from peatland fires causes, on average, around 33,100 adults and 2900 infants to die prematurely from air pollution, in the islands of Sumatra and Kalimantan alone. In addition, peatland fires cause around 4390 additional hospital admissions related to respiratory diseases, 635 thousand cases of severe asthma attacks in children, and 8.9 million lost work-days each year in these two islands. Sumatra and Kalimantan have the country’s highest concentrations of smoke from peatland fires, but there will also be health impacts of peatland fires in other Indonesian islands and in neighboring countries, in particular Singapore and Malaysia . The main source of uncertainty is in the CRF that we use, with different CRFs leading to annual premature adult mortality ranging from 19,900 to 64,800 cases. Currently, the population of both islands is relatively young. With aging of the population over time, vulnerabilities to air pollution and health effects from peatland fires will increase. Hence, it is critical that health effects from peatland fires are acknowledged and considered in decision making, providing an additional element in support of current Government of Indonesia policies to reduce drainage-based agriculture in peatlands.
Availability of data and materials
The datasets during and/or analysed during the current study available from the corresponding author on reasonable request.
Information is available online: http://ghdx.healthdata.org/gbd-results-tool
Vetrita Y, Cochrane MA. Fire frequency and related land-use and land-cover changes in Indonesia’s peatlands. Remote Sens. 2020;12(1):5.
Stockwell CE, Jayarathne T, Cochrane MA, Ryan KC, Putra EI, Saharjo BH, et al. Field measurements of trace gases and aerosols emitted by peatland fires in Central Kalimantan, Indonesia during the 2015 El Niño. Atmos Chem Phys. 2016;16(18):11711–32.
World Bank and Government of Indonesia. Pilot ecosystem account for Indonesian peatlands. Washington, D.C.: Published by the World Bank, WAVES Program; 2018.
Cahyono SA, Warsito SP, Andayani W, Darwanto DH. Faktor-Faktor Yang Mempengaruhi Kebakaran Hutan Di Indonesia Dan Implikasi Kebijakannya. Jurnal Sylva Lestari. 2015;3:103–12.
Rasyid F. Permasalahan dan Dampak Kebakaran Hutan. Widyaiswara Network Journal. 2014;1:47–59.
WHO. World Health Organization global air quality guidelines: particulate matter (PM2.5 and PM10), ozone, nitrogen dioxide, sulfur dioxide and carbon monoxide. Geneva: World Health Organization; 2021.
WHO. Health risks of air pollution in Europe – HRAPIE project. Recommendations for concentration–response functions for cost–benefit analysis of particulate matter, ozone and nitrogen dioxide. Copenhagen: WHO Regional Office for Europe; 2013.
Hoek G, Krishnan RM, Beelen R, Peters A, Ostro B, Brunekreef B, et al. Long-term air pollution exposure and cardio- respiratory mortality: a review. Environ Health. 2013;12:43.
Burnett R, Chen H, Szyszkowicz M, Fann N, Hubbell B. Global estimates of mortality associated with long-term exposure to outdoor fine particulate matter. Proc Natl Acad Sci. 2018;115(38):9592–7.
Frankenberg E, Ds MK, Thomas D. Health consequences of Forest fires in Indonesia. Demography. 2005;42(1):109–29.
Uda SK, Hein L, Atmoko D. Assessing the health impacts of peatland fires: a case study for Central Kalimantan. Indonesia Environmental Science and Pollution Research. 2019;26:31315–27.
Kiely L, Spracklen DV, Wiedinmyer C, Conibear L, Reddington CL, Arnold SR, et al. Air quality and health impacts of vegetation and peat fires in Equatorial Asia during 2004–2015. Environ Res Lett 2020;15.
Sulong NA, Latif MT, Khan F, Amil N, Ashfold MJ, Wahab MIA, et al. Source apportionment and health risk assessment among specific age groups during haze and non-haze episodes in Kuala Lumpur. Malaysia Science of The Total Environment. 2017;601–602:556–70.
Ritung S, Wahyunto NK, Sukarman H, Suparto TC. Peta Lahan Gambut Indonesia Skala 1:250.000. Bogor: Balai Besar Litbang Sumberdaya Lahan Pertanian. 2011.
MOEF. Ministry of Environment and forestry Republic of Indonesia National Forest Reference Emission Level for REDD+ in the context of decision 1/CP.16 in directorate general of climate change. The Ministry of Environment and forestry: Indonesia. 2019.
Giglio L, Boschetti L, Roy DP, Humber ML, Justice CO. The collection 6 MODIS burned area mapping algorithm and product. Remote Sens Environ. 2018;217:72–85.
Hammer M S, van Donkelaar A, Li C, Lyapustin A, Sayer AM, Hsu NC, Levy RC, et al. Global Estimates and Long-Term Trends of Fine Particulate Matter Concentrations (1998-2018). 2020.
Kaufman YJ, Tanré D, Gordon HR, Nakajima T, Lenoble J, Frouin R, et al. Passive remote sensing of tropospheric aerosol and atmospheric correction for the aerosol effect. J Geophys Res Atmos. 1997;102(D14):16815–30.
Meteorologi B. Klimatologi, dan Geofisika (BMKG) Indonesia. Informasi Kualitas Udara BMKG Edisi Spesial 2020. Jakarta: BMKG; 2020.
MOEF. Laporan Tahunan PM. Direktorat Pengendalian Pencemaran Udara (AQMS KLHK), Ministry of Environment and forestry Republic of Indonesia: Indonesia; 2020.
Field RD,, van der Werf GR, Fanin T, Fetzer EJ, Fuller R, Jethva H, Levy R, Livesey NJ, Luo M, Torres O, , Worden HM. Indonesia 2015 fire and haze. Proc Natl Acad Sci 2016; 113 (33):9204–9209.
Khan F, Hamid AH, Rabin HA, Maulud KNA, Latif MT, et al. El Niño driven haze over the southern Malaysian peninsula and Borneo. Sci Total Environ. 2020;730:13909.
MOEF. Ministry of Environment and forestry Republic of Indonesia land cover maps, available on-line for 2015. Indonesia: The Ministry of Environment and Forestry; 2015a.
Ballhorn U, Siegert F, Mason M, Limin S. Derivation of burn scar depths and estimation of carbon emissions with LIDAR in Indonesian peatlands. Proc Natl Acad Sci. 2009;106(50):21213–8.
Mulyani A, Susanti E, Dariah A, Maswar W, Agus F. Basis data karakteristik tanah gambut di Indonesia. Kementerian Pertanian Indonesia: Badan litbang Pertanian; 2012.
Lestari PF, Muthmainnah DA, Permadi. Characterization of carbonaceous compounds emitted from Indonesian surface and sub surface peat burning. Atmospheric. Pollut Res. 2020;11(26):1465–72.
Christian TJ, Kleiss B, Yokelson RJ, Crutzen P J, Hao W. M, Saharjo B. H D E Ward Comprehensive laboratory measurements of biomass-burning emissions: 1. Emissions from Indonesian, African, and other fuels. J Geophys Res 2003;108: D23.
Sk A, Yokelson RJ, Wiedinmyer C, Alvarado MJ, Reid JS, Karl T, et al. Emission factors for open and domestic biomass burning for use in atmospheric models. Atmos Chem Phys Discuss. 2010;10:27523–602.
Burnett RT, Pope CA III, Ezzati M, Olives C, Lim SS, Mehta S, et al. An integrated risk function for estimating the global burden of disease attributable to ambient fine particulate matter exposure. Environ Health Perspect. 2014;122(4):397–403.
Vodonos A, AbuAwad Y, Schwartz J. The concentration-response between long-term PM2.5 exposure and mortality; a meta-regression approach. Environ Res. 2018;166:677–89.
Badan Perencanaan Pembangunan Nasional (BAPPENAS) Indonesia. Analisis Kesenjangan Antarwilayah Jakarta:BAPPENAS 2013.
WHOGlobal Health Estimates (GHE) 2016 Summary Tables. 2018(http://www.who.int/healthinfo/global_burden_disease/en/).
UN Inter-agency Group for Child Mortality Estimation. Levels and trends in child mortality. Paris: UNICEF; 2021.
Dong TTT, Hinwood AL, Callan AC, et al. In vitro assessment of the toxicity of bushfire emissions: a review. Sci Total Environ. 2017;603–604:268–78.
Helinor J, Johnston HJ, Mueller W, Steinle S, Vardoulakis S, Tantrakarnapa K, et al. How harmful is particulate matter emitted from biomass burning? A Thailand Perspective. Curr Poll Rep. 2019;5:353–77.
Aguilera R, Thomas Corringham T, Gershunov A, Tarik BT. Wildfire smoke impacts respiratory health more than fine particles from other sources: observational evidence from Southern California. Nat Commun. 2021;12:1493.
Crouse DL, Erickson AC, Christidis T, Pinault L, van Donkelaar A, Li C, et al. Evaluating the sensitivity of PM2. 5–mortality associations to the spatial and temporal scale of exposure assessment. Epidemiology. 2020;31(2):168–76.
Heft-Neal S, Burney J, Bendavid E. Burke M robust relationship between air quality and infant mortality in Africa. Nature. 2018;559(7713):254–8.
Arceo E, Hanna R, Oliva P. Does the effect of pollution on infant mortality differ between developing and developed countries? Evidence from Mexico City Econ J. 2016;126:257–80.
Loomis D, Castillejos M, Gold DR, McDonnell W, Borja-Aburto VH. Air pollution and infant mortality in Mexico City. Epidemiol. 1999;10:118–23.
Goyal N, Karra M, Canning D. Early-life exposure to ambient fine particulate air pollution and infant mortality: pooled evidence from 43 low- and middle-income countries. Int J Epidemiol. 2019;48:1125–41.
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:189–97.
Yunus F, Antaria R, Rasmin M, Mangunnegoro H. Asthma prevalence among high school students in East Jakarta, 2001, based on ISAAC questionnaire. Med J Indonesia. 2003;12(3):178–86.
Praptiningsiha CY, Lafond KE, Wahyuningrum Y, Storms AD. Healthcare-seeking behaviors for acute respiratory illness in two communities of Java, Indonesia: a cross-sectional survey. J Epidemiol Global Health. 2016;6:77–86.
Miller B, Hurley JF. Life table methods for quantitative impact assessments in chronic mortality. J Epidemiol Commun Health. 2003;57:200–6. https://doi.org/10.1136/jech.57.3.200.
Uda SK, Hein L, Sumarga E. Towards sustainable management of Indonesian tropical peatlands. Wetl Ecol Manag. 2017;25:683–701.
Koh LP, Miettinen J, Liew SC, Ghazoul J. Remotely sensed evidence of tropical peatland conversion to oil palm. Proc Natl Acad Sci. 2011;108(12):5127–32. https://doi.org/10.1073/pnas.1018776108.
MOEF. Peta Kesatuan Hidrologi Gambut. Ministry of Environment and forestry Republic of Indonesia. Jakarta; 2015b.
Kiely C, Spracklen DV, Wiedinmyer C, Conibear L, Reddington CL, Archer-Nicholls S, et al. New estimate of particulate emissions from Indonesian peat fires in 2015. Atmos Chem Phys. 2019;19:11105–21.
Jayarathne T, Stockwell CE, Gilbert AA, Daugherty K, Cochrane MA, Ryan KC, et al. Chemical characterization of fine particulate matter emitted by peat fires in Central Kalimantan, Indonesia, during the 2015 El Niño. Atmos Chem Phys. 2018;18:2585–600.
Wooster M, Gaveau DLA, Salim M, Zhang T, Xu W, Green DC, et al. New Tropical Peatland Gas and Particulate Emissions Factors Indicate 2015 Indonesian fires released far more particulate matter (but less methane) than current inventories imply. Remote Sens. 2018;10.
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 Health. 2014;13:105.
Park M, Soo Joo H, Lee K, Jang M, Kim SD, Kim I, et al. Differential toxicities of fine particulate matters from various sources. Nature. 2018;17007.
World Bank. Are All Air Pollution Particles Equal? – How Constituents and Sources of Fine Air Pollution Particles (PM2.5) Affect Health. Washington, DC 2020.
Gaveau DLA, Salim MA, Hergoualc'h K, Locatelli B, Sloan S, Wooster M, et al. Major atmospheric emissions from peat fires in Southeast Asia during non-drought years: evidence from the 2013 Sumatran fires. Sci Rep. 2014;4:6112.
Crippa P, Castruccio S, Archer-Nicholls S, Lebron GB, Kuwata M, Thota A, Sumin S, Butt E, Wiedinmyer C, Spracklen DV. Population exposure to hazardous air quality due to the 2015 fires in equatorial Asia. Sci Rep 2016;6:37074. https://doi.org/https://doi.org/10.1038/srep37074
Ruchi S, Rajasekhar B. Indoor human exposure to size-fractionated aerosols during the 2015 southeast Asian smoke haze and assessment of exposure mitigation strategies. Environ Res Lett 2017;12(11):114026. https://doi.org/https://doi.org/10.1088/1748-9326/aa86dd.
Murdiyarso D, Hergoualc’h K, Verchot LV. Opportunities for reducing greenhouse gas emissions in tropical peatlands. Proc Natl Acad Sci. 2010;107(46):19655–60. https://doi.org/10.1073/pnas.0911966107.
Miettinen J, Shi C, Liew SC. Two decades of destruction in Southeast Asia’s peat swamp forests. Front Ecol Environ. 2011;10(3):124–8. https://doi.org/10.1890/100236.
Hooijer A, Page S, Canadell JG, Silvius M, Kwadijk J, Wösten H, et al. Current and future CO2 emissions from drained peatlands in Southeast Asia. Biogeosciences. 2010;7:1505–14.
Thorburn CC, Kull CA. Peatlands and plantations in Sumatra, Indonesia: complex realities for resource governance, rural development and climate change mitigation. Asia Pacific Viewp. 2015;56(1):153–68. https://doi.org/10.1111/apv.12045.
World Bank. The cost of fire: an economic analysis of Indonesia’s 2015 fire crisis. Indonesia sustainable landscapes knowledge. Note no. 1. Washington, D.C.: World Bank; 2016.
Normile D. Indonesia's fires are bad, but new measures prevented them from becoming worse. Commentary. Science. 2019.
Wei L, Zhao-Wu T, Lei W, Ming-Li Y, Kui L, Ling Z, et al. Analysis of factors associated with disease outcomes in hospitalized patients with 2019 novel coronavirus disease. Chin Med J. 2020.
Murdiyarso D, Lilleskov E, Kolka R. Tropical peatlands under siege: the need for evidence-based policies and strategies. Mitig Adapt Strateg Glob Chang. 2019;24:493–505.
Wösten H, Hooijer A, Siderius C, Satriadi RD, Idris A, Rieley J. Tropical peatland water management modelling of the air Hitam Laut catchment in Indonesia. International Journal of River Basin Management. 2006;4(4):233–44.
Sumarga E, Hein L, Hooijer A, Vernimmen R. Hydrological and economic effects of oil palm cultivation in Indonesian peatlands. Ecol Soc. 2016;21(2).
Sheldon TL, Sankaran C. The impact of Indonesian Forest fires on Singaporean pollution and health. Am Econ Rev. 2017;107(5):526–9.
Sayer AM, Munchak LA, Hsu NC, Levy RC, Bettenhausen C, Jeong MJ. MODIS collection 6 aerosol products: comparison between Aqua’s e-deep blue, dark target, and “merged” data sets, and usage recommendations. J Geophys Res Atmos. 2014;119(24):13,965–89.
Levy RC, Mattoo S, Munchak LA, Remer LA, Sayer AM, Patadia F, et al. The collection 6 MODIS aerosol products over land and ocean. Atmos Meas Tech. 2013;6(11):2989–3034.
Hsu NC, Lee J, Sayer AM, Kim W, Bettenhausen C, Tsay SC. VIIRS deep blue aerosol products over land: extending the EOS long-term aerosol data records. J Geophys Res Atmos. 2019;124(7):4026–53.
Gupta P, Levy RC, Mattoo S, Remer L, Munchak LA. A surface reflectance scheme for retrieving aerosol optical depth over urban surfaces in MODIS dark target retrieval algorithm. Atmos. Meas. Tech. 2016;9(7):3293–308.
Sayer AM, Hsu NC, Lee J, Kim WV, Dutcher ST. Validation, stability, and consistency of MODIS collection 6.1 and VIIRS version 1 deep blue aerosol data over land. J Geophys Res Atmos. 2019;124(8):4658–88.
Diner DJ, Beckert JC, Reilly TH, Bruegge CJ, Conel JE, Kah RA, et al. Multi-angle imaging SpectroRadiometer (MISR) instrument description and experiment overview. IEEE Trans Geosci Remote Sens. 1998;36(4):1072–87.
Martonchik JV, Kahn RA, Diner DJ. Retrieval of aerosol properties over land using MISR observations. In: Satellite Aerosol Remote Sensing over Land. Berlin: Springer; 2019. p. 267–93.
Garay MJ, Witek ML, Kahn RA, Seidel FC, Limbacher JA, Bull MA, et al. Introducing the 4.4 km spatial resolution multi-angle imaging SpectroRadiometer (MISR) aerosol product. Atmospheric. Meas Tech. 2020;13(2):593–628.
Holben BN, Eck TF, Slutsker I, Tanré D, Buis JP, Setzer A, et al. AERONET—A federated instrument network and data archive for aerosol characterization. Remote Sens Environ. 1998;66(1):1–16.
Van Donkelaar A, Martin RV, Brauer M, Hsu NC, Kahn RA, Levy RC, et al. Global estimates of fine particulate matter using a combined geophysical-statistical method with information from satellites, models, and monitors. Environ Sci Technol. 2016;50(7):3762–72.
Eck TF, Holben BN, Reid JS, Dubovik O, Smirnov A, O’Neill NT, et al. Wavelength dependence of the optical depth of biomass burning, urban, and desert dust aerosols. J Geophys Res Atmos. 1999;104(D24):31333–49.
Giles DM, Sinyuk A, Sorokin MG, Schafer JS, Smirnov A, Slutsker I, et al. Advancements in the aerosol robotic network (AERONET) version 3 database – automated near-real-time quality control algorithm with improved cloud screening for sun photometer aerosol optical depth (AOD) measurements. Atmos. Meas. Tech. 2019;12(1):169–209.
Van Donkelaar A, Martin RV, Park RJ, Estimating ground-level PM2. 5 using aerosol optical depth determined from satellite remote sensing. J Geophys Res Atmos. 2006;111(D21).
Molod A, Takacs L, Suarez M, Bacmeister J. Development of the GEOS-5 atmospheric general circulation model: evolution from MERRA to MERRA2. Geosci Model Dev. 2015;8(5):1339–56.
Li T, Zhang Y, Wang J, Xu D, Yin Z, Chen H. All-cause mortality risk associated with long-term exposure to ambient PM2.5 in China: a cohort study. The lancet. Public Health. 2018;3:e470–7.
Zanobetti A, Schwartz J. Particulate air pollution, progression, and survival after myocardial infarction. Environ Health Perspect. 2007;11:76988.
Laden F, Schwartz J, Speizer FE, Dockery DW. Reduction in fine particulate air pollution and mortality: extended follow-up of the Harvard six cities study. Am J Respir Crit Care Med. 2006;173:667–72.
Schwartz J, Coull B, Laden F, Ryan L. The effect of dose and timing of dose on the association between airborne particles and survival. Environ Health Perspect. 2008;116:64–9.
Li C, Martin RV, van Donkelaar A, Boys BL, Hammer MS, Xu JW, et al. Trends in chemical composition of global and regional population-weighted fine particulate matter estimated for 25 years. Environ Sci Technol. 2017;2017a. https://doi.org/10.1021/acs.est.7b02530.
Yin P, Brauer M, Cohen A, Burnett RT, Liu J, Liu Y, et al. Long-term fine particulate matter exposure and nonaccidental and cause-specific mortality in a large national cohort of Chinese men. Environ Health Perspect. 2017;125(11):117002.
Masripatin N. Cadangan karbon pada berbagai tipe hutan dan jenis tanaman di Indonesia. Kementerian Kehutanan, Badan Penelitian dan Pengembangan Kehutanan, Pusat Penelitian dan Pengembangan Perubahan Iklim dan Kebijakan 2011.
Khalid H, Zin ZZ, Anderson JM. Quantification of oil palm biomass and nutrient value in a mature plantation. I. Above-ground biomass. J Oil Palm Res. 1999;11(1):23–32.
Albrecht A, Kandji TS, et al. Agriculture Ecosyst Environ. 99(1–3):15–27.
Marjolein Lof, Wageningen Univeristy, assisted with the spatial data analyses.
This research has been funded through regular work programmes of the authors.,
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
A1.1 Occurrence of Peatland fires. Table A1. Burned peat and mineral lands by land cover class, 2013 to 2017 (hectare). A1.2 Analysis of monthly PM2.5 concentrations. PM2.5 concentrations are estimated using satellite Aerosol Optical Depth (AOD). AOD is a measure of the extinction (scattering and absorption) of light by aerosols in a column of air from the Earth’s surface up to space. The AOD is used to represent the abundance of total aerosol particles in an atmospheric column . To analyze AOD, we use data from three satellite instruments: twin MODIS (MODerate resolution Imaging Spectroradiometer) instruments and the MISR (Multi-angle Imaging Spectroradiometer) instrument . We retrieve AOD with two algorithms that process MODIS radiances on both the Terra and Aqua satellites: Dark Target (DT) and Deep Blue (DB). The DT retrieval algorithm , is designed to retrieve AOD over dark surfaces (e.g. vegetated land surfaces and dark soils). The DB retrieval algorithm , uses blue wavelength measurements where the surface reflectance over land is typically much lower than at longer wavelengths, allowing for the retrieval of aerosol properties over both bright and dark surfaces. This study uses the recently released collection 6.1 of the MODIS retrieved AOD products, which include spatial resolution of 10 km and several updates to the DT , and DB algorithms [71, 73]. The MISR instrument , retrieval algorithm uses the same-scene multi-angular views provided by the nine view- angles to solve for surface and top-of-atmosphere reflectance contributions, providing AOD retrievals over bright and dark land surfaces without absolute surface reflectance assumptions . Specifically, we use AOD retrieved from the recently released MISRv23 algorithm , which provides AOD at a spatial resolution of 4.4 km. Data from MODIS and MSIR is resampled at 4.4 km resolution and combined based on their relative uncertainties with Aerosol Robotic Network (AERONET) AOD ground measurements [17, 77, 78]. AERONET is a global sun photometer network established by NASA and PHOTONS (PHOtométrie pour le Traitement Opérationnel de Normalisation Satellitaire) . This study uses AOD at 550 nm from level 2 of the version 3 AERONET data . To convert the combined AOD to surface PM2.5 concentrations, we use the simulated ratio of total column AOD to surface PM2.5. The ratio of total column AOD to surface PM2.5 is a function of the factors that relate PM2.5 mass to satellite observations of AOD (e.g., aerosol size, aerosol composition, diurnal variation, relative humidity, and the vertical structure of aerosol extinction . We simulate this relationship using the GEOS-Chem chemical transport model. A full description of the GEOS-Chem simulation used can be found in . We use v11–01 of GEOS-Chem, and our simulation is driven by assimilated meteorological data from the MERRA-2 Reanalysis of the NASA Global Modeling and Assimilation Office (GMAO) . The simulation is conducted for the years 2013 to 2017 with 47 vertical layers at a spatial resolution of 0.5° × 0.625° (~ 50 km × 60 km) across all of Asia (including Indonesia). The top of lowest model layer is ~ 100 m. Anthropogenic emissions of aerosols and their precursors are provided by the MIX inventory . Biomass burning emissions are provided for individual years by the GFED4 open fire emissions inventory. A1.3 Emission factors. Table A2. Emission factors for above ground biomass (kg PM2.5 per kg dry biomass). Table A3. Above ground biomass per hectare. Table A4. Emissions from peat biomass burning. Table A5. Emissions from biomass burning, per hectare. A1.4. Concentration-Response functions. The meta-analysis study included 53 cohort studies, 39 studies from North America, eight from Europe, and six from Asia . The mean concentration of PM2.5 across the studies was 15.7 μg/m3, with higher concentrations observed in the Asian studies with a mean of 30.5 μg/m3. The authors’ analysis indicated robust association between PM2.5 and premature mortality. At the mean concentration of 15.7 μg/m3, a 10 μg/m3 increase in PM2.5 was associated with a 10.3% (95% confidence interval (CI) of 9.7 to 11.1%) increase in adult premature mortality. This risk function was used for the base case estimates since it included the greatest number of studies in total, with a few from Asia. To examine the sensitivity of the results to the choice of this study, adult mortality estimates were generated from three other cohort studies (Section 4.3). The health risk was applied to the mortality rate in the population over 30 years old (Section 4.1.3) . A second estimate of adult mortality by Burnett et al. (2018)  involved a meta-analysis of 41 cohort studies from around the world while allowing for flexibility in determining the shape of the CRF. The derived CRF was an improvement over previous risk functions used for the GBD estimates, since it added studies from countries with very high concentrations of PM2.5 such as China, and relied solely on air pollution studies (previous CRF for GBD incorporated other combustion sources of PM2.5 such as secondhand smoke and household air pollution). In addition, the CRF is restricted to mortality from non-communicable diseases (NCD) and lower respiratory illness (LRI), rather than all-causes, since these two outcomes have been closely linked to PM2.5 exposure in both mortality and morbidity studies. The resulting CRF is non-linear and fairly complex statistically. The Burnett et al. (2018) risk function, , was applied to the NCD + LRI mortality rate in the population older than 25 years (Section 4.1.3). Figure A2.1 shows the shape of the Burnett et al. (2018) function, along with the other studies used for the CRF. A third analysis was based on Crouse et al. (2020) , who utilized census data on 2.4 million Canadian adults to examine non-accidental and cause-specific mortality between 2001 and 2011. PM2.5 exposures were assigned to participants’ residence using satellite-based estimates. Of particular interest, the study compared impacts of using three different average exposure periods (1, 3, and 8 years), so it is relevant to the five-year average of concentrations due to peatland fires. For a 10 μg/m3 change in PM2.5 the associated with risk estimates for non-accidental mortality of approximately 11% (95% CI = 8, 13%), 20% (95% CI = 17, 23%) and 23% (95% CI = 20, 27%), respectively, for one-, three- and eight-years of average prior exposure. Existing studies indicate the effects of exposure changes can be experienced within one or 2 years [84,85,86], but , found larger risks when even longer exposure times were considered. For the sensitivity analysis, the risk associated with 3 years of exposure was used to approximate the impact of the five-years of exposure from the peatland fires, and the risk associated with 1 year of exposure was used for calculating the impact of the single year (2015). The Crouse et al. (2020) health risk , was applied to the mortality rate in the population over 30 years old. The fourth, and final, study used in the sensitivity analysis is based on the recommendations from WHO’s HRAPIE (2013) review , which was based primarily on an earlier meta-analysis . The estimated risk for all-cause mortality (for the population 30+) was based on eleven studies available at the time. The resulting risk was 6% (95% CI = 4, 8%). As shown in Fig. 6, the CRF from Vodonos et al. (2018)  lies in between those of the three other studies. Fig. A1. Approximate change in mortality risk vs PM2.5 concentration for different published studies. The risk estimates from two recent cohort studies conducted in China are presented for comparison with the studies used for the base case and sensitivity analysis. Li et al. (2018) examined a cohort of over 13,000 participants from the Chinese Longitudinal Healthy Longevity Survey . PM2.5 exposures were assigned to residential locations and derived from remote sensing satellite data at a 1 km grid resolution. The median concentration was 50.7 μg/m3. The estimated all-cause mortality risk for a 10 μg/m3 change in PM2.5 was 8% (95% CI = 6, 9%). Another study  (which was included in the meta-analyses of Vodonos and Burnett cited above) analysed a cohort of 190,000 men randomly selected from 45 of China’s 145 Disease Surveillance Points around the country. PM2.5 concentrations were developed by combining data from remote sensing satellites, chemical transport model and ground-based monitors. With a mean PM2.5 concentration of 43.7 μg/m3, the non-accidental mortality risk estimate for a 10 μg/m3 change in PM2.5 was 9% (95% CI = 8, 9%). A1.5. Sensitivity analysis. Table A6. Sensitivity analysis of alternative assumptions for adult and infant mortality.
About this article
Cite this article
Hein, L., Spadaro, J.V., Ostro, B. et al. The health impacts of Indonesian peatland fires. Environ Health 21, 62 (2022). https://doi.org/10.1186/s12940-022-00872-w