Electronic Theses and Dissertations (PhDs)

Permanent URI for this collection

Browse

Recent Submissions

Now showing 1 - 7 of 7
  • Item
    Trans-crustal and temporal perspectives of Palaeoproterozoic porphyry copper deposit formation (Haib deposit, Namibia)
    (University of the Witwatersrand, Johannesburg, 2023) Jogee, Bibi Ayesha
    Porphyryِcopperِdepositsِareِtheِworld’sِprimaryِsourceِofِcopperِwhichِisِanِessentialِcriticalِ metal in the green energy sector. These deposits are most often temporally found within the Phanerozoic Eon and occur principally at convergent plate boundaries. The Haib porphyry copper deposit is a rare example of a Palaeoproterozoic porphyry copper deposit. The Haib porphyry copper deposit is hosted by the rock types of the Richtersveld Subprovince. The Richtersveld Subprovince consists of an older volcano-sedimentary sequence, the Orange River Group (ORG), which was intruded by an extensive composite cogenetic suite of plutonic rocks, the Vioolsdrif Intrusive Suite (VIS). The Vioolsdrif Intrusive Suite is subdivided into several subunits, known as the Vuurdood Subsuite (mafic-ultramafic), Gaobis Diorite, Khoromus porphyry, Goodhouse Subsuite (Granodiorite) and the Ramansdrif Alkali Granite. These subsuites are variably composed of diorites, granodiorites, monzogranites, quartz- feldspar porphyry granodiorites and leucogranites as well as poorly-described clinopyroxenites and gabbros (named the Kokerboom Intrusion in this study). The age-range for the Richtersveld Subprovince is between approximately 1910 and 1865 Ma. The mineralization occurs as disseminations, blebs and veins in the feldspar porphyry andesites and quartz-feldspar porphyry granodiorites. The age of mineralization is bracketed at 1884ِ ±ِ 2ِMa. The Haib porphyry copper deposit is in close vicinity to clinopyroxene and amphibole-rich cumulates, which host novel sulphide mineralization, and present lines of evidence suggesting these are the deeper parts of the porphyry copper system. The mineralization is characterized into two main phases based on drill-core sample descriptions, petrographic observations, mineral liberation maps and isotopic analyses. Previous studies show that garnet and/or amphibole fractionation is related to the petrogenesis of porphyry copper deposits. The Haib porphyry copper deposit shows evidence for amphibole fractionation as opposed to garnet fractionation, which is supported by petrological and geochemical observations in this study. The Haib porphyry copper deposit also has limited ranges in δ65Cuِ andِ δ 56Fe signatures supporting exposure at deeper crustal levels of an arc system. Typical cumulates from deep crustal arc sections (arclogites) are coarse-grained, show cumulate textures and are composed of garnet-clinopyroxene-amphibole-iron-titanium oxides and in some cases plagioclase. The mafic-ultramafic cumulates such as the clinopyroxenites observed in the Kokerboom Intrusion, in close proximity to the Haib porphyry copper deposit, are coarse-grained, are composed of clinopyroxene-amphibole-iron-titanium oxides and consist of non-cotectic proportion of mineral phases. They also comprise a high modal abundance of pyroxene (>70 %) and show distinct cumulate textures. This study has demonstrated that the clinopyroxene and amphibole-rich cumulates which belong to the ultramafic end member of the Kokerboom Intrusion form part of the Vuurdood Subsuite. This conclusion was made through the analysis of lithological characteristics, mineral and alteration assemblages, major and trace element geochemistry and U-Pb age data. Various lines of evidence from this study indicate that the Haib porphyry copper deposit represents the deeper levels of the porphyry system. These lines of evidence include the lack of shallow alteration assemblages (propylitic and argillic alteration), the presence of early dark micaceous veins, the coarse-grained nature of the Vioolsdrif Intrusive Suite (clinopyroxenites, granodiorites, diorites, gabbros), which is indicative of slow cooling at mid-crustal batholith depths, and the presence of clinopyroxene and amphibole rich cumulates of the Kokerboom Intrusion. The clinopyroxenites are mafic, dense, and usually occur at mid- to lower-crustal levels. Geothermobarometry conducted on rocks of the Vioolsdrif Intrusive Suite yield anapproximate depth of crystallization between ~4-12 km which forms part of the mid-crustal levels of an arc system.
  • Item
    Cost-effective and novel seismic methods for mineral and coal exploration: Examples from Witwatersrand goldfields and Bushveld Complex
    (University of the Witwatersrand, Johannesburg, 2023) Sihoyiya, Mpofana; Manz, Musa S. D.
    The reflection seismic method has its origin from oil and gas exploration in the land and offshore sedimentary soft rock environment. Since the inception of the reflection seismic technology in hardrock environment, more advanced processing methods such as migration algorithms have been established to revamp the quality of hardrock seismic data for deep mineral targeting. This study shows the value of recovering and reprocessing legacy reflection seismic data using advanced processing techniques that were not available at the time of acquisition. This is achieved through different novel processing workflows that incorporate iterative static corrections, Kirchhoff pre-stack depth migration (KPreSDM), Kirchhoff pre-stack time migration (KPreSTM), and the newly developed Fresnel-volume (FV) and coherency migration (CM) techniques. This research also shows the value of novel processing of the seismic data acquired in-mine challenging environments. The legacy data from the Kaapvaal Craton in South Africa have been processed to improve the delineation of the deep-seated mineral deposits such as the gold-bearing horizons (termed reefs) in the Witwatersrand Basin and Platinum Group Element (PGE) horizons in the Bushveld Complex. Furthermore, the shallow coal seams and associated geological structures of the Karoo Supergroup in the Evander Basin, an arcuate basin characterized by the Witwatersrand Basin towards the south and the Bushveld Complex towards the north, have been delineated for future mine planning and designs. To better image and understand the geometry of the gold-bearing reefs and crosscutting geological structures in the South Rand goldfield, an 18 km long legacy two-dimensional (2D) reflection seismic line was reprocessed using today’s standard processing workflow which employed iterative static corrections. Improved structural imaging of the steeply dipping faults (mostly normal and reverse), as well as dolerite intrusions, was mostly achieved by using KPreSTM and KPreSDM, with KPreSDM providing better structural imaging than other techniques. Moreover, reprocessing of this profile using an improved velocity model and numerical simulations assisted in delineating near surface stratigraphic units and deep-seated (> 1 km depth) geological structures that are associated with the gold-bearing reef but missed by the legacy post-stack time migrated section. In the Bushveld Complex, the PGE deposits (known as platinum reefs) are delineated using a newly developed depth migration technique called coherency migration. Imaging through CM workflow shows evidence of the complex structural architecture that controls the platinum deposits in the study area. The 50 km long legacy 2D profile traverses towards the center of the complex and provides a better understanding of the tectonic evolution in the region. Additionally, magnetic data were utilised to constrain seismic interpretation and delineate the highly magnetic major geological structures such as the Chaneng structure and strongly magnetized rocks of the Rustenburg Layered Suite. The Chaneng structure, in particular, was poorly defined by the seismic data due to its steeply dipping nature. In 2020, in-mine seismic experiments were conducted to delineate the PGE-bearing horizons such as the Upper Group-2 (UG2) and Merensky Reef. Three reflection seismic profiles from these experiments were cautiously processed to attenuate the infrastructure-generated noise. The experiments were successful and delineated the Merensky Reef and UG2 mineralisations at depths between 55 m and 124 m beneath the developmental tunnel (~550 m below the surface). Six legacy 2D reflection seismic profiles acquired for gold exploration in 1986 in the Evander Basin are processed for imaging both the Witwatersrand Basin and the Bushveld Complex, as well as the coal seams of the Highveld and Witbank coalfields. The profiles were acquired perpendicular to each other providing pseudo-three-dimensional (3D) imaging of the subsurface and better delineation of the geological structures. The degree of faulting and folding that affect the mineralisation observed on the processed profiles in all sites would not have been achieved if advanced processing workflows were not implemented. Processing of the legacy data provides a more cost‐effective way to explore the mineral deposits than reacquiring new data, which could be costly and limited by surface conditions. The processing approaches used in this thesis can be used in other brownfield mining regions where the legacy data exist, and they can also be used to revamp the quality of the data acquired in noisy mining environment
  • Item
    Utilizing legacy seismics and non-seismic geophysical methods for deep mineral targeting and near-surface characterization: implications for mine development planning
    (University of the Witwatersrand, Johannesburg, 2023) Mutshafa, Ndamulelo; Manzi, Musa
    This thesis demonstrates through several case studies how the reprocessing of legacy reflection seismic data using advanced algorithms can be of value to mineral exploration and mine development, especially in hard-rock environments. The thesis also showcases how the integration of seismic and non-seismic geophysical datasets can assist in delineating near- surface geological structures (e.g., boulders and fractures) for mine planning and designs. Papers I and II demonstrate how the reprocessed legacy reflection seismic data have been used to delineate and image the world-class gold deposits such as the Black Reef (BLR) and Ventersdorp Contact Reef (VCR) of the Witwatersrand goldfields in South Africa. The two legacy reflection seismic profiles (Paper I) were acquired in 1988 for deep mineral exploration and mine planning over an area that is dominated by dolomitic outcrops that cause scattering of seismic energy at the near-surface, preventing energy propagation into the subsurface. Various migration approaches, namely, pre-stack time, pre-stack depth, and post-stack time migration were applied to test their capabilities in improving structural imaging. Reprocessing results from the pre-stack depth migration using the Kirchhoff algorithm provided the most improved subsurface images, especially the deeper targets due to its ability to honour complex lateral variations in the velocity field. In addition, Kirchhoff's pre-stack and post-stack time migration techniques improved the continuity of the near-surface reflections below the dolomitic rocks. Paper II presents the results from the recovered and processed 25.3 km long legacy seismic survey that was acquired in 1983 by the Gold Division of Anglo-American as part of the Witwatersrand goldfields exploration program. The reprocessing of the data improved the imaging of the gold-bearing horizon termed Ventersdorp Contact Reef (VCR), which is situated at depths between ~2400 and ~4100 m below the ground surface near the South Deep mine in Fochville, South Africa. The pre-stack time and phase-shift migration approaches were tested during processing, and both revealed a dipping reflection associated with the gold-bearing horizon and major steeply dipping faults that crosscut and displace the deposit. The interpretation of the results was constrained using borehole logs and surface geology. This is encouraging and motivates the use of legacy seismic data in the exploration of deep-seated targets. Papers III, IV and V present the results from the use of multi-geophysical methods (resistivity, magnetics, seismic, ground penetrating radar and multichannel analysis of surface waves) at Tharisa platinum mine to provide a comprehensive understanding of the subsurface geology by accurately delineating and locating boulders, mapping fractures and groundwater aquifers to improve Platinum Group Elements (PGEs) mining efficiency and reduce risks. The results from these integrated geophysical methods were successful in complementing each other in terms of providing a clear picture of the near-surface geological structures to help the mine plan better for future operations.
  • Item
    Constraints on the genesis of orbicular granites and sulphide mineralization in the Koperberg Suite, South Africa and the Diana’s Pool area, Zimbabwe
    (University of the Witwatersrand, Johannesburg, 2024) Dumisa, Senamile Siyaya
    This work tackles a long-standing problem in petrology, the formation of orbicular granitoids. These bodies occur as entire facies or as areas enriched in orbicules within distinct facies of plutons. Some European examples are highly prized for their decorative qualities. However, their origin is a matter of much debate going back to times when the origins of granitoids themselves were not obvious and are confused by topics such as ‘granitization’ and over-exaggeration of metasomatic effects. Here, poorly studied outcrops of orbicular rocks from the Matopos granite batholith in the Diana’s Pool area, Zimbabwe and the Koperberg Suite, South Africa are tackeld. Enzman (1953) and Garvie (1969; 1971) studied the Koperberg Suite and the Diana’s Pool orbicules, respectively and employed field observations and petrography to examine and characterize the genesis of these rocks. However, there is little evidence presented to confirm their theories on how these rocks formed as previous work lacks detailed geochemical, mineral chemistry and isotope data to support Enzman and Gravie’s conclusions on the genesis of these rocks. Furthermore, there is no connection between the metallogenesis of the sulphide mineralization and the origin of orbicular rocks, which is a unique characteristic of particular orbicule sites in the Koperberg Suite. This study focuses on four different orbicular bodies (Orbicule Koppie, Henderson South, Henderson North and Hoogkraal Lease) from the Koperberg Suite, hosted in lithologies ranging from diorite to pyroxenite compositions. The orbicules from different orbicular share similar characteristics. Coarse- grained felsic cores, fine-grained and alternating ferromagnesian and feldspathic shells, and coarse- grained to pegmatitic matrices characterize them. The orbicules are generally spherical to ellipsoidal in shape, however, some appear to be abraded and deformed (e.g., Orbicule Koppie). The compositions and grain sizes of cores and the matrix are comparable in all localities. Both the matrix and the cores are medium- to coarse-grained and dominated by plagioclase (the matrix probably in slightly lesser proportions), microcline (in variable proportions, and seemingly absent in some cores), quartz, biotite, magnetite and orthopyroxene (in the case of Hoogkraal Lese, Henderson North and Henderson South). Contrary to the cores and matrix, shells are fine-grained and exhibit polygonal textures. In addition to this, the shells are dominated by biotite and orthopyroxene. Plagioclase in the Henderson North and Henderson South orbicules is more calcic than plagioclase in Hoogkraal Lease and Orbicule Koppie orbicules. Biotite in the Henderson North and the Orbicule Koppie orbicules is more magnesian than those at Henderson South and Hoogkraal Lease orbicules. Plagioclase in the Henderson North orbicules are characterized by elevated and radiogenic initial 87Sr/86Sr ratios than the other orbicule localities. There is also a variation in 87Sr/86Sr ratios within the individual orbicular structures (cores, shells and matrices) in individual localities. The orbicules at Henderson South, Orbicule Koppie and the Jubilee Pit host a bornite-dominated sulphide assemblage where disseminated chalcopyrite and bornite grains are characterized by granular textures and chalcopyrite is replaced by magnetite and bornite. Bismuth, Ni and Se in these sulphides appear to partition into bornite while Ag, In, Cd, Sn, Mn, Ge and Co partition into chalcopyrite. Diana’s Pool orbicular samples exhibit closely packed orbicules in a granitic matrix characterized by different types of orbicules containing coarse-grained felsic cores, fine-grained and alternating ferromagnesian and feldspathic shells, and a coarse-grained to pegmatitic matrix. The orbicules are generally spherical to ellipsoidal in shape, however, some appear to be abraded and deformed. The compositions and grain sizes of cores and the matrix are comparable. Both the matrix and the cores are medium- to coarse-grained and dominated by plagioclase (the matrix probably in slightly lesser proportions), microcline (in variable proportions, and seemingly absent in some cores), quartz, biotite and accessory hornblende and magnetite. Contrary to the cores and matrix, shells are fine-grained and exhibit polygonal textures. In addition to this, the shells are dominated by biotite and magnetite; however, they do not contain hornblende. Plagioclase in cores, shells and matrices shows an almost complete overlap of An contents. Biotite composition in the shells is significantly less magnesian than in core and matrix, whose compositions overlap. Initial 87Sr/86Sr ratios from plagioclase in shells are slightly more radiogenic than in the matrix and cores. Formation of both the Koperberg Suite and Diana’s Pool orbicules has been attributed to metasomatic processes. However, modelling of the Koperberg Suite orbicules using the Magma Chamber Simulator and a variety of textural and geochemical constraints rules out a metasomatic origin. The quartz + biotite-dominated diorite (Orbicule Koppie) and the pyroxene-diorites (Hoogkraal Lease, Henderson South, Henderson North), together with more felsic anorthosite and the more mafic pyroxenites, are all produced by progressive AFC + recharge processes, where the magma was most likely anatectic melts of the country rocks that were at granulite grade. In both Diana’s Pool and the Koperberg Suite, cores are autoliths, which are plagioclase-rich, cumulate, or rim fragments reworked by new magma inputs or injections. Heterogeneous nucleation leading to the formation of orbicular shells around the cores is attributed to adiabatic decompression of magma pulses ascending in dykes leading to superheating and resorption of early solids, and volatile exsolution, inducing undercooling, supersaturation, and shell crystallization. An alternative process that triggered superheating is magma mixing (e.g., Henderson South and Henderson North). The coarse-grained matrix crystallized later, after the orbicules formed, creating the groundmass, and locking the orbicules in place. The deformation of shells and cores (e.g., Orbicule Koppie and Diana’s Pool) suggests that the orbicules continued to evolve in the presence of a melt (matrix material) until they were emplaced at their present setting. The in situ sulphide assemblages in the Koperberg Suite have been reassessed in view of recent research on sulphide behavior. The bornite-chalcopyrite assemblage is atypical of the intermediate solid solution (iss) assemblage (chalcopyrite and pyrrhotite) observed in most Cu-Ni magmatic sulphide deposits. The high concentrations of trace elements that are incompatible with monosulphide solid solution and the depletion of Ni and Co in the sulphides are consistent with the derivation of sulphides from a Cu-rich sulphide melt that separated from a Ni-rich sulphide melt prior to magma emplacement and orbicule formation. The oxidation and Cu-enrichment of a sulphide melt that formed the mineralization in the orbicules and Koperberg Suite more generally began prior to the emplacement of the silicate Koperberg Suite magmas that formed the orbicules, suggesting no link between orbicule genesis and sulphide metallogenesis.
  • Item
    Innovative surface, tunnel, and in-pit geophysical methods for mineral exploration and mine planning: case studies from the Bushveld Complex mines, South Africa
    (University of the Witwatersrand, Johannesburg, 2023) Rapetsoa, Moyagabo Kenneth
    Innovative geophysical methods were used to study the platinum group element mineralisation and their associated geological structures at Maseve and Tharisa mines, western Bushveld Complex. Four case studies are presented in this thesis that incorporate the use of in-mine or near mine geophysical methods for mineral exploration. The first one being in-mine seismic data acquired in 2020 at Maseve mine using cost-effective seismic source and sensors, followed by innovative seismic experiments acquired in 2022 at Maseve mine to evaluate the viability of using tunnel and surface experiments for mineral exploration in a noisy, logistically difficult mine environment. Thirdly, the 2021 integrated geophysical surveys conducted at Tharisa mine to image fractures that act as water pathways into the pit. Finally, integrated geophysical techniques are used to delineate boulders to enhance future mine planning and designs at Tharisa mine. The acquired geophysical data were processed using modern processing algorithms to enhance the target mineralization and complex geological structures in all the sites. In-mine reflection seismic datasets acquired in 2020 at Maseve mine proved useful as they provided optimum imaging of the economic Platinum Group Elements (PGEs) such as the Merensky Reef and Upper Group 2 chromitite layers (known as reefs). This is one of the few in-mine seismic experiments to have been conducted in South Africa for mineral exploration. In 2022, 2D reflection seismic profiles were acquired on surface above the Merensky Reef and Upper Group 2 chromitite, together with four 2D reflection seismic profiles acquired along the mine tunnel at ~ 550 m below the surface and tens of meters above known mineralisation: Merensky Reef and Upper Group 2 chromitite layer. Interpretation of the in-mine and surface seismic data were complemented by the use of 3D ray tracing numerical simulations to understand the distribution and out-of-plane reflectivity from the target mineralization. The 2022 Maseve reflection seismic data improved the imaging of geological structures and mineral deposits. The geophysical data acquired in 2021 at Tharisa mine demonstrated the importance of using near-surface integrated geophysical methods (magnetics, seismics, and electrical resistivity) with other datasets such as borehole logs and physical property measurements to understand the geophysical response of the mineral deposits. Ground magnetic data delineated a major dyke that was identified on the aeromagnetic data and geological mapping. Electrical resistivity tomography, on the other hand, identified linear low resistivity zones that differentiateiii fractured and undisturbed hard rock. Seismic methods were important for depth to bedrock imaging. Integration of geophysical methods was encouraged by the need to understand geological structures (e.g., faults, dykes, iron-rich ultramafic pegmatites, boulders) that can have impact on the efficiency, safety and costs of mining in South Africa. Moreover, this approach encourages the implementation of innovative geophysical surveys in brownfield sites for better mine design and planning, and to increase a life of mine (LoM)
  • Item
    3D seismic constraints on the strato-structural evolution of the deep-water Orange Basin, South Africa
    (University of the Witwatersrand, Johannesburg, 2023) Maduna, Nombuso Gladys; Jinna, Zubair; Manzi, Musa
    This research utilizes seismic attributes and advanced machine learning methodologies to analyse high-resolution 3D reflection seismic data from the deep-water Orange Basin, located offshore western South Africa. The primary goal is to gain valuable insights into the basin's tectonic setting, depositional environment, and hydrocarbon potential. Significant features are delineated within the basin including (1) a gravitational collapse system in the Mesozoic Late Cretaceous, (2) mass flow features in the Cenozoic, (3) natural gas and fluid escape structures, (4) a large slope-perpendicular submarine canyon cutting Oligocene strata, and (5) multiple slope-parallel, sinusoidal channel features in the Miocene. The Late Cretaceous succession exhibits a gravitational collapse system with a translational and compressional domain detaching on seaward-dipping Turonian shales. Gravitational collapse during margin uplift formed fold-and-thrust belts along the slope characterizing the compressional domain. As they are commonly linked to hydrocarbons, the compressional domain of these systems has been extensively studied, while the translational domain has been poorly constrained due to its structural complexity. In this research, the translational domain is shown to contain a mixture of extensional tectonics (normal faults) up-dip and compressional tectonics (thrusts) down-dip, with extensive oblique-slip faults cutting thrusts perpendicularly during the translation of sediment. Variance and chaos, conventional seismic attributes, were used to manually pick and interpret the >500 regional-scale faults arising from the gravitational collapse system. Fault-net, a convolutional neural network (CNN), was compared with these edge-enhancing seismic attributes for extracting faults from the seismic volume. The CNN offers several notable advantages over conventional seismic attributes, such as automation, accelerated analysis, and improved time-efficiency on large datasets. Analysing the distribution, type, and geometry of faults within the basin gave valuable insights into the potential hydrocarbon system at work. Numerous natural gas and fluid escape features are identified in the seismic volume including an elongated mud volcano, pockmarked surfaces, and polygonal faults. The stability of the evolving margin is influenced by the underlying structure of a Late Cretaceous gravitational collapse system, also referred to as a deep-water fold and thrust belt (DWFTB) system. The fault framework within provides primary migration pathways for hydrocarbons. Major seafloor slumping occurs directly above a syncline of the Late Cretaceous DWFTB system. This slumping surrounds a structurally controlled, 4.2 km long elongated mud volcano situated between the translational and compressional domains of the underlying DWFTB system. The late Campanian has the largest accumulation of hydrocarbons evidenced by (1) an anticline with a positive high amplitude anomaly situated at the intersection of the two domains, and (2) >950 pockmarks preserved on the palaeo-surface compared to the 85 pockmarks observed on the seafloor. In addition to tectonics, the onset of stratified oceanographic circulation patterns and climate played a large role in changing depositional trends since the mid-Cenozoic. The Oligocene is characterized by a ~2.3 km wide, >13 km long, slope-perpendicular canyon formed at ~30 Ma during a major sea-level fall by a turbidity current. The Miocene is characterized by a ~14 km wide zone of slope-parallel, sinusoidal channels between water depths of 1 200–1 500 m. The formation and preservation of these features during the Miocene are attributed to the erosive interaction between two distinct water currents: (1) the Antarctic Intermediate Water flowing northwards, and (2) the deep North Atlantic Deep Water bottom currents flowing southwards; and the effects of the Benguela Upwelling System and a dry climate prevailing in southwest Africa all intensifying around 11 Ma. While pre-Miocene hydrocarbons originate from Turonian and Aptian source rocks, the origin of hydrocarbons on the seafloor is likely biogenic, arising from organic-rich sediment in the Miocene
  • Item
    Application of integrated methods to assess and characterise the hydrogeology of coastal aquifers in parts of Lagos, Southwest, Nigeria
    (2020) Yusuf, Mumeen Adebayo
    The scope of this Thesis was to apply integrated methods to characterise the groundwater systems of the Lagos Coastal Basin. Like every coastal area in the world, saline intrusion has been the major challenge threatening the fresh groundwater aquifers of the study area over the last couple of decades, and thus, necessitating its assessment. Environmental isotopes, being a reliable and standard tool in hydrological investigation, was employed in combination with geophysical and hydrogeochemical methods to study the coastal aquifer systems. Geophysical probing of the subsurface revealed an alternating sequence of clay and sand, constituting the major lithological units in the study area. The basin aquifers are hosted essentially by sands and clayey sand, while the modes of aquifer occurrences are unconfined to semi-confined and confined for shallow and deep aquifers, respectively. Hydrochemical interpretation identified a surficial thin layer of fresh groundwater overlying the main zone of saline intrusion, which essentially comprises Ca-HCO3 and Ca-Mg-HCO3, Ca-Mg-HCO3 and Ca-Mg-Cl-SO4 hydrochemical facies for both dry and wet seasons, whereas the surface waters are characterised by Mg-Cl and Na-Cl water types for the lagoon and the ocean, respectively. The evaluation of the chemical processes revealed the dominance of carbonate weathering in the shallow aquifer. Hydrochemical, statistical and geochemical model analyses identified that the groundwater chemistry is significantly controlled by water–rock interaction and ion exchange processes as well as anthropogenic activities. Stable isotopes revealed precipitation as the main source of recharge into the basin aquifer systems. Analyses of the 3 H and 14 C activities were in agreement, revealing an interesting fact about the increase in the groundwater residence time from the surface through deeper depths deducible from 3H values range between 0.1 TU and 2.8TU; 0.0 TU and 0.3 TU; and 14 C age range from 4350±10 to 1050±10 years and between 12030±69 and 7400±50 years for the shallow and deep aquifers, respectively. The mean residence time was supported by the aquifer systems‟ recharge which took place in Holocene for the shallow aquifer and Late Pleistocene–early Holocene for the deep aquifers evident from the calculated ambient temperature, 18 O and 14 C plots. The hydrogeological conceptual models showed that saline incursion severely impacted the second aquifer from a depth ≥20m to 170m in the western and central parts of the study. However, the observed local saline occurrence in places <20m was attributed to groundwater overexploitation. Conclusively, the hydrological systems of the Lagos coastal basin is continually being modified by both anthropogenic and natural activities that constitute not only a major threat to the groundwater sustainability of the Lagos coastal basin but can also consume the entire study area