Abstract

GIS and remote sensing data for allowing detection of structural features, such as faults, offer opportunities for improving mapping and identifying the areas that are likely to be locations of faulting areas. Landsat ETM-7 satellite data images were used and band-5 was found as the most suitable band for lineament delineation, based on the ability to identify geological features. Four contributing factors, namely, drainage patterns, faults (previously mapped), lineaments, and lithological contacts layers, were parameters used in this study to produce a fault potential prediction map using the overlay model techniques. The potential map (fault susceptibility map) classifies the study area into five potential zones, namely, very low, low, moderate, high, and very high potential. The areas covered by moderate to the highest potential zones were considered as fault segments (fault lines) in the area. The comparison of the potential map and the published fault map by using GIS matching techniques shows that 75 fault segments (fault lines) in the potential map were not properly identified in the study area. The correlation between fault segments and faults data collected from field work stations shows that there were 39 fault segments which may represent new faults in the area being identified. The presence of these faults is not known from the literature; this leads to updating and revising of existing geological map of the study area.

1. Introduction

Faults are weakness zones in the brittle part of the lithosphere, along which movement can take place in response to induced stresses. When faults undergo displacement, depending on geological and structural conditions, strain markers can be formed on the fault surface [1]. The presence of faults in any area is based on displacement of rock layers. But also, most of faults are represented by some geological features such as drainage patterns, lineaments (linear features), and lithological contacts between rock units within the rocks of the area. The presence of faults may be indicated by these geological features (factors).

The term lineament was first introduced by [2, 3] who recognized the existence of linear geomorphic features and interpreted them as surface expressions of zones of weakness or structural displacement of the earth’s crust. Lineaments are linear features on the Earth’s surface, usually related to the subsurface phenomena. Generally, lineaments are related to large fractures and faults where their orientation and number give an idea of fracture pattern of rocks [4]. In the recent years, the lineaments have been defined as natural crustal structures that may represent a zone of structural weakness [5].

The drainage system, which develops in an area, is strictly dependent on the slope, the nature, and attitude of bedrock and on the regional and local fracture pattern [6]. Most stream networks are adapted to regional slope and geological structures, picking out the main fractures in the underlying rocks [7].

The contact between two lithologies can also appear as a linear feature. This contact may appear as a change in drainage pattern across the structural features [8] or the two units may have different spectral properties [9].

Based on the definitions of lineaments, drainage patterns (including pattern the length, spatial distributions), and the lithological contacts between different rock types, the faults could be mapped in the study area. These geological features mostly resemble a fault lines in the area. The most important features in the area are the presence of drainage lines patterns and fractures. Lineaments, drainages, lithological contacts, and previous fault lines data are important data and used in this research for fault segments mapping using GIS technique.

With the advent of remote sensing and computer technology in the geosciences, geological investigation and interpretation have entered a new era. Remote sensing technology is very efficient for collecting data. Computer technology, such as computer-based geographic information system (GIS), supplies a different method for data storage, integration, analysis, and display. The combination of remote sensing and GIS provides an optimum system for various geological investigations such as fault mapping [10].

Several studies have been carried out to enhance geological knowledge and revise existing geological maps by using optical remotely sensed data and discovered new faults that were previously unmapped [1116]. Furthermore, authors in [17, 18] conclude by suggesting that satellite images and/or geographic information systems (GIS) are useful for tectonic mapping and identifying previously unmapped faults.

The benefit of integrating different data and techniques in structural features such as faults identification and analysis made possible using remote sensing and GIS techniques is the ability to identify faults based on their characteristics. Creation and classification of various types of geological features can be carried based on the purpose of the study. Once the identification, preparation, and processes are complete, geographic information systems (GIS) functionality, such as vector and raster spatial analysis and overlay, can be employed for structural mapping and analysis using powerful software programs [19].

Generally, Yemen is located in the southwestern part of Arabian Peninsula, and it is affected by an active rifting zone of the Gulf of Eden and the Red Sea opening. This active (rifting) zone can cause the reactivation of the old fault systems and may create new faulting lines in the area. The active faults (such as, rock mass movements or sliding on the fault planes and bedrock lithologies cracking) may cause severe damages in some places, especially in urban or settlements areas. For example, bedrock cracking and rock mass sliding happened in the 2003 and 2010 in some places in Taiz city caused severe damages, and many houses, roads, and other infrastructures were destroyed.

