Abstract
GEOMATICS, NATURAL HAZARDS AND RISK 2023, VOL. 14, NO. 1, 2197104 https://doi.org/10.1080/19475705.2023.2197104 Seismicity pattern of African regions from 1964–2022: b-value and energy mapping approach Alemayehu Letamo, Kavitha B and Tezeswi TP Department of Civil Engineering, National Institute of Technology Warangal, Telangana, India ABSTRACT ARTICLE HISTORY Received 12 September 2022 The African continent is one of the few tectonic plates where Accepted 24 March 2023 very few earthquakes have been reported. Sporadic poor data have jeopardized efforts to assess seismic hazards. In light of the KEYWORDS tectonic setting, a review of the continent’s seismicity is crucial Earthquake; b-values; for predicting seismic hazards. For this study, seismicity data from energy mapping; seismicity USGS and ISC catalogs were collected from 1964–2022, in the lati- assessment; African plate; o o o o tude range of 40 S-40 N and longitude 20 W-52 E, from the sur- subduction zone face to 700 km deep. Gutenberg’s parameters as well as the spatial variation of seismic energy were then analyzed as an assessment of Africa’s seismicity. The results show that NE Africa is at a higher stress level in terms of energy release than NW Africa. The Richter and Guttenberg’s constants of a¼ 5.61, b¼ 0.59; a ¼ 6.55, b¼ 0.77; a¼ 6.44, b¼ 0.78; a¼ 7.41, b¼ 0.92; a¼ 5.94, b¼ 0.86; a¼ 5.97, b¼ 0.77 with the corresponding mag- nitude of completeness, mc, observed to be 3.7, 4.1, 4.4, 4.7, 4.8 and 4.6 were found for regions NW, NE, southern, eastern, west- ern and central Africa, respectively. The distribution of focal depths vs. magnitude revealed that most of the seismicity is inherently shallow in its nature with an average focal depth of 21.22 km. In comparison, NW experiences relatively more pro- found events with an average depth of 32.19 km. A substantial rise in focal depth in North Africa is associated with the collision boundary. This information as a constraint on the tectonic plate near convergent boundary could be useful in plate modeling. KEY POLICY HIGHLIGHTS The higher stress observed in the North Africa and East Africa continent is vital for policy makers for sustainable built environment The estimated minimum magnitude cut-off Mw 3.7 reveals that the recording facilities in African continent is now improving. An increased focal depth distribution of events in North Africa infers there is plate convergence and can be used to constrain African tectonic plate for plate modeling. CONTACT Kavitha B kavithab@nitw.ac.in 2023 The Author(s). Published by Informa UK Limited, trading as Taylor & Francis Group. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/ licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. The terms on which this article has been published allow the posting of the Accepted Manuscript in a repository by the author(s) or with their consent. 2 A. LETAMO ET AL. 1. Tectonics of African continent The seismic activity centers, such as rift zones, thrust and fold mountain belts, transform faults, and volcanic areas significantly influence Africa’s seismotectonics. The most notice- able regions of continuing tectonic deformation with significant seismic activity are the North African thrust and fold belt and the East African Rift system (Meghraoui and Pondrelli 2012; Misra and Mukherjee 2015). However, the most seismically active zones on the continent include the Cameroon Volcanic Line, the Congo Basin and the Plateau in Southern Africa (Meghraoui and IGCP-601 Working Group 2016). The Atlas orogenic system, which includes the High Atlas, Tell-Atlas, Anti-Atlas, and Sahara Atlas, is where recent deformation and fault zones are concentrated. However, they are sparse to the east of Sidra (Great Sirt) Bay, where faults are repre- sented pass through the Jebel El-Ahdar Mountains. This region of North Africa is part of the Alpine-Himalayan belt. According to estimates for the Neogene over the Middle and High Atlas Mountains, the Atlas deformation contributed significantly to the development of the plate boundary, accounting for 20–40% of the convergence between Africa and Europe (Gomez et al. 2000). Gibraltar-Sisly and the Atlas Mountains are notable collisional plate boundaries, which can be considered typical tectonic features along the Africa-Eurasia plate boundary. This area experienced mod- erate to large earthquakes, with the most notable event being the El Asnam Ms 7.3 earthquake on 10 October 1980, which was linked to a NW dipping, developing thrust fault (Philip and Meghraoui 1983; Yielding et al. 1989). It can be observed from Figure 1 that the fault system of NA where a complex fault system can be seen along the converging Africa-Eurasia plate boundary in the chain of Atlas Mountains. The Sub-Saharan African region’s most notable active continental rift system is the East African Rift System (EARS), which is also known as the East African Rift System. This divergent plate boundary, which separates the Nubian and Somalian plates, runs N-S across eastern Africa. The eastern and western branches of the EARS flank the Tanzanian craton, and the rift valleys form two significant lines. The Mozambique Channel has a third, southeast branch that finally connects to the Southwest Indian Ocean Ridge. The biggest seismicity recorded on the African Figure 1. North African fault system (Adapted from GEM Global Active Faults Database). GEOMATICS, NATURAL HAZARDS AND RISK 3 continent was caused by the western branch of the EARS, which originated around 25 million years ago (Macgregor 2015). According to elastic rebound theory (Reid 1911), earthquakes are caused by the release of accumulated energy within fault raptures of rocks during sudden move- ments of tectonic plates. The distribution and mechanism of active faults, which are the source of significant seismicity, are crucial in seismicity studies. The map of active faults was recently compiled by Styron and Pagani (2020) under the project ‘The GEM Global Active Faults Database’. Figure 1 shows the mechanism of faulting for North Africa. It demonstrates that convergent northern Africa is subject to reverse faulting, and the formation of sizable thrust systems and orogenic belts (the Atlas and Betic/Rif chains) is primarily to responsible for the continent’s seismic activity. Diverging borders, however, show typical faulting (dextral and sinistral). It can be observed from Figure 2, that the East African rift system is now experiencing typical faulting in that it is a developing divergent tectonic plate boundary. Large earth- quakes of various sizes have been recorded in historical and current times in the region around the triple intersection (Afar triangle). The area is dominated by normal faulting, according to the region’s surface geology and the focal mechanisms of earth- quakes (e.g. Kebede and Kulhanek 1991; Ayele 2002), which is in line with current worldwide analyses of active faults by Styron and Pagani (2020). Figure 2. The East African rift system with major active faults shown green (adapted from Chorowicz 2005). 4 A. LETAMO ET AL. The spatial distribution of faults within the African plate is covered in different lit- eratures (e.g. Morley 1999; Skobelev et al. 2004; Macgregor 2015; Poggi et al. 2020). 2. Introduction Earthquake risks have been significant in African countries, specifically along and near plate boundaries and away from plate boundaries in the intraplate regions. Even though the African plate is widely considered a stable seismicity region, the presence of potentially damaging earthquakes has occurred throughout history. For instance, the 1960 earthquake in Agadir, with a magnitude of 5.8 at a depth of 15 km, caused estimated fatalities of 12,000–1500 (Utsu 2002). This demonstrates that even moder- ate events can lead to significant loss of lives and properties. During the last few dec- ades, several earthquakes caused significant losses in Africa (Ambraseys and Adams 1991; Fenton and Bommer 2006). Moreover, with the increasing development of urban areas in Africa, especially in sub-Saharan Africa with modern houses (Tusting et al. 2019), there is an increasing risk of major earthquakes. Recent studies have shown that the African plate is home to many significant and damaging earthquake tremors (Table 1). In general, Africa has been considered as a sta- ble continent. Nevertheless, a large portion of modern-day Africa exhibits a variety of deformations, from compression and Orogenesis in NW Africa to compression across the Congo Basin to continental rifting in conjunction with the East African Rift System to the transition from continental rifting to active oceanic spreading in Afar, the Red Sea, and the Gulf of Aden (Craig et al. 2011). This paper examines the seismicity of the African continent by dividing the whole continent into five regions depending on the proximity of the extensional boundary for West Africa, East Africa, and Southern Africa, the collisional boundary for Northern Africa, and the intraplate region for cen- tral Africa. The study area in this study is bounded by the Mediterranean Sea in the north, the Red Sea in the North East, the Indian Ocean in the East, and the Atlantic Ocean in the West. It is divided almost equally in half by the Equator. Here, the classification of what to be stated is based on tectonic forces and sub-litho- spheric activities. For instance, North Africa is much affected by the Africa-Eurasia col- lision, whereas Eastern Africa is widely affected by the GEARS and the accompanying Table 1. Some significant earthquakes observed in Africa in recent decades. Region Country Time Magnitude Mw References Bukavu D.R. Congo 2008 6.0 Meghraoui and IGCP-601 Working Group (2016) north-western Guinea 1983 6.3 Irinyemi (2021) Guinea Zemmouri- Algeria 2003 6.8 Meghraoui and IGCP-601 Boumerdes Working Group (2016) Machaze Mozambique 2006 7.0 Fenton and Bommer (2006) Juba South Sudan 1990 7.1 Meghraoui and IGCP-601 Working Group (2016) EL Asnam Algeria 1980 7.3 Meghraoui and IGCP-601 Working Group (2016) GEOMATICS, NATURAL HAZARDS AND RISK 5 tensional stresses emanating from the elevated topography. Whilst some part of Africa (e.g. Congo) is affected by compressive stresses from ridge push from the Atlantic Ocean (Ayele 2002). On the other hand, West Africa is a relatively stable region with small far-field stresses from the ridge push in the Atlantic Ocean (Zoback 1992). In comparison, southern Africa is characterized by tensional stresses similar to Eastern Africa, with some shallow earthquakes from human-induced seismicity. Due to the interaction between seismically active Africa-Eurasia plate boundaries and faults, North Africa is characterized by active seismicity on the continent. The region covers countries such as Tunisia, Mauritania, West Sahara, Morocco, Algeria, Libya and Egypt. Along the plate boundary between Africa and Eurasia, the NW is known for recent major thrusts and strike-slip earthquakes. However, owing to two different tectonic regimes, the collision between Africa and Eurasian plates and the rifting of the Red Sea, the north-east- ern area, which includes Libya and Egypt, is characterized by uncommon big and moderate regular faulting earthquakes (Meghraoui and IGCP-601 Working Group 2016). Although numerous significant historical and contemporary earthquakes have devas- tated the area, West Africa, which spans from Nigeria to Senegal, is primarily regarded as a stable section of Africa. A notable illustration of an instance of intraplate seismicity that took place on a stable West African craton is the Guinea earthquake of 22 December 1983, with a magnitude of Mw 6.4. The portions of the CA with moderate to low seismicity zones include Angola, Chad, the Congo basin, and the Cameroun Volcanic Line (CVL). One of the world’s most prominent and least researched intracon- tinental basins, the Congo Intracratonic Basin is bordered to the east and west, respect- ively, by the Cenozoic East African Rift System (EARS) and the Atlantic passive continental margin. The primary deformation in both cases was of the thrust type, which is unusual for the African Plate. A compressive stress environment might be pro- duced by the mid-Atlantic Ridge and the EARS plate boundary pressures (Ayele 2002). A continuous review of the continent’s seismicity is essential for predicting seismic hazards within the tectonic setting. Because large earthquakes usually have a signifi- cant return period (e.g. for example Guinea, West Africa, the return period could reach 100 years for magnitudes that exceed Mw 6.0; Irinyemi 2021). Seismic studies need to be updated so that they can be used as input for seismic hazard assessment. Africa’s poorly-built infrastructures and the region’s vulnerability make it very per- tinent to understand the potential seismicity in the region. Even moderate-sized earthquakes can prove that they can cause significant loss of lives and properties, so thorough mapping is critical. This work assesses the potential seismicity of the region. The Assessment of the seismicity of Africa in terms of Gutenberg’s parameter vari- ation in the six regions of Africa for the period 1964–2022 as well as the investigation of the planar variation of seismic energy in light of seismicity assessment in the Continent are considered as the main objectives of the present study. 3. Seismicity data 3.1. Sources Earthquake data of six regions of Africa and its adjoining regions were compiled in 0 0 0 0 the range of 40 S–40 N Latitude and 20 W–52 E Longitude for magnitude range 3– 6 A. LETAMO ET AL. 7.5, depth 0–700 km, magnitude type m , M , M , M , and M . For the present study, b s w L D 93,691 earthquake events were collected from both ISC Catalogue and USGS cata- logues for 1964–2022. Based on the USGS classification, data depths up to 700 km are used to incorpor- ate shallow, intermediate, and deep earthquakes. Shallow events range from 0–70 km, intermediate ones from 70–300 km, and deep ones from 300–700 km. Classification of seismic provinces is based on active tectonics and present-day continental deform- ation (Meghraoui and IGCP-601 Working Group 2016). The geographical limits of the study area and the epicentral distribution of the earthquake events are shown in Figure 3(a) and Figure 3(b), respectively. The date of the earthquake, the epicentral coordinates, and the magnitude of the earthquake in different scales, such as moment magnitude (M ), surface-wave magni- tude (M ), body-wave magnitude (m ), duration magnitude (M ), local magnitude s b D (M ), moment magnitude (M ), depth, and time of events, are all included in the L w earthquake information collected for each event in the database. Figure 4 presents the distribution of earthquakes. 3.2. Conversion of different magnitudes into M The ISC catalog consists of heterogeneous magnitude scales reported from different agencies, thus necessitating homogenization with regression equations. In this study, empirical relations were used to convert from body wave magnitude (mb) both in Figure 3. (a) The six divisions of the study area. (b) Epicentral distribution of the declustered earth- quake events. GEOMATICS, NATURAL HAZARDS AND RISK 7 Figure 4. Distribution of earthquake events by depth. The distribution of intermediate to deep events is shown in the African plate’s northern part, while the remaining continent exhibits a shal- low depth distribution. USGS and ISC, from surface-wave magnitude (Ms) local magnitude (ML) and Duration magnitude (MD) to corresponding moment magnitude (Mw). For conversion of different magnitudes to Mw (Scordilis 2006), for global data including Africa; M ¼ 0:85ð60:04Þm þ 1:03ð60:23Þ; 3:5 m 6:2 (1) w b b Then the relation in Equation (1) can be used to estimate the moment magni- tude (M ). Conversion M M s w M ¼ 0:67ð60:05ÞM þ 2:07ð60:03Þ; 3:0 M 6:1 (2) W s S 8 A. LETAMO ET AL. Figure 5. Seismicity distribution of Africa continent in G-R. The raw data were compiled from the ISC and USGS catalogues and then declustered. M ¼ 0:99ð60:02ÞM þ 0:08ð60:13Þ; 6:2 M 8:2 (3) W s S Conversion M M L w M ¼ 0:722M þ 0:743; M > 2 , for intraplate domains (4) W L L Conversion of M M according to Akkar and Bommer (2010) D w M ¼ 0:764ðÞ 60:04 M þ 1:379ðÞ 60:2 ; 3:7 M 6:0 (5) w D D Conversion from M m to M according to Weatherill et al. (2016) sðISCÞ bðISCÞ < 0:616M þ 2:369; M 6 S S M ¼ 0:994M þ 0:1; M > 6 (6) w S S 1:084m 0:142; m < 6:5 b b Conversion from M m to Mw, according to Weatherill et al. (2016) sðNEICÞ bðNEICÞ 0:723M þ 1:798; M < 6:5 S S M ¼ 1:005M 0:026; M > 6:5 (7) w S S 1:159m 0:659; m > 6:5 b b 3.3. Declustering of catalogue Prior to computing seismicity, removing dependent events from the catalog are the first step in assessing earthquake hazards. Therefore, only shocks that trigger inde- pendent behaviour should be considered in the analysis. This can be achieved by removing dependent events (foreshocks and aftershocks) from the updated catalog by specifying a window both in the space and time domain. GEOMATICS, NATURAL HAZARDS AND RISK 9 Finding mainshocks and aftershocks is easy using windowing methods. The follow- ing shocks are classified as aftershocks for each earthquake in the catalog of magni- tude M if they happen within a specific time window, t(M), and a certain distance, d(M). When the strongest earthquake happens later in the series, foreshocks are clas- sified similarly to aftershocks, i.e. as aftershocks. The size of the most significant shock in a series determines how the time-space windows are reset. Even though there are several window sizes in the computation of declustering of catalogs, Gardner and Knopoff (1974) method is a widely used method owing to its ease of applicability. According to Gardner and Knopoff (1974), the window sizes in the space and time domain are shown below (Equation (9)). 0:032M þ2:7389 10 , if M 6:5 0:1238M þ0:983 w dðkmÞ¼ 10 tðdaysÞ¼ (9) 0:5409M 0:547 10 , if M < 6:5 Hence, in this study, ZMAP (Wiemer 2001) with the choice of Gardner and Knopoff (1974) algorithm is used to decluster our catalog. 3.4. Catalogue completeness In seismic risk analysis and prediction, the Gutenberg-Richter b-value is an essential parameter. However, sometimes the FMD could be erroneous and lead to bias due to the incompleteness of the instrumental catalogue. This could happen because some magnitudes may not be felt and remain undetected by the sensors, thus affecting the reliability of the empirical relation. Therefore, it is necessary to set the local cut-off point, M , in such a way as to maximize the reliability of the G-R relation. Furthermore, the Mc value depends on many general factors, such as the number and quality of the sensors, their geometrical configuration, and the environmental noise level. Hence, accurate computation of the b-value is crucial to the effectiveness of the method. In this study, a maximum curvature (MAXC) method proposed by Wiemer and Wyss (2000) was used, where the magnitude of completeness is obtained by maximiz- ing the best-fit equation of the frequency magnitude relation. The estimated com- pleteness magnitude of the catalogue is 4.4 which is clearly shown from Figure 5. The G-R parameters were computed in the maximum likelihood method (Aki 1965) with consideration of cut-off magnitude. 3.5. Review of G-R parameters Here, the review of some representative G-R parameters with the main findings for the six regions of Africa is presented for different time ranges. Gutenberg’s constants, a ¼ 5.75 and b¼ 1.23, were discovered by Alabi et al. (2012) in their investigation of the seismicity of southern Africa. Their findings showed that southern Africa is less seismic. The presence of a moderate earthquake, however, could be problematic in the future. The likelihood of significant earthquakes occurring in this area shortly was low. However, a significant earthquake could 10 A. LETAMO ET AL. Figure 6. Temporal distribution of cum. no. of events. happen due to the uncertainties surrounding earthquake projections, particularly in intraplate locations. Kadiri and Kijko (2021) studied the seismicity of West Africa in G-R relations. According to their finding, West Africa lies in the tectonic seismicity with a value of 4.1 and the corresponding b-values of 0.77. On the other hand, Irinyemi (2021) calcu- lated the G-R parameters to determine hazard levels in west Africa and found that the estimated hazard levels were high in the Palaeozoic area of Guinea with a uni- form b-value of 0.70 ± 0.12 and a ¼ 4.2. From 2013 up to 2016, Asefa and Ayele (2021) studied the space-time distribution of earthquakes in the East African Rift System (EARS). The average b- and a-values of their study events over this period were 1.01 and 6.5, respectively, according to the frequency- magnitude distribution (FMD) of those events. The b-value result suggests a region under relatively high stress, with earthquakes of tectonic origin predominating. While Lamessa o o o o et al. (2019), in their study, for the latitude range of (4 N-20 N) and longitude (34 E-48 E) for the period 1960 to 2016, found a value of 6.9 and b¼ 1.03. Moreover, Poggi et al. (2017) analysed the seismicity parameters of the Gulf of Aden (Near the Horn of Africa) for the time range of 1904–2013 using the MLE approach and found b¼ 1.02, a¼ 5.38. From October 1997 to December 2013, Ali (2016) examined the seismicity of Egypt and surrounding regions using data from the Egyptian National Seismological Network (ENSN) that was limited to latitude (22 N-34 N) and longitude (25 E- 36 E). His b-values ranged between 0.50 and 1.19. The Levant Aqaba with b value¼ 0.96 and the Cyprian Arc with b value¼ 1.19 have the highest b values, whereas the Suez-Cairo-Alexandria fault with b-value¼ 0.812 and the active Kalabsha fault with b value¼ 0.554 have the lowest b values. At the same time, Poggi et al. (2020) looked into the b-values as part of a study project on hazard assessment in northwestern Africa (High/Middle Atlas). As per them, the Tell Atlas’ b-value was 0.98, whereas the Middle Atlas’ was 1.1 for the G-R relation. The G-R parameters for some repre- sentative areas of Africa are presented in Table 2. GEOMATICS, NATURAL HAZARDS AND RISK 11 Table 2. Review different b-value studies for the six regions of Africa. Area Region a-value b-value References Western Rift CA – 0.773 Mavonga and Durrheim (2009) South Sudan EA 3.579 0.92 Uluka et al. (2020) Egypt NEA – 0.5–1.19 Ali (2016) High/Middle Atlas NA 3.92–4.45 1.10 Poggi et al. (2020) Rif/Tell Atlas NA 3.72–4.47 0.98 Poggi et al. (2020) EARS EA 6.5 1.01 Asefa and Ayele (2021) Southern-Ethiopian Rift EA 6.9 1.03 Lamessa et al. (2019) South Africa SA 5.75 1.23 Alabi et al. (2012) Gulf of Aden EA 5.38 1.02 Poggi et al. (2017) Mweru–South Katanga CA 4.05 0.99 Poggi et al. (2017) Western Africa WA 4.1 1.01 Kadiri & Kijko (2021) Guinea WA 4.2 1.03 Irinyemi (2021) Figure 7. Magnitude versus time distribution of the study area. 3.6. Final catalogue After removing duplicate events, 93,691 events were eligible for the final screening, which were then declustered with the Gardner and Knopoff (1974) algorithm. The remaining events were then filtered with a minimum cut-off value, and the final cata- log with 13,019 events in moment magnitude is selected for analysis. The temporal variation of the cumulative number of events in the catalogue is shown in Figure 6. From the graph, the slope increases for the years 2004–2010, indicating the seismicity rate is higher in this time interval. It can be noted form Figure 7 that a dense distribution of large-magnitude earth- quake events occurred in this time interval. 4. Seismicity parameters The relationship between the size distribution of earthquakes and magnitude was dis- covered by Gutenberg and Richter (1944) and takes the form of the power law (Equation (10)). 12 A. LETAMO ET AL. Figure 8. The relative proportion of energy intensity in J/km for the six regions of Africa. log NðmÞ¼ a bm (10) where Nis the number of earthquakes with a magnitude greater than or equal to m, and a and b are constants that describe the frequency-magnitude distribution (FMD). Parameter a describes the productivity of the volume or quantity of earthquakes within the study area, and b the slope of the spatial mapping of frequency-magnitude distribution (FMD), describes the relative size distribution of events. The value of b is generally 1 or close to 1 for seismicity of tectonic origin (Mul € ler & Jokat, 2000). A value of b greater than 1 indicates predominantly small earthquakes, while a smaller value of b indicates the dominance of more significant earthquakes concerning smaller events. Low-magnitude earthquake frequency declines, especially in seismically active regions, signify stress build-up. Lower ‘b’ readings before the great earthquakes reflect this (Gadkari and Mukherjee 2023). Moreover, high b-values in the volcanic compartments could explain the high seismicity in the region (Sanchez et al. 1995). In contrast, the regions with low b values imply the considerable differential stress of the Earth’s crust, thereby pointing toward the end of the seismic cycle (Schorlemmer et al. 2005). Researchers (e.g. Utsu, 1965;Scholz 1968) pointed out in their studies that a relation could exist between b-values and stress and/or strain for a given locality. Due to enormous implications in geodynamics and hazard assessment, correct esti- mation of Mc and Gutenberg’s constants is crucial in the GR relation. Generally, there are two approaches to computing the b values: least square fit and the max- imum likelihood methods. However, the least square method is subject to bias because of the assumption that all the data point carries the same weight and the residuals are Gaussian-distributed, which cannot be justified at higher magnitudes (Naylor et al. 2010). To avoid such inconveniences, the maximum likelihood method (MLE) proposed by Aki (1965) was used (Equation (11)). log e 0:434 (11) M ðÞ M 0:1=2 M M þ 0:05 av c av c where M is the mean value of the sample, M is the threshold above which the av c FMD is considered complete, and 0.1 is the bin size. The constant of GR relation in Equation (10), a value is determined from the rela- tion: GEOMATICS, NATURAL HAZARDS AND RISK 13 a ¼ logNðÞ M þ bM (12) c C In this work, b-value, the slope of the G-R relation, and the constant ‘a’ are deter- mined using an MLE method in log scale for absolute magnitudes only. 5. Energy mapping Energy mapping is a critical parameter in analysing stress induced during earthquake tremors in the tectonic plates and the process of corresponding seismic hazard map- ping therein. That is, it could be a measure of the potential for damage to man-made structures. Earthquakes exceeding magnitude M 3.0 are used for energy mapping since most energy budgets dissipated are from major and great earthquakes. In the analysis, the relation suggested for stable continental regions such as the African plate by Kanamori (1977) is used (Equation (13)) log M ¼ 1:5M þ 9:1 (13) 0 w To convert the moment magnitude into the radiated seismic energy, Equation (14) introduced by Kanamori (1977) was used, where the radiated seismic energy is lin- early related to the seismic moment, M E ¼ 0:5 10 M (14) s 0 where Es is the energy radiated in Joules and M is the seismic moment in Nm. Hence, from Equation (13) & Equation (14), energy radiated in Joules can be solved as in Equation (15) ð1:5M þ4:8Þ E ¼ 10 (15) The resulting sum of energy for the regions can be used to represent the stress released and accumulation of stress due to instrumental seismicity within the region from 1964–2022 time. The intensity of the energy release of African continent is shown in Table 3.To enhance the comparison, energy intensity variation is also provided in graphic repre- sentation in Figure 8 for each regions. Table 3. Intensity of seismic energy release from 1964–2022 in the African continent and adjoin- ing regions. It is illustrated in the graph form in Figure 8 for better understanding. 16 2 6 2 9 Region Energy released (J) 10 Area (Km ) 10 Energy (J/km ) NE Africa 2.79 2.87 9.74 E. Africa 3.26 5.19 6.29 NW Africa 3.26 6.34 5.14 C. Africa 0.99 6.97 1.42 S. Africa 1.41 11.23 1.25 W. Africa 0.32 6.56 0.48 14 A. LETAMO ET AL. 6. Results & discussions Figure 3(b) shows the seismicity pattern of Africa, and the adjacent area follows the triple junction branching into the E and SE parts of Africa, NE Africa, and the Gulf of Eden, with the centre emanating from the Afar region of Ethiopia. It follows the Great East Africa Rift System (GEARS). Apart from most earthquakes, epicenters look clustered in N Africa. The graph in Figure 5 shows the seismicity of the African continent. As shown in the figure, the G-R constants for the continent are a ¼ 7.25, b ¼ 0.77, and Mc ¼ 4.4. For the six tectonic provinces viz, NW Africa, NE Africa, S Africa, E Africa, W Africa, and C. Africa are a ¼ 5.612, b ¼ 0.59; a ¼ 6.55, b ¼ 0.77; a ¼ 6.44, b ¼ 0.78; a ¼ 7.408, b ¼ 0.92; a ¼ 5.94, b ¼ 0.86; a ¼ 5.969, b ¼ 0.77, respectively. The magnitude of completeness, Mc, was observed to be 3.7, 4.1, 4.4, 4.7, 4.8, and 4.6 for NWA, NEA, SA, EA, WA, and CA, respectively. It can be observed from Figure 8 that NE Africa dominates the energy release budget per unit area in the African continent with 40%, while E. Africa and NW Africa constitute 26% and 21%, respectively. The regions such as western, central and southern Africa total up 13% of the energy reservoir. Though Africa is generally con- sidered a stable continent in seismic hazard events, northern Africa and its proxy (eastern Africa) can be regarded as regions of active seismicity. This study also highlights the presence of collisional earthquakes near Eurasia- Africa plate boundaries due to subducting slabs into the mantle. A close observation of the distribution of earthquakes in depth (Figure 4) reveals a significant amount of intermediate and deep earthquakes in the N. Africa, which is attributed to the sub- duction zone (Hellenic), and the collision of Africa with the Eurasia plate. In essence, it is possible to extract information about the tectonic setting of North Africa and the extent to which the cooling slab deeps. A critical observation from the absence of intermediate to deep earthquakes in the EA and SA is that these areas are home to events of volcanic origin and tectonics of extensional basins. Due to the anomalous temperature in, particularly in EA, the crust is thinner, and the seismogenic depths are relatively shallow. Due to the sporadic nature of earthquakes in the intraplate region lying at a remote distance away from the plate boundaries, coupled with insufficient amount of data in seismic recording networks, the study requires viable solutions such as numerical methods to reinforce the efforts being done in seismic hazard assessment programs economically with less effort. 7. Conclusions Compared to other parts of the African continent, seismic stress release per area is comparatively higher in the northern region. More likely, this was caused by the seis- micity caused by collisions between plates in the Eurasia-Africa boundary zone. Furthermore, it may be inferred from the greater b-value that the population of lesser events is higher in Northern Africa. With higher a-values of 7.408 for eastern Africa and 6.55 for NE Africa, respectively, the seismicity rate is higher. While the b-value for the NW region is lower in terms of the seismicity parameter, NW Africa will likely experience substantially higher stresses because of the complicated fault systems GEOMATICS, NATURAL HAZARDS AND RISK 15 in the Atlantic range. Notably, the b-value of 0.92 and the high seismicity rate, or a- value of 7.408, demonstrating that eastern Africa is the region with active seismicity. With a b-value of roughly 1.0, Southern Africa is still significant. Despite being far from plate boundaries, the East African Rift System’s spreading rate (Pule & Saunders, 2009), the region’s elevated topography (Craig et al. 2011), and human- induced activities like mineral mining (Foulger et al. 2018) may have substantially influenced southern Africa’s seismicity. As far as the epicentral distribution is considered, the pattern follows the margins of both collisional boundaries in the north, the combination of collisional-transform boun- daries in NE Africa, and the extensional boundaries forming a tri-branching pattern in the Afar region of Ethiopia extending to Malawi in the far southern African region. The regions average focal depths range from 11.53 km to 32.19 km, with South Africa having the lowest average and northwest Africa having the highest, with NE Africa com- ing in second. N. Africa experiences a disproportionately high number of moderate and profound earthquakes compared to the rest of Africa (Figure 4). This might have been connected to subduction zone collisional boundaries. Instead, with average values of 21.22 km focal depths, the seismicity in the African lithosphere is primarily shallow in depth. In actuality, southern Africa’s average focal depth is substantially lower than the global average. Of course, human-induced activities like mineral mining could have con- tributed to southern Africa’s extremely shallow nature of seismicity. The study revealed significant findings about the tectonic setting of North Africa, most notably the influence zone of collisional boundaries near the Africa-Eurasia plate boundaries. In addition, the information extracted related to the depth of inter- mediate to deep earthquakes (Figure 4) is essential to constrain the model’s boundary condition for developing a mechanistic model of the African tectonic plate. Acknowledgements Authors are grateful for detailed review comments in two rounds provided by Prof. Soumyajit Mukherjee (IIT Bombay) and Dr. Mery Biswas (Presidency University, Kolkata) which were helpful for enhancing the quality of the manuscript. Disclosure statement The authors report no conflict of interest. Data availability statement The seismicity data used in the study are compiled from the USGS website (https://earthquake. usgs.gov/earthquakes/) accessed in 2022 GC and the ISC website (http://www.isc.ac.uk/) accessed the same year. In addition, derived data supporting the findings of this study are available from the corresponding author [Kavitha B] on request. References Akkar S, Bommer JJ. 2010. Empirical equations for the prediction of PGA, PGV, and spectral accelerations in Europe, the Mediterranean region, and the Middle East. Seismol Res Lett. 81(2):195–206. 16 A. LETAMO ET AL. Aki K. 1965. Maximum likelihood estimate of b in the formula log10N ¼ a-bm and its confi- dence limits. Bull Earthq Res. 43:237–239. Alabi AA, Akinyemi OD, Adewale A. 2012. Seismicity pattern in Southern Africa from 1986 to 2009. Earth Sci Res. 2(2):1–10. Ali SM. 2016. Statistical analysis of seismicity in Egypt and its surroundings. Arabian J Geosci. 9(1):1–16. Asefa J, Ayele A. 2021. Seismicity of the East African rift system for the period 2013 to 2016. J Afr Earth Sci. 183:104315. Ayele A. 2002. Active compressional tectonics in central Africa and implications for plate tec- tonic models: evidence from fault mechanism studies of the 1998 earthquakes in the Congo Basin. J Afr Earth Sci. 35(1):45–50. Ambraseys NN, Adams RD. 1991. Reappraisal of major African earthquakes, south of 20 N, 1900–1930. Nat Hazards. 4(4):389–419. Craig TJ, Jackson JA, Priestley K, McKenzie D. 2011. Earthquake distribution patterns in Africa: their relationship to variations in lithospheric and geological structure, and their rheological implications. Geophys J Int. 185(1):403–434. Chorowicz J. 2005. The East African rift system. J Afr Earth Sci. 43(1–3):379–410. Fenton CH, Bommer JJ. 2006. The Mw7 Machaze, Mozambique, earthquake of 23 February 2006. Seismol Res Lett. 77(4):426–439. Foulger GR, Wilson MP, Gluyas JG, Julian BR, Davies RJ. 2018. Global review of human- induced earthquakes. Earth Sci Rev. 178:438–514. Gadkari D, Mukherjee S. 2023. Assessment of efficacy of “b” value as a seismic precursor for select major seismic events. In: mukherjee S, Dasgupta S, Majumdar C, Mandal S, Dasgupta T, editors. Handbook of petroleum geoscience: exploration, characterization, and exploit- ation of hydrocarbon reservoirs. Hoboken, NJ: John Wiley & Sons Ltd; p. 1–464. Gardner JK, Knopoff L. 1974. Bulletin of the seismological society of America is the sequence of earthquakes in southern California, with aftershocks removed, Poissonian?. Bull Seismol Soc Am. 64(5):1363–1367. Gomez F, Beauchamp W, Barazangi M. 2000. Role of the Atlas Mountains (northwest Africa) within the African-Eurasian plate-boundary zone. Geol. 28(9):775–778. Gutenberg B, Richter CF. 1944. Frequency of earthquakes in California. Bull Seismol Soc Am. 34(4):185–188. Irinyemi A. 2021. Seismic hazard levels in Guinea. West Africa. 0–15. Kadiri AU, Kijko A. 2021. Seismicity and seismic hazard assessment in West Africa. J Afr Earth Sci. 183:104305. Kanamori H. 1977. The energy release in great earthquakes. J Geophys Res. 82(20):2981–2987. Kebede F, Kulhanek O. 1991. Recent seismicity of the East African rift system and its implica- tions. Phys Earth Planet Inter. 68(3–4):259–273. Lamessa G, Mammo T, K, Raghuvanshi T. 2019. Homogenized earthquake catalog and b-value mapping for Ethiopia and its adjoining regions. Geoenviron Disasters. 6(1):1–24. Macgregor D. 2015. History of the development of the East African rift system: a series of interpreted maps through time. J Afr Earth Sci. 101:232–252. Mavonga T, Durrheim RJ. 2009. Probabilistic seismic hazard assessment for the Democratic Republic of Congo and surrounding areas. S Afr J Geol. 112(3–4):329–342. Meghraoui M, IGCP-601 Working Group. 2016. The seismotectonic map of Africa. Episodes J Int Geosci. 39(1):9–18. Meghraoui M, Pondrelli S. 2012. Active faulting and transpression tectonics along the plate boundary in North Africa. Ann Geophys. 55(5):955. Misra AA, Mukherjee S. 2015. Tectonic inheritance in continental rifts and passive margins. http://link.springer.com/10.1007/978-3-319-20576-2. Morley CK. 1999. Patterns of displacement along large normal faults: implications for basin evolution and fault propagation, based on examples from East Africa. AAPG Bull. 83(4): 613–634. GEOMATICS, NATURAL HAZARDS AND RISK 17 Muller € C, Jokat W. 2000. Seismic evidence for volcanic activity discovered in central Arctic. Eos, Transactions American Geophysical Union, 81(24), 265–269. Naylor M, Orfanogiannaki M, Hart D. 2010. Exploratory data analysis: magnitude, space, and time. Community Online Resource for Statistical Seismicity. 1–42. Philip H, Meghraoui M. 1983. Structural analysis and interpretation of the surface deforma- tions of the El Asnam earthquake of October 10, 1980. Tectonics. 2(1):17–49. Poggi V, Durrheim R, Tuluka GM, Weatherill G, Gee R, Pagani M, Nyblade A, Delvaux D. 2017. Assessing seismic hazard of the East African rift: a pilot study from GEM and Africa array. Bull Earthquake Eng. 15(11):4499–4529. Poggi V, Garcia-Pelaez J, Styron R, Pagani M, Gee R. 2020. A probabilistic seismic hazard model for North Africa. Bull Earthquake Eng. 18(7):2917–2951. Pule T, Saunders I. 2009, September. Recent seismicity in the Mozambique region and its impact/effects on South Africa. In: 11th SAGA Biennial Technical Meeting and Exhibition. EAGE Publications BV; pp. cp-241. Reid HF. 1911. The elastic-rebound theory of earthquakes. Univ. Calif. Publ. Bull. Dept. Geol. 6(19):413–444. Sanchez JJ, Mccausland W, Clifford T, Sanchez JJ. 1995. Spatial variations in the frequency- magnitude distribution of earthquakes at Mount Pinatubo Volcano. Scholz CH. 1968. The frequency-magnitude relation of microfracturing in rock and its relation to earthquakes. Bull Seismol Soc Am. 58(1):399–415. Schorlemmer D, Wiemer S, Wyss M. 2005. Variations in earthquake-size distribution across different stress regimes. Nature. 437(7058):539–542. Scordilis EM. 2006. Empirical global relations converting MS and MB to moment magnitude. J Seismol. 10(2):225–236. Skobelev SF, Hanon M, Klerkx J, Govorova NN, Lukina NV, Kazmin VG. 2004. Active faults in Africa: a review. Tectonophysics. 380(3–4):131–137. Styron R, Pagani M. 2020. The GEM global active faults database. Earthquake Spectra. 36(1_ suppl):160–180. Tusting LS, Bisanzio D, Alabaster G, Cameron E, Cibulskis R, Davies M, Flaxman S, Gibson HS, Knudsen J, Mbogo C, et al. 2019. Mapping changes in housing in sub-Saharan Africa from 2000 to 2015. Nature. 568(7752):391–394. Uluka GEMAT, Ukindula JEL, Urrheim RAJD. 2020. Pure and Applied geophysics seismic hazard assessment of the Democratic Republic of Congo and environs based on the GEM– SSA catalogue and a new seismic source model. Pure Appl Geophys. 177:195–214. Utsu T. 1965. A method for determining the value of b in a formula log n ¼ a ¼ bM showing the magnitude frequency relation for earthquakes. Geophys Bull Hokkaido Univ. 13:99–103. Utsu T. 2002. A list of deadly earthquakes in the world: 1500–2000. In: Lee WK, Kanamori H, Jennings PC, Kisslinger C editors. International handbook of earthquake engineering and seismology. Amsterdam and Boston: Academic Press; p. 691–717. Weatherill GA, Pagani M, Garcia J. 2016. Exploring earthquake databases for the creation of magnitude-homogeneous catalogues: tools for application on a regional and global scale. Geophys J Int. 206(3):1652–1676. Wiemer S. 2001. A software package to analyze seismicity: ZMAP. Seismol Res Lett. 72(3): 373–382. Wiemer S, Wyss M. 2000. Minimum magnitude of completeness in earthquake catalogs: exam- ples from Alaska, the western United States, and Japan. Bull Seismol Soc Am. 90(4):859– Yielding G, Ouyed M, King GCP, Hatzfeld D. 1989. Active tectonics of the Algerian Atlas Mountains—evidence from aftershocks of the 1980 El Asnam earthquake. Geophys J Int. 99(3):761–788. Zoback ML. 1992. First and second order patterns of stresses in the lithosphere: the World Stress Map Project. J Geophys Res. 97:11,703–11,729.
Journal
Geomatics Natural Hazards and Risk
– Taylor & Francis
Published: Dec 31, 2023
Keywords: Earthquake; b-values; energy mapping; seismicity assessment; African plate; subduction zone
