Spatial Analysis of Linear Structures in the Exploration of Groundwater

The analysis of linear structures on major geological formations plays a crucial role in resource exploration in the Inner Niger Delta. Highlighting and mapping of the large lithological units were carried out using image fusion, spectral bands (RGB coding), Principal Component Analysis (PCA), and band ratio methods. The automatic extraction method of linear structures has permitted the obtaining of a structural map with 82,659 linear structures, distributed on different stratigraphic stages. The intensity study shows an accentuation in density over 12.52% of the total area, containing 22.02% of the linear structures. The density and nodes (intersections of fractures) formed by the linear structures on the different lithologies allowed to observe the behavior of the region’s aquifers in the exploration of subsoil resources. The central density, in relation to the hydrographic network of the lowlands, shows the conditioning of the flow and retention of groundwater in the region, and in-depth fluids. The node areas and high-density linear structures, have shown an ability to have rejections in deep (pores) that favor the formation of structural traps for oil resources.


Introduction
Knowledge of the geological nature of a given region, and its natural fracture networks, greatly affect the possibilities in analyzing and resolving many current problems, such as the exploration of subsoil resources and the rehabilitation of aquifers [1]. Remote sensing is used in the exploration of underground resources as a prerequisite for various studies, such as geological studies and geophysical studies (gravimetric, magnetometric, seismic, reflection, and refraction) [2,3]. It is, therefore, the research of geological structures in relief, because current technologies do not directly detect subsoil resources. By knowing these different aspects, the study of linear structures has become indispensable in resource exploration. Remote sensing, with its global view and geographical locations, is used to identify and orient linear structures. These linear structures are susceptible to reach in-depth rejections at a minimum of 2 km. The development of very high secondary porosities along the deep fault plane may favor the formation of structural traps for oil resources. Remote sensing, through its synoptic vision for studying vast geographical fields, is well adapted to the mapping of structures and lithologies in outcroppings [4,5]. Indeed, the processing of satellite data greatly contributes to the inventory and evaluation of natural resources and surface geological mapping. Remote sensing is also used for the visualization of linear structures and the lithologies of different zones. This facilitates the development and production of linear Structures maps and landfill maps, which are classified according to different lithologies.
The main objective of this work was to develop a map of linear structures (fractures, faults, etc.) and geological outcroppings (litho-structural map) from Landsat 7 images of the Inner Delta of Niger. The aim was to create an inventory of the linear structures and to trace the contours of the geological formations in order to help the relevant authorities to make decisions regarding underground resource exploration of the region.

Data and Study Area
Located in the center of Mali, the study area is approximately limited from west to east by the coordinates −6.20° and −1.50°, and, from south to north, by the coordinates 13.32° and 17.34° ( Figure  1). With an area of 131,886.607 km 2 , it is a depressed north-northwest dip zone [6], located in the transition zone between the Precambrian tabular of the Mandingo Plateau in the south, and the Precambrian folded of the Gourma in the north. And at northwest are accumulated the formations of continental terminal/Intercalary. The underlying primary formations include southwest/south-east oriented folds and west/east to northwest/south-east fractures that occur on the outcrops of Goundam Hills, the Bandiagara sandy plateaus, and outcrop sandstone to the south of Debo Lake. The geology of the shallow delta (flood zone) is mainly marked by continental or desert sedimentary deposits, comprising, from bottom to top, the following formations: The continental intercalary, continental terminal, and the quaternary [7]. The entire area is covered by a mosaic of 14 Landsat TM scenes, ranging from path 198 row 48 to path 196 row 51, of 2010. The images were taken from the database of the United States Geological Survey (USGS), and were acquired from their website. These raster images are meshed images forming a block of cells, and each of them contains information concerning the occupied space through the mean of the spectral signature of the objects in place. These images have a spatial resolution of 30 m and the appropriate bands for their spectral enhancement facilitate significant discrimination of the fracture networks. ENVI 4.5, Geomatica 2017, Arc Map 10.1, and RockWorks15 were used for preprocessing, processing and analyses.

