{'fraction_non_alphanumeric': 0.03027137338742985, 'fraction_numerical': 0.030975923811374845, 'mean_word_length': 4.871897289586305, 'pattern_counts': {'":': 0, '<': 0, '<?xml version=': 0, '>': 3, 'https://': 1, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 0, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Fire and deforestation dynamics in Amazonia (1973–2014)

Abstract Consistent long‐term estimates of fire emissions are important to understand the changing role of fire in the global carbon cycle and to assess the relative importance of humans and climate in shaping fire regimes. However, there is limited information on fire emissions from before the satellite era. We show that in the Amazon region, including the Arc of Deforestation and Bolivia, visibility observations derived from weather stations could explain 61% of the variability in satellite‐based estimates of bottom‐up fire emissions since 1997 and 42% of the variability in satellite‐based estimates of total column carbon monoxide concentrations since 2001. This enabled us to reconstruct the fire history of this region since 1973 when visibility information became available. Our estimates indicate that until 1987 relatively few fires occurred in this region and that fire emissions increased rapidly over the 1990s. We found that this pattern agreed reasonably well with forest loss data sets, indicating that although natural fires may occur here, deforestation and degradation were the main cause of fires. Compared to fire emissions estimates based on Food and Agricultural Organization's Global Forest and Resources Assessment data, our estimates were substantially lower up to the 1990s, after which they were more in line. These visibility‐based fire emissions data set can help constrain dynamic global vegetation models and atmospheric models with a better representation of the complex fire regime in this region.


Introduction
In their natural state, tropical rainforests rarely burn, because of moist conditions underneath the canopy and because dry lightning is rare [Cochrane, 2003]. Human settlements leading to, for example, agriculture and road development in general lead to increased ignitions [Moran, 1993]. One way is the slash-and-burning method where small-scale patches of forest are burned and used for several years for agricultural purposes. At a later stage, large-scale deforestation aided by fire occurred, where forest is removed permanently for industrial-scale crop production. Currently, fire-driven deforestation is the main source of carbon emissions in Amazonia [Morton et al., 2006[Morton et al., , 2008Cochrane, 2009;van der Werf et al., 2010]. It is often assumed that fire emissions have increased over time here, but the timing and pattern of this increase is not well known. Anthropogenic fire regimes in the rainforests make the ecosystem more vulnerable for fire [Nepstad et al., 1999] and may increase the occurrence of fire-prone conditions due to the effects of emitted aerosols on precipitation formation [Bevan et al., 2009;Gonçalves et al., 2015].
The Amazon region has been deforested by about 15% from 1976 to 2010 and deforestation is one of the drivers of changing energy and water cycles in the southern and eastern portions of the Amazon basin [Davidson et al., 2012;Aragão et al., 2014]. Satellite-based forest loss data sets show substantial interannual variability (IAV) over the 1990s and 2000s and an overall decrease after 2005 as a result of reduced deforestation rates in Brazil due to conservation policies, namely, stricter law enforcements and expansion of the governmental protection of the Amazon area [Nepstad et al., 2009;Hansen et al., 2013]. Although forest loss has decreased in Brazil, neighboring Argentina, Paraguay, and Bolivia show an increase over the 2000s [Chen et al., 2013;van Marle et al., 2016]. Nowadays, deforestation in the Amazon often goes hand in hand with fire, emitting particulate matter into the atmosphere, which reduces air quality and can affect human health [Johnston et al., 2012;Marlier et al., 2012]. Reddington et al. [2015] showed that the reduction in Brazilian fire emissions related to the decrease in deforestation fires prevented 400 to 1700 premature adult deaths annually. VAN  To understand how fires influence and interact with the Earth system and specifically the atmosphere, quantitative information on emissions and a breakdown of emissions into different sources is required. These can be derived using satellite-observed burned area or active fires potentially in combination with biogeochemical modeling, as done in, for example, the Global Fire Emissions Database (GFED) or the Global Fire Assimilation System [van der Werf et al., 2010;Kaiser et al., 2012]. Alternatively, atmospheric concentrations of trace gases and aerosols, partly influenced by fire emissions, can also be observed by satellites directly. The Measurements Of Pollution In The Troposphere (MOPITT) sensor is often used for this purpose and observes carbon monoxide (CO) from space. It has shown to be related to South American biomass burning emissions [Edwards et al., 2006]. Other examples include the Microwave Limb Sounder that can observe variations in CO and ozone with increased CO emissions over the Amazon region for 2004, 2005, 2007[Livesey et al., 2013 and the aerosol index observed by the various sensors.
However, these data sets only provide data for the satellite era, which starts for the sensors or data sets mentioned above in 1997 or later. Information about fire emissions in the historic record before that time is mainly based on country-level statistics. A widely used metric is the deforestation and afforestation inventory of the Food and Agricultural Organizationˈs (FAO) Global Forest and Resources Assessment (FRA) [Food and Agriculture Organization of the United Nations, 2012], where country-scale estimates of forest loss and gain area are converted to carbon emissions via a bookkeeping model [Houghton, 2003]. While several proxies exist that provide information before the satellite period (e.g., charcoal-based records from lake-sediment and levoglucosan from ice core records), for the Amazon region only charcoal [Power et al., 2010;Marlon et al., 2013] and historical [Mouillot and Field, 2005] records have been constructed. The latter suggests that burned area or fire in general has increased slowly over the 1960s and 1970s and increased more strongly over the 1980s and 1990s. This data set lacks information about interannual variability (IAV), which limits our ability to understand the contributions of different possible drivers of fire activity, for example, land use and climatic factors.
In this study we used horizontal visibility data observed by the World Meteorological Organization (WMO) weather stations throughout the Amazon region to reconstruct a fire history for this region since 1973. Horizontal visibility are human observations using landmarks with known distances during the day and using point light sources at night [World Meteorological Organization, 1996]. Horizontal visibility is influenced by several sources including dust, air pollution, haze, fog, and precipitation. Fires also have a strong impact on visibility as shown in Figure 1: an image of Bolivian fires taken from the Moderate Resolution Imaging Spectroradiometer (MODIS) during the 2010 fire season. Visibility data sets are inherently more subjective than more conventional weather station records but have shown to be useful in a wide range of studies constraining regional air quality changes [Doyle and Dorling, 2002;Che et al., 2007;Chang et al., 2009], desert dust emissions [Mahowald et al., 2007], and global trends in atmospheric haze [Li et al., 2016].
Our goal was to determine the extent to which visibility observations could be used to constrain fire emissions in the main deforestation region of the Amazon. This has shown to be the case for peat fires in Indonesia [Heil and Goldammer, 2001;Field et al., 2009Field et al., , 2016. The fire regime in Indonesia is strongly linked to anthropogenic use of fire for deforestation and peatland cultivation and rainfall patterns. New estimates of fire emissions and their temporal pattern in Amazonia may help to understand longer-term fire dynamics, which we do here in the context of estimates of deforestation. We also anticipate that the record can be used to evaluate models that require information about fire and deforestation dynamics in this region, including biogeochemical, atmospheric, and Earth system models, especially those with prognostic fire models.

Horizontal Visibility Observations
We obtained raw, synoptic weather station records from the NOAA National Climatic Data Center Integrated Surface Database (https://catalog.data.gov/dataset/integrated-surface-global-hourly-data), which is composed of data from WMO-level surface stations provided by national meteorological agencies. While 55 WMO stations in our study region [7°S-17°S, 48°W-68°W] were available, many were excluded after applying our selection criteria as described below. Out of 55 stations 18 were found suitable and are depicted in Figure 2 and cover the areas with the highest total particulate matter (TPM) emissions based on the annually averaged Global Fire Emissions Database version 4 with small fires (GFED4s) over 1997-2014 (updated version of van der Werf et al. [2010]). Data filtering and station selection was based on three steps.
1. The individual observations were filtered and adjusted as in Husar et al. [2000] where observations that contained flags for rain or wet fog (relative humidity (RH) > 75%) were excluded, as were observations with RH > 90% or precipitation > 2.5 mm over the day. Furthermore, records with missing dew point or Global Biogeochemical Cycles 10.1002/2016GB005445 air temperature data were removed, because these are needed to calculate RH. Visibility records of 0 m were replaced by a set minimum of 50 m, the next lowest reported value, necessary for calculating extinction coefficients (described below), and we used a uniform maximum value of 20 km. Finally, to exclude years with missing landmarks or not enough observations to account for seasonality, we removed years with a variance less than 5 km. 2. Observations are often done at irregular time intervals, which vary per day, year, and station. For consistency we therefore aggregated the individual observations to one daily average. 3. The average daily observations were subsequently manually inspected for temporal inconsistencies and sudden jumps in the signal, which could be the result of, for example, a change in landmarks being used to estimate visibility. Horizontal visibility can be influenced not only by fire but also by various other sources, such as dust and fog, and is influenced by the observersˈ interpretation. This subjectivity introduces a random component in the signal. Also, when observations were too fragmented (for example, some stations were lacking observations during the fire season or only had continuous observations shorter than 1 year with gaps of several years) or when the maximum distance recorded was smaller than 12 km, stations were excluded [Husar et al., 2000].
To focus on low-visibility events and to correct for limitations of the human eye and imperfections of the landmarks, the observations in meters were then converted to the more often used extinction coefficient (B ext ) in km À1 using the Koschmeider relationship (equation (1)), where HV is the horizontal visibility in km and the Koshmeider constant (K) was set to a value of 1.9 [Husar et al., 2000;Schichtel et al., 2001].
We constructed our monthly record from 1973 to 2014, the period over which most stations had data albeit with gaps of various lengths of time ( Figure 3). We took the median B ext over the selected stations in the region [7°S-17°S, 48°W-68°W], the region with the highest fire emissions over 1997-2014 according to GFED4s. We also split this region into eastern and western halves to identify possible differences in fire activity and relationships to satellite data within the larger Amazon region. We compared the GFED TPM emissions with the aggregated B ext for the whole study region ( Figure 2) to understand how well they covaried and to get an indication of how B ext relates to TPM emissions. The comparison of B ext values with fire emissions was done for the dry season, thus the fire season, without taking atmospheric transport into account.
We also compared the visibility data set to total column CO observations retrieved from the MOPITT satellite as an independent constraint on fire emissions. CO is released from incomplete combustion of vegetation and fossil fuels, and fires substantially enhance the CO concentration over background conditions. CO observations from MOPITT have shown to be able to capture IAV in emissions in regions which are heavily impacted by biomass burning [Edwards et al., 2006;Logan et al., 2008;Mieville et al., 2010;Liu et al., 2011]. We used total column observations of CO for the months June-October from 2000 to 2014 in our comparison.
In situ surface measurements of particulate matter mass concentrations with an aerodynamic diameter smaller than 10 μm (PM 10 ) were also available for August 1992 to February 2005 in Alta Floresta (Brazil, State of Mato Grosso, 9.871°S, 56.104°W, 277 m above sea level), which provided a nonsatellite record for comparison, and extending further back than GFED or MOPITT [Artaxo et al., 1994;Maenhaut et al., 2002;Martin et al., 2010]. The particles sampled in Teflon filters were analyzed by gravimetry with filter explosion times ranging from 1 to 10 days. Summed uncertainties were estimated to be 15% of PM 10 mass concentrations, based on uncertainties in calibration of the mass flowmeters (5%) and elemental and mass analysis (10%, including filter handling).
In our approach, we assumed that potential changes in wind patterns are of second-order importance and have thus not been taken into account. This is supported by visual inspection of surface wind patterns derived from the National Centers for Environmental Prediction/National Center for Atmospheric Research (NCEP/NCAR) reanalysis [Kalnay et al., 1996]. From June to October the general wind pattern was relatively constant over our study period, with light winds from the east bending southward when reaching the Andes. Future work, using an atmospheric transport model and observed wind patterns, could lower the uncertainty arising from not taking variability in atmospheric transport into account.

The Role of Deforestation as a Driver of Fire Emissions
To understand the role of deforestation as a driver of fires and potentially also to yield information on how deforestation may have changed over time, we compared the visibility-based emissions with two forest loss data sets over the Brazilian Amazon and one data set of carbon emissions related to deforestation. We only used data sets which cover a relatively long time period from at least the 1990s onward.
First, we compared our results to net forest loss based on passive-microwave satellite observations of vegetation optical depth (VOD), which has shown to be sensitive to changes in biomass density [Liu et al., 2015]. The VOD-based deforestation data set we used here [van Marle et al., 2016] provides annual estimates of net forest loss, which is the net result of deforestation, forest degradation, and regrowth within a grid cell, on a 0.25°× 0.25°resolution for all of South America from 1990 to 2010. Interannual VOD changes over the region agreed reasonably with the Global Forest Change data set derived from Landsat, during the overlapping period since 2001 [Hansen et al., 2013]. VOD data for 1991 were missing, because the eruption of Mount Pinatubo influenced the passive microwave signal.
The second data set we compared our results with is the Program for Deforestation Assessment in the Brazilian Legal Amazon with Satellite Imagery (PRODES), which is developed by the Brazilian National Institute for Space Research (INPE). PRODES estimates annual deforestation events of primary forest larger than 6.25 ha based on a multidata approach utilizing Landsat data, supplemented by data from the China-Brazil Earth Resource Satellite and from the UK-DCM2 from the Disaster Monitoring Constellation International Imaging [Shimabukuro et al., 1998]. PRODES provides data from 1989 onward, although no data are available for 1992. In this study we used the summed total of the provinces of Acre, Mato Grosso, Rondônia, Pará, and Tocantins, the provinces within the Amazonia Legal closest to our study region (Figure 2), from 1989 to 2014. PRODES does not include Peru and Bolivia, which were part of our study region.
Finally, we compared our results, including decadal changes, with the updated country-scale carbon emission estimates associated with deforestation fires based on the bookkeeping method of Houghton [2003] using deforestation rates reported by the FRA 2015. For further details on the collection and validation of national data for FRA2015 see MacDicken [2015]. Deforestation is defined there as the conversion from forest to some other land cover. The rates are based on changes in forest area over a period of 5 years. In Houghton [2003] annual changes in cropland and pastures are used to provide year-to-year differences in deforestation rates. Deforestation is then assigned to forest types (i.e., biomass densities), and a fraction (35%) of that biomass is assumed to be burned [Houghton et al., 1991;Houghton, 2003]. The updated data set is not published yet. For comparisons with the western box and for the whole region we used the total of Bolivian and Brazilian country estimates. For the eastern box only the Brazilian estimates were used.
For the comparison on a decadal scale we converted the visibility-based results to carbon emissions by converting the B ext to TPM emissions based on their linear relationship found per 10°× 10°box and subsequently converting the TPM emissions to carbon emissions based on the emissions factors used within GFED (drawn largely from Akagi et al. [2011]). The converted visibility-based carbon emissions were then compared with the summed FAO-based carbon emissions for Bolivia and Brazil.
Forest loss is not occurring at the same location every year and spatial variability in deforestation and fire occurrence could thus impact the visibility signal. We calculated the mean weighted distance from pixels with VOD-based forest loss to the WMO stations which revealed that the average distance was relatively stable with 550 AE 43 km (1σ) between deforestation and observations for the western box and 614 AE 39 km (1σ) for the eastern box without a clear increasing or decreasing trend over time. We therefore assumed in this study that the effect of a changing deforestation front did not influence B ext substantially, but we do acknowledge that this is a source of uncertainty, especially for the 1970s and 1980s.

Results
4.1. Visibility as Fire Indicator Over the South American Burning Region Figure 4a shows the visibility-based B ext and GFED TPM fire emissions for the whole region ( Figure 2). These visibility-based results are averaged over the two boxes, because B ext behaved slightly different for the two grid boxes, see section 4.2. Over this large region, the B ext signal explained 47% of the variability in GFED TPM (  (1) (Table 1). Years where MOPITT CO shows relatively high values are 2000, 2005, and 2007, and the visibility-based estimate is relatively high in 2001. The correlation improves when we take only the western region into account as will be explained later.
Before the fire satellite era, beginning in 1997, visibility-based estimates indicate that fire emissions were relatively low until 1987 and 1988. Those 2 years were the first relatively high-fire years in the B ext record, followed by a few low-fire years after which fire emissions started to increase during the 1990s, with 1995 as peak year over the 1990s (Figure 4). After 1995 fire emissions remained elevated above the pre-1987 values, but with substantial IAV. During the 2000s there is, in general, a decreasing trend toward the end of the record although 2010 was the highest fire year in our data set.

Subregional Visibility-Fire Relationships
To illustrate how visibility relates to local aerosol load, Figure 5 shows the visibility from the Alta Floresta station marked with "AF" in Figure 2. Seasonality is evident in most years, with higher B ext values during the dry fire season and lower values during the wet season. We compared the visibility signal from Alta Floresta with ground-based measurements of PM 10 at the same station (15% uncertainty based on summed estimated uncertainties for sampling method and analysis), and the IAV is captured reasonably well with Pearson r 2 = 0.84 on a monthly time step ( Figure 5) and r 2 = 0.79 for annual summed values.
Separating the analysis between western ( Figure 6) and eastern (Figure 7) parts of the large burning regions, and for the dry season only, led to different correlations between B ext , TPM and CO (Table 1). Both the western and eastern boxes had a higher agreement for the dry months only, June to October (Figures 6b and 7b), with r 2 = 0.61 for the western box and r 2 = 0.40 for the eastern box, compared to all monthly observations (Figures 6a and 7a). The 5 months during the dry season account for 98% and 95% of all GFED TPM emissions for the western and eastern boxes, respectively. For the dry season, monthly averaged B ext values were converted to TPM emissions according to equations (2) and (3) for the western and eastern boxes, respectively. The negative intercept in the equations is due to nonfire impacts on visibility.
We tested how sensitive changing the slopes and intercepts were because the initial best fit resulted in negative values in the 1970s, which is not plausible. To circumvent this issue, we increased the intercepts with steps of 0.5E11 until it started to negatively impact the correlations (r 2 = 0.61 and r 2 = 0.40 for the western and eastern boxes, respectively) and mean squared errors and used the corresponding slopes and intercepts as used in equations (2) and (3) in our calculations back in time.
The western box (Figures 6d and 6e) shows equal or higher agreement compared to the total region for both GFED TPM and MOPITT CO, explaining 47% and 32%, respectively (Table 1). The IAV for the western box is similar to the IAV of the total region. This is not the case for the eastern box (Figures 7d and 7e). Although the monthly aggregated visibility observations do show a linear relationship (r 2 = 0.40) on a monthly time step (Figures 7a-7c), there are no significant relationships with GFED TPM emissions or MOPITT CO  In the discussion we touch upon potential reasons behind these differences between the boxes.

Deforestation as Driver of IAV in Fire Activity
We finally compared visibility-based B ext with forest loss data sets to understand to what degree the fire signal may be explained by deforestation. For this comparison we mainly focused on the whole study region. Figure 8a shows that annual VOD-based net forest loss in the study region explained 31% of the variability in the B ext signal ( Table 2). The high-fire years in 1995 and 2010 are captured in both data sets, whereas VOD shows higher deforestation rates from 2003 to 2006. Figure 8b shows that PRODES-based deforestation for the provinces of Acre, Mato Grosso, Pará, Rondônia, and Tocantins cannot explain the B ext signal. The IAV between both data sets partly agrees, especially during the 1990s, including the peak in 1995. The poor correlation is for a large part due to disagreement between PRODES and the visibility-based results on emissions from 2009 onward, including the peak in high-fire year 2010.
We compared the annual averaged B ext with the FAO-based carbon emissions related to deforestation from an updated version of Houghton [2003]. Correlations with B ext were poor (Figure 8c and Table 2). One of the most striking discrepancies is that the FAO-based carbon emissions related to forest fires from the 1970s to the early 1990s were higher than the visibility-based emissions. In addition, there are key differences in IAV (as was the case when comparing against the two deforestation data sets). For example, 1996 and 2000, years with little visibility-based emissions, show higher values for the FAO-based results. On the other hand, 2010 was a high-fire year, according to the visibility-based data set, which did not show up in the FAO-based data set.
Although the eastern region did not show any significant relations for the B ext signal with GFED TPM and MOPITT CO, VOD-based net forest loss explained 33%, whereas FAO-based carbon emissions did explain 18% of the B ext IAV for the eastern region. For the western region both VOD and Houghton explained, respectively, 23% and 9% of the B ext signal (Table 2).
We converted the visibility-based B ext signal to TPM emissions based on the linear relationships with GFED TPM mentioned before for the different regions and to carbon emissions using emission factors from GFED. We then compared the visibility-based carbon emissions with the summed FAO-based carbon emissions related to deforestation based on Houghton [2003] for Brazil and Bolivia at a decadal scale. The footprint of the two data sets differs because the FAO-based values are from all of Brazil and Bolivia, but our study region (Figure 2) covered the region with the highest fire emissions from 1997 to 2014 (Figure 2). Our estimates are therefore conservative. From the 1970s to the 1980s, the FAO-based data set and the visibility-based emissions show opposing trends, with a large increase (+747%) for the visibility-based carbon emissions and a decrease (À32%) for FAO-based deforestation-linked fire emissions (Figure 9). The large relative increase in the visibility data set is mostly because of low fire emissions in the 1970s. For both data sets the 1990s had the highest fire emissions. Both data sets also show a decrease from the 1990s to the 2000s, although the FAO-based emissions decreased more strongly (À41%) than the visibility-based emissions (À8%).

Visibility as Fire Indicator
On a local level, the good agreement (r 2 = 0.84 with monthly time step and r 2 = 0.79 with annual data) between visibility-based TPM and in situ observed PM 10 at Alta Floresta station ( Figure 5) indicated that visibility is strongly linked to aerosol loads and can be used, with limitations, as a proxy for fire emissions. This has been shown earlier to be the case as well in Indonesia [Field et al., 2009]. This also holds on larger scales; for the western box the pattern of the visibility-based signal agreed reasonable with the GFED-derived TPM emissions (r 2 = 0.61) and MOPITT CO concentrations (r 2 = 0.42) for the overlapping time period. For the whole Figure 8. Annual estimates of fire emissions and forest loss for the whole study region (Figure 2). In all plots dark green indicates the annual averaged visibility-based B ext , which are compared with (a) VOD-based net forest loss (1990-2010), (b) summed PRODES deforestation for the provinces of Acre, Mato Grosso, Rondônia, Pará, andTocantins (1989-2014, and averaged FAO emissions attributed to deforestation fires for Brazil and Bolivia. region (the western and eastern boxes combined) results were less good; IAV of B ext explained 47% of GFED TPM emissions and only 18% of MOPITT CO concentrations. The unexplained variability is likely due to measurement errors in the visibility signal and issues related to its spatial inhomogeneous network, the conversion to B ext , other sources that impact B ext , and the effects of atmospheric transport and chemistry. The error in atmospheric transport may be larger for the eastern box because variability in wind patterns is larger than in the western box, and smoke gets more concentrated in the western box because of the proximity to the Andes (Figure 1). Given the overall reasonable results, our findings may be suitable to fulfill the need for more detailed information on land use and land cover change over the tropics before the satellite era [Willcock et al., 2016].
Annual averaged B ext values, and thus fire emissions, were relatively low over the first part of the time series from 1973 until the first high-fire years in 1987 and 1988. These two relatively high-fire years were followed by some low-fire years, after which fire emissions started to increase over the 1990s, with 1995 being the peak year over the 1990s, which was also a high-fire year according to fire season severity predictions based on sea surface temperature anomalies [Chen et al., 2011]. After 1995 fire emissions were higher than the pre-1987 values. During the 2000s there was a general decrease, although 2010 was the highest fire year in our data set. This downward trend is also reported by the Global Forest Change project and satellite-observed burned area, suggesting that fires and deforestation are linked over that time period [Hansen et al., 2013;Fanin and van der Werf, 2015].
The first high-fire years 1987 and 1988 correspond to two El Niño years, which are associated with dryer conditions over the burning regions [Chen et al., 2011], although earlier in the 1970s no increase in fire emissions was observed during El Niño years. This suggests that fire in the region increased due to conversion of forest to other land uses, similar to what was observed in Indonesia [Field et al., 2009]. A key difference with Indonesia is the stronger control of El Niño-Southern Oscillation (ENSO) over drought and thus fire emissions. At a decadal scale the trend of increasing fires from 1970 to the 1990s corresponds with the results based on country statistics of Mouillot and Field [2005] for South America. However, the timing of the increase is different. Mouillot and Field [2005] suggested that the strongest increase was during the 1980s, while our data indicated that this was a decade later in the 1990s. There is thus apparently a difference between what was reported by countries and what is actually seen in the data, assuming that the average distance of deforestation to the stations before the 1990s showed equally little variability as after the 1990s. Another reason could be that fire was used to a smaller extent in the early part of our analysis to convert forest to other land uses.
The correlations between B ext and other data sets we found for the whole region were negatively impacted by the eastern box, which showed for both the GFED TPM and MOPITT CO comparisons lower agreement than the western box (Table 1). The most likely reason for this is the geographical location of the western box adjacent to the Andes. The prevailing wind during the dry season is from the northeast [Killeen et al., 2007], and thus, the smoke will be concentrated in front of the Andes mountain range (Figure 1). This implies that the visibility signal of the western box is probably representative for a much larger region. We tested this hypothesis by comparing the visibility signal for the western box only with the GFED TPM emissions for the whole region (i.e., the summed values for the two 10°× 10°from Figure 2, 7°S-17°S, 48°W-68°W). We found that 65% of the GFED TPM time series could be explained by B ext , which is 18% more than the comparison of the visibility signal for the western region, with the GFED TPM emissions for the western region (r 2 = 0.47, Table 1). Comparing MOPITT CO concentrations and VOD net forest loss over the whole region with the B ext for the western box only also revealed that B ext could explain more of the interannual variability with an increase of 8% and 7% for correlations between B ext with MOPITT CO concentrations and B ext with VOD-based net forest loss, respectively. This also explains why the combined visibility observations of the eastern and western boxes still explained 47% of GFED TPM emissions and 18% of MOPITT CO concentrations (Table 1 and Figure 4). Future research could investigate this in more detail including regional transport and atmospheric chemistry modeling.
For the western box, the overlapping time period of visibility with GFED (1997GFED ( -2014 showed for most years an agreement in IAV, including the high-fire years in 1999in , 2004in . However, in 2001in , 2005 the GFED and visibility-based fire signal disagreed, potentially due to uncertainties in both data sets or changes in atmospheric transport. A potential error in GFED is related to assigning the burned area to the wrong fuel classes. For example, if most of the burned area occurs in open landscapes but is assigned to deforestation fires, then emissions are overestimated. Visibility observations are also prone to errors, due to the nature of this human-observed proxy and our approach of combining the different records. Furthermore, the WMO stations are not evenly distributed over the region, which could have the effect that fires that occurred close to a station will have a relatively large effect on the average signal and that some fires may be missed, because there is no station located near or downwind from those fires.

Deforestation as Driver of Fire Activity
To answer the question to what extent deforestation can explain trends and IAV in fire emissions, we compared B ext with several deforestation data sets. For the whole study region VOD-based forest loss explained 31% of the fire signal, which suggests that deforestation and fires were coupled to some degree in this region for the overlapping time period from 1990 to 2014. PRODES deforestation explained little of the fire signal. We used the summed values of Acre, Mato Grosso, Rondônia, Pará, and Tocantins for PRODES. These are the provinces within the Legal Amazon closest to the two boxes, and results did not improve when we used PRODES data for the whole Legal Amazon or decreased our subset to a smaller number of provinces. One potential reason is that PRODES does not capture forest degradation, which caused high fire emissions in, for example, the 2010 dry year [Fanin and van der Werf, 2015]. Furthermore, we assumed that the state statistics for the five Brazilian provinces used were representative for the western box, because we expected the source regions to be largely from these provinces as a result of the prevailing wind direction. To some degree the comparison is not fair because PRODES is limited to Brazil, while the visibility signal potentially will also be affected by dynamics in Bolivia and Peru. Future research could include specific source regions, thus providing more insight into the actual area that is represented by the visibility signal.
Although correlations were higher, the comparisons with VOD-based net forest loss also suffered from limitations in the data sets. This forest loss data set captures mainly large-scale forest losses and in areas with many small-scale degradation events, VOD-based forest loss tends to overestimate, which could explain why the VOD-based interannual variability provides overall much higher values in 2003 to 2006 compared to the visibility-based fire emissions. Also, regrowth within a pixel is observed which can offset some of the signal [van Marle et al., 2016].
Overall, the comparison with the VOD-based deforestation data set provided evidence that deforestation is what at least partly drives these fires (r 2 = 0.23 and r 2 = 0.33 for the western box and eastern box, respectively, and r 2 = 0.31 for the whole region), although fires are a complex mix of deforestation, degradation, and natural causes. After deforestation, the land cover or use may change leading initially to an increase in burned area to maintain an open landscape or remove leftover debris. This is not captured in the deforestation data but will impact the visibility signal although fuel loads are, in general, lower in these types of fires [van Leeuwen et al., 2014]. Overall, it is difficult to explain the fire pattern without taking deforestation into account supporting Reddington et al. [2015], who concluded that deforestation fires dominate regional air quality  Morton et al. [2008], who showed that most fires in the Amazon during the 2000s were related to conversion of forest to agricultural land. We therefore cautiously conclude that the forest to other land use conversion started slowly at the end of the 1980s and increased strongly during the 1990s with accompanying deforestation fires. We are cautious, because the drivers of deforestation could have changed over time, influenced by changes in infrastructure, logging, and governmental subsidies [Moran, 1993;Fearnside, 2005;Rudel, 2007;Cochrane, 2009], and not all deforestation includes fire. An analysis of the GFED time series showed that the ratio between deforestation and nondeforestation fire emissions is relatively stable though, with on average 55% of fire emissions stemming from deforestation fires (1997)(1998)(1999)(2000)(2001)(2002)(2003)(2004)(2005)(2006)(2007)(2008)(2009)(2010)(2011)(2012)(2013)(2014)(2015).
Finally, the FAO-based deforestation emissions could explain 9% of the fire signal in the western box and 18% of the fire signal in the eastern box over 1973-2014. There was less agreement between B ext and FAO than between B ext and VOD, and this is probably the result of what VOD and FAO observe and how IAV is captured. VOD observes the net change in biomass, whereas FAO FRA only observes the conversion from forest to another land use type. For example, changes in biomass due to understory fires or from fires that do not lead to land use change may be detected by VOD, but not by PRODES and FAO; a potential explanation of why the peak in 2010 is missed. Furthermore, the IAV of carbon emissions in the bookkeeping method of Houghton [2003] is based on 5 year estimates of FAO FRA in combination with annual changes in cropland and pastures, thus not fully capturing IAV. The third reason for differences is the scale on which deforestation is observed. For most comparisons we found a better agreement with visibility and other data sets for the western box, compared to the eastern box. However, the FAO-based values compared better for the eastern box. Since the deforestation estimates are on a country level and not separated by province, this could be the result of discontinuity between Brazilian and Bolivian data. Although the two boxes we selected include the region with most fires over 1997-2014, this does not mean that all IAV for Brazil and Bolivia is captured.
This potential mismatch in footprint obviously also impacted the results on a decadal scale. The biggest difference between the FAO-based and our visibility-based fire emission estimates was in the early part of our time series from the 1970s until the first half of the 1990s. The visibility data indicated a gradual increase from very low to moderate emissions, while FAO-based emissions changed from relatively high to moderate. FAO data are based on the country totals of Brazil and Bolivia and therefore also captured carbon emissions related to deforestation occurring outside the Amazon region. Both data sets do agree on decadal averages for the 1990s and 2000s.

Conclusions
We have reconstructed fire emissions with an annual time step from 1973 to 2014 based on visibility observations from weather stations across the Amazon. Our study region encompassed the Arc of Deforestation and Bolivia, and we compared and constrained our results with satellite-derived fire emission estimates for the overlapping 1997-2014 period and MOPITT CO concentrations over the 2000-2014 period. These comparisons indicated that, with limitations, visibility can be a useful proxy for fire emissions in the Amazon. In general, the signal improves towards the western part of the Amazon where smoke gets more concentrated, because the Andes blocks the overall eastern winds.
The visibility signal showed that there was relatively little fire activity until 1987, after which fires increased and plateaued in the 1990s and 2000s before leveling off more recently. Substantial interannual variability has been present since the late 1980s and was partly related to ENSO, and comparisons with deforestation data sets showed that fires and deforestation were linked here although not all deforestation data sets were suitable for direct comparison with the visibility signal. Our estimates point towards lower (fire-driven) deforestation in the 1970s than found in FAO-based data sets. In our work we have assumed that the variability in B ext was directly proportional to emissions, omitting changes due to variability in atmospheric transport as well as changes due to spatial variability in the location of fire and deforestation events. Initial estimates indicated that variability in wind patterns was small as was the average distance of fire to the stations. While the errors due to this seem relatively small, a more extensive and spatially explicit study including a regional transport model is needed to better understand what impact progressing fire and deforestation fronts had on the observations. Global Biogeochemical Cycles