Thermal Perturbations beneath the Incipient Okavango Rift Zone, Northwest Botswana

We used aeromagnetic and gravity data to investigate the thermal structure beneath the incipient Okavango Rift Zone (ORZ) in northwestern Botswana in order to understand its role in strain localization during rift initiation. We used three-dimensional (3-D) inversion of aeromagnetic data to estimate the Curie Point Depth (CPD) and heat flow under the rift and surrounding basement. We also used two-dimensional (2-D) power-density spectrum analysis of gravity data to estimate the Moho depth. Our results reveal shallow CPD values (8-15 km) and high heat flow (60-90 mW m-2) beneath a ∼60 km wide NE-trending zone coincident with major rift-related border faults and the boundary between Proterozoic orogenic belts. This is accompanied by thin crust ( \u3c 30 km) in the northeastern and southwestern parts of the ORZ. Within the Precambrian basement areas, the CPD values are deeper (16-30 km) and the heat flow estimates are lower (30-50 mW m-2), corresponding to thicker crust (∼40-50 km). We interpret the thermal structure under the ORZ as due to upward migration of hot mantle fluids through the lithospheric column that utilized the presence of Precambrian lithospheric shear zones as conduits. These fluids weaken the crust, enhancing rift nucleation. Our interpretation is supported by 2-D forward modeling of gravity data suggesting the presence of a wedge of altered lithospheric mantle centered beneath the ORZ. If our interpretation is correct, it may result in a potential paradigm shift in which strain localization at continental rift initiation could be achieved through fluid-assisted lithospheric weakening without asthenospheric involvement

Thermal Perturbations beneath the Incipient Okavango Rift Zone, Northwest Botswana

We used aeromagnetic and gravity data to investigate the thermal structure beneath the incipient Okavango Rift Zone (ORZ) in northwestern Botswana in order to understand its role in strain localization during rift initiation. We used three-dimensional (3-D) inversion of aeromagnetic data to estimate the Curie Point Depth (CPD) and heat flow under the rift and surrounding basement. We also used two-dimensional (2-D) power-density spectrum analysis of gravity data to estimate the Moho depth. Our results reveal shallow CPD values (8-15 km) and high heat flow (60-90 mW m-2) beneath a ∼60 km wide NE-trending zone coincident with major rift-related border faults and the boundary between Proterozoic orogenic belts. This is accompanied by thin crust ( \u3c 30 km) in the northeastern and southwestern parts of the ORZ. Within the Precambrian basement areas, the CPD values are deeper (16-30 km) and the heat flow estimates are lower (30-50 mW m-2), corresponding to thicker crust (∼40-50 km). We interpret the thermal structure under the ORZ as due to upward migration of hot mantle fluids through the lithospheric column that utilized the presence of Precambrian lithospheric shear zones as conduits. These fluids weaken the crust, enhancing rift nucleation. Our interpretation is supported by 2-D forward modeling of gravity data suggesting the presence of a wedge of altered lithospheric mantle centered beneath the ORZ. If our interpretation is correct, it may result in a potential paradigm shift in which strain localization at continental rift initiation could be achieved through fluid-assisted lithospheric weakening without asthenospheric involvement

Late Pleistocene hydrological settings at world heritage Tsodilo Hills (NW Kalahari, Botswana), a site of ancient human occupation

Based on genetic studies, a Pleistocene Kalahari “palaeo-wetland”, which spanned the region of the Okavango Delta and the Makgadikgadi Basin, was recently considered the geographic origin of evolutionary modern humans. It was proposed that subsequent out-of-homeland migration was induced by climate shifts. The Tsodilo Hills, which are in relative proximity to the Okavango Delta, represent a site of ancient human occupation since at least 100 ka. Local hydrological dynamics were predominately controlled by climate variability and are archived in the sediments of Palaeolake Tsodilo. This study seeks to better understand the Late Pleistocene environments of the ancient Tsodilo people with a focus on palaeo-hydrological settings, which played a major role for their livelihoods. Our multidisciplinary approach included different remote sensing and geophysical methods, comprehensive application of differential GPS, and sedimentological analyses concentrating on the lake beds. Four palaeo-shorelines could be identified, three of which indicate highstands during which the Tamacha palaeo-river drained Palaeolake Tsodilo towards the Okavango Panhandle. Two highstands during MIS 3b and LGM are related to periods of largely increased fish consumption by humans as has been documented by archaeologists. The palaeolake was likely most extended about 100 ka ago or earlier, when it covered ca. 70 km2 and was 16 m deep. A single (neo-)tectonic fault could be detected. We assume that the Tamacha palaeo-river was a gateway for ancient humans to reach the Tsodilo Hills from these palaeo-wetlands. The people took advantage of the Tsodilo Hills as shelter from weather hazards and as a natural fortress against predators and elephants. Geologically, the Tsodilo Hills were comparatively calm. They represented a relatively safe haven where the social behaviour of early modern humans could evolve to a higher complexity, which relates to the fundamental question when and where modern human behaviour began