Methodology
The adopted methodology in this work is based on two points. It begins with the application of remote sensing methods for lithological and structural mapping, and is followed by the characterization and analysis of linear structures and geologic formations. Figure 2 shows the elaboration and analysis of the litho-structural map. The entire area is covered by a mosaic of 14 Landsat TM scenes, ranging from path 198 row 48 to path 196 row 51, of 2010. The images were taken from the database of the United States Geological Survey (USGS), and were acquired from their website. These raster images are meshed images forming a block of cells, and each of them contains information concerning the occupied space through the mean of the spectral signature of the objects in place. These images have a spatial resolution of 30 m and the appropriate bands for their spectral enhancement facilitate significant discrimination of the fracture networks. ENVI 4.5, Geomatica 2017, Arc Map 10.1, and RockWorks15 were used for preprocessing, processing and analyses.

Methodology
The adopted methodology in this work is based on two points. It begins with the application of remote sensing methods for lithological and structural mapping, and is followed by the characterization and analysis of linear structures and geologic formations. Figure 2 shows the elaboration and analysis of the litho-structural map. Enhanced images make it possible to extract maximum lithological information. Several methods of investigation were applied in the discrimination of lithological units [8,9]; namely, spectral bands (enhancement), Principal Component Analysis (PCA), ratio bands (TM4/TM5, TM4/TM3 and TM4/TM7). In all cases, the application fusion using RGB coding was used first to merge the multispectral images of the visible (TM1, 2, 3) and then with the infrared images (TM4, 5, 7). The images resulting from these different processing methods were used for the realization of additive trichromatic synthesis. From the bands of the Principal Component (PC), a clear colored composition permits to discriminate the different lithological units, especially the PC band for automatic extraction of linear structures. The automatic extracting system of material is equipped with morphological processing operators, in which the directional criterion is directly integrated in order to automatically identify anisotropic structuring elements (rectilinear, curvilinear segments) [10,11]. Thus, the proposed extraction method adopted the following characteristics from Table 1 for the sought linear structures. Then, a split line at the vertices was applied for each of the vector files of the observed geological structures. Additionally, certain structures escaping automatic extraction were traced manually. The attribute information (coordinates of the two vertices of each linear structure) permitted the generation of statistics of the frequencies and lengths of the fractures. Density analysis of the linear structures of the lithology formations were the object of understanding of the subsoil state of groundwater and fluids. Enhanced images make it possible to extract maximum lithological information. Several methods of investigation were applied in the discrimination of lithological units [8,9]; namely, spectral bands (enhancement), Principal Component Analysis (PCA), ratio bands (TM4/TM5, TM4/TM3 and TM4/TM7). In all cases, the application fusion using RGB coding was used first to merge the multispectral images of the visible (TM1, 2, 3) and then with the infrared images (TM4, 5, 7). The images resulting from these different processing methods were used for the realization of additive trichromatic synthesis. From the bands of the Principal Component (PC), a clear colored composition permits to discriminate the different lithological units, especially the PC band for automatic extraction of linear structures. The automatic extracting system of material is equipped with morphological processing operators, in which the directional criterion is directly integrated in order to automatically identify anisotropic structuring elements (rectilinear, curvilinear segments) [10,11]. Thus, the proposed extraction method adopted the following characteristics from Table 1 for the sought linear structures. Then, a split line at the vertices was applied for each of the vector files of the observed geological structures. Additionally, certain structures escaping automatic extraction were traced manually. The attribute information (coordinates of the two vertices of each linear structure) permitted the generation of statistics of the frequencies and lengths of the fractures. Density analysis of the linear structures of the lithology formations were the object of understanding of the subsoil state of groundwater and fluids.

