From spatio-temporal data to chronological networks: An application to wildfire analysis

12/04/2018 ∙ by Didier A. Vega-Oliveros, et al. ∙ Universidade de São Paulo Unife 0

Network theory has established itself as an important tool for complex systems analysis and pattern recognition. In the context of spatiotemporal data analysis, correlation networks are used in the vast majority of works. However, the Pearson correlation coefficient captures only linear relationships and does not properly capture recurrent events. This missed information is very important for temporal pattern recognition. In this work, we propose a chronological network construction process that is capable of capturing various events. Similar to the previous methods, we divide the area of study into grid cells and represent them by nodes. In our approach, links are established if two consecutive events occur in two different nodes. Our method is computationally efficient, adaptable to different time windows and can be applied to any spatiotemporal data set. As a proof-of-concept, we evaluated the proposed approach by constructing chronological networks from the MODIS dataset for fire events in the Amazon basin. We explored two data analytic approaches: one static and another temporal. The results showed some activity patterns on the fire events and a displacement phenomenon over the year. The validity of the analyses in this application indicates that our data modeling approach is very promising for geographical systems mining.



There are no comments yet.


page 4

page 5

page 6

page 7

This week in AI

Get the week's most popular data science and artificial intelligence research sent straight to your inbox every Saturday.

1. Introduction

Every second, a vast amount of spatio-temporal data are produced in the whole world. Examples include phone calls, live concerts, crimes, and car accidents. The data collected from these events are usually used to understand the behavior and to get much other insightful information from the complex systems they are part of. In this sense, Network Sciences have provided several resources for this purpose. The network representation allows the study of the interactions and dynamics between the small components from the represented complex systems in a unified way.

A common approach to represent geographical data into networks is based on the construction process of linking nodes according to their correlation coefficients, which are calculated from the underlying time series for each point of the spatial grid. This approach has been successfully used in a wide range of scientific areas. For instance, in Earth Sciences, networks have been used to analyze global climate (Zhou et al., 2015), to predict El-Niño, and explore its impact around the world (Tsonis and Swanson, 2008; Meng et al., 2018; Fan et al., 2017). In Bioinformatics, networks were applied to study gene expression (Farkas et al., 2003). In Finance, it has been used to understand and find the dynamics of financial markets (Bialonski et al., 2011).

Network science has helped to identify valuable information in many domains. However, several questions have raised around the limitations and applications of correlation networks. For example, what is the minimum time series length to consider the correlation? Is this long-length enough to find statistically significant correlations? What is the minimum correlation threshold to connect two nodes? Which temporal patterns can and cannot be captured by this construction process? Notwithstanding having answers to all of these questions, the correlation-based networks are not appropriate to many real-world spatio-temporal systems. For instance, if we only have short-length time series, how can we get a precise correlation coefficient between those time series? How can we measure how much the system changed between short time intervals?

To tackle the above questions, we propose an alternative approach to overcome those constraints: The chronological network construction. Similar to the previous works, our method represents grid cells of a geographical region by nodes. However, we connect them in a different way. A link is created between nodes if two consecutive events occur between them. If those events are in the same grid point, the nodes form a self-loop. Although being easy to build and use, very few studies have explored the potential of this network construction process to model and analyze spatio-temporal data sets. In this paper, we apply our new method to study the Amazon Basin fire activity. The Amazon, a vast region in South-America, is covered by rain-forest that contains a colossal biodiversity (Vieira et al., 2008). Unfortunately, over the recent decades, the region has been directly affected by deforestation, which is influenced by several factors like urban growth, farming areas, cattle industry, roadworks, fires, among others (Vieira et al., 2008; Zemp et al., 2017). The Amazon Basin is the largest drainage basin in the world, discharging about 209,000 cubic meters per second to the ocean (Vieira et al., 2008; Zemp et al., 2017)

. With that quantity of water flowing through the Amazon region, it seems that fire propagation in large areas, as occurs in other parts of the world, is an unlikely event. However, looking at fire data sets collected from satellite, wildfires in the Amazon are very dense and active in different regions. Although those fires events are not directly related, as a complex system, they share conditions that increase the probability of happening in some areas, even when they are far between each other.