The previous works (data) such as geological and structural maps are very old data, showing a few of fault lines in the study area. For this purpose, mapping of faults may lead to updating these data. In this paper, the using of remote sensing (data processing) and GIS (modeling) techniques could be helpful tools for mapping of faults, where the different geological features (factors) are taken in the considerations. The results of the work could be considered from the government or the people before or during the different construction projects, to increase the safety, especially in the urban areas.

The purpose of this study was to test the integration between remote sensing and geographic information system (GIS) for detecting the fault segments (fault lines) over the study area and to investigate the ability of this method in giving real results compared to the previous data with respect to the ground truth (field work).

The study area located in the southwest of Yemen has an area of 176 km2. It includes Taiz city and extends between Rubayi and Al-Hawban areas (Figure 1).

According to GSY [20] classification, the rocks in the study area are divided into two groups as follows.(1)Tertiary (granite body intrusive).(2)Tertiary (volcanic rocks such as basalt, rhyolite, dacite, ignimbrite, and ashflow deposits).

The geological map of the study area (Figure 2) was digitized from the geological sheet map of Taiz with scale 1 : 250,000.

2. Materials and Methods

Various data were used for fault segment susceptibility mapping (potential fault zones mapping) in the study area. Various image processing and enhancement techniques were applied on different remote sensing data including the Landsat ETM-7 (Enhanced Thematic Mapper) satellite images and SRTM (The Shuttle Radar Topography Mission) to get maximum information by using different approaches of extraction methods. The other data used in this study were topographic maps, geological maps, and field data.

During conducting this study, many different software packages were used since there was no single software that would process all steps in the analyses. Therefore, the diverse software used for the analysis in the current study is PCI Geomatica (version 9.1), ERDAS (version 8.4), and Ilwis software (version 3.3).

2.1. Data Processing

According to this paper the linear features delineation was based on decision of the most appropriate data, such as drainage patterns, lineament, fault (previously mapped), and lithological contact layers for mapping of fault segments using GIS. Many of researchers such as [2126] have utilized lineaments, drainage patterns, faults, and lithology factors as important factors to measure and map the susceptibility of geological hazard. Therefore, the input attributes of the structural features of an area should be studied as a first consideration.

The previous data about the study area including geological, structural, and topographic maps are digitized. To prepare these maps, first, the maps are converted into a digital format by using the scanner. There are some data which were extracted from data in digital formats, such as Landsat ETM satellite images and SRTM.

Generally, four digital layers (maps) such as drainage patterns, fault (previously mapped), lineament, and lithological contact were prepared and converted into secondary data such as drainage buffer, fault buffer, lineament buffer, and lithological contacts buffer layers.

These layers were converted into slicing layers with five classes. Each class was determined as 30 m. This is because these layers were generated from different data (ETM, SRTM, and previous maps) with different spatial resolutions and scales. Due to this, these layers should be enhanced and resampled (rescaled) to be suitable in resolution for GIS model. The resolution was selected as 30 m for this process based on the ability of satellite image of 30 m spatial resolution to identify the geological features in the area. Then, slicing layers were converted into weight layers (thematic layers). These procedures are illustrated in Figure 3.

Arbitrary classifications are still common; however, the main classification approaches are ranking, natural breaks, equal interval classes, equal area classes, and mean value and standard deviation intervals [25, 2729]. Hence, the classification determines the spatial distribution of buffering zones or susceptibility (classifying susceptibility) based on equal area classes.

Classes (five classes of each weight layer) and their weight values are given in Table 1. Each class has a weight to express the contribution to the occurrence of fault segments. The weight value was assigned to be between 1 to 5 (Table 1). These weight layers (maps) were input into GIS model. The resultant map (susceptibility map) was classified into different zones, very low, low, moderate, high, and very high fault susceptibility zones or equal area classes using ranking approach.

2.1.1. Lineament Factor

Lineaments are related to large structural fractures [4] and it may represent zones of weakness [22]. Also, the presence of tectonic structures such as linear features could be important factor in geological hazard occurrences [30, 31].

Filtering is useful in image enhancement. An edge-enhancing filter can be used to highlight any changes of gradient within the image features, such as structural lines [17, 28, 3237]. They generated lineament maps and determined several significant structural features.