Lithological Mapping
Observation of the tones, textures, and arrangement of the various thematic objects contained in the image revealed the geological formations, as compared with pre-existing geological data of the area of interest. Figure 3 shows a color composition produced from the spectral bands, resulting from the application of image fusion techniques; namely, spectral enhancement (TM2, TM4, TM6). The three main areas of color represent the relevant regions for the analysis of this image. In the central part of the study area can be clearly observed as a domain of degraded texture. This texture extended from the south to the northeast and was more widened to the south. The quaternary of dominant sandy lithology showed variable dispersed granules, gravels, clay sands and sandy clays in the pools of ponds and lakes [12,13]. Unlike the sandy formations in areas adjacent to floodplains, clay-dominated components in lowlands, ponds, and lakes tended to limit water exchange between flood waters and the underlying formations of the continental terminal, and play a fundamental role in the hydrology of the delta [14,15]. To the west of the study area, there were some small areas with coarse textures. This is the continental terminal, with "sandy" as the dominant lithology. Often, local lenticular layers, with layers of interstratified clays, laterite horizons of hard cuirass, or gravel, are also interstratified. Marls and lacustrine clays are formed locally in the lithological series. The thickness of the continental terminal is quite variable. The eastern part of the image shows a field in texture that is less smooth and less coarse, and is light brown in color. With a rounded form, it was identified as the formation of calcareous and quartzite of lenticular forms, which, together, correspond to Precambrian A-folded.
Band ratios based on the notion of reflectance reduced the topography effects and enhanced the contrast reflections between mineral surfaces [16][17][18]. Ratios TM4/TM3, TM4/TM5, and TM4/TM7 were realized, and the colored compositions generated from the neo-bands allowed lithological identification ( Figure 4). A less-rough texture of blue-green was observed in the southwest part, with an extension towards the north. The identified lithology corresponded to limo-sandy alluvium (ancient alluvium). identification ( Figure 4). A less-rough texture of blue-green was observed in the southwest part, with an extension towards the north. The identified lithology corresponded to limo-sandy alluvium (ancient alluvium). Principal Component Analysis (PCA) consists of defining new channels that summarize the information contained in an image in multispectral space [19,20]. This method aims to maximize (statistically) the amount of information (or variance) of original data in a restricted number of components [11,21]. From the bands of Principal Components (PCs), a clear colored composition discriminated the different lithological units, and those same bands are specially used for the automatic extraction of linear structures ( Figure 5). The smooth textured area was mostly located in the north and was the pink and blue color. The generated lithology was a mixture of dune and sand. A scratched structure, dark blue, was located to the east. These longitudinal textures on the eastern coast were identified as metamorphic schists. These metamorphic schists are part of the Precambrian A-folded metamorphic. Two different textures were identified in the south-east, the first a greenish color, and the second in violet crossed by yellow. The first occupied more space to the south and was also spread across some areas to the north-east of the study area. The existing lithology is schistose sandstone. The second was observed to the south in certain small areas representing granules of sandstone as its lithology. Considering the arrangement of these two lithological formations, they belong to the same stratigraphic family of the Precambrian A-Tabular [22]. The first (schistose sandstone) is Precambrian A-Tabular of higher formations, and the second (sandstone of granules) is Precambrian A-Tabular of medium formations. Principal Component Analysis (PCA) consists of defining new channels that summarize the information contained in an image in multispectral space [19,20]. This method aims to maximize (statistically) the amount of information (or variance) of original data in a restricted number of components [11,21]. From the bands of Principal Components (PCs), a clear colored composition discriminated the different lithological units, and those same bands are specially used for the automatic extraction of linear structures ( Figure 5). The smooth textured area was mostly located in the north and was the pink and blue color. The generated lithology was a mixture of dune and sand. A scratched structure, dark blue, was located to the east. These longitudinal textures on the eastern coast were identified as metamorphic schists. These metamorphic schists are part of the Precambrian A-folded metamorphic. Two different textures were identified in the south-east, the first a greenish color, and the second in violet crossed by yellow. The first occupied more space to the south and was also spread across some areas to the north-east of the study area. The existing lithology is schistose sandstone. The second was observed to the south in certain small areas representing granules of sandstone as its lithology. Considering the arrangement of these two lithological formations, they belong to the same stratigraphic family of the Precambrian A-Tabular [22]. The first (schistose sandstone) is Precambrian A-Tabular of higher formations, and the second (sandstone of granules) is Precambrian A-Tabular of medium formations.