Wildfire has a considerable impact on human life. This environmental process is responsible for vegetation composition changes (Brando et al., 2014), emission of gases and particles into the atmosphere (Gatti et al., 2014), and damage to properties (Dey and Schweitzer, 2018). Furthermore, climate change has been changing the frequency and severity of fire (Rogers, 2018). The many factors that influence fire activity make this complex system hard to forecast. In recent years, complex networks emerged as a powerful tool to study systems like this. Since fires dynamics is still poorly studied using network science, we believe that this application domain suits very well for our method. In summary, our goal in this paper is to use our network construction method to model and study wildfires focusing on the advantages and limitations of our method.

2. Related work

Concerning spatio-temporal events, Abe and Suzuki (Abe and Suzuki, 2006) employed a similar approach called sequential networks, for describing and finding patterns on the earthquake network of United States of America and Japan. They also analyzed the topological properties on growing networks. Zemp et al. (Zemp et al., 2017) model the moisture recycling process in South America and developed a similar framework to event-based networks with weighted nodes and directed links. More recently, Ferreira et al. (Ferreira et al., 2018) reported a variation the work in  (Abe and Suzuki, 2006) introducing window times and new rules to link nodes. As results, the authors obtained a better understanding of long-range seismic activities. The before works are interesting due to the simplicity to build the data model network and to find connections between different regions, even from the world (Ferreira et al., 2018). Therefore, the following sections will introduce the proposed methods to explore the advantages of network characterization based on the chronological approach, with potential applications in several fields.

3. Background

Let consider a geographical data with a set of spatial points. They can be represented by the network , where is the set of spatial located points called as nodes or vertices , and the set of edges or links denoting some similarity or distance relationship. The adjacency matrix is the mathematical entity of the network, where means there exists a connection between nodes and , and , otherwise. The links can be directed, indicating the initial and target node of the relationship, or undirected where both have the same relationship. Also, the elements of can be weighted, meaning the similarity strength between the nodes.

The degree of connectivity of node , called as , is the number of links incident on . When the links are weighted, the alternative is to calculate the strength of the nodes, which is the sum of the weights of incident links of each node. In the case of directed networks, is the sum of the degrees of input (links that reach the node) and output (links that leave the node). The degree distribution of a network is the probability of randomly select a node with degree . The level of disorder or heterogeneity of nodes connections is obtained with the entropy of the degree distribution, calculated by the normalized Shannon entropy (Wang et al., 2006), i.e.,


with . The more heterogeneous the distribution the greater the entropy, resulting in with a uniform and when all vertices have the same degree. The entropy of a network is related to the robustness and level of resilience (Wang et al., 2006).

The network decomposition into shells or K-Core (Seidman, 1983) (K or core of order K) is the maximum subset of nodes that have at least degree K and the K-Core is the highest-order core they belong to. Formally, a node is in a core of order K belongs to the K-Core but not the (K+1)-Core decomposition (Seidman, 1983; Kitsak et al., 2010). Vertices with the highest coreness are the most central. On the other hand, communities are sets of densely interconnected vertices and sparsely connected with the rest of the network (Newman, 2010). Nodes that belong to the same community, in general, share common properties and perform similar roles. Therefore, the division of a network into communities helps to understand their topological structure (structural and functional properties) and its dynamical processes, obtaining relevant information and features to the network domain. Several methods have been reported to detect communities on networks (Fortunato and Hric, 2016). Two methods adopted here are the agglomerative and optimization fastgreedy algorithm (Newman, 2004), and the fast modularity optimization method (Lambiotte et al., 2014).

4. Exploring spatiotemporal patterns in complex networks

In a geographical complex system we have the set of spatial-located points producing some signal or information over time. This information can be represented as chronological events and employed for discovering global, local and intermediate patterns of the studied system. We presented our approach of grouping this points in grid cells and constructing an event-based characteristic network. The construction of this complex network, associated with the geographical system, can be conducted for the whole period of time or in fixed/dynamic intervals following a multigraph or temporal networks technique. Next, we present our event-based modeling network and two data analytic approaches with results and interpretations related to the fire-event Geographical system.

4.1. Data

