Mapping Main Structures and Related Mineralization of the Arabian Shield (Saudi Arabia) Using Sharp Edge Detector of Transformed Gravity Data

Saudi Arabia covers most of the Arabian Peninsula and is characterized by tectonic regimes ranging from Precambrian to Recent. Using gravity data to produce the lateral boundaries of subsurface density bodies, and edge detection of potential field data, a new subsurface structural map was created to decipher the structural framework controls on the distribution of gold deposits in Saudi Arabia. Moreover, we detected the relationships between major structures and mineral accumulations, thereby simultaneously solving the problem of edge detectors over complex tectonic patterns for both deeper and shallower origins. Analytic signal (ASg), theta map (TM), TDX, and softsign function (SF) filters were applied to gravity data of Saudi Arabia. The results unveil low connectivity along the Najd fault system (NFS) with depth, except perhaps for the central zones along each segment. The central zones are the location of significant gold mineralization, i.e., Fawarah, Gariat Avala, Hamdah, and Ghadarah. Moreover, major fault zones parallel to the Red Sea extend northward from the south, and their connectivity increases with depth and controls numerous gold mines, i.e., Jadmah, Wadi Bidah, Mamilah, and Wadi Leif. These fault zones intersect the NFS in the Midyan Terrane at the northern part of the AS, and their conjugation is suggested to be favorable for gold mineralization. The SF maps revealed the boundary between the Arabian Shield and Arabian Shelf, which comprises major shear zones, implying that most known mineralization sites are linked to post-accretionary structures and are not limited to the Najd fault system (NFS).

Sedimentology and Evolution of the Fluvial-Deltaic System: A Modern Depositional Model Analog from the Red Sea Coastal Region, Saudi Arabia

Siliciclastic sediments represent important lithological unit of the Red Sea coastal plain. Their subsurface equivalents are important targets of groundwater aquifer and hydrocarbon reservoirs in the region. The lithofacies of the modern fluvial deltaic system has several distinct geomorphic units and sub-environments such as alluvial, fluvial, delta plain, aeolian, intertidal, coastal sabkha and eustuarine sediments. This study intends to characterize the lithofacies and the depositional environments and to produce an integrated facies model for this modern fluvial-deltaic system. The study might provide a valuable modern analog to several important subsurface Neogene formations that act as important hydrocarbon reservoirs and groundwater aquifers. The study integrates information and data obtained from landsats, maps and detailed field observation and measurements of facies analysis of the fluvial and deltaic along traveses from the Arabian Shield to the Red Sea coast. The lithofacies sediment analysis revealed four main lithofacies associations namely lithofacies A,B,C ad D. Lithoacies Associations A, which represents the oldest unit is dominated by coarse gravel with minor sands facies. While the lithofacies B is dominated byfine gravel and sand lithofacies, occasionally pebbly, vary from horizontal, planar to massive sands with minor laminated to massive silts and mud facies. The lithofacies in A and B show lateral proximal to distal variation as well as characteristic vertical stacking patterns. The Facies Association A and B indicates a change in fluvial depositional styles from gravelly alluvial fans to gravelly sandy fluvial systems. The lithofacies association C represents the recent fluvial system which consists of minor gravel lag deposits associated maily with various sand lithofacies of planner, horizontal and massive sand associated with massive and limainted sand and mud lithofacies. The lithofacies Association D is dominated with Barchan sand dunes local interfigger with muddy iinterdunes and sand sheets. Lithofacies D occupies rather more distal geomporphic position of the fluvial deltaic system that is adjace to coastal sabkha. The lithofacies associations described here document the evolution and development of the coastal plain sediments through space and time under various autocyclic and allocyclic controls. This included the tectonics and structural development associated with the Red Sea rifting and opening since the Oligocene – Miocene time. Others controls include the evolution of the Arabian shield (provenance) and the coastal plain through space and time as controlled by tectonics, sediment supply, climate and locally by autocyclic environmental This study might be beneficial for understanding the controls and stratigraphic evolution of the Red Sea region and will be of great value for reservoir and aquifer characterization, development and management. This modern analog model can also help in providing geological baseline information that would be beneficial for understanding similar ancient fluvial deltaic sediments. The study might provide guides and leads to understand the subsurface facies, stratigraphic architecture and heterogeneity of any potential groundwater aquifers and hydrocarbon reservoirs.

