Data in Brief 17 (2018) 588-603,
Contents lists available at ScienceDirect
Data in Brief
ELSEVIER journal homepage: www.elsevier.com/locate/di
Data Article
Evaluation of O < M < 8 earthquake data sets in ®
African - Asian region during 1966-2015 | Soe
Theophilus Aanuoluwa Adagunodo*™, Sebastian Lining”,
Adekunle Michael Adeleke ‘, Julius Oluwasegun Omidiora“,
Ahzegbobor Philips Aizebeokhai*, Kehinde David Oyeyemi*,
Olaide Sakiru Hammed *
* Department of Physics, Covenant University, Ota, Ogun State, Nigeria
® Institute for Hydrography, Geoecology and Climate Sciences, HauptstraSe 47, 6315 Ageri, Switzerland
* Department of Physics, University of Morin, loin, Kwara State, Nigeria
* Department of Languages and General Studies, Covenant University, Ota, Ogun State, Nigeria
Department of Physics, Federal University of Oye Ekit, Oye Ekiti, Ekii State, Nigeria
ARTICLE INFO ABSTRACT
“ticle history This article evaluates the occurrence of 0 < M = 8 earthquake data
Received 14 August 2017 sets for the period of 50 years (that is, January 1, 1966 to December
Received in revised form 31, 2015) in African and Western Asia region. It is bounded by
30 November 2017
Accepted 16 January 2018
Available online 31 January 2018
latitude 40° S to 40° N and longitude 30° W to 60° E with the focal
depth of 0-700km. Seventy seven thousand, six hundred and
ninety-six data points were presented for the analysis. The data
Keywords: used were extracted from earthquake catalog of Advanced National
African plates Seismic system via htp://quake-geo.berkeley.edujenss), an official
aban plates ‘website of the Northern California Earthquake Data Centre, USA.
Each datum comprised the earthquake occurrence date, time of the
Gutenberg-Richer's model
finch newrlnetwork earthquake occurrence, epicenter’s coordinates, focal depth and
Farthquake magnitudes ‘magnitude. The Gutenberg-Richter’s relationship being the longest
Focal depth observed empirical relationship in seismology, analysis of variance
Seismic events and time series were used to analyze the seismicity of the study
Seismograps area Annual distributions of earthquake occurrence based on
magnitude variations with the limit 0 < M < 8 were presented,
‘The two constants a and b in the Gutenberg-Richter's equation,
‘magnitude of completeness (MC) adjusted R-Square and F-value for
the period of 1966-1975, 1976-1985, 1986-1995, 1996-2005, 2006-
2015, and the entire period of investigation ranging from 1966 to
Tectonic stress
* Corresponding author,
E-mail addresses:
[email protected], taadagunodo@gmailcom,
theophilus adagunodo@covenantuniversity edu ng (TA. Adagunodo)
https: doi org/10.1016]j.dib.2018.01049
2352-3409}0 2018 The Authors, Published by Elsevier Inc. This is an open access article under the CC BY license
(http: jJereativecommons orgilicenses/by/4.0).TA. Adagumodo et al / Data in Brief 17 (2018) 588-603 589
2015 were determined so as to investigate the variations of these
parameters on earthquake occurrence over time. The histograms of
earthquake occurrence against magnitude of earthquakes for the
selected years (1966-1975, 1976-1985, 1986-1995, 1996-2005,
2006-2015, and 1966-2015), and the decadal frequency distribu-
tions of earthquake occurrence were also plotted. The focal depth
occurrence for each magnitude bins (0-0.9, 1-19, 2-29, 3-3.9, 4
49, 5-59, 6-6.9, 7-79, 8-8.) were grouped into shallow, inter
mediate, and deep depths ranging from 0 to 70, 71 to 300, and 301
to 700 km as being used in seismology. The neural network analysis
was also applied to the magnitude of the earthquake. The network
uses a time series magnitude data as input with the output being
the magnitude of the following day. tthe nature of the earthquakes
time series is stochastic, modeling and prediction is possible. The
earthquake data sets presented in this article can further be adopted
in the study of seismicity pattern, b-value using series of models,
earthquake prediction and variations of earthquake parameters on
African andjor Arabian plates. When this approach is integrated
with other techniques), it can provide insights to stability of African
lithospehric plates especially the coastal region of Africa.
© 2018 The Authors. Published by Elsevier Inc. This is an open access
article under the CC BY license
(http: /ereativecommonsorgilicenses/by/4.0/).
Specifications Table
Subject area
More specific sub-
ject area
‘Type of data
How data was
acquired
Data format
Experimental
factors
Experimental
features
Data source
location
Data accessibility
Computational Geophysics
Earthquake
Table and figure
The seismic events were recorded by the seismographs of the Northern California
Farthquake Data Centre, USA.
Raw and processed
The data were extracted from the earthquake catalog of Advanced National
Seismic system.
Computational analysis of earthquake parameters for the period of 50 years
(1966-2015) using Microsoft Excel, SPSS and MATLAB R2013a software.
The data were obtained for 0 < M < 8 earthquake latitude 40° S to 40° N and
longitude 30° W to 60° E, focal depth distribution from 0 to 700 km for the
period of January 1, 1966 to December 31, 2015. There were 77,696 data points in
all,
The data sets are with this article. It is also available on hittp://quake.geo.ber
keley.edu/cnss).
Value of the data
Can be used to study the seismicity pattern in African and/or Western Asia region.
Can be used for b-value estimation using integrated models in African ~ Western Asia seismology.
Can be used to study the effect of earthquake occurrence on African and/or Arabian lithospheric
plates.
Can be used to estimate the time scale dependence of earthquake parameters in subregions of
Africa (Northern, Central, Western, Southern and Eastern Africa) (Fig. 1) and Middle East.
Can be used to forecast the earthquake occurrence in African and/or Western Asia region.590 TA, Adagunodo etal. / Data in Brief 17 (2018) 588-603
Can be integrated with other computational approach for earthquake interpretation.
Can be used to further explain the stability of African lithospheric plates.
For educational purposes on seismically active zones in African ~ Asian region.
Can be correlated with other earthquake data for seismic activity studies in coastal region of Africa
and Middle East.
Can be employed in the study of seismic activities around the equator when integrated with other
techniques such as aeromagnetic data and geographic information system approach.
© It can provide insights to further exploration of aseismic zones being affected by tremors in Africa
especially Nigeria.
eee
1. Data
The data in this article contains the record of earthquake occurrence in African - Western Asia
region. The seismic events were recorded by the seismographs of the Northern California Earthquake
Data Centre, USA. The data were obtained for the 0 < M < 8 magnitude between latitude 40° S to
40° N and longitude 30° W to 60° E (Fig. 2), focal depth distribution from 0 to 700 km for the period of
January 1, 1966 to December 31, 2015. There were 77, 696 data points in all. Each datum comprised
UN Subregions of Africa
Northern Africa
[1] Western Atrica
Gi Central Africa
[7] Eastern Africa
Il Southern Africa
Fig. 1. Map of Africa showing her five subregions.TA. Adagumodo et al / Data in Brief 17 (2018) 588-603 591
i Turkey
bertonear Syria on ASIA
PS 0 mie east
Ue cog 38, / We
o ‘Qatar— 7
Atlantic Ocean
30"
@ worldatlas |
= . OE
30 w isw r ISE We 45° E OGraphichtaps.com
Fig. 2. The map of Africa and Asia showing the coverage of the dataset.
Indian Ocean
the earthquake occurrence date, time of the earthquake occurrence, epicenter’s coordinates, focal
depth and magnitude.
‘An earthquake is caused by a sudden slip along a fault zone. It has been recognized as one of the
most destructive of all natural hazards which can severely destroy the entire vicinity in seconds
without an explicit warning [1-3]. Evaluation of earthquake parameters such as magnitude, focal
depth and frequency is the fundamental in the study of earthquake pattern and its prediction. These
parameters are essential in seismology and serve as reference point to the applied theoreticians. In
the study of earthquake, one of the most determined parameters is b-value which varied from 0.2 to
2.0, and generally found around 1. This is the measure of stress on the lithospheric plates, because
lower b-values indicate that the stress is optimum in the investigated region. Generally, very low b-
values are found in case of immediate aftershocks and higher values are found in case of swarm.
There are two mostly used methods in estimation of b-values: least square and maximum likelihood
methods [4]. In this data article, least square method which is based on Gutenberg-Richter’s (GR)
relationship has been adopted to evaluate the data sets of 0 < M < 8 magnitude in African ~ Arabian
region to determine the decadal variations of seismicity levels (a-values) and tectonic character (b-
values) for the period of 50 years (1966-2015). Globally, GR equation has been applied to earthquake
data for the estimation of b-values and related parameters, but few reports from African continent are
available in the literature [5] ~ a gap that is essential to be bridged in the study of seismic activities in
Africa and Western Asia region.592 TA, Adagunodo etal. / Data in Brief 17 (2018) 588-603
The GR law has remained one of the oldest empirical relationships that are still relevant in seis
mology till date. The relationship is based on power scaling relationship, which relates the frequency
and the magnitude of earthquake together in order to predict the degree of stress on the lithospheric
plate in a region. The GR equation [6] is presented in Eq. (1).
LogioN =a-bM a
where N is the cumulative number of earthquakes of magnitude > M,a characterizes the seismicity level
of a region, which represents the M > 0 earthquake. b defines the tectonic character, which is a function
of the accumulated stress of a region. In addition, a and b are constants that vary in space and time.
The descriptive analysis has also been found useful in the evaluation of earthquake occurrence in a
region. This ranged from description of earthquake occurrence by plotting the graphs of frequencies
of earthquakes against their coordinates [5], number of earthquakes against its magnitudes [7,8],
cumulative number of earthquake against its magnitude [7|, and measure of central tendencies [9]. It
has been reported that study of previous and present activities of earthquake pattern is vital in
prevention of lives and properties from earthquake destructions [9].
Furthermore, in earthquake predictions, several phenomena have been considered by researchers.
The considered parameters are electromagnetic fields, seismicity pattern, unusual cloud and weather
parameters, unusual emanation of hydrogen and radon gases from the subsurface (e.g. groundwater
or soil), animal behaviours [10], and unbalancing level in surface and groundwater [3]. The most
unsolved issues in seismology, that is, the time, location and magnitude of the impending earthquake
are the major aim of earthquake prediction which can further be improved on via the approach
presented at the latter part of this data article.
The neural network developed in this article uses only time series magnitude data as input with the
output being the magnitude of the following day. Time series is defined as a sequence of values
documented in chronological order over time. Occurrence of previous events may be extremely valuable
in prediction of its behaviours in the future. As reported by [11] that, ‘if given a set of past values, it is not
possible to predict future values with reliability, the time series is said to be chaotic’. However, if the
nature of the earthquakes time series is stochastic, modeling and prediction is possible. The available
data in this article can be adapted by the seismologists in understanding, modeling and prediction of
earthquake occurrence in African ~ Arabian region. Furthermore, this analysis can be integrated with
other computational approach for better earthquake interpretation. Similar computational analyses to
solve other challenges in Man’s environment have been presented in [12-14].
2. Study area, Tectonic Settings and its Geology
African-Arabian or Western Asia region constitutes all the countries presented in Fig. 2. The study
area is bounded by latitude 40° $ to 40° N and longitude 30° W to 60° E. The African plate has recently
been reported as the third largest plate [5]. It is bounded by a total area of about 60 million square
kilometer, with about half of it being covered by land. African plate is composed of old Cratonic units
and growth of younger Crust, which represent a period > 2.5 billion years of oceanic and continental
crust growth [15]. The African plate is a significant tectonic plate bestriding the equator and the prime
meridian. It encapsulates larger percent of the African continent, as well as oceanic Crust which
reclines between series of oceanic and continental ridges. The Arabian plate is a minor tectonic plate
that falls on the eastern and northern hemispheres. It is one of the three continental plates (the
Arabian, Indian, and African plates) that have been moving northward in the recent geological record,
and colliding with the Eurasian plate. The African-Arabian region is composed of five tectonic plates:
Madagascar, Arabia, Seychelles, Nubia and Somali as presented in Fig. 3. The historical record showed
that African tectonic setting was constituted by the breakup of Gondwana in 200 Ma (Mega-annum).
This resulted to the interaction of Nubia with Eurasia along the former northern margin [16]. During
160-117 Ma, Madagascar separated from southeastern Africa and rifted to its present location. During
the Oligocene (that is, between Eocene and Miocene), the Neotethys Sea (previously located between
Nubia and Eurasia) closed through subduction as the two plate collided [17]. The Arabian plate got
separated from Africa about 25 Ma ago. This separation led to the closure of the Neotethys Sea, withTA. Adagumodo et al / Data in Brief 17 (2018) 588-603 593
oo. 2 2 0 0 0 oe Oo wm Ow
Fig. 3. Aftican-Arabian tectonic plates Adapted from [17|
the succeeding rifting which lead to the formation of Red Sea. From 10 to 60 Ma, the Somali plate
began to rift of from African plate. Logatchev et al. [18] predicted another Sea and a new continent
between Somali and Africa in the next 1 and 10 Ma respectively. Madagascar and Seychelles (Plateau)
plates are microplates within the Somali plate. About 84 Ma, a spreading ridge formed a new location
in the Indian Ocean, from which the Mascarene Basin was formed. Further rifting between Seychelles
and India at the Tertiary or Cretaceous boundary resulted to the hot spot magmatism, which further
sedimented to produce the carbonate shelves on the microplate [19]
Generally, Arabian and African continent are made up of a PreCambrian basement of crystalline
meta-sedimentary, igneous and meta-igneous rocks (Fig. 4). This crystalline basement is overlain by
series of geological settings ranging from volcanic and sedimentary sequences to unconsolidated
Cenozoic sediments (20). African continent is made up of primary units known as Cratons, which are
the aforementioned sediments or weathered rocks overlying the crystalline basement.These Cratons
are predominantly granitic series, gneisses, and low-grade greenstone belts [20].
3. Experimental design, materials and methods
‘The magnitude of an earthquake is determined based on the information received by the seis-
‘mograph. The Richter magnitude involves measuring the amplitude of the largest recorded wave at a
specific distance from the seismic source. The magnitude of earthquake and its implications are
presented in Table 1.
The annual distributions of seismic activities based on the magnitude of earthquakes are pre-
sented in Table 2. The 0 < M < 0.9 earthquakes showed total events of 154,1 < M < 1.9 earth-
quakes showed total events of 2347,2 < M < 2.9 earthquakes showed total events of 17640, 3 < M
< 3.9 earthquakes showed total events of 33010, 4 < M < 4.9 earthquakes showed total events of
20922, 5 < M < 59 earthquakes showed total events of 3388,6 < M < 69 earthquakes showed
total events of 216,7 < M < 7.9 earthquakes showed total events of 18, and 8 < M < 8.9 earth-
quakes being the least recorded event occurred once in 1969. Table 2 revealed that African-Arabian594 TA. Adagunodo et al. / Data in Brief 17 (2018) 588-603
Bie is) PLATE TOPICS
Kaapreat
(Craton
gseag
ue
{
i
oH
y
|
Table 1
‘The Richter magnitude and its implications.
Magnitude Implications
ot ‘Cannot be felt, but itis detectable by seismograph
2 ‘Smallest quake to be felt. Hangling objects may swing
3 People near the epicenter feel the quake. Comparable to vibra-
tions of a passing truck
4 Causes damage around the epicenter. It isthe same as a small
fission bomb
5 “The weak buildings around the epicenter are damaged
6 Causes greater damage around the epicenter
7 Causes serious damage. Capable to create energy that will heat up
4 country, It can be felt globally
8 Causes major destruction and death
9 Rare, but can cause unbelievable damage or total destructionTA. Adagumodo et al / Data in Brief 17 (2018) 588-603 595,
Table 2
Yearly distribution of seismic activities according to the magnitude of earthquakes.
Year! 009 «1-19 «2-290 3:39 4490 5-59 6697-79889
Mag.
1966 0 0 0 9 230 4 4 0 0
1967 0 0 0 5 179 80, 5 0 0
1968 ° 1 oO B 209 5 4 1 0
1969 0 o 0 9 241 33 9 0 1
1970 ° o oO 3 268 64 3 1 °
1971 ° o 0 6 223 59 3 1 0
1972 ° ° 0 15 182 37 5 0 0
1973 0 0 0 38 218 62 3 0 0
1974 0 0 1 41 183 46, 3 0 0
1975 ° o 2 106 266 92 4 1 °
1976 ° o 7 18 255 59 3 1 °
197 ° o o 35 262 20, 5 0 0
1978 ° o 1 78 216 44 3 1 °
1979 0 0 1 7 27 R 4 1 0
1980 0 0 1 155 327 54 6 1 0
1981 0 0 u 302 480 55 7 0 0
1982 0 0 9 254 33 54 2 0 0
1983 0 0 10 25 500 96 3 0 0
1984 0 o 7 259 333 54 0 0 0
1985 ° ° B 222 343 45 5 0 0
1986 ° ° 0 270 359) 46 2 0 0
1987 ° 1 16 222 337 40 0 0 0
1988 0 o 2 ‘400 397 4 1 oO 0
1989 0 1 184 801 385, n 5 0 0
1990 2 8 505, 843, 537 70 3 3 0
1991 1 “4 459 838 344 44 0 oO 0
1992 2 6 360 845, 363 55 3 1 0
1993 8 18 2097 114 409) 82 2 0 0
1994 0 19 1998 1st 324 60. 6 o 0
1995 ° 25 957 ms 355, 56 4 1 0
1996 3 2 527 716 509) 61 7 0 0
1997 ° 9 286 366 567 57 3 1 0
1998 ° 234 102 652 418 a7 10 0 0
1999 ° 132 ou 959 416 44 3 0 °
2000 1 213 676. 73 381 55 2 1 °
2001 ° 166 a2 1276 448, 50 4 0 °
2002 0 248 407 280485, 5 8 o 0
2003 80 n9 1197 1794 438 68 9 oO 0
2004 at 433 121 2165 636 7 4 0 0
2005 ° 15 849 2861 763 80 5 0 0
2006, 0 3 737 2m 706 96, 6 1 0
2007 1 10 682 3306737 96 6 0 °
2008 ° o 785 3314 670 101 " oO 0
2009 0 0 0 10 365 98 3 0 0
2010 0 ° 8 531 105 3 0 0
2011 ° o 0 10 696 85 6 1 0
2012 2 o oO 2 583 20, a o 0
2013 3 o 1 25 580. 69. 6 o 0
2014 ° o 0 28 803 78 5 oO °
2015 0 o 4 2 619 67 3 1 0
Total 154 2347 17640-33010» 208223388216 8 1
magnitude of earthquakes fluctuates between 3 < M < 59, Seismic events were recorded yearly
within this range, which varied from as low as 2 events in 2012 for 3 < M < 3.9 to very high events
of 3314 in 2008 for the same range of magnitude.
‘The data sets were further explored by constructing the histograms of frequency of earthquake
occurrence against the magnitude for the period of 50 years (1966-2015) (Fig. 5a), and 10-year596 TA, Adagunodo etal. / Data in Brief 17 (2018) 588-603
é
b=
= |
i j
t. i
f
cm _
Fig. 5. a. Earthquake occurrence - magnitude distribution from 1966 to 2015 b. Earthquake occurrence - magnitude di
tribution from 1966 to 1975 e. Earthquake occurrence - magnitude distribution from 1976 to 1985 d. Earthquake occurrence ~
‘magnitude distribution from 1986 to 1995 e. Earthquake occurrence - magnitude distribution from 1996 to 2005 f. Earthquake
occurrence ~ magnitude distribution from 2006 to 2015,TA. Adagumodo et al / Data in Brief 17 (2018) 588-603 597
i on ov
i~
aH |
Fig. 6. Decadal frequency distribution of earthquake occurrence in all the data sets.
Table 3
‘Magnitude ~ focal depth relationship of the earthquake occurrence,
Magnitude ‘Shallow depth Intermediate depth Deep depth
range (0-70 km) (71-300 km) (301-700 km)
0-09) 154 0 °
119 2343 4 °
229 ara 165 1
339 32,510 478 2
449 20.232 699 38
559 3281 60 °
669 ait 5 0
7-19 18 0 °
889 1 0 °
Total 76223 412 6
interval: 1966-1975 (Fig. 5b), 1976-1985 (Fig. 5c), 1986-1995 (Fig. 5d), 1996-2005 (Fig. 5e), and
2006-2015 (Fig. 5f). The histogram of Fig. Sa (1966-2015) revealed that events of 3 =< M <
3.9 earthquakes were most frequent, with the total events of 33010. In Fig. 5b and c (1966-1975 and
1976-1985), the most frequent magnitude is 4 < M < 4.9, with the total events of 2199 and 3302
respectively. However, Fig. 5d-f (1986-1995, 1996-2005, 2006-2015) showed that the most frequent
magnitude is 3 < M < 39, with the total events of 8062, 13392, and 9496 respectively. The decadal
frequency distribution of earthquake occurrence plotted in Fig. 6 showed that 1996-2005 is the most
active decade of seismic events, with 30283 being the total earthquake occurrence in this decade.
The focal depths of the earthquake data sets were also analyzed. The focus of an earthquake is the
actual point where the rocks break in the subsurface. The depths of the focus are categorized as
shallow (0-70 km), intermediate (71-300 km), and deep depth (301-700 km below the earth surface)
respectively [8]. Most earthquakes occurred at shallow depth, with total events of 76223, while deep
depth had a total record of 61 times, with the magnitude ranging from 2 to 4.9. The focal depth
distributions of earthquakes in the African-Arabian region during 1966-2015 are presented in Table 3.
The frequency-magnitude distributions of earthquake occurrences for the period of 50 years
(1966-2015) and 10-year interval covering the entire investigated period were produced through the
GR relation (Eq. (1)). This was achieved by plotting the graph of cumulative number of earthquakes
against their respective magnitudes. The graphs were fitted, with a linear fitting. The equation from
the linear fitting represents the GR relation, where the slope of the graph stands for the b-value.8
TA. Adagunodo et al. / Data in Brief 17 (2018) 588-603
ae d
8 g° a92042
Z° 3 bao. = 0.03
€ g Me=1.80
2 8 Ag.R square = 087
: 5
z zg
& 3
Ze ze
5 5
b° e
z z-
4 2
& é
a. § Adj R square =087
% 5
z a?
z z
3 q
5 5
g g
inquaes Mogae of eatnquates
c f
ge f- assaat oie
3 3 b=062 = 003
é a Mc=1.90
3 i ‘Ad, R square = 083,
i e
3 3
3 5
Zo ze
g 8
Magatuce of earthquakes Magnus of eatmquakes
Fig. 7. a. LogN - magnitude distribution From 1966 to 2015 b. LogN - magnitude distribution from 1966 to 1975 e. Log N -
magnitude distribution from 1976 to 1985 d. Log N - magnitude distribution from 1986 to 1995 e. LogN - magnitude dis-
tribution from 1996 to 2005 f. Log N - magnitude distribution from 2006 to 2015,
However, the a-value (seismicity level) was estimated by substituting the known parameters into the
GR equation. The magnitude of completeness (Mc) of the earthquake catalogue was also determined
for each period. The Mc is referred to as the threshold magnitude which is the magnitude above
which all earthquakes were recorded.
The GR plots for the period of 1966-2015, 1966-1975, 1976-1985, 1986-1995, 1996-2005, 2006-
2015 were presented as Fig. 7a-f, with the a-value, b-value, Mc value and Adjusted R-square (Adj.TA. Adagumodo et al / Data in Brief 17 (2018) 588-603 599
R-square) value of each analysis being displayed on each graph. The Adj. R-square is a corrected
goodness of fit (model accuracy). It is calculated by dividing the residual mean square of error by the
total mean square error. The complete data sets during the period of 1966-2015 revealed that the b-
value of the study area is 0.61, a-value is 6.01, Mc is 1.85, and Adj, R-square is 0.84 respectively. The
Table 4
ANOVA results and summary of the parameters from frequency-magnitude curves.
Year Source DF Sum of squares Mean of square F-value Prob > F a-value p-value Mc Adj. R-square
1966-2015 Model 1 16530 16530 4073 0 601 061185 O84
Error 79 20.32 038
Total 80 19562
1966-1975 Model 1 8791 8701 26762 0 440 045 200 0.77
Error 79 25.95 033,
Total 80 113.86
1976-1985 Medel 1 11800 118.00 a0 0 478 052-195 080
Error 79 28.87 037
Total 80 14686
1986-1995 Model 1 17146 146 52775 0 539 062 1.80 087
Error 79 25.67 032
Total 80 19713
1996-2005 Model 1 19504 19504 53507 0 560 065 180 087
Error 79 2727 035
Total 80 21231
2006-2015 Model 1 16818 16818, 40509 0 543 062 1.90 083
Error 79 32.80 042
Total 80 20098
§
da
‘J
> AU LLL ay JUL
* # + + z ¢ % ¥ *
Fig. 8. Neural network training of error autocorrelation, Epoch 291600 TA, Adagunodo etal. / Data in Brief 17 (2018) 588-603
analysis across the five (5) decades showed that the b-value fluctuates from 0.45 (during the first
decade) to 0.65 (during the fourth decade), a-value from 4.40 (during the first decade) to 5.60 (during
the fourth decade). The Mc varied from 1.80 to 2,00.
Statistical analysis involving the Analysis of Variance (ANOVA) was carried out on the magnitude
of the earthquakes. This analysis is used to determine the difference in the mean of sets of data or
groups. The ANOVA uses F-tests to statistically analyze or test the equality of means. The test was
named after Sir Ronald Fisher. The F-test is the ratio of two variances (measure of dispersion), being
the square of standard deviation. The ANOVA results of the data sets covering the period of 1966-
2015, and five decades (1966-1975, 1976-1985, 1986-1995, 1996-2005, and 2006-2015) are pre-
sented in Table 4, The ANOVA results showed that F-values recorded for the periods of 1966-2015
(being the whole coverage year for the analysis) is 430.73, 1966-1975 is 267.62, 1976-1985 is 322.91,
1986-1995 is 527.75, 1996-2005 is 535.97, and 2006-2015 is 405.09 respectively.
sa ‘Conan btween put 1 ancrr1= Target? - Output
I lll,
Contain
us
Fig. 9. Neural network training of input-error cross-correlation, Epoch 291
a Ever stogram wih 20 88
PEPSt ET EE a
Eom» Targets - Opts
Fig. 10. Neural network training of error histogram with 20 bins, Epoch 291Output ~= 0.26"Target + 2.4
Output ~= 0.26*Target + 2.4
TA. Adagunodo et al. / Data in Brief 17 (2018) 588-603,
Best TrtnngPertermancels 886 atopoch 200
21 poche
Fig. 11. Neural network training of performance 064566 at epoch 200 and iteration 20.
Training: R=0.51201
Test: R=0.5042
=
a
+
3
2
gi
Fr
&
°
4
3
ce}
4 6
Target
All: R=O.51082
Target
Fig, 12, Regression analysis of the dataset with R-square of
05 in all cases.602 TA, Adagunodo etal. / Data in Brief 17 (2018) 588-603
Response of Oupit Een fr Tne Svs
Fig. 13. Time series response for the dataset, Epoch 291,
The earthquake forecasting approach from historic seismic data is prevalent nowadays [11], which
was further used to evaluate the earthquake data sets. This technique employing a dynamic neural
network was further used to explore the earthquake data sets (the neural network function and the
source code were attached as Supplementary files). This approach is good at time series prediction. With
the total magnitude of 77, 696, graphical user interfaces and command-line functions were used to
produce the codes and figures for earthquake prediction in this article. The dataset was selected and the
problem to solve was defined in the MATLAB. The network was trained in order to fit a time series data
set. The beauty of using neural network time series tool (ntstool) is that it is capable to solve three
problems differently. The first task which was adopted in this analysis is to predict future values from
previous values y(¢) and past values from a second time series x(t) using a nonlinear autoregressive with
exogenous input (NARX). The second task is to only predict future values from the past values of such
time series using a nonlinear autoregressive (NAR) input. The third task is to predict future values from
the previous values without having knowledge of previous values using input-output model instead of
NARX. The magnitudes of the earthquake were imported, validated and tested such that 70% of the data
were used for the training while 15% each will be used for validation and testing of the dataset. The
graphical interfaces of the treated data were presented from Figs. 8-13.
The qualitative and quantitative approach employed in this article can be beneficial in the study of
stability of the African lithospheric plates. The data sets can be reexamined to estimate the time scale
dependence of earthquake parameters in subregions of Africa andjor Middle East. The approach pre-
sented in this article can provide insights to researchers on further explorations of aseismic zones being
affected by tremors in Africa such as Nigeria [21-23]. However, despite the challenges being faced in
earthquake predictions, it has been noted that the beauty of neural network is to predict the next major
seismic event [1. Analysis of the dataset presented in this article can be used to forecast the earthquake
occurrence in African - Western Asia region. If accurate forecast is achieved in this region, it would be
beneficial for the masses since danger of loss of lives and properties would be reduced.
Acknowledgements
Dataset in this article were accessed through the Northern California Earthquake Data Center
(NCEDC), doi: 10.7932/NCEDC. We also acknowledge the Centre for Research, Innovation and Dis-
covery, Covenant University, Ota, Nigeria for partial sponsorship.TA. Adagumodo et al / Data in Brief 17 (2018) 588-603 603,
‘Transparency document. Supplementary material
Supplementary data associated with this article can be found in the online version at htt)
org/10.1016/j.dib.2018.01.049,
Appendix A. Supplementary material
Supplementary data associated with this article can be found in the online version at http://dx.doi.
org/10.1016/j.dib.2018.01.049.
References
[1] TA. Adagunodo, LA. Sunmonu, Earthquake: a terrifying of all natural phenomena, . Adv. Biol. Basic Res. 1 (1) (2015) 4-11.
[2] EFlorido, LL Aznarte, A. Morales-Estebam, F. Martinez-Alvare2, Earthquake magnitude prediction based on artificial neural
networks: a survey, Croat. Oper. Res. Rev. 7 (2016) 159-169,
[3] M. Moustra, M. Avraamides, C. Christodoulou, Artificial neural networks for earthquake prediction using time series
magnitude data or seismic electric signals, Expert Syst. Appl. 38 (2011) 15032-15039,
14] W. Marzocchi, . Sandri, A review and new insights on the estimation of the b-value and its uncertainty, Ann. Geophys. 46,
(6) (2003) 1271-1282.
15] OS. Hammed, MO, Awoyemi, WIN. Igboama, C.O. Badmus, UC. Essien, Pattern of seismicity associated with the African,
lithospheric plate, Phys. Sci. Int. . 10 (3) (2016) 1-11,
{6] 8. Gutenberg, C. Richter, Magnitude and energy of earthquakes, Ann. Geofis. 9 (1956) 115,
[7] 0S. Hammed, M.0. Awoyemi, C.0. Badmus, 0.0. Sani, Interdependence an variations of earthquake parameters on Aftican
lithosperic plate using Guterberg and Richter relations, | Geogr. Environ. Earth Sci, Int. 4 (4) (2016) 1-12.
[8] OS. Hammed, 0. Popoola, AA. Adetoyinbo, M.O. Awoyemi, CO. Badmus, 0.8. Ohwo, Focst depth, magnitude, and fre-
‘quency distribution of earthquakes along oceanic trenches, Eartha, Sc. 26 (2) (2013) 75-82.
19] AA. Alabi, 01D. Akinyemi, A. Adewale, Seismicity pattern in southern Africa from 1986 to 2009, Earth Sci. Res. 2 (2) (2013)
01-10,
[10] RA. Grant, T. Haliday, Predicting the unpredictable; evidence of pre-seismic anticipatory behaviour in the common toad,
J Zool, 281 (4) (2010) 263-271.
[11] A. Morates-Esteban, F. Martinez-Alvarez, A. Troncoso, JL Justo, C Rubio-Escudero, Pattern recognition to forecast seismic
time series, Expert. Syst. Appl. 37 (2010) 8333-8342.
[12] TA. Adagunodo, LA. Sunmonu, M.A. Adabanija, EA. Suleiman, O.A. Odetunmibi, Geoexploration of Radioelement’s
datasets in a flood plain of Crystalline bedrock, Data Brief 15 (2017) 809-820. http://dx-do\.org/ 101016)}.