In a global scale, the fire event activity is collected mainly through satellite instruments like the Moderate Resolution Imaging Spectroradiometer (MODIS) and Visible Infrared Imaging Radiometer Suite (VIIRS). The MODIS runs in both, Aqua and Terra satellites, operated by the National Aeronautics and Space Administration (NASA), while the VIIRS works at Suomi-NPP satellite, also by the NASA and the National Oceanic and Atmospheric Administration (NOAA). In this research, we employ the MODIS data due to the more than 15 years of available fire event around the world and the data confidence for the continuous refinements and calibrations performed in the MODIS system.

Here, we perform the research using data from the last version (C6) of MODIS. The time interval is between 01 January 2003 and 31 January 2018 and the region under study is a portion of the Amazon basin, that is located between longitude W, W and latitude S, N. From the MODIS data, we consider the UTC date and hour, both satellites Terra and Aqua, the geographic coordinates, and the detection confidence of the fire event. The total number of occurrences is 1,684,600, after a filtering process considering detection confidence above 70%. It is relevant to notice that the satellites sequentially scan the earth surface, i.e., the capture is by resolution points. Therefore, there are not parallel events at the same time in the data set.

4.2. Event-based Characteristic Network

Our network characterization process for spatio-temporal events is based on these three steps:

  1. Grid-division: A geographical region under consideration is divided in a grid. Each grid cell is represented as a node in the network.

  2. Time length: The network data modeling can be defined for specific periods of time, e.g., the whole time period, fixed or dynamic intervals.

  3. Links construction: From the data set, two successive events create a link between the grid cells where they are located.

In Figure 1, we show a representation of the network construction method and the tackled problem. Every event has a different timestamp, even that wildfires in different geographic areas can occur at the same time (Figure 1.A). The before is because fire events in the MODIS data are sequential. However, in the case of parallel occurrences, the construction process can proceed branching the connections. The spatio-temporal data is represented in a grid-division, and the events are linked as they occur (Figure 1.B). Then, in Figure 1.C, we have the chronological network representation of the data. Our approach has low code complexity and linear computational cost , with being the number of processed geographical events. Furthermore, the faster network construction allows a real-time data streams ingestion. To make reproducibility easier, we share the implementation of our method online111Code available at

Figure 1. Illustrative example of a set of spatio-temporal events transformed into a network by following the chronological characterization.

Some points concerning this method are raised here. First, what is the optimal grid-division? It depends on the quantity of data and the spatial distribution, in order to avoid dense networks – where almost every node is connected to each other (a feasible solution is to increase the grid-division or to create time slices), and also to avoid sparse networks (a solution could be to decrease the grid-division or to increase the time length). Second, what happens if we have self-loops or multiple links between two nodes? For instance, Figure 1.C we have a self-loop in node . After the linking process, the network can be simplified by removing self-loops, the directions and multiple links. Each type of network can reveal different properties of the characterized system. However, following different networks analysis can lead to similar results (see discussion in Section 5).

According to the given considerations, we analyze two different network approaches to show the contributions of our chronological method: (I) Building a single network with multiple directed links and no self-loops for the whole period (Section 4.3). Thus, we analyze the long-term properties of the system. (II) Building unweighted and undirected temporal networks from the whole period. Thus, we explore the geographical system evolution over the 15 years (Section 4.4). The best parameter combination (grid size and time length) was found by performing a sensibility analysis, seeking to maintain a consistent links density after the simplification process. As illustrated in Figure 2, an optimal grid-division for the Amazon basin area is , and for the temporal division, periods of seven days. The reason is that at this points we achieve a link density plateau with almost of the original edges after the simplification.

Figure 2. Sensibility analysis by measuring the links density for grid size division and days length values. The darker the color, the higher the links density.
Figure 3. Network representing the whole data set. The small dots are the nodes (grid cells) and red lines represent the links (consecutive events). The background colors represent different land uses. The line transparency is inversely proportional to the link weight. This network is characterized by six regions with fire activity: (A) east Colombia (B) Roraíma state in Brazil, (C) Park of Tumucumaque in north Pará state in Brazil, (D) Amapá state in Brazil, (E) Amazon river banks, and (F) southeast Amazon forest region.

4.3. Single Network Approach