In the recent years, the lineaments have been defined as natural crustal structures that may represent a zone of structural weakness [38]. Landsat ETM-7 satellite data were used and the first step was to select the band that should be used for lineament extraction [39]. Visual inspection of the individual bands was carried out, based on the ability to identify features, and band-5 (1.55–1.75 μm) (SWIR) was selected as shown in Figure 4. And it was stretched linearly to output range from 0 to 255 (Figure 5), because it is the least affected band by the scattering, through the travel path of the long wavelength in the atmosphere. Therefore, it shows a good contrast and better display of geological features compared to other bands.

The second step was to select the filter type. For this purpose, different types of filters are tested such as 3 by 3 sobel kernel filter, 5 by 5 edge kernel filter, and 7 by 7 edge kernel filter. The 5 by 5 and 7 by 7 edge enhancement filters give thicker and less linear areas after threshold application; sobel results in thinner and more linear zones that give the best result compared with other two filters.

The sobel filter was the best which convolved over band-5 in four principal directions and the filtered image was added back to the original image. The sobel filters in four principal directions are given in Table 2.

Then, the four enhanced filtered images combined together to produce one threshold filtered image (Figure 6), and this image was converted to binary image by applying adaptive thresholding. This adaptive threshold used to create binary images is controlled by two parameters, and the suitable values of adaptive thresholding parameters are windows size of 3, and this value is equal to thinning of pixels and the mean offset of 3. The final binary image of the study area is shown in Figure 7.

The lineaments detected during the interpretation process were digitized directly on the image on the screen, and the final lineament map was recorded and stored in vector files (Figure 8).

There are important points which need to be mentioned and should be followed in relation to the above procedure for lineament mapping. First, care is needed in the interpretation of the result of such analysis to exclude man-made features; therefore, the final map was screened to remove lineaments related to cultural features such as roads, canals, and field plantation boundaries by comparing it with the topographic maps and geological map of the area. All procedures mentioned above have been followed during trace lineaments from binary image.

The lineament map was prepared and converted into lineament buffer layer and then converted into slicing layer with five classes. Finally, this layer was converted into weights layer (Figure 9). The weight value was assigned to be between 1 and 5 (Table 2). And, this layer was ready to be input into GIS modeling.

2.1.2. Drainage Patterns Factor

The drainage pattern is apparently being controlled by structure and lithology in the study area. The lithologic variation has given a rise to different drainage patterns. For example, radial and dendritic drainage patterns are developed over granitic rocks. Moreover, the most important feature in the area is the presence of drainage lines patterns and fault lines. It is clearly to see that there is a good relationship between these fault lines and drainage pattern system distribution especially with third and fourth river orders in the area.

Geological features are any alignment of features on satellite images such as the various types recognized including topographic, drainage, vegetative, and color alignments. Digital elevation models (DEMs) are very useful in aiding the classification of landforms for many purposes such as recognition of drainage patterns. Extraction of topographic feature information from DEMs has become increasingly popular in structural analysis [40]. Digital elevation models (DEMs) data were used to trace tectonic features and mapping geologically and topographically defined structures in many areas [41]. The Shuttle Radar Topography Mission (SRTM) data were also used to trace drainage pattern system and mapping structures features [42]. The extraction of negative system lines (fault lines) from DEM data mostly resembles drainage lines patterns [43].

SRTM radar images with 90 m resolution were used to obtain contour line data. This contour data was used to create the DEM with 30 m resolution of the study area. The interpolation and resampling processing techniques were used to convert the contour line data (polygon) into 30 m grid layer and the grid layer was converted into a DEM with 30 m spatial resolutions (Figure 10).

This DEMs data was used for automatic extraction of drainage patterns over the study area using drainage basin from elevation data function (DRAIN) of PCI Geomatica software version 9.1. The drainage basin from elevation data function (DRAIN) locates the drainage networks (river or stream) within the DEM. The output data channel will appear as lines representing the rivers and streams [44].

The result of the drainage pattern map using automatic extraction technique is shown in Figure 11. This technique helps us to extract the drainage patterns over the settlement area, and the topographic maps were used as reference data.

This map was converted into vector layer. The drainage pattern layer (vector layer) was converted into drainage buffer layer and then converted into slicing layer with five classes. This map was converted into weights layer (Figure 12). The weight value was assigned to be between 1 and 5 (Table 1). This layer was ready to be input into GIS modeling.

