# “El Niño” events increase seasonal wildfire danger in Ecuador

Abstract

Weather is a key driver of catastrophic wildfire seasons. On a global scale, different phases of the El Niño-Southern Oscillation, “El Niño” or “La Niña”, modulate weather and therefore wildfire activity. Ecuador is one the most biodiverse countries in the world and wildfires produce severe impacts to its ecosystems. This investigation explores for the first time the relationship between wildfire weather and El Niño-Southern Oscillation in the tropical Andes. Wildfire weather variability has been quantified using seasonal metrics of the McArthur Forest Fire Danger Index. The calculation of this index requires weather station data that were only available since 1997. Therefore, the availability of data was extended using the Twentieth Century Reanalysis Project to cover the period 1963-2010. Bushfire weather in the seasons April-May-June and July-August-September show significant correlations with the NINO3.4 index (r=0.37,p=0.003, n=47 and r=0.4, p=0.001, n=47 respectively). Additionally, a Chi-squared test showed that every extreme wildfire season in this region was linked to an ‘El Niño’ event. Our results demonstrate that ‘El Niño” events bring hot and dry conditions conducive to wildfires in the antecedent—April-May-June—and simultaneous—July-August-September—wildfire seasons in Ecuador.

# Introduction

Wildfires are a global phenomenon that produces severe impacts (Luke 1978, Conard 1997, Morton 2003). Catastrophic wildfire events claim lives (Haynes 2008, Haynes 2010), degrade the environment (Shakesby 1993, Stephens 2004, Lane 2010) and destroy building infrastructure (Morton 2003, McAneney 2009, Crompton 2010). Wildfires also contribute to biodiversity loss (Kodandapani 2008, Pastro 2011) and affect the climate (Conard 1997). They are part of the natural evolution of vegetation ecosystems (HaynesBradstock 2012). However, human actions amplify the impact over people, nature and infrastructure (Pausas 2011). In fact, wildfires have a strong anthropogenic ignition pattern. In addition, population growth in rural areas increase human vulnerability to this hazard. The threat of wildfires have encouraged extensive research on the behavior of its precursors.

Weather is one of the key components that contributes to wildfire activity (Powell 1983, Williams 1999, Mills 2005, Haines 1988, Mccaw 2007, Dowdy 2012, Engel 2013). Wildfire weather refers to the meteorological influence over wildfires over a period of time —usually within days—. Common variables associated with wildfire weather are temperature, relative humidity, wind speed and precipitation (BoM 2009). Each variable exerts influence in the liberation of water content in vegetation. This process—Evapotranspiration—controls the availability of fuels to be burnt. Most of the research conducted on wildfires and climate focuses on this daily approach (Roads 2005). Yet, the study of wildfire predictability on a seasonal scale is also important. A key factor to explore this perspective is El Niño Southern-Oscillation (ENSO) (Swetnam 1990, Simard 1985, Kitzberger 2002, Beckage 2003, Wooster 2012).

ENSO produces negative rainfall anomalies in some wildfire-prone regions (Swetnam 1990, Simard 1985, Kitzberger 2002, Beckage 2003, Wooster 2012). The intensity of these anomalies depends on the regional influence of ENSO over climate. These rainfall deficits occur over the months before—and during–the wildfire season. ENSO’s influence may even increase by the superposition of the Inter-decadal Pacific Oscillation (IPO) (Verdon 2004). ENSO also has an influence over ignition patterns. In some regions it fosters the occurrence of lightning strikes (Beckage 2003). These facts encourage the investigation of new studies of ENSO and wildfire weather around the world.

The west coast of South America is one of the regions most affected by ENSO. The coastal areas of Peru and Ecuador experience massive rainfall episodes during strong El Niño events (Rodbell 1999). Yet, there has not been enough research on its influence over other natural hazards. Wildfires occur all over South America (Manta 2008). However, studies about wildfire weather in the northern Andean region are scarce. This study aims to describe the seasonal variability of wildfire weather in Ecuador. Also, we investigate if there is a link between this type of weather and ENSO in this region. This is important because wildfires occur every year in this country producing severe impacts (Ministerio del Ambiente 2013, Secretaria de Ambiente 2013). A better understanding of wildfire weather variability and its precursors can provide useful information. These insights could be of special interest to emergency response services and planning authorities.

This article comprises six sections. In the first section we provide background and justify the importance of this study. Section two describes the Ecuadorian Andean region as our study area. Section three presents details about the datasets employed in this research. Section four explains the different methods applied. Section five shows the results of our investigation. Finally, section six presents a discussion of the results and concluding remarks.

## Study area

### Geography

Ecuador is located in the north-west side of South-America (see Figure 1). The country has a diverse geography. It comprises a continental territory between the latitudes 01o 28’ N and 05o 02’ S and longitudes 75o 11’ W and 81o 04’ W (Instituto Geográfico Militar (IGM) 2013). The Galapagos Islands are also part of the country. These islands are approximately 1000 Km from the continent. The national territory has a total area of 256.370 Km2 (Instituto Geográfico Militar (IGM) 2013). The Andes mountains divide the continental Ecuador into three areas: the Coast, Andean and Amazon regions. Each of these regions has diverse geographical features, climate, and ecosystems. The highlands are the most wildfire-prone region because of its specific natural characteristics.

