Environmental Research Letters LETTER • OPEN ACCESS Smoke pollution disrupted biodiversity during the 2015 El Niño fires in Southeast Asia To cite this article: Benjamin P Y-H Lee et al 2017 Environ. Res. Lett. 12 094022 View the article online for updates and enhancements. Related content - Singapore’s willingness to pay for mitigation of transboundary forest-fire haze from Indonesia Yuan Lin, Lahiru S Wijedasa and Ryan A Chisholm - Impact of the June 2013 Riau province Sumatera smoke haze event on regional air pollution Sheila Dewi Ayu Kusumaningtyas and Edvin Aldrian - Public health impacts of the severe haze in Equatorial Asia in September–October 2015: demonstration of a new framework for informing fire management strategies to reduce downwind smoke exposure Shannon N Koplitz, Loretta J Mickley, Miriam E Marlier et al. This content was downloaded from IP address 97.113.180.199 on 25/02/2018 at 20:56 Environ. Res. Lett. 12 (2017) 094022 https://doi.org/10.1088/1748-9326/aa87ed LETTER OPEN ACCESS Smoke pollution disrupted biodiversity during the 2015 El Niño fires in Southeast Asia RECEIVED 17 April 2017 Benjamin P Y-H Lee1,2,3 , Zoe G Davies1 and Matthew J Struebig1,3 REVISED 1 17 August 2017 ACCEPTED FOR PUBLICATION 2 23 August 2017 3 Durrell Institute of Conservation and Ecology (DICE), School of Anthropology and Conservation, University of Kent, Canterbury, CT2 7NR, United Kingdom Centre for Urban Greenery and Ecology, National Parks Board, 1E Cluny Road, S259601, Singapore Both authors contributed equally to this work. Authors to whom any correspondence should be addressed. PUBLISHED 20 September 2017 E-mail: benjamin_lee@nparks.gov.sg and m.j.struebig@kent.ac.uk Keywords: acoustic indices, haze, soundscape ecology, forest fire, Singapore, environmental change, bioacoustics Original content from this work may be used under the terms of the Creative Commons Attribution 3.0 licence. Any further distribution of this work must maintain attribution to the author(s) and the title of the work, journal citation and DOI. Supplementary material for this article is available online Abstract Forest and peatland fires during the 2015 El Nĩno drought were amongst the worst on record in Southeast Asia. They were a major contributor of carbon emissions across the region, with the associated smoke-induced haze causing an air pollution crisis that affected millions of people. We present evidence of air pollution impacts on biodiversity. Using daily acoustic recordings in central Singapore, we monitored the dawn chorus before, during and after the haze event. We demonstrate that levels of ecological community acoustic activity dropped dramatically during the haze, and that this decline was significantly associated with levels of air pollution considered ‘unhealthy’ to the human population. Acoustic disruption was apparent across four common indices of soundscape activity, with only a partial recovery to pre-haze levels observed 16 weeks after the smoke had dissipated. These impacts on ecological communities were likely to be even more severe closer to the fires, where air pollution levels were reported to be 15-fold greater than those recorded in Singapore. Our results indicate that large-scale air pollution crises may have hitherto underestimated and potentially far-reaching impacts on biodiversity, especially in parts of the world prone to extensive forest fires. Introduction Air pollution is increasing globally (World Health Organization 2016). While the human health, stress and mortality impacts associated with exposure to air pollutants, such as sulphur dioxide and particulate matter, are widely acknowledged as significant (World Health Organization 2016, Doherty et al 2009), the effects of air pollution on wildlife are much less understood. Previous studies indicate that animals respond negatively to pollution in similar ways to humans in terms of oxidative stress (Isaksson 2010), but it remains unclear whether these impacts translate to lowering of activity or extirpation for particular animal species. This is especially the case for adverse pollution episodes. Between September and October 2015, much of maritime Southeast Asia was covered with a smokeinduced haze that originated from fires burning in © 2017 The Author(s). Published by IOP Publishing Ltd Indonesia’s forests and peatlands. Although uncontrolled fires occur annually in the region, those in 2015 were exacerbated by a prolonged drought caused by the El Nino–Southern ̃ Oscillation and Indian Ocean Dipole (Tacconi 2016, Huijnen et al 2016). An average of 11.3 Tg CO2 was emitted each day during September and October 2015, a figure that exceeded the fossil fuel CO2 emissions (8.9 Tg CO2 per day) across the European Union (Huijnen et al 2016). The transboundary haze event was so serious that it forced airports and schools to shut in parts of the region, at an estimated economic cost of ∼US$16.124 billion for Indonesia alone (World Bank 2016). Throughout the two months, the resulting air pollution regularly reached ‘unhealthy’ or ‘very unhealthy’ levels for humans, as measured by the Pollutant Standards Index (PSI) in Indonesia, Malaysia and Singapore (Tacconi 2016). Harrison et al (2016) reviewed the potential impacts of the fire and resultant haze on Kalimantan’s Environ. Res. Lett. 12 (2017) 094022 Table 1. Summary of the acoustic monitoring sample periods used to assess the impact of the 2015 haze event on ecological communities in central Singapore. Sample periods Dates Number of days Pollutant Standards Index (PSI) range Before During 3 week after 16 weeks after 11 January−1 February 2015 24 September−28 October 2015 16 November−14 December 2015 29 February−14 March 2016 22 21 20 15 31−68 97−267 28−53 37−55 biodiversity, speculating that the fauna of peat-swamp forests declined due to habitat loss, toxic particulate inhalation, drought and much reduced sunlight, and increased river acidity. To further examine impacts on biodiversity we used the 2015 haze event to examine ‘soundscape’ data from a wildlife-monitoring programme, which recorded acoustic activity each morning on an overpass in Singapore’s central forest block. Singapore is located approximately 300 km away from the locations of the fires, with the prevailing winds blowing the smoke over from Indonesia into the citystate. Acoustic recordings were made for two hours each morning after dawn for three weeks during the haze, and were compared to data available from eight months before the haze began, as well as three and 16 weeks after the smoke had dissipated. For each morning, we calculated acoustic activity using four indices that measure different aspects of the soundscape. The basic premise underpinning the application of acoustic indices is that ecological communities with more vocal species are characterised by greater acoustic diversity, and that biodiversity and ecosystem functioning are positively correlated with acoustic diversity (Gage et al 2001, Ji et al 2007, Pijanowski et al 2011). These validated metrics derived from acoustic monitoring can therefore be used to detect periods of environmental change (Sueur et al 2014). Materials and methods Monitoring of air pollution We obtained air pollution data, measured by the Pollutant Standards Index (PSI), from the National Environment Agency (2016) of Singapore. Hourly PSI values are available for five regions on the island (north, south, east, west and central). We used PSI values from central Singapore to best reflect the location of the overpass. The PSI is derived from the average concentrations of six common air pollutants: (i) fine particulate matter (PM2.5 ); (ii) particulate matter (PM10 ); (iii) sulphur dioxide; (iv) carbon monoxide; (v) ozone; and, (vi) nitrogen dioxide. PSI values up to 50 are considered ‘good’ in relation to the human population; 51−100 ‘moderate’; 101−200 ‘unhealthy’; 201−300 ‘very unhealthy’; and, >300 as ‘hazardous’. During the peak of the haze event, PSI values reached 267, and residents were advised by the Singapore Government to avoid prolonged or strenuous outdoor physical exertion. We compiled PSI values for our four acoustic monitoring sample periods in relation to the haze event: ‘before’; ‘during’; ‘3 weeks after’; and, ‘16 weeks after’ (table 1). These periods were identified according to distinct shifts in PSI values and public health advisory announcements on air quality, as well as availability of quality recordings. As an equatorial island, Singapore is largely aseasonal, although it tends to experience two monsoons: December to March, and June to September. By sampling 16 weeks after the event, we were able to quantify acoustic activity mid-way through the northeast monsoon in a similar period sampled before the onset of haze. We used the 0800 hrs PSI value for each sampling day as this represented the midpoint of our two-hour morning recordings. Study system Acoustic data were collected as part of a wildlifemonitoring programme on the ‘EcoLink’ wildlife overpass (1.357◦ N, 103.784◦ E) in central Singapore. The overpass is 62 m in length and 50 m in width, Acoustic recording, processing and statistical analyses and serves to re-connect two tropical lowland rain- We deployed a Song Meter SM2BAT+ unit (Wildlife forest reserves, Bukit Timah Nature Reserve (163 ha) Acoustics Inc., USA) at the centre of the overpass with a and Central Catchment Nature Reserve (c. 2000 ha); SMX-II weatherproof microphone attached to a 1.8 m protected areas that have been separated by the Bukit etre aluminium pole. The SMX-II is an omnidirecTimah Expressway for 30 years. Construction and tional and audio broadband microphone recording planting of native vegetation on the bridge were com- frequencies of up to 48 000 Hz, and with an estipleted in December 2013. The likely fauna contributing mated range of ca. 60 m radius. Our acoustic data to the soundscape include birds and Orthoptera. At therefore covered the overpass, as well as ca. 30 m of least 55 bird species have been identified from the forest either side. The recordings comprised both ambisite via mist-netting (Yifei Chung, unpublished data), ent sounds from the ecological community, as well and Nisitrus sp. crickets and Conocephalus sp. katydids as non-biological phenomena, including those arishave been encountered (Ming Kai Tan, unpublished ing from human activities (Pijanowski et al 2011). data). In addition, opportunistic camera trapping has The pole extended the microphone above vegetation detected macaques (Macaca fascicularis) and other to provide a clear recording space without obstrucsmall mammals using the overpass. tion from foliage or branches. The microphone was 2 Environ. Res. Lett. 12 (2017) 094022 checked and tested regularly, with the foam windscreen changed after sustained periods of rain. The Song Meter unit was configured to record in an uncompressed.wav format at a sampling rate of 24 000 Hz and 16 bits to target audible sound (i.e. up to ca. 11 025 Hz, the highest frequency of bird song in the region). A high pass-filter was applied at 180 Hz to attenuate excessive low-frequency sounds so that bird vocalisations could be recorded adequately. The device was placed in a locked box, from which the microphone extended, and secured to a tree trunk. Recordings were scheduled for between 0700 and 0900 hrs each day to capture the dawn chorus comprising bird vocalisations and insect calls. We checked all recordings within each file to assess levels of rain. Any recordings from mornings with heavy rain were subsequently excluded from the study, as bird and insect acoustic activity declines in rainy conditions (Robbins 1981) and the detection range of the microphone is reduced when wet. One morning was removed as a result. Each two-hour recording was split into 12 files of 10 min duration using WavePad v6.37 (NCH Software), and acoustic indices calculated for each 10 min segment. Sensitivity analyses using segments of various lengths (10, 30 and 60 min) for all four soundscape indices revealed that 10 min segments gave sufficient amount of variation across the duration of recordings for statistical comparisons to be made (supplementary information figure S1 available at stacks.iop.org/ERL/12/094022/mmedia). Calculation of acoustic indices We generated four acoustic indices for the Singapore soundscape recordings: • The Bioacoustic (BA) index is designed to assess the relative abundance and composition of bird communities, and is a function of the acoustic level and the number of frequency bands used by animals (Boelman et al 2007). The index is derived from the area under each curve in all frequency bands, and so the area values are a function of both activity level and the number of frequency bands used by vocalising species, but in particular, avifauna. Here, calculations targeted the ecological frequency range of 2000−11 025 Hz (anthropogenic sound is typically <2000 Hz). • Acoustic Complexity (AC) directly quantifies acoustic intricacy by calculating the variability of recording intensities, despite the presence of constant anthrophony (human-generated noise) in the background (Pieretti et al 2011, Farina and Pieretti 2014). The index correlates well with the number of bird calls in a community (Pieretti et al 2011), has been used to describe avian soundscapes (Farina et al 2011), to explore the association between these soundscapes and vegetation complexity (Farina and Pieretti 2014), and to investigate the influence of traffic noise on the singing of a bird community (Pieretti and Farina 2013). 3 • Acoustic Diversity (AD) is typically used as a direct proxy for species diversity, since it quantifies the number of unique sounds at unique frequencies, and these different sounds are largely emitted by different taxa (Villanueva-Rivera et al 2011). It is based on the Shannon index (Villanueva-Rivera et al 2011), and is derived by dividing the frequency spectrogram into bins and quantifying the proportion of sounds in each bin above a volume threshold (Sueur et al 2014, Villanueva-Rivera et al 2011). We applied this index to 1000 Hz bins. • The Normalised Difference Soundscape Index (NDSI) reflects the relative contribution of ecological (biophonic) and anthrophonic sounds to the soundscape. The calculation involves segregating the spectral profile into two main frequency bands comprising anthrophonic (1000−2000 Hz) and ecological (2000−11 025 Hz) components (Ji et al 2007), before calculating the ratio of the two. The value therefore ranges from −1 to +1, with low and high values indicating the prevalence of anthrophony or biophony respectively. Acoustic data were calculated in R version 3.2.4 using the ‘multiple_sounds’ function in the package ‘soundecology’ (Villanueva-Rivera and Pijanowski 2013), as well as the packages ‘seewave’, ‘ineq’, ‘tuneR’ and ‘vegan’ (using signal settings in S1 table). Mean values of each index were calculated across the twelve 10 min sound files recorded each morning, thereby producing a single average value per morning. The distribution of acoustic data across each sample period was checked for homoscedasticity (Levene’s test), visualised using the R package ‘ggplot2’ to generate violin plots, and the residuals checked for normality. The variance in mean acoustic data for each index was unequal between sample periods. Kruskal-Wallis tests were thus used to determine if the indices varied significantly between sample periods, and pairwise Mann-Whitney tests were applied post-hoc to determine where differences occurred. A standard Bonferroni correction was applied to the Mann-Whitney tests to control for Type 1 error due to multiple testing. General linear models were then used to determine whether variation in each acoustic index across the entire study period could be predicted by PSI. Since samples collected after the haze were clearly not independent from those collected before or during the event, we report regressions restricted to the ‘before’ and ‘during’ data. However, our findings remain consistent even when all data are included in the models. Results Recordings were available for a total of 78 mornings between January 2015 and March 2016. Our analyses reveal that daily averages for all four acoustic indices corresponded clearly and significantly to the onset and, Environ. Res. Lett. 12 (2017) 094022 Pollutant Standards Index Bioacoustic Index, BA 90 Acoustic Complexity, AC b b 0 50 100 150 200 250 R2 = 0.54 p < 0.001 70 60 50 H = 55.46 p < 0.001 40 a b c ac R2 = 0.48 p < 0.001 8600 8500 8400 H = 25.10 p < 0.001 8300 2.2 a b b c R2 = 0.61 p < 0.001 2.0 1.8 1.6 H = 52.17 p < 0.001 1.4 1.00 Normalized Difference Soundscape Index, NDSI a 80 8700 Acoustic Diversity, AD a a b b b R2 = 0.48 p < 0.001 0.75 0.50 0.25 0.00 H = 35.80 p < 0.001 –0.25 BEFORE DURING 3 WEEKS AFTER 16 WEEKS AFTER 0 50 100 150 200 250 Pollutant Standards Index unhealthy good moderate very unhealthy Figure 1. Daily variation in dawn chorus acoustic indices, and their association with air pollution levels, in Singapore before, during and after the 2015 El Nino ̃ forest and peatland fire smoke-induced haze event. Left panels: violin plots show variation (median, range, kernel density; 25th–75th percentiles) in daily measures of dawn chorus acoustic activity ‘before’ (22 sample days), ‘during’ (21 sample days) and ‘after’ (20 and 15 sample days at 3 and 16 weeks respectively). Different letters denote significant differences between sample periods, based on Kruskal-Wallis and post-hoc Mann-Whitney tests. Right panels: regressions (+/− 95% CI, indicated by shaded area around line of best fit) of daily average acoustic indices against Pollutant Standards Index (PSI) values. Due to data non-independence issues associated with the two ‘after’ sample periods, the regressions are calculated from data ‘before’ and ‘during’ the haze event. Dashed lines show the public health advisory PSI thresholds. to some degree, the passing of the air pollution crisis (figure 1). The greatest acoustic activity was pre-haze. Subsequently, average levels of acoustic activity across all indices declined substantially (by up to 37.5%) during the haze event, a disruption that lasted for at least 4 three weeks. A near complete recovery to pre-haze levels was evident after 16 weeks for the acoustic complexity and bioacoustic indices, two well-established indicators of bird vocalisation activity (Pijanowski et al 2011, Towsey et al 2014). However, acoustic diversity Environ. Res. Lett. 12 (2017) 094022 and normalized difference soundscape index, indices that reflect broader patterns of ecological activity relative to anthropogenic sounds (Gage et al 2001, Sueur et al 2008, Pijanowski et al 2011, Fuller et al 2015), remained low after 16 weeks. This suggests that some components of the ecological community continued to be absent or torpid for at least four months after the smoke dissipated. All acoustic indices correlated negatively and significantly with the average daily PSI readings for central Singapore, and were particularly low when air pollution reached levels deemed unhealthy for humans (>100 PSI; figure 1). Discussion We provide a unique evaluation of the effects of the 2015 Southeast Asia forest and peatland fires on biodiversity, demonstrating that the transboundary haze crisis had a clear negative impact on the ecological community in Singapore. To the best of our knowledge, no other study has investigated the impacts of air pollution on biodiversity in this way. Previous studies have focused on the autecological responses of particular species to air pollution (e.g. great tits, Parus major (Gorissen et al 2005), yet acoustic indices give us a valuable and, importantly, rapid insight into the entire audible ecological community, which is not restricted to individual species. Both the bioacoustic and acoustic complexity indices we used capture mainly bird vocalisations, with the former index being a reliable indicator of avian abundance (Boelman et al 2007) and the latter serving as a proxy for singing activity (Pieretti et al 2011). Our findings therefore imply that while bird abundance and associated activity were substantially impacted by the pollution event, both recovered within four months post-haze (figure 1). Conversely, the acoustic diversity index and normalised difference soundscape index showed little sign of recovery over this time period. In biodiversity assessments, acoustic diversity is renowned for being a reliable surrogate for overall (vocal) faunal richness, because it quantifies the number of unique sounds made at different acoustic frequencies, which are typically made by different species (VillanuevaRivera et al 2011, Depraetere et al 2012). Low acoustic diversity values therefore indicate that sound is being recorded at fewer frequencies across all sound emitting taxa. The drop in NDSI, an index characteristically high and stable in undisturbed environments (Gage et al 2001), also reflects this trend. Potential causes and mechanisms of acoustic activity change There are several possible explanations for our findings, which may not be mutually exclusive. First, the drop in acoustic activity might signal a temporary cessation in ecological activity during the dawn chorus. For example, in temperate regions the singing 5 behaviour of male great tits is known to be reduced and less complex at sites near to heavy metal air pollution sources (Gorissen et al 2005). In nearby Borneo, gibbons are reported to sing less during periods of poor air quality (Cheyne 2008). However, this seems unlikely to be the sole explanation of the soundscape patterns observed, given that acoustic recovery back to pre-haze levels was only demonstrated by two of the four indices. Second, smoke-induced air pollution may have caused direct mortality in key components of the ecological community. Although understudied, terrestrial vertebrates are predicted to suffer from the same direct effects of air pollution as humans, including respiratory diseases, hypoxia, irritation of the eyes and skin, increased stress and death (Isaksson 2010, Kimura et al 1988). During the 1997 El Nino, ̃ 90% of people in Indonesian villages exposed to similar levels of pollution as Singapore in 2015 had severe respiratory problems from inhaling toxic smoke (Kunii et al 2002). Third, mortality might have arisen indirectly via reduced fitness or foraging success, since airborne particles drastically increases atmospheric attenuation of light and sound (Hardy et al 2001, Pye 1971), and thus compromises the ability of animals to forage. This is particularly likely to be the case for flying and site-dependent species such as birds. Fourth, haze could reduce the abundance/size of prey species or affect plant phenology (see Harrison et al 2016), thereby driving bottom-up trophic cascades (Zvereva and Kozlov 2010). Fifth, there is scope for some seasonal variation in ecological activity to be occurring in the dataset although, given the highly aseasonal nature of Singapore and the scheduling of our sample periods, this is unlikely to have contributed to the observed trends. Finally, it is possible that our findings reflect a postponement of ecological activity to outside of the two-hour window we recorded at dawn. Delaying activity to sub-optimal times in this way would still be expected to have long-term repercussions for individual fitness and phenology (McNamara et al 2011). Nevertheless, this hypothesis is not well supported by anecdotal observations during the haze (i.e. forests were uncharacteristically quiet). Conclusion While we can only speculate about the precise mechanisms underpinning the reduction in acoustic activity during the 2015 haze crisis, the trends observed in Singapore could have far reaching implications for biodiversity if realised across pollution-stricken Asia. Ecological disruption across all taxa is likely to be greatest nearer to smoke sources. For example, PSI readings exceeded 2000 in Indonesian Borneo during the 2015 fires (Chan 2015), which is 15-times greater than the levels recorded in Singapore. Without a major concerted international effort to tackle the causes of forest and peatland fires, annual air pollution events Environ. Res. Lett. 12 (2017) 094022 are expected to continue and worsen in Southeast Asia (Tacconi 2016, Huijnen et al 2016). Indeed, high levels of transboundary air pollution are not restricted to drought years and, in fact, occur almost annually (Gaveau et al 2014). The long-term prognosis for the region’s imperilled biodiversity is therefore poor, particularly in burnt areas, which already support fewer taxa and depauperate wildlife populations compared to forests elsewhere (e.g. Cleary et al 2006). Further research into the effects of air pollution on tropical biodiversity is thus warranted. The funding and human capacity required to effectively document and monitor biodiversity in tropical countries, in addition to significant logistical challenges on-the-ground, are major barriers to large-scale implementation of traditional ecological field surveys (Harrison et al 2012). There is therefore a well justified need for rapid biological monitoring techniques, especially in situations in which health and safety is a concern. The capacity of acoustic techniques to yield information on species identity and turnover is currently still limited. However, as technological capabilities improve, it is highly probable that key vocalising taxa, such as those of conservation concern, will be identifiable. Such acoustic techniques therefore offer a practical and relatively cheap alternative to traditional wildlife monitoring techniques, could be safer to implement in fire/smoke affected areas, and provide much-needed data in near real-time. Moreover, acoustic indices are representative of entire ecological communities, rather than just a select few species that are indicators of environmental change. Just as with standard measures of species diversity, these proxies have the potential to be up-scaled over large landscape areas using spatial imagery (Boelman et al 2007, Pekin et al 2012), which is becoming increasingly available and affordable for tropical countries. Information on the spatial and temporal extent of acoustic disruption, and how well this corresponds to observed changes in ecological communities, will be fundamental to understanding the full environmental implications of future air pollution crises. Acknowledgments This research was supported by the Wildlife Reserves Singapore Conservation Fund and National Parks Board Singapore. We thank Yifei Chung and Joey Gan for field assistance, Yifei Chung and Ming Kai Tan for discussions, and James Cugley for help in devising a protocol for the processing of acoustic data. BPYH Lee was supported by a Ministry of National Development EDGE Scholarship. ORCID iDS Benjamin P Y-H Lee 9952-1011 https://orcid.org/0000-0002- 6 https://orcid.org/0000-0003-0767Zoe G Davies 1467 https://orcid.org/0000-0003Matthew J Struebig 2058-8502 References Boelman N T, Asner G P, Hart P J and Martin R E 2007 Multi-trophic invasion resistance in Hawaii: bioacoustics, field surveys, and airborne remote sensing Ecol. Appl. 17 2137–44 Chan F 2015 Indonesia set to evacuate kids in worst haze hit areas The Straits Times (www.straitstimes.com/asia/se-asia/ indonesia-set-to-evacuate-kids-in-worst-haze-hit-areas) Cheyne S M 2008 Effects of meteorology, astronomical variables, location and human disturbance on the singing apes: hylobates albibarbis Am. J. Primatol. 70 386–92 Cleary D F R, Fauvelot C, Genner M J, Menken S B J and Mooers A 2006 Parallel responses of species and genetic diversity to El Nino ̃ southern oscillation-induced environmental destruction Ecol. Lett. 9 304–10 Depraetere M, Pavoine S, Jiguet F, Gasc A, Duvail S and Sueur J 2012 Monitoring animal diversity using acoustic indices: implementation in a temperate woodland Ecol. Ind. 13 46–54 Doherty R M et al 2009 Current and future-climate and air pollution-mediated impacts on human health Env. Health 8 S1–8 Farina A, Pieretti N and Piccioli L 2011 The soundscape methodology for long-term bird monitoring: a Mediterranean Europe case-study Ecol. Inform. 6 354–63 Farina A and Pieretti N 2014 Sonic environment and vegetation structure: a methodological approach for a soundscape analysis of a Mediterranean maqui Ecol. Inform. 21 120–32 Fuller S, Axel A C, Tucker D and Gage S H 2015 Connecting soundscape to landscape: which acoustic index best describes landscape configuration? Ecol. Ind. 58 207–15 Gage S H, Napoletano B M and Cooper M C 2001 Assessment of ecosystem biodiversity by acoustic diversity indices J. Acoust. Soc. Am. 109 2430 Gaveau D L et al 2014 Major atmospheric emissions from peat fires in Southeast Asia during non-drought years: evidence from the 2013 Sumatran fires Sci. Rep. 4 06112 Gorissen L, Snoeijs T, Van Duyse E and Eens M 2005 Heavy metal pollution affects dawn singing behaviour in a small passerine bird Oecologia 145 504–9 Hardy C C, Ottmar R D, Peterson J L, Core J E and Seamon P 2001 Smoke management guide for prescribed and wildland fire (www.fs.fed.us/pnw/pubs/journals/pnw_2001_ottmar001.pdf) Harrison M E, Boonman A, Cheyne S M, Husson S J, Marchant N C and Struebig M J 2012 Biodiversity monitoring protocols for REDD+: can a one-size-fits-all approach really work? Trop. Conserv. Sci. 5 1–11 Harrison M E, Capilla B R, Thornton S A, Cattau M E and Page S E 2016 Impacts of the 2015 fire season on peat-swamp forest biodiversity in Indonesian Borneo Peatlands in Harmony—Agriculture, Industry and Nature: Oral Presentations: Proc. 15th Int. Peat Congress (Sarawak, Malaysia, 15–19 August 2016) (Kuching: International Peat Society and Malaysian Peat Society) pp 713–17 Huijnen V et al 2016 Fire carbon emissions over maritime Southeast Asia in 2015 largest since 1997 Sci. Rep. 6 26886 Isaksson C 2010 Pollution and its impact on wild animals: a meta-analysis on oxidative stress EcoHealth 7 342–50 Ji W, Joo W, Napoletano B, Qi J, Gage S H and Biswas S 2007 Soundscape characteristics of an environment a new ecological indicator of ecosystem health Wetland and Water Resource Modeling and Assessment: A Watershed Perspective 1st edn ed W Ji (Boca Raton, FL: CRC Press) pp 201–11 Environ. Res. Lett. 12 (2017) 094022 Kimura R, Traber L D, Herndon D N, Linares H A, Lubbesmeyer H J and Traber D L 1988 Increasing duration of smoke exposure induces more severe lung injury in sheep J. Appl. Physiol. 64 1107–13 Kunii O et al 2002 The 1997 haze disaster in Indonesia its air quality and health effects Arch. Environ. Health 5 16–22 McNamara J M, Barta Z, Klaassen M and Bauer S 2011 Cues and the optimal timing of activities under environmental changes Ecol. Lett. 14 1183–90 National Environment Agency 2016 Historical PSI readings (www.haze.gov.sg/haze-updates/historical-psi-readings) Pekin B K, Jung J, Villanueva-Rivera L J, Pijanowski B C and Ahumada J A 2012 Modeling acoustic diversity using soundscape recordings and LIDAR-derived metrics of vertical forest structure in a neotropical rainforest Landsc. Ecol. 27 1513–22 Pieretti N, Farina A and Morri D 2011 A new methodology to infer the singing activity of an avian community: the acoustic complexity index (ACI) Ecol. Ind. 11 868–73 Pieretti N and Farina A 2013 Application of a recently introduced index for acoustic complexity to an avian soundscape with traffic noise J. Acoust. Soc. Am. 134 891–900 Pijanowski B C, Farina A, Gage S H, Dumyahn S L and Krause B L 2011 What is soundscape ecology? An introduction and overview of an emerging new science Landsc. Ecol. 26 213–32 Pijanowski B C et al 2011 Soundscape ecology: the science of sound in the landscape Bioscience 61 203–16 Pye J D 1971 Bats and fog Nature 229 572–4 Robbins C S 1981 Bird activity levels related to weather Estimating Numbers of Terrestrial Birds 1st edn ed C J Ralph and J M Scott (Chicago, IL: Cooper Ornithological Society) pp 301–10 7 Robson C 2008 A Field Guide to The Birds of Southeast Asia (London: New Holland Publishers) Sueur J, Farina A, Gasc A, Pieretti N and Pavoine S 2014 Acoustic indices for biodiversity assessment and landscape investigation Acta Acust. United Acust. 100 772–81 Sueur J, Pavoine S, Hamerlynck O and Duvail S 2008 Rapid acoustic survey for biodiversity appraisal PloS One 3 e4065 Tacconi L 2016 Preventing fires and haze in Southeast Asia Nat. Clim. Change 6 640–3 Towsey M, Wimmer J, Williamson I and Roe P 2014 The use of acoustic indices to determine avian species richness in audio-recordings of the environment Ecol. Inform. 21 110–9 Villanueva-Rivera L J, Pijanowski B C, Doucette J and Pekin B 2011 A primer of acoustic analysis for landscape ecologists Landsc. Ecol. 26 1233–46 Villanueva-Rivera L J and Pijanowski B C 2013 R package soundecology World Health Organization 2016 WHO global urban ambient air pollution database (www.who.int/phe/health_topics/ outdoorair/databases/cities/en/) World Bank 2016 The cost of fire: an economic analysis of Indonesia’s 2015 fire crisis (http://documents.worldbank. org/curated/en/776101467990969768/pdf/103668-BRICost-of-Fires-Knowledge-Note-PUBLIC-ADD-NEWSERIES-Indonesia-Sustainable-Landscapes-KnowledgeNote.pdf) Zvereva E and Kozlov M 2010 Responses of terrestrial arthropods to air pollution: a meta-analysis Environ. Sci. Pollut. Res. Int. 17 297–311