Structural Mapping
PCA images also allow the highlighting of linear structure networks, and, possibly, the geological accidents associated with it. Numerous regional linear structures have, thus, been automatically extracted from the geological formations of the PCA image. This leads to the formation of main corridors by the linear structure densities located from the center to the north and to the southeast. Figure 6 illustrates a linear structural map of the interior delta region. This mode of extraction, based on a few parameters, was the object of vectorization of an important network of linear structures, ranging from a few meters to several kilometers. This map resulted from the exploitation of raw and derived images, which led to the delimitation of geological formations [23]. This map will be used to generate the litho-structural map. Based on this knowledge, the characterization of geological formations and hydrological analyses were carried out.

Study of the Intensity of Linear Structures
The realization of iso-value maps of fracture density distribution using stratigraphic layers permits better observation of this relationship. A global analysis of this map showed that highdensity areas occupy less than a quarter of the studied surface. These high-density fractures occupy

Structural Mapping
PCA images also allow the highlighting of linear structure networks, and, possibly, the geological accidents associated with it. Numerous regional linear structures have, thus, been automatically extracted from the geological formations of the PCA image. This leads to the formation of main corridors by the linear structure densities located from the center to the north and to the southeast. Figure 6 illustrates a linear structural map of the interior delta region.

Structural Mapping
PCA images also allow the highlighting of linear structure networks, and, possibly, the geological accidents associated with it. Numerous regional linear structures have, thus, been automatically extracted from the geological formations of the PCA image. This leads to the formation of main corridors by the linear structure densities located from the center to the north and to the southeast. Figure 6 illustrates a linear structural map of the interior delta region. This mode of extraction, based on a few parameters, was the object of vectorization of an important network of linear structures, ranging from a few meters to several kilometers. This map resulted from the exploitation of raw and derived images, which led to the delimitation of geological formations [23]. This map will be used to generate the litho-structural map. Based on this knowledge, the characterization of geological formations and hydrological analyses were carried out.

Study of the Intensity of Linear Structures
The realization of iso-value maps of fracture density distribution using stratigraphic layers permits better observation of this relationship. A global analysis of this map showed that highdensity areas occupy less than a quarter of the studied surface. These high-density fractures occupy This mode of extraction, based on a few parameters, was the object of vectorization of an important network of linear structures, ranging from a few meters to several kilometers. This map resulted from the exploitation of raw and derived images, which led to the delimitation of geological formations [23]. This map will be used to generate the litho-structural map. Based on this knowledge, the characterization of geological formations and hydrological analyses were carried out.

Study of the Intensity of Linear Structures
The realization of iso-value maps of fracture density distribution using stratigraphic layers permits better observation of this relationship. A global analysis of this map showed that high-density areas occupy less than a quarter of the studied surface. These high-density fractures occupy 16,507.84 km 2 (12.52%) of the total area that covers the study area (Figure 7a). High-densities form two large blocks, located in the north and center areas, while still occupying the central hinge of the study area with extensions in the center/east and the east/south directions. A medium density zone in terms of the number of fractures clearly appeared along the regional axis Djene-Mopti-Bandjagara-Douentza. Low-density areas surrounded high-density areas and are located on South/east, east/center and center/north directions.
From a geological point of view, the high fracturing intensities corresponded mainly to the fields of recent alluvium deposits. An area with a high structural tendency, identified using remote sensing in the south of Tombouctou (in dunes and sandy areas) was also a high-density area of intersection of fractures (Figure 7b). The area of high densities corresponded to an underground corridor, which is expressed by the existence of low-lands, housing the Niger River and its tributaries, and by a large number of linear structures. The least-structured parts were identified as: On the one hand, precambrian A-folded metamorphic, the sand dunes and the sandy zones in the Gourma-rharous; and, on the other hand, in ancient alluvium, the continental terminal, and in one part of the recent alluvium deposits at Niono.  16,507.84 km 2 (12.52%) of the total area that covers the study area (Figure 7a). High-densities form two large blocks, located in the north and center areas, while still occupying the central hinge of the study area with extensions in the center/east and the east/south directions. A medium density zone in terms of the number of fractures clearly appeared along the regional axis Djene-Mopti-Bandjagara-Douentza. Low-density areas surrounded high-density areas and are located on South/east, east/center and center/north directions. From a geological point of view, the high fracturing intensities corresponded mainly to the fields of recent alluvium deposits. An area with a high structural tendency, identified using remote sensing in the south of Tombouctou (in dunes and sandy areas) was also a high-density area of intersection of fractures (Figure 7b). The area of high densities corresponded to an underground corridor, which is expressed by the existence of low-lands, housing the Niger River and its tributaries, and by a large number of linear structures. The least-structured parts were identified as: On the one hand, precambrian A-folded metamorphic, the sand dunes and the sandy zones in the Gourma-rharous; and, on the other hand, in ancient alluvium, the continental terminal, and in one part of the recent alluvium deposits at Niono.

