Several studies show how the aftershocks observed on time scales of minutes to months after a mainshock can be triggered by static and dynamic stresses imparted by the stress change and the mainshock waves’ passage  .
Branching structures  , associated with earthquakes and their development over time, have proven to be a helpful instrument to obtain crucial information about the occurrence of future events, the interaction between the earthquakes they consist of and the local faults’ critical and hazard state.
The purpose of this study is to describe the analysis method of seismic sequences in order to locate the nodes of higher order seismic rods in the upper and lower branching structure. Further, by analyzing the seismic sequence and using the nodes in the upper branching structure, we will propose a novel method for identifying the zone within which a future strong earthquake’s epicenter may fall.
Through the analysis of the seismic rods in the upper branching structure, we will also introduce a sequence seismic hazard classification in order to define the warning signs during the seismicity’s time evolution of the analyzed area.
Finally, the joint analysis of the seismic rods in the upper branching structure allows us to obtain relevant information on possible interactions between different earthquakes on local and global scale, probably due to the effects of the changes caused by static and dynamic stresses.
The analysis of the seismic sequence structure is performed by using: 1) a datasets spatial coverage involving parts of a region or several countries (from a few tens to thousands of km), 2) a range of magnitude values of 2.0 - 10 m, 3), a range of depths of 1 - 50 km.
2.1. Upper and Lower Branching Structure
The sequences of higher order seismic rods in the upper and lower branching structures are patterns that allow the graphical evaluation of the earthquake preparation process’s magnitude and therefore contribute to outline the past and future seismicity evolution of an area. These are identified through the consideration that after an energy shock the magnitude values gradually decrease and increase in the upper and lower structure, respectively, up to the Trigger Point (TP) by following a hierarchization process.
The entire preparation time of the upper branching structure can be divided into three stages: 1) energy accumulation, 2) trigger, 3) energy release  .
Figure 1 shows the patterns of the upper and lower branching structures: as it can be seen, the first developmental stage (the largest) in the upper branching structure starts to develop after a strong earthquake (point source), while the subsequent ones are related to the decrease in the magnitude values over time. In the lower branching structure, however, a trend to increase the magnitude values over time can be observed.
Furthermore, we note a decrease in the recurrence rates or the number of events in the final portion of the branching structures and a common trigger point (TP) determining the transition from the energy accumulation to the energy release stage. The energy release stage develops within a shorter time frame compared to the accumulation stage. The distance between the TP and the mainshock is the warning time frame of an energetic seismic event. Usually, a first warning signal results from the convergence of a third order lower branch with a second order higher one (point A). For example, in the case of the Samoan earthquake on 29 September 2009, a second abnormal
Figure 1. Energy accumulation phase’s developmental stages in the upper and lower branching structure. The vertical axis shows the magnitude values, while the horizontal axis displays the number of events recorded in the period analyzed.
warning sign resulted from the convergence of a lower fourth order branch with a higher second order one (point B).
2.2. Upper Branching Structure Hazardousness
The classification process provides not only information on the upper branching structure’s developmental state, but also a more consistent hazard classification related to the upper branching structure, which can be represented with colors.
The branching structure’s hazard level is classified according to the seismic branches magnitude  , by establishing seven hazard classes ranging from extremely high to low. The different hazard levels were encoded with colors that let us know, in a simple and quick fashion, the critical level reached by the upper branching structure that is being developed in an area either on local or global scale.
We chose the magnitude as a criterion for the division into hazard classes in order to identify the branching structure hazardousness at any point in its evolution.
The branching structures classification process involves several steps, including the following:
・ Upper branching structure detection;
・ Identification of the magnitude values of the seismic branches nodes it consists of and of the energy release stage;
・ Determining the hazard class of each stage of the energy accumulation and release phases;
・ Determining the hazard class of the energy release last stage.
As it uses seven hazard classes, this method can be successfully adopted to determine the appropriate warnings during the energy release stage, which can lead to a strong earthquake.
Table 1 lists the different hazard classes of the upper branching structure per magnitude range, while Figure 2 shows the hazard classes, determined according to the magnitude values of the energy accumulation and release stages during the Samoan earthquake on 29 September 2009.
Initially, we note a gradual reduction of the upper branching structure’s hazard level up to the TP, thereafter the hazard rises until it reaches an extremely high value coinciding with the mainshock. The midpoint of the fifth order branch, which defines the alleged value of the expected magnitude  , enables us to estimate the energy release stage’s minimum hazard, which in this case was estimated to be high.
2.3. Epicentral Area Identification-Aimed Procedure
This procedure, which can be of help in identifying the epicentral area of a strong earthquake, consists in plotting on a map the nodes’ epicenters in the upper branching structure along with the trigger signals. Figure 3 illustrates the position of the nodes’ epicenters in the Japan’s upper branching structure (Table 2), whose release stage was characterized by two foreshocks measuring 5.7 Mw and 7.3 Mw magnitudes respectively, recorded on 06/12/2010 and 09/03/2011 and by the earthquake occurred on 11 March 2011, with a magnitude of 9.0 Mw.
If we observe carefully the map, we will see how the epicenters tend to group over time depending on their type, i.e. nodes’ epicenters in the energy accumulation and energy release stages and TP-related ones. One of the epicenters of the shocks which the TP or the energy release stage consists of is almost always close to the mainshock.
Therefore, in order to predict the mainshock epicenter, it is considered more appropriate to use the clusters including the shocks epicenters which the TP or the energy
Table 1. Hazard classes.
Table 2. Coordinates and magnitude of nodes and warning sign in the upper branching structure of Japan earthquake, 2011.
Figure 2. Samoan earthquake. Branching structure and hazard classes.
Figure 3. Japan’s seismic sequence. Yellow circles: energy accumulation nodes’ epicenters; red circles: TP epicenters; green circles: energy release stage epicenters; red star: mainshock epicenter.
Figure 4. Earthquake’s seismic sequence, Japan, 2011, and upper branching structure.
release stage consists of.
Figure 4 shows the seismic sequence and the upper branching structure that preceded the earthquake on 11 March 2011.
2.4. Global Scale Branching Structures
Several earthquakes do not occur as independent events: just as faults’ mutual interaction causes a change in the static stress, so too earthquakes can interfere through the dynamic variation of the stress field. The co-seismic stress change, caused by the wave passage, can lead to a seismicity increase even in some portions of the areas far from the mainshock, while in the neighbouring areas we can detect a seismicity quiescence resulting from the stress drop. A large number of studies describe the interference between earthquakes  -  .
After a major seismic event, we observe the development of aftershocks, which generally may occur over a short period of time, due to the dynamic stresses associated with the passage of seismic waves produced by the main shock. Medium-long term aftershocks may develop, which are associated with changes in the static stress    .
In general, dynamic stresses can trigger long distance aftershocks, even more than 10,000 km, where the stress caused by the passage of the mainshock seismic wave, affects the crustal deformation by disrupting the crossed medium’s elastic properties change    . Static stresses have a greater influence in the regions closest to the mainshock    .
The stress changes induced by static stresses, can anticipate or delay a series of aftershocks depending on the sign associated to the stress variation   , while the trigger produced by dynamic stresses, linked to the medium’s elastic variation, is capable of generating seismic events following the main one  .
The branching structure’s long-term development  can provide important clues related to future events’ occurrence. In addition, the analysis of the branching structure suggests that the triggering of other earthquakes can occur after a mainshock, in areas where faults may be in a state of critical stress.
Below we analyze on global scale the upper branching structure of two recent strong earthquakes in order to obtain information on the interaction between the most energetic earthquakes and the energy accumulation and release stage in the upper branching structure. From our studies we infer that the delays due to the co-seismic stress changes’ application depend on the branching structure developmental stage.
We note a delays decrease in the final portion of the branching structure and during the energy release stage, while the different time-scales observed in the higher order branches most likely reflect different tectonic loading rates in the regions.
Figure 5 shows the branching structures of the Ecuador earthquake. The energy accumulation stage that developed after the 6.9 Mw earthquake occurred on 03 April 2016 (point source) consists of eight seismic events featuring a decreasing magnitude over time and delays between them ranging from a few hours to a maximum of 3 days (Table 3). The TP formed on 13/04/2016 with a delay of about 3 days compared to the mainshock. The energy release stage was “Progressive earthquakes”-type and consisted of 4 foreshocks of increasing magnitude and one mainshock  . The foreshocks straight interpolation line (red dashed line) has enabled us to make a preliminary estimate of the mainshock magnitude, while the minimum magnitude, calculated from the
Figure 5. Ecuador earthquake. Branching structures and hazard classes.
Table 3. Coordinates and magnitude of the nodes and the warning signal in the 2016 Ecuador earthquake’s upper branching structure.
midpoint of the first seismic branch, is estimated to be 6.8 Mw (black dashed line). The trigger pattern preceding the mainshock was “Symmetrical triangle”-type  .
The hazard level of the upper branching structure shows a decreasing trend from high to moderate while after the TP it reaches very high values just before the mainshock. The achieved results make us assume a remote interaction between different faults stressed by the seismic waves generated by the most energetic branching structure’s earthquakes.
Figure 6 shows the location of strongest earthquakes epicenters in the energy accumulation (yellow stars) and release (green and red stars) stage and of the three shocks the trigger point consists of (red-brown stars).
It appears that the epicenter of one of the three TP’s shocks is located near the mainshock, while the energy accumulation and release stages have developed far from the epicenter.
A similar developmental pattern is shown by the Nepal earthquake’s branching structures (Figure 7).
The energy accumulation stage that developed after the earthquake on 17 April 2015 (point source) consisted of six earthquakes with decreasing magnitude over time and with delays ranging from a few hours to a maximum of 3 days (Table 4). The TP formed on 23/04/2015 with a delay of about 2 days compared to the mainshock. The energy release stage, consisting of 4 foreshocks of increasing magnitude and one mainshock, was “Progressive earthquakes”-type  . The foreshock straight interpolation
Figure 6. Ecuador earthquake. Earthquakes epicenters of the upper branching structure and TP.
Figure 7. Earthquake in Nepal. Branching structures and hazard classes.
Table 4. Coordinates and magnitude of the nodes and the warning signal in the upper branching structure of Nepal 2015 earthquake.
line (red dashed line) allowed us to calculate with sufficient accuracy mainshock’s, while the minimum magnitude, which was estimated to be 6.45 Mw, was calculated from the midpoint of the first seismic branch (black dashed line). The breaking pattern preceding the mainshock was “Falling wedge”-type  .
The upper branching structure hazardousness shows a decreasing trend from high to moderate while after the TP it reaches very high values just before the mainshock.
In Figure 8 we can note how the epicenter of one of the three TP shocks (red-brown stars) is located close to the mainshock (red star), while the energy accumulation (yellow stars) and release (green stars) stages have developed away from the epicenter.
2.5. Examples of Branching Structures
The examples shown in Figures 9-12 display the main branches in the upper and lower branching structures that have developed before some strong earthquakes as well as the branching structure’s hazard classes developed over time.
The branches in the two branching structures, identified with the long-term seismic sequence hierarchization method, highlight how they tend to converge over time to generate the TP, which in some cases is followed by various order warning signs  . The TP may be followed by “Flash earthquake”-type (Figure 9 and Figure 10)  or “Progressive earthquakes”-type (Figure 11 and Figure 12) energy release stages.
The upper branching structure hazard level tends to decrease from the point source up to TP, then sharply increase up to the main event.
The alternation of the various seismic branches in the upper and lower branching structures show a repetitive development that can be used to obtain information on the seismic sequences’ developmental state. In fact, it is clear that the formation of a certain order branch in the upper branching structure is subsequent to the formation of a
Figure 8. Earthquake in Nepal. Earthquakes’ epicenters in the upper branching structure and TP.
Figure 9. Ecuador earthquake. Branching structures and hazard classes.
branch of the same order in the lower one.
For example, Figure 9, relative to the Ecuador earthquake’s seismic sequence, shows how, after the formation of a fifth order seismic branch in the lower branching structure, a seismic branch of the same order has formed in the upper branching structure. Then a new fifth order seismic branch in the lower branching structure has formed, followed, in time, by the formation of a fourth order rod in the upper structure (seismic branches alternating process).
Usually, this developmental process ends with the convergence of a third order seismic branch in the lower branching structure with a second order seismic branch in the upper branching structure, followed by the TP formation. The repetition of seismic branches of the same order during the branching structures developmental process is suggestive of a seismic delay in the ongoing energy accumulation stage.
Table 5 shows the formation order of seismic branches in upper and lower branching structures and the seismic sequence’s TP of the Ecuador’s earthquake.
The upper branching structure’s hazard level from the point source onward varies from very high to moderately high, while from the TP on, the hazard level returns to very high values. The energy release stage was “Flash earthquake”-type  , with a
Figure 10. Earthquake in Nepal. Branching structures and hazard classes.
delay between the TP and the mainshock approximately of 44 days.
A similar trend is shown in the seismic sequence preceding the Nepal earthquake, on 25 April 2015 (Figure 10) where branching structures started from a fourth order branch and ended with the convergence of the third order branches with a second order one in the upper branching structure. The trigger point was formed on 14 May 2013 while the warning signs preceding the mainshock, generated by a DB-3SE microstructure, appeared on 22 January 2015  . The delay between the TP and the mainshock was about 23 months, whereas that between the trigger point and the mainshock was approximately 3 months. The upper branching structure’s hazard level after the TP decreased from moderate to low up, then returned to high values. The energy release stage was “Flash earthquake”-type.
Figure 11 shows the branching structure of the L’Aquila earthquake, occurred on 6 April 2009. We note a first upper branching structure with a point source that has formed on 15/08/1998, consisting of 5 developmental stages and an early convergence of the lower branching structure. The TP formed on 10 July 2007 was followed by a 4.0 ML magnitude seismic event that generated a second, abnormal fourth order branching structure, with a TP formed on 15/02/2009. The energy release stage was “Progressive
Figure 11. L’Aquila’s earthquake. Branching structures and hazard classes.
earthquakes”-type, with the last trigger point generated on 02/04/2009. The delays between the first and the second TP compared to the mainshock were approximately 21 months (duration of the whole energy release stage) and 46 days, respectively, while the delay between the last TP and the mainshock was around 4 days.
The upper branching structure hazard level shows a decreasing trend from high to low up to the first TP, while during the energy release stage, it increased from moderate to potentially moderate-high if we consider the trendline joining the earthquakes the energy release stage.
Figure 12 shows the branching structures of the earthquake occurred in Emilia on 20 May 2012. The point source from which a fourth order upper branching structure has developed, has formed on 15/10/1996, while the TP was generated on 21/09/1997.
The energy release stage that followed was “Progressive earthquakes”-type with the last strongest earthquake, featuring a magnitude of 4.5 Md, recorded on 18/06/2000. This event is a second point source that generated a second third order upper branching structure and a TP that has formed on 22/02/2008.
The energy release stage consisted of two first order foreshocks of 4.5 Mw and 4.9 Mw magnitude respectively and a second order foreshock of 4.0 Md magnitude, recorded on 19/05/2012. The second order foreshock was preceded by a warning signal that was generated on 01/03/2012. The delay between the second TP and the mainshock
Figure 12. Earthquake in Emilia. Branching structures and hazard classes.
Table 5. Development of branching structures.
was approximately 4 years and three months (duration of the whole energy release stage), the delay between the last TP on 01/03/2012 and the mainshock was about 80 days, while the delay between the second order foreshock and the mainshock was one day.
The upper branching structure hazard level remained moderate/moderate low until the mainshock and then became moderately high.
By using the nodes of upper and lower branched structures’ seismic branches, we have proposed new approaches to define a classification of the hazard levels reached by the upper branched structure during its development. This classification plays a crucial role in determining the hazardousness of the energy release phase that is triggered following the TP formation.
In fact, the upper and lower branched structures, in the last part of their development, show a convergence of branches of the same order with a lower order one, which is an important warning signal that usually precedes a strong earthquake and determines a warning sign of the sequence’s hazardousness and its development.
Analyzing the upper branched structure’s temporal and spatial changes both on global and local scale, it appears that they are adequate enough to assume a long-dis- tance interaction between the various earthquakes.
These interactions should, however, be linked both to the static stress changes and the passage of seismic waves from the most energetic earthquakes. These can generate transient stresses in the earth’s crust, thus triggering earthquakes far from the original earthquake source.
We also believe that recurrence rates that can be observed in the branching structure’s first order seismic branches most likely reflect the different tectonic loading rates of the seismic regions concerned.
Conversely, in the last stages, the decrease in the recurrence rates is probably due to a very high tectonic load, where small shocks can trigger the energy release stage.
 Pollitz, F.F. and Johnston, M.J.S. (2006) Direct Test of Static Stress versus Dynamic Stress Triggering of Aftershocks. Geophys?cal Research Letters, 33, L15318. http://dx.doi.org/10.1029/2006gl026764 https://profile.usgs.gov/myscience/upload_folder/ci2010Nov2222030742871152.pdf
 Riga, G. and Balocchi, P. (2016) Short-Term Earthquake Forecast with the Seismic Sequence Hierarchization Method. Open Journal of Earthquake Research, 5, 79-96.
 Harris, R.A. (1998) Introduction to Special Section: Stress Triggers, Stress Shadows, and Implications for Seismic Hazard. Journal of Geophysical Research: Solid Earth, 103, 24347-24358. http://dx.doi.org/10.1029/98jb01576
 King, G.C.P. and Cocco, M. (2000) Fault Interaction by Elastic Stress Change: New Clues from Earthquake Sequences. Advances in Geophysics, 44, 1-38. http://dx.doi.org/10.1016/S0065-2687(00)80006-0
 Stein, R.S., Barka, A.A. and Dieterich, J.H. (1997) Progressive Failure on the North Anatolian Fault Since 1939 by Earthquake Stress Triggering. Geophysical Journal International, 128, 594-604. http://dx.doi.org/10.1111/j.1365-246X.1997.tb05321.x
 Brodsky, E.E. and van der Elst, N.J. (2014) The Uses of Dynamic Earthquake Triggering. Annual Review of Earth and Planetary Sciences, 42, 317-339. http://dx.doi.org/10.1146/annurev-earth-060313-054648
 Freed, A.M. (2005) Earthquake Triggering by Static, Dynamic, and Postseismic Stress Transfer. Annual Review of Earth and Planetary Sciences, 33, 335-367.
 Delorey, A.A., Chao, K., Obara, K. and Johnson, P.A. (2015) Cascading Elastic Perturbation in Japan Due to the 2012 Mw 8.6 Indian Ocean Earthquake. Science Advances, 1, e1500468. http://dx.doi.org/10.1126/sciadv.1500468
 Lin, J. and Stein, R.S. (2004) Stress Triggering in Thrust and Subduction Earthquakes and Stress Interaction between the Southern San Andreas and Nearby Thrust and Strike-Slip Faults. Journal of Geophysical Research: Solid Earth, 109, B02303. http://dx.doi.org/10.1029/2003jb002607
 Belardinelli, M.E., Bizzarri, A. and Cocco, M. (2002) Earthquake Triggering by Static and Dynamic Stress Changes. Journal of Geophysical Research: Solid Earth, 108, 2135.
 Riga, G. and Balocchi, P. (2016) Seismic Sequence Structure and Earthquakes Triggering Patterns. Open Journal of Earthquake Research, 5, 20-34. http://dx.doi.org/10.4236/ojer.2016.51003 http://www.scirp.org/Journal/PaperInformation.aspx?PaperID=63538
 Riga, G. and Balocchi, P. (2016) How to Predict Earthquakes with Microsequences and Reversed Phase Repetitive Patterns. Open Journal of Earthquake Research, 5, 153-164.