Thank you for visiting nature.com. You are using a browser version with limited support for CSS. To obtain the best experience, we recommend you use a more up to date browser (or turn off compatibility mode in Internet Explorer). In the meantime, to ensure continued support, we are displaying the site without styles and JavaScript.
Advertisement
Scientific Reports volume 13, Article number: 21569 (2023)
Metrics details
The Bir El-Qash area, located in the Central Eastern Desert of Egypt, is characterized by a diverse range of igneous, metamorphic, and sedimentary rocks with ages spanning from the Late Proterozoic to Quaternary. Integration of remote sensing with aeromagnetic data was conducted to generate surface and subsurface structural lineaments. Shaded relief from digital elevation models, principal component analysis of Landsat-8 data, and ALOS/PALSAR images were utilized to create lineament maps. Airborne magnetic data were employed to reveal subsurface characterizations. The study area has undergone various tectonic activities, resulting in complex structures. Multiple fault trends and fractures were identified, including the NW–SE (Red Sea-Gulf of Suez) trend, the NE–SW trending Syrian arc trend, the N–S trending East African trend and the WNW–ESE trend. By analyzing the tectonic features of the Bir El-Qash area, this study provides insights into the geological history and evolution of the Eastern Desert of Egypt.
Remote sensing and geophysical data integration have been integral to accurately mapping lineaments. By utilizing remote sensing technologies such as digital elevation models and Landsat-8 data, alongside geophysical data such as airborne magnetic data, researchers can gain a comprehensive understanding of both surface and subsurface tectonic features. This integrated approach has proven to be a powerful tool in geological exploration and analysis1.
Lineament mapping is an important aspect of geological studies as it helps in identifying and understanding the structural elements of the subsurface. In recent years, remote sensing technologies have emerged as a valuable tool for lineament mapping. In particular, shaded relief maps were developed from the digital elevation models (DEMs) used for identifying the lineaments. DEMs enhanced the lineaments by using eight hill shade maps with different sun angles. Combining the first group of hill shade maps (East maps) into the final map also combines the second group of hill shade maps (West maps) into the final map. Landsat-8 data and ALOS/PALSAR images have been utilized, and the principal component analysis technique has been applied to both data to create lineament maps.
The main objective of an aeromagnetic survey is to identify the subsurface geological setting based on variations in the earth’s magnetic field that appear as anomalies as a result of the contrast in the magnetic susceptibility between the basement complex and the sedimentary cap2. The magnetic method is used in different applications, such as locating faults, folds, shear zones, contact, and determining geological structures that may play a role in mineral exploration. Magnetic maps of the area are then used to identify tectonic settings and structures within the sediments through various advanced interpretation methods like filtering, power spectrum analysis, tilt derivatives, Euler deconvolution, and source parameter imaging3,4.
The area under investigation is located in the central part of the Eastern Desert and lies between longitudes 33° 14′ 19.47″ E and 33° 50′ 37.65″ E and latitudes 25° 33′ 40″ N and 25° 54′ 41.11″ N (Fig. 1). The study area has many natural resources that allow sustainable development plans; this makes the area suitable for comparative analysis of available data. The current study aims to exploit remote sensing and airborne magnetic data to determine the structural lineaments that may have a role in mineral exploration in the study area.
Location map of the Bir El-Qash area, Central Eastern Desert, Egypt (using Arc GIS 10.8 and ENVI 5.3 software).
The area was intensively investigated by many workers:5,6,7,8,9,10,11,12,13,14,15,16. The regional rock units and the observed structures investigated in the study area are shown (Figs. 2, 3). The lithostratigraphy of the study area is divided into:
Geological map of the Bir El-Qash Area, Central Eastern Desert, Egypt after7,23 (using Arc GIS 10.4).
(A) Geologic fault trends and fractures derived from the geologic map of Bir El-Qash area, Central Eastern Desert, Egypt7 and (B) Rose diagram of geologic lineaments of the study area.
A) Late Proterozoic rocks include the (i) Ophiolite group, which is made up of serpentine talc carbonate and other related rocks; Undifferentiated metavolcanics that are intermediate to acidic; Shelf sediments that have been changed through metamorphosis; and Pyroclastics. The ophiolitic rocks appear as a long belt stretching NW–SE with generally moderate to high relief 14; (ii) Hammamat Clastics are essentially unmetamorphosed conglomerates, greywackes, sandstone, and siltstone. The basins include deformed molasse sequences (Hammamat Group) with thicknesses of several thousand meters. More than twenty-five Hammamat basin sites have been discovered in Egypt’s Eastern Desert. The (CED) basins of Wadi Hammamat, Zeidun, El Qash, Arak, Zeidun, Um Esh-Um Seleimat, Queih, El Mayah, Kareim, and Igla are the largest and most well-studied examples11,12,17,18,19,20; (iii) Younger granitic rocks: are represented by weakly deformed calk-alkaline granitic rocks that are composed of pink coarse-grained biotite perthite syenogranite, and overlain by pink grey, coarse- to medium-grained equigranular biotite monzogranite and alkaline generally undeformed granitic to alkali-feldspar granitic rocks; (iv) Natash volcanics, which are composed of undeformed volcanic rocks that range from basic to acidic alkaline in composition; (v) Dokhan volcanics; these rocks mostly unmetamorphosed that are well known as purple colored rocks (These rocks are named as imperial porphyry and used as antique in the Romanian age). These rocks are represented by rhyodacite, and plagioclase-phyric andesite subvolcanics21; (vi) Post-Hammamat rocks, which are primarily represented by felsite and quartz veins that overlain by a series of trachyte plugs and sheets.
B) The Cretaceous rocks in the examined area are subdivided into three units, arranged geochronologically as following from the oldest to the youngest: (i) Taref Formation: is composed of red, brown, medium- to coarse-grained ferruginous cross-bedded sandstone and pebble beds; (ii) Quseir Formation: is made up of coastal varicolored shale, siltstone, and flaggy sandstone, with mixed marine and freshwater fossils of gastropods, pelecypods, plants, and vertebrates; (iii) Dawi Formation: is built up of three phosphate horizons separated by beds of marls, shale, and Oyster limestone with flint. The lower horizon of phosphate is known as the Hamadat or Abu Shigeila horizon, which is followed by a section of variegated shales of non-marine to marginal marine origin. The middle horizon is the best-developed and most consistent of the phosphate horizons that are exploited in the Duwi and Hamadat mines. The phosphatic rock is dark in color and has silicified phosphatic nodules. The upper phosphate horizon named Atshan horizon is explored in the Atshan and an-Nakhil mines in the Quseir area; (iv) Lower part of Dakhla Formation is composed of dark gray shallow marine marl and shale with limestone intercalations.
C) Cenozoic rocks are composed of the following units, from bottom to top: (i) Upper part of the Dakhla Formation; (ii) Tarawan Chalk Formation; this Formation is made up of marl and marly limestones bearing Gryphaea vesicularis; (iii) Esna Formation; is made up of grey laminated shales of about 47–50 m thick in Gabal Duwi and Safaga consequently. The formation ranges in age from the late Paleocene to early Eocene times; (iv) Thebes Formation, which consists of fossiliferous limestone (Oyester beds) with flint and chert concretions interbed. The Thebes Formation reveals the following three-stage depositional history22: First stage: basin-wide pelagic carbonate deposition with thin, shallow-water facies near the basin margin; second stage: gradual shallowing and establishment of widespread benthonic communities during intermittent pelagic deposition; simultaneous with tectonic uplift of block faults in the east, which produced many small carbonate platforms separated by deeper basins and third stage: abrupt lowering of sea level, the establishment of shallow-water organic buildups within the basin, and exposure and erosion of uplifted fault blocks to the east; v) Pliocene deposits (Issawiya Formation): consist of fluviatile siltstone, sandstone, and claystone; vi) Quaternary deposits, raised beaches, alluvial fans, and sand covered by wadi deposits.
Structurally, the investigated area was subjected to different tectonic movements, giving rise to some complex structures12; (i) The area was developed through three compression deformation phases, followed by a fourth phase of extension and wrenching. D1: includes SW-ward and their accompanied NW–SE imbricate, D2: is represented by N–S, D3: is manifested by NE–SW trending, and D4 is of transtensional nature and is represented by several episodes of normal and strike-slip faulting; (ii) The normal faults have three principle trends, NE–SW, NW–SE, and N–S, which controlled the intrusion of post-Hammamat felsites, post-granite dykes, and post-tectonic granites, respectively; (iii) the strike-slip faults include NNW-dextral and NNE-sinistral faults dissecting every previously developed structure.
The main objective of this study is to extract surface geologic structures within the study area. A range of remote sensing data and techniques were employed to evaluate the most efficient method for achieving this goal. Specifically, Digital Elevation Models (DEM), Landsat-8, and ALOS/PALSAR data were utilized.
The SRTM project has produced digital elevation models (DEMs) across a large portion of the globe, from 56° south to 60° north. This has resulted in the creation of the most detailed and comprehensive high-resolution topographic database available on Earth24. The SRTM DEM data can be downloaded from http://earthexplorer.usgs.gov.
To distinguish topographic features from the DEM (Fig. 4), the research area’s maximum elevation is found in the northeastern, eastern, and southwestern regions, while the lowest altitude is found in the western parts.
Digital Elevation Model (DEM) of Bir El-Qash area, Central Eastern Desert, Egypt.
This study utilizes data from the OLI and TIRS instruments on the Landsat-8 sensor. Only seven bands of the OLI data were employed, which measure reflected radiation in the VNIR and SWIR regions (bands 2–7) with a spatial resolution of 30 m and the panchromatic band 8 with a spatial resolution of 15 m (Table 1). The deep blue visible channel (band 1) and shortwave infrared channel (band 9) were not utilized as they are primarily designed for water resources, coastal zone investigation, and detecting cirrus clouds, respectively. The thermal bands of the Landsat-8 (TIR), which collect emitted radiation in two bands (10 and 11), were not used in this study. Two cloud-free images were obtained for the same date from the US Geological Survey website (http://earthexplorer.usgs.gov) and are Level 1 Terrain-corrected (L1T) products. The images were geometrically corrected, and a mosaic of VNIR-SWIR data was constructed and preprocessed using ENVI 5.1 software.
Pre-processing is an essential step in this study as it minimizes sensor, solar, atmospheric, topographic effects and distortion for surface reflectance analysis25. The FLAASH algorithm was utilized to correct the atmospheric effects and convert radiance data to reflectance in the Landsat-8 OLI data. The process was carried out by utilizing the algorithm on radiance data that had been radiometrically calibrated in the format of band interleaved by line (BIL)26.
In this research, the Phased Array Type L-Band Synthetic Aperture Radar (PALSAR) Fine Beam Single (FBS) Polarization Mode (HH) was also employed due to its unique radar wavelength (23.62 cm with 1.27 GHz) (Table 1). SAR images are more effective than optical images in identifying subsurface features as they have a longer wavelength, which allows for subsurface penetration27,28. The ability of waves to penetrate a surface layer is related to the amount of scattering loss present in the layer. Factors that affect this scattering loss include the grain size being smaller than 10% of the wavelength of the incident wave and the moisture content being less than 1%. The skin depth of the radar L-Band was calculated at ~ 1.5–2 m in non-vegetated dry sand areas, which allows for the detection of buried (invisible) structural features.
In this study, Digital Elevation Models (DEMs), the principal component band (PC1) of Landsat-8 imagery, and the PALSAR band were utilized to extract lineaments. To identify linear topographical features within the Digital Elevation Model (DEM), eight shaded relief images were generated by illuminating the image from eight distinct directions.
The primary technique commonly employed for remote sensing data involving multivariate statistics and dimensionality reduction is known as principal component analysis (PCA). PCA is utilized to generate uncorrelated bands, separating noise elements and reducing the spectral complexity of the data. In order to get the lines of the study area, principal component analysis (PCA) was used to change the VNIR and SWIR bands of Landsat-8 (OLI).
The magnetic measurements were carried out in 1982 and recorded in Sheet No. 65. The collaboration between the Egyptian General Petroleum Corporation, the Egyptian Geological Survey and Mining Authority, and the Aero-Service Division of the Western Geophysical Company of America provides aeromagnetic data in the research region. This partnership carried out an airborne magnetic survey over a large portion of Egypt’s Eastern Desert, as well as certain portions of the central Western Desert, in 1982 to provide data to help in the identification and assessment of the region’s minerals, petroleum, and groundwater resources29. The traverse lines were spaced 1.5 km apart and ran in the N45°E direction. The tie lines were roughly 10 km apart and perpendicular to the traverse lines. The separation between stations was 92.65 m, and the aircraft’s average speed ranged from 222 to 314 km per hour, with a mean terrain clearance of 120 m. The aeromagnetic data were subjected to different processing techniques, including reduction to magnetic pole (RTP), 2D average power spectrum, regional and residual separation, vertical and tilt derivative, 3D Euler deconvolution, Enhanced horizontal gradient amplitude, and interpreted structural lineaments.
The RTP technique is a valuable method used to eliminate the influence of magnetic inclination on magnetic data. By transforming inclined magnetic data to a hypothetical vertical field, RTP eliminated anomaly asymmetry caused by inclination and enabled the precise localization of anomalies immediately above their causative sources. According to30, the field reduced to the pole at a fixed point above the measurement plane in the frequency domain is given by:
where, I = Geomagnetic inclination, Iɑ = Inclination for amplitude correction (never less than I), and D = geomagnetic declination, θ = wavenumber direction, i = imaginary component.
Reduction to the pole has an amplitude component (the sin term) and a phase component [the icos (I).cos (D + θ) component]. An amplitude inclination of 90º causes only the phase component to be applied to the data (no amplitude correction), and a value of 0º (zero) causes phase and amplitude components to be applied over the entire range31.
The application of spectral analysis to gravity and magnetic data has been widely utilized during the last few years to determine the depth of specific geological characteristics, such as the magnetic basement32. The power spectrum of magnetic data is just the product of the two-dimensional Fourier transform of the map and its complex conjugation. It is a function of wavelength in both the x and y directions. While traditional depth estimation methods examine individual sources, spectral analysis techniques can provide depth information on networks of anomalies, such as collections of bodies. Since the depth of a single anomaly is less significant, it might not produce results as precise as an average depth. The spectrum analysis method is typically used because it is more accurate than other depth computation methods at determining the average depth of the basement. To calculate these depths, one can utilize the slope of a line that has been fitted to any linear section of the curve using equation33.
where, h = depth S = slope of the log (energy) spectrum. These estimates can be used as a rough guide to the depth of the gravity or the magnetic source.
Filtering magnetic data is a necessary step before analyzing and interpreting the RTP map. The main objective of the filtering procedures is to prepare the data set and present it in a way that makes it simpler to understand the significance of magnetic anomalies in relation to their geological origins34. A method called regional-residual separation can be applied to observed magnetic data to separate it into two parts: shallow sources that are local in nature and regional sources. This technique allows for a better understanding of the magnetic data by separating it into distinct components. In the present study, the Fast Fourier Transform (FFT) was applied to the total magnetic intensity data, using the Oasis Montaj package31 to explore the frequency content of these data. The process of regional-residual separation was conducted by applying a Gaussian filter with a cutoff wavelength of 15 km in order to extract the magnetic signature. The Gaussian filter possesses the capability to do either low-pass or high-pass filtering, based on the specific parameters employed during the calculations.
The vertical derivative (vertical gradient) is a good method for resolving anomalies over individual structures in total magnetic intensity data and, more importantly, suppressing the regional content of the data35. It also makes anomalies smaller in width and matches the causative body more closely. Vertical derivative transforms are intended to facilitate the interpretation of reduced-to-pole magnetic data. They are enhancement techniques that amplify the shorter wavelength anomalies relative to those with longer wavelengths. The nth-order vertical derivatives for a given potential field (x, y, z) are calculated as follows:
where F refers to the Fourier transform operator, F−1 refers to the inverse Fourier transform operator, and
where kx and ky are the wavenumbers in the x and y directions, respectively. Clearly, multiplying the transformed potential by k to any power will magnify short-wavelength features of the potential field typically associated with near-surface sources while attenuating long-wavelength components.
The tilt derivative method (TDR) is employed to enhance the edges of magnetic sources. The TDR filter was suggested and developed by35,36,37. The vertical derivative is divided by the total horizontal derivative to calculate this filter38 as follows:
The tilt angle responses can range from positive values over the source to zero at or near the edge to negative values outside the body. When attempting to determine the relative contrast in magnetization, this sign variation is particularly important because of its ability to provide more details39. The TDR approach brings out characteristics that are near the surface in magnetic fields.
Euler deconvolution is a data enhancement technique for estimating the location and depth of a magnetic anomaly source. It relates the magnetic field and its gradient components to the location of the anomaly source with the degree of homogeneity expressed as a structural index, and it is a suitable method for delineating anomalies caused by isolated and multiple sources40. Euler deconvolution is expressed in Eq. 6.
Applying Euler’s expression to profile or line-oriented data (2D source), the x-coordinate is a measure of the distance along the profile, and the y-coordinate is set to zero along the entire profile. Equation 6 is then written in the form of Eq. 7
where (x0, z0) is the position of a 2D magnetic source whose total field T is detected at (x, z). The total field has a regional value of B, and n is a measure of the fall-off rate of the magnetic field. n is directly related to the source slope, is referred to as the structural index, and depends on the geometry of the source40.
Edge enhancement techniques are widely employed to accentuate geologic structures resulting from magnetic field variations41,42,43. A lot of these advanced methodologies for this purpose is extensively documented in existing literature44,45.
In this research endeavor, a cutting edge filter based on the derivatives of the horizontal gradient of magnetic data was introduced by46 to delineate the lateral boundaries of geologic lineaments. This filter, known as the Edge Enhancement based on Horizontal Gradient Amplitudes (EHGA) operator, can be formally represented as follows
where, Ʀ is real part, p is a constant greater or equal to 2, (left(frac{partial THG}{partial y}right)+left(frac{partial THG}{partial y}right)+left(frac{partial THG}{partial z}right)) directions of the total horizontal gradient (THG) is approximated by:
Lineaments are tonal or topographical features on the land that reflect zones of structural weakness47. The Central Eastern Desert has been intensively investigated by many workers21,48. The identification of lineaments has been thoroughly explored by numerous researchers49,50,51. The identification of linear features from remotely sensed data is carried out by applying two main techniques: (1) manual digitizing techniques47. (2) Algorithms and software of computers52. In this research, Digital Elevation Models (DEM), the first principal component band (PC1) of Landsat-8, and the PALSAR band were utilized to extract lineaments.
To identify linear topographical characteristics in the Digital Elevation Model (DEM), eight shaded relief images were created by illuminating the image from eight different directions. By using 0 azimuth of solar, we construct the first image of shaded relief. The light directions in the last seven shaded relief images were varied by 45°, 90°, 135°, 180°, 225°, 270°, and 315°.
The input parameters’ values determine the number and size of the extracted lineaments. These values are represented by optional digits in the PCI Geomatica software’s LINE module. The module’s algorithm comprises three steps: edge detection, thresholding, and curve extraction.
However, the LINE module itself is responsible for extracting lineaments from an image and transforming these linear features into vector form, employing six optional parameters. This step involves combining four shaded relief images into a single image53 (Fig. 5A, B). This stage entails the creation of four shaded relief images through GIS techniques, achieved by superimposing four shaded relief images to produce a composite image with varying illumination directions (0°, 45°, 90°, and 135°) (Fig. 6A) and another image with multiple lighting directions (180°, 225°, 270°, and 315°) (Fig. 6C). The Geomatica 10.3 software is used to combine the images for automated lineament extraction across the study area. The generated lineaments from both images with different sun angles (0°, 45°, 90°, 135°, 180°, 225°, 270°, and 315°) show an NW–SE trend (Fig. 6B, D).
(A) Combination of shaded relief images of the Bir El-Qash region, Central Eastern Desert, Egypt, with sun angles of 0°, 45°, 90°, and 135°. (B) Shaded relief image generated by integrating four shaded relief images of the Bir El-Qash region in Egypt’s Central Eastern Desert with sun angles of 180°, 225°, 270°, and 315°.
(A) Lineament map generated by integrating shaded relief images with sun angles of 0°, 45°, 90°, and 135° and (B) Rose diagram of the lineaments of the Bir El-Qash region, Centr of Eastern Desert, Egypt. (C) Lineament created by integrating of shaded relief images with sun angles of 180°, 225°, 270°, and 315° and (D) Rose diagram of the automated lineaments of Bir El-Qash region, Central of Eastern Desert, Egypt.
PCA is a widely used multivariate statistical method. The principal component analysis was utilized in this study to separate noise components and reduce the spectral dimensionality of the data, resulting in uncorrelated output bands54,55,56. The PCA technique was applied to OLI data, and the first principal component band (PC1) was chosen for lineament extraction in the study area (Fig. 7A). The results of the principal component analysis (PCA) showed that the first band had a significantly higher variance (96.71%) compared to the other bands (Table 2). Therefore, the PC1 band was selected for the automatic extraction of lineaments in the study area (Fig. 8A).
(A) Principal Component Analysis (PC1) image of the study area. (B) Grayscale image of AlOS/PALSAR SAR data.
(A) Automatic lineament map of PCA1 image of the study area. (B) Rose diagram of automatic Lineament map. (C) Automatic lineament map of PALSAR image of the study area. (D) Rose diagram of automatic Lineament map.
The lineament analysis of the study area indicates that it has undergone significant tectonic activity, resulting in extensive faulting with varying orientations and magnitudes. Notably, the WNW–ESE and NW–SE trending Najd fault system stands out as a prominent feature (Fig. 8B).
The ability to outline structural lineaments in radar data has been improved by applying four directional filters (0º, 45°, 90º, and 135°) to the backscatter PALSAR data, which enhances linear features for the extraction process of key structural trends (NE–SW, E–W, and NW–SE, respectively). Additionally, four DEM-derived shaded relief images were extracted from the ALOS/PALSAR DEM which have 12.5 m spatial resolution with four different directions (0º, 45°, 90º, and 135°) to aid in the lineament extraction process with the ALOS/PALSAR data (Fig. 7B). It has been established that when this approach is used together with optical images, it can be trusted as a dependable way to gather information about lineaments and faults57. The study area was analyzed using an ALOS/PALSAR image to extract lineaments automatically. Automated lineament extraction was performed using algorithms and computer software e.g.58,59,60,61,62. In this research, the PCI Geomatica software was utilized to automatically extract lineaments from the PALSAR image, resulting in vector segments that depict linear features. (Fig. 8C).
The lineament analysis carried out in the research area has provided valuable insights into the tectonic processes that have shaped the geological features of the region. The findings reveal that the area has undergone complex tectonic activity, resulting in extensive faulting with diverse orientations and magnitudes. The identified lineaments include prominent trends such as the WNW–ESE oriented and NW–SE (Red Sea-Gulf of Suez) trend, as well as less conspicuous trends such as the NE-SW oriented Gulf of Aqaba trend. The contradiction arises from the ALOSE/PALSER image, which reveals a WNW–ESE trend. This inconsistency with other study data can be attributed to the unique ability of ALOSE/PALSER to unveil hidden features beneath natural phenomena, employing high-texture lineaments through backscattering radiations. Additionally, it’s noteworthy that remote sensing data primarily focuses on surface structures, in contrast to magnetic data, which is concerned with both shallow and deep structures. The interpretation of these lineaments allows for a better understanding of the geological history and evolution of the study area, which may have implications for geological exploration and natural resource management (Fig. 8D).
The examination of a magnetic map (Fig. 9a) for a specific study area. The map reveals two primary magnetic zones within this area. The first zone, found in both the southern and northern regions of the study area, displays a significant magnetic anomaly with a peak value of 42 nT. This high magnetic anomaly is characterized by high-amplitude variations and is attributed to the presence of high-rise basement outcrops in these areas. The high magnetic values likely result from geological features or materials in the Earth’s crust, such as rocks with a strong magnetic susceptibility. The second zone, in contrast, is situated in the central and eastern parts of the study area and is associated with a low magnetic anomaly. This low magnetic anomaly reaches a minimum value of -115 nT. The geological composition of this area is described as sedimentary and acidic rocks. It is also noted that this region trends in an NW–SE direction.
(A) Total Magnetic Intensity Map (TMI), and (B) Reduced to northern magnetic pole (RTP) map of the of Bir El-Qash area, Central Eastern Desert, Egypt.
The total intensity magnetic map (Fig. 9A) was processed by applying 2-D wave number filtration to adapt it to the research area’s northern magnetic pole using its inclination and declination. The application of the International Geomagnetic Reference Field (IGRF) technique indicates that the total magnetic intensity is 42,425 gammas in the study area and the angle of the declination and inclination are 39.5 degrees and 2 degrees, respectively. The RTP—magnetic map is created using these parameters to produce the reduced to the north magnetic pole (Fig. 9B). The reduction to the northern magnetic pole has been used to overcome the inclination problems and to locate the magnetic anomalies directly above the causative sources.
Because the inclination effects on the magnetic field at this point have been removed, a general look at this RTP map in comparison to the original total aeromagnetic intensity map shows that the magnetic anomalies appear to have moved northward. The magnetic anomalies on the reduced to the northern magnetic pole map (Fig. 9B) are slightly different from those on the total magnetic intensity (TMI) map (Fig. 9A). In the study area, a magnetic anomaly with a maximum amplitude value of 84 nT has been observed. This anomaly is characterized by positive magnetic anomalies, which are most pronounced in the southern part of the study area. They extend towards the central region and consist of multiple irregular and elongated closures. Furthermore, there’s a notable positive magnetic anomaly in the northwest part of the study area, which is elongated and runs in a northwest-southeast direction. Geologically, this area exhibits a diverse range of rock types, including the Quseir Formation, Taref Formation, Hammamat clastics, granites, wadi deposits, metasediments, and ophiolitic serpentinite.
The significant magnetic anomalies might be caused by the existence of subsurface basic intrusions with high magnetic content. The research region exhibits the presence of felsite dykes along Wadi El-Qash in N–S and NE–SW directions that penetrate the basement rocks and molasse deposits and there are felsite sills situated in the northern region of the Zeidun basin. The intrusion of felsite results in the formation of an alteration zone inside the molasse sediments. The Arak basin is affected by NW-trending dykes composed of porphyritic microgranite16.
The reduced-to-the-pole magnetic amplitude map demonstrates a notable negative magnetic anomaly with a minimum value of −155 nT. This negative anomaly is primarily located in the western and northeastern sections of the study area. The geological characteristics of these regions include a diverse array of rock types, such as the Taref Formation, Quseir Formation, Duwi Formation, Dakhla Formation, Tarawan Formation, wadi deposits, metavolcanic, and metasediments.
The primary orientation of both positive and negative magnetic anomalies on the map follows two main trends: northwest–southeast (NW–SE) and north–south (N–S). This alignment of magnetic anomalies suggests the presence of structural features within the study area.
The analysis of the two-dimensional power spectrum curve has revealed the presence of two distinct linear segments (Fig. 10). Each of these segments corresponds to specific characteristics of the sources of the observed phenomena.
The first segment is associated with deep and/or wide sources. It is characterized by low frequencies, indicating long wavelengths. In other words, the sources contributing to this segment are located at greater depths within the Earth’s subsurface.
The second segment, on the other hand, is attributed to shallow sources. It exhibits high frequencies, representing short wavelengths. Shallow sources are typically closer to the Earth’s surface.
Power spectrum technique shows the depth of magnetic sources of Bir El-Qash area, Central Eastern Desert, Egypt.
The deep magnetic sources are characterized by wavenumbers ranging from 0.00 to 0.152 cycles per k-unit, and the average depth for regional sources is approximately 3 km. Meanwhile, the shallow sources are associated with wavenumbers ranging from 0.152 to 1.30 cycles per k-unit at the first level of residual, and the average depth for residual sources is approximately 0.58 km.
After removing the residual effects from the RTP map, the regional magnetic component map (Fig. 11A) at the assigned interface is generated. Magnetic anomalies in the region can be separated into two groups. The first group is located in the southern and northwestern parts of the study area. These anomalies have a high range of magnetic values, between 0 and 43 nT. The second group is located in the northeastern and western parts of the study area. These anomalies have a low range of magnetic values, between 0 nT and −124 nT. The main trend of the positive and negative anomalies is NW–SE and N–S.
(A) Regional anomaly map and, (B) Residual anomaly map of airborne magnetic data of Bir El-Qash area, Central Eastern Desert, Egypt.
Geophysicists use residual maps to identify local features that are often obscured by broader features in the field63. The residual magnetic anomaly map (Fig. 11B) shows the magnetic field distribution at a shallow depth after removing the regional influence. As a result, the map shows a high frequency of both high and low magnetic anomalies with a limited area of coverage and variations in amplitude and direction. These variations indicate different anomalies, sources, depths, compositions, and structure configurations. The residual magnetic anomaly map (Fig. 11B) shows a pattern of alternating negative and positive magnetic zones, with a general direction of NW–SE, and N–S. This variation in the predominant structural trends (NW–SE and N–S) is attributed to a variety of tectonic processes within the investigated region12.
In data analysis, the first vertical derivative is frequently employed to find edges and geologic limits64,65. To improve the shallow sources of geologic features of the data, the vertical derivative is commonly applied to RTP magnetic data. The first or second vertical derivative is often used to emphasize gradients near the edges of shallow magnetic sources. This allows them to be used to find the boundaries of magnetic bodies and identify sources at shallow depths66.
The analysis of the first vertical derivative technique of airborne magnetic data is shown in Fig. 12A. The NW–SE trending anomaly system is located exactly in the area of an NW–SE striking fault system confined in the major fault lineaments. As a result, the long, narrow, and linear negative magnetic anomalies in particular parts of the area have a strong spatial relationship with fault structures. In other parts of the study area, the direction of elongated and semi-round positive and negative magnetic anomalies aligns significantly with the primary structural trend that affects the research region.
(A) First order anomaly map of vertical derivative technique and, (B) Tilt Derivative (TDR) technique of Bir El-Qash area, Central Eastern Desert, Egypt.
It is feasible to discern numerous unique lineaments in the TDR anomaly map that are very difficult to notice in the total magnetic anomaly field map. The identification of magnetic lineaments is depicted in the TDR anomaly map (Fig. 12B), where the TDR’s zero contour is shown in the black line on the tilt derivative map. This map reflects a series of structural lineaments, the NN–SE and NS trends, which represent the main structural trends in the study area and are associated with the Red Sea-Gulf of Suez trend and the East African trend, respectively.
To create a map that shows the locations and depths of magnetic sources in a two-dimensional grid, 3D Euler deconvolution is used67. Any geologic model (faults, magnetic contacts, dykes, sills, etc.) can be used with the Euler Deconvolution method, and it is unaffected by geomagnetic inclination and declination or magnetic remanence68. The depth of a magnetic source is crucial for analyzing underlying structures using geological and geophysical methods55,69,70. Only when the depth uncertainty of the calculated depth falls beneath a predetermined threshold and the solution’s position falls within a predetermined range of the data window’s center is the solution deemed to be valid71.
Utilizing the RTP magnetic grid of the research region, the depth estimate by 3D Euler deconvolution has been applied31. The available RTP anomaly map was subjected to this method using the structural indices (0, and 1.0) and window size equal to 20, as well as a tolerance value of 5.
The simple model of a magnetic field produced by a geological contact or fault structure, concentrated along the borders of the anomalies, is given a structural index of 0. The calculated depth estimates with a structural index of 0 are displayed in (Fig. 13A). According to the Euler deconvolution approach, the basement structural depths in the investigation area range from less than 400 m to more than 2800 m. Additionally, the research shows that the lineament patterns are mainly oriented in the NNW–SSE, NE–SW, and N–S directions, which represents the prevailing structural pattern and tectonic trend of the studied area.
(A) Structural index of (0) and (B) Structural index of (1) based on the Euler deconvolution technique of Bir El-Qash area, Central Eastern Desert, Egypt.
The magnetic dyke or sill solution for structural sources (dykes and sills) is built using RTP magnetic data (Fig. 13B). Using the structural index model (SI = 1), this model provides solutions for the positions and depths of the dykes and sills that intruded the various basement rocks at great depths (Fig. 13B). The depths obtained range from 400 to 2800 m. The analysis of these solutions reveals that the lineament patterns are generally oriented in the NNW–SSE, NW–SE, and N–S directions. These dykes or sills are related to felsite dykes that penetrate the basement rocks and molasse deposits and are situated along Wadi El-Qash and along the northern region of the Zeidun basin16.
The EHGA technique can generate well defined edges for magnetic sources. The peaks of the EHGA are located directly over the source borders, and this filter generates more reliable responses. The EHGA filter is applied to the aeromagnetic data and clearly outlines and refines the lineaments of the studies region.
The Enhanced horizontal gradient amplitude map (Fig. 14) shows alternating between positive (high) and negative (low) values. The positive values reach up to 0.79 rad, while the negative values reach to -1rad. The areas with high amplitude suggest stronger magnetic gradients. This could be indicative of geological feature such as faults, changes in rock type or mineral deposits that exhibit significant contrasts in magnetic propereties. The analysis and interprtation of EHGA shows that the study area is highly deformed and the main structural trending the NW–SE, NNW–SSE and N–S directions.
Enhanced horizontal gradient amplitude of Bir El-Qash area, Central Eastern Desert, Egypt.
Finding linkages between the data and the region’s subterranean structure is a key aspect of interpreting potential field data in geology. In several areas of geology and geophysics, structural difficulties are frequently found through analyzing structural patterns. A significant relationship between the pattern, strength, and direction of magnetic anomaly trends was discovered66. This is due to the fact that the presence and strength of magnetic contrasts in the rocks have a significant impact on the visibility of faults on magnetic maps.
The region under investigation is primarily composed of igneous and sedimentary rocks that have been impacted by various types, trends, and positions of regional and local geological structures. The interpreted magnetic structure maps of the study area were constructed using surface geologic information, the interpretation of magnetic maps, and the research results of the structural trend analysis. Several magnetic methods were used to extract the subsurface structures in the study area, such as tilt derivative structural elements (Fig. 15A), shallow-seated structural elements (Fig. 15C), and deep-seated structural elements (Fig. 15E).
(A) Structural elements resulting from tilt derivative map, (B): rose diagram of structural elements of tilt derivative map, (C) shallow structural elements, (D) rose diagram of shallow structural elements and, (E) deep-seated structural elements and, (F) rose diagram of deep-seated structural elements of the study area.
The analysis of these diagrams reveals three primary structural trends of varying degrees of strength and duration. According to magnetic data, these three trends (NW, NE, and N-S) are the most significant structural trends affecting the area under study. As illustrated below, each of the major trend categories is briefly discussed in descending order. The trend analysis of the interpreted magnetic structural maps is used to explore their relationship in the present area and obtain an idea about the main features of its structural pattern. Rose diagrams (Fig. 15B, D, F) were built to varying degrees from available data to aid in the understanding of structural trends affecting the study area. These diagrams reveal two major structural trends: the Red Sea-Gulf of Suez trend (NW–SE) and the Syrian arc trend (NE-SW). Meanwhile, the Meridional or East African (N–S) trend appears as a minor trend. Each of the major and minor trend categories is described briefly, as shown below.
The magnetic anomaly maps reveal the NE-SW trend as the dominant trend. The analysis of the rose diagrams of the regional, residual, and tilt derivative structure maps (Fig. 15) indicates that this trend has similar azimuth and intensity significance and is considered a strong trend on the deep- and shallow structural elements maps (Fig. 15). This trend is known as the Aualitic trend and is linked to rift tectonics in the Gulf of Aqaba73. He also noted that it is less prominent near the stable shelf and characterizes several of Egypt’s fold systems in the middle (Syrian arc system). The Jordan-Dead Sea-Aqaba rift has been active as a strike-slip fault since the Oligocene period and has evolved into a graben since the late Pliocene74.
It is the most prevalent one in the area at all depths presented by the geologic and processed maps. In addition to its surfacial significance, the NW–SE trend is a significant deep-seated tectonic trend. This trend appears to be a dominant trend of both deep and shallow structures in the Nile Valley, Sinai, and the Gulf of Suez75. He further said that the large extent of this direction, as well as its occurrence in different depth zones of the crust, give the idea that it is associated with an old tectonic movement that has lately been reactivated, most likely with the movement created in the Aqaba direction. The peaks of the NW trend were identified as a major trend in the shallow-seated structural elements map (Fig. 15C); however, this trend was identified as a moderate tilt derivative structure map (Fig. 15A) and a deep-seated structural elements map (Fig. 15E).
This trend encompasses all trends that run from N 11.25°E to N 11.25° W and is known as the “East-African trend73. The N–S trend is specifically highlighted on the rose diagram of deep-seated structural elements (Fig. 15F) indicating its significance in the deeper parts of the study area. However, it’s worth noting that this north–south trend is not prominently observed in the shallow depths (Fig. 15B, D) It is a prominent Precambrian direction as well as a tensional late Tertiary direction75 and is considered the oldest fracture system in the Eastern Desert63. The Nubia trend, which is found in the basement rocks of the Eastern Desert, Sinai, and inliers in southern Egypt, is thought to be the northern extension of the Mozambique belt76. The (N–S) trend is also visible as a main trend on the surface structure map.
This paper deals with the evaluation of the surface and subsurface structure of Bir El-Qash Area, Central Eastern Desert, Egypt, from DEM and aeromagnetic data. A digital elevation model (DEM), Landsat-8 Operational Land Imager (OLI) data, and Phased Array Type L-Band Synthetic Aperture Radar (PALSAR) images are the three types of remote sensing data that are used to find surface lineaments. Meanwhile, the subsurface lineaments were detected using aeromagnetic data. The results from the analysis of remote sensing data and field geology revealed four main trends, namely NW–SE, WNW–ESE, NE–SE, and N–S. The results demonstrate that the ALOS/PALSAR SAR image is the most effective data source for extracting surface structures. The striking differences in the results of the previous structural map of the study area arise WNW–ESE trend revealed only in the ALOSE/PALSER image. The divergence from prior research outcomes can be attributed to the unique capacity of ALOSE/PALSER to expose hidden features within natural phenomena, utilizing high-texture lineaments through backscattering radiations. Additionally, it is noteworthy to highlight that remote sensing data primarily focuses on surface formations, in contrast to magnetic data, which encompasses both shallow and deep formations.
Initially, the total intensity aeromagnetic map of the study area was transformed to the northern magnetic pole (RTP) to identify the magnetic anomalies’ location above their underlying sources. Regional-residual separation techniques were applied to the RTP magnetic map to differentiate between shallow anomalies (structures) and deeper ones. This separation was performed in the frequency domain using the Fast Fourier transformation, which is implemented by analyzing the energy spectrum of the magnetic maps. Through power spectral analysis, the RTP anomaly map of the study area was separated into residual and regional magnetic anomalies.
Spectral analysis, and Euler deconvolution are all automatic source depth estimation techniques that may be used to map the basement depth. The analysis of spectral data from the research area’s aeromagnetic data showed two major magnetic anomaly sources depths, deep magnetic sources, and shallow magnetic sources. The deep source anomalies are located at a depth of 3km, which corresponds to the magnetic basement depth. At a depth of 0.58km. Based on the interpretation of the EHGA, residual, regional, and TDR techniques, the study area reveals pronounced deformation, with predominant structural orientations primarily observed along the NW–SE, NNW–SSE, and N–S.
In summary, the study has identified three major structural trends within the investigated area through the examination magnetic and remote sensing data. The geological framework of the region under consideration is marked by prominent trends, namely, the Red Sea-Gulf of Suez trend, the Syrian Arc trend, and the WNW–ESE trend. These identified trends are considered to be the most dominant and influential tectonic features in the study area. They have played a crucial role in shaping the geological characteristics of the region, including mineral distributions.
In contrast, the Meridional or East African (N–S) trend appears to be of minor significance compared to the Red Sea-Gulf of Suez and Syrian Arc trends. While it may still have an impact on the area’s geological features, it is not as prominent or influential as the other two major trends.
The data presented in this study are available on request from the corresponding author.
Eldosouky, A. M., El-Qassas, R. A. Y., Pour, A. B., Mohamed, H. & Sekandari, M. Integration of ASTER satellite imagery and 3D inversion of aeromagnetic data for deep mineral exploration. Adv. Space Res. 68, 3641–3662. https://doi.org/10.1016/j.asr.2021.07.016 (2021).
Article ADS CAS Google Scholar
Reynolds, J. An Introduction to Applied and Environmental Geophysics, Wiley, Chester, 2nd edition, 710p (1997).
Khalil, A. et al. Subsurface structural characterization as deduced from potential field data—West Beni Suef, Western Desert, Egypt. Arab. J. Geosci. 15, 1627. https://doi.org/10.1007/s12517-022-10900-1 (2022).
Article Google Scholar
Lasheen, E. S. R., Mohamed, W. H., Ene, A., Awad, H. A. & Azer, M. K. Implementation of petrographical and aeromagnetic data to determine depth and structural trend of Homrit Waggat Area, Central Eastern Desert, Egypt. Appl. Sci. 12, 8782. https://doi.org/10.3390/app12178782 (2022).
Article CAS Google Scholar
Dixon, T. H. Evolution of Continental Crust in the Late Precambrian Egyptian Shield. California University Press, San Diego, USA (1979)
Akaad, M. K., Noweir, A. M., Geology and lithostratigraphy of the Arabian Desert orogenic belt of Egypt between latitudes 25◦35′ and 26◦30′N. In A. M. S. Al Shanti, (ed.), Evolution and Mineralization of the Arabian–Nubian Shield, vol. 4. Pergamon Press, Oxford, pp. 127–135 (1980).
Conoco. Geological map of Egypt, Scale (1:500,000), NH 35 SE Bahariya Sheet (1987).
Akaad, M. K., Abu El Ela, A. M. On the clasts of the Atud conglomerates: a reconsideration. M.E.R.C. Ain-Shams University, Egypt. Earth Sci. 4, 1–21 (1990).
El-Gaby, S., List, F. K., Tehrani, R. The basement complex of the Eastern Desert and sinai. In Said, R. (Ed.), The Geology of Egypt. A.A. Balkema, Rotterdam, pp. 175–184, 734p (1990).
Akaad, M. K., Abu El Ela, A. M. Geology and petrochemistry of the Muweilih submarine metabasalts, qift-quseir region, Eastern Desert, Egypt. Egypt. J. Geol. 40(1), 321–349 (1996).
Abd El-Wahed, M. A. Thrusting and transpressional shearing in the Pan-African nappe southwest El-Sibai core complex, Central Eastern Desert, Egypt. J. Afr. Earth (2008).
Abd El-Wahed, M. A. The role of the Najd Fault System in the tectonic evolution of the Hammamat molasse sediments, Eastern Desert, Egypt. Arab. J. Geosci. 3, 1–26 https://doi.org/10.1007/s12517-008-0030-0 (2010).
Khalil, K. I. & El-Shazly, A. K. Petrological and geochemical characteristics of Egyptian banded iron formations: Review and new data from Wadi Kareem. Geochem. Explor. Environ. Anal. 12, 105–126 (2012).
Article CAS Google Scholar
El Bahariya, G. Classification of the Neoproterozoic ophiolites of the Central Eastern Desert, Egypt based on field geological characteristics and mode of occurrence. Arab. J. Geosci. 11(12) (2018).
Abd El-Wahed, M. A. Hamdy, M. M., Ali, K. A., Ayad, A. E., Zeh, A., Zoheir, B., Ren, M., Koutsovitis, P. Reworked Neoproterozoic sub-arc lithosphere beneath the Eastern Desert (Egypt): Evidence from geochemistry and geochronology of metabasites and granitoid gneisses, Gabal Um Gunud area. Int. Geol. Rev. (2023). https:// doi.org/https://doi.org/10.1080/00206814.2023.2207305.
Kassem, O. M., Abd El-Wahed, M., Abudeif, A. M., Abdelhameed, M. Microstructural and strain studies of neoproterozoic rocks encapsulating hammamat molasse deposits in theWadi El-Qash area, Central Eastern Desert, Egypt. J. Afr. Earth Sci. 205, 105009 (2023). https://doi.org/10.1016/j.jafrearsci.2023.105009.
Fowler, A. R. and El-Kalioubi, B. Gravitational collapse origin of shear zones, foliations and linear structures in the Neoproterozoic cover nappes, Eastern Desert, Egypt. J. Afr. Earth Sci. 38, 23–40 (2004).
Abdeen, M. M. & Greiling, R. O. A quantitative structural study of late Pan-African compressional deformation in the Central Eastern Desert (Egypt) during Gondwana assembly. Gondwana Res. 8(4), 457–471 (2005).
Article ADS Google Scholar
Fowler, A. R. & Osman, A. F. Gneiss-cored interference dome associated with two phases of late Pan-African thrusting in the central Eastern Desert, Egypt. Precambrian Res. 108, 17–43 (2001).
Article ADS CAS Google Scholar
Fowler, A. R. & Osman, A. F. Sedimentation and inversion history of three molasse basins of the western Central Eastern Desert of Egypt: implications for the tectonic significance of Hammamat basins. Gondwana Res. 23, 1511–1534 (2013).
Article ADS Google Scholar
Badawy, M. B. M. Application of remote sensing and geographic information systems in geomorphological studies: Safaga-El-Quseir Area, Red Sea Coast, Egypt as an Example. Dissertation zur Erlangung des Grades “Doktor der Naturwissenschaften” im Promotionsfach Geographie am Fachbereich Chemie, Pharmazie, Geowissenschaften der Johannes Gutenberg-Universitt Mainz (2008).
Snavely, P. O., Garrison, R. E. & Meguid, A. A. Stratigraphy and regional depositional history of the thebes formation (lower eocene), Egypt. Ann. Geol. Surv. Egypt 9, 344–362 (1979).
Google Scholar
Abdelkareem, M. & El-Baz, F. Analyses of optical images and radar data reveal structural features and predict groundwater accumulations in the central Eastern Desert of Egypt. Arab J Geosci 8, 2653–2666. https://doi.org/10.1007/s12517-014-1434-7 (2015).
Article CAS Google Scholar
Nikolakopoulos, K. G., Kamaratakis, E. K., Chrysoulakis, N. SRTM vs ASTER elevation products. Comparison for two regions in crete, Greece. Int. J. Remote Sens., 27(21), 4819–4838 (2006).
Young, N. E. et al. A survival guide to Landsat preprocessing. Ecology. 98(4), 920–932 (2017).
Article PubMed Google Scholar
Cooley, T., Anderson, G., Felde, G., Hoke, M., Ratkowski, A., Chetwynd, J., Gardner, J., Adler-Golden, S., Matthew, M., Berk, A., Bernstein, L., Acharya, P., Miller, D. and Lewis, P. FLAASH, a MODTRAN4-based atmospheric correction algorithm, its application and validation. IEEE International Geoscience and Remote Sensing Symposium (2002).
Elachi, C., Granger, J. Spaceborne imaging radars probe `in depth’. IEEE Spectrum, 19(11), pp.24–29 (1982).
McCauley, J. et al. Subsurface valleys and geoarcheology of the Eastern Sahara Revealed by Shuttle Radar. Science 218(4576), 1004–1020 (1982).
Article ADS CAS PubMed Google Scholar
Aero-Service. Final operational report of airborne magnetic/radiation survey in the Eastern Desert, Egypt. For the Egyptian General Petroleum Corporation (EGPC) and the Egyptian Geological Survey and Mining Authority (EGSMA), Aero-Service Division, Houston, Texas, USA, Six Volumes (1984).
Blakely, R. J. Potential theory in gravity and magnetic applications (Cambridge University Press, 1995).
Book Google Scholar
Geosoft Program Oasis Montaj. Geosoft Mapping and Processing System. (Geosoft Inc., Toronto, Canada. Scientific Research Publishing An Academic Publisher Copyright 2006–2017 Scientific Research Publishing Inc (2010).
Garcia-Abdeslem, J. & Ness, G. E. Inversion of the power spectrum from magnetic anomalies. Geophysics 59, 391–401 (1994).
Article ADS Google Scholar
Spector, A. Application of aeromagnetic data for porphyry copper exploration in area of volcanic cover. In The 45th Annual International Meeting of the Soci. of Exp. Geophysicists, Colorado, USA (1975).
Bird, D. Interpreting magnetic data: geophysical corner AAPG Explorer (1997).
Fairhead, J. D. (2009). Potential Field Methods for Oil and Mineral Exploration GETECH/ Univ. of Leeds, pp. 290.
Salem, A., Williams, S., Fairhead, D., Smith, R. & Ravat, D. Interpretation of magnetic data using tilt-angle derivatives. Geophysics 73, L1–L10 (2008).
Article ADS Google Scholar
Salem, A., Williams, S., Fairhead, J., Ravat, D. & Smith, R. Tilt-depth method: a simple depth estimation method using first-order magnetic derivatives. Lead Edge 26, 1502–1505 (2007).
Article Google Scholar
Verduzco, B. New insights into magnetic derivatives for structural mapping. Lead. Edge 23, 116–119 (2004).
Article Google Scholar
Cooper, G. R. J., Cowan, D. R. Enhancing potential field data using filters based on the local phase. Comput. Geosci. 32(10), 1585–1591 (2006).
El Dawi, M. G., Tianyou, L., Hui, S. & Dapeng, L. Depth estimation of 2-D magnetic anomalous sources by using Euler deconvolution method. Am. J. Appl. Sci. 1(3), 209–214 (2004).
Article Google Scholar
Pham, L. T., Le-Huy, M., Oksum, E. & Do, T. D. Determination of maximum tilt angle from analytic signal amplitude of magnetic data by the curvature-based method. Vietnam J. Earth Sci. 40(4), 354–366 (2018).
Article Google Scholar
Pham, L. T., Oksum, E., Do, T. D., Nguyen, D. V., Eldosouky, A. M. On the performance of phase-based filters for enhancing lateral boundaries of magnetic and gravity sources: a case study of the Seattle Uplift. Arab. J. Geosci. 14(2). https://doi.org/10.1007/s12517-021-06511-x (2021a).
Pham, L. T. et al. An improved approach for detecting ridge locations to interpret the potential field data for more accurate structural mapping: a case study from Vredefort dome area (South Africa). J. Afr. Earth Sci. 175, 104099. https://doi.org/10.1016/j.jafrearsci.2020.104099 (2021).
Article Google Scholar
Eldosouky, A. M. & Mohamed, H. Edge detection of aeromagnetic data as effective tools for structural imaging at Shilman area, South Eastern Desert, Egypt. Arab. J. Geosci. 14, 13. https://doi.org/10.1007/s12517-020-06251-4 (2021).
Article Google Scholar
Ekwok, S. E., Eldosouky, A. M., Achadu, O. I. M., Akpan, A. E., Pham, L. T., Abdelrahman, K., & Fnais, M. S. Application of the enhanced horizontal gradient amplitude (EHGA) filter in mapping of geological structures involving magnetic data in southeast Nigeria. J. King Saud Univ. – Sci. 34(8), 102288, https://doi.org/10.1016/j.jksus.2022.102288 (2022).
Pham, L. T., Oksum, E., Gómez-Ortiz, D. & Do, T. D. MagB_inv: A high performance Matlab program for estimating the magnetic basement relief by inverting magnetic anomalies. Comput. Geosci. 134, 104347 (2020).
Article Google Scholar
Williams, R. S. Geological applications, In. Colwell, R. N. (Editors). “Manual of Remote Sensing”, American Society of Photogrammetry, Falls Church, pp. 1667- 1951(1983).
Abou El-Maaty, M. A., Taman, Z. I., Abdel Wahab, W. Characterization of the Atalla Talc-Carbonate rocks, central Eastern Desert, Egypt. Egypt. J. Geol. 56(1), 67–83 (2012).
Koike, K., Nagano, S. & Ohmi, M. Lineament analysis of satellite images using a Segment Tracing Algorithm (STA). Comput. Geosci. 21, 1091–1104 (1995).
Article ADS Google Scholar
Masoud, A., Koike, K. Morphotectonics inferred from the analysis of topographic lineaments auto-detected from DEMs: Application and validation for the Sinai Peninsula, Egypt. Tectonophysics 510, 291–308 (2011).
Obiegbu, O. L., Laake, A., Brabham, P. An estimation of regional geologic structures from the geomorphology to characterize the Wadi systems, southeast Sinai, Egypt, Published by; American Association of Petroleum Geologists and Society of Exploration Geophysics. Current issue; volume 5, number 1 (2016).
Baumgartner, A., Ateger, C., Mayer, H., Eckstein, W., Ebner, H. Automatic road extraction based on multi-scale, grouping and context. Photogrammetric Engineering and Remote Sensing V 65, 777–785 (1999).
Thannoun, R. Automatic extraction and geospatial analysis of lineaments and their tectonic significance in some areas of Northern Iraq using remote sensing techniques and GIS. https://doi.org/10.13140/RG.2.2.20851.99363 (2013).
Cheng, Q., Jing, L. & Panahi, A. Principal component analysis with optimum order sample correlation coefficient for image enhancement. Int. J. Remote Sens. 27, 3387–3401 (2006).
Article Google Scholar
Gupta, R. P. Remote Sensing Geology, 3rd ed.; Springer: Berlin, pp. 180–190, 235–240, 332–336 (2017).
Schowengerdt, R. A. Remote Sensing: Models and Methods for Image Processing, 3rd ed. Elsevier: Burlington, MA, USA; pp. 193–202, 229–243 (2007).
Ganas, A., Pavlides, S. & Karastathis, V. Dem-based morphometry of range-front escarpments in Attica, Central Greece, and its relation to fault slip rates. Geomorphology 65(3–4), 301–319. https://doi.org/10.1016/j.geomorph.2004.09.006 (2005).
Article ADS Google Scholar
Kim, T., Park, S., Kim, M., Jeong, S. & Kim, K. Tracking road centerlines from high resolution remote sensing images by least squares correlation matching. Photogrammetric Eng. Remote Sens. 70, 1417–1422 (2004).
Article Google Scholar
Hung, L. Q., Batelaan, O., Smedt, F. Lineament extraction and analysis, comparison of LANDSAT ETM and ASTER imagery. Case study: Suoimuoi tropical karst catchment, Vietnam. In Proceedings of SPIE, V. 5983, 12p. Remote Sensing Symposium. 2, 724–726 (2005).
Abdullah, A., Akhir, J. M. & Abdullah, I. A comparison of Landsat TM and SPOT data for lineament mapping in Hulu Lepar Area, Pahang, Malaysia. Eur. J. Sci. Res. 34(3), 406–415 (2009).
Google Scholar
Abdullah, A., Akhir, J. M. & Abdullah, I. Automatic mapping of lineaments using shaded relief images derived from digital elevation model (DEMs) in the Maran-Sungi Lembing area, Malaysia. Electron. J. Geotech. Eng. 15, 949–956 (2010).
Google Scholar
Zoheir, B. & Emam, A. Field and ASTER imagery data for the setting of gold mineralization in Western Allaqi-Heiani belt, Egypt: A case study from the Haimur deposit. J. Afr. Earth Sci. 99, 150–164 (2014).
Article ADS CAS Google Scholar
Ammar A. A., Meleik, M. L. and Fouad K. M., Tectonic analysis of a sample area, Central Eastern Desert, Egypt, applying aeroradiometric and aeromagnetic survey data. Bull. Faculty Earth. Sci. King Abdulaziz Univ. Jeddah, KSA, 6, 459–482 (1983).
Clarke, K. G. Optimum second-derivative and downward continuation filters. Geophysics 34, 424–437(1969).
Henderson, R. G., Zietz, I. The computation of second vertical derivatives of geomagnetic fields. Geophysics 14 (1949).
Dobrin, M. B., Savit, C. H. Introduction to Geophysical Prospecting, 4th Edition, McGraw-Hill Book Co., 867 (1990).
Reid, A. B., Allsop, J. M., Granser, H., Millet, A. J. & Somerton, I. W. Magnetic interpretation in three dimensions using Euler deconvolution. Geophysics 55, 80–91 (1990).
Article ADS Google Scholar
Thompson, D. T. EULDPH: A technique for making computerassisted depth estimates from magnetic data. Geophysics 47, 31–37 (1982).
Article ADS Google Scholar
Naudy, H. Automatic determination of depth on aeromagnetic profiles. Geophysics 36, 717–722. https://doi.org/10.1190/1.1440207 (1971).
Article ADS Google Scholar
Bournas, N., Galdeano, A., Hamoudi, M. & Baker, H. Interpretation of the aeromagnetic map of Eastern Hoggar (Algeria) using the Euler deconvolution, analytic signal and local wavenumber methods. J. Afr. Earth Sc. 37, 191–205 (2003).
Article ADS Google Scholar
Whitehead, N., Musselman, C. Montaj gravity/magnetic interpretation: Processing, analysis and visualization system for 3D inversion of potential field data for oasis montaj v6.3, Geosoft Incorporated, 85 Richmond St. W., Toronto, Ontario, M5H 2C9, Canada (2008).
Hall, D. H. Magnetic and tectonic regionalization in Texada Island. British Columbia. Geophysics 29, 4 (1964).
Article Google Scholar
Said, R. Geology of Egypt. Elsevier Publ. Co., Amsterdam and New York, pp. 293–319 (1962).
Neev, D. Tectonic Evaluation of the Middle East and Levantine Basin. Easternmost Mediterranean Geology, printed in USA (1975).
Tealeb, A. Predominant Tectonic Trends in North Egypt. Academy of Scientific Research & Technology, Helwan Institute of Astronomy & Geophysics, Bulletin No. 189 (1979).
Meshref, W. M. Tectonic framework. In The Geology of Egypt (Chapter VII), Ed. R. Said, Balkema, Rotterdam, pp. 113–155 (1990).
Download references
Open access funding provided by The Science, Technology & Innovation Funding Authority (STDF) in cooperation with The Egyptian Knowledge Bank (EKB). The authors declare that no funds, grants, or other support were received during the preparation of this manuscript.
Geology Department, Faculty of Science, Al-Azhar University, Cairo, Egypt
Waheed H. Mohamed & Mahmoud H. Elyaseer
Egyptian Mineral Resources Authority, PO. Box 11517, Abbassiya, Cairo, Egypt
Mohamed Elsadek M. Sabra
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
Waheed H. Mohamed: Analyzed and interpreted the magnetic data, contributed regents, materials, analysis tools or data and wrote the paper. Mahmoud H. Elyaseer and Mohamed Elsadek M. Sabra are analyzed and interpreted the remote sensing data, contributed regents, materials, analysis tools or data and wrote the paper.
Correspondence to Waheed H. Mohamed.
The authors declare no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
Reprints and Permissions
Mohamed, W.H., Elyaseer, M.H. & Sabra, M.E.M. Structural lineament analysis of the Bir El-Qash area, Central Eastern Desert, Egypt, using integrated remote sensing and aeromagnetic data. Sci Rep 13, 21569 (2023). https://doi.org/10.1038/s41598-023-48660-x
Download citation
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41598-023-48660-x
Anyone you share the following link with will be able to read this content:
Sorry, a shareable link is not currently available for this article.
Provided by the Springer Nature SharedIt content-sharing initiative
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.
Advertisement
© 2023 Springer Nature Limited
Sign up for the Nature Briefing newsletter — what matters in science, free to your inbox daily.