The seismological station of Abéché located in eastern Chad operates in a small network of short-period sensors and one three-component sensor. The geographical position of the station in the heart of Africa allows it to record most seismic events occurring in the world. Comparative studies to assess the station’s recording capacity were carried out from 1999 to 2003. These studies consisted of comparing the data from the Chad station with those from the REB bulletin (Review Event Bullen). The result for earthquakes of magnitude greater than or equal to 4 showed that approximately 80% of the teleseisms listed on the REB bullen are seen by the station in Chad, which reflects its high signal reception capacity. At the local level, the signals recorded are few in number and of very low magnitudes (on the order of 2 to 3), mostly from quarry shootings. The earthquake of May 2025 that occurred in the Borkou region (Figure 1) is the first of its kind to be felt by the population. The seismic calm that prevailed in the region after the earthquake shows that this is an isolated event. Although not destructive, the unexpected nature of the earthquake requires a thorough study to determine its real causes.

Figure 1. Map of Chad shown the Borkou Region.

The extract from the geological map of Chad ([14]) shows that the Chadian territory has mountainous areas distributed across the regions of Borkou, Ouaddaï, Wadi Fira, Guéra, and Mayo Kebi. These areas correspond to outcrops of Precambrian rocks. In the far north of the country lies the largest mountain range, the Tibesti, composed of a Precambrian basement covered by primary sandstone and then by tabular intercalated continental sandstone, all fractured and covered by Tertiary and Quaternary volcanic formations. The landscape features volcanoes such as Emi Koussi, which reaches 3415 m in altitude (the highest point in Chad), Pic Toussidé at 3315 m, and high plateaus called “Tarso” with altitudes between 2000 and 3000 m. The Emi Koussi volcano, whose last eruption is estimated to have occurred between 2.4 and 1.33 million years ago, remains an extinct volcano. To the northeast, one encounters Erdis, which are sandstone plateaus with ruin-like features, ranging in altitude from 500 to 800 meters. From a geophysical point of view, the studies in Chad started in 1954 with a team from ORSTOM. The team highlighted some large gravimetric routes ([6]), and short online electrical soundings. The electrical and seismic refraction prospecting methods practiced by the General Geophysical Company (C.G.G), have been oriented towards the search for water. The gravimetry work carried out subsequently highlighted:

the large trench that crosses all of southern Chad (Doba, Baké-Birao, Salamat) and which is located on the extension of the northern coast of the Gulf of Guinea.the large transversal Sao-Tome, Cameroon known as the “line of Cameroon” which would extend to Tibesti by a detour to the Erdis basin (Ounianga Kebir) ([5]; [3]).

In seismology, research in Chad began with the exploitation of data from Abeche station since 2000. Since this period, no major earthquakes have been recorded in Chad. One can nevertheless note the existence of small signals of very low magnitudes, mostly from quarry shootings but also from natural earthquakes sometimes felt like the one in 2015 in the Waddifira region and 2021 in the Hadjerlamis region.

Many studies have now been carried out in order to understand and interpret all the phases present on regional seismograms. Classically, we record the Pn, Pg, Sn, Sg and Lg phases (Figure 2). These phases are diffracted and refracted by the different heterogeneities they encounter along their propagation paths and thus form the secondary phases or the coda of the seismogram.

Figure 2. Sismogram of an earthquake in Sud Soudan located 300 km from station.

While classical wave trains are often well reproduced, the secondary phases that owe their existence to the heterogeneity of the propagation medium are generally more difficult to model ([1]). This difficulty comes on the one hand from the limited knowledge that we have of the propagation medium and on the other hand, from the misunderstanding that we have of the interactions between the wave field and the different heterogeneities of the medium. Indeed, although we now know relatively precisely the evolution of speed as a function of depth, lateral variations in speed are more difficult to assess. To this is added the difficulty of implementing digital techniques able to taking into account the three-dimensional propagation effects of the wave field ([11]).

A plane wave propagating through a seismological antenna is characterized by an oscillation frequency f 0 and a wave number vector:

where So designates the amplitude of the signal, k → 0 is the wave vector that indicates the direction of wave propagation, f 0 is the signal frequency, and r → ( x , y , z ) is the vector position.

The parameters of this wave field will be fully defined when we have known the frequency f 0 and the wave number k₀ that govern its propagation. One will then be able to estimate the direction of the wave which represents its azimuth θ 0 measured in degrees and its apparent propagation speed v a .

The real velocity v of the wave is given by:

Figure 3. Average background noise of Chad station.

The Abeche station consists of a mini network of short-period sensors and a three-component sensor. The geometry of the mini-network consists of a succession of nested triangles suitable for locating seismic events at local, regional and telesismic distances. The short-period seismometers used are of the ZM500 type, sensitive to the vertical component of the seismic wave. The stations are located on well-rooted granite blocks, which ensures good coupling with the ground, good protection against atmospheric disturbances. In addition, their distance from the urban area and the sea allows to reduce as much as possible the background noise of industrial origin and the permanent effect of the marine swell that spreads over very long distances. Figure 3 is a comparison of the average background noise of the stations with the NLNM (New Low Noise Model) background noise model of Peterson ([12]). This figure shows a good correlation with the model thus reflecting the quality of the signals recorded at the station. We see that between 0 and 1 Hz, all the stations have the same noise spectrum. This low-frequency noise is due to the combined action of the oceans and the atmosphere (oceanic swell). Beyond 1 Hertz, the local noise specific to each station is due to human activity and to the climatic environment.