In this section, we use our method (Sec. 4.2) to construct a single network and study the historical data (Sec. 4.1). The goal is to use this network to characterize the whole data set and search for patterns. We start by constructing a directed and weighted network. The link weights represent the frequency of fire occurrences between two nodes on the data set, and the direction is the chronology flow. Then, we simplify it into an undirected and weighted network without self-loops, by adding the weights (in- and out-links) between two nodes, i.e., . Figure 3 illustrates the resulting network that can be divided into six regions of fire activity. For the sake of simplicity, we will refer to this network as “historical network”.

Figure 4. Strength and degree (inset) cumulative distributions. Both distributions have decay sharper than a power-law.
Figure 5. Nodes strength for the historical network. The sizes of the nodes (red dots) are proportional to the strength.
Figure 6. Community structure found with the fast greedy algorithm (weighted) (Newman, 2004). The colors represent the four communities achieved by the highest modularity. In our network construction method, a community represents a region of frequent fire activity that occur in a period of time.
Figure 7. Fire events over time.

Our first analysis consists on studying two simple network properties: node degree and node strength. The degree of a node here represents the number of different areas whose fire events occur consecutively with the grid cell of . The strength measures the frequency that a fire event occurs consecutively between all the neighbors of . Since the degree does not account for frequency, the strength measure seems to bring more information about the networks generated by our method. In Fig 4 we present the cumulative strength and degree (inset) distributions for the historical network. Both distributions have a decay faster than a power-law, what suggests the network is not scale-free. These distributions show that many nodes have low degree and low strength. Conversely, the network have just a few nodes with high degree and/or strength. It is important to note that a high degree does not necessarily imply a high strength. But, in fact, the weights and degrees of the historical network are highly correlated (, p-value ).

Figure 8. Normalized entropy of the weekly temporal networks

Since the historical network have just a few nodes with high strength, one next step is to analyze the geographical locations of these nodes. In Figure 5 we illustrate the node strengths for the historical network. The regions with high frequency of consecutive events are similar to those ones found in Figure 3. However, here it is possible find what are the most frequent sub-regions. In fact, these sub-regions are related to the land use (background colors in Figure 5). As pointed by a recent research (Curtis et al., 2018), the land activity in these regions in the last years are mainly characterized by commodity driven deforestation (agriculture, mining, etc) and shifting agriculture (that is later abandoned).

4.4. Temporal network approach

In the previous section, we explore an approach of constructing a single network from the entire historical period. We observe general patterns related to some macro-regions of fire activity, recognized by a community detection algorithm. In this section, we face the problem form the point of view of multiple layers or temporal networks approach. We generate networks per weeks from the fire events data by employing our proposed construction method. In this way, a set of layers, or temporal networks, represents the spatial fire activities and patterns during each week. The set of networks is constructed in consecutive intervals of days, beginning from Jan. 01, 2003 to Jan. 24, 2018. Formally, with layers, where is the network of the first days, the next days, and so on.

We start the temporal analysis considering intervals of days and the frequency of fire events in the grid cells. At first glance in Figure 7(a), we can assume there is a pattern of fire season in the Amazon basin. However, this pattern is not entirely clear over the years and the start and end seasons are not so well defined. As an example, the marked interval from Aug. 25, to Oct. 06, 2010 in Figure 7(a) seems to be at the end of the fire season of 2010. However, in Figure 7(b) we observe there was a high activity of fire events in the basin during this interval of time. Motivated by the lack of precision of the direct frequency approach, we analyze the spatio-temporal fire event data considering weekly temporal networks. For this purpose, we employ the proposed event-occurrence construction method to mine the activity patterns in global (networks), local (nodes) and intermediate (communities) scale.

(a) Fire frequency
(b) Degree centrality
(c) K-Core centrality
Figure 9. (Color online) Comparison of the temporal series values of each node (grid cells) of the studied area per weeks: (a) the fire frequency, (b) the node degree and (c) the K-Core value in the respective weekly temporal network. Higher the frequency or centrality, redder the cell.

In a global scale analysis, we explore some general measures to characterize the temporal networks. We calculate the normalized entropy (Eq. 1) for each network from , as shown in Figure 8, in which higher the entropy, higher the fire events activity and more heterogeneous the network. We can observe a clear pattern of fire activity in the studied region, starting at the beginning of winter and finishing in the middle of summer, which is, in general, the predominant dry season. Different to Figure 7(a) in the marked interval, the fire activity is in the middle of the peak in Figure 8. This result shows the advantages of considering the proposed data modeling method and network mining for obtaining extra knowledge.