2.1.3. Lithological Contact Factor

Lithological contacts were found in the field, such as bedding planes and lithological contacts between different rock types. The lithological contact layer of the study area (Figure 13) was digitized from the geological sheet map of Taiz with scale 1 : 250,000. The lithological contact layer (vector layer) was converted into lithological buffer layer and then converted into slicing layer with five classes. And this map was converted into weights layer (Figure 14). The weight value was assigned to be between 1 and 5 (Table 1). This map was ready to be input into GIS modeling.

2.1.4. Fault Factor

Geological structures are patterns such as faults, folds, joints, and others found in the rocks resulting from natural deformation [45]. The fault map of the study area (Figure 15) was digitized from the geological and the hydrogeological sheet maps of Taiz area with scales 1 : 250,000 and 1 : 50,000. The first sheet was prepared by geological survey of Yemen (GSY) and the second was prepared by National Water Resources Authority, Dar El-Yemen Hydro consultant [46]. This layer (fault layer) was converted into fault buffer layer and then converted into slicing layer with five classes. And this map was converted into weights layer (Figure 16). The weight value was assigned to be between 1 and 5 (Table 1). This layer was ready to be input into GIS modeling.

2.2. GIS Model

The major goal of this study is to produce the fault susceptibility and the fault segment maps in order to get information about the possible locations of new faults existing in the area.

Generally, four digital weight layers such as drainage patterns, fault (previously mapped), lineament, and lithological contact layers were prepared from different data. The first step of the GIS modeling involves the input data (weight layers), and the second step includes selection of the suitable equation for this model to determine the susceptibility of faulting segments in the study area.

In this study, the overlay layers model (spatial multilayers analysis) was used preformed on these thematic maps with different weights by combining and aggregating these maps together with considering all different parameters. The fault susceptibility model produced with a modified linear equation of drastic model [47], using the following equation: where, DRNW is the drainage weight layer, FLTW the fault weight layer (previously mapped), LINW the weight layer, and LCTW the lithological contact weight layer.

Researchers such as [4850] classified the susceptibility into five potential classes, very low (VL), low (L), medium (M), high (H), and very high (VH).

The resultant map (susceptibility map) was classified into different zones, including very low, low, moderate, high, and very high fault susceptibility zones or equal area classes using ranking values.

2.3. Evaluation of the Susceptibility Map

The above-mentioned technique was used to produce the fault segment susceptibility map from the different data. The result obtained from this technique needs to be evaluated in some manner. The first evaluation was made with respect to the fault map (previous fault map) using GIS overlay technique to determine where the fault segments and true faults are matched over the study area.

During fieldwork, structural data (strike/dip of faults) were collected from 75 stations in the study area. These data were used for the fault segment map evaluation. These stations were located in the fault segments map to make whether there is any correlation between the fault segments and the faults collected from the field.

3. Results and Discussion

The final susceptibility map of faulting (potential fault zones) was produced as a result of the integration of all parameters of weight layers (drainage, fault (previously mapped), lineament, and lithological contact). This map was classified into five zones, namely, very low potential zone (VL), low potential zone (L), medium potential zone (M), high potential zone (H), and very high potential zone (VH), as shown in Figure 17.

The potential zones in the fault susceptibility map show that the very low potential zone has a total area of 22.87%, and the low potential zone has a total area of 51.89%. More than half of the study area are dominated by very low and low potential zones. And this area does not reflect faults (faults in the previous map), while the 22.40% of moderate, 2.73% of high, and 0.11% of very high potential zones areas almost cover all the faults.

The evaluation of the susceptibility map has been computed by using GIS overlay capability techniques with the fault map of the study area to determine where the potential zones and faults are matched. In order to evaluate our susceptibility map for faulting, a comparison with the known structural fault lines has been done. The result of this technique shows that all of the faults were matched with the moderate to high and very high potential zones (continuous areas from moderate to high and very high).

There were some of moderate to high and very high potential areas in the susceptibility map that could not match any fault lines in the output map of the fault map. These fault segments identified by this technique may be surface reflections of the new fault locations in the study area which were not identified by previous works. This means that the use of this applied model worked successfully to identify linear features (e.g., fault lines) in the area of interest. The results suggest an additional 75 of new fault segment lines, which were not recognized by the fault map over the study area (Figure 18).

