Possible Influence of Brittle Tectonics on the Main Road Network Built in the Central African Environment

1. Introduction

Roads and corridors where construction interacts with both the environment and ecosystems, have been means of communication and exchange between nations, peoples and civilisations since antiquity. And to ensure their sustainable use, they are the subject of numerous studies (sociological, scientific, technical, anthropological and many others). However, the elements (bridges, pavements, tunnels, scuppers, embankments) of a road or corridor built sometimes at very high cost (whose average price per kilometre is between 336,944 $ to 4,930,885 $ in Central Africa [1,2,3,4]) are threatened by faulting and landslides.

These movements are one of the most destructive natural phenomena of infrastructure in the world and Africa is not spared. All over the world in general and in Central Africa in particular, the phenomenon is not to be neglected, as Figure 1 and Figure 2 show that the road network is intercepted by faults and fractures, causing significant damages not only to houses, but also to road infrastructure and human lives, and therefore to socio-economic development. Although not widely reported in Central Africa, these movements are very recurrent and produce scouring of culverts, collapse of structures (bridges, scuppers, etc.) [5], and landslides (Figure 1) [6,7,8], slope failure, as studied in the work of Biswas et al [9] and Mbouombouo et al [10] in which it is established that slopes become unstable for slope values greater than 6.17⁰ in sandy-clay formations and for slope values greater than 10⁰ in gneissic formations respectively. Figure 1, based on the literature review, shows that roads are generally subject to rockfalls, collapses (longitudinal and transverse) and shears (longitudinal and transverse), which very often precede collapses. However, geological information does not always indicate the presence of several faults in the study area (Figure 2), yet according to the work of Aretouyap et al, Gweth et al, and Nguemhe et al, [8,11,12] there would be an enormous number of them. In the field, scientists need to know about fracturing to better orient the road network but due to many hazards such as vegetation cover, sedimentation, these faults and fractures are often hidden. The structural geology of an area has a significant influence on landslide occurrence and one way of integrating structural information into landslide risk assessment is through lineament mapping [13,14,15,16,17,18]. Lineaments can be defined as the surface expression of deep geological structures, usually corresponding to structural features such as geological contacts, foliations, schistosities, fold hinges, faults and fractures [19], which are structural lineaments. They may also correspond to man-made features (roads, pipelines, railways etc.) which are anthropogenic lineaments. Many studies have shown the importance of structural lineaments in land movement [8,20,21]. In addition, it is important to incorporate lineament density as it has proven to be effective in several fields of study, including geology, hydrogeology, geotechnics etc., in that it can highlight the most fractured areas, and therefore highlight areas with a particular anomaly [22,23,24].

Lineament extraction from geospatial data generally involves two approaches namely, the manual approach and the automatic approach, based on visual interpretation and the use of computer algorithms respectively [25]. In the manual approach, the extraction of lineaments is mainly influenced by the user’s experience while the automatic approach mainly depends on the performance of the software used and the data provided in the image used [26]. Since the introduction of Landsat MSS in 1980, remote sensing has been widely used in lineament studies [25]. Indeed, remote sensing is a set of techniques used to acquire information about an area of the Earth’s surface, without any physical contact with it. It allows the observation of very large and inaccessible areas in a short period of time, and therefore has a good ability to detect lithology, landforms and lineaments on imagery [27]. Several types of remote sensing images such as Landsat, IRS, LISS, ASTER, aerial photos and SPOT have been used for land movement assessment [25].