Concerning the local scale, we process the information contained in the grid-cells or nodes in the case of networks. The traditional approach is to analyze the activation frequency of the grid cells over time. Then, statistical measures like mean, moments, percentiles, are used to understand the cells distribution. On the other hand, we can calculate some centrality measures over the network for finding the role of each node in the dynamic evolution. Several topological measures describe the relevance of nodes according to structural and dynamical properties 

(Lü et al., 2016). Many have been employed in climatological problems, like in forecasting and prediction (Fan et al., 2017; Meng et al., 2018), disaster risk, and management (Santos et al., 2019), among others. In particular, the degree and K-Core have great prominence in the area. For instance, K-Core centrality presents a good agreement for identifying the essential nodes in different dynamics, like the most influential spreaders in diffusion process (Kitsak et al., 2010; Vega-Oliveros et al., 2015), the best target nodes for vaccination or marketing campaigns (Hébert-Dufresne et al., 2013; Vega-Oliveros et al., 2017), and other domains.

4.4.1. Centrality Measures over time

Figure 9 shows the temporal values per week of each node or grid cell of the Amazon basin. As defined in Section 4.2, we have a total of grid cells, where node is the leftmost cell at the bottom of the grid. We can observe a faint seasonal pattern in the frequency of fire events per weeks for each cell over time (Figure 9(a)). In particular, for id cells over (y-axis), it is difficult to see the frequency values, leading to wrongly understand that there is a non-significant fire dynamic in the central and north regions. Besides, for x-axis from week to , it seems the seasonal pattern tend to vanish. Opposite, when considering the temporal networks, the degree and K-Core centralities (Figure 9(b) and (c) respectively, with the darker color representing the highest centrality) capture well the seasonal fire pattern of each year. These centrality measures enable to characterize better the role of each node in the fire dynamic, e.g., nodes of the north region (above id y-axis) have a more clear fire activity in Figures 9(b) and (c) than Figure 9(a).

The K-Core measure revealed to be particularly interesting in the context of fire activities. Nodes with higher K-Core centrality are most central in the fire activity, i.e., they are the centroids of the fire events. Nodes in the south (from to ), are the central focus of fire activity in the basin. Also, Figure 9(c) illustrates a lag pattern in the starting and finishing fire activity between the southern regions (bottom in the Figure) and northern regions (top in the Figure). This result indicates a dynamic of movement or displacement of fire activities in the Amazon basin throughout the year. The fire displacement may occur due to favorable environmental conditions and land use in the regions. The fire season starts at the southern areas and advances to the northern areas through the year.

4.4.2. The Centrality-series Similarity Network

We define the centrality series of a node as the consecutive measurements of its centrality over time. Formally, for and the centrality value , with the network where is calculated, the time series is the set of sequential centrality values among the temporal networks, i.e., , and the number of temporal networks. describes the importance of the node in the geographical system according to its centrality evolution. Thus, we aim to understand how the nodes are related and are similar over time. The first step is to calculate the similarity among all the time series of centralities. Then, we generate the centrality-series similarity matrix , defined as follow:


Here, we adopt the similarity function as the Pearson correlation and the centrality as the K-Core measure. The full matrix brings the similarity of temporal K-Core series between pairs of nodes. For the sake of simplicity, we will refer to this centrality-series similarity matrix as “CSS matrix”. Next, we employ the CSS matrix of K-Core for obtaining the centrality-series similarity (CSS) network. The constructed network allows finding local and dense structures according to the temporal centrality similarity of the nodes. For this goal, we employ the well known NN construction method, where a hub in the data space is a hub in the NN network (Berton et al., 2018). The NN method consists in connect each node with its most similar or nearest neighbors. There exists a link between and if , where is the set formed by the nearest neighbors of , and the main diagonal of . The larger the , the denser the network.

Figure 10. (Color online) Community structure of the geographic fire events in the Amazon basin. The size of the grid cells represents the overall intensity of fire activity over the time. In the background, deforestation (2013-2017) is shown in red, with yellow point for 2018. Data from Instituto Nacional de Pesquisas Espaciais  (INPE, 2018).
Figure 11. Average monthly precipitation for the detected communities. Big symbols show the interval of intense dry period.