Statistical Analysis of Linear Structures
The fracture map obtained after various processing stages had approximately 82,869 linear structures of varying sizes, at an interval of 0.03 km; 9.069 km, with a mean of 0.648 km, a median of 0.524 km and a standard deviation of 0.410 km. Beyond the mean length of 0.648 km, are listed the majority of the linear structures which represent 64.06%. 84.35% of linear structures are less than 1 km, 99.98% are less than 5 km, which means only 15.63% of the linear structures have a length between 1 km and 5 km. Therefore, a low density of structures was located at the edges of the study area. This proved that high-density areas were suitable for groundwater circulation, including the central-north corridor with its contribution of watercourses. The rejection of these fractures, at great depths, shows a great influence on the retention of fluids.
Linear structures were grouped according to their orientation into 12 classes by an angular crescent of 15 degrees. The frequency and density of the fractures, per orientation class, were calculated. The directional fracturing rosettes, expressed in number and in length, are shown in Figure 8. The distribution of the fracturing on the directional rosette was almost homogeneous. Frequencies oscillated between 3% and 8.25%, and no family of fractures exceeded 9% in frequency. However, some families stood out, with frequencies close to 9%. These fragments were N60-75, N45-60, and N270-285 with frequencies between 5.5% and 9%, and N0-15, N15-30, N285-300 with

Statistical Analysis of Linear Structures
The fracture map obtained after various processing stages had approximately 82,869 linear structures of varying sizes, at an interval of 0.03 km; 9.069 km, with a mean of 0.648 km, a median of 0.524 km and a standard deviation of 0.410 km. Beyond the mean length of 0.648 km, are listed the majority of the linear structures which represent 64.06%. 84.35% of linear structures are less than 1 km, 99.98% are less than 5 km, which means only 15.63% of the linear structures have a length between 1 km and 5 km. Therefore, a low density of structures was located at the edges of the study area. This proved that high-density areas were suitable for groundwater circulation, including the central-north corridor with its contribution of watercourses. The rejection of these fractures, at great depths, shows a great influence on the retention of fluids.
Linear structures were grouped according to their orientation into 12 classes by an angular crescent of 15 degrees. The frequency and density of the fractures, per orientation class, were calculated. The directional fracturing rosettes, expressed in number and in length, are shown in Figure 8. The distribution of the fracturing on the directional rosette was almost homogeneous. Frequencies oscillated between 3% and 8.25%, and no family of fractures exceeded 9% in frequency.
However, some families stood out, with frequencies close to 9%. These fragments were N60-75, N45-60, and N270-285 with frequencies between 5.5% and 9%, and N0-15, N15-30, N285-300 with frequencies between 4% and 5.5%. Apart from these families, the other frequencies were less than 4%. The analyses showed a strong similarity between length and frequency distributions of fractures. This result shows that the main directional classes of identified fractures appeared to be the longest. In other words, the most abundant fractures in frequency are also the longest. frequencies between 4% and 5.5%. Apart from these families, the other frequencies were less than 4%. The analyses showed a strong similarity between length and frequency distributions of fractures. This result shows that the main directional classes of identified fractures appeared to be the longest.
In other words, the most abundant fractures in frequency are also the longest.