The Andes mountains are the most important geographical feature of Ecuador—and South-America— (Insel 2010). This mountain range crosses Ecuador from North to South. The Ecuadorian Andes comprises two flanks of mountains: the eastern and western ’Cordilleras’. Between these two flanks—that include active volcanoes—there are several inter-Andean valleys. Many of these Andean plateaus are over 3000 m above sea level. The steep elevations have an important effect on local climate. The lapse rate create the conditions for different Bio-climatic zones to exit. This enhances biodiversity. Yet, the original vegetation of the Andean valleys has almost disappeared. Most of the endemic vegetation was replaced by the Australian specie Eucalyptus globulus around 1860 (Ministerio de Ambiente del Ecuador 2012). Other introduced species are the Pinus radiata (from California) and the Pinus patula (from Mexico). This occurred for commercial reasons (Anchaluisa 2013) with great ecological impact. Further, these species are prone to wildfires in the Andean dry season.

### Climate

The tropical Andes climate has several drivers. The Amazon forest, oceanic currents, and the topography are among the main drivers (Martínez 2011). The Amazon forest—through Evapo-transpiration—produces a great amount of water vapor. This water precipitates along the western Cordillera because of an orographic precipitation process. This precipitation is also produced in the eastern Cordillera. Air masses advected from the Pacific Ocean produce this effect. The influence of the Pacific Ocean in the climate variability is important. The Pacific Decadal Oscillation (PDO) is the main driver of decadal climate variability in this region (Martínez 2011). On the other hand, the inter-annual variability depends on the El Niño Southern Oscillation (ENSO) (Vuille 2000, Villacis 2003, Martínez 2011). These influences create several climate regions within the Andes.

Most of the Ecuadorian Andes—between 1500 to 3000 m—have a “Semi-Humid Mesothermal” climate (Pourrut 1995). Mean temperatures in this region range from 8 oc to 20 oc. Maximum temperatures span from 22 oC and 30 oC. Minimum temperatures range from -4 oc to 5 oC. Air masses from the Pacific Ocean and the Amazon region create three seasons (Pourrut 1995). Two wet seasons during the periods February-May and October-November. A dry season spans from June to September. Total annual precipitation ranges from 800 mm to 1500 mm along the Andes. Relative humidity varies from 65 to 85 percent.

### Fire

Wildfires mostly occur in three geographical regions of Ecuador. The Andean region has large extensions of introduced Eucalyptus forests (Anchaluisa 2013). Thus, the availability of fuels makes it the most wildfire-prone region. The fire season in the Ecuadorian Andes spans from July to November (Secretaria de Ambiente 2013). Yet, the most critical months are July, August and September (Estacio 2012, Secretaria de Ambiente 2013). The coast and the Galapagos Islands experience wildfire episodes from January to May (Secretaria de Ambiente 2013).

Wildfires produce severe impacts in Ecuador. Yet, information about these impacts is scarce and disperse. Wildfires produced 21.570 Ha of area burnt in Ecuador in 2012 (Ministerio del Ambiente 2013). The Metropolitan District of Quito (MDQ) is particularly vulnerable to wildfires. The vulnerability is high because 15.4 percent of the national population lives here. Additionally, vegetation covers 60.46 percent of its territory (Secretaria de Ambiente 2013). Wildfires destroyed 2700 Ha of the MDQ vegetation in 2009 (Estacio 2012). In 2012, the area burnt by wildfires increased to 4.882,16 Ha (Secretaria de Ambiente 2013). This wildfire season cost the MDQ more than 50 million dollars (Secretaria de Ambiente 2013).

Unfortunately, wildfires in Ecuador are anthropogenic in nature (Ministerio del Ambiente 2013). In fact, the MDQ reported that humans cause 95 percent of wildfires in the city (Secretaria de Ambiente 2013). Burning forests to convert them to agricultural use is the main cause of wildfires in the country (Ministerio del Ambiente 2013, Secretaria de Ambiente 2013, Rodas 2015). Forcing land use changes for urbanization purposes is another reason people start wildfires (Secretaria de Ambiente 2013). Rural communities burn agricultural waste, which is a traditional practice that fosters wildfire occurrences (Estacio 2012). Rural people also burn regular waste because of the lack of waste management services. (Secretaria de Ambiente 2013). Finally, wildfires occur because of the damaging acts of arsonists (Ministerio del Ambiente 2013, Secretaria de Ambiente 2013).

# Data and methods

## Data

### Weather station

This study required weather station (WS) data from the Andean region in Ecuador. A limitation with this dataset was that only six stations had enough records for this study. The selected period of analysis was 1997-2012. This period had almost complete daily information. The variables for this research are air temperature (Tmax), relative humidity (RH), wind speed (W) and precipitation (P). Each of the stations provided data on an eight-hourly basis. The hour selected for the computations was 1 p.m. local time (GMT 18:00). Table 1 describes each of the stations.

### Reanalysis