Spectral Angle Mapping and AI Methods Applied in Automatic Identification of Placer Deposit Magnetite Using Multispectral Camera Mounted on UAV

The use of drones in mining environments is one way in which data pertaining to the state of a site in various industries can be remotely collected. This paper proposes a combined system that employs a 6-bands multispectral image capturing camera mounted on an Unmanned Aerial Vehicle (UAV) drone, Spectral Angle Mapping (SAM), as well as Artificial Intelligence (AI). Depth possessing multispectral data were captured at different flight elevations. This was in an attempt to find the best elevation where remote identification of magnetite iron sands via the UAV drone specialized in collecting spectral information at a minimum accuracy of +/- 16 nm was possible. Data were analyzed via SAM to deduce the cosine similarity thresholds at each elevation. Using these thresholds, AI algorithms specialized in classifying imagery data were trained and tested to find the best performing model at classifying magnetite iron sand. Considering the post flight logs, the spatial area coverage of 338 m(2), a global classification accuracy of 99.7%, as well the per-class precision of 99.4%, the 20 m flight elevation outputs presented the best performance ratios overall. Thus, the positive outputs of this study suggest viability in a variety of mining and mineral engineering practices

Paleolimnological features of a mega-lake phase in the Makgadikgadi Basin (Kalahari, Botswana) during Marine Isotope Stage 5 inferred from diatoms

The Makgadikgadi–Okavango–Zambezi basin (MOZB) is a structural depression in the south-western branch of the East African Rift System of the northern and middle Kalahari, central southern Africa. In the present day, the mainly dry subbasins of the MOZB are part of a long-lived lacustrine system that has likely existed since Early Pleistocene and from which an extant freshwater fish radiation emerged seeding all major river systems of southern Africa. During hydrologically favourable periods the subbasins were connected as a single mega-lake termed Lake Palaeo-Makgadikgadi. Previous geomorphological studies and OSL dates have provided evidence for repeated mega-lake periods since approximately 300 ka. The environmental and climatic implications of such large scale late Quaternary lake-level fluctuations are controversial, with the duration of mega-lake phases poorly constrained. Here, we present the first evidence for a Marine Isotope Stage (MIS) 5 mega-lake period (about 935–940 m a.s.l.) reconstructed from a diatom-rich, 30-cm-thick lacustrine sediment section, exposed close to a palaeo-shoreline of the Makgadikgadi Basin. Based upon the environmental setting and in comparison with sedimentation rates of other similar lake environments, we tentatively estimated that the highstand lasted approximately 1 ka during MIS 5d–b. The 30-cm section was sampled in 0.5-cm steps. Diatom species diversity ranges from 19 to 30 through the section. The dominant species are Pseudostaurosira brevistriata, Rhopalodia gibberula, Cyclotella meneghiniana and Epithemia sorex. The total of 60 sediment samples provide us with a record at decadal to bi-decadal resolution. Based on diatom assemblages and their oxygen isotope composition (δ18O) we infer an alkaline and mostly oligohaline lake with shallow water conditions prevailing in MIS 5, and is potentially analogous to a Heinrich event. The climate over southern Africa during MIS 5 has been considered very arid but the hydromorphological context of our sediment section indicates that we captured a mega-lake period providing evidence that short-term excursions to significantly higher humidity existed. A hydrologically more favourable environment during MIS 5 than formerly presumed is in line with the early human occupation of the Kalahari

Partitioning and distribution of silver in sediment-hosted Cu-Ag deposits: Evidence from the Ghanzi-Chobe Belt portion of the Kalahari Copper Belt

