Technical note: Low meteorological influence found in 2019 Amazonia fires
The sudden increase in Amazon fires early in the 2019 fire season made global headlines. While it has been heavily speculated that the fires were caused by deliberate human ignitions or human-induced landscape changes, there have also been suggestions that meteorological conditions could have played a role. Here, we ask two questions: were the 2019 fires in the Amazon unprecedented in the historical record, and did the meteorological conditions contribute to the increased burning? To answer this, we take advantage of a recently developed modelling framework which optimises a simple fire model against observations of burnt area and whose outputs are described as probability densities. This allowed us to test the probability of the 2019 fire season occurring due to meteorological conditions alone. The observations show that the burnt area was higher than in previous years in regions where there is already substantial deforestation activity in the Amazon. Overall, 11 % of the area recorded the highest early season (June–August) burnt area since the start of our observational record, with areas in Brazil's central arc of deforestation recording the highest ever monthly burnt area in August. However, areas outside of the regions of widespread deforestation show less burnt area than the historical average, and the optimised model shows that this low burnt area would have extended over much of the eastern Amazon region, including in Brazil's central arc of deforestation with high fire occurrence in 2019. We show that there is a 9 % likelihood of the observed August fires being caused by meteorological conditions alone, decreasing to 6 %–7 % along the agricultural–humid forest interface in Brazil's central states and 8 % in Paraguay and Bolivia dry forests. Our results suggest that changes in land use, cover or management are the likely drivers of the substantial increase in the 2019 early fire season burnt area, especially in Brazil. Burnt area for September in the arc of deforestation had a 14 %–26 % probability of being caused by meteorological conditions, potentially coinciding with a shift in fire-related policy from South American governments.
South American fires made global headlines in August 2019, with the largest increase in fire activity seen in nearly 10 years (INPE, 2019; Lizundia-Loiola et al., 2020). Of the roughly 100 000 fires burning by the end of the month, around half were in the Amazon rainforest region (de Oliveira Andrade, 2019; INPE, 2019). While fires in drier savannah regions of South America such as the Cerrado are more common, fires in the rainforest are not a natural occurrence and are rarely ignited without human intervention (Aldersley et al., 2011). As such, fires in humid, tropical regions where the vegetation is not adapted to frequent burning (Kelley, 2014; Staver et al., 2020; Zeppel et al., 2015), have much higher tree mortality rates (Brando et al., 2014; Cochrane and Schulze, 1999; Pellegrini et al., 2017). As a result, an estimated 906 000 ha of the Amazon biome was lost to fires in 2019 (Butler, 2017). The amount of carbon and trace gas emission was also a major concern given the high biomass of the areas being burnt, and smoke from these fires reached cities as far as São Paulo more than 2700 km away (Lovejoy and Nobre, 2019). Usually, small-scale fires in Amazonia are associated with deliberate but localised deforestation, although in dry years, there is more risk of these fires escaping into much larger areas (Aragão et al., 2018). Hence the substantial increase in fires in 2019 sparked much debate about whether the level of burning was unprecedented, whether increased burning was driven by a drier than normal fire season and whether raised levels of direct deforestation played a role (Arruda et al., 2019; Escobar, 2019).
The Amazon has not been the only place with recent unusual and high fire activity, with large-scale fire events worldwide in the last couple of years including in the Arctic, Mediterranean, Australia, the UK and the US. In November 2018 over 80 people were killed in the Camp Fire in Paradise, California, the most destructive in California's history, with the Camp, Woolsey and Carr fires together costing an estimated USD 27 billion in damages (Nauslar et al., 2019). California 2020 fire season is set to be even more damaging (Anon, 2020), though the true extent and costs are still increasing at the time of writing. Hundreds of fires burnt in Siberia and Alaska throughout the 2019 and 2020 Northern Hemisphere summers, releasing over 150 Mt of CO2 in 2019 and 244 Mt in 2020 into the atmosphere (Witze, 2020). Of major concern was the release of large quantities of black carbon with the potential to further accelerate local arctic ice melt (Patel, 2019). The UK saw some burning, including a peatland fire in north-east Sutherland that doubled Scotland's carbon emission for 6 d in May 2019 (Wiltshire et al., 2019). Between September 2019 and February 2020, fires across south-eastern Australian burnt around 18.6 × 106 ha, destroyed over 5900 buildings and killed at least 34 people (Boer et al., 2020; RFS, 2019; Sanderson and Fisher, 2020). Unusual fire events such as these are expected to increase in frequency in the future from both changes in climate and socio-economic pressures on the landscape (Fonseca et al., 2019; Jones et al., 2020). Given the concerns raised and the extent to which much of these fire events captured the attention of the public and press in recent months, in the aftermath, it is important to look at these events objectively. In particular, it is essential to determine whether they were unusual in the context of the historical record and, if so, what might be new and emerging drivers.
There are many ways to assess drivers of historical fire events. Some studies simply correlate individual drivers with burnt area in isolation (Andela et al., 2017; Van Der Werf et al., 2008). However, these do not consider the complex interaction of multiple drivers regarding fire and are therefore unable to go beyond a loose attribution of a particular forcing to fire, which can easily be confused as causations due to the number of drivers. Fire danger indices (FDIs), which can capture simultaneous drivers, are useful for calculating the level of risk of a fire spreading and becoming severe in a particular area (de Groot et al., 2015). FDIs have been adapted to assess recent and future trends in climate regarding fire weather (Burton et al., 2018; Jolly et al., 2015) and attribute increases in fire risk to anthropogenic changes in climate (van Oldenborgh et al., 2020). These metrics thereby provide rapid policy-relevant information for fire management (De Groot et al., 2010; Perry et al., 2020). However, FDIs by themselves do not account for fuel or ignitions, differentiating them from fire observations such as burnt area, which makes them an unsuitable tool for assessing fire in the holistic context of weather, fuel dynamics, ignitions, and human land and fire management (Kelley and Harrison, 2014). Fire-enabled land surface models (LSMs) can, however, account for these drivers (Kelley and Harrison, 2014; Lasslop et al., 2016; Prentice et al., 2011b) to simulate a physical, observable measure of fire regimes, such as burnt area or number of fires. Some LSM fire schemes achieve this by modelling fuel moisture using FDIs (Lenihan et al., 1998; Rabin et al., 2017; Venevsky et al., 2002). However, most LSMs have been developed to study long-term, often decadal-timescale carbon dynamics and Earth system feedbacks and therefore often fail to reproduce year-to-year patterns of fire with the required accuracy to determine the causes of individual fire seasons (Andela et al., 2017; Hantson et al., 2016, 2020). Embedding fire within a complex vegetation model system also prevents rapid observation–model fusion, as iterative optimisation techniques are too computationally expensive and instabilities arise from non-linear responses of fire to simulated vegetation and fuel dynamics. Many large-scale vegetation-modelling projects, therefore, simulate up to a “present day” that can be several months or years out of date (Friedlingstein et al., 2019; Hantson et al., 2020). This lack of predictive capability has led to calls for simulation frameworks that fuse statistical representations of fire drivers with modelling techniques that consider such interactions (Fisher and Koven, 2020; Forkel et al., 2017; Krawchuk and Moritz, 2014; Sanderson and Fisher, 2020; Tollefson, 2018; Williams and Abatzoglou, 2016).
Kelley et al. (2019) developed a methodology which addresses this gap by coupling the same representation of processes found in simple fire-enabled LSMs (Rabin et al., 2017) using a Bayesian inference framework. The main advantage of this system is that it can assess the contribution of different fire drivers directly from observations and track uncertainty in the model parameters and the models' ability to reproduce observations. We apply this methodology here, using monthly meteorological conditions and burnt area (BA) observations to constrain and drive the model, thus capturing interannual variability within the context of changing meteorological conditions. We use this framework to answer the specific question: did the meteorological conditions contribute to the Amazonia fires of 2019?
The modelling protocol and optimisation framework largely follow Kelley et al. (2019), where a more detailed description may be found. Here, we summarise and outline further refinements. Monthly BA is modelled as a product of limitations imposed by four controls: (1) fuel availability; (2) moisture in live and dead fuel; (3) anthropogenic and natural ignition; (4) both active suppression and landscape fragmentation effects from human land use (Table S1 in the Supplement). Each control is calculated as a linear combination of its respective drivers. The impact each control has on fire is represented by a logistic curve describing the maximum allowed burnt area based on that control, and the product of all four limitations is used to determine BA.
We made several changes to the previous modelling protocol in order to utilise near-real-time meteorological and fire variables so that we can produce relevant results that closely follow the fire event and to better describe the conditional probability distribution between parameter samples and burnt area observations. We used the MODIS Collection 6 MCD64A1 burnt area product (Giglio et al., 2018) as our target dataset and replaced actual over potential evapotranspiration in the moisture control with soil moisture (Table S1). We also used both the top 10 cm and 10–200 cm soil moisture (Kalnay et al., 1996) as independent moisture drivers in order to capture the impact of previous drought years on deep water availability for live fuel. As near-real-time wet-day information is also not available, we replaced wet days in the calculation of dead fuel drying potential (Kelley et al., 2014) with a proxy for wet days (WDs), using Global Precipitation Climatology Project (GPCP) precipitation (Adler et al., 2003) (pr) based on (Prentice et al., 2011a)
where wd is an optimised parameter.
All variables were resampled and, where necessary, interpolated to a monthly time step as per Kelley et al. (2019). All driving variables were provided on a resolution of 2.5∘ except land use, provided at 0.5∘. We, therefore, choose to regrid all datasets to a resolution of 2.5∘, as interpolating to a finer resolution would provide no new information about the meteorological drivers tested. MCD64A1, soil moisture and equilibrium fuel moisture content were processed using the “rgdal” (Bivand et al., 2016) and “raster” (Hijmans and van Etten, 2014) packages in R (R Core Team, 2015). For MODIS vegetation continuous field (VCF) fractional covers (Dimiceli et al., 2015), tiles were merged and resampled to the model grid using the “gdal” package (GDAL/OGR contributors, 2018). Land use, population density, precipitation, humidity, temperature and lightning were processed using the iris package (Met Office, 2013) with Python version 3 (Python Software Foundation, https://www.python.org/, last access: 27 March 2019).
The model was optimised against MCD64A1 burnt area (Giglio et al., 2018) for the period July 2002 to June 2018, which was the common period for all datasets (Table S2) over South America, south of 13∘ N. We used the same Bayesian inference technique as per Kelley et al. (2019). Bayes' theorem states that the likelihood of the values of the set, β, which contain our 24 unknown parameters (i.e. the 21 parameters from Kelley et al. (2019), wd from Eq. (1) and error term parameters P0, σ in Eq. 3) and our known model inputs, given a set of observations Ys is proportional to the prior probability distribution of β (P(β)) multiplied by the conditional probability of Ys given β:
As 41.47 % of the burnt area observations are zero and the remaining are normally distributed under logit transformation (Fig. A1 in Appendix A), we, therefore, defined the likelihood, P(Ys|β) using a zero-inflated normal distribution on the logit transformed burnt area, as opposed to a simple normal distribution as used in Kelley et al. (2019). This better described the observational to the simulated difference in burnt area during times of very low or very high burning. Our zero-inflation likelihood term is therefore described as
where i represents an individual data point, yi is the burnt area observations, N is the observation sample size and .
The posterior solution was inferred for all model parameters using a Metropolis–Hastings Markov chain Monte Carlo (MCMC) step with the PyMC3 Python package (Salvatier et al., 2016), running 10 chains each over 10 000 iterations. We used all of the 44 750 grid cells on our 2.5∘ grid and monthly time step for 16 years in our assimilation procedure. This is a departure from Kelley et al. (2019), where only 10 % of grid cells were used, as our sample size was much smaller and we did not face the same computational demand. Due to our sample size, our likelihood dominates over our priors, and as with Kelley et al. (2019), priors predominantly were employed to set physically plausible bounds on our parameters.
Once optimised, the model was then run from January 2002–December 2019, and so the trained model was in a predictive mode for 2019. Due to data availability at the time of writing, July 2017–June 2018 land cover, land use and population density were recycled for July 2018 onwards (Table S2). We sampled 100 parameter ensemble members from the last 5000 iterations of each of the 10 chains, providing us with 1000 ensemble members to estimate the models' posterior solution to Eq. (2). Sampling was performed using the iris package (Met Office, 2013) with Python version 3 (Python Software Foundation, https://www.python.org/). The posterior solution provides an estimate of the burnt area based on the parameter uncertainty of our model, corresponding to the yellow areas in time series in Figs. 1 and 2. The mean burnt area for a particular parameter combination ( was obtained from
was evaluated using the implementation of the fire Model Intercomparison Project (fireMIP) benchmarking metrics (Hantson et al., 2020; Kelley et al., 2013; Rabin et al., 2017) as per Kelley et al. (2019). We also performed additional benchmarking metrics of the models' ability to reproduce seasonality and inter-annual variability of fire. See Figs. S1–S3 in the Supplement and the model evaluation Supplement for validation methods and results.
We chose five regions (marked A–E in Fig. 1 and Table 1; see Fig. A2 for locations) to represent forest areas already under pressure from deforestation. Regions A–C form a transect (west to east) across the agricultural–humid tropical forest interface in Brazil's arc of deforestation, often associated with deforestation (Fig. S4), whereas the D and E regions are found in agricultural regions embedded in savanna and grassland regions that experience frequent burning:
Acre, southern Amazonas state and Brazilian–Peruvian border
Rondônia and northern Mato Grosso, Brazil
Maranhão and Piau in coastal deforestation regions
Brazilian, Bolivian and Paraguayan border.
We also assessed an overall area of active deforestation (AAD) in the Amazon region (Fig. A2). This area is defined as the parts of South American southern tropics with significant decreasing tree cover trends, as seen in VCF (Dimiceli et al., 2015) and significant increasing agricultural fractions in the HYDEv3.1 dataset (Klein Goldewijk et al., 2010). Trend analysis used the same technique described in Kelley et al. (2019), where we took significance as p<0.05 on the linear trend for each month in the year on logit transformed cover variables, using the green–brown R package (Forkel et al., 2013, 2015). AAD was additionally assessed over three sub-region areas, primarily to evaluate the models' historic performance and assess the increase in 2019 fires across the humidity gradient:
area of active deforestation
the southern end of agricultural–humid tropical forest interface in Brazil's central states, often associated with the arc of deforestation in Brazil's central states
drier savanna and woodland in Cerrado and Caatinga in the eastern basin where land has already been heavily converted to agriculture
southern seasonally dry–deciduous Chiquitano and Gran Chaco forests, mainly along the Amazon and La Plata watersheds.
We assessed the probability of 2019 fire activity being explained by information provided to the model in three ways (Table 1):
The likelihood of observed monthly burnt area based on the information provided to our model: in the predictive model, the probability of a burnt area y (where y can be outside training data Ys, as is the case for our year 2019 analysis) being explained by our model (Pred(y) – full model uncertainty, or model error, in tan areas on time series in Figs. 1 and 2) is proportional to the probability of y given a parameter set, β, weighted by P(β| Ys):
Where the observed burnt area, y, falls within the model's full posterior, (D(y)) is then the sum of all probabilities greater than y:
When suggests y is towards the extremes of the posterior. As our posterior solution is not normally distributed, observations can fall at the extremes (i.e. when there is no burnt area, y=0 and by definition D(y)=1) and still have a high likelihood (i.e. if P(Ys=0| β) in Eq. (3) is much greater than 0; see Fig. B2 as an example). We, therefore, define the significance of D(y) as the probability of y occurring by chance (pv(y)) from the sum of all probabilities below P(y) (Fig. B2):
where is the set of solutions to P(yi)=P(y).
Whenever D(y) and pv(y) are low, this indicates burning is significantly higher than expected by the model in that month.
The likelihood that the burnt area would have been higher than the annual average, i.e. the fraction of the model's full posterior being greater than the model's annual average climatological posterior (the point where the vertical lines cross 1 in the right-hand columns in Figs. 1 and 2): a climatological burnt area “clim” for a given month, m, in the year (i.e. January, February, etc.) can be calculated from the convolution of each year's posterior solutions, βyr,m. Note that it is the model inputs, incorporated in βm, and not the model parameters that vary with time:
The probability of an anomaly A in a given year, yr, for month m, is, therefore
The likelihood of a year having a higher anomalous A is the sum of probabilities of <A:
And the likelihood of the year an average burnt area is given by .
The likelihood of the observed anomalous year occurring is given by
where ym,yr is the burnt area for the month, m, and year, yr, in question and is the climatological average of that month
3.1 Burnt area from 2001–2018
The highest historic burnt areas are found in the Savanna regions of tropical South America (Fig. 3), though some burning still occurs in forested areas, particularly in areas which have experienced an increase in agriculture and decrease in tree cover since 2002 (Fig. S4). The model reproduces this spatial pattern, and the model's full posterior encompasses the full range of burnt areas (Fig. 3, and benchmarking in the Supplement). Burnt area starts to increase in May and dies out in October throughout most of the area of active deforestation, although it can start as late as July in more humid areas and can continue through to December in drier Savanna (Fig. S3). The bulk of the burnt area occurs in August and September. September typically sees the highest burnt area in central Brazil, whereas fire peaks in August around Bolivia and Paraguay (Figs. 1E and S3). Our model reproduces this seasonal pattern in burning across all regions (see benchmarking in the Supplement), including onset and peak (Fig. S3). As our model maintains constant human ignitions and suppression throughout the year, this suggests that the seasonal pattern can be largely reproduced from meteorological variations, although a slight increase in uncertainty in early fire season burning could point to increased human ignitions not captured in the model (Fig. S3).
Unusually high levels of burning occurred in 2004 in the Bolivian–Paraguayan dry forest (red line in Fig. 1E), in 2005 in the eastern arc of deforestation (Fig. 1A) and Paraguay dry forest (Fig. 1E), in 2007 in monsoonal coastal forests (Fig. 1D), and in 2010 in Bolivia and Paraguay dry forests (Fig. 1D and E). The 2005 and 2010 burnings have previously been associated with droughts driven by a tropical North Atlantic warming anomaly (Marengo and Espinoza, 2016). The model reproduced the spatial pattern of this increased burning in 2005 and 2010 (Figs. C1 and C2 in Appendix C). In our different regions, observed levels of burning fall within our model's posterior (Figs. B1 and B2 in Appendix B), although they are at the higher end, with a high value of expectation (height of the posterior curve in Figs. B1 and B2) and high p value (blue shaded area, Figs. B1 and B2). This, along with the model's high Spearman's rank performance (Figs. S1 and S2) suggests that the model is able to capture the interannual variations driven by meteorological conditions. Deforestation rates in 2004/05 were high (Marengo et al., 2018), and an increase in fire activity in 2007 has also been linked to deforestation across the Amazon (Morton et al., 2008). Additionally, in the early part of our observational record, much of the region has been shown to be less coupled to meteorological drivers and more heavily influenced by human fire and land management (Aragão et al., 2018). This is reflected by the improved performance of the model, which depends solely on changing population density and land use cover and not on changes in landscape management, during this later period in the AAD (Fig. 2F, 2011 onwards), particularly in areas dominated by agriculture (Fig. 2H).
On the whole, the framework's posterior is better able to encompass extremes in observations in humid regions (Fig. 2G vs. H and I), particularly across the Brazilian arc of deforestation (Fig. 1A–C vs. D and E). Overall, 19 out of 204 months up to 2018 for the AAD (∼ 9 %) fall outside the 90 % confidence interval (tan in Fig. 2F), suggesting that the frameworks posterior accurately describes the occurrence of more extreme months for the region as a whole. That only 13 months out of (204 months × 5 regions) 1020 months (∼ 1 %) fall outside the posterior for smaller regions (Fig. 1) suggests that the posterior is wider than expected. Our assessments of mismatch between observations and model for these regions will, therefore, likely be conservative, particularly for humid regions B and C, with no months prior to 2019 falling outside the 90 % confidence interval.
3.2 Burnt area in 2019 in context
The year 2019 burnt area during the early fire season (defined as June to August) was higher than the 2002–2019 average in areas of recent historical deforestation, despite a lower than average burnt area over much of the rest of the continent (Fig. 3). The AAD as a whole saw the third highest levels of burning in the fire season (Fig. 2F, August), behind 2007 and 2010. Overall, 11 % of the AAD, particularly around the central region of Brazil's arc of deforestation, experienced more burning in August than in any previous year since our observation record started in 2002 (Fig. 1B and C). Despite burnt area returning to normal levels in September across most of the AAD (Fig. 4), burning remained high in humid forest areas (Fig. 2G). Burning also remained higher than average along the border between Brazil, Bolivia and Paraguay (Figs. 1E, 2I and 4). This meant that, while the burnt area was higher than usual in 2019, it was not exceptionally higher over the entire fire season (June–September) for the entire AAD, though individual regions still stand out as having much higher burning than in any previous year, particularly in Brazil's central states (Fig. 1B and C).
3.3 Climatic conditions in 2019
The model shows with high confidence that most of eastern Amazonia should have, in fact, experienced less fire than normal for June–August when accounting for 2019 meteorological conditions. This expected low fire rate included areas in the Brazilian humid-forest–agricultural interface in the AAD that saw higher than annual average burning in observations (Figs. 1C and 3). The western Amazon shows an increase in fire compared to the annual average (Fig. 3). The observed burnt area, however, still exceeds the model in all our regions in 2019 (Fig. 1) except region D in the already heavily converted agricultural land near the Brazilian coast (Fig. S4). Only 0 %–2 % of the area of the AAD showed unprecedented high burning in the model, compared to the 11 % in observations. Observed burnt area in August falls outside the full model's posterior (at 90 % confidence interval) for the AAD (Fig. 2F), with only a 9 % likelihood of being explained by the model (Table 1). This is particularly prominent in more humid areas, with a 10 % likelihood in humid forests (Fig. 2G and Table 1, G) – tied as lowest likelihood in the observational record with September 2005, compared to 17 % likelihood in seasonally dry forests regions (Fig. 2I) and 17 % the driest, savanna areas (Fig. 2H) where observations tend to fall outside the model posterior more often (see Sect. 3.1). Regions B, C in central Brazil, and E on the Brazil, Bolivia and Paraguay border are even less likely to be explained by the model (7 % for B, 6 % for C and 8 % for E; Table 1), despite all previous months falling within the full model's posterior confidence range in these regions, except for August and September in the 2005 in region A, although it is more likely that burnt area regions A and D at either end of Brazil's arc of deforestation could be explained by the model (18 % and 20 %, respectively).
The observed anomaly for August 2019 is higher than the model across all regions except D. This is particularly prominent in regions B and C, where observations show that burnt area was 196 % and 138 % greater than the August average (Table 1), whereas the model suggests that meteorological conditions alone should have resulted in a fire season with a 16 %–22 % increase (based on 5 %–95 % parameter uncertainty range for parameter uncertainty) in burnt area in B and 2 % reduction to 4 % increase for C compared to the August average, with only a 57 % and 53 % chance of a greater burnt area than the average for B and C, respectively. The likely occurrence of the observed anomaly was 7 % and 10 % for B and C, respectively (Table 1) – much greater than any previous year (Fig. 1B and C, August column).
The higher observed anomaly vs. the model extends over much of the AAD (Fig. 2, August column, red points). The model suggests a 4 %–6 % reduction for the AAD, with a 49 % probability of greater than the annual average burnt area (Table 1). By comparison, the observed burnt area was 45 % greater than the annual average, with a 20 % likelihood. Again, the observed anomaly seems to be least likely in more humid regions. For our humid area, G, the model suggests a small (10 %–14 %) increase in burnt area, with only a 12 % probability of a 107 % increase seen in observations, whereas the 5 % observed reduction in drier savanna regions in H seems to be in line with the model (at 67 % likelihood).
By September, there was less disparity between observations and models. In regions B and C, for example, observations had a 14 % and 26 % probability of being explained by the model and a 66 % and 25 % likelihood of the anomalous year as seen in the observations. Across the AAD, the more modest 16 % increase in the observations had a likelihood of 40 %.
The observed spatial pattern of burnt area in June–August 2019 shows that unprecedented burning was only seen in Brazilian regions normally associated with deforestation. Our modelling framework demonstrates that, based on meteorological conditions alone, reduced burning seen across the rest of tropical South America should have extended into these regions. Specifically, our analysis suggests that there is only a 9 % probability that the levels of burning in the early fire season would have been caused by 2019 meteorological conditions or natural ignitions alone (time series in Figs. 1 and 2). Eastern areas normally associated with deforestation did show expected levels of burning, but in the western and central parts of the arc of deforestation and Bolivia and Paraguay dry forests, burning was much higher. Here there is a 6 %–8 % of such high levels of burning compared to the background rate (Fig. 1, August column), with areas where agriculture meets more humid forest seeing the most unusual levels of burning. As our model's posterior reflects the levels of burning in previous dry years, we can eliminate drier conditions as a possible driver of increased 2019 fires. We also account for, and therefore eliminate, longer-term drier conditions through deep soil moisture as a possible driver. The cause of increased burning in 2019 is therefore either a driver left static in the model for 2019 or a process not considered. Because of the non-availability of near-real-time data, drivers held unchanged at 2018 values for 2019 are tree cover, land use and human population. The only plausible way tree cover could have substantially changed is through increased deforestation rates (Zhang et al., 2015). Thereby changes in drivers not accounted for in 2019 would only have caused increased burning through direct human manipulation of the landscape rather than the particular meteorological features of that year.
Improved descriptions of evolving changes in human fire and landscape interactions over time may also be required to capture direct human-driven changes in burnt area. This is likely to include changes in demography or human behaviour. For example, we currently account for the impact of a changing population on fire starts and suppression but not for how fire ignitions per person change over time. An evolving policy could have also been the cause of the unusual fire activity. It should also be noted that observed fire activity returned to expected levels given meteorological conditions in September over most of the deforestation region (Figs. 1 and 4). This reduction could be after the June–August fires received international media coverage, triggering efforts in combating fires from South American governments (BBC news, 2019; NASA, 2019).
In this study, we have used a novel Bayesian modelling approach, which tests the likelihood of observed extremes in fire against inferred historical relationships by tracking uncertainties in modelling fire in the land surface. Our framework provides a rapid assessment of whether there was any influence of meteorological conditions across the Amazon that exacerbated fire levels in 2019.
The model predicts a lower burnt area than we see in the observations for Amazonia during June–August 2019, indicating that from observed meteorological data alone, we would not expect 2019 to be a high-fire year. This result points to socio-economic factors having a strong role in the high recorded fire activity. Specifically, we conclude that it is likely (> 90 % probability) based on past relationships between burnt area and meteorological conditions that the weather conditions did not trigger the increase in burning in Amazonia during the early fire season in 2019. This result holds over the entire area of active deforestation and furthermore is likely (93 %) in central Amazonia.
The model code and Bayesian inference framework used to support the findings of this study are archived at https://doi.org/10.5281/zenodo.4298760 (Kelley et al., 2020a). Model output is archived at https://doi.org/10.5281/zenodo.4298902 (Kelley et al., 2020b). Driving data availability is listed in Table S1.
The supplement related to this article is available online at: https://doi.org/10.5194/bg-18-787-2021-supplement.
DIK, CB, MAJB and CH developed the modelling framework; DIK, RW designed the Bayesian inference framework; DIK, MAJB and ND collated and regridded input data. DIK, CB, CH performed the analysis. CB wrote the first draft of the paper with input from DIK, CH and ND. All authors contributed to the final paper.
The authors declare that they have no conflict of interest.
GPCP precipitation and NCEP reanalysis data were provided by the NOAA/OAR/ESRL PSD, Boulder, Colorado, USA, from their website at https://www.esrl.noaa.gov/psd/ (last access: 20 January 2020).
The contribution by Douglas I. Kelley and Megan A. J. Brown was supported by the UK Natural Environment Research Council through The UK Earth System Modelling Project (UKESM, grant no. NE/N017951/1). Chris Huntingford was supported by the NERC CEH National Capability Fund. Ning Dong is supported by the Australia Research Council (DP170103410) and Macquarie university research seeding grant to Melanie Zeppel. Megan A. J. Brown was additionally supported by the Science and Technology Facilities Council (grant no. ST/P006760/1) through the DISCnet Centre for Doctoral Training.
This paper was edited by Sönke Zaehle and reviewed by Sönke Zaehle and one anonymous referee.
Adler, R. F., Huffman, G. J., Chang, A., Ferraro, R., Xie, P.-P., Janowiak, J., Rudolf, B., Schneider, U., Curtis, S., Bolvin, D., Gruber, A., Susskind, J., Arkin, P., and Nelkin, E.: The Version-2 Global Precipitation Climatology Project (GPCP) Monthly Precipitation Analysis (1979–Present), J. Hydrometeorol., 4, 1147–1167, 2003.
Aldersley, A., Murray, S. J., and Cornell, S. E.: Global and regional analysis of climate and human drivers of wildfire, Sci. Total Environ., 409, 3472–3481, 2011.
Andela, N., Morton, D. C., Giglio, L., Chen, Y., van der Werf, G. R., Kasibhatla, P. S., DeFries, R. S., Collatz, G. J., Hantson, S., Kloster, S., Bachelet, D., Forrest, M., Lasslop, G., Li, F., Mangeon, S., Melton, J. R., Yue, C., and Randerson, J. T.: A human-driven decline in global burned area, Science, 356, 1356–1362, 2017.
Anon: 2020 Incident Archive, CAL FIRE, available at: https://www.fire.ca.gov/incidents/2020/ (last access: 10 April 2019), 2020.
Aragão, L. E. O. C., Anderson, L. O., Fonseca, M. G., Rosan, T. M., Vedovato, L. B., Wagner, F. H., Silva, C. V. J., Silva Junior, C. H. L., Arai, E., Aguiar, A. P., Barlow, J., Berenguer, E., Deeter, M. N., Domingues, L. G., Gatti, L., Gloor, M., Malhi, Y., Marengo, J. A., Miller, J. B., Phillips, O. L., and Saatchi, S.: 21st Century drought-related fires counteract the decline of Amazon deforestation carbon emissions, Nat. Commun., 9, 536, https://doi.org/10.1038/s41467-017-02771-y, 2018.
Arruda, D., Candido, H. G., and Fonseca, R.: Amazon fires threaten Brazil's agribusiness, Science, 365, 1387, https://doi.org/10.1126/science.aaz2198, 2019.
BBC news: Amazon fires: Brazil bans land clearance blazes for 60 days, BBC news, available at: https://www.bbc.co.uk/news/world-latin-america-49507405 (last access: 3 November 2020), 2019.
Bivand, R., Keitt, T., and Rowlingson, B.: rgdal: Bindings for the Geospatial Data Abstraction Library, available at: http://CRAN.R-project.org/package=rgdal, last access: 10 November 2016.
Boer, M. M., de Dios, V. R., and Bradstock, R. A.: Unprecedented burn area of Australian mega forest fires, Nat. Clim. Change, 10, 171–172, 2020.
Brando, P. M., Balch, J. K., and Nepstad, D. C.: Abrupt increases in Amazonian tree mortality due to drought–fire interactions, P. Natl. Acad. Sci. USA, 111, 6347–6352, https://doi.org/10.1073/pnas.1305499111, 2014.
Burton, C., Betts, R. A., and Jones, C. D.: Will fire danger be reduced by using Solar Radiation Management to limit global warming to 1.5 ∘C compared to 2.0 ∘C?, Geophys. Res. Lett., 45, 3644–3652, https://doi.org/10.1002/2018GL077848, 2018.
Butler, R.: Calculating deforestation figures for the Amazon, Rio de Janeiro, available at: https://rainforests.mongabay.com/amazon/deforestation_calculations.html (last access: 29 February 2020), 2017.
Cochrane, M. A. and Schulze, M. D.: Fire as a Recurrent Event in Tropical Forests of the Eastern Amazon: Effects on Forest Structure, Biomass, and Species Composition1, Biotropica, 31, 2–16, 1999.
De Groot, W. J., Goldammer, J. G., Justice, C. O., and Lynham, T. J.: Implementing a global early warning system for wildland fire, in: Proc. Virchow. Pirquet. Med. Soc., University of Coimbra, Coimbra, Portugal, 15–18, 2010.
de Groot, W. J., Wotton, B. M., and Flannigan, M. D.: Chapter 11 – Wildland Fire Danger Rating and Early Warning Systems, in: Wildfire Hazards, Risks and Disasters, edited by: Shroder, J. F. and Paton, D., Elsevier, Oxford, 207–228, 2015.
de Oliveira Andrade, R.: Alarming surge in Amazon fires prompts global outcry, Nature, https://doi.org/10.1038/d41586-019-02537-0, 2019.
Dimiceli, C., Carroll, M., Sohlberg, R., Kim, D. H., Kelly, M., and Townshend, J. R. G.: MOD44B MODIS/Terra Vegetation Continuous Fields Yearly L3 Global 250 m SIN Grid V006, NASA EOSDIS Land Processes DAAC, https://doi.org/10.5067/MODIS/MOD44B.006, 2015.
Escobar, H.: Amazon fires clearly linked to deforestation, scientists say, Science, 365, 853, https://doi.org/10.1126/science.365.6456.853, 2019.
Fisher, R. A. and Koven, C. D.: Perspectives on the future of Land Surface Models and the challenges of representing complex terrestrial systems, J. Adv. Model. Earth Sy., 12, e2018MS001453, https://doi.org/10.1029/2018MS001453, 2020.
Fonseca, M. G., Alves, L. M., Aguiar, A. P. D., Arai, E., Anderson, L. O., Rosan, T. M., Shimabukuro, Y. E., and Aragão, L. E. O. C.: Effects of climate and land-use change scenarios on fire probability during the 21st century in the Brazilian Amazon, Glob. Change Biol., 25, 2931–2946, 2019.
Forkel, M., Carvalhais, N., Verbesselt, J., Mahecha, M. D., Neigh, C., and Reichstein, M.: Trend Change Detection in NDVI Time Series: Effects of Inter-Annual Variability and Methodology, Remote Sensing, 5, 2113–2144, 2013.
Forkel, M., Migliavacca, M., Thonicke, K., Reichstein, M., Schaphoff, S., Weber, U., and Carvalhais, N.: Co-dominant water control on global inter-annual variability and trends in land surface phenology and greenness, Glob. Change Biol., 21, 3414–3435, 2015.
Forkel, M., Dorigo, W., Lasslop, G., Teubner, I., Chuvieco, E., and Thonicke, K.: A data-driven approach to identify controls on global fire activity from satellite and climate observations (SOFIA V1), Geosci. Model Dev., 10, 4443–4476, https://doi.org/10.5194/gmd-10-4443-2017, 2017.
Friedlingstein, P., Jones, M. W., O'Sullivan, M., Andrew, R. M., Hauck, J., Peters, G. P., Peters, W., Pongratz, J., Sitch, S., Le Quéré, C., Bakker, D. C. E., Canadell, J. G., Ciais, P., Jackson, R. B., Anthoni, P., Barbero, L., Bastos, A., Bastrikov, V., Becker, M., Bopp, L., Buitenhuis, E., Chandra, N., Chevallier, F., Chini, L. P., Currie, K. I., Feely, R. A., Gehlen, M., Gilfillan, D., Gkritzalis, T., Goll, D. S., Gruber, N., Gutekunst, S., Harris, I., Haverd, V., Houghton, R. A., Hurtt, G., Ilyina, T., Jain, A. K., Joetzjer, E., Kaplan, J. O., Kato, E., Klein Goldewijk, K., Korsbakken, J. I., Landschützer, P., Lauvset, S. K., Lefèvre, N., Lenton, A., Lienert, S., Lombardozzi, D., Marland, G., McGuire, P. C., Melton, J. R., Metzl, N., Munro, D. R., Nabel, J. E. M. S., Nakaoka, S.-I., Neill, C., Omar, A. M., Ono, T., Peregon, A., Pierrot, D., Poulter, B., Rehder, G., Resplandy, L., Robertson, E., Rödenbeck, C., Séférian, R., Schwinger, J., Smith, N., Tans, P. P., Tian, H., Tilbrook, B., Tubiello, F. N., van der Werf, G. R., Wiltshire, A. J., and Zaehle, S.: Global Carbon Budget 2019, Earth Syst. Sci. Data, 11, 1783–1838, https://doi.org/10.5194/essd-11-1783-2019, 2019.
GDAL/OGR contributors: GDAL/OGR Geospatial Data Abstraction software Library, Open Source Geospatial Foundation, Version 3.2, available at: http://gdal.org (last access: March 2020), 2018.
Giglio, L., Boschetti, L., Roy, D. P., Humber, M. L., and Justice, C. O.: The Collection 6 MODIS burned area mapping algorithm and product, Remote Sens. Environ., 217, 72–85, 2018.
Hantson, S., Arneth, A., Harrison, S. P., Kelley, D. I., Prentice, I. C., Rabin, S. S., Archibald, S., Mouillot, F., Arnold, S. R., Artaxo, P., Bachelet, D., Ciais, P., Forrest, M., Friedlingstein, P., Hickler, T., Kaplan, J. O., Kloster, S., Knorr, W., Lasslop, G., Li, F., Mangeon, S., Melton, J. R., Meyn, A., Sitch, S., Spessa, A., van der Werf, G. R., Voulgarakis, A., and Yue, C.: The status and challenge of global fire modelling, Biogeosciences, 13, 3359–3375, https://doi.org/10.5194/bg-13-3359-2016, 2016.
Hantson, S., Kelley, D. I., Arneth, A., Harrison, S. P., Archibald, S., Bachelet, D., Forrest, M., Hickler, T., Lasslop, G., Li, F., Mangeon, S., Melton, J. R., Nieradzik, L., Rabin, S. S., Prentice, I. C., Sheehan, T., Sitch, S., Teckentrup, L., Voulgarakis, A., and Yue, C.: Quantitative assessment of fire and vegetation properties in simulations with fire-enabled vegetation models from the Fire Model Intercomparison Project, Geosci. Model Dev., 13, 3299–3318, https://doi.org/10.5194/gmd-13-3299-2020, 2020.
Hijmans, R. J. and van Etten, J.: raster: Geographic data analysis and modeling, R package version, available at: http://CRAN.R-project.org/package=raster (last access: 15 March 2018), 2014.
INPE: INPE Fire Count Data, available at: http://queimadas.dgi.inpe.br/queimadas/portal-static/estatisticas_paises/, last access: 29 February 2019.
Jolly, W. M., Cochrane, M. A., Freeborn, P. H., Holden, Z. A., Brown, T. J., Williamson, G. J., and Bowman, D. M. J. S.: Climate-induced variations in global wildfire danger from 1979 to 2013, Nat. Commun., 6, 7537, https://doi.org/10.1038/ncomms8537, 2015.
Jones, M. W., Smith, A., Betts, R., Canadell, J. G., Prentice, I. C., and Le Queìreì, C.: Climate change increases the risk of wildfires, ScienceBrief Review, available at: https://ueaeprints.uea.ac.uk/id/eprint/77982/, last access: 29 February 2020.
Kalnay, E., Kanamitsu, M., Kistler, R., Collins, W., Deaven, D., Gandin, L., Iredell, M., Saha, S., White, G., Woollen, J., Zhu, Y., Chelliah, M., Ebisuzaki, W., Higgins, W., Janowiak, J., Mo, K. C., Ropelewski, C., Wang, J., Leetmaa, A., Reynolds, R., Jenne, R., and Joseph, D.: The NCEP/NCAR 40-Year Reanalysis Project, B. Am. Meteorol. Soc., 77, 437–472, 1996.
Kelley, D. I.: Modelling Australian fire regimes, PhD thesis, Macquarie, University, Ryde, NSW, Australia, 2014.
Kelley, D. I. and Harrison, S. P.: Enhanced Australian carbon sink despite increased wildfire during the 21st century, Environ. Res. Lett., 9, 104015, https://doi.org/10.1088/1748-9326/9/10/104015, 2014.
Kelley, D. I., Prentice, I. C., Harrison, S. P., Wang, H., Simard, M., Fisher, J. B., and Willis, K. O.: A comprehensive benchmarking system for evaluating global vegetation models, Biogeosciences, 10, 3313–3340, https://doi.org/10.5194/bg-10-3313-2013, 2013.
Kelley, D. I., Harrison, S. P., and Prentice, I. C.: Improved simulation of fire–vegetation interactions in the Land surface Processes and eXchanges dynamic global vegetation model (LPX-Mv1), Geosci. Model Dev., 7, 2411–2433, https://doi.org/10.5194/gmd-7-2411-2014, 2014.
Kelley, D. I., Bistinas, I., Whitley, R., Burton, C., Marthews, T. R., and Dong, N.: How contemporary bioclimatic and human controls change global fire regimes, Nat. Clim. Change, 9, 690–696, 2019.
Kelley, D. I., Whitley, R., Brown, M. A. J., and Burton, C.: Scripts used in the submission of “Low Climatic Influence found in 2019 Amazonia Fires”, Zenodo, https://doi.org/10.5281/zenodo.4298760, 2020a.
Kelley, D. I., Whitley, R., Burton, C., Brown, M. A. J., Huntingford, C., and Dong, N.: ConFire Model input/output for South America, Zenodo, https://doi.org/10.5281/zenodo.4298902, 2020b.
Klein Goldewijk, K., Goldewijk, K. K., Beusen, A., Van Drecht, G., and De Vos, M.: The HYDE 3.1 spatially explicit database of human-induced global land-use change over the past 12,000 years, Global Ecol. Biogeogr., 20, 73–86, 2010.
Krawchuk, M. A. and Moritz, M. A.: Burning issues: statistical analyses of global fire data to inform assessments of environmental change, Environmetrics, 25, 472–481, 2014.
Lasslop, G., Brovkin, V., Reick, C. H., Bathiany, S., and Kloster, S.: Multiple stable states of tree cover in a global land surface model due to a fire-vegetation feedback, Geophys. Res. Lett., 43, 6324–6331, 2016.
Lenihan, J. M., Daly, C., Bachelet, D., and Neilson, R. P.: Simulating broad-scale fire severity in a dynamic global vegetation model, Northwest Sci., 72, 91–101, 1998.
Lizundia-Loiola, J., Pettinari, M. L., and Chuvieco, E.: Temporal Anomalies in Burned Area Trends: Satellite Estimations of the Amazonian 2019 Fire Crisis, Remote Sens.-Basel, 12, 151, https://doi.org/10.3390/rs12010151, 2020.
Lovejoy, T. E. and Nobre, C.: Amazon tipping point: Last chance for action, Sci. Adv., 5, eaba2949, https://doi.org/10.1126/sciadv.aba2949, 2019.
Marengo, J. A. and Espinoza, J. C.: Extreme seasonal droughts and floods in Amazonia: causes, trends and impacts: EXTREMES IN AMAZONIA, Int. J. Climatol., 36, 1033–1050, 2016.
Marengo, J. A., Souza, C. A., Thonicke, K., Burton, C., Halladay, K., Betts, R. A., and Soares, W. R.: Changes in climate and land use over the Amazon Region: current and future variability and trends, Front Earth Sci. Chin., 6, 228, https://doi.org/10.3389/feart.2018.00228, 2018.
Met Office: Iris: A Python library for analysing and visualising meteorological and oceanographic data sets, v1.2 edn., Exeter, Devon, available at: http://scitools.org.uk/ (last access: 8 October 2019), 2013.
Morton, D. C., Defries, R. S., Randerson, J. T., Giglio, L., Schroeder, W., and Van Der Werf, G. R.: Agricultural intensification increases deforestation fire activity in Amazonia, Glob. Change Biol., 14, 2262–2275, 2008.
NASA: Uptick in Amazon Fire Activity in 2019, NASA Earth Observatory, available at: https://www.earthobservatory.nasa.gov/images/145498/uptick-in-amazon-fire-activity-in-2019 (last access: 3 November 2020), 2019.
Nauslar, N. J., Brown, T. J., McEvoy, D. J., and Lareau, N. P.: Record-setting 2018 California wildfires, in: State of the Climate 2018, edited by: Blunden, J. A. D. S. A., Special Supplement to the Bulletin of the American Meteorological Society Vol. 100, No. 9, September 2019, p. 195, 2019.
Patel, K.: Arctic Fires Fill the Skies with Soot, NASA Earth Observatory, available at: https://earthobservatory.nasa.gov/images/145380/arctic-fires-fill-the-skies-with-soot, last access: 1 December 2019.
Pellegrini, A. F. A., Anderegg, W. R. L., Paine, C. E. T., Hoffmann, W. A., Kartzinel, T., Rabin, S. S., Sheil, D., Franco, A. C., and Pacala, S. W.: Convergence of bark investment according to fire and climate structures ecosystem vulnerability to future change, Ecol. Lett., 20, 307–316, 2017.
Perry, J. J., Cook, G. D., Graham, E., Meyer, C. P. M., Murphy, H. T., and VanDerWal, J.: Regional seasonality of fire size and fire weather conditions across Australia's northern savanna, Int. J. Wildland Fire, 29, 1–10, 2020.
Prentice, I. C., Harrison, S. P., and Bartlein, P. J.: Global vegetation and terrestrial carbon cycle changes after the last ice age, New Phytol., 189, 988–998, 2011a.
Prentice, I. C., Kelley, D. I., Foster, P. N., Friedlingstein, P., Harrison, S. P., and Bartlein, P. J.: Modeling fire and the terrestrial carbon balance, Global Biogeochem. Cy., 25, GB3005, https://doi.org/10.1029/2010gb003906, 2011b.
Rabin, S. S., Melton, J. R., Lasslop, G., Bachelet, D., Forrest, M., Hantson, S., Kaplan, J. O., Li, F., Mangeon, S., Ward, D. S., Yue, C., Arora, V. K., Hickler, T., Kloster, S., Knorr, W., Nieradzik, L., Spessa, A., Folberth, G. A., Sheehan, T., Voulgarakis, A., Kelley, D. I., Prentice, I. C., Sitch, S., Harrison, S., and Arneth, A.: The Fire Modeling Intercomparison Project (FireMIP), phase 1: experimental and analytical protocols with detailed model descriptions, Geosci. Model Dev., 10, 1175–1197, https://doi.org/10.5194/gmd-10-1175-2017, 2017.
R Core Team: R: A Language and Environment for Statistical Computing, available at: http://www.R-project.org/, last access: 16 April 2015.
RFS: New South Wales Rural Fire Service, available at: https://www.rfs.nsw.gov.au/fire-information/fires-near-me, last access: 29 February 2019.
Salvatier, J., Wiecki, T. V., and Fonnesbeck, C.: Probabilistic programming in Python using PyMC3, PeerJ. Comput. Sci., 2, e55, https://doi.org/10.7717/peerj-cs.55, 2016.
Sanderson, B. M. and Fisher, R. A.: A fiery wake-up call for climate science, Nat. Clim. Change, 10, 175–177, 2020.
Staver, A. C., Brando, P. M., Barlow, J., Morton, D. C., Paine, C. E. T., Malhi, Y., Araujo Murakami, A., and Del Aguila Pasquel, J.: Thinner bark increases sensitivity of wetter Amazonian tropical forests to fire, Ecol. Lett., 23, 99–106, 2020.
Tollefson, J.: ECOLOGY Huge wildfires defy explanation, Nature, 561, 16–17, 2018.
Van Der Werf, G. R., Randerson, J. T., Giglio, L., Gobron, N., and Dolman, A. J.: Climate controls on the variability of fires in the tropics and subtropics, Global Biogeochem. Cy., 22, GB3028, https://doi.org/10.1029/2007GB003122, 2008.
van Oldenborgh, G. J., Krikken, F., Lewis, S., Leach, N. J., Lehner, F., Saunders, K. R., van Weele, M., Haustein, K., Li, S., Wallom, D., Sparrow, S., Arrighi, J., Singh, R. P., van Aalst, M. K., Philip, S. Y., Vautard, R., and Otto, F. E. L.: Attribution of the Australian bushfire risk to anthropogenic climate change, Nat. Hazards Earth Syst. Sci. Discuss. [preprint], https://doi.org/10.5194/nhess-2020-69, in review, 2020.
Venevsky, S., Thonicke, K., Sitch, S., and Cramer, W.: Simulating fire regimes in human-dominated ecosystems: Iberian Peninsula case study, Glob. Change Biol., 8, 984–998, 2002.
Williams, A. P. and Abatzoglou, J. T.: Recent Advances and Remaining Uncertainties in Resolving Past and Future Climate Effects on Global Fire Activity, Current Climate Change Reports, 2, 1–14, 2016.
Wiltshire, J., Hekman, J., and Milan, B. F.: Carbon loss and economic impacts of a peatland wildfire in north-east Sutherland, Scotland, in: Report on impact of wildfire on the peatland and released carbon from the peatland, 12–17 May 2019, Ricardo Energy & Environment for WWF-UK, 2019.
Witze, A.: The Arctic is burning like never before – and that's bad news for climate change, Nature, 585, 336–337, 2020.
Zeppel, M. J. B., Harrison, S. P., Adams, H. D., Kelley, D. I., Li, G., Tissue, D. T., Dawson, T. E., Fensham, R., Medlyn, B. E., Palmer, A., West, A. G., and McDowell, N. G.: Drought and resprouting plants, New Phytol., 206, 583–589, 2015.
Zhang, K., de Almeida Castanho, A. D., Galbraith, D. R., Moghim, S., Levine, N. M., Bras, R. L., Coe, M. T., Costa, M. H., Malhi, Y., Longo, M., Knox, R. G., McKnight, S., Wang, J., and Moorcroft, P. R.: The fate of Amazonian ecosystems over the coming century arising from changes in climate, atmospheric CO2, and land use, Glob. Change Biol., 21, 2569–2587, https://doi.org/10.1111/gcb.12903, 2015.