Regional In-Situ Stress Prediction in Frontier Exploration and Development Areas: Insights from the First-Ever 3D Geomechanical Model of the Arabian Plate

In this paper, we present results from the first-ever 3D geomechanical model that supports pre-drill prediction of regional in-situ stresses throughout the Arabian Plate. The results can be used in various applications in the petroleum industry such as fault slip-tendency analysis, hydraulic fracture stimulation design, wellbore stability analysis and underground carbon storage. The Arabian tectonic plate originated by rifting of NE Africa to form the Red Sea and the Gulfs of Aden and Aqaba. The continental rifting was followed by the formation of collisional zones with eastern Turkey, Eurasia and the Indo-Australian Plate, which resulted in the formation of the Eastern Anatolian fault system, the fold-thrust belts of Zagros and Makran, and the Owen fracture zone. This present-day plate tectonic framework, and the ongoing movement of the Arabian continental lithosphere, exert a first-order control on the of in-situ stresses within its sedimentary basins. Using data from published studies, we developed a 3D finite element of the Arabian lithospheric plate that takes into account interaction between the complex 3D plate geometry and present-day plate boundary velocities, on elastic stress accumulation in the Arabian crust. The model geometry captures the first-order topographic features of the Arabian plate such as the Arabian shield, the Zagros Mountains and sedimentary thickness variations throughout the tectonic plate. The model results provide useful insights into the variations in in-situ stresses in sediments and crystalline basement throughout Arabia. The interaction between forces from different plate boundaries results in a complex transitional stress state (thrust/strike-slip or normal/strike-slip) in the interior regions of the plate such that the regional tectonic stress regime at any point may not be reconciled directly with the anticipated Andersonian stress regimes at the closest plate boundary. In the sedimentary basin east of the Arabian shield, the azimuths of the maximum principal compressive stresses change from ENE in southeast to ~N-S in northern portions of the plate. The shape of the plate boundary, particularly along the collisional boundaries, plays a prominent in controlling both the magnitude and orientations of the principal stresses. In addition, the geometry of the Arabian shield in western KSA and variations in the sedimentary basin thickness, cause significant local stress perturbations over 10 – 100 km length scales in different regions of the plate. The model results can provide quantitative constraints on relative magnitudes of principal stresses and horizontal stress anisotropy, both of which are critical inputs for various subsurface applications such as mechanical earth model (MEM) and subsequently wellbore stability analysis (WSA). The calibrated model results can potentially reduce uncertainties in input stress parameters for MEM and WSA and offer improvements over traditional in-situ stress estimation techniques.

Nb–Ta Behaviour during Magma-to-Pegmatite Transformation Process: Record from Zircon Megacrysts in Pegmatite

Due to the absence of early magma records in pegmatites, it is difficult to investigate the behavior of Nb and Ta during the transformation from magma to pegmatite melt. Zircon megacrysts in an NYF-type (Nb-Y-HREE-F) pegmatite from the Arabian Shield could be divided into three phases from core to margin. The Phase Ι zircon in the core of the zircon megacrysts had typical magma oscillatory zonation with ∑REE content from 300 to 400 ppm, Th/U ratios of less than 0.1 and Nb/Ta ratios of less than 1.0. Phase ΙΙ zircon had oscillatory zonation and was enriched with LREEs mostly with Th/U ratios of 0.1–0.2 and Nb/Ta ratios of 1.0–3.0. Phase ΙΙΙ unzoned zircon had the highest ∑REE content, from 8000 to 15,000 ppm, with Th/U ratios higher than 3.0 and Nb/Ta ratios higher than 5.0. The Hf-O isotopic composition was similar in the different phases of zircon with initial 176Hf/177Hf ratios of 0.28258–0.28277, εHf(t) values from 8.0 to 12.0 and δ18OVSMOW from +4.0‰ to +5.0‰. Zircon megacrysts in the NYF-type pegmatite from the Arabian Shield record the transformation from magma to pegmatite melt. Similar Hf-O isotopic compositions mean a closed magmatic system without contamination by external melt, rock or fluid. The proposed modeling shows that magma with low Nb and Ta concentrations and Nb/Ta ratios could evolve into residual pegmatite melt with a high Nb content and superchondrite Nb/Ta ratio during several stages of melt extraction and fractional crystallization of Ti-rich minerals, such as rutile and titanite. The Nb/Ta ratio can be used as an effective indicator of the transformation process from magma to pegmatite melt.