The geological maps available [28,29] show that there are nevertheless points where the road and fault networks meet in Central Africa. However, these maps do not provide a sufficient understanding of all the accidents such as landslides, shearing and collapses that occur on the aforementioned road network. Furthermore, Figure 2 shows that the fault network intercepts the road network in only 16 sectors, but in view of the many accidents recorded on these roads, it is clear that this map cannot be used to justify these numerous accidents. Landsat satellite images are often marred by a lot of cloud, which does not always allow good observation of the target [30,31]. In this work, the idea is to see whether by combining Landsat 8 images (for the 28 scenes used with cloud percentages of 0%, 0.31%, 3.24%, 6.02% and 46%) and DEM SRTM, we can achieve better mapping of accident-prone areas on these roads, and plan a sustainable layout for future infrastructure and, above all, better prevention and monitoring of disasters on the current road network. So, this work examines whether Landsat 8 OLI (Operational Land Imager) with a high cloud level and SRTM (Shuttle Radar Topographic Mapper) data can be used in both equatorial and coastal Central Africa zones to produce relevant mapping, locating fracture and landslide zones, in order to contribute not only to better road layout at lower cost and in a relatively short time, but also to better and sustainable prevention of geological disasters that may occur on these road network.

2. Study area

2.1. Geological setting

The Precambrian basement complex of Central Africa which recorded the evolution of the crust from the Mesoarchean to the Neoproterozoic era [32,33] comprises two major units, namely the Pan-African Fold Belt and the Congo Craton (Ntem Complex).

In the study area, the Pan-African Fold Belt is subdivided into three domains, namely the Southern Domain, the Central Domain and the Northern Domain [34] (Figure 2).

(1)

The southern domain is a huge allochthonous nappe unit bounded to the south by the Congo Craton, to the west by the Kribi-Campo fault and extending eastwards into the Central African Republic in the Bolé and Gbaya series [35]. It comprises high-grade metamorphic and igneous rock units, low to medium grade shales of sedimentary origin and metamorphosed Neoproterozoic metasediments under high pressure metamorphism [36]; (2) the central domain, which is bounded by the Sanaga Fault in the south and the Tcholliré-Banyo Fault in the north, includes the Adamaoua Fault and other anastomosing faults, constituting the N70E central Cameroonian shear zone system [37].

This domain consists of Archean to Paleoproterozoic high-grade gneisses intruded by broad syntectonic Neoproterozoic plutonic rocks of calc-alkaline affinity with high potassium (K) content, with emplacement ages around 600 Ma [44,45,46]; (3) the northern domain located west of the Tchollire-Banyo fault and extending along the western border of Cameroon consists of: high-grade metamorphic complex, comprising rocks of volcano-sedimentary and sedimentary origin due to the presence of orthogneiss with predominant dioritic content; metasedimentary rocks from the Neoproterozoic schist belts; and predominantly calc-alkaline Pan-African granitoids deposited between 660 and 580 Ma [47]. The volcanic rocks found there have tholeiitic and alkaline affinities. This area has experienced at least two episodes related to the Pan-African plutonism, the oldest of which, dated between 640 and 620 Ma, corresponds to pre- to syntectonic calc-alkaline granitoids and the youngest (580 Ma) corresponds to late-tectonic granitoids [48,49]. Three main successive tectonic events have been identified in this area by Ngako et al [50]: Crustal thickening (630-620 Ma); left lateral distortion movements (613-585 Ma), marked by the northern area shear zones and the Sanaga shear zone; right lateral distortion movements (585-540 Ma), marked by the northern area dextral shear zones and the central area shear zone.

The Ntem Complex representing the northern part of the Congo Craton is subdivided into three units: the Ntem Unit, the Nyong Unit and the Ayna Unit.

The Ntem Unit consists mainly of younger intrusive complexes, banded series and greenstone belts. The intrusive complexes include the tonalitic, trondjemite and granodiorite suites constituting the TTG (Tonalite-Trondjemite-Granodiorite) suite [52], and the charnokitic-magmatic suite [53]. The banded series consists of granitic gneisses, charnokitic gneisses, enderbergites and leptynites [54]. The banded iron formation (BIF), metagraywackes, pyroxenites, sillimanite paragneiss and garnet amphibolites constitute the greenstone belts [55].

The Nyong unit is a high-grade gneiss unit, dominated by biotite-hornblende gneisses locally occurring as grey TTG, charnockites, garnet-amphibole-pyroxenites and banded iron formations (BIF). The magmatic rocks found here (granodiorites, augen metadiorites, syenites) are represented by a SW-NE trending group of small intrusions [56].