Moreover, the fault segments obtained from this model were shorter in length compared with faults that were observed in the field. For this reason, this paper attempted to determine the accuracy of this technique for fault segments mapping and the relationship between these segments and ground truth (faults) in terms of the length and the spatial distribution in the study area.

Based on the pattern of fault segments including the length, spatial distributions, and visual interpretation with a topographic map, a drainage pattern map, and three dimensions of DEM and Landsat ETM, the true length of fault segments could be recognized from this analysis and could be implemented in order to assess the extraction accuracy of every true fault length in the area.

Validation of result was made with the geological truth (field work) in the study area. During fieldwork, many faults were identified in different field work outcrops. Strike/dip readings of faults have been collected during the field work from different rock outcrops in the area. Based on the field data, these linear features identified on Landsat imagery and DEM are interpreted to be faults (Figure 19).

The fieldwork stations were located in the fault segments map to make whether there is any correlation between the fault segments (faults lines identified by this technique) and the faults data collected from field work stations. There was good correlation between a few fault segments lines and faults, where 39 locations came in the same or nearby the location of fault segments (faults lines). Therefore, these 39 fault segments (faults lines) were identified as new faults (Figure 19). The structural reading (strike/dip) of faults from 93 fieldwork stations are given in Table 3.

Three fault systems were identified. These are the E-W, NW-SE, and NE-SW fault systems [46]. The area was subjected to intensive tectonic activity which resulted in a series of grabens, horsts, and tilted fault blocks. The tectonic activity deformed the area before the Tertiary volcanic activity, especially through the NE-SW faults which run parallel to the rift system of the Gulf of Aden, and it is considered the oldest fractures in the area. The NW-SE faults are run parallel to the opening of the Red Sea and occurred in the area during the Tertiary volcanic activity. The N-S and E-W fault system are the youngest major faults in the area. Smaller faults (as well as joints) also follow the same trends as the major faults. However, more intensive faulting is in the NW-SF direction [46].

Conjugate faults completely defined the principal stress direction. Strike and dip of the faults were plotted in equal-area projection of Schmidt net to obtain and calculate the directions of principal stresses which acted on those faults.

Totally 10 strike-dip readings of the conjugate faults from different rock units of the area were analyzed with the Stereonet program (Stero32 Version 0.9). The stereographic plot was prepared (Figure 20). According to the result of the conjugate faults analysis, compressional principal stress () is measured to be in the NE-SW direction.

This analysis indicates that the N-S and E-W orientated fault sets are sets of shear fractures which could be typical for strike slip faults. The NW-SE and NE-SW orientated sets are sets of tensional and extensional fractures, which could be typical for normal and reverse faults, respectively. Different field observations have also indicated that many damages were found or occurred along the NW-SE faults.

4. Conclusions

This study shows that remote sensing data and GIS technique are powerful tools in identifying geological structures such as fault segments in the urban area such as the Taiz city, Yemen. The use of different types of data such as drainage patterns, faults (previously mapped), lineament, and lithological contacts with spatial resolution of 30 m was combined through a process of integration and index overlay modeling technique for producing the susceptibility map (potential fault zones) of fault segments in the study area. GIS spatial overlay technique was used to determine the spatial relationships of all the criteria (factors) and subcriteria (classes) within layers (maps) to classify and map the potential area of fault zones of the study area. The fault segments map shows that there is good correlation and distribution between fault segments and faults (previous fault lines) in the study area. There are 39 fault segments which may represent new faults in the area being identified. The presence of these faults is not known from the literature. These faults could be active due to the rift system of the Gulf of Eden and opening Red Sea. Active or reactivated faults can cause enormous damage to settlement or urban (houses or construction project) areas by rock mass movement along the fault surfaces or bedrock cracks, damage which refers to loss of life or injury, or damage to roads. Furthermore, geological survey of Yemen is requested to update the maps and create new ones using new data and techniques and advise the people about the proper sitting of houses and restricting settlement activities in danger areas. GIS models can then yield more accurate estimates of the fault mapping using more data such as geophysical data.

Acknowledgment

The authors would like to thank the Geology Department, Faculty of Applied Science, Taiz University, Taiz, Yemen for the support to this research and staff for their valuable helps and comments.