Systems of Aquifers
The study area has five deep aquifer systems that correspond to main stratigraphic stages. They were classified into three categories, according to the type of groundwater deposit: • Crack-type aquifers are semi-continuous or discontinuous, depending on the density, extension, and degree of interconnection of cracking networks, which affect the host rock and its hydraulic relationship with surface layers in the overlay [24,25].

• Cracked Aquifers
Sedimentary, metamorphic, and Precambrian formations are characterized by very low intrinsic permeability. Water resources are almost exclusively associated with secondary permeability, originating, on the one hand, from cracking in the deep part of these formations, and, on the other hand, from alterations of their upper parts [26]. -

Aquifer of the Precambrian Tabular
With an area of 15,374.8 km 2 (Table 2), the Precambrian Tabular is the most important cracked aquifer in the region, and the most exploited, with 313 boreholes on a total of 970 boreholes of the

Systems of Aquifers
The study area has five deep aquifer systems that correspond to main stratigraphic stages. They were classified into three categories, according to the type of groundwater deposit: • Crack-type aquifers are semi-continuous or discontinuous, depending on the density, extension, and degree of interconnection of cracking networks, which affect the host rock and its hydraulic relationship with surface layers in the overlay [24,25]. Cracked aquifers are represented by crystalline and sedimentary formations of the Precambrian.
• Generalized aquifers are associated with low or unconsolidated formations with intergranular porosities, mainly consisting of deposits of continental origin, accumulated in large sedimentary basins from the Secondary to the quaternary.
• These deep aquifer systems are covered by superficial aquifers located in lateritic alteration formations on the surface of plateaus or in the alluvial and colluvial deposits of the plains and valley bottoms. Depending on the thickness of the upper layer, and the local geomorphological and rainfall conditions, these superficial aquifers can be semi-continuous and in hydraulic connection with deep or discontinuous aquifers.

Cracked Aquifers
Sedimentary, metamorphic, and Precambrian formations are characterized by very low intrinsic permeability. Water resources are almost exclusively associated with secondary permeability, originating, on the one hand, from cracking in the deep part of these formations, and, on the other hand, from alterations of their upper parts [26].

-Aquifer of the Precambrian Tabular
With an area of 15,374.8 km 2 (Table 2), the Precambrian Tabular is the most important cracked aquifer in the region, and the most exploited, with 313 boreholes on a total of 970 boreholes of the study area (Figure 9). This subdivided aquifer corresponded to the sandstone strip occupying the southeastern part of the study area in the Sudano-Sahelian climatic zone. With the exception of the Dogon Plateau, where sandstones are flush, the Precambrian is almost completely covered by lateritic formations. The Precambrian is made up of monoclinal griseo-schistous formations, which were characterized by an important development of secondary permeability, with an average main water flow of 48 m, average alteration thicknesses of 16.6 m, and an average static level of 48 m/ground.

-Aquifer of the Precambrian folded Metamorphic
The Precambrian folded covers a total area of 11,558.71 km 2 ( Table 2). The aquifer of the Precambrian folded metamorphic has a width of 4061.68 km 2 , and differs from the Precambrian Tabular due to of its hydrogeological characteristics, with Gourma-Rharous beings the most extensive unit. -

Aquifer of the Precambrian folded Metamorphic
The Precambrian folded covers a total area of 11,558.71 km 2 ( Table 2). The aquifer of the Precambrian folded metamorphic has a width of 4061.68 km 2 , and differs from the Precambrian Tabular due to of its hydrogeological characteristics, with Gourma-Rharous beings the most extensive unit. •

Generalized Aquifers
These aquifers are characterized by permeability of intergranular types and by a continuous tablecloth covering of the major part of the study area (Sahelian and desert region), making up a total area of more than 59,883.81 km 2 (45.41% of the study area). They are essentially constituted by detrital formations of continental origin, which accumulated since the secondary period in the Niger Basin. -