Formed by crystallophyllian rocks, greenstone belts and intrusive granitoids [28], the Ayna unit comprises N120-N140E trending foliations superimposed by antepan-African folds and blastomylomites, and overlain by the Pan-African nappe [29].

The Ntem complex has experienced two major periods of deformation. The first period of deformation is marked by successive diapiric emplacement of the TTG and Mesoarchean charnockites, characterised by isoclinal stretching and folding, vertical lineation and foliation [52]; and the second one includes the development of N-S to NE-SW trending sinister shear zones, partial melting of the greenstone belt and TTG suite rocks, followed by granite formation [57].

2.2. Roads

The study area extends geographically between latitudes 1°38’5 "N and 13°4’50.7 "N and longitudes 8°30’50.6 "E and 16°12’11.8 "E and has a surface area of 475442 km². There are 18 national roads in total, but taking into account the slip roads we can count 25, including N1, N1A, N2, N2A, N3, N4, N5, N6, N6A, N7, N8, N9, N10, N11, N12, N13, N13A, N14, N15, N15A, N16, N17, N17A, N17B and N18 (Figure 3a–c).

N1, N1A, N6 and N14 ensure trade between Cameroon and Nigeria; N1, N12, N13 and N13A connects Cameroon to Chad; N2 and N17B link Cameroon to Gabon; N2A and N7 join Cameroon to Equatorial Guinea; N9 connects Cameroon to Congo; N1 and N10 connect Cameroon to the Central African Republic. N3 is the main road linking the port area (from the Atlantic Ocean) to all the other countries in Central Africa. The other national roads are intermediate roads leading to N1A, N1, N2, N2A, N6, N7, N9, N10, N12, N13, N13A, N14 and N17B.

3. Material and methods

3.1. Data

To visualise the geological features on the Earth’s surface, remote sensing methods were used for their ability to cover large areas in a short time.

3.1.1. Landsat 8 OLI/TIRS

Although the method used is broad-spectrum, it was important to select the sharpest scenes from those acquired by the Landsat 8 satellite since its launch into orbit.

In short, 28 Landsat 8 scenes acquired have been geometrically corrected in World Geodetic System (WGS-84, Zone 32N and 33N) and converted in Universal Transverse Mercator (UTM) coordinate systems, were used for this work and are listed in Table 1. The use of these images is justified by their spectral characteristics and their spatial resolution (30m), which are suitable for structural analysis. Launched in 2013, Landsat 8 is an American Earth observation satellite. The eighth satellite in the Landsat programme and the seventh to reach orbit successfully, it comprises two sensors, including OLI and TIRS (Thermal Infrared Sensor). The OLI sensor is characterised by 9 bands, including 4 in the visible range (0.43 - 0.67µm), 1 in the near infrared range (0.85 - 0.88 µm), 2 in the mid-infrared range (1.57 - 2.29 µm), 1 band in the Cirrus range (0.50 - 0.68 µm) and 1 panchromatic band (0.50 - 0.68 µm). All bands of this sensor have a resolution of 30 m except the panchromatic band which has a resolution of 15 m. The TIRS sensor consists of two bands, all in the Thermal Infrared range (10.60 - 12.51 µm) and each with a resolution of 100 m. In this work only some of the OLI bands were used.

3.1.2. SRTM (Shuttle Radar Topography Mission)

SRTM data of 30m resolution were also used in this work. 63 SRTM scenes covering the entire study area were used. Information on the SRTM can be found in Farr et Kobrick [60].

Both datasets used in this work were downloaded from the website https://earthexplorer.usgs.gov [61].

3.1.3. Field data

Field campaigns were carried out during the dry season from 2019 to 2023 using compasses, GPS, distance meters and altimeters to identify fault zones and areas of possible landslides [10], particularly on the N3 national road. This road was chosen as the pilot zone because it is the main road linking the port area (from the Atlantic Ocean) to all the other countries in Central Africa and belongs to both the equatorial zone and the coastal zone, characterised by a high level of cloud cover. Longitudinal and transverse faults more than 5 km long were investigated, as they have a very large lengthwise extension and are therefore difficult to modify by human activity.