This research also used reanalysis data from the Twentieth Century Reanalysis Project (20CR) (Compo 2011). The 20CR relies on a data assimilation algorithm of surface pressure values. Sea-surface temperature and sea ice distributions are boundary conditions in this model. The 20CR yields an ensemble mean of realizations. This produces the best possible representation of the state of the atmosphere. This dataset yields six-hourly data on 24 vertical levels. The 1000 hPa-level data is the first vertical level. Thus, it is the most likely to represent surface values. The hour selected for calculations with this data was 4 p.m. local time (GMT 06:00). This is the time closest to the hour adopted for most wildfire potential computations (3 p.m.) with the Australian Forest Fire Danger Index (FFDI)—to be explained next—. The 20CR has data with a resolution of 2o x 2o. Thus, the extracted data is a grid in the Ecuadorian Andean region.

This dataset presents limitations common to all reanalysis products. There is greater uncertainty in regions were there are less observations. This is a particular problem in Africa, some regions in Asia and South America. There is less reliability in the first half of the century due to lack of observations. Although there have been studies that reflect this data represents well climate variability during this period (Ashcroft 2014). Yet, studies for South America are not available. Additionally, some variables are better represented than others.

### ENSO Indices

A common ENSO metric is the Southern Oscillation Index (SOI). The SOI represents the difference in normalized surface pressure between Taihiti and Darwin, Australia. We adopted the SOI’s calculation method proposed by the National Oceanic and Atmospheric Administration (NOAA 2015). The index reflects the changes in air-mass fluxes in the tropical Pacific Ocean. These fluxes, are a manifestation of changes in sea surface temperature (SST). The SOI has an open ended scale with positive and negative values. A negative SOI value represents above-average warm SST conditions in the eastern tropical Pacific Ocean. This represents an ’El Niño’ phase of the phenomenon. El Niño brings wet conditions to the coasts of southern Ecuador and northern Peru. The opposite pattern occurs with a positive SOI value, representing a ’La Niña’ phase.

Another ENSO index is the Oceanic Niño Index (ONI) (NOAA 2015a). The ONI measures monthly SST values in the tropical Pacific Ocean. The measurements take place in the region from 5 oN to 5 oS latitude and 120 oW to 170 oW longitude. This is the ’Niño 3.4’ region. The index requires the calculation of 3-month running means and its climatology. The ONI is the difference between these two values. An ’El Niño’ condition exits if there is at least an anomaly of +0.5 oC in this computation. An anomaly of -0.5 oC indicates a ’La Niña’ phase.

## Methods

### Data Pre-processing

The computation of wildfire potential usually requires continuous daily data. Thus, only six stations in the Ecuadorian Andean region were suitable for this study. The selected stations are representative of the ecosystems that are prone to wildfire. The available records in these stations spanned the period 1997-2012. Yet, the period of analysis is short. Thus, reanalysis data was bias-corrected to extend the availability of records.

The bias-correction was necessary for several reasons. First, 1000 hPa reanalysis data does not represent surface level values. Second, there was a 3-hour difference between the time of the two dataset’s values. Third, reanalysis data in the highlands may differ in quality to sea-level outputs. Thus, we plotted probability density functions (PDFs) of the reanalysis and weather station datasets. This allowed to determine the differences between their data distributions. Applying linear scaling techniques was a suitable approach to bias-correct the required variables.

The aim of the study was to estimate wildfire potential for the entire Ecuadorian Andean region. Thus, a calculation of daily averages for each wildfire weather variable yielded regional values. With the reanalysis dataset the procedure is straightforward. These data is spatially and temporally continuous. However, the weather stations data had some missing days. Therefore, calculating averages over daily anomalies made the results less sensitive to missing data.

### Forest Fire Danger Index (FFDI) calculation

Ecuador currently does not have any scientifically designed wildfire potential index. Therefore, an alternative approach was to use an index suitable for this region. The FFDI is an Australian wildfire potential index (Noble 1980). This index uses daily values of temperature, relative humidity and wind speed. The index also incorporates a drought factor. This factor represents the availability of fuels. This is measured indirectly using soil moisture indexes such us the Keetch–Byram drought index (Keetch JJ 1968). The FFDI equation is expressed as follows:

$\label{eqn:FFDI} FFDI = 2e^{-0.45+0.987*Ln(DF)-0.0345RH+0.0338T+0.0234v}$

The variables represented are drought factor (DF), relative humidity (RH), temperature (T) and wind speed (v). This empirical index is valid for the southeastern forest ecosystem of Australia. In this wildfire-prone ecosystem, the eucalyptus is the dominant vegetation. On the other side of the world, Latin-America—including Ecuador—experienced a massive introduction of this specie (Anchaluisa 2013). Thus, is reasonable to expect that this index may represent the Ecuadorian wildfire potential.

The application of this index with weather station and reanalysis data yielded daily values. With weather station data the calculation comprised the period 1997-2012. The bias-corrected reanalysis data allowed to extend the method for the period 1963-2012. For each year, a sum of the daily FFDI values provided a seasonal metric of wildfire potential. This cumulative FFDI corresponds to the wildfire season July-August-September (JAS). Finally, we calculated the 90th percentile on the seasonal FFDI time series. This determined which were extreme wildfire seasons during the period 1963-2012.

