Downloaded from orbit.dtu.dk on: Oct 04, 2021
Investigating groundwater salinity in the Machile-Zambezi Basin (Zambia) with hydrogeophysical methods
Chongo, Mkhuzo
Publication date: 2015
Document Version Publisher's PDF, also known as Version of record
Link back to DTU Orbit
Citation (APA): Chongo, M. (2015). Investigating groundwater salinity in the Machile-Zambezi Basin (Zambia) with hydrogeophysical methods. Technical University of Denmark, DTU Environment.
General rights Copyright and moral rights for the publications made accessible in the public portal are retained by the authors and/or other copyright owners and it is a condition of accessing publications that users recognise and abide by the legal requirements associated with these rights.
Users may download and print one copy of any publication from the public portal for the purpose of private study or research. You may not further distribute the material or use it for any profit-making activity or commercial gain You may freely distribute the URL identifying the publication in the public portal
If you believe that this document breaches copyright please contact us providing details, and we will remove access to the work immediately and investigate your claim.
Investigating groundwater salinity in the Machile-Zambezi Basin (Zambia) with hydrogeophysical methods
Mkhuzo Chongo
PhD Thesis June 2015
Investigating groundwater salinity in the Machile-Zambezi Basin (Zambia) with hydrogeophysical methods
Mkhuzo Chongo
PhD Thesis June 2015
DTU Environment Department of Environmental Engineering Technical University of Denmark
Mkhuzo Chongo
Investigating groundwater salinity in the Machile-Zambezi Basin (Zambia) with hydrogeophysical methods
PhD Thesis, June 2015
The synopsis part of this thesis is available as a pdf-file for download from the DTU research database ORBIT: http://www.orbit.dtu.dk
Address: DTU Environment Department of Environmental Engineering Technical University of Denmark Miljoevej, building 113 2800 Kgs. Lyngby Denmark
Phone reception: +45 4525 1600 Fax: +45 4593 2850
Homepage: http://www.env.dtu.dk E-mail: [email protected]
Printed by: Vester Kopi June 2015
Cover: Torben Dolin Preface This PhD thesis is a presentation of the work conducted under the supervision of Associate Professor Peter Bauer-Gottwein between 1st December 2011 and 2nd March 2015 at the Technical University of Denmark (DTU) with extensive field work in the Machile-Zambezi Basin in south-western Zambia. The work was funded by the Danish Ministry of Foreign Affairs (DANIDA) under the programme for capacity building for the Zambian Water Sector - Phase II. Co supervisors included Anders Vest Christiansen, PhD (Associate Professor, Department of Geosciences, Aarhus University); Imasiku A. Nyambe, PhD (Professor, Geology Department, University of Zambia); and Flemming Larsen, PhD (Head of Department, Geological Survey of Denmark and Greenland). The first part of the thesis is a synopsis that puts into context the main findings of the PhD project. It is based on the second part which gives a detailed account of the work done in the form of three scientific papers referred to as papers I, II, and III. The titles of the three papers are:
I Chongo, M., Christiansen, A.V., Tembo, A., Banda, K.E., Nyambe, I.A., Larsen, F., and Bauer-Gottwein, P. 2014. Airborne and ground based transient electromagnetic mapping of groundwater salinity in the Machile-Zambezi Basin, south-western Zambia. Accepted for publica- tion in Near Surface Geophysics Journal (in press).
II Chongo, M., Christiansen, A.V., Fiandaca, G., Nyambe, I.A., Larsen, F., and Bauer-Gottwein, P .2015. Mapping localized freshwater anomalies in the brackish Paleo-Lake sediments of the Machile- Zambezi Basin with transient electromagnetic sounding, geoelectrical imaging and induced polarization. Submitted to Journal of Applied Geophysics.
III Chongo, M., Nyambe, I.A., Larsen, F., and Bauer-Gottwein, P, 2015. Coupled hydrogeophysical inversion using a groundwater flow and transport model, and transient electromagnetic and direct current elec- trical resistivity data. Manuscript.
iii
iv Acknowledgements I wish to express my thanks and gratitude to all institutions and individuals that have been instrumental to the successful execution and completion of my PhD studies. It is not possible for me to list each and every one of them due to constraints of space on the paper except for a very select few; nevertheless I am truly grateful to all. I would like to thank the Danish Government for funding this research through its Ministry of Foreign Affairs and to the Zambian Government through the Department of Water Affairs under Phase II of the programme to build capacity in the Zambian Water Sector. I would also like to thank my principle supervisor Peter Bauer Gottwein (PhD) for the coaching, encouragement and guidance throughout this period. I am equally grateful to my co- supervisors Anders Vest Christiansen (PhD), Imasiku Anayawa Nyambe (PhD) and Flemming Larsen (PhD). A note of acknowledgement and thanks is also extended to the entire departmental and PhD administration at DTU-Environment and also to the Water Resources Engineering Section faculty members and colleagues. I would also like to appreciate the significant technical support from the Hydrogeophysics group at the Department of Geosciences, Aarhus University with respect to processing and inversion of various geophysical data. Finally, I would like to thank friends and colleagues who supported me by their companionship during my field work in Sesheke and Kazungula districts in Zambia. This also applies to the management and staff at the University of Zambia School of Mines. Last but not the least; I would like to thank my family for only they know how much water has gone under the bridge and the Lord God Almighty for all his goodness and tender mercies.
v Summary The importance of knowing the current state of groundwater resources cannot be over emphasized in rural areas in the developing countries with limited water resources. As such, innovative geophysical techniques are now part of the norm for quick and effective characterization of groundwater resources worldwide. This thesis presents the application of geo-electrical and electromagnetic methods for the investigation of groundwater salinity in the Machile-Zambezi Basin in south western Zambia, southern central Africa. Aerial and ground based transient electromagnetic measurenments were used to map the spatial distribution of apparent electrical resistivity on a regional scale in order to obtain a regional overview of groundwater salinization based on electrical resistivity correlation. Furthermore, ground based transient electromagnetic soundings and direct current and induced polarization measurements were used to investigate on a local scale, indications of surface water/ groundwater exchange from electrical resistivity anomalies coincident with alluvial fans and flood plains as deduced from the aerial electrical resistivity result. New and innovative geophysical data inversion schemes were also developed and used to gain a better explanation of the data collected. These include a new scheme for the joint inversion of direct current and induced polarization data, and transient electromagnetic data; and a new coupled hydrogeophysical inversion setup to allow for the first time the joint use of direct current and transient electromagnetic data in one optimization.
The result from the regional mapping with transient electromagnetic measurenments showed a spatial distribution of electrical resistivity that indicated block faulting in the Machile-Zambezi Basin. Saline groundwater was found to occur predominantly in the low lying graben areas that are essentially an extension of the Palaeo Lake Makgadikgadi system into south western Zambia. In addition, surface water from the Zambezi River was found to interact with saline aquifers to such an extent that the surficial physical form of alluvial fans and flood plains was visible in the spatial distribution of electrical resistivity from the aerial survey up to a depth of about 40 m.
Interpretation of direct current and induced polarization, and transient electromagnetic data using the new joint inversion scheme revealed a fresh water lens overlying the saline aquifer at Kasaya in Kazungula District, Zambia. The freshwater lens appeared to be in hydraulic contact with the
vi Zambezi River where it was thickest (60 m) and had the highest electrical resistivity values (about 200 Ωm) which steadily declined to about 30 Ωm whereas the thickness reduced to around 22 m at the end of the 6 600 m long transect line measured perpendicular to the Zambezi River towards the North. The distribution of chargeability along the Kasaya transect line was found to be correlated with the distribution of electrical resistivity thus giving a strong indication of the intrusion of fresh surface water into a pre-existing saline aquifer. It is postulated that the intrusion of fresh surface water in to the saline aquifer was driven by evapotranspiration. Finally, the new coupled hydrogeophysical inversion approach resulted in sharp estimates of hydrogeological model parameters. This was for a coupled flow and solute transport model setup for the Kasaya transect under the forcing of evapotranspiration. Performance of the coupled hydrogeophysical inversion was better with the inclusion of direct current data in comparison to the use of transient electromagnetic data alone.
The broader implications of these findings is that groundwater salinization in the Machile Zambezi Basin is now known to be strongly influenced by the tectonism of the Palaeo Lake Makgadikgadi system and is therefore not expected to increase over time. Rather, surface water tends to interact with the saline aquifers in places to create freshwater lenses that are an important source of clean drinking water. Therefore, the findings of this thesis will need to be augmented with data and further research from other geoscience disciplines such as surface water hydrology, geochemistry, petrology and meteorology in order to come up with sustainable water resources management practices that are applicable to arid and semi-arid sedimentary basins in general and the Machile-Zambezi Basin in particular within the broader context of the Kalahari Basin of southern Africa. This will make it possible to project the effects of climate variability and anthropogenic activities on the water resources with a view of disaster preparedness and mitigation for which this part of the world is particularly vulnerable.
vii Dansk sammenfatning Vigtigheden af et kendskab forekomsten af grundvandsressourcens kan ikke overvurderes i landområder i den tredje verden med begrænset vandressour- cer. I forbindelse med beskrivelse forekomsten af vandressourcer, er anven- delsen af nye, hurtige og effektive geofysiske metoder et vigtigt element. I denne PhD afhandling præsenteres anvendelse af geoelektriske og elektro- magnetiske metoder til undersøgelse saltforekomster i grundvandet i Machi- le-Zambezi bassinet i Zambia, i det sydlige Afrika. Den rumlige fordeling af jordlagenes tilsyneladende elektriske modstand er målt med anvendelse af luftbåren og landbaseret transient elektromagnetiske målinger med henblik på kortlægning af forekomsten af salt grundvand. Derudover er landbaseret elek- tromagnetiske sonderinger og elektrisk målinger, sammen med induceret po- larisations, blevet anvendt på lokal skala med det formål at påvise overflade- grundvandsinteraktion. Nye og innovative metoder til inversion af geofysiske data er blevet udviklet, og anvendt til at opnå en bedre tolkning og anvendel- se af indsamlede data. Dette omfatter blandt andet en ny metode for sammen- hængende inversion af elektriske data, induceret polarisations data og transi- ent elektromagnetiske data; samt en ny koblet hydrogeofysisk inversionsme- tode, der for første gang gør det muligt at foretage en samlet anvendelse af elektriske og elektromagnetiske data i en optimeringsrutine. Den regionale kortlægning med den transiente elektromagnetiske metode vi- ser en rumlig fordeling af elektriske modstande, der indikerer blokforkastning i Machile-Zambezi bassinet. Salt grundvand forekommer hovedsagelig i lavt- liggende dele af grabenstrukturen, der er en forlængelse af Palaeo Makgadik- gadi søen ind i den sydlige del af Zambia. Det er blevet påvist, at interaktio- nen af overfladevand fra Zambezi floden med salt grundvand er så omfatten- de, at den kan registreres til 40 meters dybde med luftbårne opmålte elektri- ske modstand i alluvialfane og flodplan sedimenterne. Tolkningen af elektriske målinger, induceret polarisationsmålinger og elek- tromagnetiske data, med anvendelse af den nyudviklede koblerede inversi- onsmetode, viste, at fersk grundvand ligger over salt grundvand ved Kasaya i Kazaungula distriktet, Zambia. Hvor ferskvandslinsen er tykkest (60 m), ser den ud til at være i hydraulisk kontakt med Zambezi floden, og have maksi- male elektriske modstand (200 Ωm), som gradvist aftager til 30 Ωm, hvor tykkelsen af ferskvandslaget er reduceret til 22 m for enden af et 6600 m langt transekt, opmålt fra Zambezi floden og vinkelret mod nord. Fordelingen af sedimenterne kapacitet til elektrisk opladning langs Kasaya transektet kor-
viii relerer med dets fordeling af elektrisk modstand, og viser derved en stærk indikation på intrusion af fersk overfladevand ind i salt grundvand, hvor salt- dannelsen skyldes evapotranspiration. Endelig har den nyudviklede, koblede hydrogeofysiske inversionsmetode re- sulteret i en pålidelig bestemmelse af hydrogeologiske model parametre. Det- te blev påvist i en koblet strømnings- og transportmodel for Kansaya transek- tet, med klimatisk betinget forøgelse af evapotranspirationen. Med anvendel- se af den koblede hydrogeofysiske inversionsrutine, blev resultaterne forbed- ret med indarbejdning af elektriske modstandsdata, sammenlignet med an- vendelsen af udelukkende transient elektromagnetiske data. Den overordnede betydning af disse resultater er den, at forekomsten af salt grundvand i høj grad er påvirket af tidligere tektoniske forhold i området, hvor Palaeo Makgadikgadi søens sedimenter forekommer, er det må derfor forventes, at saltforekomsterne ikke fremover vil blive mere udbredt i områ- det. Tværtimod synes interaktionen med overfaldevand at skabe en linse af ferskvand, der har stor betydning med henblik på etablering af drikkevands- forsyninger. Derfor bør resultaterne fra dette studie uddybes med nye data og yderligere forskning fra andre geovidenskabelig discipliner, såsom hydrologi, geokemi, petrologi og meteorologi for at tilvejebringe generelle, bæredygtige løsninger til god praksis til forvaltning af vandressourcer i aride og semi- aride områder, og især i Machile-Zambezi Bassinet og i Kalakari Bassinet i det sydlige Afrika. Dette vil gøre det muligt at belyse effekterne af klimafor- andringer og andre antropogene påvirkninger af vandressourcen med det for- mål at skabe beredskab og forebyggelse i forhold til sygdomsbekæmpelse for en befolkning, der især i denne del af verden, er sårbar.
ix
x Table of contents Preface ...... iii Acknowledgements ...... v Summary ...... vi Dansk sammenfatning ...... viii Table of contents ...... xi 1 Introduction ...... 1 1.1 Motivation ...... 1 1.2 State of the art ...... 1 1.3 Research objectives ...... 2 1.4 Structure of the thesis ...... 3 2 Background and context ...... 5 2.1 TEM for groundwater applications ...... 5 2.2 DCIP for groundwater applications ...... 6 2.3 Hydrogeophysical inversion ...... 6 2.4 Hydrodynamic modelling ...... 7 3 The study area ...... 9 3.1 Geomorphological setting ...... 9 3.2 Geological setting ...... 10 3.3 Hydrological setting ...... 11 3.4 Socio-economic aspects ...... 11 3.5 Wildlife ...... 12 3.6 Data archive ...... 13 4 Materials and methods ...... 15 4.1 Airborne and ground based transient electromagnetic measurements ...... 15 4.2 Direct current and induced polarization measurements ...... 16 4.3 Processing and inversion of TEM and DCIP data ...... 18 4.3.1 Data processing ...... 18 4.3.2 Data inversion ...... 18 4.4 Coupled hydrogeophysical inversion ...... 19 5 Main findings...... 21 5.1 Regional electrical resistivity variations ...... 21 5.2 Localised freshwater/ groundwater interaction ...... 22 5.3 Freshwater intrusion into saline aquifer ...... 23 5.4 Utility of electrical and electromagnetic methods ...... 24 5.5 Coupled hydrogeophysical inversion with combined DC and TEM data ..... 24 6 Discussion ...... 27 7 Conclusion ...... 31 8 Outlook and Perspectives ...... 33 9 References ...... 35 10 APPENDIX: Inventory of electronic material ...... 43 11 Papers ...... 45
xi xii 1 Introduction The Earth is a very watery planet considering that about 70% of it is covered by water in addition to occurrences in the atmosphere and subsurface. How- ever, with respect to human habitat and environmental concerns, only a very small percentage is available as renewable freshwater reserves in the form of surface water and shallow groundwater Postel et al. (1996). The renewability of these fresh water resources (and hence their availability for continued sus- tained use) is very sensitive to climatic variability and anthropogenic activi- ties (Srinivasan et al., 2012; Cassardo, 2014; Collet et al., 2014). Therefore in order to mitigate against the adverse effects of climate change or human so- cio-economic development on fresh water reserves such as groundwater, it is imperative that the condition of the respective water resources are constantly assessed and the impacts of the various stresses (such as reduced precipita- tion, pollution or abstraction) on the water resources themselves are quanti- fied. This would then form a basis upon which decision makers can make judgements about the best choices that would ensure sustainability of water resource utilization (Kinzelbach et al., 2003). 1.1 Motivation Arguably, nowhere is this more important and critical than in arid and semi- arid regions where much of the hydrological cycle oscillates between ex- tremes of precipitation (or flooding) and drought (Buytaert et al., 2012). As a result the dependence on groundwater is increased in such regions (Kinzelbach et al., 2003) like in the Kalahari Basin of southern Africa (McCarthy and Haddon, 2005; McCarthy, 2013). Thus for countries like An- gola, Botswana, Namibia and Zambia that have got significant territory and populations within the Kalahari Basin, knowledge of the properties and char- acteristics of the Kalahari aquifer system is absolutely crucial. Where availa- ble, such knowledge could then be used for important planning decisions such as the placement of water points in relation to human settlements or the zoning of new settlements based on the quality and quantity of groundwater reserves. 1.2 State of the art Various tools and methods have been developed and are still being developed that enable a relatively rapid and cost effective means of characterization of groundwater resources (Kinzelbach et al., 2003). These include ground based and airborne transient electromagnetics (TEM)(Auken et al., 2003;
1 Christiansen, 2003; Yan et al., 2009; Ezersky et al., 2011; Xue et al., 2012), direct current and induced polarization (DCIP)(Dahlin, 2001; Dahlin et al., 2002; Dahlin and Zhou, 2006; Loke et al., 2013), and hydrodynamic model- ling(Rosbjerg and Madsen, 2005; Bauer et al., 2006a; Bauer et al., 2006b). The TEM methods are applied around the world on both regional and local scales and have been shown to be useful for investigation of regional trends on groundwater quality variations and detailed local scale site characterisa- tion (Melloul and Goldenberg, 1997; Bauer-Gottwein et al., 2010; Chongo et al., 2011; Ezersky et al., 2011; Xue et al., 2012; Herckenrath et al., 2013b; Juanah et al., 2013; Podgorski, 2014). On the other hand, DCIP techniques are more suited for detailed local scale investigations since they can only be effectively deployed on ground based systems. However, recent innovations in data acquisition systems has given rise to DCIP instrumentation that com- bines the benefits of high spatial density, cost effectiveness and the integra- tion of the measurement of both direct current and induced polarization in one field setting (Aristodemou and Thomas-Betts, 2000; Dahlin and Zhou, 2006; Poulsen et al., 2010; Gazoty et al., 2012; Loke et al., 2013). Therefore it is increasingly becoming common to design hydrogeological investigations that incorporate multiple types of geophysical data which can then be incor- porated within the context of a hydrodynamic model (Ferré et al., 2009; Hinnell et al., 2010). This then has resulted in a powerful means of evaluating aquifer system behaviour under various stresses such as evapotranspiration and anthropogenic factors like contamination and water abstraction. 1.3 Research objectives The objectives of the research presented in this thesis were to: Gain a better understanding of the occurrence of saline groundwater in the Machile-Zambezi Basin in south-western Zambia related to propagation of the Palaeo Lake Makgadikgadi system through the Okavango-Linyati fault system into Zambia; Make a detailed local scale evaluation, at Kasaya in south western Zambia, of high electrical resistivity anomalies (> 100 Ωm) overlying a generally low electrical resistivity (< 13 Ωm) subsurface based on local scale DCIP and TEM measurements using a new DCIP-TEM joint inversion scheme. Areas with such high electrical resistivity anomalies were identified based on the results of the regional scale TEM mapping and appeared to be relat- ed to surface water groundwater exchange; and
2 Further develop a coupled hydro geophysical inversion framework for ap- plication to a river-aquifer solute transport problem under evapotranspira- tion to enable calibration with both DC and TEM data for the first time. 1.4 Structure of the thesis As mentioned above this thesis is divided into two parts. The first part is a synopsis that provides a background, context and overview of the research conducted at DTU Environment and extensive field campaigns in Kazungula and Sesheke districts of south western Zambia since December 2011 to date. The synopsis therefore comprises an introduction which briefly outlines the purpose and significance of the research (Chapter 1); a background and con- text which puts the thesis in the context of recent scientific developments and state of the art (Chapter 2); a description of the study area (Chapter 3); and a short overview of the methods used (Chapter 4). The main findings are pre- sented in Chapter 5 whereas the discussion and conclusion are presented in chapters 6 and 7 respectively. The second part of the thesis comprises three papers which are a fulfilment of the three overall research objectives stated above respectively. To the best of our knowledge, this study has mapped for the first time the Palaeo Lake Makgadikgadi sediments in the Machile-Zambezi Basin using regional scale ground based and airborne TEM measurements. It has also demonstrated that TEM is effective for mapping groundwater quality (salinity) variations in a sedimentary setting such as that of the Machile-Zambezi Basin (Paper I). Re- gional scale TEM and local scale DCIP and TEM surveys give evidence of significant groundwater surface water interaction in the Machile-Zambezi Basin. Furthermore, the distribution of electrical resistivity and induced po- larization at a study site at Kasaya in Kazungula District, Zambia indicatives recent intrusion of surface water into a pre-existing saline aquifer (Paper II). The intrusion of fresh surface water under evapotranspiration was evaluated using a new CHI that includes both DCIP and TEM data. The CHI showed the evapotranspiration does seem to be the main climatic excitation driving the surface water/ groundwater exchange given the simplifications and con- sequent limitations in the coupled flow and salinity transport model that was used. The CHI approach was thus shown to have tremendous potential as an effective and cost saving tool for hydrological model validation and geophys- ical data interpretation (Paper III).
3
4 2 Background and context Hydrogeophysics is a relatively young area of research often classified as a sub discipline of Hydrology even though it has its origins in the fields of mineral and petroleum exploration (Rubin and Hubbard, 2006; Auken et al., 2009a). The increasing interest in the use of geophysical methods as part of hydrological investigations probably stems from the need for a better under- standing of how hydrological systems particularly groundwater resources re- spond to stresses such as climatic variability, pollution and long term abstrac- tion. Furthermore, classical methods of groundwater investigation such as pumping tests, slug tests, core sample analysis and geophysical borehole log- ging can only be conducted as point measurements within a well or a bore- hole as and as a consequence will usually not be able to adequately capture lateral or spatial variations of hydrological parameters (Rubin and Hubbard, 2006). This lack of spatial data density and variability can potentially be ad- dressed by the incorporation of appropriate geophysical methods into hydro- logical investigations. Recent research efforts have therefore been directed at finding ways and means of using geophysical methods to provide accurate and reliable esti- mates of hydrological parameters such as hydraulic conductivity and porosity in addition to describing hydrological processes like groundwater flow and solute transport (Christiansen, 2010; Dafflon et al., 2010; Pollock and Cirpka, 2010; Christiansen et al., 2011; Busch et al., 2013; Herckenrath et al., 2013b). This has led to rapid developments in geophysical instrumentation, data handling and inversion. Transient electromagnetics (TEM) and direct current and induced polarization (DCIP), which are the main focus of this thesis, are two of the most commonly used geophysical methods for water resource or hydrogeological applications that have been extensively devel- oped and improved in recent years. 2.1 TEM for groundwater applications TEM systems which are typically deployed as either ground based or airborne have seen the advent of instrumentation that is equally suited for traditional geological exploration as well as for hydrogeophysical applications using a high moment for deeper information and a low moment for shallow infor- mation. Furthermore, older instrumentation that is not equipped with the low moment-high moment concept can equally be used for hydrogeophysical ap- plications provided it is calibrated with a just a few accurate ground based
5 TEM soundings from the survey area (Podgorski et al., 2013). Thus TEM has been successfully used for hydrogeological applications such as regional scale salinity mapping (Guerin et al., 2001; Chongo et al., 2011; Ezersky et al., 2011; Viezzoli et al., 2012)and aquifer characterization (Auken et al., 2003; Auken et al., 2009b; Siemon et al., 2009) and local scale salty water – freshwater interaction studies(Albouy et al., 2001; Bauer-Gottwein et al., 2010; Nenna et al., 2013). 2.2 DCIP for groundwater applications DC and IP have benefited from recent innovations that now make it possible to measure DC and IP in one field setup (hence the term DCIP) using multi- channel instruments that not only greatly reduce the amount of time spent in the field but can at the same time increase the amount of data collected per field setting (Dahlin, 2001; Dahlin and Zhou, 2006; Loke et al., 2013; Singha et al., 2014). Examples of hydrogeological considerations to which DCIP has been applied include groundwater flow and solute transport, groundwater salinization (including salty water freshwater interaction in coastal and inland areas), surface water/ groundwater exchange and groundwater contamination or pollution at landfill or other sites (Aristodemou and Thomas-Betts, 2000; Nassir et al., 2000; Binley et al., 2002; Muller et al., 2005; Bauer et al., 2006c; Muller et al., 2010; Vaudelet et al., 2011; Fernandes et al., 2012; Gazoty et al., 2012; Singha et al., 2014; Van Dam et al., 2014). 2.3 Hydrogeophysical inversion CHI is one of three approaches that can be used to inform groundwater mod- els with geophysical data. The other two are sequential hydrogeophysical inversion (SHI) in which a groundwater model is setup, parameterised and calibrated based on prior and independently inverted geophysical models; and joint hydrogeophysical inversion (JHI) in which geophysical and hydrologi- cal model parameters are inverted at the same time and related through petro physical relations (Herckenrath et al., 2013a). Thus SHI follows the tradi- tional approach with which geophysical data has been used for hydrological purposes whereas CHI and JHI are relatively recent developments that run the geophysical and hydrological inversions simultaneously with exchange of information taking place both ways. Nevertheless they differ in the sense that CHI minimises the difference between hydrological state variables (e.g. so- lute concentration) and instrument responses through a petro physical trans- formation of the state variables and subsequent forward response calculation. JHI however inverts for both hydrological model parameters (e.g. hydraulic
6 conductivity) and geophysical model parameters (e.g. electrical resistivity) at the same time; with the hydrological and geophysical parameters being con- strained by an appropriate petro physical relationship such as the log-linear relationship between hydraulic conductivity and electrical resistivity. 2.4 Hydrodynamic modelling Hydrodynamic modelling as applied to water resources in general and groundwater in particular provides an avenue through which scientific theo- ries about hydrological phenomenon and the interplay between different hy- drological processes can be tested, validated or disproved. Thus in practice hydrodynamic modelling formulates required hydrological processes in a ge- neric mathematical/ numerical construct which can then be applied to specific hydrological conceptualizations and parameterizations (Rosbjerg and Madsen, 2005; Todini, 2007). Examples of standard mathematical/ numerical constructs applicable to groundwater include MODFLOW (Harbaugh, 2005 ; Hanson et al., 2014) for groundwater flow modelling, MT3DMS for reactive and solute transport modelling (Zheng and Wang, 1999) and SEWAT (Langevin et al., 2007) for coupled flow and reactive transport modelling. Thus any given model application is at best a parsimonious simplification of natural systems which explains to some degree observed processes and phe- nomenon. The performance of the model will therefore largely depend on the chosen conceptualization (e.g. layering of the subsurface and stresses in a groundwater model) and parametrizations (e.g. aquifer properties such as hy- draulic conductivity and porosity). However there is often a great deal of un- certainty about some or all aspects of the conceptualization and model pa- rameters. A standard approach of dealing with this uncertainty is through the use of in situ data to calibrate and validate the model (Rosbjerg and Madsen, 2005). Quite often the data is very sparsely distributed (e.g. point measuren- ments of hydraulic head and solute concentration from a limited number of boreholes in a wide geographic area) and because of this, model calibration efforts can often be unsatisfactory. However, geophysical measurenments overcome the limitations on data density since they can be more easily and cost effectively deployed on both local and regional scales. Thus the incorpo- ration of geophysical data into classical and innovative model calibrations should lead to improvements in model parameter estimations and reduced uncertainty (Hinnell et al., 2010).
7
8 3 The study area The Machile-Zambezi Basin straddles Kazungula and Sesheke districts in the Southern and Western provinces of Zambia. It was hydrologically defined for purposes of this research as the catchment area with outlet on the Zambezi River at 17o 50.4’ S and 25o38’E but not exceeding the areal extent defined by latitude 16o - 17o 54’ S and longitude 24o 13’ – 26o 22’ E and within the Zambian territory. Thus it covers an area of about 26, 260 km2 which in- cludes the northern extent of the Palaeo Lake Makgadikgadi Basin (Okavan- go-Linyati depression) and the catchment areas of three major streams called Loanja, Machile and Ngwezi. The location of the Machile Zambezi Basin is shown in Fig. 1.
N
22°30’0"E 24°0’0"E 25°30’0"E 27°0’0"E WE !
! !
! S
! !
!
!
! ! ! 1:5,000,000
!
! ! !
!
! ! ! R.Zambezi
!
! R.Cuando
! !
!
! ! !
!
! ! ! !
! ZAMBIA ! !
16°30’0"S !
!
!
!
! ! ! ! !
!
!
! ! ! ! !
! !
!
! ! ! ! ! !
!
! ANGOLA ! !
!
! !
! !
! ! !
! !
! !
! !
!
! ! ! !
!
! !
! ! ! !
!
! ! !
!
! !
!
! ! !
F.Sesheke! !
! !
! !
! !
! ! ! !
! ! !
! !
! ! ! !
! ! !
! ! !
! !
! ! ! ! ! ! ! !
!
!
! ! ! ! !
! ! ! ! !
!
! !
! ! ! !
! !
! !
! !
! R.Zambezi
R.Okavango !
! ! !
! ! !
! ! !
! ! !
! !
R. ! !
! ! !
! ! R.Chobe!
! !
! ! ! ! !
!
! ! ! ! ! !
! ! !
! ! ! !
! ! ! !
! ! !
! ! ! ! ! !
!
!
! !
! !
! ! ! !
! ! !
! ! ! ! !
! !
!
! ! ! ! !
! ! !
! ! !
! ! ! ! !
! F.Mambova
! ! ! ! !
! !
! ! ! !
!
!
! ! ! ! !
! ! ! !
! ! !
! !
! !
! !
!
!
!
! ! ! NAMIBIA ! !
! !
!
! ! !
! ! !
! ! !
! ! ! ! !
! ! ! ! !
! ! ! ! ! !
! ! !
! ! ! !
! !
! !
! ! ! !
! ! !
! !
! Linyati
! !
! ! ! ! !
! !
!
! ! ! ! ! !
!
! ! ! ! ! ! !
! !
! ! ! ! !
! ! !
! ! ! ! ! ! !
!
! !
! ! ! ! !
! !
! !
! !
! ! ! ! ! ! !
! ! ! ! !!!!! !
18°0’0"S !
F.Chobe !
!
! !
!
! ! ! ! ! !
!
! ! ! !
F.Linyati! ! !
!
! ! !
! !
! !
!
! ! !
!
!
!
!
! !
! !
!
!
! !
!
!
! !
!
! !
! LEGEND
! ! !
!
! !
! !
!
! ! ! Stream !
BOTSWANA ! !
!
19°30’0"S ! ! ! !
! ! !
! !
! Fault
! !
! ! ! !
! ! ! !
! F.Thamalakane ! !
F.Gumare ! Makgadikgadi Basin
!
!
!
!
! !
!
! !
! ! ! ! Machile-Zambezi Basin ! ! ! !
!
!
! !! ! ! !!
!
! !
!
!!! !! Internationl Boundary !!
! !
!
!
!
!
! !
!
!
! !
! !
! !
!
! !
21°0’0"S
! !
! 020 40 80 120 160 ! ! ! Km 22°30’0"E 24°0’0"E 25°30’0"E 27°0’0"E Fig. 1: Location of the Machile-Zambezi in Zambia, southern Africa in relation to the Lake Palaeo Makgadikgadi Basin. 3.1 Geomorphological setting The landforms and features of the Machile-Zambezi Basin are characteristic of a fluvial sedimentary environment with an intermittent - ephemeral river network that is prone to extensive flooding. This river network is entrenched onto a topography that is broadly divided into upland and low land areas. The upland areas or hilly belt with altitude between 950 – 1150 m above
9 mean sea level surround the low lying areas (935 -950 m above mean sea lev- el) in form of an arc that stretches from the north western extent of the basin through the Machile-Kafue Watershed in the north east until the south eastern extent of the Basin. The Zambezi River forms the northwest to southeast boundary of the Basin into which many of the streams drain whereas others terminate into alluvial fans before reaching the Zambezi River. A topographic map of the Machile-Zambezi Basin which also depicts the surface geology is shown in Fig. 2.
N 200,000 250,000 300,000 350,000 400,000 450,000 WE
8,250,000 S ! ! ! ! ! !
! ! ! ! !
! !
! UTM 35S ! ! !
! ! !
! ! ! !
! ! ! ! ! ! ! ! !
! ! ! ! ! !
! ! ! ! !
!
! ! 1:2,000,000
! ! ! ! ! ! !
! ! ! ! !
! !
! ! ! ! ! ! ! !
8,200,000 8,200,000 ! !
! !
! ! ! ! ! !
! !
! ! !
! ! ! ! ! ! ! ! ! ! ! ! ! ! ! !
! !
! !
! ! ! ! ! LEGEND ! ! ! ! ! ! ! !
! ! ! ! ! ! ! !
! Town
! ! ! ! !
! ! !
8,150,000 8,150,000 !
!
! !
! !
! ! Kafue-Machile Watershed !
! ! ! ! !
! !!! ! ! !
! !
! ! !!! Machile-Zambezi Basin
! ! ! ! !
! ! ! !
! ! ! !
!
!
! Surface Geology
! ! ! ! ! ! ! !
! ! !
! Quartenary !
8,100,000 8,100,000
!
!
! ! ! ! !
! ! Cenozoic
! ! !
! ! !
Sesheke!
! !
! !
!
! ! !
!
! !
! !
!
!
! Mesozoic Igneous !
! !
!
! !
!
! !
!
! ! ! ! ! !
! ! ! !
! ! !
! !
! ! Pre Cambrian
!
!
!
!
!
! !
!
!
! ! !
! ! !
Kazungula!
! ! !
8,050,000 8,050,000
! ! Relief (m.a.m.s.l) !
! !
! !
!
! ! !
!
! !
! High : 1300 Low : 850
02040608010
8,000,000 200,000 250,000 300,000 350,000 400,000 450,000 Km
Fig. 2: Topographic map of the Machile-Zambezi Basin also depicting the surface geology. 3.2 Geological setting The low lying areas occupy a position that appears to be a continuity in the form of a topographic depression that extends all the way down and across the Zambezi River then into the Okavango-Linyati depression (McCarthy, 2013). This depression, which is a remnant of the Palaeo Lake Makgadikgadi, is controlled by the Okavango-Linyati fault system whose genesis is the propagation of the East African Rift System through the Tanganyika-Mweru- Upemba-Kabopmbo Gorge axis (Moore et al., 2012) . The tectonic activity responsible for the creation of the Okavango-Linyati Depression is said to be still ongoing and has been responsible for evolution of landforms and drain-
10 age patterns in the entire Kalahari Basin (Moore et al., 2012). The low lying areas are thus underlain by sediments that have in-filled a graben structure often referred to as the Machili Graben or the Machili Basin (Main et al., 2008; Moore et al., 2012; McCarthy, 2013). The stratigraphy of the Machile-Zambezi Basin is a variant of the general stratigraphic setting of the entire Kalahari Basin and comprises a Kalahari Supergroup sand member as the uppermost unit. This in turn is typically un- derlain by a Kalahari Supergoup sandstone member or in its absence any se- quence of stratification ranging from Karoo to Basement Complex. The blan- ket of sand varies in thickness from place to place and is generally thinnest on the outer northwest and southeast fringes of the basin particularly close to the Zambezi River where outcrops of Basalt and Basement Complex rocks respectively are clearly visible. A map depicting the surface geology of the Machile-Zambezi Basin is shown in Fig. 2. 3.3 Hydrological setting As has been mentioned above, the Machile-Zambezi Basin has three main streams of which the Machile and Ngwezi flow directly into the Zambezi River whereas the Loanja terminates in a seasonally flooded alluvial fan ad- jacent to the Zambezi. The upper reaches of the stream network originate in the forested hilly belt and are characterised by narrow flood plains and inter- mittent to perennial flow. On the other hand, the lower reaches traverse the low lying areas as they head towards the Zambezi River and are characterised by very wide flood plains and an ephemeral flow regime. Annual precipitation is around 808 mm whereas the actual and potential evapotranspiration are in the order of 616 and 1718 mm per annum (JICA/MEWD-B, 1995). Groundwater recharge has been estimated at 5 % of annual precipitation or around 40 mm per annum. The groundwater regime is predominantly characterised by unconsolidated sand and sandstone aquifer units that have a moderate yield of 3 – 10 l/s and hydraulic conductivities in the order of 1.05 m/day (Chenov, 1978; JICA/MEWD-D, 1995). 3.4 Socio-economic aspects The Machile-Zambezi Basin is a vast rural area with sparse population whose main centres are the towns of Kazungula and Sesheke. The main economic activity is trading with nearby towns like Livingstone in Zambia, Kasane in Botswana and Katima Mulilo in Namibia. Agricultural activity is mostly for subsistence (although some commercial agriculture does take place in the
11 Kazungla area) and involves staple crop production and animal husbandry. A thriving commercial timber industry used to exist in the area but this has now been reduced to small scale operators due to depletion of the timber stock over time. 3.5 Wildlife There exists a great variety of plant and animal species in the Machile- Zambezi Basin. The vast majority of the animal species are confined to the game management areas which are a buffer zone between the Kafue National Park in the northern fringes of the basin and human settlements in the central and southern areas. The Machile-Zambezi Basin also is home to the Machile and Simungoma Important Bird areas given the unique and diverse bird spe- cies that live there. Furthermore a new wildlife sanctuary has been estab- lished in the Machile Flats (the Machile Graben area) called Simalaha Com- munity Conservancy (SCC). It is an important wildlife sanctuary in the Ka- vango-Zambezi trans-frontier wildlife conservation area that seeks to re- establish wildlife populations and their migratory routes between the Chobe National Park in Botswana and the Kafue National Park in Zambia. SCC also aims at improving the economic outlook of the local people through tourism and wildlife development (Peace-Parks, 2015).
Fig. 3: Sable Antelope grazing in the dry season on the southern fringes of the Kafue Na- tional park interfacing with the northern areas of the Machile-Zambezi Basin.
12 The main vegetation types are deciduous tall canopy trees and flood plain and dambo grass. The major tree species include Mopane (Veenendaal et al., 2008; Smit, 2014) and Baikiaea plurijuga locally known as Mukusi or Zam- bezi Teak (Childes and Walker, 1987; Richer, 2008). The Mukusi is promi- nent in the hilly belt whereas the Mopane thrives in the low lying areas of the Basin that are close to the Zambezi River. The tall flood plain and dambo types of grass are primarily found in flood plains and dambos (local scale seasonally flooded pans) (Fanshawe, 2010 ).
(a) (b)
Fig. 4: (a) Typical vegetative cover during the rainy season in the low lying areas of the Machile-Zambezi Basin in south-western Zambia. (b) Vegetative cover in a typical flood plain in the Machile-Zambezi Basin in south western Zambia during the rainy season. In the background is the tree line denoting the beginning of the forested area shown in (a). 3.6 Data archive During the course of the research reported in this thesis, a large amount of old and new data was gathered and processed. This includes old ground based TEM data from Chongo et al. (2011) and new DCIP and TEM data (aerial and ground based) from this work. Furthermore, a large amount of new data was created in the form of intermediary TEM and DCIP data files required for inversion in addition to the final inverse models. A significant number of C++ and Matlab scripts were also produced or modified in order to support and implement the innovative inversion schemes herein described. Further- more, geospatial data for the study area was created and managed in two geo- graphic information systems. These were: ArcGIS (ESRI, 2015) for general spatial data processing and management; and Aarhus Workbench (HGG, 2014) for processing and inversion of aerial TEM data and spatial visualization of aerial and ground based TEM in- verse models.
13 This archive comprising raw data, processed data, programming scripts and geospatial databases is contained in an electronic appendix as out lined in Chapter 9 (Appendix). It may be accessed by contacting:
DTU Environment Technical University of Denmark Miljøvej, Building 113 2800 Kgs. Lyngby Denmark [email protected]
14 4 Materials and methods 4.1 Airborne and ground based transient electromagnetic measurements The TEM (Nabighian, 1988; Nabighian, 1991; Danielsen et al., 2003; Kirsch, 2006) method was used on both local and regional scale for mapping of elec- trical resistivity variations in the Machile-Zambezi Basin. The regional scale TEM was deployed on both airborne and ground based systems for mapping of the entire basin whereas the local scale TEM was conducted using only ground based TEM equipment on a 6.6 km transect line on the northern bank of the Zambezi River at Kasaya in south western Zambia. The instrumenta- tion that was used for the airborne TEM measurements was the VTEM sys- tem from GEOTECH (GEOTECH, 2011) whereas the instrumentation for the ground based TEM comprised WalkTEM (ABEM(b), 2014) and ProTEM equipment. However, data from the ProTEM equipment is not discussed in this thesis. This is because more data was collected with the WalkTEM in- strument. In addition, at locations where measurements were collected with both instruments, similar responses were obtained with slightly better meas- urements from the WalkTEM. The regional TEM survey was designed with two objectives in mind. These were:
To cover as much area of the Machile-Zambezi Basin as possible given the constraints on time and resources; and To capture as much geological variations as possible. Therefore with the constraint of 1000 line kilometers for the airborne meas- urements, the aerial survey was designed with two sets of four flight paths. One set was oriented northeast to southwest whereas another was oriented southeast to northwest. In order to extend the coverage of tem measurements beyond the 1000 line kilometers of the aerial survey, the regional ground based survey was de- signed to conduct spatially distributed TEM measurements at a spatial density of approximately 25 km2 taking into account the survey objectives mentioned above. The ground based survey was conducted only on the eastern half of the Machile-Zambezi Basin because spatially distributed TEM sounding on the western half of the basin were available from Chongo et al. (2011). The
15 locations of the flight paths and all TEM soundings considered in this thesis are shown in Fig. 5. Further details about the regional TEM mapping includ- ing calibration, processing and inversion of the airborne TEM data are given in Paper I. 200,000 250,000 300,000 350,000 400,000 N
8,250,000 8,250,000 WE
S UTM 35S O le X hi ac 1:2,000,000 .M 8,200,000 R 8,200,000 O O iX z O e b O lo O O O O O u .M RX X OO X X X O O O O X OX O XX X O O O X 8,150,000 lo 8,150,000 O fu X hi O X ic X X R O .S O X R . L X O o a O X n j O LEGEND O a X X X O O O X zi O OO e O e w Xl g Town O O O i N O h X . O c R OO O a X 8,100,000 O O X X 8,100,000 O OXXX O .M X O R Stream O O XXO O X O a O R.Za X Oy O mb XXO a ez OX s X X i a X VTEM Flight Path K . X Sesheke R XX XX WalkTEM 2013 X X 8,050,000 X Kazungula 8,050,000 R.Z X Chongo et. al (2011) ProTEM ambezi O Machile Zambezi Basin
025507510012.5
8,000,000 200,000 250,000 300,000 350,000 400,000 8,000,000 Km Fig. 5: Map showing the location of VTEM flight paths and ground based TEM soundings in the Machile-Zambezi Basin. The ground based TEM soundings coincident with the flight paths were used to calibrate the VTEM signal for amplitude and time shift error. The local scale TEM measurements were collected together with DCIP meas- urements on a transect line at Kasaya in the Machile Zambezi Basin as de- scribed in Section 4.2. Sixty-four TEM soundings were thus collected on the transect line at an approximate spacing of 100 m. 4.2 Direct current and induced polarization measurements DCIP (Binley and Kemna, 2006)measurements were conducted on a 6 600 m transect line across the Simalaha flood plain at Kasaya in the Machile- Zambezi Basin. The Terrameter LS (which is an innovative data acquisition system for self-potential, electrical resistivity and induced polarization) (ABEM(a), 2012) was used for the DCIP measurements in a roll along con- tinuous vertical electrical sounding (CVES) setup (Loke, 1999; Loke et al.,
16 2013) with 5 m electrode spacing. The gradient array (Dahlin and Zhou, 2006) protocol was used in order to take advantage of the multi-channel ca- pabilities of the instrument. This resulted in significant reduction in the amount of time spent in the field but with data quality as good as traditional electrode configurations such as the Wenner and Schlumberger arrays (Dahlin and Zhou, 2004). Location of the Kasaya transect within the Machile-Zambezi Basin is shown in Fig. 6. The site was chosen based on electrical resistivity anomalies from the aerial TEM survey that were indica- tive of surface water groundwater exchange. Another site that showed similar traits is the Loanja alluvial fan. However this is not treated in this thesis.
24°40’0"E 24°50’0"E 25°0’0"E 25°10’0"E N WE
S 1:400,000
ile ch a M . R R . Lo a n ja
17°22’30"S Lipumpu Village X Kasaya School X R a .Zam y be a zi s R.Z a LEGEND a K
m . be zi R X Borehole Simalaha Plain
17°33’0"S Kasaya Transect VTEM flight paths 24°40’0"E 24°50’0"E 25°0’0"E 02.5 5 10 15 20 Km Fig. 6: Location of the Kasaya Transect across the Simalaha Plan in the Machile-Zambezi Basin. Also note the termination of the Loanja River into an alluvial fan which was another area investigated but is not reported for purposes of this thesis. More information about the Kasaya Transect including details about how the data collected was processed in an innovative DCIP-TEM inversion scheme can be found in Paper II.
17 4.3 Processing and inversion of TEM and DCIP data 4.3.1 Data processing The electromagnetic data from the aerial survey was presented in the form of a geo-database by GEOTECH (GEOTECH, 2011). This was imported into the Aarhus Workbench (HGG, 2014) for processing and inversion. The aim of the processing was to remove data affected by coupling and systematic noise and to supress random noise by averaging of the electromagnetic data into soundings. The processing also comprised filtering and averaging of ge- ographic position and altitude data and removal of soundings exceeding a roll or pith angle of 25 degrees. This was followed by visual inspection and edit- ing where required (HGG, 2011). Calibration parameters (Sørensen et al., 2011; Podgorski et al., 2013) for correction of amplitude and time shift errors were added to the data at the processing stage. The calibration parameters (time shift and amplitude shift factor) were derived from calibration of the VTEM signal with corresponding ground based WalkTEM soundings (Paper I). Ground based TEM data was processed for removal of data points affected by noise and coupling using SiTEM-SEMDI (HGG, 2001) and sometimes manually with a text editor after download from the WalkTEM instrument. Unlike the airborne electromagnetic data, the ground based TEM data was processed on an individual sounding basis because the data was collected one sounding at a time at different locations. Finally DCIP data was also imported into the Workbench for automatic processing followed by visual inspection and manual editing (Paper II). 4.3.2 Data inversion Following processing, the airborne electromagnetic data was inverted in the Aarhus Workbench using the laterally constrained inversion (LCI) (Auken et al., 2005) approach on a flight line by flight line basis. This was followed by spatially constrained inversion(Christiansen, 2008; Viezzoli et al., 2009) of the whole aerial dataset in one inversion job. Ground based TEM data col- lected for the regional survey was inverted separately for each site using Aar- husInv (Auken et al., 2014). However the inverse models were imported into the Aarhus Workbench for synchronization with the aerial TEM inversion result. This facilitated the creation of mean horizontal electrical resistivity maps at various depth intervals for the Machile-Zambezi Basin (Paper I).
18 DCIP and TEM data collected for the Kasaya transect was inverted using AarhusInv under various schemes. These included separate LCIs of DCIP and TEM data; mutually and laterally constrained inversion (MCI)(Christiansen et al., 2007) of DC and TEM data; and finally an innovative MCI of DCIP and TEM data. Details about the inversion of DCIP and TEM data on the Ka- saya Transect are given in Paper II. 4.4 Coupled hydrogeophysical inversion Coupled hydro-geophysical inversion was performed using the Gauss- Marquardt-Levenberg parameter estimation approach (Doherty, 2005) for a simple coupled flow and solute transport model using TEM and DC data. The model was setup for a case of freshwater intrusion into a saline aquifer under evapotranspiration with a fixed river head boundary. SEAWAT (Langevin et al., 2007) was used to implement this hydro-dynamic model which was first calibrated using TEM data and then with an innovative combination of TEM and DC data. The PEST (Doherty, 2005) implementation of the Gauss- Marquardt-Levenberg algorithm (Marquardt, 1963) was used for the minimi- zation of square differences between transformations of simulated SEAWAT concentrations and apparent electrical resistivity from TEM and DC meas- urements. Thus Archie’s law (Archie, 1941) was used to transform the simu- lated concentrations into simulated electrical resistivities. The electrical resis- tivities were then sampled at locations corresponding to the location of re- spective geophysical data and used for calculation of forward instrument re- sponses or apparent resistivities using AarhusInv (Auken et al., 2014). The optimization was thus performed on the SEAWAT model parameters (e.g. hydraulic conductivity and porosity) and petro-physical parameters (e.g. Archie’s exponent and cementation factor) but evaluated based on the sum of squared differences between the calculated forward responses and measured geophysical data. Further details about the calibration of the SEAWAT model with TEM and DC data are given in Paper III.
19
20 5 Main findings 5.1 Regional electrical resistivity variations Regional aerial and ground TEM electrical resistivity variations in the Machile-Zambezi Basin reflected the tectonic setting of the Okavango- Linyati fault system and Palaeo Lake Makgadikgadi extension into south- western Zambia. A regional graben-horst system with a northeast to south- west strike was outlined. The horst regions were characterized by resistivities values greater than 100 Ωm on the eastern and western flanges of the basin. However, the graben region was characterized by an area than can be consid- ered to be a full graben in the low lying southern central portions of the basin and an expanse that dips into the full graben from the north east to the south- west. The full graben was found to be characterized by electrical resistivity values less than 13 Ωm whereas the dipping expanse had electrical resistivity values ranging between 10-100 Ωm. The low electrical resistivity values in the graben are attributed to the Palaeo Lake Makgadikgadi sediments. A map depicting the spatial distribution of electrical resistivity values for the deep sub subsurface(80-100 m) is shown in Fig. 7. 200,000 250,000 300,000 350,000 400,000 450,000 N
WE
8,250,000 S
UTM 35S ile ch 1:2,000,000 a i 8,200,000 M z 8,200,000 . e R e b il lo h u c .M a R .M a j R R n lo a .M 8,150,000 u 8,150,000 o if a h L c . c i
R h .S
i l R e LEGEND zi e w g Town 8,100,000 .N 8,100,000 R Sesheke R . Stream R Za .Za zi m mbe be zi VTEM Flight Path Kazungula
8,050,000 R 8,050,000 .Z am 010 20406080 bez i Km
8,000,000 200,000 250,000 300,000 350,000 400,000 450,000 Fig. 7: Mean horizontal electrical resistivity distribution at depth interval 80-100 m for the Machile-Zambezi Basin from aerial and ground TEM measurements.
21 5.2 Localised freshwater/ groundwater interaction The drainage network of the Machile-Zambezi Basin is such that streams originate in the high elevation areas on the north-western to north-eastern peripheries of the basin. As they flow southwardly towards the Zambezi riv- er, they cross the low relief full graben area. At this stage, some of the streams form very wide flood plains whereas others terminate into alluvial fans. At several locations within the basin, the flood plains and alluvial fans make an imprint onto the electrical resistivity distribution in the shallow sub- surface (0-40 m) by anomalously high resistivities (greater than 13.2 Ωm) that are coincident with these surface water features in an otherwise low elec- trical resistivity subsurface. This was indicative of intermediate to local scale surface water groundwater exchange. A satellite image showing the coinci- dence of an alluvial fan and flood plain with anomalously high electrical re- sistivity values in the Machile-Zambezi Basin is shown in Fig. 8.
240,000 250,000 260,000 270,000 280,000 290,000 300,000 310,000 N R .L WE oa nj a S
8,080,000
ile h ac .M R
8,070,000 R.Zambezi a y i a ez s gw a N
R.Z R. K a . mbe
z R (a) i
8,060,000 R .L o an ja
ile ch a .M R R.Zam a bezi y i a ez s gw a N
R.Z R. K a .
mbe z R (b) i 240,000 250,000 260,000 270,000 280,000 290,000 300,000 310,000 05102.5 UTM 35S Km Fig. 8: (a.) Alluvial fan and flood plain features for the Loanja River and the Kasaya River as it enters the Zambezi River respectively. (b.) Overlay of 0-20 m electrical resistivity map depicting the imprint of the alluvial fan and flood plain features onto the shallow sub- surface.
22 5.3 Freshwater intrusion into saline aquifer Local scale evaluation of the anomalous superficial electrical resistivities mentioned in Section 5.2 using DCIP and TEM measurements at Kasaya, de- picted a layer with variable high end electrical resistivity distributions and thickness. This layer exhibited electrical resistivity variations from about 200 Ωm at the beginning of the transect (at the edge of the Zambezi River) to about 30 Ωm at the end of the transect line (total distance 6,600 m). It in turn was overlain by a thin low resistivity layer (less than 12.6 Ωm) and underlain by the low electrical resistivity (3.59 Ωm) subsurface and was 60 m thick at the beginning and 22 m thick at the end. In addition, the spatial distribution of chargeability from joint inversion of DCIP and TEM data was more or less correlated to that of the electrical resistivity. This strong correlation between induced polarization and electrical resistivity is a strong indicator of freshwa- ter intrusion into a pre-existing saline environment. The electrical resistivity and chargeability cross sections for the Kasaya Transect are shown in Fig. 9.
Elevation [m]
Elevation [m]
Fig. 9: (a.) Electrical resistivity cross section and (b.) chargeability cross section from DCIP-TEM joint inversion for DCIP and TEM data collected at Kasaya, Southern Prov- ince, Zambia. The mechanism governing the intrusion of freshwater from the Zambezi Riv- er into the salty aquifer at Kasaya is unknown. However, evapotranspiration seems to be a key process. Evaluation of the evapotranspiration hypothesis conducted using an innovative coupled hydro-geophysical inversion of a SEAWAT model with DC and TEM data added credence to the hypothesis. Nevertheless, processes such as seasonal recharge/ flooding, infiltration and percolation, and through flow could possibly have had an impact on the mod- el outcome and could have helped to obtain a better explanation for the DCIP and TEM data.
23 5.4 Utility of electrical and electromagnetic methods TEM and DCIP are effective methods that can be used for quick characteriza- tion of groundwater resources particularly in sedimentary terrain. Major ben- efits can be drawn from innovations in instrumentation and data processing or inversion. Thus aerial TEM systems originally designed for mineral explora- tion such as the VTEM system can now be used for hydrogeological applica- tions provided they are calibrated with accurate and precise ground based measurements. Similarly, innovative data processing schemes such as DCIP- TEM joint inversion can be used to take maximum advantage of recent and innovative data acquisition systems that enable the simultaneous collection of electrical resistivity and induced polarization measurements. Furthermore, TEM and DC data can now be jointly used to inform or calibrate coupled flow and solute transport models to such an extent that reliable estimates of parameters can be achieved as within sharp confidence intervals. 5.5 Coupled hydrogeophysical inversion with combined DC and TEM data An implementation framework for the calibration of hydrodynamic models (SEAWAT freshwater intrusion model in this case) was successfully devel- oped. Model output in form of solute concentrations were transformed through a series of steps into simulated instrument responses and used in the parameter estimation process. Thus the Gauss-Marquardt-Levenberg algo- rithm implemented in PEST was used to minimise the difference between the simulated instrument responses and the measured TEM and DC data. This implementation approach has demonstrated the benefits of using multiple da- tasets in terms of both type and quantity by the improved sensitivities and reduced uncertainty in parameters and parameter estimates. In addition, the developed framework id flexible enough to incorporate any data type (hydro- logical or geophysical) for as long as it can be related to model variables (e.g. solute concentration and hydraulic head) using either direct or empirical rela- tionships). Qualitative evaluation of the CHI for a simple coupled flow and solute transport problem yielded a simulated electrical resistivity profile from the distribution of solute concentrations that is a close approximation of the elec- trical resistivity distribution obtained from joint inversion of DC and TEM data.
24 Simulated electrical resistivity x-section from SEAWAT model Wm 950 1000 900 100 850 10 800 1 Elevation [m] (a) 0 1000 2000 3000 4000 5000 6000 Profile Coordinate [m]
Resistivity cross section from Joint DC-TEM Inversion Wm 950 1000 900 100 850 10 800 1 Elevation [m] (b) 0 1000 2000 3000 4000 5000 6000 Profile Coordinate [m] Fig. 10: (a.) Simulated electrical resistivity cross section generated from solute concentra- tion output of a SEAWAT excited by evapotranspiration at Kasaya, south-western Zambia. (b.) Electrical resistivity cross section from joint inversion of DC and TEM data. It is envisaged that the main limitations to the CHI is the simplicity of the hydrodynamic model that was used. Improved process descriptions and pa- rameterizations in the model would lead to a closer match between the CHI resistivity cross section and the DC-TEM joint inversion cross section.
25
26 6 Discussion Electrical and electromagnetic methods are an effective means of evaluating groundwater resources particularly with respect to groundwater salinity varia- tions and aquifer delineation in sedimentary basins. This is because of the sharp electrical resistivity contrasts between fresh groundwater and salty groundwater, and between aquifers and the background geology (Nobes, 1996; Albouy et al., 2001). Thus using a combination of both aerial and ground based deployment, regional trends in groundwater salinization can be cost effectively mapped. A typical regional survey would therefore begin with a ground based reconnaissance survey at appropriate sample spacing, followed by a more targeted aerial survey to gain finer resolution of regional features of interest as detected by the preliminary ground survey. In this way a balance can be struck between the need for high data density, reduced un- certainty in model interpretations, spatial coverage and costs associated with aerial surveys (Sapia et al., 2014). Furthermore, ground based measurements are more amenable to deployed for finer scale investigations such as site characterization and they can also be used for quality control checks and cal- ibration of aerial measurements (Podgorski et al., 2013). Thus geophysical measurenments offer the possibilities of determining hydrological parameters for a broad range of scales(Rubin and Hubbard, 2006) particularly when used in conjunction with traditional means of hydrological exploration such as core recovery, pumping tests and groundwater sampling (Danielsen et al., 2003; Auken et al., 2009b; Ruggeri et al., 2014). However there are a number of significant challenges that hinder the use or incorporation of geophysical methods into hydrological investigations. These include terrain/geology specific suitability or appropriateness; requirement for sophisticated interpretation algorithms; general association of geophysical methods with mineral and petroleum industries; very high initial costs of in- strumentation; and measurement of other physical earth properties than the hydrological parameters of interest (Rubin and Hubbard, 2006; Kirsch, 2009). Consequently, DCIP and TEM measurenments are not suitable for use in all geological environments. They are very effective in sedimentary systems be- cause of the relatively high electrical resistivity contrasts between targets of interest and their host materials (e.g. clay lenses embedded in sand or fresh water lenses overlying salty aquifers) (Nassir et al., 2000; Guerin et al., 2001; Bauer et al., 2006c; Gazoty et al., 2012). This is in contrast to igneous or metamorphic geology in which groundwater occurs predominantly in the su- perficial weathered zone and in deeper fractures (Raju and Reddy, 1998)
27 which may not exhibit a sufficiently high electrical resistivity contrast to have a meaningful distinction between aquifers and their host environment. This is in addition to the high electrical resistivity nature of such environ- ments which typically results in very low signal to noise ratios. Furthermore, geophysical datasets that are collected for surveys beyond a certain size (e.g. aerial surveys of more than a few 100 km2 or ground based transects of more than 1 km) often demand excessive computer memory and calculation speeds that cannot be met by ordinary computers. As a result, a significant invest- ment in appropriate computer technology and expertise is often required if DCIP and TEM data are to be interpreted and used in a meaningful way. This comes with the requirement for complicated data handling software and in- version algorithms that are often proprietary as a result of many years of ex- pensive research that go into their production. The fact that the same geo- physical principles and methods as are used for hydrological applications can also be used for mineral and petroleum exploration does not help matters ei- ther. The mineral and petroleum industries are in general more willing and able to pay for geophysical instrumentation, data interpretation algorithms and the innovations thereof for as long as they can be shown to improve their operations and productivity. Therefore the laws of supply and demand dictate that the water and environmental sectors are at the losing end of the geophys- ical supply chain and would thus rather use alternative methods that fit their funding profiles in order to answer the pertinent hydrological questions. Many would go on to further argue that after all, geophysical methods do not measure the properties that are of most interest to hydrologists. For example, DCIP and TEM methods measure the chargeability or the electrical resistivity of the subsurface whereas hydrologists are interested in the aquifer properties such as porosity and hydraulic conductivity. There is no direct relationship between most geophysical parameters and hydrological parameters of interest except for the use of empirical relations such as Archie’s Law and hydrogeo- physical inversion. Notwithstanding, the advantages of incorporating geophysical measurements into hydrological investigations are increasingly being recognised. A case in point is the characterisation of regional and local scale subsurface features in the Machile Zambezi Basin using DCIP and TEM measurenments. Standard aerial and ground based TEM measurenments were used to for the first time visualise the extent of the propagation of the Okavango-Linyati Fault System (OLFS) into the Machile-Zambezi Basin in south-western Zambia. The OLFS is considered to be an important structural control on the extent of groundwa-
28 ter salinization in the Machile-Zambezi Basin as it relates to the spatial extent of the Palaeo Lake Makgadikgadi sediments. Thus a block faulting system was inferred in the Machile-Zambezi Basin whose main feature is a northeast to southwest oriented graben structure which is deepest at its south western extent. Salinized groundwater was found to occur predominantly in this gra- ben structure particularly at the deepest end in the southwest (Paper I). One implication of this is that saline groundwater was created at the time of depo- sition of the Palaeo Lake Makgadikgadi sediments and that as a consequence, further salinization of groundwater is not taking place (Banda et al., 2014). Indeed as indicated from the aerial TEM data, significant surface water\ groundwater exchange based on high electrical resistivity values that were found to be coincident with surface water features appears to be taking place, leading to the creation of freshwater lenses that are an important source of clean drinking water. A combination of DCIP and TEM measurements indi- cated that the high electrical resistivity anomalies were as a result of intrusion of fresh surface water into a pre-existing saline aquifer (Paper II). However a more detailed understanding of the surface water / groundwater exchange in- dicated by the geophysical data is still needed. This can only be achieved by conducting further research that will integrate geochemical, hydrological, petrological and meteorological data around the new knowledge that has been brought out by the geophysics as closely as possible under hydrogeophysical inversion approaches similar to those presented in Paper III. This will of no doubt be useful to water resource managers and planning authorities who have to weigh societal and environmental interests with respect to scarce freshwater resources particularly where decisions have to be made with the aid of groundwater models. Thus better estimates of pertinent hydrological parameters and processes can result from new formulations that incorporate as much geophysical data as is available in addition to data from other geo- science disciplines such as petrology, geochemistry, hydrology and meteorol- ogy. Furthermore, model outcomes will be more reliable and better able to predict system behaviour for as long as the model faithfully represents the natural system under consideration and this requires an integrated multi- disciplinary approach.
29
30 7 Conclusion The main conclusions that have been arrived at with respect to the research reported in this thesis are as follows:
Saline groundwater in the Machile-Zambezi Basin in south-western Zam- bia was found to occur in a northeast to southwest oriented graben struc- ture related to propagation of the Okavango-Linyati fault system stretching from Botswana in the south into south-western Zambia. The graben struc- ture is deepest at its south-western extent where the observed low electrical resistivity values have been attributed to Palaeo Lake Makgadikgadi sedi- ments and associated saline groundwater. This marks the first time that the northern extent of the Paleo Lake Makgadikgadi system (including the as- sociated fault structures) extending into Zambia has been mapped using TEM methods;
Considerable surface water/ groundwater exchange takes place in the Machile-Zambezi Basin as indicated by the coincidence of high electrical resistivity anomalies (> 100 Ωm) in a low electrical resistivity (< 13 Ωm) geological setting with surface water features such as alluvial fans and flood plains. One such feature, the Simalaha Flood Pain, was investigated at a local scale at Kasaya in Kazungula District, Zambia, using a combina- tion of DCIP and TEM measurenments. Evaluation of the measurenments under a new DCIP-TEM inversion scheme revealed a freshwater lens with high electrical resistivity values (30-200 Ωm) overlying a Saline aquifer with electrical resistivity values 3.59 Ωm and below. The correlation of electrical resistivity and chargeability from the DCIP-TEM inversion result indicates intrusion of fresh surface water into a pre-existing saline envi- ronment; and
An innovative coupled hydro geophysical inversion framework was devel- oped to enable the calibration of a coupled flow and transport model with both DC and TEM data for the first time. When tested for a river-aquifer exchange problem under evapotranspiration for the Kasaya Transect in Kazungula District, Zambia, the CHI resulted in sharp parameter confi- dence intervals and reasonable estimates of parameter values. The estimat- ed parameters also exhibited strong sensitivity with respect to model out- comes.
31 Finally, the practical relevance of the findings is that there is now new knowledge on a regional scale about the spatial distribution of saline groundwater in the Machile-Zambezi Basin including local scale surface water/ groundwater exchange in places. This knowledge will be useful in planning for new water points or human settlements taking into account the location of fresh and saline aquifers. Furthermore, these findings have an important impact on the management of the water resources not only of the Machile-Zambezi Basin, but of the entire Kalahari Basin in Southern Afri- ca and other similar semi-arid sedimentary settings around the world. This is because standard transient electromagnetic and geoelectrical (direct cur- rent and induced polarization) methods have been shown to be effective means of groundwater characterization in a semi-arid sedimentary envi- ronment. In addition, the findings have also shown that the data require- ment for groundwater models which are an important water resources management tool and decision support can be effectively complimented by using geophysical data.
32 8 Outlook and Perspectives The hydrogeophysical characterization of groundwater salinization reported in this thesis is an important outcome for water resources management in arid and semi-arid environments particularly in the Machile-Zambezi Basin. However is does not represent the whole picture from the hydrogeological point of view because by and large the geophysical data was not integrated with other data such as geochemical (e.g. prevailing groundwater types and their evolution), petrological (e.g. laboratory determinations of rock and sed- iment properties like porosity and cementation) and hydrogeological (e.g. hydraulic conductivity and hydraulic head) data due to a lack of it. The few such data available in the study area were mostly not coincident with the ge- ophysical data. Therefore there is still a strong need to conduct further geo- chemical, petrophysical, sedimentological, hydrological and hydrogeological research to generate data that can help to strengthen and or validate the out- comes from the hydrogeophysical investigation reported in this thesis. Furthermore, it is highly likely that given the nature of the propagation of the East African Rift Valley system into southern African, there could be other graben structures in the western parts of Zambia similar to the Machile gra- ben that are hidden by the extensive Kalahari sand cover but nevertheless harbour salinized groundwater as a result of palaeo lake systems. One such palaeo lake system that could be investigated in the future is the area around the Barotse Flood plains which has been designated as Palaeo Lake Bulozi by Moore et al. (2012). The Barotse flood plains are an important wetland on the course of the Zambezi River and are protected under the Ramsar Conven- tion on conservation of wetlands of international importance. Thus in order to come up with measures that mitigate against the adverse effects of climate variability such as floods and droughts, and to develop water resource man- agement and development practices that maintain the heritage status of the Barotse Flood plains whilst sustaining the requirements of the local people, the need for ongoing and integrated multi-disciplinary water resource related research cannot be over emphasised. Such research will need to integrate sur- face water and groundwater evaluation using multiple data from remote sens- ing and airborne geophysics to ground based geophysics and sedimentologi- cal, petrological, groundwater and geochemical sampling.
33 A recommended research strategy with respect to the Kalahari System of western Zambia would therefore be to:
Use satellite data for a preliminary description of the hydrology of the Zambezi River system and design of a network of permanent hydrological monitoring stations;
Conduct preliminary ground based TEM measurements in a 25 x 25 km grid in order to gain a regional perspective of ground water salinization based on electrical resistivity correlation. Geochemical sediment sampling could also be incorporated into the preliminary TEM survey in order to have an initial basis for relating geological mineralization to electrical re- sistivity variations;
Design and develop a ground water monitoring network based on the pre- liminary TEM and geochemical sampling survey mentioned above with a rough spatial density of one monitoring borehole per 500 km2. The hydro- geological, petrological and geochemical data derived from these bore- holes could then be integrated with the geophysical data using various geo- statistic and inversion approaches to interpolate and extrapolate parameters of interest. It is recommended that the monitoring boreholes should have a minimum depth of 500 m in order to capture all relevant geological his- tory and not only that of the relatively recent Kalahari Supergroup;
Design and conduct site specific airborne TEM surveys based on areas of interest identified by the preliminary ground based TEM and geochemical sampling surveys mentioned above, followed by detailed local scale sur- veys that integrate different geoscience disciplines in order to gain a com- prehensive picture about the ecosystem and water resources particularly with regards to the impact of climatic variability and anthropogenic activi- ties; and
Create an electronic and online database that integrates local and interna- tional multi-disciplinary data. This database should go together with coor- dinated environmental and water resources related research across the in- ternational boundaries of the southern African countries with significant territory in the Kalahari Basin.
34 9 References ABEM(A) 2012. Terrameter LS User's Guide. Stockholm: ABEM Geophysics.
ABEM(B) 2014. WalkTEM User's Guide. Stockholm: ABEM Geophysics.
ALBOUY, Y., ANDRIEUX, P., RAKOTONDRASOA, G., RITZ, M., DESCLOITRES, M., JOIN, J. L. & RASOLOMANANA, E. 2001. Mapping coastal aquifers by joint inversion of DC and TEM soundings - Three case histories. Ground Water, 39, 87-97.
ARCHIE, G. E. 1941. The Electrical Resistivity Log as an Aid in Determining Some Reservoir Characteristics. SPE Reprint Series, 9-16.
ARISTODEMOU, E. & THOMAS-BETTS, A. 2000. DC resistivity and induced polarisation investigations at a waste disposal site and its environments. Journal of Applied Geophysics, 44, 275-302.
AUKEN, E., CHRISTIANSEN, A. V., JACOBSEN, B. H., FOGED, N. & SORENSEN, K. I. 2005. Piecewise 1D laterally constrained inversion of resistivity data. Geophysical Prospecting, 53, 497-506.
AUKEN, E., CHRISTIANSEN, A. V., KIRKEGAARD, C., FIANDACA, G., SCHAMPER, C., BEHROOZMAND, A. A., BINLEY, A., NIELSEN, E., EFFERSO, F., CHRISTENSEN, N. B. I., SORENSEN, K., FOGED, N. & VIGNOLI, G. 2014. An overview of a highly versatile forward and stable inverse algorithm for airborne, ground-based and borehole electromagnetic and electric data. Exploration Geophysics.
AUKEN, E., GUE´RIN, R., MARSILY, G. D. & SAILHAC, P. 2009a. Hydrogeophysics - Foreword. Comptes Rendus Geoscience.
AUKEN, E., JØRGENSEN, F. & SØRENSEN, K. I. 2003. Large Scale TEM Investigations for Groundwater. Exploration Geophysics, 34, 88–194.
AUKEN, E., SORENSEN, K., LYKKE-ANDERSEN, H., BAKKER, M., BOSCH, A., GUNNINK, J., BINOT, F., GABRIEL, G., GRINAT, M., RUMPEL, H.-M., STEUER, A., WIEDERHOLD, H., WONIK, T., CHRISTENSEN, P.-F., FRIBORG, R., GULDAGER, H., THOMSEN, S., CHRISTENSEN, B., HINSBY, K., JORGENSEN, F., BALLING, I. M., NYEGAARD, P., SEIFERT, D., SORMENBORG, T., CHRISTENSEN, S., KIRSCH, R., SCHEER, W., CHRISTENSEN, J. F., JOHNSEN, R., PEDERSEN, J., KROEGER, J., ZARTH, M., REHLI, H.-J., ROETTGER, B., SIEMON, B., PETERSEN, K., KJAERSTRUP, M., MOSE, K.-M., ERFURT, P., SANDERSEN, P., JOKUMSEN, V., NIELSEN, S. O. & GRP, B. W. 2009b. Buried Quaternary valleys - a geophysical approach. ZEITSCHRIFT DER DEUTSCHEN GESELLSCHAFT FUR GEOWISSENSCHAFTEN, 160, 237-247.
BANDA, K. E., JAKOBSENC, R., GOTTWEIN, P. B.-., MURRAYD, A. S., NYAMBE, I. & LARSEN, F. 2014. The Lake Palaeo-Makgadikgadi in western Zambia: its formation and role in producing recent saline groundwater. In preparation. Geological Survey of Denmark and Greenland.
35 BAUER-GOTTWEIN, P., GONDWE, B. N., CHRISTIANSEN, L., HERCKENRATH, D., KGOTLHANG, L. & ZIMMERMANN, S. 2010. Hydrogeophysical exploration of three-dimensional salinity anomalies with the time-domain electromagnetic method (TDEM). Journal of Hydrology, 380, 318-329.
BAUER, P., GUMBRICHT, T. & KINZELBACH, W. 2006a. A regional coupled surface water/groundwater model of the Okavango Delta, Botswana. Water Resources Research, 42.
BAUER, P., HELD, R. J., ZIMMERMANN, S., LINN, F. & KINZELBACH, W. 2006b. Coupled flow and salinity transport modelling in semi-arid environments: The Shashe River Valley, Botswana. Journal of Hydrology, 316, 163-183.
BAUER, P., SUPPER, R., ZIMMERMANN, S. & KINZELBACH, W. 2006c. Geoelectrical Imaging of Groundwater Salinisation in the Okavango Delta, Botswana. Journal of Applied Geophysics, 60, 126-141.
BINLEY, A., CASSIANI, G., MIDDLETON, R. & WINSHIP, P. 2002. Vadose zone flow model parameterisation using cross-borehole radar and resistivity imaging. Journal of Hydrology, 267, 147-159.
BINLEY, A. & KEMNA, A. 2006. DC Resistivity and Induced Polarization Methods. In: RUBIN, Y. & HUBBARD, S. S. (eds.) Hydrogeophysics. Springer.
BUSCH, S., WEIHERMÜLLER, L., HUISMAN, J. A., STEELMAN, C. M., ENDRES, A. L., VEREECKEN, H. & VAN DER KRUK, J. 2013. Coupled hydrogeophysical inversion of time-lapse surface GPR data to estimate hydraulic properties of a layered subsurface. Water Resources Research, 49, 8480-8494.
BUYTAERT, W., FRIESEN, J., LIEBE, J. & LUDWIG, R. 2012. Assessment and Management of Water Resources in Developing, Semi-arid and Arid Regions. Water Resources Management, 26, 841-844.
CASSARDO, C. 2014. Global warming and water sustainability. Science and the Future, 2.
CHENOV, C. D. 1978. Groundwater Resources Inventory of Zambia. UNESCO/NORAD Water Resources Project. Lusaka: National Council for Scientific Research.
CHILDES, S. L. & WALKER, B. H. 1987. ECOLOGY AND DYNAMICS OF THE WOODY VEGETATION ON THE KALAHARI SANDS IN HWANGE-NATIONAL- PARK, ZIMBABWE. Vegetatio, 72, 111-128.
CHONGO, M., WIBROE, J., STAAL-THOMSEN, K., MOSES, M., NYAMBE, I. A., LARSEN, F. & BAUER-GOTTWEIN, P. 2011. The use of Time Domain Electromagnetic method and Continuous Vertical Electrical Sounding to map groundwater salinity in the Barotse sub-basin, Zambia. Physics and Chemistry of the Earth, 36, 798-805.
CHRISTIANSEN, A. V. 2003. Application of Airborne TEM methods in Denmark and Layered 2D Inversion of Resistivity Data. PhD, University of Aarhus.
36 CHRISTIANSEN, A. V. 2008. A25 Delaunay Triangulation-Based Spatially Constrained Inversion for Quasi 3D Modeling of Airborne TEM Data.
CHRISTIANSEN, A. V., AUKEN, E., FOGED, N. & SORENSEN, K. I. 2007. Mutually and laterally constrained inversion of CVES and TEM data: a case study. Near Surface Geophysics, 5, 115-123.
CHRISTIANSEN, L. 2010. Using ground-based time-lapse gravity observations for hydrological model calibration.
CHRISTIANSEN, L., BINNING, P. J., ROSBJERG, D., ANDERSEN, O. B. & BAUER- GOTTWEIN, P. 2011. Using time-lapse gravity for groundwater model calibration: An application to alluvial aquifer storage. Water Resources Research, 47, W06503.
COLLET, L., RUELLAND, D., BORRELL-ESTUPINA, V. & SERVAT, E. 2014. Assessing the long-term impact of climatic variability and human activities on the water resources of a meso-scale Mediterranean catchment. Hydrological Sciences Journal-Journal Des Sciences Hydrologiques, 59, 1457-1469.
DAFFLON, B., IRVING, J. & HOLLIGER, K. 2010. Calibration of high-resolution geophysical data with tracer test measurements to improve hydrological predictions. Advances in Water Resources, 33, 55-68.
DAHLIN, T. 2001. The development of DC resistivity imaging techniques. Computers and Geosciences, 27, 1019-1029.
DAHLIN, T., LEROUX, V. & NISSEN, J. 2002. Measuring techniques in induced polarisation imaging. Journal of Applied Geophysics, 50, 279-298.
DAHLIN, T. & ZHOU, B. 2004. A numerical comparison of 2D resistivity imaging with 10 electrode arrays. Geophysical Prospecting, 52, 379-398.
DAHLIN, T. & ZHOU, B. 2006. Multiple-gradient array measurements for multichannel 2D resistivity imaging. Near Surface Geophysics, 4, 113-123.
DANIELSEN, J. E., AUKEN, E., JØRGENSEN, F., SØNDERGAARD, V. & SØRENSEN, K. I. 2003. The application of the transient electromagnetic method in hydrogeophysical surveys. Journal of Applied Geophysics, 53, 181-198.
DOHERTY, J. 2005. PEST Model-Independent Parameter Estimation User Manual
5th ed. Australia: Watermark Numerical Computing.
ESRI. 2015. ArcGIS for Desktop [Online]. California, USA: Environmental Systems Research Institue. Available: http://www.esri.com/software/arcgis/arcgis-for-desktop.
EZERSKY, M., LEGCHENKO, A., AL-ZOUBI, A., LEVI, E., AKKAWI, E. & CHALIKAKIS, K. 2011. TEM study of the geoelectrical structure and groundwater salinity of the Nahal Hever sinkhole site, Dead Sea shore, Israel. Journal of Applied Geophysics, 75, 99-112.
37 FANSHAWE, D. B. 2010 VEGETATION DESCRIPTIONS OF THE UPPER ZAMBEZI DISTRICTS OF ZAMBIA. In: TIMBERLAKE, J. R. & BINGHAM, M. G. (eds.) Occasional Publications in Biodiversity No. 22 ed. Bulawayo, Zimbabwe: Biodiversity Foundation for Africa.
FERNANDES, F., CELLIGOI, A., SILVA, S. M. C. P., DALL'ANTÔNIA, L. H., TEIXEIRA, R. S. & LOPES, D. D. 2012. Geophysical technique and groundwater monitoring to detect leachate contamination in the surrounding area of a landfill - Londrina (PR - Brazil). Journal of Environmental Management, 113, 481-487.
FERRÉ, T., BENTLEY, L., BINLEY, A., LINDE, N., KEMNA, A., SINGHA, K., HOLLIGER, K., HUISMAN, J. A. & MINSLEY, B. 2009. Critical Steps for the Continuing Advancement of Hydrogeophysics. Eos, Transactions American Geophysical Union, 90, 200-200.
GAZOTY, A., FIANDACA, G., PEDERSEN, J., AUKEN, E., CHRISTIANSEN, A. V. & PEDERSEN, J. K. 2012. Application of time domain induced polarization to the mapping of lithotypes in a landfill site. Hydrology and Earth System Sciences, 16, 1793-1804.
GEOTECH 2011. Survey and logistcs report on a helicopter borne versatile time domain electromagnetic survey on the Zambezi River Basin Kazungula Zambia for Ministry of Energy and Water Development (Republic of Zambia). West Indies: GEOTECH AIRBORNE LIMITED.
GUERIN, R., DESCLOITRES, M., COUDRAIN, A., TALBI, A. & GALLAIRE, R. 2001. Geophysical surveys for identifying saline groundwater in the semi-arid region of the central Altiplano, Bolivia. Hydrological Processes, 15, 3287-3301.
HANSON, R. T., BOYCE, S. E., SCHMID, W., HUGHES, J. D., MEHL, S. M., LEAKE, S. A., III, T. M. & NISWONGER, R. G. 2014. One-Water Hydrologic Flow Model (MODFLOW-OWHM). Reston, Virginia: U.S. Geological Survey.
HARBAUGH, A. W. 2005 MODFLOW-2005, The U.S. Geological Survey Modular Ground-Water Model—the Ground-Water Flow Process. U.S. Geological Survey Techniques and Methods 6-A16. Reston, Virginia: U.S. Geological Survey.
HERCKENRATH, D., FIANDACA, G., AUKEN, E. & BAUER-GOTTWEIN, P. 2013a. Sequential and joint hydrogeophysical inversion using a field-scale groundwater model with ERT and TDEM data. Hydrol. Earth Syst. Sci., 17, 4043-4060.
HERCKENRATH, D., ODLUM, N., NENNA, V., KNIGHT, R., AUKEN, E. & BAUER- GOTTWEIN, P. 2013b. Calibrating a Salt Water Intrusion Model with Time-Domain Electromagnetic Data. Ground Water, 51, 385-397.
HGG 2001. Getting Started with SiTEM and SEMDI. Aarhus, Denmark University of Aarhus.
HGG 2011. Guideline and standards for Skytem measurements, processing and inversion. Aarhus, Denmark: Department of Earth Sciences, Aarhus University.
38 HGG. 2014. Aarhus Workbench [Online]. Aarhus University: Hydrogeophysics Group. Available: http://hgg.au.dk/software/aarhus-workbench/ [Accessed 30-06-2014 2014].
HINNELL, A. C., FERRE, T. P. A., VRUGT, J. A., HUISMAN, J. A., MOYSEY, S., RINGS, J. & KOWALSKY, M. B. 2010. Improved extraction of hydrologic information from geophysical data through coupled hydrogeophysical inversion. Water Resources Research, 46.
JICA/MEWD-B 1995. National Water Resources Master Plan (the Study of the). Supporting Report [B], Meteorology. Lusaka: Ministry of Energy and Water Development/Japanese International Cooperation Agency.
JICA/MEWD-D 1995. National Water Resources Master Plan (The Study on the). Suppoting Report [D], Hydrogeology. LUSAKA: Ministry of Energy and Water Development/Japanese International Cooperation Agency.
JUANAH, M. S. E., IBRAHIM, S., SULAIMAN, W. N. A. & LATIF, P. A. 2013. Groundwater resources assessment using integrated geophysical techniques in the southwestern region of Peninsular Malaysia. ARABIAN JOURNAL OF GEOSCIENCES, 6, 4129-4144.
KINZELBACH, W., BAUER, P., SIEGFRIED, T. & BRUNNER, P. 2003. Sustainable groundwater management - problems and scientific tools. Episodes, 26, 279-284.
KIRSCH, R. (ed.) 2006. Groundwater Geophysics, A Tool for Hydrogeology, Flintbek, Germany: Springer.
KIRSCH, R. 2009. Groundwater geophysics : a tool for hydrogeology.
LANGEVIN, C. D., DANIEL T. THORNE, J., DAUSMAN, A. M., SUKOP, M. C. & GUO, W. 2007. SEAWAT Version 4: A Computer Program for Simulation of Multi- Species Solute and Heat Transport:Techniques and Methods Book 6, Chapter A22. 39: U.S. Geological Survey.
LOKE, M. H. 1999. Electrical Imaging Surveys for Environmental and Engineering Studies, A Practical Guide to 2-D and 3-D Surveys. Malaysia.
LOKE, M. H., CHAMBERS, J. E., RUCKER, D. F., KURAS, O. & WILKINSON, P. B. 2013. Recent developments in the direct-current geoelectrical imaging method. Journal of Applied Geophysics, 95, 135-156.
MAIN, M. P. L., MOORE, A. E., WILLIAMS, H. B. & COTTERILL, F. P. D. 2008. The Zambezi River. Large Rivers: Geomorphology and Management, 311-332.
MARQUARDT, D. W. 1963. AN ALGORITHM FOR LEAST-SQUARES ESTIMATION OF NONLINEAR PARAMETERS. JOURNAL OF THE SOCIETY FOR INDUSTRIAL AND APPLIED MATHEMATICS, 11, 431-441.
MCCARTHY, T. S. 2013. The okavango delta and its place in the geomorphological evolution of Southern Africa. South African Journal of Geology, 116, 3-54.
39 MCCARTHY, T. S. & HADDON, I. G. 2005. The Mesozoic-Cenozoic Interior Sag Basins of Central Africa: The Late Cretaceous-Cenozoic Kalahari and Okavango Basins. Journal of African Earth Sciences, 43, 316-333.
MELLOUL, A. J. & GOLDENBERG, L. C. 1997. Monitoring of seawater intrusion in coastal aquifers: Basics and local concerns. Journal of Environmental Management, 51, 73-86.
MOORE, A. E., COTTERILL, F. P. D. & ECKARDT, F. D. 2012. THE EVOLUTION AND AGES OF MAKGADIKGADI PALAEO-LAKES: CONSILIENT EVIDENCE FROM KALAHARI DRAINAGE EVOLUTION SOUTH-CENTRAL AFRICA. South African Journal of Geology, 115, 385-413.
MULLER, K., VANDERBORGHT, J., ENGLERT, A., KEMNA, A., HUISMAN, J. A., RINGS, J. & VEREECKEN, H. 2010. Imaging and characterization of solute transport during two tracer tests in a shallow aquifer using electrical resistivity tomography and multilevel groundwater samplers. Water Resources Research, Water Resour. Res, 46.
MULLER, K., VANDERBORGHT, J., ENGLERT, A., KEMNA, A. & VEREECKEN, H. 2005. Characterization of transport processes in a heterogeneous aquifer using electrical resistivity tomography (ERT). In: THOMSON, N. R. (ed.) Bringing Groundwater Quality Research to the Watershed Scale.
NABIGHIAN, M. N. 1988. Electromagnetic Methods in Applied Geophysics Vol. 1.
NABIGHIAN, M. N. 1991. Electromagnetic Methods in Applied Geophysics Vol. 2. Electromagnetic Methods in Applied Geophysics Vol. 2.
NASSIR, S. S. A., LOKE, M. H., LEE, C. Y. & NAWAWI, M. N. M. 2000. Salt-water intrusion mapping by geoelectrical imaging surveys. Geophysical Prospecting, 48, 647-661.
NENNA, V., HERCKENRATH, D., KNIGHT, R., ODLUM, N. & MCPHEE, D. 2013. Application and evaluation of electromagnetic methods for imaging saltwater intrusion in coastal aquifers: Seaside Groundwater Basin, California. Geophysics, 78.
NOBES, D. C. 1996. Troubled waters: Environmental applications of electrical and electromagnetic methods. Surveys in Geophysics, 17, 393-454.
PEACE-PARKS. 2015. SIMALAHA COMMUNITY CONSERVANCY IN THE KAZA TFCA [Online]. Available: http://www.peaceparks.org/programme.php?pid=25&mid=1120.
PODGORSKI, J., AUKEN, E., SCHAMPER, C., VEST CHRISTIANSEN, A., KALSCHEUER, T. & GREEN, A. 2013. Processing and inversion of commercial helicopter time-domain electromagnetic data for environmental assessments and geologic and hydrologic mapping. Geophysics, 78, E149-E159.
PODGORSKI, J. E. 2014. HYDROGEOLOGICAL INVESTIGATION OF THE OKAVANGO DELTA, BOTSWANA USING HELICOPTER TEM AND GROUND-BASED GEOPHYSICAL METHODS. PhD, ETH ZURICH.
40 POLLOCK, D. & CIRPKA, O. A. 2010. Fully coupled hydrogeophysical inversion of synthetic salt tracer experiments. Water Resources Research, 46, W07501.
POSTEL, S. L., DAILY, G. C. & EHRLICH, P. R. 1996. Human appropriation of renewable fresh water. Science, 271, 785-788.
POULSEN, S. E., RASMUSSEN, K. R., CHRISTENSEN, N. B. & CHRISTENSEN, S. 2010. Evaluating the salinity distribution of a shallow coastal aquifer by vertical multielectrode profiling (Denmark). Hydrogeology Journal, 18, 161-171.
RAJU, N. J. & REDDY, T. V. K. 1998. Fracture pattern and electrical resistivity studies for groundwater exploration. Environmental Geology, 34, 175-182.
RICHER, R. A. 2008. Leaf phenology and carbon dynamics in six leguminous trees. African Journal of Ecology, 46, 88-95.
ROSBJERG, D. & MADSEN, H. 2005. Concepts of hydrological modelling. Encyclopedia of Hydrological Sciences.
RUBIN, Y. & HUBBARD, S. S. 2006. Hydrogeophysics, Springer.
RUGGERI, P., GLOAGUEN, E., LEFEBVRE, R., IRVING, J. & HOLLIGER, K. 2014. Integration of hydrological and geophysical data beyond the local scale: Application of Bayesian sequential simulation to field data from the Saint-Lambert-de-Lauzon site, Quebec, Canada. Journal of Hydrology, 514, 271-280.
SAPIA, V., VIEZZOLI, A., JORGENSEN, F., OLDENBORGER, G. A. & MARCHETTI, M. 2014. The Impact on Geological and Hydrogeological Mapping Results of Moving from Ground to Airborne TEM. Journal of Environmental and Engineering Geophysics, 19, 53-66.
SIEMON, B., CHRISTIANSEN, A. V. & AUKEN, E. 2009. A review of helicopter-borne electromagnetic methods for groundwater exploration. Near Surface Geophysics, 7, 629-646.
SINGHA, K., DAY-LEWIS, F. D., JOHNSON, T. & SLATER, L. D. 2014. Advances in interpretation of subsurface processes with time-lapse electrical imaging. Hydrological Processes, n/a-n/a.
SMIT, N. 2014. Response of Colophospermum mopane to different intensities of tree thinning in the Mopane Bushveld of southern Africa. African Journal of Range & Forage Science, 31, 173-177.
SØRENSEN, K., AUKEN, E., CHRISTIANSEN, A. V., FOGED, N., PEDERSEN, J. & SCHAMPER, C. 2011. Refinement of the national TEM reference model at Lyngby. Aarhus: Aarhus University.
SRINIVASAN, V., LAMBIN, E. F., GORELICK, S. M., THOMPSON, B. H. & ROZELLE, S. 2012. The nature and causes of the global water crisis: Syndromes from a meta-analysis of coupled human-water studies. Water Resources Research, 48.
41 TODINI, E. 2007. Hydrological catchment modelling: past, present and future. Hydrology and Earth System Sciences, 11, 468-482.
VAN DAM, R. L., EUSTICE, B. P., HYNDMAN, D. W., WOOD, W. W. & SIMMONS, C. T. 2014. Electrical imaging and fluid modeling of convective fingering in a shallow water-table aquifer. Water Resources Research, 50, 954-968.
VAUDELET, P., SCHMUTZ, M., PESSEL, M., FRANCESCHI, M., GUERIN, R., ATTEIA, O., BLONDEL, A., NGOMSEU, C., GALAUP, S., REJIBA, F. & BEGASSAT, P. 2011. Mapping of contaminant plumes with geoelectrical methods. A case study in urban context. Journal of Applied Geophysics, 75, 738-751.
VEENENDAAL, E. M., MANTLANA, K. B., PAMMENTER, N. W., WEBER, P., HUNTSMAN-MAPILA, P. & LLOYD, J. 2008. Growth form and seasonal variation in leaf gas exchange of Colophospermum mopane savanna trees in northwest Botswana. Tree Physiology, 28, 417-424.
VIEZZOLI, A., AUKEN, E. & MUNDAY, T. 2009. Spatially constrained inversion for quasi 3D modelling of airborne electromagnetic data – an application for environmental assessment in the Lower Murray Region of South Australia. Exploration Geophysics, 40, 173-183.
VIEZZOLI, A., MUNDAY, T. & COOPER, Y. L. 2012. Airborne electromagnetics for groundwater salinity mapping: case studies of coastal and inland salinisation from around the world. BOLLETTINO DI GEOFISICA TEORICA ED APPLICATA, 53, 581- 599.
XUE, G.-Q., BAI, C.-Y., YAN, S., GREENHALGH, S., LI, M.-F. & ZHOU, N.-N. 2012. Deep sounding TEM investigation method based on a modified fixed central-loop system. Journal of Applied Geophysics, 76, 23-32.
YAN, S., CHEN, M. S. & SHI, X. X. 2009. Transient electromagnetic sounding using a 5 m square loop. Exploration Geophysics, 40, 193-196.
ZHENG, C. & WANG, P. P. 1999. A modular three-dimensional multispecies transport model for simulation of advection, dispersion and chemical reactions of contaminants in groundwater systems. Tuscaloosa, Alabama 35487-0338: University of Alabama.
42 10 Appendix Inventory of electronic material Electronic materials associated with the production of this thesis are listed below. They could not be placed in a CD as is usually the practice due to the requirement for large storage space (63.9 gigabytes). The containing folder is called Chongo_Thesis_Electronic_Materials and maybe obtained by contact- ing the following: DTU Environment Technical University of Denmark Miljøvej, Building 113 2800 Kgs. Lyngby Denmark [email protected]
The contents of the electronic appendix are as follows: Raw data items Geotech db/dt aerial data Kasaya transect DCIP data Kasaya transect differential GPS files Kasaya transect TEM data Kasaya CHI setup Machile-Zambezi ground based TEM (raw and processed)
Processed data items Machile-Zambezi Aarhus workbench workspace for aerial and ground based TEM data Kasaya transect Aarhus workbench processed DCIP data files Kasaya transect TEM data (.usf and .tem files) Kasaya transect DCIP data (RES2DINV format and electrode coordi- nates) DCIP_TEM_Joint (Folder containing the DCIP-TEM joint inversion scheme) Joint_DC_TEM_Inversion (Folder containing the joint DC-TEM inver- sion scheme)
43 Databases Sesheke GIS Endnote references
Scripts C++ CHI Scripts Matlab CHI Scripts
Documents Thesis Papers I-III
44 11 Papers
The scientific papers that form the basis of this thesis are:
I Mkhuzo Chongo, Anders Vest Christiansen, Alice Tembo, Kawawa E. Banda, Imasiku A Nyambe, Flemming Larsen, Peter Bauer-Gottwein, 2015. Airborne and ground based transient electromagnetic mapping of groundwater salinity in the Machile-Zambezi Basin, south-western Zam- bia. Accepted for publication in Near Surface Geophysics Journal (in press).
II Mkhuzo Chongo, Anders Vest Christiansen, Gianluca Fiandaca, Imasiku A Nyambe, Flemming Larsen and Peter Bauer-Gottwein, 2015. Mapping localized freshwater anomalies in the brackish Paleo-Lake sediments of the Machile-Zambezi Basin with transient electromagnetic sounding, geoelectrical imaging and induced polarization. Submitted to Journal of Applied Geophysics.
III Mkhuzo Chongo., Nyambe, I.A., Larsen, F., and Peter Bauer-Gottwein, 2015. Coupled hydrogeophysical inversion using a groundwater flow and transport model, and transient electromagnetic and direct current electrical resistivity data.
45
I
Airborne and ground based transient electromag- netic mapping of groundwater salinity in the Machile-Zambezi Basin, south-western Zambia
Chongo, M., Christiansen, A.V., Tembo, A., Banda, K.E., Nyambe, I.A., Larsen, F., Bauer-Gottwein, P.
Accepted for publication in Near Surface Geophysics Journal (Manuscript ID: nsg-2014-1142.R1)
Airborne and ground based transient electromagnetic mapping of groundwater salinity in the Machile-Zambezi Basin, south-western Zambia.
Mkhuzo Chongo1,*, Anders Vest Christiansen2, Alice Tembo3, Kawawa E. Banda1,4, Imasiku A Nyambe4, Flemming Larsen5, Peter Bauer- Gottwein1
1Technical University of Denmark, Department of Environmental Engi- neering, Miljøvej, Building 113, 2800 Kgs. Lyngby, Denmark 2Aarhus University, Department of Geoscience, C.F. Møllers Alle 4, 8000-Århus C, Denmark 3University of Zambia, School of Mines, Integrated Water Resources Management Centre, Great East Road Campus, P.O Box 32379, Lu- saka, Zambia 4University of Zambia, School of Mines, Department of Geology, Great East Road Campus, P.O Box 32379, Lusaka, Zambia 5 Geological Survey of Denmark and Greenland (GEUS), Geochemis- try, Øster Volgade 10, 1350 København K, Denmark
Accepted for publication in Near Surface Geophysics Journal (Manu- script ID: nsg-2014-1142.R1)
*Email: [email protected]; [email protected].
Abstract The geological and morphological evolution of the Kalahari Basin of South- ern Africa has given rise to a complex hydrogeological regime that is affected by water quality issues. Among these concerns is the occurrence of saline groundwater. Airborne and ground based electromagnetic surveying is an ef- ficient tool for mapping groundwater quality variations and has been used extensively to explore the Kalahari sediments e.g. in Botswana and Namibia. Recently, airborne and ground based mapping of groundwater salinity was conducted in the Machile-Zambezi Basin, south-western Zambia using the VTEM and WalkTEM systems respectively, incorporating earlier ground based ProTEM 47D measurements. The data were inverted using the lateral- ly constrained inversion technique followed by a separate spatially con- strained inversion scheme. WalkTEM data was inverted as ordinary single site 1D inversions. The regional electrical resistivity signature of the Machile-Zambezi Basin was found to be characterised by high elevation (1,000-1,050 m amsl), high electrical resistivity (above 100 Ωm) areas that form the western and eastern boundaries of a low resistivity (below 13 Ωm) valley that extends south-westwards into the Makgadikgadi salt pans. The electrical resistivity distribution is indicative of a full graben related to the Okavango-Linyati Fault system as a result of propagation of the East African Rift Valley System into Southern Africa. The saline lacustrine sediments in- filling the Machile Graben are responsible for the low formation resistivity (below13 Ωm) and high salinity (above 7,000 µS/cm) observed in the groundwater and are probably related to the complex evolutionary history of palaeo Lake Makgadikgadi.
I - 1
1 Introduction The Kalahari sediments contain important groundwater resources throughout Southern Africa which however are constrained by a very variable water quality. Much effort has been made to characterize the sediments and groundwater resources of the Kalahari Basin (Worthington, 1977; Arad, 1984; Thomas and Shaw, 1993; Thomas and Shaw, 2002; McCarthy and Haddon, 2005; McCarthy, 2013), a basin which stretches from South Africa in the south into the Democratic Republic of Congo over a distance of 2200 km (McCarthy and Haddon, 2005) to form an extensive sea of sand (Thomas and Shaw, 1990). Geological events related to the breakup of Gondwana have had a major influence on the geological evolution and geomorphology of Southern Africa in general and the Kalahari Basin in particular (McCarthy, 2013). This has been characterised by a series of tectonic, volcanic, fluvial and aeolian processes over geological time scales. At the base of the simpli- fied Kalahari stratigraphy is conglomerate sitting on pre-Kalahari formations. The conglomerate is in turn overlain by sandstone which in turn is overlain by the unconsolidated sands (Money, 1972; McCarthy and Haddon, 2005). The simplified stratigraphy exhibits local variations in terms of the presence or absence of lithological units. For example in places, the conglomerate is absent and sandstone overlain by aeolian sands rests directly on the bedrock whereas in others mudstones and siltstones overly the conglomerate before coarsening upwards into sandstone (McCarthy, 2013). Fig. 1 depicts the sur- face geology (un-consolidated sands) in the Okavango-Makgadikgadi- Machile area as being predominantly Cenozioic. One way of addressing water quality concerns is to use the transient electro- magnetic (TEM) method to identify water quality variations in potential aqui- fers (Melloul and Goldenberg, 1997; Land et al., 2004; Ezersky et al., 2011). In the last few decades TEM has been developed and used extensively around the globe to identify fresh water aquifers, delineate salt water-freshwater in- terfaces and map the spatial extent of salty groundwater (Guerin et al., 2001; Auken et al., 2003; Danielsen et al., 2007; Levi et al., 2008; Boucher et al., 2009; Bauer-Gottwein et al., 2010; Herckenrath, 2012; Nenna et al., 2013). Large-scale airborne geophysical mapping in the Okavango Delta, conducted by the Government of Botswana, has played a major role in the understanding of the groundwater system of the Kalahari Basin (Campbell et al., 2006; Kgotlhang, 2008; Podgorski et al., 2013a).
I - 2
Fig. 1: Contour 945 m above mean sea level (amsl) corresponding to the palaeo Lake Makgadikgadi shoreline in southern Africa. In the background is the topography and Ce- nozoic surface geology corresponding to the extent of the unconsolidated Kalahari sedi- ment. Also depicted are various faults in the Okavango-Machile depression here collective- ly referred to as the Okavango-Linyati Fault System (Persits et al., 2002; Jarvis et al., 2008; Moore et al., 2012; McCarthy, 2013; USGS, 2014) and the outline of the Machile- Zambezi Basin.
Airborne and follow-up ground based TEM data were used in these studies to map potential fresh water aquifers on the basis of their relatively high for- mation electrical resistivities. Such aquifers appeared as sinuous zones of moderate resistivity in an otherwise low resistivity environment. This reflect- ed the wide spread distribution of saline groundwater and to a lesser extent clay rich units inter-bedded with the sands. The TEM data was used by Campbell et al. (2006) to set the lateral extent of potential freshwater aquifers as well as the fresh-saline water interface. Earlier work by Sattel and Kgotlhang (2004) in the Boteti area of Botswana used airborne TEM data collected with the TEMPEST system to map geological features that are con- sistent with the occurrence of fresh groundwater in the Lake Palaeo Makgadikgadi system (Burrough et al., 2009). From this they noted that freshwater occurred as recharge lenses above saline groundwater. Their cor- relation of TEM derived resistivity depth profiles with borehole data showed
I - 3
that the resistivity of the Kalahari Super-group lithology is defined by the level of clay content, amount of saturation and pore water salinity. Using a numerical modelling approach, Bauer et al. (2006) reproduced the salinity distribution in and around the Shashe River Valley in Botswana by modelling transpiration as a function of groundwater salinity. The model results were in agreement with airborne and ground based TEM data (including water level data and geochemical evaluation). The Airborne TEM data clearly outlined fresh groundwater lenses resulting from infiltration of fresh surface water along stream channels into a stagnant saline aquifer in the interfluves. In ad- dition, profiling with ground based TEM across the Shashe River valley sug- gested the thickness of the freshwater lenses to be in the order of 60 m, and to be superficially embedded in a highly saline aquifer (formation electrical re- sistivity between 0.5 – 5 Ωm) with strong formation resistivity correlation to the total dissolved solids of the groundwater given the saturated conditions and relatively uniform sandy lithology (Bauer et al., 2006). In order to define the occurrence of salty groundwater in the Sesheke area in Zambia, Chongo et al. (2011) conducted single site ground based TEM measurements on a re- gional scale and showed that salty groundwater mainly occurs in a topograph- ic depression east of Sesheke Town which is related to faulting associated with the East African Rift Valley System (EARS) (McCarthy, 2013). They further hypothesised that the origin of the salty groundwater was related to evapo-concentration of salts in inter-dune deposits and evaporation of the palaeo Lake Makgadikgadi. The general layered earth model from the TEM measurements indicated an unconfined freshwater aquifer with electrical re- sistivities in excess of 70 Ωm underlain by salty groundwater with resistivi- ties below 35 Ωm. Podgorski et al. (2013a) describe the Okavango Delta as an alluvial mega-fan and used airborne TEM data constrained by ground based TEM, seismic data and borehole records to provide additional evidence for an ancient large scale alluvial fan complex underlying the Makgadikgadi sediments. Their airborne TEM data was a reprocessing of a commercially obtained VTEM dataset from the Okavango Delta affected by systematic er- rors and with large uncertainties regarding key system parameters. This was the same type of VTEM system used for the airborne TEM survey in the Machile-Zambezi Basin for this study and a number of steps were taken to retrieve the best possible results; given that the data collected was prone to contamination of early time gates by transmitter currents, noise in late time gates, and amplitude shifts between adjacent flight lines. Most importantly, accurate ground based TEM data were used to calibrate the VTEM data in order to eliminate absolute amplitude inaccuracies and timing errors between
I - 4
receiver and transmitter instrumentation (Podgorski et al., 2013b). Contami- nated early time and noisy late time gates were semi-automatically eliminated for each record. Late Palaeogene uplift along the Okavango-Kalahari-Zimbabwe axis and as- sociated subsidence of the Kalahari Basin disrupted the south-westerly flow of the Palaeo Chambeshi-Kafue-Upper Zambezi (PCKU) into the Indian Ocean via the proto Limpopo River. This Palaeogene tectonism also resulted in the damming of the Cubango-Cuito (proto Okavango River system) and Cuando rivers to form a closed fluvial/ lacustrine Kalahari Basin (Moore et al., 2012). Headward erosion of the Mid-Zambezi led to the capture and di- version of the PCKU system to the Indian Ocean via the lower Zambezi (Main et al., 2008). This was marked by a significant decrease in sedimenta- tion in the Kalahari Basin between the Pliocene to Early Pleistocene (Moore et al., 2012). Secondary faulting related to the south-westerly development of the EARS through the Mweru-Tanganyika-Kabompo Gorge axis resulted in the Okavango-Linyati fault system responsible for formation of the Machile Graben bound to the north-west by the Sesheke fault and to the south-east by the Chobe/Mambova Fault (Fig. 1) (Moore et al., 2012; McCar- thy, 2013). Early Pleistocene uplift of the Chobe Horst and related subsidence of the Machile Graben re-directed the flow of the PCKU system into the en- dorheic Kalahari Basin to result in the highest lake stand - palaeo Lake De- ception. This was later reduced to the 945 m level corresponding to the pal- aeo Lake Makgadikgadi following loss of the palaeo Chambeshi River flow as a result of capture by the Luapula River initiated by 2.2-1.8 Ma uplift of the Congo-Zambezi Watershed (McCarthy, 2013). The palaeo Lake Makgadikgadi was maintained by flow from the Cubango- Cuito, Cuando and Upper Zambezi -proto Kafue rivers and associated climat- ic feedback (Moore et al., 2012). At this point the proto Kafue River still drained into the Upper Zambezi via the Machile Flats. Early to Mid- Pleistocene Uplift along the Machile-Kafue Watershed finally severed the link between the Upper Zambezi and the Proto Kafue with consequent con- traction of the palaeo Lake Makgadikgadi to the 936 m shoreline. Incision of the Chobe Horst initiated during high lake stands by over topping maintained some flow in the Middle-Zambezi which in turn sustained riverine erosion and lead to the final breach at Mambova (Fig. 1) to establish the modern course of the Zambezi River in the Mid-Pleistocene. In addition, tectonic events re-organised the Okavango and Cuando river systems resulting in re- duced flow into the Palaeo Lake Makgadikgadi system which resulted in
I - 5
eventual drying up leaving behind remnants such as the Makgadikgadi Salt Pans (Main et al., 2008; Burrough et al., 2009; Moore et al., 2012; McCarthy, 2013). Thus endorheic conditions were initiated in the Kalahari Basin in the Late Palaeogene (approximately 25 Ma) with subsequent formation of palaeo Lake Makgadikgadi by Early Pleistocene (2.2 to 1.8 Ma). Early Mid- Pleistocene (more than 500 Ka) isolation of the proto-Kafue from the Upper Zambezi initiated desiccation of palaeo Lake Makgadikgadi until its final ex- tinction in the Mid-Late Pleistocene (100 Ka) to present (Moore et al., 2012). The objectives of this study were to map for the first time the paleo Lake Makgadikgadi sediments in the Machile flats of Zambia using airborne TEM methods; demonstrate the utility of TEM in the Machile-Zambezi geological setting for groundwater quality mapping; extend the regional coverage and spatial resolution of the work by Chongo et al. (2011) to encompass the entire Machile-Zambezi Basin using a combination of airborne and ground based TEM measurements; and to come up with a consistent conceptual explanation for the observed electrical resistivity variations. This was necessary in order to shed more light on and gain a better understanding of the Palaeo Lake Makgadikgadi system that extends into Zambia and in order to understand the occurrence of saline groundwater and how it relates to the geology. From a management perspective this was also important because it will lead to im- proved access to clean drinking water in rural and remote areas though a bet- ter understanding of the groundwater regime and quicker, efficient and cost effective means of evaluating groundwater for borehole placement.
2 Materials and methods 2.1 Description of the study area The Machile-Zambezi basin is located between latitude 16o - 17o 54’ S and longitude 24o 13’ – 26o 22’ E. It encompasses river basins for the three main northerly tributaries of the Zambezi namely, Loanja, Machile and Ngwezi, and a number of other smaller river basins. The outlet to the Machile- Zambezi Basin is defined at 17o 50.4’ S and 25o38’E (Fig. 2 ) and is there- fore hydrologically defined but encompasses the northern tip of the Makgadikgadi Basin that extends into Zambia which is called the Machili Basin by Moore and Larkin (2001) (Fig. 2). The dominant surface geology is unconsolidated sediments of Cenozoic Era with outcrops of basaltic and basement rocks in places particularly in the eastern areas where the presence of basement complex rocks is significant. The topography is characterised by
I - 6
a valley that is widest in the south-central areas and narrows towards the north-east with elevations ranging between 800 -1200 m amsl (Fig. 2).
Fig. 2: The Study area map depicting aerial TEM survey flight lines and ground based TEM points. The surface geology is predominantly unconsolidated Cenozoic sediments with significant outcrop of the basement complex in the eastern region after Persits et al. (2002). Elevation data from Jarvis et al. (2008). The outline of the Machile Basin after Moore and Larkin (2001) is also contrasted with the outline of the Machile-Zambezi Basin. 2.2 Transient electromagnetic method Based on the principles of electromagnetism, secondary currents can be in- duced in the subsurface by a time-varying electromagnetic field. The second- ary magnetic fields generated by the secondary currents can then be measured by appropriate electromagnetic receivers (Nabighian, 1988). The two main types of electromagnetic prospecting are frequency-domain electromagnetic (FEM) methods and time-domain (or transient) electromagnetic (TEM) meth- ods (Nabighian, 1988; Nabighian, 1991; Kirsch, 2006). TEM methods were used for this research. These comprised both ground based and airborne measurements. Airborne TEM measurements were conducted with the heli- copter-borne versatile time-domain electromagnetic (VTEM) system devel- oped by GEOTECH Limited who were contracted by the Ministry of Energy and Water Development (MEWD) of the Republic of Zambia to conduct the
I - 7
survey (GEOTECH, 2011). Ground-based TEM measurements were carried out with the WalkTEM system (ABEM, 2014) developed by Aarhus Univer- sity and ABEM Instrument AB. In both cases a central loop configuration was used. The VTEM system had a nominal terrain clearance of 48 m and nominal flight speed of 90 km/h and a total of 1000 line kilometres were flown. This comprised two sets of four parallel lines as shown in Fig. 2. One set was southwest to northeast trending and the other northwest to southeast trending. In both instances the line separation was 2000 m. The 24 ground- based TEM soundings were conducted on the eastern side of the Machile River. These, together with the 66 TEM soundings conducted by Chongo et al. (2011) and the 27 soundings collected for calibration form the ground based TEM dataset for this paper (Fig. 2). The total number of TEM sound- ings collected for the Machile-Zambezi Basin was therefore 117 soundings over 26, 260 km2. Field schematics, transmitter waveform and system re- sponse of both the airborne and ground-based deployment are illustrated in Fig. 3. A central loop configuration was used for the VTEM setup with concentric transmitter and receiver coils. The outer loop in form of a 12 sided polygon acted as the transmitter loop with diameter 26 m and 4 turns transmitting a peak current of 240 A. The inner loop had a diameter of 1.2 m and 100 turns and acted as the receiver coil (GEOTECH, 2011) (Fig. 3a). The WalkTEM system also used the central loop configuration with a single turn 40 m by 40 m transmitter loop and a 0.5 m by 0.5 m receiver coil in the centre with 20 internal turns(ABEM, 2014) (Fig. 3b). The VTEM transmitter current wave- form is 4.7 ms long with a peak at -1.646 ms. Receiver gates commence shortly after transmitter turn off at 0 ms with gate widths progressively in- creasing at later times. Gate centre times start at 21 µs and end at 10667 µs with a total of 45 gates (Fig. 3c). The WalkTEM system had two wave- forms; a high-moment waveform transmitting a 10 ms 8 A pulse and a low- moment waveform of amplitude 1.0 A and duration 10 ms. Turn off time was 5.9 µs and 10.4 µs for low and high moments, respectively. The low-moment response was measured over 23 gates with the first gate centre time at 3.6 µs and the last gate centre time at 364 µs. The high moment response was meas- ured over 37 time gates with the first gate centre time position at 3.6 µs and the last at 8850 µs. The low moment is optimised for early time gate meas- urements and thus near surface geological information whereas the high mo- ment is optimised for later time gate measurements and thus deeper surface geological information (Fig. 3d).
I - 8
Fig. 3: Field Schematics, transmitter waveform and system response for the WalkTEM and VTEM systems: (a) VTEM Setup; (b) WalkTEM Setup; (c) VTEM waveform and system response; and (d) WalkTEM waveform and system response. 2.3 Aerial survey The instrumentation for the airborne survey was mounted on an AS 350 B3 helicopter equipped with a Terra TRA 3000/ TRI40 radar altimeter and a NoVatel wide area augmentation system (WAAS) enabled OEM4-G2-3151W GPS navigation system with a sampling interval of 0.2 s. The TEM transmit- ter and receiver loops were carried as a sling-load beneath the helicopter with sampling interval at 0.1s. The VTEM decay sampling scheme was configured for 45 time gates. The first gate centre time was at 21 µs and the last gate centre time was at 10.667 ms (GEOTECH, 2011). The transmitter current waveform and receiver time gates are illustrated in Fig. 3c. Electromagnetic data collected for the survey was presented as db/dt (time derivative of mag- netic flux through the receiver coil). The data were imported into the Aarhus Workbench (HGG, 2014) where it was processed and inverted. The pro- cessing performed in the Workbench was on data that was already pre- processed by GEOTECH (2011) using an un-documented three stage digital
I - 9
filtering process to reject major broadband electromagnetic impulses or spherics caused by lightening events and to reduce system noise. The signal to noise ratio was further improved by the application of a low pass linear digital filter but the characteristics of this filter and its effects on the data were not documented by GEOTECH (2011). Thus the uncertainties about how the digital filtering was conducted by GEOTECH (2011) in addition to transmitter-receiver timing errors and amplitude shifts necessitated calibra- tion of the airborne TEM dataset with accurate ground based TEM measure- ments (Podgorski et al., 2013b) in order to obtain a more reliable inversion and interpretation. Interpretation of resistivity data obtained from TEM measurements is described in more detail in the Discussion section. 2.4 Ground survey The ground based TEM dataset comprised data from both ProTEM and WalkTEM equipment. The ProTEM equipment was set to measure at three different repetition rates under 20 gate mode. The repetition rates were desig- nated as u(237.5 Hz), v(62.5Hz) and H(25Hz) and the current was set at 3 A throughout with a standard square transmitter waveform characteristic of the ProTEM equipment. Thus for each location a complete sounding comprised one measurement for each of the three repetition rates. The gate times for the u repetition rate were optimised for early times whereas those for the H repe- tition rate were optimised for late times. The repetition rate v had an overlap between the gate times for the u and H repetition rates. The WalkTEM equipment used the concept of low moment and high moment. For each sounding, one measurement was made with the current setting at 1.0 A (low moment) and another with the current setting at 8.0 A (high moment). Fig. 3 (d) illustrates the concept of low and high moments used by the WalkTEM equipment. The low moment is optimised for early times and the high mo- ment is optimised for late times. Characteristics of the transmitter waveforms and receiver gate times were all controlled through user defined script files. The location of the ground based TEM soundings for the Machile-Zambezi Basin are shown in Fig. 2. Each sounding was inverted individually using AarhusInv (Auken and Christiansen, 2007; Auken et al., 2014) and the best fit model with residuals below 1 that made geologic and physical sense was selected as the most appropriate at that location. 2.5 Inversion methodology Inversion of the airborne dataset from the time derivative of the secondary electromagnetic field (db/dt) passing through the receiver coil to resistivity
I - 10
layered earth models was done using two separate but similar techniques. These are laterally constrained inversion (LCI) (Auken et al., 2005) and spa- tially constrained inversion (SCI) (Viezzoli et al., 2008). Using the LCI technique, each of the eight flight lines was inverted individually segment by segment. A typical flight line comprised two or three segments and with the SCI technique, all flight lines were inverted together as one inversion job. In both cases tight constraints were set on a smooth or minimum structure start- ing model given the relatively homogenous sedimentary terrain. Thus a 29 fixed-layer starting model was setup such that the mean apparent resistivity for each sounding was used as the initial resistivity for each layer. No a priori information was added. The thickness of the first layer was 7 m and that of the last layer was 60 m; i.e. the thicknesses increased progressively with depth. Vertical and lateral reference constraints on resistivity values were set to 2.5 and 1.3, respectively. The reference distance for the reference con- straints was set to 25 m with the with the power law factor (Christiansen et al., 2007) set to 0.5. 2.6 Calibration of the airborne dataset High quality ground based TEM soundings that were coincident with the air- borne TEM data within a distance of 50 m were used for calibration, and in total these were 14 soundings. At each of the 14 sites, a forward response based on the layered earth model from inversion of the ground based TEM sounding was compared with the corresponding VTEM sounding at the same location. A manual calibration was then conducted by first shifting the VTEM sounding curve vertically upwards to get the level right and then hori- zontally to get the time shift right. The amount of vertical and horizontal shift comprised the calibration parameters designated as shift factor and time shift, respectively. The procedure and considerations used for processing and cali- bration of the airborne dataset are similar to those of Podgorski et al. (2013b) who report a final time shift of 30 µs and a final shift factor of 1.44. For this study the time shift was determined as -44.929 ± 10.05µs, and the shift factor was determined as 1.071 ± 0.112 based on averaging of the calibration results at 14 sites. The difference in sign between the time shift factor reported by Podgorski et al. (2013b) and the one determined in this study probably has to do with the amount of system drift experienced by the respective VTEM sys- tems used each survey. The calibration was thus completed by subtracting 44.929 µs from all time gates and multiplying all raw data values (db/dt) with 1.071.
I - 11
A qualitative assessment of the effect of calibration on the inversion was conducted by comparing the inversion result from an LCI with un-calibrated data (Fig. 4a) to that of a similar LCI with calibrated data (Fig. 4b) and an- other also with calibrated data but using SCI (Fig. 4c).
Fig. 4: Comparison of horizontal resistivity thematic maps over the Loanja Alluvial Fan at depth interval 0 to 20 m for a 29 layer smooth inversion using (a) LCI with un-calibrated data; (b) LCI with calibrated data; and (c) SCI with calibrated data. 3 Results The results from the airborne and ground based TEM survey are presented below as maps and cross sections. These are extractions from the Machile- Zambezi geophysical database produced and maintained in the Aarhus Work- bench (HGG, 2014). The mean horizontal electrical resistivity map (Fig. 5) at depth interval 0 – 20 m was superimposed over the topography in the Machile-Zambezi Basin and depicts low electrical resistivity values (1 – 10 Ωm) confined to the low lying south-central region. At higher elevation par- ticularly in the north-eastern and north-western regions, the electrical resis- tivity values are very high in the range of 1000 Ωm. In between and sur- rounding the high and low electrical resistivity blocks are areas of moderate values (10 to 100 Ωm).
I - 12
Fig. 5: Mean horizontal electrical resistivity variations at depth interval 0 – 20 m.
The mean horizontal electrical resistivity map at 80-100 m (Fig. 6) is similar to that at 0 to 20 m (Fig. 5 above) except that the south-central low resistivity region is smaller in spatial extent and roughly coincides with the 945 m con- tour. In addition, the western and eastern areas with high resistivities are larger and elongated from south-west to north-east. North of the area with low electrical resistivity oriented towards the north east, is an elongated mod- erate resistivity region which together with the low resistivity block appear to be in valley structure between the western and eastern electrical resistivity blocks.
I - 13
Fig. 6: Mean horizontal electrical resistivity variations at depth interval 80 – 100 m.
The electrical resistivity structure along the northwest to southeast profile (Fig. 7a) indicates a low lying area with low electrical resistivity in the centre bound by an area with high electrical resistivity to the west (left hand side) and another area with high electrical resistivity to the east (right hand side) (Fig. 7a). The low electrical resistivities are interrupted in the shallow subsur- face by moderate resistivity anomalies. Along the southwest to northeast pro- file line a low resistivity or high conductivity region shaped in form of a pot with a handle extends from south (left hand side) to north (right hand side) with the pot handle pinching out towards the north. Underlying the pot handle is the moderate resistivity region which crops out in the northern areas. The main low resistivity pot is interrupted in the southern shallow subsurface by moderate anomaly resistivity values ranging from 10 – 100 Ωm (Fig. 7b).
I - 14
(a.)
(b.)
Fig. 7: (a) - Northwest to southeast cross section along profile line NW_SE (Fig. 5). (b)-Southwest to northeast cross section along profile line SW-NE (Fig. 5). The depth of investigation is the top of the faded part of the electrical resistivity cross sections.
I - 15
Evaluation of borehole log data from various boreholes in the Sesheke area formed a basis of inferring groundwater salinity from the electrical resistivity mapped by TEM. An example borehole log at borehole RV31 (borehole number 1 in Table 1 below) is shown in Fig. 8 whereas the location of the borehole logs considered are shown in Fig. 5 and Fig. 6 above.
Natural Gamma Fluid Long [Counts/s] Conductivity[µS/cm] Conductivity[mS/cm] 0
5
10
15
Depth [m]
20
25
30 020406080 0 200 400 10024 10 10 Fig. 8: Example borehole log from borehole RV31
From the borehole logs, formation factors (Tsallis et al., 1992; Neretnieks et al., 2001) were determined using the pore water resistivity (which was meas- ured as fluid conductivity with the borehole probe in µS/cm) and formation resistivity below the water table (which was measured as long conductivity in mS/m with an induction probe). The construction of most boreholes imple- mented a screen interval from the water table up to 3m above the bottom of the borehole. A single formation factor was thus assigned to each borehole by calculating the average formation factor from individual entries in the bore- hole log and the nominal formation resistivity around each borehole was de- fined as the product of the average pore water resistivity and the borehole formation factor. The formation factors and formation electrical resistivity values calculated for each borehole are shown in Table 1.
I - 16
Table 1: Pore water and formation resistivities, and formation factors derived from borehole logging in the Machile-Zambezi Basin. Locations of the boreholes are shown on figures 2, 5 and 6.
Location (UTM 35S) Pore Water Formation Resistivity Formation Nominal Formation Borehole Resistivity Threshold on non-saline General Lithology Factor Resistivity (Ωm) Northing (m) Easting (m) (Ωm) water (Ωm)
RV_31 8098183.40 237931.96 26.86 6.68 179.48 95.44 Sand
Clayey Sand and sand- RV_08 8099759.76 262479.77 17.31 5.74 99.31 81.97 stone
RV_29 8076339.01 211449.36 21.76 1.50 32.75 21.50 Sandstone/ Basalt
RV_12_02 8137425.81 299048.19 21.43 1.52 32.52 21.68 Sand/ Sandstone
RV_36 8068825.60 288165.00 15.72 1.54 24.19 21.98 Sand/ Sandstone
RV_01 8070051.34 231262.64 4.50 1.94 8.74 27.74 Sandy Clay
RV_26 8096066.40 289747.90 2.15 2.64 5.66 37.66 Clayey sand
I - 17
Also shown in Table 1 are formation resistivity thresholds of ground water salinity below which the groundwater is considered salty and above which it is considered fresh water or non-saline water. The formation resistivity thresh holds were calculated by multiplying the respective formation factor with the threshold for non-saline water which was taken as 700 µS/cm or 14.3 Ωm (FAO, 2014). Based on Table 2 the average formation electrical resistivity threshold for the Machile-Zambezi Basin is estimated as 44.0 ± 31.3 Ωm im- plying that in areas with very high clay content (i.e. high surface electrical conductivity), electrical resistivity values below 12.6 Ωm are indicative of salty groundwater whereas in areas with very little clay content the threshold is at 75.3 Ωm. This high degree of variability of formation factors is probably a consequence of the different geological domains from which the respective borehole logs were sampled but the values are typical of what can be ex- pected in sedimentary terrain comprising mostly clay, sand and sandstone (Salem, 2001). Similarly the average formation factor for the Machile- Zambezi Basin is 3.08 ± 2.19.
4 Discussion
Overall, the resistivity structure mapped by TEM in the Machile-Zambezi Basin Fig. 6 is characterised by three southwest to northeast trending areas of distinct electrical resistivity. These are: The westerly high electrical resistivity area with resistivity values greater than 100 Ωm; The middle area with low to moderate electrical resistivity (with values less than 100 Ωm) which is further subdivided into the south central area with low electrical resistivity and the central, south west to north east trending area with moderate electrical resistivity area; and The easterly high electrical resistivity area with resistivity values greater than 100 Ωm. The high resistivity structure of the westerly block is attributed to a typical Kalahari Basin stratigraphy (McCarthy and Haddon, 2005) of unconsolidated sands underlain by sandstone or basalt. A borehole profile at Munyeula (Kameyama, 2003) (Borehole 3 in Table 2) reveals a 12 m sequence of sand underlain by a 58 m succession of sandstone, which in turn overlies weath- ered and fractured Batoka Basalt. The electrical conductivity recorded in the
I - 18
Table 2: Selected borehole records in Machile-Zambezi Basin (Kameyama, 2003; Bäumle et al., 2007)
Borehole 1 Borehole 2 Borehole 3
Name Mbanga Name Kasaya Name Munyeula
Latitude 17.45241667 Latitude 17.45497222 Latitude 17.06072222
Longitude 24.96958333 Longitude 25.00333333 Longitude 24.66130556
EC (µS/cm) 7250 EC (µS/cm) 12180 EC (µS/cm) 263
Depth (m) Lithology Depth (m) Lithology Depth (m) Lithology
1 Fine sand 6 Clay 12 Fine sand
2 Clay 12 Fine sand 18 Sandstone
Clayed Fine Clayed fine 4 16 40 Sandstone sand sand Clayed fine 10 Fine sand 34 52 Sandstone sand Clayed Fine 14 50 Fine sand 58 Basalt sand Clayed fine 16 Clay 56 70 Basalt sand
22 Medium sand 60 Fine sand
Clayed Fine Clayed fine 26 62 sand sand
66 Fine sand 82 Fine sand
Borehole 4 Borehole 5 Borehole 6
Name Siamundele Name Sekute School Name Kamenyani
Latitude -17.17558 Latitude -17.65472 Latitude 17.06680556
Longitude 25.83619 Longitude 25.65889 Longitude 25.16697222
EC (µS/cm) - EC (µS/cm) - EC (µS/cm) 1932
Depth (m) Lithology Depth (m) Lithology Depth (m) Lithology
2 Top soil 6 Sand with clay 4 Fine sand
28 Granite 18 Sandstone 6 Clay
38 Granite 49 Basalt 22 Fine sand + silt
66 Granite 36 Sandstone
55 Quartz 68 Sandstone
I - 19
groundwater from this borehole is 263 µS/cm. This shows consistency be- tween the formation resistivity measured by the TEM methods and the intrin- sic high resistivity of the unconsolidated sands, sandstone, basalt and pore water at this location. The borehole record at Munyeula also shows the water table to be 4.4 m implying that a 4.4 m sequence of dry Kalahari sand sits on top of the water table and would typically register as a high resistivity layer with resistivity values greater than 500 Ωm. Below the water table, the resis- tivity values would be moderated down wards to below 100 Ωm until solid bed rock is encountered at which point they would also register high resistivi- ty values in the 1000 Ωm order of magnitude. Similarly, the surface geology of the easterly block is defined in the northeast by a northeast to southwest trending basement complex (Pre Cambrian geology) region and in the south- west by what appears to be a basin filled with Kalahari sediments and bound to the northeast and southwest by Mesozoic igneous horsts (refer to Fig. 2 above) herein referred to as the Sekute Basin. A borehole record at Siamun- dele (Bäumle et al., 2007) (Borehole 4 in Table 2) indicates a granitic for- mation consistent with the high formation resistivity detected by the TEM methods. In addition, a borehole record at Sekute School (Bäumle et al., 2007) (Borehole 5 in Table 2) indicates that the stratigraphy of the Sekute Basin is characterised by sand and clay (6 m) on top of sandstone (18 m) which in turn is underlain by basalt, a lithology log similar to the borehole profile at Munyeula in the easterly block and a similar relationship to the formation resistivity measured by the TEM methods is thus inferred. A summary of selected borehole records in the Sesheke and Kazungula area of south western Zambia in the 3 resistivity blocks of the Machile-Zambezi Ba- sin is shown in Table 2. The electrical resistivity structure of the middle area as mentioned above has two sub divisions. The first one is the south and central low electrical resis- tivity area that has an extent closely following the 945 m contour synony- mous with the Palaeo Lake Makgadikgadi shoreline (Fig. 5 above). However, the spatial extent of the low electrical resistivity area is much larger at shal- lower depths (Fig. 4 above) and extends north-easterly beyond the 945 m contour by about 32 km with a vertical rise of about 20 m. The low resistivity block therefore overlies and thins out into the southwest to northeast trending moderate resistivity block northeast of the 945 m contour. The moderate re- sistivity block is abruptly truncated in the vicinity of the 945 m contour and has a gentle rise towards the northeast until it overlies the thinning out low electrical resistivity area (Fig. 6b above). Existence of the low resistivity
I - 20
block can be attributed to sediments related to the palaeo Lake Makgadikgadi with evidence of this correlation arising from the borehole profiles at Kasaya and Mbanga (Table 2). These depict alternating and mixed sequences of fine sand and clay with very high pore-water salinity; 7,250 and 12,180 µS/cm at Mbanga and Kasaya (boreholes 1 and 2 in Table 2). Milzow et al. (2009) ob- served a similar setting in the Okavango Delta and attributed the alternating and mixed sequences of clay and sand to alternating cycles of fluvial and la- custrine depositional setting. The work by Banda et al. (Technical University of Denmark/University of Zambia, Unpublished data, 2014) used borehole logging, sediment core characterization and geochemical analysis to postulate that palaeo Lake Makgadikgadi sediments were responsible for groundwater salinity in the Machile Basin (after Moore and Larkin (2001)) which is con- gruent with the low electrical resistivity area. Banda et al. (Technical Univer- sity of Denmark/University of Zambia, Unpublished data, 2014) thus describe a lithological setting of alternating and mixed sequences of silty-clay, clay and sand attributed to fluvial and lacustrine deposition associated with vary- ing climatic conditions containing saline groundwater entrapped at the time of deposition of the lake sediments. In addition, the fact that the low electri- cal resistivity area closely follows the 945 m contour at the depth interval 80 -100 m is further supporting evidence of palaeo Lake Makgadikgadi sedi- ments in the Machile-Zambezi Basin. However, this is in contrast to the ex- tension of the low resistivity area beyond the 945 m contour towards the northeast for depth interval 0 to 20 m (refer to Fig. 5 above). A possible ex- planation for this is likely related to uplift along the Kafue-Machile Water- shed and associated closure of the palaeo-Kafue/Upper Zambezi link (McCar- thy, 2013). Alternatively this could also be as a result of shallow ground wa- ter related to the Palaeo Lake Makgadikgadi which was subjected to evapo- concentration. According to Moore et al. (2012) the palaeo Chambeshi River, which was a combination of the proto Kafue and proto Chambesh river systems, used to flow in to the Makgadikgadi Basin through the Machile Flats, where it would join the Upper Zambezi River. This south westerly flow of the palaeo Chambeshi was associated with deposition of Kalahari super group sediments (Moore and Larkin, 2001) which correlated with the Barotse sand member (Money, 1972). Thus the deposition of sandstone is believed to have occurred in braided streams (McCarthy and Haddon, 2005) indicative of low slope ter- rain which could have easily been below or at the level of the 945 m contour by the time the sands were fully deposited. Lake sediments would have then
I - 21
been deposited in a shallow stretch up to the present northern extent of the low resistivity block followed by uplift along the Kafue-Machile Watershed (Moore et al., 2012). This up-lift could have given rise to the gentle north easterly slope responsible for occurrence of part of the low resistivity block above the 945 m contour. However, Banda et al. (Technical University of Denmark/University of Zambia, Unpublished data, 2014) postulate that this was probably the result of evapo-concentration of shallow groundwater from the palaeo lake system. A borehole record at Kamenyani (Borehole 6 in Table 2) reveals an uppermost 4 m dine sand, overlying a 6 m sequence of clay, which again is underlain by a 22 m sequence of silt and sand suggestive of a lacustrine setting, overlying is a 4m layer of sand, probably as a result of re- cent deposition by aeolian processes. The silt and sand sequence is underlain by sandstone. One implication of this interpretation is that the moderate resis- tivity area could have had resistivity values as low as in the low electrical resistivity area. However uplift along the Kafue-Machile Watershed resulted in a net groundwater flow towards the southwest resulting in flushing of the system with fresh groundwater resulting in moderated electrical conductivity values such as the 1,932 µS/cm observed at Kamenyani compared to the 12,180 µS/cm observed at Kamenyani.
5 Conclusion The resistivity structure of the Machile-Zambezi Basin is characterised by three distinct southwest to northeast trending electrical resistivity areas. These are:
The high electrical resistivity westerly area with resistivity values greater than 100 Ωm associated with a sand, sandstone and Basalt downward se- quence; The middle low to moderate electrical resistivity area with values less than 100 Ωm. This is further subdivided into the south and central low electrical resistivity area (with values less than 13 Ωm) and the south west to north east trending moderate electrical resistivity area (with values between 10 – 100 Ωm). The low electrical resistivity area is congruent with the northern tip of the palaeo Lake Makgadikgadi sediments extending into south- western Zambia and is filled with saline groundwater trapped at the time of deposition of the sediments ; and
I - 22
The easterly high resistivity block with resistivity values greater than 100 Ωm associated with Basement Complex and Kalahari stratigraphic set- tings. The southwest to northeast orientation of the resistivity blocks is a clear indi- cation of continuity of the Okavango-Linyati Fault System and associated palaeo Lake Makgadikgadi into south-western Zambia. The combination of ground based and airborne TEM methods was thus effec- tive in mapping the regional electrical resistivity structure of the Machile- Zambezi Basin from which groundwater salinity variations could be inferred in addition to the regional tectonic structure or geological fault system. Cali- bration of VTEM data with accurate ground based WalkTEM data ensured a strong agreement between the airborne TEM inverted using an SCI scheme and ground based data inverted as single site 1D inversions.
6 Acknowledgements Gratitude is given to the governments of the Kingdom of Denmark and the Republic of Zambia for funding this research under the capacity building ini- tiative for the Zambian Water Sector through the University of Zambia Inte- grated Water Resources Management Centre. Thanks also go to the Hydro- geophysics group of Aarhus University for assistance in processing and in- version of the airborne dataset. Special thanks also go to colleagues and friends who assisted in field work especially Kebby Kapika (District Water Officer-Sesheke) and Mweemba Sinkombo (District Water Officer-Nalolo) who was then a Post-graduate diploma student at UNZA. Logistical support by Ms. Ingrid Mugamya Kawesha of UNZA IWRM Centre is also acknowl- edged.
7 References
ABEM 2014. WalkTEM User's Guide. Stockholm: ABEM Geophysics.
ARAD, A. 1984. RELATIONSHIP OF SALINITY OF GROUNDWATER TO RE- CHARGE IN THE SOUTHERN KALAHARI DESERT. Journal of Hydrology, 71, 225-238.
AUKEN, E. & CHRISTIANSEN, A. V. 2007. Manual for the Inversion Programme em1dinv. In: GROUP, H. (ed.) version 5 ed. Aarhus Denmark: Hydrogeophysics Group, University of Aarhus.
I - 23
AUKEN, E., CHRISTIANSEN, A. V., JACOBSEN, B. H., FOGED, N. & SORENSEN, K. I. 2005. Piecewise 1D laterally constrained inversion of resistivity data. Geophysical Prospecting, 53, 497-506.
AUKEN, E., CHRISTIANSEN, A. V., KIRKEGAARD, C., FIANDACA, G., SCHAMP- ER, C., BEHROOZMAND, A. A., BINLEY, A., NIELSEN, E., EFFERSO, F., CHRISTENSEN, N. B. I., SORENSEN, K., FOGED, N. & VIGNOLI, G. 2014. An overview of a highly versatile forward and stable inverse algorithm for airborne, ground-based and borehole electromagnetic and electric data. Exploration Geophysics.
AUKEN, E., JØRGENSEN, F. & SØRENSEN, K. I. 2003. Large Scale TEM Investigations for Groundwater. Exploration Geophysics, 34, 88–194.
BAUER-GOTTWEIN, P., GONDWE, B. N., CHRISTIANSEN, L., HERCKENRATH, D., KGOTLHANG, L. & ZIMMERMANN, S. 2010. Hydrogeophysical exploration of three-dimensional salinity anomalies with the time-domain electromagnetic method (TDEM). Journal of Hydrology, 380, 318-329.
BAUER, P., HELD, R. J., ZIMMERMANN, S., LINN, F. & KINZELBACH, W. 2006. Coupled flow and salinity transport modelling in semi-arid environments: The Shashe River Valley, Botswana. Journal of Hydrology, 316, 163-183.
BÄUMLE, R., NEUKUM, C., NKHOMA, J. & SILEMBO, O. 2007. Groundwater Re- sources of Southern Province (Phase 1) ANNEX 3 GReSP Borehole Data. In: MINIS- TRY OF ENERGY AND WATER DEVELOPMENT, D. O. W. A. (ed.). Lusaka: De- partment of Water Affairs.
BOUCHER, M., FAVREAU, G., DESCLOITRES, M., VOUILLAMOZ, J.-M., MAS- SUEL, S., NAZOUMOU, Y., CAPPELAERE, B. & LEGCHENKO, A. 2009. Contri- bution of geophysical surveys to groundwater modelling of a porous aquifer in semiar- id Niger: An overview. Comptes Rendus Geoscience, 341, 800-809.
BURROUGH, S. L., THOMAS, D. S. G. & BAILEY, R. M. 2009. Mega-Lake in the Kala- hari: A Late Pleistocene record of the Palaeo Lake Makgadikgadi system. Quaternary Science Reviews, 28, 1392-1411.
CAMPBELL, G., JOHNSON, S., BAKAYA, T., KUMAR, H. & NSATSI, J. 2006. Air- borne geophysical mapping of aquifer water quality and structural controls in the Lower Okavango Delta, Botswana. South African Journal of Geology, 109, 475-494.
CHONGO, M., WIBROE, J., STAAL-THOMSEN, K., MOSES, M., NYAMBE, I. A., LARSEN, F. & BAUER-GOTTWEIN, P. 2011. The use of Time Domain Electromag- netic method and Continuous Vertical Electrical Sounding to map groundwater salini- ty in the Barotse sub-basin, Zambia. Physics and Chemistry of the Earth, 36, 798-805.
CHRISTIANSEN, A. V., AUKEN, E., FOGED, N. & SORENSEN, K. I. 2007. Mutually and laterally constrained inversion of CVES and TEM data: a case study. Near Surface Geophysics, 5, 115-123.
I - 24
DANIELSEN, J. E., DAHLIN, T., OWEN, R., MANGEYA, P. & AUKEN, E. 2007. Geo- physical and Hydrogeologic Investigation of Groundwater in the Karoo Stratigraphic Sequence at Sawmills in Nothern Matebeleland, Zimbabwe: A Case History. Hydro- geology Journal, 15, 945-960.
EZERSKY, M., LEGCHENKO, A., AL-ZOUBI, A., LEVI, E., AKKAWI, E. & CHALI- KAKIS, K. 2011. TEM study of the geoelectrical structure and groundwater salinity of the Nahal Hever sinkhole site, Dead Sea shore, Israel. Journal of Applied Geophysics, 75, 99-112.
FAO. 2014. The use of saline waters for crop production [Online]. Food and Agricultural Organisation, Natural Resources Management and Environment Department. Availa- ble: http://www.fao.org/docrep/t0667e/t0667e05.htm#classification of saline waters 09/09/2014].
GEOTECH 2011. Survey and logistics report on a helicopter borne versatile time domain electromagnetic survey on the Zambezi River Basin Kazungula Zambia for Ministry of Energy and Water Development (Republic of Zambia). West Indies: GEOTECH AIR- BORNE LIMITED.
GUERIN, R., DESCLOITRES, M., COUDRAIN, A., TALBI, A. & GALLAIRE, R. 2001. Geophysical surveys for identifying saline groundwater in the semi-arid region of the central Altiplano, Bolivia. Hydrological Processes, 15, 3287-3301.
HERCKENRATH, D. 2012. Informing groundwater models with near-surface geophysical data.
HGG. 2014. Aarhus Workbench [Online]. Aarhus University: Hydrogeophysics Group. Available: http://hgg.au.dk/software/aarhus-workbench/ [Accessed 30-06-2014 2014].
JARVIS, A., REUTER, H. I., NELSON, A. & GUEVARA, E. 2008. Hole-filled seamless SRTM data V4 [Online]. International Centre for Tropical Agriculture (CIAT). Available: http://srtm.csi.cgiar.org [Accessed 22/06/2014.
KAMEYAMA, N. 2003. Completion Report: JICA borehole drilling project-Sesheke. Lu- saka: Ministry of Energy and Water Development.
KGOTLHANG, L. P. 2008. Application of airborne geophysics in large scale hydrological mapping; Okavango Delta, Botswana. PhD Doctoral and Habilitation Thesis, ETH Zurich.
KIRSCH, R. (ed.) 2006. Groundwater Geophysics, A Tool for Hydrogeology, Flintbek, Germany: Springer.
LAND, L. A., LAUTIER, J. C., WILSON, N. C., CHIANESE, G. & WEBB, S. 2004. Geo- physical monitoring and evaluation of coastal plain aquifers. Ground Water, 42, 59- 67.
LEVI, E., GOLDMAN, M., HADAD, A. & GVIRTZMAN, H. 2008. Spatial delineation of groundwater salinity using deep time domain electromagnetic geophysical measure- ments: A feasibility study. Water Resources Research, 44.
I - 25
MAIN, M. P. L., MOORE, A. E., WILLIAMS, H. B. & COTTERILL, F. P. D. 2008. The Zambezi River. Large Rivers: Geomorphology and Management, 311-332.
MCCARTHY, T. S. 2013. The Okavango delta and its place in the geomorphological evo- lution of Southern Africa. South African Journal of Geology, 116, 3-54.
MCCARTHY, T. S. & HADDON, I. G. 2005. The Mesozoic-Cenozoic Interior Sag Basins of Central Africa: The Late Cretaceous-Cenozoic Kalahari and Okavango Basins. Journal of African Earth Sciences, 43, 316-333.
MELLOUL, A. J. & GOLDENBERG, L. C. 1997. Monitoring of seawater intrusion in coastal aquifers: Basics and local concerns. Journal of Environmental Management, 51, 73-86.
MILZOW, C., KGOTLHANG, L., BAUER-GOTTWEIN, P., MEIER, P. & KINZEL- BACH, W. 2009. Regional review: the hydrology of the Okavango Delta, Botswana- processes, data and modelling. Hydrogeology Journal, 17, 1297-1328.
MONEY, N. J. 1972. An Outline of the Geology of Western Zambia. In: DRYSDALL, A. R. & JOHNSON, R. L. (eds.) Records of the Geological Survey Volume 12. Lusaka: Geological Survey Department.
MOORE, A. E., COTTERILL, F. P. D. & ECKARDT, F. D. 2012. THE EVOLUTION AND AGES OF MAKGADIKGADI PALAEO-LAKES: CONSILIENT EVIDENCE FROM KALAHARI DRAINAGE EVOLUTION SOUTH-CENTRAL AFRICA. South African Journal of Geology, 115, 385-413.
MOORE, A. E. & LARKIN, P. A. 2001. Drainage evolution in south-central Africa since the breakup of Gondwana. South African Journal of Geology, 104, 47-68.
NABIGHIAN, M. N. 1988. Electromagnetic Methods in Applied Geophysics Vol. 1.
NABIGHIAN, M. N. 1991. Electromagnetic Methods in Applied Geophysics Vol. 2. Elec- tromagnetic Methods in Applied Geophysics Vol. 2.
NENNA, V., HERCKENRATH, D., KNIGHT, R., ODLUM, N. & MCPHEE, D. 2013. Application and evaluation of electromagnetic methods for imaging saltwater intru- sion in coastal aquifers: Seaside Groundwater Basin, California. Geophysics, 78.
NERETNIEKS, I., OHLSSON, Y. & LÖFGREN, M. 2001. Formation factor determina- tions by in-situ resistivity logging. Materials Research Society Symposium - Proceed- ings, 663, 1197-1205.
PERSITS, F., AHLBRANDT, T., TUTTLE, M., CHARPENTIER, R., BROWNFIELD, M. & TAKAHASHI, K. 2002. MAP SHOWING GEOLOGY, OIL AND GAS FIELDS AND GEOLOGIC PROVINCES OF AFRICA, VER. 2.0 [Online]. Denver: United States Geological Survey. Available: http://pubs.usgs.gov/of/1997/ofr-97-470/OF97- 470A/ [Accessed 19/07/2014.
I - 26
PODGORSKI, J., AUKEN, E., SCHAMPER, C., VEST CHRISTIANSEN, A., KALSCHEUER, T. & GREEN, A. 2013b. Processing and inversion of commercial helicopter time-domain electromagnetic data for environmental assessments and geo- logic and hydrologic mapping. Geophysics, 78, E149-E159.
PODGORSKI, J. E., GREEN, A. G., KGOTLHANG, L., KINZELBACH, W. K. H., KALSCHEUER, T., AUKEN, E. & NGWISANYI, T. 2013a. Paleo-megalake and paleo-megafan in southern Africa. Geology, 41.
SALEM, H. S. 2001. Determination of porosity, formation resistivity factor, Archie ce- mentation factor, and pore geometry factor for a glacial aquifer. ENERGY SOURCES, 23, 589-596.
SATTEL, D. & KGOTLHANG, L. 2004. Groundwater Exploration with AEM in the Boteti area, Botswana. Exploration Geophysics, 35, 147-156.
THOMAS, D. S. G. & SHAW, P. A. 1990. The deposition and development of the Kalaha- ri Group sediments, Central Southern Africa. Journal of African Earth Sciences (and the Middle East), 10, 187-197.
THOMAS, D. S. G. & SHAW, P. A. 1993. THE EVOLUTION AND CHARACTERIS- TICS OF THE KALAHARI, SOUTHERN AFRICA. Journal of Arid Environments, 25, 97-108.
THOMAS, D. S. G. & SHAW, P. A. 2002. Late Quaternary environmental change in cen- tral southern Africa: new data, synthesis, issues and prospects. Quaternary Science Reviews, 21, 783-797.
TSALLIS, C., CURADO, E. M. F., DESOUZA, M. D., ELIAS, V. L., BETTINI, C., SCU- TA, M. S. & BEER, R. 1992. GENERALIZED ARCHIE LAW - APPLICATION TO PETROLEUM RESERVOIRS. Physica A, 191, 277-283.
USGS. 2014. Shuttle Radar Topography Mission (SRTM) [Online]. College Park, Mary- land: Global Land Cover Facility, University of Maryland. Available: ftp://ftp.glcf.umd.edu/glcf/SRTM/Degree_Tiles/ [Accessed 20/06/2014.
VIEZZOLI, A., CHRISTIANSEN, A. V., AUKEN, E. & SORENSEN, K. 2008. Quasi-3D modelling of airborne TEM data by spatially constrained inversion. Geophysics, 73.
WORTHINGTON, P. F. 1977. GEOPHYSICAL INVESTIGATIONS OF GROUNDWA- TER RESOURCES IN THE KALAHARI BASIN. Geophysics, 42.
I - 27 II
Mapping localized freshwater anomalies in the brackish Paleo-Lake sediments of the Machile- Zambezi Basin with transient electromagnetic sounding, geoelectrical imaging and induced polarization
Chongo, M., Christiansen, A.V., Fiandaca, G., Nyambe, I.A., Larsen, F., and Bauer-Gottwein, P.
Submitted to Journal of Applied Geophysics (Manuscript No: - AP- PGEO3721)
Mapping localized freshwater anomalies in the brackish Paleo-Lake sediments of the Machile-Zambezi Basin with transient electromagnetic sounding, geoelectrical imaging and induced polarization
Mkhuzo Chongo1,*, Anders Vest Christiansen2, Gianluca Fiandaca2, Imasiku A Nyambe3, Flemming Larsen4 and Peter Bauer-Gottwein1
1Technical University of Denmark, Department of Environmental Engi- neering, Miljøvej, Building 113, 2800 Kgs. Lyngby, Denmark 2Aarhus University, Department of Geoscience, C.F. Møllers Alle 4, 8000-Århus C, Denmark 3University of Zambia, School of Mines, Department of Geology, Great East Road Campus, P.O Box 32379, Lusaka, Zambia 4 Geological Survey of Denmark and Greenland (GEUS), Geochemis- try, Øster Volgade 10, 1350 København K, Denmark
Submitted to Journal of Applied Geophysics (Manuscript No: - AP- PGEO3721)
*Email: [email protected]; [email protected].
Abstract A recent airborne TEM survey in the Machile-Zambezi Basin of south west- ern Zambia revealed high electrical resistivity anomalies (around 100 Ωm) in a low electrical resistivity (below 13 Ωm) background. The near surface (0- 40 m depth range) electrical resistivity distribution of these anomalies ap- peared to be coincident with superficial features related to surface water such as alluvial fans and flood plains. This paper describes the application of tran- sient electromagnetic soundings and multiple electrode profiling using geo- electrics and time domain induced polarization to evaluate a freshwater lens across a flood plain on the northern bank of the Zambezi River at Kasaya in south western Zambia. Coincident TEM and CVES measurements were con- ducted across the Simalaha Plain from the edge of the Zambezi River up to 6.6 km inland. The resulting TEM, direct current and induced polarization data sets were inverted using a new mutually and laterally constrained joint inversion scheme. The resulting inverse model sections depict a freshwater lens sitting on top of a regional saline aquifer. The fresh water lens is about 60 m thick at the boundary with the Zambezi River and gradually thins out and deteriorates in water quality further inland. It is postulated that the freshwater lens originated as a result of interaction between the Zambezi Riv- er and the salty aquifer under evapotranspiration. Similar high electrical re- sistivity features were also associated with other surface water features locat- ed in the airborne survey area.
II - 1 1 Introduction The interaction between surface water and groundwater has been studied ex- tensively around the world (Winter, 1999; Sophocleous, 2002; Westbrook et al., 2005; Milosevic et al., 2012; Shanafield and Cook, 2014; Zhou et al., 2014) using different approaches, and increasingly geophysical methods are being incorporated into such studies. For example, Bauer et al. (2006) de- scribed the process of salt accumulation on islands within the Okavango Del- ta, related to the interaction between surface water and groundwater under evapo-concentration using a combination of electrical resistivity tomography (ERT) and hydrodynamic modelling. Sonkamble et al. (2014) evaluated the extent of aquifer pollution from industrial effluent across the flood plain of the Palar River at Ambur Town (India) using 1D and 2D geo-electrics corre- lated with in-situ water quality data and ground penetrating radar. They found that 60.4% of the study area (particularly in the shallow aquifers underlying the flood plain and valley fills) corresponding to a spatial extent of 33.4 km2 was polluted. Their criterion for pollution threshold was based on a formation electrical resistivity of less than 20 Ωm representing pore water electrical conductivity in the order of 15,345µS/cm. Shalem et al. (2014) studied the interaction of the Alexander River with groundwater as it cuts its way across a mostly sandy Quaternary coastal aquifer on the eastern coast of the Medi- terranean Sea. They accomplished this using a combination of 3 continuous vertical electrical sounding (CVES) transects of not more than a few 100 m each; single site transient electromagnetic (TEM) measurements; in-situ elec- trical conductivity, temperature and sediment sampling during borehole drill- ing; groundwater and stream level measurements; salinity and depth profiles along the river bed; and evaluation of regional water levels and chloride con- centrations from databases of the Hydrological survey of Israel. From this they observed that sea encroachment resulted in a seasonal increase of river water salinity from which they could infer the extent to which the river inter- acted with the underlying aquifer. Aquifer salinization due to the river was found to be limited to several meters from the river bed and less than 5 m be- low ground level within 500 m from the coastline. However, about 1,800 m from the coast, the river salinized the aquifer more than 50 m away towards the north with fresh water occurring 120 m away at depths of 20-25 m. Thus geophysical measurements are well suited for gathering data at high spatial resolution and an overall assessment strategy using a combination of different geophysical methods can be advantageous (Brodie et al., 2007). In this regard, TEM (Harthill, 1976; Nabighian, 1991; Danielsen et al., 2003;
II - 2 Xue et al., 2012), direct current geo-electrics (DC) (Loke, 1999; Dahlin, 2001; Aizebeokhai, 2010; Loke et al., 2013) and induced polarization (IP) (Bertin and Loeb, 1969; Dahlin et al., 2002; Titov et al., 2002; Fiandaca et al., 2012; Fiandaca et al., 2013) techniques are quite suitable for environmen- tal and hydro-geological investigations particularly in sedimentary terrain. Traditionally, TEM, DC and IP techniques have been deployed separately even for investigations at the same study site (Nassir et al., 2000; Guerin et al., 2001; Bauer et al., 2006; Ezersky et al., 2011; Vaudelet et al., 2011) alt- hough it is now common to have instrumentation that measures both DC and IP in the same field setup (Aristodemou and Thomas-Betts, 2000; Marescot et al., 2008). As a result, different types of datasets are quite often generated for the same physical or environmental phenomenon by inverting each type of dataset individually. Nevertheless, major benefits can be derived from joint inversion of different types of data that observe the same phenomenon and can lead to more accurate interpretations. Thus many studies have successful- ly used one form of joint inversion or another such as DC-TEM (Albouy et al., 2001; Christiansen et al., 2007; Danielsen et al., 2007), and MRS-TEM (Behroozmand et al., 2012; Vilhelmsen et al., 2014). Examples of DCIP joint inversions are scarce in the literature with the normal practice being to inde- pendently invert the DC and IP data either as separate inversion jobs or in one inversion job but without any of the datasets influencing the other during the inversion process. Furthermore, joint DCIP-TEM inversions have not been reported in the literature before. This paper therefore presents a first case study of joint inversion of DCIP-TEM data. The focus of this paper is on local scale electrical resistivity anomalies de- rived from interpreting high resolution regional scale TEM data in terms of surface water/ groundwater interaction in the Machile-Zambezi Basin. The objectives were to describe the occurrence of high electrical resistivity anom- alies in the low electrical resistivity background environment of the Machile- Zambezi Basin; conduct high resolution local scale TEM and direct current- induced polarization (DCIP) CVES measurements along a transect cutting across an area exhibiting electrical resistivity anomalies; evaluate the bene- fits of joint inversion of the local scale TEM and DCIP data in comparison to separate inversions; and to evaluate the inverse resistivity section in terms of surface water groundwater interaction taking place at the local site.
II - 3 2 Materials and methods 2.1 Study site The study area is in the southern central low lying areas of the Machile- Zambezi Basin on the northern banks of the Zambezi River. The area is drained by three main tributaries of the Zambezi namely Loanja, Machile (or Kasaya) and Ngwezi. The Machile and Ngwezi streams cut across flood plains before entering the Zambezi but the Loanja forms an inland delta or alluvial fan instead. The Loanja alluvial fan and the Simalaha flood plain (bound by Kasaya River to the west and the Zambezi River to the south) were the two local areas of interest for this study. However, the combined TEM/CVES transect is reported only for the Simalaha flood plain (Fig. 1).
Fig. 1: The study area depicting Loanja alluvial fan, Simalaha flood plain, Kasaya Tran- sect and TEM flight paths. Satellite image courtesy of ESRI (2014). Note that the Simalaha Profile is shorter than the Kasaya Transect since there is no more airborne electromagnetic data outside the flight lines, whereas the CVES and TEM data on the Kasaya Transect ex- tend beyond the flight path where the Simalaha Profile ends up to the Zambezi River. 2.2 Data collection and pre-processing Airborne data was conducted along 8 flight lines totalling 1000 line kilome- tres using the VTEM system (GEOTECH, 2011). Four of the flight lines were oriented northeast to southwest whereas the other 4 were oriented from
II - 4 northwest to southeast (Fig. 1 in Section 2.1). Details about the airborne sur- vey and about the processing, inversion and interpretation of the collected TEM data are given in Chongo et al. (in press, 2014). Cross sections of the airborne TEM data along the Loanja and Simalaha profiles (Fig. 1 in Section 2.1) are shown in Fig. 2 (a.) and (b.) respectively. These depict superficial electrical resistivity anomalies in an otherwise low electrical resistivity (sa- line environment) and were the basis of the detailed local scale study con- ducted on the Kasaya transect presented in this paper.
(a.)
(b.)
Fig. 2: (a) Electrical resistivity cross section along Loanja Profile (Fig. 1 in Section 2.1) from the airborne TEM data. (b) Interpolated electrical resistivity cross section along Simalaha Profile (Fig. 1 in Section 2.1) from the airborne TEM data. Note that Loanja and Simalaha profiles are not drawn to scale nor are they the same length since Loanja Profile (about 106 km long) is along a flight line whereas Simalaha Profile (about 5.7 km long) cuts across flight lines and as a result has a more limited data extent.
The detailed local scale geophysical investigation conducted across the Sima- laha Plain at Kasaya (Fig. 1 in Section 2.1) comprised:-
6.6 km of CVES (Loke, 1999; Nassir et al., 2000; Loke et al., 2013) meas- urements at 5 m electrode spacing using the gradient array (Dahlin and Zhou, 2006) with 25,003 data points. The Terrameter LS (ABEM(a), 2012) was used for the CVES to measure both direct current electrical resistivity (DC) (Loke et al., 2013) and time domain induced polarization (IP)(Johnson, 1984) hence the term DCIP to denote the combination of DC and IP measurements in a roll along setup; and 64 single site TEM (Kirsch, 2006) soundings using the Aarhus Universi- ty/ABEM WalkTEM system; and another set of 64 central loop TEM soundings using the Geonics ProTEM 47D instrument (Geonics, 2006) at
II - 5 the same positions as the WalkTEM soundings. Thus the total number of TEM soundings along the transect line was 128 spaced at approximately 100 m along the 6.6 km transect line per pair of WalkTEM/ ProTEM 47D soundings. However data from the ProTEM instrument was not used for this paper. The DCIP data was pre-processed by removing all data points with negative electrical resistivity and data variations greater than 1.5 %. The data that was removed this way represented only 3.3 % of the original data i.e. 851 filtered out measurements from a total of 25, 857 DCIP measurements. The data was then imported into the Aarhus Workbench with the IP data gated into 10 channels. Data processing in the Workbench comprised semi-automatic re- moval of bad IP data by setting a maximum slope change for the IP decay curves followed by visual inspection of the DC and IP data points along the profile and consequent disabling of the outliers. The DCIP noise model was set to 1.03 uniform standard deviation (USTD) on DC and 1.15 USTD on IP whereas the threshold on voltage was set to 2.0 mV. For the WalkTEM data we used data from 77.6 µs to 2.84 ms, focusing on the deep information only. 2.3 Inversion methodology DC and TEM data were inverted separately using the 1D laterally constrained inversion (LCI) (Auken et al., 2005) scheme. Subsequently, a joint inversion using the mutually and laterally constrained inversion scheme of Christiansen et al. (2007) was conducted on the DC and TEM data as a single inversion. This was then extended to include IP data using the Cole-Cole model setup (Fiandaca et al., 2012; Gazoty et al., 2012b) so that the final inversion was a joint inversion of DCIP and TEM data. Thus the DCIP and TEM model pa- rameters being modelled comprised (intrinsic chargeability (M0), frequency dependence constant (c), time constant (τ), formation electrical resistivity (ρ) and layer thicknesses). The inversion algorithm AarhusInv (Auken et al., 2014) was used for all inversions presented in this paper. As mentioned in Section 2.2, the interval of TEM soundings along the Ka- saya transect was approximately every 100 m whereas the DCIP data was col- lected with 5 m gradient array electrode spacing. The lateral constraints on the TEM models were setup such that each TEM model was constrained only to the adjacent TEM model on either side along the transect line. Similarly, each DCIP model was constrained only to the adjacent DCIP models. Treat- ment of mutual TEM/DCIP constraints is explained below.
II - 6 The reference lateral constraint on electrical resistivities was set to 0.3 and scaled according to: ∗ (1)
where Ci = lateral constraint on resistivity [dimensionless fraction]; Cr is the reference constraint [dimensionless fraction]; d is the distance between re- spective models [m]; and dr is the reference distance which was set to 10 m. Furthermore, the reference constraint on depths was set to 1 m and scaled according to depth so that the deeper layers had relatively tighter constraints. The constraint values mentioned above can be considered as medium for the resistivity values and tight for the depths. These were used because they were found to give a reasonable inversion result using trial and error. Lastly, mutual constraints were applied between the TEM and DCIP models using depths of layers to set the constraint width and scaled according to the power law given above (equation 1). Thus, the deeper layers had wider and tighter constraints between TEM and DCIP models because the constraints were only applied if the distance between respective TEM and DCIP models was less than or equal to the layer depth. The reference constraint between TEM and DCIP models was set to 0.1 but the reference depth was the same as for the lateral constraints (Christiansen et al., 2007). 2.4 Petro-physical considerations The petro physical relation for the Kasaya area was estimated using the fol- lowing equations (Rhoades et al., 1989; Mualem and Friedman, 1991; Kirsch, 2009), ∗ (2) ∅