3.2. Methodology

3.2.1. Preprocessing

In this step, after mosaicing all these scenes, radiometric processing, haze reduction and histogram equalisation were applied in order to improve the contrast of the satellite image; noise reduction was also applied, with the aim of filtering out the noise from the image. In addition, Pansharpening, which merges high resolution panchromatic and low resolution multispectral images to create a single high resolution colour image, was applied to reduce the resolution of the satellite image from 30m to 15m.

3.2.2. Processing

• Principal component analysis (PCA)

PCA is an analytical method used to compress the information contained in the original bands into new bands called principal components, by eliminating data redundancy [58,59,60]. It has been used in radiometric studies by several authors [62,63,64,65].

In this study, it was applied to the combination of bands 7, 6, 2.

• Directional filtering

The filtering operation consists in modifying the value of a pixel according to that of its neighbours [66,67,68]. The aim of this method is to highlight structural discontinuities and facilitate the discrimination of lineaments in all possible directions and to bring out the maximum number of lineaments. The filters used in this work are Sobel directional filters applied in four directions: N-S, E-W, NE-SW, SE-NW with convolution matrices of size 7×7 [11,69]. Several studies have shown the effectiveness of combining different image processing techniques to improve image contrast, thus better highlighting geological discontinuities and faults [70,71]. In this respect, PCA is applied here in combination with the Sobel filtering method. This task was carried out using ERDAS Imagine software.

• Shaded relief

Shaded relief thematic maps are a visual representation of the terrain consisting of grey values stored in raster images, derived from the Digital Elevation Model (DEM). Shaded relief is commonly used in the extraction of lineaments from the SRTM image [72]. For this purpose, some parameters have to be applied to the input SRTM data, such as changing the virtual azimuth of the sun while maintaining its elevation marker in order to create a shaded relief [73]. This will result in four shaded relief images that will be output files with an azimuth of 0°, 45°, 90° and 135°.

This task was carried out with ArcGIS software using the spatial analyst tool hillshade.

The boundaries between shaded and unshaded areas could indicate the presence of lineaments [74].

• Manual and automatic lineament extraction

Automatic extraction was applied to the principal component of the filtered Landsat image and to the shaded relief of the SRTM image in the four preferred directions. This was done with the PCI Geomatica software via the LINE module algorithm using the optimal parameter values presented in Table 2. Similarly, manual extraction using ArcGIS software was applied to the filtered Landsat 8 principal component and SRTM shaded relief, in the four preferred directions NS, NE-SW, NW-SE and EW. The application of the different processing methods on the Landsat 8 and SRTM images led to the lineament map of the study area. The main steps that led to the extraction of the lineaments are summarised in Figure 4. Lineaments indicating either faults or fractures are considered as faults if they result in a variation in height of at least 30 m on the DEM, and fractures if they do not.

• Rose diagram

The conventional method is to produce the directional rosettes proportional to the cumulative length of the lineaments by 10⁰ orientation class [73,75,76]. This task was performed using Rockworks software.

As lineament extraction can generate many artefacts especially by automatic methods [77], any lineament corresponding to anthropogenic features was identified and removed, with the exception of linear forms similar to roads, given the purpose of the study. Lineaments digitised from the geological map, DEM and slope were also useful for validation.

• Slope

Derived from the DEM (Figure 5), the slope map (Figure 6) expresses the change in elevation, where high values correspond to steep slopes and low values to flat areas. High values of slope correspond to abrupt and uneven changes, and are often key indicators of the presence of lineaments [78,79]. A slope is declared a landslide zone if it is greater than 6⁰ and its height greater than 30 m, even though it is true that, for slopes less than 30 m high, we can still observe small falls of particles torn off not only from the rock, but also from the small sedimentary cover above.

4. Results and discussion

4.1. Results

4.1.1. Lineaments obtained