### Investigation of ENSO influence over wildfire seasons

We used several approaches to investigate ENSO’s link over wildfire seasons. Firstly, it was important to know the influence of ENSO in the months leading to the wildfire season JAS. Thus, we applied linear correlations of ENSO’s indexes (seasonal averages) with FFDI values. These averages comprised three-months seasons starting from October-November-December (OND) of the previous year. The correlation also included the SOI and Niño 3.4 values for the concurrent JAS season.

Secondly, we explored the relationship of each wildfire variable with ENSO. This comprised exploring the link of concurrent FFDI variables (JAS) with antecedent (AMJ) and concurrent ENSO indexes values—SOI and Niño 3.4—. The computation of linear correlations allowed to establish these relationships.

A third approach compared the extreme wildfire seasons with ENSO’s categories of intensity. We used the categories proposed by the Australian Bureau of Meteorology. Finally, we expanded this analysis using a $$\chi^2$$ test. In this test, we used categories of seasonal FFDI and ENSO strength. The categories for seasonal FFDI were extreme (FFDI >90th percentile), medium (50th percentile<FFDI<=90th percentile), and low (FFD<=50 percentile). The numbers of ENSO years were added within three categories: very strong El Niño, Weak-Moderate-Strong El Niño, and La Niña. The strength of ENSO years corresponds to the Australian Bureau of Meteorology classification available at http://www.bom.gov.au/climate/enso/enlist/.

# Results

Figure 2 shows PDFs for each wildfire weather variable. These results represent the deviations between original WS and 20CR data. All 20CR PDFs have deviations from the WS data. The daily deviations from the mean including all seasons are 8.9 oC (Tmax), 11 (RH), 1.6 m/s (W) and 4.5 mm (P).

Figure 3 presents time series of the same variables. The values correspond to seasonal averages for Tmax, RH and W during JAS. P corresponds to a seasonal sum. The figure shows Pearson’s linear correlation coefficients between 20CR and WS data for 16 seasons. The results were statistically significant only for W (r=0.52,p=0.04).

The bias-correction process comprised the application of linear scaling techniques for every wildfire weather variable. However, the bias-correction of all variables did not improve the seasonal FFDI calculation (not shown). Only the temperature bias-correction significantly improved the FFDI results. Thus, for this study we only show the bias-correction results for this variable (see Figure 4).

Figure 4 shows PDFs of original and bias-corrected Tmax data. The application of additive linear scaling successfully centers the mean of the 20CR. Additionally, Figure 5 presents original and bias-corrected results of the seasonal time series for Tmax. After the bias-correction, 20CR maximum temperature data has a comparable magnitude with WS data. The linear scaling did not affect the correlation coefficients between the two types of data.

Figure 6 presents the results of the FFDI calculation with original and bias-corrected data. The results show PDFs of WS and bias-corrected 20CR data. The bias-correction increases the similarity between the distribution of data.

Figure 7 presents seasonal (JAS) time series of FFDI. The figure displays three time series. FFDI calculated with WS data comprises the period 1997-2012. Two more time series show FFDI calculated with reanalysis data. The calculation used original (20CR) and bias-corrected (20CR(bc)) reanalysis data. These time series span the period 1963-2012. The FFDI time series calculated with original reanalysis data over-estimates seasonal wildfire potential. The overestimation is—on average—of 65 units during the period 1997-2012. These data—WS and 20CR—have a moderate linear correlation (r=0.53,p=0.04). FFDI calculated with bias-corrected reanalysis data under-estimates the FFDI computed with observed values. The magnitude of the underestimation is less than the initial over-estimation (47 units). The bias-correction does not significantly affect the linear correlation between WS and 20CR(bc) data (r=0.53,p=0.05). These results encouraged us to extend the FFDI calculation for the period 1963-2012.

Having an extended FFDI dataset allowed us to investigate the influence of ENSO over seasonal wildfire variability in Ecuador. Table 2 shows linear correlation’s ’R’ coefficients of seasonal ENSO metrics with FFDI. The table shows how the evolution of ENSO in previous—and the concurrent season—is associated with wildfire potential. Previous seasons start one year earlier on October-November-December (OND), and continue in three-months intervals including the wildfire season July-August-September (JAS). These correlations show that the influence of ENSO is associated to wildfire potential in the preceding (AMJ) and concurrent (JAS) season. ENSO does not seem to have a direct relationship further back in time. This is valid using the SOI and Niño 3.4 indices.

These results encouraged further analysis into the concurrent (JAS) and antecedent (AMJ) seasons. Table 3 presents the results of linear correlations applied to every wildfire variable in the concurrent season, with ENSO indices in the antecedent and concurrent season. The relationships among FFDI variables are physically consistent. The R coefficients for these variables with ENSO indices are moderate are statistically significant at the 5 percent level. This holds valid for all variables with the exception of precipitation.

