Preliminary Seismic Hazard Assessment in West Africa Based on Incomplete Seismic Catalogs ()
1. Introduction
West Africa spans the area between latitude 0˚N to15˚N and longitude 20˚W to 10˚E. It corresponds to a vast complex that includes the West African Craton, the Pan-African mobile belt and the Mauritanides fold belt (Figure 1).
The edification of West Africa has developed grabens, horsts, and tectonic lineaments which are superimposed on older structures. The tectonics of the Western Africa plate have revealed several continental or oceanic fracture zones. These fracture zones are often interpreted as Atlantic passive margins, the site of seismic strain rates and active deformation [1] [2] [3]. The most seismically active regions are located along the transform faults [4] [5]. The Atlantic Ocean has extended along these faults, especially in its southern part where there is junction of the oceanic fracture zones with the West African continental margin in the Gulf of Guinea (GOG) and continental faults [6]. The GOG region consists of several tectonic structures which are linked to the pan-African Proterozoic orogeny [7] [8]. The seismic activity inland would be related to the reactivation of Paleozoic and Proterozoic tectonic structures of West Africa in the Quaternary [9] [10] [11]. Kutu et al. [12] show the continental active tectonics seems to be linked with the oceanic transform faults in the Gulf of Guinea. Thus, several continental faults may generate large earthquakes (Mw ≥ 6.0) and can be associated with the Romanche, Chain and Saint Paul transform faults [4].
Bertil [13] studied the intraplate seismicity for the period 1967 to 1991. He used data from Lamto seismic networks in Ivory Coast and neighboring countries. According to his results, earthquakes occur at lithospheric preexisting zones
Figure 1. Geological map of West Africa.
of weakness: continental faults systems and oceanic fracture zones. Although several major historical and recent earthquakes were recorded in the region, the seismicity of West Africa is infrequent according to the geological context. Using the geological and geophysical features, Meghraoui et al. [14] defined from the tectonic regime and seismicity background as well as the geographical situation that West Africa is a stable seismotectonic province. This intraplate area is generally considered as aseismic or having low seismicity. However, the largest magnitude recorded is 6.8 in Ghana. Several other major earthquakes of magnitude 4 to 6 are felt in West Africa and must be taken into consideration for seismic hazard analysis.
The purpose of this work is to assess the seismic hazard level in this province from an incomplete seismic catalog compiled from the available data gathered from previous works and international data centres during the period 1615 to 2021. The seismic parameters are estimated using the computer code describing the probabilistic approaches of Kijko and Sellevoll [15] [16] and Kijko et al. [17].
2. Geological Setting
The geological history of West Africa begins with the Archaean era (about 2.5 billion years old). Its different deformation processes correspond to main geodynamic events that affected this region during successive orogenic cycles. West Africa is characterized by formation of craton and its borders, western (Atlantic passive margin) and eastern (area of the active intracontinental rift system of Central and West Africa).
The West African craton has remained stable since the Liberian orogeny. It was partly reworked during Eburnean and Pan-African cycles [18] [19] [20].
The western and northern edges of the craton contain any evidence of significant deformations undergone during the Variscan orogeny [21]. The Pan-African suture in the East was reactivated [22] [23]. At the east end of the craton, the Apto-Albian corresponds to an important geodynamic event marked by the opening of the Equatorial Atlantic Ocean [24] which was characterized by an initiation of oblique rifting in the GOG [25]. This opening accompanied by the “abortion” of the Benue Trough [26], had repercussions in terms of deformation in the north-eastern part of the craton [27] [28] [29]. This episode of deformation is known as the “Austrian phase”, and is underlined by tectonic structures that follow North/South directions along the Pan-African suture northwest of the Tuareg Shield [22]. The Cenomano-Turonian bears witness to the most important transgression of the Phanerozoic [30]. West Africa was then partially covered by the sea which was moving southwards [31] [32] [33]. The series associated with this transgression are preserved in the eastern part of the craton [34] [35]. Its extension towards the west on the craton is still unknown.
The Archean and Paleoproterozoic domains of craton are located in two broad shields (the Reguibat Shield to North and the Man-Leo Shield to the south). They are separated by the vast intracratonic Precambrian basin of Taoudeni with its sedimentary cover, as well as mobile belts of late Precambrian age, locally affected by a Paleozoic imprint. At the extreme north of the craton is the Anti-Atlas (thrust and fold belt), which represents the northern end of the West African Craton. Just below the Anti-Atlas, is the Tindouf basin whose sedimentary filling is predominantly Paleozoic, with a basal cover of upper Proterozoic age. The western end of the Reguibat Shield consists of the Mauritanides, which are a portion of the Variscan chain in north-west Africa [21] [36]. This complex shows tectonic sheets with E-SE vergence, which extends the fold belt of the Zemmour and then of the Anti-Atlas moving northwards, and whose age of emplacement is estimated between 270 and 330 Ma [21] [36] [37].
West Africa is a seismotectonic area yet poorly known in terms of the current faulting activity, crustal deformation, and their geodynamic causes. However, this region is also seismically active (Figure 2). The geological context is propitious to the presence of active faults that can generate destructive earthquakes [4] [5] [11]. A hazard assessment is therefore required in the area.
3. Data Analysis and Methodology
The historical and instrumental regional seismicity data used in our study was obtained from the seismic catalogs prepared during the studies of intraplate seismicity in West Africa from 1615-1991 [4] [13] and the international data
centres such as IRIS, USGS, BGS and ISC. Ambraseys and Adams [4] listed several historical and instrumental events between 1615 and 1970 without taking into account those that were doubtful. Bertil seismic catalog (1967-1991) was built entirely from the recordings of networks of seismic stations in Ivory Coast and neighboring countries [13]. These two seismic catalogs were supplemented by those of the international data centres (ISC, USGS, BGS and IRIS) to cover a period from 1615 to 2021. The resulting catalog has been homogenized to moment magnitude (Mw). According to the approximations carried out by Musson [38], we also adopted that macroseismic magnitude values quoted by Ambraseys and Adams [4] can be taken as roughly equivalent to Mw. Local magnitude values of Bertil catalog (ML) were converted to moment magnitude through the relation derived by Grünthal et al. [39]:
. (1)
Finally, this seismicity catalog was divided into incomplete historical data and complete observations (Figure 3). The incomplete historical part of catalog covers a period from 1615 to 1795 while the complete part can be subdivided into smaller complete subcatalogs each one with a different time period.
Thus, the earthquake hazard parameters of West Africa are estimated from an algorithm developed by Kijko and implemented in Matlab. The applied procedure provides a maximum likelihood estimate of earthquake hazard parameters, together with their uncertainties, by combining historical and instrumental data. In the description of this analytical method, the moment magnitude (Mw) is reported as m; we have chosen to do this because it will be easier to include it in the formulae than to use the conventional notation Mw.
The theoretical background assume that the historical part contains
largest seismic events, each with a magnitude
(
) such as
and
,
being overall minimum magnitude of the whole catalog. The time span of the historical part of the catalog
can be expressed as the sum of the time intervals
between the historical events. The complete part can be divided into s subcatalogs (s being the number of subcatalogs), each one complete has
earthquakes, each with a magnitude
(
;
Figure 3. Earthquakes magnitude distribution of seismicity of West Africa.
) such as
,
is threshold magnitude of subcatalog with a time span
. The maximum likelihood procedure require that the temporal and magnitude distributions of seismicity for a specified source are respectively described by the Poisson process [40] and the frequency-magnitude Gutenberg-Richter relationship [41].
, (2)
where a and b are parameters, m is the magnitude and N denotes the number of earthquakes of magnitude greater than or equal to the level of completeness magnitude
. The probability that a total of n earthquakes will be observed during the specified time interval t within a given area by the formula:
,
(3)
refers to the Poisson distribution parameter and describes the area-characteristic, mean seismic activity rate of earthquakes with magnitudes greater than or equal to
. The standard seismic hazard analysis procedure requires that earthquake events are independent [42] then the whole catalog must be declustered.
Considering both models (2) and (3), the classic probability density function (PDF) of earthquake magnitude defined by Aki [43]:
for
, (4)
while the cumulative distribution function (CDF) of earthquake magnitude is also given by Aki [43]:
for
, (5)
where the magnitudem is a continuous variable which may take any value greater or equal to than the completeness level
and
. The most seismic-hazard assessment procedures do not integrate the temporal variations of seismicity [40] [44] [45] because it has been found that seismicity is temporally stationary [46]. However, large temporal changes of small to moderate magnitude seismicity are often observed in stable intraplate regions [47] [48]. The variations in b values do occur regionally, these variations are often correlated with the regional features of seismic activity and may be significant [17]. In probabilistic seismic hazard analysis (PSHA) studies, the b value is crucial for the accuracy of the final evaluation results. The probabilistic model described by Kijko et al. [17] use Bayesian formalism to account the uncertainty in the parameters λ and b by considering them as random variables [49] as well as the uncertainty associated with the earthquake-occurrence models. The seismicity parameters are distributed according to the gamma distribution. Thus, the probability density function (PDF) of gamma distribution is defined as follows [50]:
with
. (6)
In Equation (6),
denotes gamma function [50]:
,
(7)
where p and q are related to the mean
(
) and the variance
(
) of the distribution.
The coefficient of variation expresses the uncertainty related to the variation of a given parameter relative to its mean is given by:
The Bayesian CDF of earthquake magnitudes derived from Equations (5) and (6) will take the form of the Bayesian Exponential-Gamma distribution, as described by [51]:
. (8)
In which
,
with
and
being respectively the mean value and the standard deviation of the parameter
.
being a normalizing coefficient of the form:
(9)
The poisson distribution (3) combined with the gamma distribution (6) will become the Poisson-gamma Bayesian distribution:
, (10)
where
and
are the parameters of gamma distribution (6),
the mean value of the activity rate
.
As
and
Equations (8) and (10) may respectively be rewritten in the form [17]:
(11)
(12)
In this case, the normalizing coefficient takes the form:
, (13)
and (
,
) are calculated by applying the maximum likelihood procedure.
Using Equations (11) and (13), the Bayesian PDF of earthquake magnitudes are defined as [17]:
(14)
The maximum likelihood estimators of the parameters
and
respectively denoted
and
are the values of
and
that maximizes the likelihood function
for a given maximum area-characteristic earthquake magnitude
. Thus, maximization of the likelihood function is obtained by
solving the system of two equations
and
, in which
. The uncertainties of the parameters
and
expressed as a formal estimate of the variance of
are given from the equations describing the approximate variance-covariance matrix of vector
such as [52]:
, (15)
in which derivatives are calculated at the point
and
. Equation (15) provides a good approximation of the event variance-covariance matrix for sufficiently largen. It should be noted that the formulation (15) does not provide estimation of the standard error of the estimator of
denoted
.
The likelihood function
calculated from all available data is the product likelihood functions, based on the historical and complete parts of the catalog, as
(16)
In Equation (16), the likelihood function
of the parameters
and
for the historical events is [17]:
,
in which
which denotes the PDF of the largest earthquake magnitude within the time period
takes the form [17]:
,
with
, the mean activity rate for earthquakes with magnitudes
and
as the mean activity rate corresponding to the magnitude value
. The functions
and
are respectively the CDF and PDF of earthquake magnitudes obtained in Equations (11) and (14). According to Equation (16), the likelihood function
based on all s complete subcatalogs is defined as [17]:
, with (
)
where
,
and
,
and
are the likelihood functions of each complete subcatalog i (
) where
is the mean earthquake activity rate, corresponding to the magnitude level of completeness
.
is defined in Equation (11).
The maximum possible magnitude
of earthquake occurring within a specified source is calculated using Bayesian formalism [53]:
. (17)
In Equation (17),
denotes the complementary incomplete gamma function [50], p and q are defined as the Bayesian distribution parameters,
and
with
being a normalizing coefficient defined in Equation (9). The variance of the Bayesian estimator
corresponding to its uncertainty can be approximated as [53]:
(18)
The procedure also accepts time gaps (periods during there is no seismic record) and takes into account the uncertainty in the magnitude of seismic events and treats them as random errors following a Gaussian distribution with zero mean and standard deviation σ [54]. Thus, the seismicity parameters are calculated simultaneously by an iterative scheme in the seismic code software developed by Kijko.
4. Results
The outlined technique was applied to West Africa to estimate the earthquake hazard parameters. Table 1 lists the incomplete data of the largest historical events while the complete part of the catalog was subdivided into four subcatalogs, each with a different level of completeness (Table 2).
For each subcatalog, the level of completeness was determined using the maximum curvature technique [55]. The magnitude of completeness
estimate is based on the departure of the Gutenberg–Richter frequency–magnitude
Table 1. A part of historical catalog of West Africa.
Table 2. A summary of the complete parts of the West Africa Catalog.
distribution (Equation (2)) from its linear trend. Table 2 summarizes the dates, completeness magnitude
and the number of events above
as well as the standard error of earthquake magnitude determination in each subcatalog. The frequency–magnitude plots for each subcatalog are presented in Figure 4.
It was assumed that earthquake magnitudes for the incomplete part of the catalog were determined with a standard error equal to 0.3 magnitude units and those of the complete part are assumed to have a magnitude standard error of 0.2 or 0.1 (Table 2). The temporal variation and the effects of swarms included in the estimation process of parameters
and
, were equal to 0.5. The possibility of including any additional information for the case of poor quality seismic catalog is required in the estimation of average values of seismic hazard parameters for West Africa. Kijko and Graham [56] suggest therefore that the value of
and its standard deviation
obtained from the maximization of likelihood function (16) are roughly considered to be the same as for another area of similar seismogenic features. Thus, the chosen arbitrary values of
and standard error of
are 1.0 and 0.1 respectively because the b value of a stable continental region is generally close to 1.0 [57] [58] [59]. The value of
parameter will also help to stabilize the results.
The obtained estimates of b value, mean activity rate occurrence
, and area-characteristic maximum possible magnitude
for West Africa, together with their standard deviations are given in Table 3.
The seismic hazard curves are also plotted; Figure 5 show the return period of magnitude event
, Figure 6 describes the annual probability of exceeding a given magnitude at least once and finally Figure 7, the probabilities of
Figure 4. Graphs showing the frequency-magnitude distributions of each seismic subcatalog of West Africa. The magnitude bin is 0.2. The level of completeness represents the magnitude at which the curve deviates from its linear trend.
Figure 5. The mean return periods of magnitude events from West Africa.
Figure 6. The probability that a given magnitude will be exceeded at least once in any year.
Figure 7. Graphs showing the probability of exceedance of earthquake events from West Africa for the next 25, 50 and 100 years.
Table 3. Earthquake hazard parameters of West Africa.
exceedance of earthquake events for the next 25, 50 and 100 years. Each graph also provides the calculated level of confidence for the calculated values. The main results show that there is about 50% of probability that at least one event will exceed the magnitude Mw = 6.5 in the next 50 years and the return period of earthquake event with moment magnitude Mw ≥ 6.5 is at least one century.
5. Discussions
Seismic hazard assessment in West Africa is rather difficult. The informations about the historical seismicity are very incomplete. Instrumental monitoring of seismicity in West Africa started from 1914 onwards [60] with the installation of seismic networks for different local studies [13] [61] [62] [63] [64]. Nevertheless, very few data are available and some testimonies were scarce and imprecise. The first formal seismological bulletins of West Africa for a significant period was published by Ambraseys and Adams [4]. They listed the largest earthquakes with a rough moment magnitude estimate (Mw) recorded between 1615 and 1970. It was supplemented by the database built from Lamto seismic network of Ivory Coast. Bertil catalog is assumed to be the most complete of the region during 26 years. However, the location of earthquake events are not well constrained for epicentral distances greater than 600 km and the most of earthquake events with Mw ≤ 2.0 are near to the Lamto seismic network. The seismic events of Mw ≥ 3.0 have been detected in neighboring countries and GOG [13]. In spite of these constraints and the choice done by Musson [13] to do not include them in the seismicity of West Africa, we considered these earthquake events and extended the time period (1615-2021) in order to obtain an accurate seismic hazard assessment of West Africa. We can explain it simply because the work carried out by Bertil [13] remains unpublished. In addition, seismicity monitoring of West Africa is very intermittent over long period and there is no available formal seismological bulletin published recently to improve knowledge of low seismicity in this moment. Thus, the technique of Kijko et al. [17] which use the incomplete historical data mixed with recent complete observations, is adequated to apply to the West Africa dataset.
The seismicity parameters calculated give a b value (Table 3) of about 0.83. This b value approaching unity supports that the seismicity of West Africa is associated with the tectonic areas. According to Agrawal [65], that tectonic earthquakes are characterized by the b value from 0.5 to 1.5 and are more frequently around 1.0 for seismically active regions. The b value depends on structural heterogeneity and stress distributions in space [66] [67]. In a stable intraplate area, the b value tend to be lower than those for plate boundary zones, which are usually around 1. One possible explanation for this observation is that within the plate, the rocks are not as intensively fractured [46]. Normally, there is a greater abundance of fractures in a plate boundary which could lead to a larger ratio of smaller events to larger events (large b values). Freeth [68] considered boundary plate tectonic forces to be responsible for continental tectonics in West Africa and suggested therefore that extensional stresses in the plate boundary could make some continental faults active. According to Talwani [69], Intraplate Earthquakes (IPEs) occur within the continental zone when the magnitude of local stresses become comparable with the regional stress (boundary plate). In addition, the strain rate is relatively very low in the western side of the African plate [70] [71]. Then, IPEs are associated with larger stress drops therefore large earthquakes should have larger return period. Thus, the hazard curves obtained (Figure 7) show that earthquake event with magnitude greater than 6.5 can probably occur every century in West Africa.
The maximum observed magnitude (6.8) in the West Africa catalog fall into the southern part of Ghana. This magnitude had previously been used to assess magnitude of the largest possible earthquake for Ghana [72]. They used the deterministic method based on empirical formula of Gupta [73]. Normally, the maximum credible magnitude is assumed to occur at the minimum possible distance from the site, without providing any indication of the likelihood of its occurrence during a specified exposure period. In addition, it is very difficult to identify individual active fault zones because of their rarity and the efforts to study them have been limited. In this case, various tectonic and fault rupture parameters such as length, area and dislocation, cannot be available. Thus, regional maximum magnitude can only be based on the seismicity, and geological, geophysical and structural features available in the West Africa. The seismological information focuses on the maximum observed magnitude and statistical analysis of the available seismic catalog while the geological information is used to identify different tectonic features. Here, we are considered West Africa as a seismogenic province so we are not need to subdivide into seismic sources. According to Gupta [73], the probabilistic approaches are able to give the results more reliable. Thereby, we are chosen to apply the probabilistic procedure developed by Kijko et al. [17] to calculate the maximum regional magnitude in West Africa because it is used for the areas with both low and high seismicity, including cases where the catalogs are incomplete. This procedure also integrates the effects of all the earthquakes of all possible magnitudes and at all significant distances from the site during a specified exposure time, taking into account the associated uncertainties and the uncertainty of earthquake-occurrence models. In PSHA,
plays an important role in assessment of long return periods [74] therefore its value should not be overestimated.
6. Conclusion
We are focused our study on the estimation of earthquake hazard parameters of West Africa. The sismotectonic analysis in the region showed the occurrence of strong earthquakes, some of which are associated with offshore faults closer to the coastline [11] [12] [14]. Unfortunately, West Africa is not well covered by seismic stations permanently operating and the earthquake catalogs are very incomplete. In this preliminary study, the probabilistic approach developed by Kijko et al. [17] has been successfully tested. The maximum regional magnitude is around 7 and hazard curves are plotted for return periods of 25, 50 and 100 years. As a result, we observed the return period of the largest earthquake is around one century but we also have 50% of chance to have an earthquake of magnitude 6 occur at least once in West Africa in the next 25 years. These results we lead to take into account the main seismic source zones of West Africa most affected by large earthquakes in the estimation of the level seismic hazard in this part of Africa.
Acknowledgements
This work was carried out at the Institute for Geological and Mining Research, Cameroon. The authors sincerely thank Professor A. Kijko for giving us access to its seismic code software and we also thank the Department of Physics of the University of Yaounde I for their prompt collaboration.