Aquifer of the Continental Terminal/Quaternary
This aquifer covers the vast alluvial plains that extend on either side of the Niger River. With a total area of around 54,746.69 km 2 ( Table 2), it is mainly located in the climatic and sub-desert zones of the sahalian area. The hydrogeological characteristics of this aquifer are related to the presence of permanent surface waters and the extension of the flood zones covered by the Niger River [25]. The continental terminal of the aquifer is composed of a mix of clays, sandy clays, and fine and coarse sands, often in lenticular and locally-gresified layers. Quartz and clay gravel or concretionary lateritic horizons are frequently intercalated in different parts of continental formations. -

Aquifer of the Continental Intercalary/Continental Terminal
This aquifer occupies certain portions in the western part of the study area, located in Niono, Niafunke, Tenekoun and Gorma-rharous, covering a total area of nearly 5137.12 km 2 , i.e., 3.90%. These ancient continental series are covered by quaternary dunes and, locally, by Holocene lacustrine deposits exposed in inter-dune valleys. Due to the depressed piezometry, which characterizes this aquifer system, the continental terminal and quaternary formations are generally de-watered and water resources are, therefore, mainly contained in the aquifer of the intercalary Continental [27]. These aquifer formations consist of local coarse sands and sandstone, with intercalations of more-orless sandy clays. The water levels are deep with an average measurement of 45 m and a maximum depth of 60 m.

Generalized Aquifers
These aquifers are characterized by permeability of intergranular types and by a continuous tablecloth covering of the major part of the study area (Sahelian and desert region), making up a total area of more than 59,883.81 km 2 (45.41% of the study area). They are essentially constituted by detrital formations of continental origin, which accumulated since the secondary period in the Niger Basin.

-Aquifer of the Continental Terminal/Quaternary
This aquifer covers the vast alluvial plains that extend on either side of the Niger River. With a total area of around 54,746.69 km 2 ( Table 2), it is mainly located in the climatic and sub-desert zones of the sahalian area. The hydrogeological characteristics of this aquifer are related to the presence of permanent surface waters and the extension of the flood zones covered by the Niger River [25]. The continental terminal of the aquifer is composed of a mix of clays, sandy clays, and fine and coarse sands, often in lenticular and locally-gresified layers. Quartz and clay gravel or concretionary lateritic horizons are frequently intercalated in different parts of continental formations.
-Aquifer of the Continental Intercalary/Continental Terminal This aquifer occupies certain portions in the western part of the study area, located in Niono, Niafunke, Tenekoun and Gorma-rharous, covering a total area of nearly 5137.12 km 2 , i.e., 3.90%.
These ancient continental series are covered by quaternary dunes and, locally, by Holocene lacustrine deposits exposed in inter-dune valleys. Due to the depressed piezometry, which characterizes this aquifer system, the continental terminal and quaternary formations are generally de-watered and water resources are, therefore, mainly contained in the aquifer of the intercalary Continental [27]. These aquifer formations consist of local coarse sands and sandstone, with intercalations of more-or-less sandy clays. The water levels are deep with an average measurement of 45 m and a maximum depth of 60 m.

Superficial Aquifers
Superficial Aquifers are recent formations, generally unconsolidated, of varying ages and origins. They play an important hydrogeological role as they constitute the first reservoir intercepting infiltration water generated by precipitation and runoff. The average thickness of the superficial aquifer is about 10 m. This aquifer is concentrated on the layers of superficial alteration (eluvial and colluvial) of other formations, and can be more-or-less thick and permeable, according to the nature of the original rock and fracture intensities. The presence of superficial alteration layers may be important in the vicinity of rivers because these rivers often settle in zones with a substratum, but also because of better feeding conditions.