4.4.3. Clustering Patterns and Heterogeneous Regions of Fire Events

We generate the CSS network using the NN method with parameter . Then, we calculate the community structure of the CSS network by the modular community detection method (Lambiotte et al., 2014) (Figure 10). We obtain significant communities, disregarding all the grid cells where the MODIS never reported a fire event, i.e., K-Core equal to zero in all the period. The subregions were analyzed concerning the spatio-temporal patterns and the fire activity. In this way, our method satisfactorily recognizes subregions that present dynamical differences in dry periods and intensity of fire activity.

We compared the subregions concerning the average monthly precipitation of the CMAP data products from NOAA/ OAR/ ESRL PSD 222avaliable at

. This is a global precipitation data set of 30-years monthly analysis based on gauge observations, satellite estimates, and numerical models 

(Huffman et al., 1997). The precipitation curves of each subregion are depicted in Figure 11. The dry season starts with values lower than mm/day, and are marked in the curves. We find the following patterns concerning the precipitation: (i) Subregions 4-7 have dry seasons occurring approximately in the same period: 4 and 5 have less precipitation than 6 and 7, especially in the dry season. (ii) Subregions 2 and 3 have the longest and more severe dry season than 4-7; 1 has a similar pattern to 2-3, but less precipitation in total (values in the rainy season are lower). (iii) Subregions 8-11 have a rainy season starting later and better defined (fewer months with high precipitation). Consequently, the dry season ends later: The rainy season of 9 also starts later but its dry season is a little shorter than 8,10,11. (iv) Subregion 12 presents the least intense rainy season (lower maximum values) and starts later. Consequently, the dry season is also relatively long and least intense. Subregion 12 has different period of the dry season than all other regions.

Finally, about vegetation cover and land use, the subregions have the following characterization: (i) Subregions 1-4 are related to land use related to pastures in old areas of deforestation, with subregion 3 relatively old. Subregion 2 has land use more related to agriculture, pasture, and areas that were forests. (ii) Subregion 5 is related to a savanna region in Bolivia. Due to the high fire activity detected in 5, there must be land use with agriculture. (iii) Subregions 6-8,10 are related to deforestation in areas surrounding reserves and roadsides with more strong gradients of forests and recently deforested areas. Subregion 8 should be related to more recent deforestation, with pasture and agriculture. Subregion 6 is related to land use with deforestation near roads. (iv) Subregions 9 and 11 are areas near the Amazon River. They are more related to recent settlements where there must be deforestation, agriculture, and pasture. (v) Subregion 12 is a savanna region surrounded by forest with deforestation and land use. It is a combination of natural propensity to fire with pasture maintenance.

5. Final remarks and Future works

In this paper, we explored the fire activity in the Amazon basin as a complex system. We presented a data modeling method for constructing chronological networks (Sec. 4.2) and two graph mining approaches (Sec. 4.3 and 4.4). We observed that these approaches captured important insights hidden in the fire data that other network-based models are incapable of highlighting. As results, we found subregions – or communities – with particular differences in the temporal (by its dry season periods), activity (by the intensity of fire events) and space (by the land-use and location) conditions for the occurrence of fire. Moreover, a pattern of fire-activity displacement over the year, starting from the southeastern subregions, spreading in the south and rising to the northwest of the basin was also identified. Also, although both approaches follow different ways to analyze the network, they recognize the same set of most active cells (Fig. 4 and Fig. 10). In summary, our method allows the study of spatio-temporal data sets from a different perspective.

For future works, we intend to employ and propose advanced techniques that could lead to more revelations on this kind of spatio-temporal problems. Evaluating more sophisticated rules, like linking nodes by considering a spatial threshold, is a natural path of this work for understanding how the new rules could affect the network structure. Another important point is that the same fixed time window does not necessarily space the events. Adding a memory-mechanism that avoid the time slicing approach, but allows to measure changes over time, could be a new strategy to understand temporal changes in complex systems.