This investigation also explored the relationship between ENSO and extreme wildfire seasons. Our definition of extreme wildfire seasons filtered years with cumulative FFDI above the 90th percentile for the period 1963-212. Table 4 shows that the most extreme seasons occurred during the years 1982, 1991, 1997, 2009 and 2012. Each season was compared with historical categorizations of strength of ENSO events. With the exception of the year 2012, extreme wildfire weather occurred during El Niño events.

Table 4 expands on the categorical analysis of seasonal wildfire potential and ENSO strength with a $$\chi^2$$ test. The table shows that most very strong El Niño events—two out of four years—have a connection with extreme wildfire seasons. Within the range of Weak-Moderate-Strong El Niño events, most years—13 out 22—occur during low wildfire potential—seasonal FFDI<=50th percentile—. However, two years belong to the extreme wildfire weather category. Finally, most La Niña years occur during medium and low seasonal FFDI—23 out 24 years—.

# Discussion

This study has investigated the variability of seasonal wildfire potential in Ecuador. The analysis included the examination of the relationship of wildfire weather and ENSO phases. Using a seasonal metric of wildfire potential—cumulative FFDI during the season July-August-September—, this study detected that five extreme seasons occurred during the period 1963-2012. An extreme season was defined as any year with a seasonal FFDI greater than the 90th percentile value. Most of these seasons—four out of five—occurred during diverse categories of El Niño strength.

The El Niño events in the Tropical Andes produce a below-average precipitation pattern (Garreaud 2009). Yet, this response is not uniform in the Ecuadorian Andes (Villacis 2003). No relationship exists between ENSO and precipitation during the wet season (October-May) (Villacis 2003). Yet, there is a consistent trend with below-average precipitation during the dry season (June-August) (Villacis 2003). This may occur because precipitation has some Atlantic influence as well (Vuille 2000). Additionally, ENSO can largely explain the inter-annual variability of temperature (Vuille 2000). El Niño events are associated with higher temperatures in most of the Ecuadorian Andes (Vuille 2000). However, these studies relied on a relatively short period—approximately 30 years—of weather station data. Moreover, they include only some wildfire variables —temperature and precipitation—. In spite of this, they show a consistent agreement on the effect of El Niño over environmental conditions conducive to wildfires. This agrees with our overall results. Nevertheless, this should be carefully analyzed considering the use of reanalysis data in this study.

The general limitations of reanalysis data were discussed in the data and methods section. The results shown in Table 3 demonstrate that the precipitation produced by ENSO’s teleconnections is not well represented by the 20CR data. This variable should be strongly influenced by ENSO. Yet, the linear correlations do not show any link. This is an unexpected result especially for the dry season. The rest of the wildfire variables do show a statistically significant relationship with ENSO. In spite of the acknowledged biases in the data, the relationship is physically consistent. Yet, the coarse resolution of the reanalysis data brings uncertainty to the results. The quantification of these uncertainties is beyond the scope of this study. This could be accomplished examining each 20CR realization. An alternative approach could be the use of other reanalysis products. Although, the limitations may be similar. Another source of uncertainty is the selected wildfire index.

The FFDI is an empirically designed index calibrated for southeast Australia. To our knowledge, this index has not being tested outside of Australia. The ecosystems of southeast Australia and the Ecuadorian Andes share the wildfire-prone Eucaliptus globulus vegetation. Yet, the particular characteristics of the Tropical Andean ecosystems may produce imprecise wildfire potential results. The computation of alternative indexes could provide a source of validation. The Canadian Forest Fire Weather Index (FWI) (Van Wagner 1987) has been used in several countries with positive results (Dowdy 2009). Although computationally more complex to calculate, it could be interesting to test how well it represents fire weather in Ecuador compared to the FFDI. A future endeavor might be as well the design of an Ecuadorian Fire Weather Index System. Finally, another approach to strengthen this analysis would be the investigation of ENSO’influence over wildfire activity rather than fire weather.

# Concluding remarks

Roboust long-term wildfire activity records do not exist in Ecuador. The MDQ has started to produce more detailed information about wildfire impacts based on satellite data starting from the year 2009 (Estacio 2012). However, the rest of the country lacks this information that would be valuable to evaluate ENSO’s influence over wildfires. For example, basic fire statistics such as the number of fire occurrences and area burnt per year would provide additional metrics to evaluate the impact of El Niño events over wildfires in Ecuador. Further investigation of ENSO’s influence over wildfires in Ecuador could enhance our capacity to anticipate these events.

Several studies agree that there is potential to use ENSO to forecast wildfire seasons (Simard 1985, Nicholls 2007, Wooster 2012, Harris 2013). In fact, a predictive index has already been proposed for the Norhtern Patagonia region in Argentina (Kitzberger 2002). However, as acknowledged by the author, the use of statistical-empirical methods to investigate seasonal wildfire predictability has limitations. For example, the assumption of stable ENSO phase’s relationships over a long period of time. Another important assumption to consider is the fixed timing of ENSO arrival. In spite of this, simple indices based on the relationship of ENSO and wildfire weather—or activity—can provide valuable information for decision makers. This information will facilitate timely decisions that can minimize economic impacts, protect the environment, and potentially save lives.