Discussion
The application of processing and preprocessing techniques resulted in a radiometric enhancement of images, making them more expressive and finer for mapping structural and lithological factors. Based on the stratigraphic boundaries of the existing geological map, the images were enhanced to delineate the different geological formations. The colored composition, produced by the fusion of the images of the different parts of the electromagnetic spectrum, improved the differences in the spectral properties of various lithological formations, which facilitate their mapping. The analysis and interpretation of all these derived images resulted in the updating of the map of certain geological formations of the interior delta region of Niger. This update affected all geological formations, including central limits (gravel and laterites), north/east and south/west. Some formations are imposed by a high density of linear structures, extended in the south/east, east/center and center/north directions. Similarly, some northwest boundaries are defined by the stratigraphic formation of the terminal/quaternary continental. Few linear structures were extracted at this regional level, from Niono to Goundam. The updating also concerned the boundaries of the Precambrian folded in the East.
The processing of satellite images through combination of ratios of TM4/TM3, TM4/TM5, and TM4/TM7, highlights the biggest part of the ancient alluvium formation in Niono, in the southwest of the study area. This is the case for Precambrian folded formations, in term of structural tendency (weak structural tendency), located in these formations at the edge of the study area.
Processing of the images using the PCA technique followed by color composition permitted to complete the update of the geological map of the inner delta region of Niger. The use of 6 bands of Landsat TM from the total of 7 bands significantly contributed to improvement of the quality of geological discrimination. The improved contrast of the histogram had a positive impact on the images. In fact, the fusion of high-resolution spectral bands generated rich images, of which analysis and interpretation facilitated the identification of geological formations. Despite the richness of the derived images, the differences between certain geological surfaces are sometimes difficult to see. Indeed, in each of these surfaces, there appeared zones in which the variation in color tones was very poor. The contrast remained very low despite image enhancement. Precambrian A-Tabular-Medium formations and Precambrian A-Tabular-Higher formations encountered this difficulty. This is also the case in some formations of cracked aquifers and generalized aquifers, respectively, the metamorphic schist and dunes and sandy areas. These obstacles in identification were overcome using remote sensing and GIS techniques. In the north of the study area has been identified a strong density zone of nodes (intersection of the linear structures), which allows strongly to orientate studies of fluids exploration (oil, gas) in the region. The 672 node areas (8332.19 km 2 ) in the dunes and sandy area units were also assumed to have deep rejects (pores), which condition the flow and retention capacities of fluids (liquid or gas). The litho-structural ( Figure 10) mapping, resulting from remote sensing, fixes the lithological boundaries that remain globally similar to those defined in existing maps.

Conclusions
The results of this work confirmed interest in the use of numerical imaging in lithological and structural mapping. Indeed, the investigation of lithology in the Delta Interior of Niger using a combination of different methods of digital processing, namely spectral enhancement, Principal Component Analysis, and ratio bands, led to discrimination of different stratigraphic units of the interior delta. PCA neo-bands derived from the elimination of raw image noise (of 95%) also highlighted the linear structures for better automatic extraction. This extraction resulted in the elaboration of a detailed map of the linear structures of the inner delta of Niger after validation of approximately 82,869 fractures. Thus, analysis of the linear structures through their densities and their nodes, in the center and to the north, shows conditioning of mobilization of underground waters, and, more deeply, the flow and retention of fluids. These different processing methods led to the litho-structural mapping of the region, which proposes an update of the old geological map and the creation of a new structural map of the region.

Conclusions
The results of this work confirmed interest in the use of numerical imaging in lithological and structural mapping. Indeed, the investigation of lithology in the Delta Interior of Niger using a combination of different methods of digital processing, namely spectral enhancement, Principal Component Analysis, and ratio bands, led to discrimination of different stratigraphic units of the interior delta. PCA neo-bands derived from the elimination of raw image noise (of 95%) also highlighted the linear structures for better automatic extraction. This extraction resulted in the elaboration of a detailed map of the linear structures of the inner delta of Niger after validation of approximately 82,869 fractures. Thus, analysis of the linear structures through their densities and their nodes, in the center and to the north, shows conditioning of mobilization of underground waters, and, more deeply, the flow and retention of fluids. These different processing methods led to the litho-structural mapping of the region, which proposes an update of the old geological map and the creation of a new structural map of the region.