This research is supported by the Sponsor Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP) Rl under Grant No.: Grant #3 and the Sponsor German Research Council (DFG-GRTK) Rl Grant No.: Grant #3. D.A.V.O acknowledges FAPESP (Grants Grant #3, Grant #3, and Grant #3). L. N. F. thanks FAPESP under Grant No.: Grant #3 and Sponsor Brazilian Higher Education Funding Council (CAPES) Rl scholarship.


  • (1)
  • Abe and Suzuki (2006) S. Abe and N. Suzuki. 2006. Complex-network description of seismicity. Nonlinear Processes in Geophysics 13, 2 (may 2006), 145–150.
  • Berton et al. (2018) L. Berton, A. de Andrade Lopes, and D. A. Vega-Oliveros. 2018.

    A Comparison of Graph Construction Methods for Semi-Supervised Learning. In

    2018 International Joint Conference on Neural Networks (IJCNN)

    . 1–8.
  • Bialonski et al. (2011) Stephan Bialonski, Martin Wendler, and Klaus Lehnertz. 2011. Unraveling Spurious Properties of Interaction Networks with Tailored Random Networks. PLoS ONE 6, 8 (aug 2011), e22826.
  • Brando et al. (2014) P. M. Brando, J. K Balch, D. C Nepstad, D. C Morton, F. E. Putz, M. T Coe, D. Silvério, M. N. Macedo, E. A. Davidson, C. C. Nóbrega, A. Alencar, and B. S. Soares-Filho. 2014. Abrupt increases in Amazonian tree mortality due to drought-fire interactions. Proceedings of the National Academy of Sciences of the United States of America 111, 17 (apr 2014), 6347–52.
  • Curtis et al. (2018) Philip G. Curtis, Christy M. Slay, Nancy L. Harris, Alexandra Tyukavina, and Matthew C. Hansen. 2018. Classifying drivers of global forest loss. Science 361, 6407 (2018), 1108–1111.
  • Dey and Schweitzer (2018) Daniel Dey and Callie Schweitzer. 2018. A Review on the Dynamics of Prescribed Fire, Tree Mortality, and Injury in Managing Oak Natural Communities to Minimize Economic Loss in North America. Forests 9, 8 (jul 2018), 461.
  • Fan et al. (2017) Jingfang Fan, Jun Meng, Yosef Ashkenazy, Shlomo Havlin, and Hans Joachim Schellnhuber. 2017. Network analysis reveals strongly localized impacts of El Niño. Proceedings of the National Academy of Sciences of the United States of America 114, 29 (jul 2017), 7543–7548.
  • Farkas et al. (2003) I. Farkas, H. Jeong, T. Vicsek, A.-L. Barabási, and Z.N. Oltvai. 2003. The topology of the transcription regulatory network in the yeast, Saccharomyces cerevisiae. Physica A: Statistical Mechanics and its Applications 318, 3-4 (feb 2003), 601–612.
  • Ferreira et al. (2018) Douglas Ferreira, Jennifer Ribeiro, Andrés Papa, and Ronaldo Menezes. 2018. Towards evidence of long-range correlations in shallow seismic activities. EPL 121, 5 (mar 2018), 58003.
  • Fortunato and Hric (2016) Santo Fortunato and Darko Hric. 2016. Community detection in networks: A user guide. Physics Reports 659 (2016), 1–44.
  • Gatti et al. (2014) L. V. Gatti, M. Gloor, J. B. Miller, C. E. Doughty, Y. Malhi, L. G. Domingues, L. S. Basso, A. Martinewski, C. S. C. Correia, V. F. Borges, S. Freitas, R. Braz, L. O. Anderson, H. Rocha, J. Grace, O. L. Phillips, and J. Lloyd. 2014. Drought sensitivity of Amazonian carbon balance revealed by atmospheric measurements. Nature 506, 7486 (feb 2014), 76–80.
  • Hébert-Dufresne et al. (2013) Laurent Hébert-Dufresne, Antoine Allard, Jean-Gabriel Young, and Louis J Dubé. 2013. Global efficiency of local immunization on complex networks. Scientific reports 3 (jan 2013), 2171.
  • Huffman et al. (1997) George J. Huffman, Robert F. Adler, Philip Arkin, Alfred Chang, Ralph Ferraro, Arnold Gruber, John Janowiak, Alan McNab, Bruno Rudolf, and Udo Schneider. 1997. The Global Precipitation Climatology Project (GPCP) Combined Precipitation Dataset. Bulletin of the American Meteorological Society 78, 1 (1997), 5–20.
  • INPE (2018) INPE. 2018. Projeto PRODES: Monitoramento da Floresta Amazônica Brasileira por Satélite. Instituto Nacional de Pesquisas Espaciais, São José dos Campos, Brazil, retrieved from (2018).
  • Kitsak et al. (2010) M. Kitsak, L.K. Gallos, S. Havlin, F. Liljeros, L. Muchnik, H.E. Stanley, and A. Makse. 2010. Identification of Influential Spreaders in Complex Networks. Nature Physics 6, 11 (2010), 888–893.
  • Lambiotte et al. (2014) R. Lambiotte, J. Delvenne, and M. Barahona. 2014. Random Walks, Markov Processes and the Multiscale Modular Organization of Complex Networks. IEEE Transactions on Network Science and Engineering 1, 2 (July 2014), 76–90.
  • Lü et al. (2016) Linyuan Lü, Duanbing Chen, Xiao-Long Ren, Qian-Ming Zhang, Yi-Cheng Zhang, and Tao Zhou. 2016. Vital nodes identification in complex networks. Physics Reports 650 (2016), 1–63.
  • Meng et al. (2018) Jun Meng, Jingfang Fan, Yosef Ashkenazy, Armin Bunde, and Shlomo Havlin. 2018. Forecasting the magnitude and onset of El Niño based on climate network. New Journal of Physics 20, 4 (apr 2018), 043036.
  • Newman (2010) Mark Newman. 2010. Networks: An Introduction. Oxford University Press, Inc., New York, NY, USA.
  • Newman (2004) M E J Newman. 2004. Fast algorithm for detecting community structure in networks. Physical Review E 69, 3 (2004), 66133.
  • Rogers (2018) Adam Rogers. 2018. The Only Thing Fire Scientists Are Sure of: This Will Get Worse. WIRED (Aug, 1, 2018), retrieved from
  • Santos et al. (2019) Leonardo B. L. Santos, Luciana R. Londe, Tiago Carvalho, Daniel S. Menasché, and Didier A. Vega-Oliveros. 2019.

    About interfaces between Machine Learning, Complex Networks, Survivability Analysis and Disaster Risk Reduction

    Springer International Publishing, Chapter In press, 1–32.
  • Seidman (1983) S.B. Seidman. 1983. Network structure and minimum degree. Social Networks 5, 3 (1983).
  • Tsonis and Swanson (2008) Anastasios A. Tsonis and Kyle L. Swanson. 2008. Topology and predictability of El Niño and la Niña Networks. Physical Review Letters 100, 22 (jun 2008), 228502.
  • Vega-Oliveros et al. (2015) D. Vega-Oliveros, L. Berton, A. Lopes, and F. Rodrigues. 2015. Influence Maximization Based on the Least Influential Spreaders. In SocInf 2015, co-located with IJCAI 2015, Vol. 1398. 3–8.
  • Vega-Oliveros et al. (2017) Didier A Vega-Oliveros, Luciano da F Costa, and Francisco A Rodrigues. 2017. Rumor propagation with heterogeneous transmission in social networks. Journal of Statistical Mechanics: Theory and Experiment 2017, 2 (2017), 023401.
  • Vieira et al. (2008) ICG. Vieira, PM. Toledo, JMC. Silva, and H. Higuchi. 2008. Deforestation and threats to the biodiversity of Amazonia. Brazilian Journal of Biology 68, 4 suppl (nov 2008), 949–956.
  • Wang et al. (2006) Bing Wang, Huanwen Tang, Chonghui Guo, and Zhilong Xiu. 2006. Entropy optimization of scale-free networks’ robustness to random failures. Physica A: Statistical Mechanics and its Applications 363, 2 (2006), 591–596.
  • Zemp et al. (2017) D. C. Zemp, C.-F. Schleussner, H. M. J. Barbosa, and A. Rammig. 2017. Deforestation effects on Amazon forest resilience. Geophysical Research Letters 44, 12 (jun 2017), 6182–6190.
  • Zhou et al. (2015) Dong Zhou, Avi Gozolchiani, Yosef Ashkenazy, and Shlomo Havlin. 2015. Teleconnection Paths via Climate Network Direct Link Detection. Physical Review Letters 115, 26 (dec 2015), 268501.