Both manual and automatic extraction methods were used to map the lineaments (Figure 7a). The statistical analysis of the lineaments obtained in this work shows 2321 (of which 1052 manual and 1269 automatic) and 3142 (of which 1874 manual and 1268 automatic) lineaments extracted from Landsat 8 and SRTM images respectively. In order to facilitate the comparison between the lineaments obtained and the faults of the geological map, in the following work, only lineaments of at least 5 km in length were taken into account. The lineaments obtained by the manual approach have lengths ranging from 5 to 65 km with an average length of 11 km for Landsat 8 images, while the lineaments extracted from SRTM images have lengths ranging from 5 to 356 km with an average length of 20 km. For the automatic extraction, the lineament lengths vary from 5 to 35 km with an average value of 6 km, and from 5 to 30 km with an average value of 6 km for Landsat 8 and SRTM images respectively. The final lineament map was generated by integrating all automatically and manually extracted lineaments for each dataset (Figure 7a); thus 5463 lineaments were identified. After removing repetitive lineaments, 5222 lineaments were retained, constituting the final lineament map; their lengths ranged from 5 to 356 km with directions N160 and N33 respectively, and a mean of 12 km and a standard deviation of 14 km (Figure 7b). Lineaments shorter than 12 km are the most numerous and constitute a percentage of 73%. The statistical analysis of the lineaments has been investigated by several authors in order to study the geometry of the lineament network and to identify the dominant directions in the area.

In this work, two types of directional rosettes were performed, the length-based rosette (Figure 7c) and the number-based rosette (Figure 7d) of the lineaments. The rose diagram obtained as a function of the length of the lineaments shows a predominance of the N20-N30, N100-N110, N30-N40 and N110-N120 direction classes. The N20-N30 direction would therefore be the main direction of deformation, and the N110-N120 direction, which is perpendicular to it, would be the direction of major stresses.

4.1.2. Altimetry

The relief map (Figure 5) obtained from the DEM is divided into three classes, namely low altitudes ranging from 0 to 490 m in the coastal regions and the north, medium altitudes ranging from 490 to 944 m in the centre and south, and high altitudes ranging from 944 to 4029 m.

4.1.3. Slope

The resulting slope map has also 3 classes, low slope class ranging from 0 to 6; medium slope class ranging from 6 to 17 and high slope class ranging from 17 to 76. The high slope areas are concentrated to the west of the study area, slightly to the north and to the southwest. Areas of medium slope are found further south, east and west, and areas of low slope to the north and east.

4.2. Discussion

4.2.1. validation of road map

The road maps proposed by [10] (Figure 3a) and [58] (Figure 3b) have shortcomings, which has led us to propose a new map (Figure 3c), showing more details than previous ones, while at the same time bringing out 100% of the road network illustrated by the previous two maps.

First validation of lineament map

4.2.2. Correlation of lineaments with the geological map

Figure 8 shows the lineaments superimposed on the geological faults. Among the 20 faults of the geological map, 16 were found, which is a percentage of 80%. By overlaying, for an initial validation, the network of lineaments obtained having more than 5 km with that of the known faults, it appears 4 faults F1, F9, F17 and F18 (Figure 8) which could not be identified by the present work. Field validation of faults F1, F9, F17 and F18 was not possible due to difficult access conditions, so a new field campaign along the N3 was carried out to validate them. Table 3 shows the lineaments obtained and the geological faults to which they correspond and overlap. Figure 7 shows that the lineaments are highly concentrated along the Cameroon Volcanic Line (CVL), the Foumban Shear Zone (FSZ), the Southern Cameroon Shear Zone (SCSZ), and in the southern part of the Central Cameroon Shear Zone (CCSZ), in the northern part of the Kribi-Campo Shear Zone (KCSZ) and in the eastern part of the Congo Craton. In the far north, the concentration of lineaments is low, which would certainly be due to the presence of sand and alluvium covering the ground. Overlaying the lineaments on the geological map shows that the fracture density in the Congo Craton is 0.13 while that in the mobile zone is 0.14. This shows that the Congo Craton, which was thought to be stable, appears to be highly fractured, especially in the western and central flanks. It is understandable why earthquakes [80] are felt in the areas contained in this unit.