Figure 2: Probability density functions (PDFs) of original data (without bias-correction). The PDFs show data for four wildfire variables: Maximum Temperature (Tmax), Relative Humidity (RH), Wind Speed (W) and Precipitation (P). The period of analysis is 1997-2012.

Figure 3: Seasonal July-August-September (JAS) time series of wildfire weather variables in Ecuador for the period 1997-2012. Original Tweentieth Century Reanalysis (20CR) and weather station (WS) data correspond to regional averages in the Ecuadorean Andean region. Pearson’s linear correlation coefficients for 16 seasons are: a) Maximum temperature (Tmax)(r=0.36,p=0.17), b) Relative humidity (RH) (r=0.14,p=0.59), c) Wind speed (W) (r=0.52,p=0.04), d) Precipitation (P) (r=0.33, p=0.21).

Figure 4: Probability density functions (PDFs) of original and bias-corrected maximum temperature (Tmax) in the Ecuadorian Andean region. The figure shows weather station (WS), Tweentieth Century Reanalysis (20CR) and 20CR data with bias-correction (20CR(bc)). These data comprises the seasons July-August-September (JAS) for the period 1997-2012.

Figure 5: Time series of seasonal averages of original and bias-corrected maximum temperature in the Ecuadorian Andean region. The figure shows weather station (WS) and Tweentieh Century Reanalysis (20CR) data. These data comprises the seasons July-August-September (JAS) for the period 1997-2012. Panel a) and b) linear correlations statistics are r=0.36,p=0.17,n=16. The bias-correction technique applied was additive linear scaling.

Figure 6: Probability density functions (PDFs) of original and bias-corrected Forest Fire Danger Index (FFDI) data in the Ecuadorian Andean region. The figure shows FFDI calculated with weather station (WS), Tweentieth Century Reanalysis (20CR) and 20CR data with bias-correction (20CR(bc)). These data comprises the seasons July-August-September (JAS) for the period 1997-2012.

Figure 7: Time series of Forest Fire Danger Index (FFDI) seasonal values in the Ecuadorian Andean region. The figure shows FFDI calculated with weather station (WS) and Tweentieh Century Reanalysis (20CR) data. These data comprises the seasons July-August-September (JAS) for the period 1997-2012 for WS, and 1963-2012 for 20CR and 20CR(bc) data. Original WS and 20CR linear correlation statistics for the period 1997-2012 are r=0.53,p=0.04,n=16. These statistics for original WS and bias corrected 20CR(bc) data are r=0.51,p=0.05,n=16.

Table 1: Selected weather stations in the Ecuadorian Andean region.

Table 2: Pearson’s R coefficients between Southern Oscillation Index (SOI), Niño 3.4 seasonal averages and bias-corrected Forest Fire Danger Index (FFDI) values. ENSO indices correspond to the seasons October-November-December (OND)—previous year—, January-February-March (JFM), April-May-June (AMJ) and July-August-September (JAS). FFDI values are the cumulative sum of daily values during the season JAS. The linear correlations correspond to the period 1963-2012.

Table 3: R linear coefficients of concurrent (JAS) wildfire weather variables with antecedent (AMJ) and concurrent ENSO indexes. Wildfire variables correspond to original 20CR data with the exception of maximum temperature (Tmax) which was bias-corrected. The period of analysis is 1964-2012.

Table 4: Extreme wildfire potential in Ecuador and associated ENSO types. A calculation of the 90th percentile of seasonal (JAS) FFDI from 1963 to 2012 yielded five extreme seasons. A strong El Niño event is common is extreme seasons. The classification of ENSO types belongs to the Australian Bureau of Meteorology available at http://www.bom.gov.au/climate/enso/enlist/

Table 5: Chi-squared test for categories of wildfire seasons—seasonal FFDI—and ENSO events in Ecuador for the period 1963-2012 (50 season in total).

### References

1. R Luke, A McArthur. Bushfires in Australia. (1978).

2. Susan G Conard, Galina A Ivanova. Conard, Ivanova - 1997 - Wildfire in Russian Boreal ForestsâPotential Impacts of Fire Regime Characteristics on Emissions and Global Carbon Balance Estimates.pdf. 98 (1997).

3. Douglas C Morton, Megan E Roessing, Ann E Camp, Mary L Tyrrell. Assessing the environmental, social, and economic impacts of wildfire. GISF Research Paper 001.. 54 (2003).

4. Katharine Haynes, Amalie Tibbits, Lucinda Coates, Gaminda Ganewatta, John Handmer, John Mcaneney, Risk Frontiers, Community Safety, Bushfire Cooperative. 100 years of Australian civilian , bushfire fatalities : exploring the trends in relation to the ’ stay or go policy ’. (2008).

5. Katharine Haynes, John Handmer, John McAneney, Amalie Tibbits, Lucinda Coates. Australian bushfire fatalities 1900-2008: exploring trends in relation to the ‘Prepare, stay and defend or leave early’ policy. Environmental Science & Policy 13, 185–194 Elsevier Ltd, 2010. Link