The Khoemacau Mining (KCM) Cu-Ag deposits in northwestern Botswana are comprised of several mineralised zones, including from south to north: the Banana Zone (made up of the Chalcocite Zone, New Discovery, North Limb, South Limb Definition, and the North East Fold), Mango, Zone 5, Zone 5 North, Boseto, Zeta NE, and Zone 6. The KCM Cu-Ag deposits are characterised by two styles of mineralisation that include disseminated and structurally controlled styles of mineralisation, both typified by the scarcity of Ag-bearing minerals. Laser ab lation inductively coupled mass spectrometry (LA-ICP-MS) and electron probe microanalysis (EPMA) techniques were employed to (i) track the sources of Ag by assessing the Ag budget of the sulphide phases, (ii) assess the partitioning and distribution of Ag in coexisting sulphide phases, and (iii) outline the mode of occurrence of Ag in the host sulphides. Bornite, chalcocite, wittichenite, tetrahedrite-like mineral, covellite, and galena were observed to be Ag-endowed. Bornite and tetrahedrite-like mineral are the principal hosts of Ag and contain, on average, up to 2090 and 4107 ppm Ag, respectively. But when bornite coexists with chalcocite, Ag preferentially partitions into chalcocite (mean Ag concentrations range from 423 to 1050 ppm) than in bornite (mean Ag concentrations range from 201 to 309 ppm). In the chalcocite-wittichenite assemblage, Ag partitions equally in these phases with mean concentrations of 1890 and 1854 ppm, respectively. Galena from the North East Fold contains a significant amount of Ag (mean 477 ppm) than that one from Zone 5, which is Ag-poor (mean up to 87 ppm). Unlike the spatially associated sulphides, chalcopyrite, however, is a poor host for Ag (mean concentration 1.3-104 ppm). The LA-ICP-MS data demonstrated that Ag is essentially bound in the crystal lattice of bornite, chalcocite, wittichenite, tetrahedrite-like mineral, covellite, and locally galena, as shown by the homogeneous distribution of Ag in the elemental maps, small to medium standard deviations (averaging 2-44% of the mean) and relatively flat time-resolved depth profiles. However, large variations in the Ag concentration between laser spots (averaging 51-70% of the mean) from the same sample, as observed in some bornite II grains from the New Discovery and the North East Fold mineralised zones, may indicate the presence of Ag-rich micro to nanoparticles in bornite II. The initial mineralising fluid is thought to have been saturated in Ag and precipitated Ag-rich phases in the mineralised zones southwest of the KCM Property, including the Chalcocite Zone, New Discovery, North Limb, South Limb Definition, and the North East Fold. Subsequently, it became Ag-undersaturated with time as it precipitated Ag-poor phases in the mineralised zones northeast of the KCM Project (such as Zone 5). The positive correlations (R > 0.60) displayed between Ag and As, Se, Sb, Hg, Te, and Bi in the Cu-(Fe)-sulphides demonstrate that these elements could be used as pathfinder elements in sediment-hosted Cu-Ag deposit exploration

Geothermal Resource Mapping in Northern Botswana Inferred from Three-Dimensional Magnetotelluric Inversion

A set of magnetotelluric (MT) data collected from 28 stations at the Kasane Hot Spring in northern Botswana was used to derive a 3-D model of the electrical resistivity distribution around deep geologic structures that we associated with geothermal resources. The dimensionality analysis results revealed that 3-D modeling should be used to determine the electrical structure. The resistivity models revealed a conductive layer with an average thickness of 200 m, representing the overlying sediments of the Proterozoic volcanic rocks. The thick high-resistivity zone (>100 Ωm) below the conductive layer can be associated with Mesoproterozoic bedrock. The MT measurements in this area show a tube-shaped conductive anomaly that could serve as a fluid pathway feeding the hot spring. A fracture-controlled meteoric fluid circulation presumably determines the existence of the Kasane Hot Spring system

Hydrogeological conceptual model of large and complex sedimentary aquifer systems – Central Kalahari Basin (Botswana)

Successful groundwater resources evaluation and management is nowadays typically undertaken using distributed numerical groundwater flow models. Such models largely rely on hydrogeological conceptual models. The conceptual models summarize hydrogeological knowledge of an area to be modelled and thereby providing a framework for numerical model design. In this study, an efficient data integration method for developing hydrogeological conceptual model of the large and hydrogeologically-complex, Central Kalahari Basin (CKB) aquifer system, was undertaken. In that process, suitability of 3-D geological modelling with RockWorks code in iterative combination with standard GIS (ArcGIS) was tested. As a result, six hydrostratigraphic units were identified, their heads and related flow system interdependencies evaluated and hydraulic properties attached. A characteristic feature of the CKB is a thick unsaturated Kalahari Sand Unit (KSU), that restricts the erratic recharge input to <1 mm yr −1 in the centre to about 5–10 mm yr −1 in the eastern fringe. The analysis of the spatial distribution of topological surfaces of the hydrostratigraphic units and hydraulic heads of the aquifers, allowed to identify three flow systems of the three aquifers, Lebung, Ecca and Ghanzi, all three having similar radially-concentric regional groundwater flow patterns directed towards discharge area of Makgadikgadi Pans. That pattern similarity is likely due to various hydraulic interconnections, direct or through aquitard leakages, and also due to the presence of the overlying unconfined, surficial KSU, hydraulically connected with all the three aquifers, redistributing recharge into them. The proposed 3-D geological modelling with RockWorks, turned to be vital and efficient in developing hydrogeological conceptual model of a large and complex multi-layered aquifer systems. Its strength is in simplicity of operation, in conjunctive, iterative use with other software such as standard GIS and in flexibility to interface with numerical groundwater model. As a result of conceptual modelling, fully 3-d, 6 layer numerical model, with shallow, variably-saturated, unconfined layer is finally recommended as a transition from conceptual into numerical model of the CKB