4.2.3. Correlation of the lineaments with some previous geophysical work

In addition to the geological faults found by [47], this work has highlighted the fault announced by [81] at the level of Akono-Mbalmayo but which until now has remained little known probably because of the particularly difficult access conditions. Among the faults highlighted by [8], 4 of the 5 confirmed faults and 7 of the 8 assumed faults were found. Only 3 of the confirmed faults intercept the main road network. In addition, the lineaments show a strong concentration around Lakes Nyos and Manengouba, as observed by [11,82].

Likewise, the directional rosette obtained as a function of the greatest number of lineaments showed a predominance of the N40-N50, N30-N40, N20-N30 and N0-N10 direction classes (Figure 7d). The N40-N50 direction was found in the work of [8], relating to the occurrences of fault and landslide in the West Cameroon region, identifying the N45-N50 direction as the predominant orientation of lineaments.

Similarly, the study carried out by [83], on mapping of major tectonic lineaments across Cameroon also revealed that lineaments are mainly oriented in the N45 direction. The N0-N10 direction obtained in this work was also identified by [83]. The N20-N30 direction was highlighted in the work of [11], on fracture comparison models in CVL, where the N20 and N30 directions were obtained for Landsat 8 and Landsat 7 images respectively.

4.2.4. Correlation of the lineaments with the elevation map

The correlations between Figure 7a and Figure 5 show a relief demarcation observed between Kribi and Ebolowa corresponding to the Kribi-Campo Shear Zone (KCSZ) fault family found in Figure 2 and Figure 7a. Between Bafoussam and Bertoua, the combination of Figure 7a and Figure 5 also shows a fault family in the Foumban Shear Zone (FSZ) and that the cliffs of Dschang and Ngaoundéré lie in zones of high elevation and high fracture density. Furthermore, the correlation between Figure 7a and Figure 5 shows that, of the 373 lineaments obtained, 255 are fractures, 115 are normal faults and 3 are fractures in one of their portions and normal faults in their extension. The lineaments categorised as fractures in Figure 8 can lead to cracks along the road like those observed in Figure 1c, 1e and those categorised as normal faults can lead to subsidence and block collapse (Figure 1f).

Second validation of lineament map

Field investigations for N3 showed that all faults and fractures (Table 4) with a longitudinal extension greater than 5 km in both coastal and forest zones (for the 28 scenes used with cloud percentages of 0%, 0.31%, 3.24%, 6.02% and 46%) can be mapped using the method described in the flow chart (Figure 4).

Similarly, a fault found on the map from satellite data may correspond to several fault systems, as in zones FZ5N3, FZ6N3, FZ7N3 and FZ8N3 (Figure 9, Table 4). The 03 faults F11N3, F12N3 and F13N3 (Figure 9) described in the geological documents [11] and mapped in this study were confirmed by the field work, as were the 26 other new faults. In addition, the zones of fractures (20) and slopes (18) highlighted in the course of this work and not all observable in the existing geological documents were also validated by this study. However, the accuracy of the location of fault zones varied from 0 to 100 m. In addition, the approach used in this work made it possible to map all the areas around which landslides may occur. But unlike the location of fault zones, the deviation in the location of landslide zones varied from 0 to 500 m on average in relation to the road. The 4 faults F1, F9, F17 and F18 not found previously could be due to the fact that their difference in height is less than 30 m (which corresponds to the resolution of the DEM used), or else they could be contacts between different geological formations, covered by a sedimentary layer. All these results show that the use of Landsat 8 and SRTM data can be effective for the sustainable layout of roads and motorways in Central Africa and throughout the world, while making a relevant contribution to reduce the cost of their construction, which is still very high in Central Africa (between 336,944 $ to 4,930,885 $ [1,2,3,4]).

4.2.5. Correlation of lineaments with the slope map

By overlaying Figure 7a and Figure 6, we observe a strong concentration of lineaments that correlate with areas of steep slopes showing that roads can be subject to numerous rockfall phenomena. Thus, Figure 6 allows to understand and justify the landslides on the national roads N1, N4, N5, N6, at the level of some localities [8,84,85].