6. Ra Shakesby, Cda Coelho, Ad Ferreira, Jp Terry, Rpd Walsh. Wildfire Impacts on Soil-Erosion and Hydrology in Wet Mediterranean Forest, Portugal. International Journal of Wildland Fire 3, 95 (1993). Link

7. Scott L. Stephens, Thomas Meixner, Mark Poth, Bruce McGurk, Dale Payne. Prescribed fire, soils, and stream water chemistry in a watershed in the Lake Tahoe Basin, California. International Journal of Wildland Fire 13, 27–35 (2004). Link

8. P.N.J. Lane, P.M. Feikema, C.B. Sherwin, M.C. Peel, a.C. Freebairn. Modelling the long term water yield impact of wildfire and other forest disturbance in Eucalypt forests. Environmental Modelling & Software 25, 467–478 Elsevier Ltd, 2010. Link

9. John McAneney, Keping Chen, Andy Pitman. 100-years of Australian bushfire property losses: is the risk significant and is it increasing?. Journal of environmental management 90, 2819–22 (2009). Link

10. Ryan P. Crompton, K. John McAneney, Keping Chen, Roger a. Pielke, Katharine Haynes. Influence of Location, Population, and Climate on Building Damage and Fatalities due to Australian Bushfire: 1925â2009. Weather, Climate, and Society 2, 300–310 (2010). Link

11. Narendran Kodandapani, Mark a. Cochrane, R. Sukumar. A comparative analysis of spatial, temporal, and ecological characteristics of forest fires in seasonally dry tropical ecosystems in the Western Ghats, India. Forest Ecology and Management 256, 607–617 (2008). Link

12. Louise a. Pastro, Christopher R. Dickman, Mike Letnic. Burning for biodiversity or burning biodiversity? Prescribed burn vs. wildfire impacts on plants, lizards, and mammals. Ecological Applications 21, 3238–3253 (2011). Link

13. R. A. HaynesBradstock, A. M. Gill, R. J.Bushfire1 Williams. Flammable Australia : Fire Regimes, Biodiversity and Ecosystems in a Changing World. 345 CSIRO Publishing, 2012.

14. Juli G. Pausas, Santiago Fernández-Muñoz. Fire regime changes in the Western Mediterranean Basin: from fuel-limited to drought-driven fire regime. Climatic Change 110, 215–226 (2011). Link

15. F. a. Powell. Bushfire Weather. Weatherwise 36, 126–131 (1983). Link

16. Allyson A J Williams, David J Karoly. Extreme fire weather in Australia and the impact of the El Nino-Southern Oscillation. Australian Meteorological Magazine 48, 15–22 (1999).

17. Graham A Mills. A re-examination of the synoptic and mesoscale meteorology of Ash Wednesday 1983. 54, 35–55 (2005).

18. Donald A Haines, Forest Service. A lower atmosphere severity index for wild fires. National Weather Digest 13, 23–27 (1988).

19. L Mccaw, P Marchetti, G Elliott, G Reader. Bushfire weather climatology of the Haines Index in southwestern Australia. 56, 75–80 (2007).

20. Andrew J. Dowdy, Graham a. Mills. Atmospheric and Fuel Moisture Characteristics Associated with Lightning-Attributed Fires. Journal of Applied Meteorology and Climatology 51, 2025–2037 (2012). Link

21. Chermelle B. Engel, Todd P. Lane, Michael J. Reeder, Michael Rezny. The meteorology of Black Saturday. Quarterly Journal of the Royal Meteorological Society 139, 585–599 (2013). Link

22. BoM. Bushfire Weather. (2009). Link

23. J. Roads, F. Fujioka, S. Chen, R. Burgan. Seasonal fire danger forecasts for the USA. International Journal of Wildland Fire 14, 1–18 (2005). Link

24. Thomas Swetnam, Julio Betancourt. Fire-Southern Oscillation Relations in the Southwestern United States. Science 249, 1017–1020 (1990). Link

25. Albert Simard, Donald Haines, William Main. Reltions between El Nino/Southern Oscillation anomalies and wildland fire activity in the United States. Agricultural and Forst Meterology 36, 93–104 (1985). Link

26. Thomas Kitzberger. ENSO as a forewarning tool of regional fire occurrence in northern Patagonia, Argentina. International Journal of Wildland Fire 11, 33–39 (2002). Link

27. Brian Beckage, William J. Platt, Matthew G. Slocum, Bob Panko. Influence of the El NinÌo Southern Oscillation on Fire Regimes in the Florida Everglades. Ecology 84, 3124–3130 (2003). Link

28. M. J. Wooster, G. L W Perry, a. Zoumas. Fire, drought and El Niño relationships on Borneo (Southeast Asia) in the pre-MODIS era (1980-2000). Biogeosciences 9, 317–340 (2012). Link

29. Danielle Verdon, Anthony Kiem, Stewart Franks. Multi-decadal variability of forest fire risk â eastern Australia. 165–171 (2004).

30. D T Rodbell, G Seltzer, D M Anderson, M B Abbott, D B Enfield, J H Newman. A 15,000 year record of ENSO-driven alluviation in Southwestern Ecuador. Science 283, 516–520 (1999).

