Seismotectonic analysis of the Anker area, Kunene region north-western Namibia

Recently, a swarm of seismic events occurred within the Kunene region (NW of Namibia) with magnitude values going up to ML 4.6. Some of the events alarmed the Anker community and have become of national interest. In order to understand the source of this seismic activity, the Geological Survey of Namibia (GSN) and the Council for Geoscience (CGS) of South Africa installed a temporary network of 10 seismic stations in the Anker area to monitor seismicity for a period of 3 months from 20 June to 24 September 2018. More than 1600 microseismic events were recorded and located within the boundaries of the array, which was roughly 50 km × 50 km. The stations were installed in a 50 km × 50 km radius to have an accurate overlap of the network coverage and also by the spatial distribution of the known seismicity in the study area. The magnitude of the events ranged from ML − 1.1 to 3.6 and occurred at depths of less than 20 km. To improve the location of earthquakes, a new velocity model was developed using high-quality seismic dataset recorded by the temporary network. The seismic data was used to determine a minimum 1-D velocity model, which was used for the relocation of the earthquakes using the double difference method. The relocated solutions provided information to explain the seismotectonics of the region and establish the source of the current seismic swarm. Lineaments were mapped using magnetic data to identify the sources of the seismicity. Focal mechanism analysis of some of the events indicates the area is highly faulted with predominantly normal faulting.

. Of late, a number of these events have occurred in areas that have not shown any sign of stress transfer (e.g. Midzi et al. 2018;Calais et al. 2016). Of note is the 2006 Machaze earthquake which occurred in Mozambique and ruptured a 60 km surface fault and the 2017 Botswana earthquake which occurred in an area that had no history of seismic occurrence. These events pose a major challenge to the estimation of seismic hazard for a region (Labidi et al. 2021). It has been suggested that major earthquakes in SCR regions may not only depend on a long-term tectonic loading under a regional, essentially time invariant, tectonic stress field, but instead they can also be determined by small transient stress changes in a crust close to failure equilibrium (Calais et al. 2016). The background tectonic stress field in the lithosphere has been observed to be effectively modulated by very small stress changes and strain rates ~ 100 Pa/year (100 times less than plate motion rates). Their driving mechanism may also be attributed to external, non-tectonic origin thereby triggering some localised earthquakes (Luttrell et al. 2007;Luttrell and Sandwell, 2010). Regardless of the specific transient stress change that brings a fault to failure, the resulting rupture mechanism will however be consistent with the regional tectonic stress field as was observed in many other studies (e.g. Manzunzu et al. 2017;Midzi et al. 2018).
Seismicity in Namibia is generally low with exceptional moderate magnitude events that occur occasionally (Mangongolo and Hutchins, 2008). The recent occurrence of five moderate earthquakes recorded by the Namibian Seismological Network with magnitudes (4.0 < M L < 4.6) in May 2018 in the Anker area, north-western Namibia raised concerns on the possible source of seismic events, whether the events could be precursors to a major earthquake, and their impact on the community. However, the region has limited both seismic stations and detailed geological information to identify the causative source of these events. At the time of the study, there was no functional seismic station in the Kunene region and generally, there is insufficient earthquake records in Namibia to determine the source of the seismic swarm. The events were felt widely around the epicentral region and were recorded by a network of permanent stations monitored by the Geological Survey of Namibia (GSN) as well as other regional stations.
Seismic swarms have been observed to occur differently in terms of the extent and intensity of occurrence (Llenos et al. 2009;Hill et al. 1975;Smith et al. 2004;Hainzl and Ogata, 2005). The occurrence of seismic swarms is often associated with response to fault creep, fluid or magmatic events (Hill et al. 1975;Smith et al. 2004;Hainzl and Ogata, 2005;McGuire et al. 2005;, and has been exploited for evaluating earthquake triggering models in magmatically dominated environments (e.g. Dieterich et al. 2000;Toda et al. 2002). Analysis of seismic swarms requires a dedicated local high-density network capable of picking up even the smallest events that are due to micro stress transfer. In many cases, a general extent of the causative source is known and the network configuration is such that it covers the extent of the source. However, in regions like Anker where faults are not well mapped and their activity is not known, it is difficult to associate seismic swarms with specific faults. In 2018, the GSN and the Council for Geoscience (CGS) installed a temporary network of seismic stations to monitor the seismic swarm to help model the seismotectonics of the area. This knowledge is useful in seismic hazard studies as well as for decisionmaking. In this study, we present accurately located seismic events of the Anker seismic swarm recorded by a temporary network of seismic stations. Focal mechanism solutions are determined for selected good quality events. The results were used in modelling the seismotectonics of the region. In comparison to events of similar nature in other regions, we discuss the results of this study including the fault rupture geometry in the stable continental region (SCR) context, bearing in mind the possible seismic hazard implications.

Geology of the study region
Large areas of the Namibian geology exposed onshore are associated with the Late Proterozoic Pan-African orogenic cycle (Miller, 1983;Porada, 1989;Gray et al. 2006). The origin of orogenic systems of northern Namibia is connected to the the break-up of Gondwana during the Late Proterozoic and Early Phanerozoic (Prave, 1996;Trompette, 1997). It forms part of the much larger Pan-African orogenic system that surrounds and transects the African continent (Porada, 1989). The conclusions of many studies on geodynamic models have been proposed for the intracontinental branch of the Pan-African Damara Orogen which assumes subduction processes of either continental or ocean floor. The study area is situated within the triple junction and Neoproterozoic rocks of the Pan-African Kaoko and Damara mobile belts as well as the Congo Craton ( Fig. 1) mainly dominates its geology. Due to the collisional suture between the Congo and the Kalahari cratons that is associated with rifting which the intracontinental Damara Belt represents, makes the geology as well as the faulting system complex and hence tectonic stress changes might also be taking place although at a very slow rate typical of intracontinental regions. The mobile belts evolved through successive phases of intracontinental rifting, spreading and subduction as well as deformation that is characterised by open folding which decreases in intensity towards the north and east (Miller, 2008). The intracontinental Damara belt represents the collisional suture between the Congo and Kalahari cratons and chronicles a complete history from rifting and associated rift sedimentation and volcanism, via an oceanic stage and marine sedimentation through to convergence and final collision (Miller, 1983;Porada, 1989;Gray et al. 2006). Major  Gray et al., 2006) structures show predominantly two trends associated with the Damara (SW-NE trend) and the Kaoko mobile belts (NNW-SSE trend) many of which have been reactivated (Miller, 1983;Wanke, et al. 2000). The Kaoko mobile belt is a sinistral transpressional belt that is characterised by prominent sets of NNW-SSE striking, westerly dipping normal faults, with a down to the west sense of displacement (Durr and Dingeldey 1996;Maloof 2000;Passchier et al. 2002;Goscombe et al. 2003a,b), whereas the Damara belt is a divergent orogen that formed during high-angle convergence between the Congo and Kalahari cratons (Miller, 1983;Porada et al. 1983). The junction between the southern Kaoko belt and the Damara belt is the distinctive Ugab Zone with complex fold interference (Coward, 1983;Porada et al. 1983;Maloof, 2000;Passchier et al. 2002;Goscombe et al. 2004). The Congo Craton is composed of the Kamanjab Inlier, in the Kunene region, the Epupa complex and Grootfontein Inlier (Fig. 1), the three units are composed of metasedimentary sequence and volcanic intrusive granites and gneiss (Frets 1969;Porada, 1974;Becker 2006a, b, c;Miller, 2008).

