The Southern African Magnetotelluric Experiment (SAMTEX) involved the collection of data at over 700 sites in Archean to Proterozoic southern Africa, spanning features including the Kalahari Craton, Bushveld Complex and voluminous kimberlites. Here, we present the first 3D inversions of the full SAMTEX dataset. In this paper, we focus on assessing the robustness of the 3D models by comparing two different inversion codes, jif3D and ModEM, and two different subsets of the data, one containing all acceptable data and the other containing a smaller selection of undistorted, high-quality data. Results show that the main conductive and resistive features are imaged by all inversions, including deep resistive features in the central Kaapvaal Craton and southern Congo Craton and a lithospheric-scale conductor beneath the Bushveld Complex. Despite this, differences exist between the jif3D and ModEM inverse models that derive mainly from the differences in regularization between the models, with jif3D producing models that are very smooth laterally and with depth, while ModEM produces models with more discrete conductive and resistive features. Analysis of the differences between these two inversions can provide a good indication of the model resolution. More minor differences are apparent between models run with different subsets of data, with the models containing all acceptable data featuring higher wavelength conductivity variations than those run with fewer stations but also demonstrating poorer data fit.

Mantle viscosity controls a variety of geodynamic processes such as glacial isostatic adjustment (GIA), but it is poorly constrained because it cannot be measured directly from geophysical measurements. Here we develop a method that calculates viscosity using empirical viscosity flow laws coupled with mantle parameters (temperature and water content) inferred from seismic and magnetotelluric (MT) observations. We find that combining geophysical constraints allows us to place significantly tighter bounds on viscosity estimates compared to using seismic or MT observations alone. In particular, electrical conductivity inferred from MT data can determine whether upper mantle minerals are hydrated, which is important for viscosity reduction. Additionally, we show that rock composition should be considered when estimating viscosity from geophysical data because composition directly affects seismic velocity and electrical conductivity. Therefore, unknown composition increases uncertainty in temperature and water content, and makes viscosity more uncertain. Furthermore, calculations that assume pure thermal control of seismic velocity may misinterpret compositional variations as temperature, producing erroneous interpretations of mantle temperature and viscosity. Stress and grain size also affect the viscosity and its associated uncertainty, particularly via their controls on deformation regime. Dislocation creep is associated with larger viscosity uncertainties than diffusion creep. Overall, mantle viscosity can be estimated best when both seismic and MT data are available and the mantle composition, grain size and stress can be estimated. Collecting additional MT data probably offers the greatest opportunity to improve geodynamic or GIA models that rely on viscosity estimates.

The prominent seismic low-velocity zone (LVZ) in the oceanic low-viscosity asthenosphere is approximately coincident with a zone of high seismic attenuation (or low seismic Q). Small grain sizes in the asthenosphere could link these seismic and rheological properties as small grain sizes reduce viscosity and also lower seismic velocity and seismic Q. Because grain-size is reduced by rock deformation, the asthenosphere’s seismic properties can place constraints on asthenospheric deformation or flow. To determine dominant flow patterns, we develop a selfconsistent analytical 1-D channel flow model that accounts for upper mantle rheology and its dependence on flow-modified grain-sizes, water content and melt fraction, both for flow driven by surface plate motions (Couette flow) and/or by horizontal pressure gradients (Poiseuille flow). From our flow models, Couette flow dominates if the upper mantle is dry, and plug flow (a Poiseuille flow for power law rheology) if it is wet. A plug flow configuration spanning the upper 670 km of the mantle best explains the low seismic Q zone in the asthenosphere, which can be attributed to significant grain-size reduction due to extensive shearing across the asthenosphere. Below the asthenosphere, high water content and minimal shear deformation promote large grain sizes and high seismic Q. We suggest that asthenospheric low-Q and LVZ can be largely explained by grain-size variations associated with plug flow in the wet upper mantle.

The tectonic history of Southern Africa includes Archean formation of cratons, multiple episodes of subduction and rifting and some of the world’s most significant magmatic events. These processes left behind a compositional trail that can be observed in xenoliths and measured by geophysical methods. The abundance of kimberlites in southern Africa makes it an ideal place to test and calibrate mantle geophysical interpretations that can then be applied to less well-constrained regions. Magnetotellurics (MT) is a particularly useful tool for understanding tectonic history because electrical conductivity is sensitive to temperature, bulk composition, accessory minerals and rock fabric. We produced three-dimensional MT models of the southern African mantle taken from the SAMTEX MT dataset, mapped the properties of $\sim36000$ garnet xenocrysts from Group I kimberlites, and compared the results. We found that depleted regions of the mantle are uniformly associated with high electrical resistivities. The conductivity of fertile regions is more complex and depends on the specific tectonic and metasomatic history of the region, including the compositions of metasomatic fluids or melts and the emplacement of metasomatic minerals. The mantle beneath the $\sim 2.05$ Ga Bushveld Complex is highly conductive, probably caused by magmas flowing along a lithospheric weakness zone and precipitating interconnected, conductive accessory minerals such as graphite and sulfides. Kimberlites tend to be emplaced near the edges of the cratons where the mantle below 100 km depth is not highly resistive. Kimberlites avoid strong mantle conductors, suggesting a systematic relationship between their emplacement and mantle composition.