31. Maria Manta, Patricio Sanhueza. Wildfires in South America. Crisis|Response 6 (2008).

32. MAE Ministerio del Ambiente. Plan de contingencia de incendios forestales. (2013).

33. DMQ Secretaria de Ambiente. Los incendios forestales en el Distrito Metropolitano de Quito evaluación Los incendios forestales en el Distrito Metropolitano de Quito evaluación 2012. 1–57 (2013).

34. Instituto Geográfico Militar (IGM). Contextos Históricos y Políticos Generales. 350 In Atlas Geográfico de la República del Ecuador. (2013).

35. Nadja Insel, Christopher J. Poulsen, Todd a. Ehlers. Influence of the Andes Mountains on South American moisture transport, convection, and precipitation. Climate Dynamics 35, 1477–1492 (2010). Link

36. Ministerio de Ambiente del Ecuador. Sistema de clasificación de los ecosistemas del Ecuador continental. Subsecretaría de Patrimonio Natural (a). 143 (2012). Link

37. Stephanie Anchaluisa. Efecto del fuego sobre la estructura, microclima y funciones ecosistémicas de los bosques introducidos de eucalipto (Eucalyptus globulus,; Mytaceae) en el Distrito Metropolitano de Quito. (2013). Link

38. Rodney Martínez, Daniel Ruiz, Marcos Andrade, Luis Blacutt, Daniel Pabon, Ena Jaimes, Gloria León, Marcos Villacís, Juan Quintana, Edgard Montealegre, Christian Euscátegui. Synthesis of the climate of the tropical Andes. 97–109 In Climate change and biodiversity in the tropical Andes. (2011).

39. Mathias Vuille, Raymond S. Bradley, Frank Keimig. Climate Variability in the Andes of Ecuador and Its Relation to Tropical Pacific and Atlantic Sea Surface Temperature Anomalies. Journal of Climate 13, 2520–2535 American Meteorological Society, 2000. Link

40. M Villacis, J D Taupin, E Heredia, E Palacios, J Palacios. Variabilité climatique dans la sierra équatorienne en relation avec le phénomène ENSO. Hydrology of the Mediterranean and Semiarid Regions (Proceedings of an international symposium held at Montpellier, April 2003). IAHS Publ. no. 278 202–208 (2003).

41. Pierre Pourrut, O Rovere, I Romo, H Villacrés. Clima del Ecuador. 13–26 In El agua en el Ecuador. (1995). Link

42. J Estacio, N Narvaez. Incendios forestales en el Distrito Metropolitano de Quito (DMQ): conocimiento e intervencion publica del reisgo. (2012). Link

43. Melissa Rodas. Análisis del comportamiento de incendios forestales en la cuenca del río Paute mediante sensores remotos. (2015). Link

44. G. P. Compo, J. S. Whitaker, P. D. Sardeshmukh, N. Matsui, R. J. Allan, X. Yin, B. E. Gleason, R. S. Vose, G. Rutledge, P. Bessemoulin, S. Brönnimann, M. Brunet, R. I. Crouthamel, a. N. Grant, P. Y. Groisman, P. D. Jones, M. C. Kruk, a. C. Kruger, G. J. Marshall, M. Maugeri, H. Y. Mok, Ø. Nordli, T. F. Ross, R. M. Trigo, X. L. Wang, S. D. Woodruff, S. J. Worley. The Twentieth Century Reanalysis Project. Quarterly Journal of the Royal Meteorological Society 137, 1–28 (2011). Link

45. Linden Ashcroft, David John Karoly, Joëlle Gergis. Southeastern Australian climate variability 1860-2009: a multivariate analysis. International Journal of Climatology 34, 1928–1944 (2014). Link

46. NOAA. Southern Oscillation Index. (2015). Link

47. NOAA. Equatorial Pacific Sea Surface Temperatures. (2015). Link

48. I. R. Noble, a. M. Gill, G. a. V. Bary. McArthur’s fire-danger meters expressed as equations. Austral Ecology 5, 201–203 (1980). Link

49. Byram GM Keetch JJ. A Drought Index for Forest Fire Control. 32 (1968).

50. R. D. Garreaud. The Andes climate and weather. Advances in Geosciences 22, 3–11 (2009). Link

51. C. E Van Wagner. Development and Structure of the Canadian Forest FireWeather Index System. (1987).

52. Andrew J Dowdy, Graham A Mills, Klara Finkele, William De Groot. Australian fire weather as represented by the McArthur Forest Fire Danger Index and the Canadian Forest Fire Weather Index Australian fire weather as represented by the McArthur Forest Fire Danger Index and the Canadian Forest Fire Weather Index. (2009).

53. Neville Nicholls, Christopher Lucas. Interannual variations of area burnt in Tasmanian bushfires: Relationships with climate and predictability. International Journal of Wildland Fire 16, 540–546 (2007). Link

54. S. Harris, N. Nicholls, N. Tapper. Forecasting fire activity in Victoria , Australia , using antecedent climate variables and ENSO indices. International Journal of Wildland Fire (2013). Link