The Kunene earthquake swarm
Seismicity in Namibia is generally low with a few isolated events within known seismogenic zones as delineated by Mangongolo and Hutchins (2008). In their findings, the Kamanjab region is not mentioned and this may be due to lack of seismic data at the time. However, of late, the region has shown some increase in seismicity that has caused a scare to the general population. In order to fully understand a seismic sequence, a network of seismic stations that can locate as many events as possible is needed close to the epicentral area, using the right configuration to constrain the location and focal mechanism solutions. This is because the accurate location of seismic events is controlled by the network geometry, available phases, arrival-time reading, suitable knowledge on the crustal structure and related velocity model (Pavlis, 1986;Gomberg et al. 1990;Waldhauser and Ellsworth, 2000). The Anker community reported an increase in frequency of seismic events that were felt since the beginning of 2018. In order to understand the source of these seismic swarms, the GSN and the CGS of South Africa installed a temporary network of 10 seismic stations to monitor the events for 3 months from 20 June to 24 September 2018 (Fig. 2). The seismometers used during the experiment were Trillium Compact 120s sensors.
More than 1600 microseismic events were recorded and located within the boundaries of the array roughly 50 km × 50 km (Fig. 3). The recorded data were analysed using the SEISAN software (Ottemöller et al. 2018) to determine the source parameters. The local magnitude scale developed by Saunders et al. (2013) for South Africa was used because it is the same scale that is currently being used by the GSN for magnitude determination in Namibia. The magnitude of the events ranges from − 1.1 to 3.6 M L with depth reaching up to 15 km. The map (Fig. 3) shows two clusters of seismicity with a general directional trend of NW-SE although there is a slight distinction between the upper and the lower clusters. The upper cluster is trending NNW-SSE and the lower cluster is trending in a more NW-SE and it is diffusely distributed (Fig. 4). The upper seismic cluster follows the direction of regional tectonic structures of the Kaoko mobile belt. The majority of lineaments from the magnetic data ( Fig. 4) indicate lineaments in the same direction, which may be an indication that these deep sitting faults maybe unstable hence resulting in an increase in microseismicity in the study area. Due to the limited number of seismic events recorded in Namibia, there is no locally developed velocity model for use in earthquake location. Rather, routine location at the GSN uses the Midzi et al. (2010) model developed for South Africa. However, for this study, because of the localised nature of the hypocentral distances for the analysed earthquakes, it was decided to recompute a new local 1-D velocity model in order to improve the earthquake locations. The computation of the local velocity model was conducted through a travel time inversion process that is implemented in the VELEST software (Kissling et al. 1994) using the data collected from the temporary network of stations. The program performs a simultaneous inversion for both velocity and hypocenter parameters by calculating a minimum 1-D velocity model which represents the least squares solution of the travel time residuals (Ellsworth, 1977;Kissling et al. 1994Kissling et al. , 1995. A selected high-quality seismic dataset was used for the velocity model inversion and to ensure accuracy in the results, only events located by at least eight stations were selected for processing. A modified version of the model by Midzi et al. (2010), which is currently being used in the analysis of data by the GSN, was used as the starting model. The model was modified by changing the thicknesses of the layers in the upper crust to about 2 km, whilst those in the lower crust were set at about 4 to 5 km. This was done to try to resolve optimum layer thicknesses, since the inversion is only done for the velocities and not for layer thickness. Smaller layer thickness values were selected for the upper crust because layer velocities vary more than in the lower crust. The starting model was iteratively refined to establish a reliable model geometry and velocity intervals in preparation for an updated a priori 1-D velocity model which was then used as the initial model for the final inversion process. In the final inversion process, travel time data were jointly inverted to obtain a 1-D velocity model (Fig. 5). This solution was then used for relocation of the analysed seismic events as well as for focal mechanism computation for selected earthquakes in the Anker region.

Relocation of the recorded seismic events
To better understand the seismotectonic model for the Anker region, the earthquakes recorded by the temporary network were relocated using the double difference method (Waldhauser and Ellsworth, 2000). This method minimises the residuals between observed and theoretical travel time differences for pairs of earthquakes at each station whilst linking together all observed event-station pairs. Absolute P-and S-wave travel time differences were used, although the location method can also incorporate the cross-correlation of differential travel time measurements. Only wellrecorded events with at least eight stations were relocated using the newly developed velocity model.
The spatial distribution of the original and relocated earthquakes in the region with depths of events is shown in Fig. 6. The main change in the distribution of relocated events is the disappearance of the two separated clusters of events. Rather, there is now one large cluster of events collocated with the original southern cluster. The clustering appears to have increased (Fig. 6), implying improved locations of events. Three distinct clusters are observed in the depth profile (Fig. 6b) all dipping in a westerly direction. The shallower part is to the eastern side and has depth between 0 and 20 km. The western deeper cluster has depth between 10 and 25 km.
After the relocation process, the earthquakes can be observed to follow a similar trend observed in Fig. 3 of NW-SE trend although they seem to cluster around the Anker region. There are no detailed structural mapping studies done for the study area to corroborate our observation; however, the regional tectonic setting within this region shows typical trends obtained in this study. The NW-SE trend is associated with the Kaoko mobile belt, a prominent geological structure in northern Namibia and known seismic zone (Mangongolo and Hutchins, 2008).

Focal mechanism computation
Focal mechanism solutions for some of the cluster earthquakes were determined using the FOCMEC program which implements the first motion polarity method (Snoke et al. 1984). Phases were picked using the SEI-SAN software (Ottemöller et al. 2018). The basic idea is that the polarity of the first P-wave arrival varies between seismic stations at different directions from an earthquake. The first motions define the four quadrants, two compressional and two dilatational, and the division between the quadrants occurs along the main fault plane and an auxiliary plane that is perpendicular to the main fault plane. Each event selected had at least eight clear polarities. The fault plane solutions vary in type of faulting but normal faulting appears to dominate ( Fig. 7 and Table 1). However, many of the solutions had a left-lateral strike-slip component to them. In a few cases, reverse faulting mechanisms were also observed. This suggests that the area is highly faulted with faults oriented in all directions as shown in Fig. 7.

Seismotectonic analysis
The seismotectonic context of northern Namibia is typically of a SCR characterised by low rate of crustal deformation and temporally and spatially diffuse distributed seismicity (Mangongolo and Hutchins, 2008). This makes it difficult to associate any recorded seismicity to faults. Though no known fault could be clearly linked to the Anker swarm, it is possible that the events are linked to the regional triple junction due to the Kaoko and Damara mobile belts. The junction is characterised by faults oriented in both NW-SE and NE-SW directions. Seismic swarms have been associated with either transient stress changes, fluid migration or magmatic activities Fig. 7 Focal mechanism solutions obtained in this study and associated with tectonic lineaments (yellow lines) in the study area. Blue beach balls represent normal faulting, red represent reverse faulting and black beach balls are for strike-slip faulting  (Hill et al. 1975;Dieterich et al. 2000;. However, in this area, neither of those activities is evident. Results from the focal mechanism analysis indicate that regional tectonic stresses are predominant although it was not possible to associate with localised structures and lineaments. The current seismic swarm may be either a result of stress transfer within a set of faults that are rupturing in response to heterogeneous stress fields or the presence of pre-existing crustal weaknesses within which deformation is localised and possibly triggered by loading, unloading or non-tectonic sources (e.g. Calais et al. 2016;Manzunzu et al. 2017). This can accommodate localised seismic swarms as seen to occur in the region. Although the Kaoko and Damara mobile belts have not shown significant movement since Precambrian, the occurrence of the seismic swarm showing localised movement of the faults. Another possible explanation for the occurrence of different faulting mechanisms at the same place and time (Fig. 7) is that events might be taking place on non-optimally oriented faults or are caused by block faulting. The Anker area lies at the confluence of two tectonic regimes, the Damara belt that trends NE-SW and the Kaoko belt trending NW-SE. The presence of steep shear zones in this region that are typically inherited from Pan-African orogenic belts in this part of Southern Africa (Wanke et al. 2000) may aid in transient stress changes. However, intraplate regions often lack surface rupture linked to an individual earthquake making it hard to identify active faults. Intraplate regions deform slowly and it is difficult to identify surface geological features related to the seismic foci with any sort of precision, and therefore, only a few events can be correlated to the surface geology, making them hard to study (Gangopadhyay and Talwani, 2003). Further studies on fault interaction and static stress transfer need to be conducted in order to understand the occurrence of this seismic swarm.

Conclusions
More than 1600 microseismic events were recorded using a temporary network of 10 stations deployed for a period of 3 months. The magnitude of these events ranges from − 1.1 to 3.6 M L which occurred down to a depth of 15 km. The locations of these events were used to try and associate the seismicity to the local tectonic structures. In this study, a new average crustal velocity model was computed for the Anker local area and used for relocating the earthquakes. The new velocity model is based on extensive good quality local earthquake data which resulted into significant improvements in swarm earthquake locations. Relocation of the earthquakes indicates a NW-SE trend in seismicity which is consistent with the orientation of the Kaoko mobile belt. However, the study region lies at the junction of two mobile belts and focal mechanism solutions suggest that rupture mechanisms can be associated with both mobile belts. A mixture of focal mechanism solutions was obtained from this study indicating the complexity of the localised faulting in the area. Seismicity might be linked to reactivation of pre-stressed faults within the region due to changes in stress conditions. The seismotectonics of the Anker region might be under the influence of regional tectonic structures. Future work need to be undertaken in order to find best solution of the source of these seismic swarms in the study area. for Geoscience for their immense contribution in data acquisition. The authors are also grateful for the assistance given by Mr. Ewereth Muvangua for his contribution on geological information and Ms. Kombada Mhopjeni in proofreading our manuscript. Funding The research was funded by the Government of the Republic of Namibia (Ministry of Mines and Energy) and the Government of the Republic of South Africa (Council of Geoscience) in a form of a joint research project.
Data availability Data can be shared with any interested parties upon request.
Code availability Not applicable.

Conflict of interest
The authors declare no competing interests.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http:// creat iveco mmons. org/ licen ses/ by/4. 0/.