Figure 4 represents the signal trace in the three Components (ABC) and in the one component station (TIRE2) and Figure 5 represents the time-frequency distribution of seismic energies.

Figure 4. The signal recorded at ABC and TIRE stations: Z component.

Figure 5. The time-frequency representation of the signal at T1R2 station.

We then calculated the polarization of the Pn, Pg and Sg waves recorded at station ABC.

The Waves polarization is estimated statistically from the data recorded by each of the North-South, East-West and vertical Z components ([9]). For this, the covariance matrix of the data is constructed:

where x_ij is the ith component of the ground motion, N is the number of points used on each component. T is the transpose of the matrix X.

The x_ij covariance matrix is a real and symmetric 3 × 3 matrix whose main diagonal elements represent the self-covariance of the components and the other elements represent the cross-covariance; it is written:

The polarization ellipsoid is estimated by finding the main axes of the ellipsoid by calculating the eigenvalues and eigenvectors of S using the equation

where I is the 3 × 3 identity matrix, and represent the eigenvalues and eigenvectors, respectively.

Then we determine:

In the horizontal plane, the azimuth of the wave by the north-south and east-west components:

The u_ij represent the director cosine of eigenvectors. The u i j sign function is introduced to resolve the ambiguity of 180˚ by taking the positive vertical component.

In the vertical plane, the angle of incidence between the vertical and the horizontal plane:

Figure 6 represents the signals of Z, N and E components of the station ABC filtered in [0.5 - 4.0] Hz frequency band. On this figure, the red brackets represent the time range used for polarization computation of each phase. On the bottom, we have the waveform polarization in horizontal plane giving the back-azimuth for Pg wave (Figure 7), Pn wave (Figure 8) and Sg wave (Figure 9). The two P phases are polarized in longitudinal plane, and Sg wave is polarized in transverse plane. Yellow marker indicates theoretical back-azimuth and red arrow represents measured back-azimuth.

Figure 6. Signals of Z, N and E components of the station ABC.

Figure 7. Polarization in horizontal plane giving the Azimut of Pg phase.

Figure 8. Polarization in horizontal plane giving the Azimut of Pn phase.

Figure 9. Polarization in vertical plane giving the Incidence angle of Sg phase.

Locationsettings

The localization is done in two steps: The first step consists to plot the arrival times at different stations. The second step consists of adjusting the parameters by hand, adjusting arrival times in order to obtain the best possible RMS. About the depth, this earthquake takes place in the upper part of the continental crust and therefore its depth is estimated between (2 - 3) km and (10 - 15) km. The localization is done by the ONYX software provided by the (LDG) Laboratory of Detection and Geophysics of Bruyères-Le-Châtel in Paris. The settings used in the localization software are listed in Table 1 below.

Table 1. Calculation parameters.

Figure 10. Arrival times at different stations.

Figure 10 shows the stations that contributed to the localization.

Only two stations in the Abéché array (TIR2 and ABC) have contributed to this location. The other stations are in the LDG network with which we collaborate. The result of the localization is summarized in Table 2 below:

Table 2. Location parameters.

The Lg waves were first identified by [13]. They are interpreted on the one hand as resulting from the superposition of upper surface wave modes (Rayleigh and Love waves) propagating inside the continental crust and on the other hand being the envelope of S-waves reflected on the base of the crust and whose all energy has trapped within [10], [4]. Lg waves propagate mainly on continental paths. Their velocity is close to the average velocity of S-wave in the upper part of the continental crust. It’s worth about 3.6 km/s. The earthquake studied have generated Lg waves recorded by two stations located in the north part of the Chad station. These are the IDL stations located beyond the Mediterranean and the EIL station located in Israel. Lg waves are very well observed on the EIL path while they are blocked on the IDL path that crossed the Mediterranean, Figure 11 and Figure 12. The Lg/P ratio is about 6 on the IDL patch, whereas this ratio is estimated at 1/10 on the EIL patch.

Figure 11. Lg waves recording.

Figure 12. Path of Lg waves from the EIL and IDL stations.

In the literature, we have not received any information on the occurrence of a major destructive earthquake in Chad. However, in recent years more or less significant earthquakes have been recorded, thus showing seismic activity in the in Chad. The current earthquake in the Borkou region is located on a fracture line highlighted by satellite, topographic and seismic data. This line consists of the spectacular lineaments on the edge of the mountain ranges oriented SE-NW and passing through the south of the Tibesti massif. This large lineament is the succession of faults that extend from the Great Lakes to the Hoggar crossing the Tibesti. These faults, locally observed in Algeria and Kenya, correspond to all the faults that have affected the Pan-African chains. They can replay and cause earthquakes. However, the non-reproducible nature of the event shows that it is an isolated earthquake, the kind we encounter but not often in stable continental regions ([2]). The homogeneity of the continental crust made it possible to record Lg waves at station EIL located at a very long distance. On the other hand, these waves are blocked on the path that crossed the Mediterranean. This observation is in perfect agreement with the mode of propagation of Lg waves which are strongly attenuated by an oceanic crust. Indeed, the station in Chad has been operating since 2012 and from then until today, only three significant earthquakes have been recorded in Chad: on June 23, 2015 in northeastern Chad, on March 20, 2021 in the Hadjer Lamis, and on May 25, 2025 in Borkou.