4.2.6. Correlation of the slope map with the geological map

Correlating Figure 2 and Figure 6, it can be observed that the road sections N1, N2, N3, N4, N6, N7, N10, N11, N12, N13, N15, N15A located on gneissic formations could be prone to slope failures if they have slopes greater than 10⁰ according to the work of [10] (Figure 10). In addition, the sections of national roads N1, N3, N5, N6, N7, N12, N16 located on sandstone and clay formations could experience slope failures if they are located on slopes greater than 6.17⁰ according to [9].

4.2.7. Correlation of the slope map with the road network

The Figure 6 shows that in several places on N2, N4, N5, N6, N8, N16, and N17 there may be a tendency to observe slope failures as was the case of [84]. These particle falls are particularly heavy on the N6 road linking Cameroon to Nigeria, in a NW-SE and NE-SW direction.

4.2.8. Correlation of lineaments with the road network

The superposition of Figure 3 and Figure 7a shows that 373 fault and fracture lineaments intercept the road network from north to south and from east to west (Figure 11a), showing that tectonic activity can have consequences on the road network, unlike the geological map, where only 14 faults intercept the road network (Figure 2), hence the validity of the study.

Their lengths vary from 5 to 356 km with an average of 28 km and a standard deviation of 38 km (Figure 11b), indicating a wide dispersion of lengths. The rose diagrams plotted against the frequency (Figure 11c) and length (Figure 11d) of the lineaments shows that the lineaments interfering with the road network are predominantly oriented in the N60-N70 and N50-N60 directions respectively. Among these lineaments, 18 run along the road network, 2 overlap and 2 intercept. The results obtained show that there are not only 14 faults (Figure 2) intercepting the road network, but that there are still others that have created and could create shears and collapses on the road network, as shown in Table 5.

Compared with the geological map (Figure 2), this new lineaments map (Figure 7) makes it easier to predict accidents linked to tectonic activity on the road network. The statistical analysis made it possible to obtain a correlation coefficient of 0.071, thus testifying to the existence of a correlation between the road network and the lineaments in the study area. It can thus be observed for the lineaments intercepting the road network, accidents such as the road shear, as well as a transverse subsidence of one of the two blocks on either side of the lineament (as in the Figure 1d, 1f). Lineaments overlayed on the road network are likely to cause cracks in the road and longitudinal subsidence of one of the blocks on either side of the lineament (Figure 1c, 1e). Lineaments along the road network can create accidents similar to those intercepting it. Table 5 shows that the road network in Central Africa may be subject to 113 zones of possible rockfall or landslide, 117 of possible transverse collapse zones, 2 of possible longitudinal collapse zones, 271 zones of possible transverse shear and 15 of possible longitudinal shear [8].

5. Conclusion

This work studies the impact of geodynamic activity on road corridors in Central Africa, with a view to monitor the present road network in a sustainable manner while planning a better route for future corridors.

To achieve this, 28 scenes with cloud percentages of 0%, 0.31%, 3.24%, 6.02% and 46% of Landsat 8 OLI, and SRTM, coupled with field data, were used. Various image processing techniques were then applied to highlight the linear structures. These included PCA, spatial filters and shaded relief, all of which proved effective in detecting and extracting lineaments. The result of the road map shows 29 main roads and 6 secondary roads. While the lineament map shows 5222 lineaments with increased fracturing densities along the Cameroon Volcanic Line, around the FSZ, CCSZ, SCSZ, Dschang and Ngaoundéré cliffs and some crater lakes; the longest following the N20-N30 orientation and the most numerous following the N40-N50. Among the linear structures identified, 376 interfere with the road network, of which 256 are fractures, 117 correspond to normal faults and 3 to fractures in one of their portions and normal faults in their extension. These results are useful for sustainable planning, management, monitoring and extension of roads worldwide, but especially in Central Africa where the road problem is still high due to the lack of large-scale information at the time of their design and construction.