Skip to content

To report a serious incident to DPC, as required by Section 85(2) of the Petroleum and Geothermal Energy Act 2000, licensees should contact the serious incident report no: (08) 8463 6666.

High-density magnetotelluric mapping across uranium deposits highlights lithospheric controls on deposit location

Stephan Thiel1, 2, Paul Soeffky2, Lars Krieger2, Klaus Regenauer-Lieb3, Jared Peacock4, Graham Heinson2 and Kate Robertson1, 2
1 Geological Survey of South Australia, Department of the Premier and Cabinet
2 School of Physical Sciences, The University of Adelaide

3 Faculty of Engineering, University of New South Wales
4 US Geological Survey

Download this article as a PDF (4.6 MB); cite as MESA Journal 85, pages 8–11


Figure 1 MT stations (black triangles, this study; white triangles, Peacock et al. 2013) over topography. Coloured circles denote earthquake epicentres and magnitudes across the survey area although depths are poorly constrained. Reprinted from Thiel et al. (2016; fig. 1) with permission from the American Geophysical Union.
Figure 1 MT stations (black triangles, this study; white triangles, Peacock et al. 2013) over topography. Coloured circles denote earthquake epicentres and magnitudes across the survey area although depths are poorly constrained. Reprinted from Thiel et al. (2016; fig. 1) with permission from the American Geophysical Union.

With the Australian Lithospheric Architecture Magnetotelluric Project (AusLAMP) deployment across South Australia nearing completion (Thiel and Roberson 2017), the subsequent 3D inversion modelling of the magnetotelluric (MT) data and integration with other geophysical, geochemical and geological data will provide a new fundamental precompetitive dataset for the state. It will highlight areas of enhanced conductivity in the crust, which were shown by Thiel et al. (2017) to be a potential proxy for fertile crust, one characteristic for establishing mineral deposits.

However, with the site spacing of ~55 km for the AusLAMP grid, the resolution in the upper crust is not sufficient to confidently map the high conductivity zones extending into the surface cover sequence. Closely spaced broadband MT surveys are required to obtain greater resolution across the areas of interest identified by AusLAMP. There are already some examples across the state where this approach has been used to great benefit across known deposits such as Olympic Dam (e.g. Heinson, Direen and Gill 2006). This article summarises the results of a recently published closely spaced (1 km site spacing) broadband MT survey across the Beverley uranium deposit, located east of the northern Flinders Ranges (Fig. 1; Thiel et al. 2016). The results provide another example of how the crustal structure appears to exert a control on the location of mineral deposits.


Figure 2 Two-dimensional resistivity models of the MT profiles obtained using the 2D OCCAM inversion of deGroot-Hedlin and Constable (1990).
Figure 2 Two-dimensional resistivity models of the MT profiles obtained using the 2D OCCAM inversion of deGroot-Hedlin and Constable (1990). The models show increased conductivity in the lower crust at around 30 km (C2). Localised near-vertical zones of higher conductivity extend to the brittle–ductile transition (BDT) at around 10 km depth (C1). Lower conductivities extend to the surface along faults at 45° angle. Earthquake locations in the vicinity of the profile were projected along strike onto the nearest MT profile (black dots). Hypocentres are poorly constrained, and we assign a generic 15 km error. There is generally higher uncertainty in the E–W direction of epicentres due to the seismometer network geometry. Reprinted from Thiel et al. (2016; fig. 2) with permission from the American Geophysical Union.

MT is a passive electromagnetic technique which involves the recording of naturally occurring magnetic and electric fields at the surface of the earth to obtain information on the electrical resistivity structure of the subsurface. During two field campaigns, MT measurements were recorded along two 35 km, ESE–WNW-oriented parallel profiles spaced 7 km apart (Fig. 1). The northern transect comprised 36 broadband (0.004–1,000 s) stations with 13 of these also recording long-period (10–10,000 s); the southern transect comprised 59 broadband stations. An additional 22 broadband MT sites along a NNE–SSW profile were used from an array of stations centred at the Paralana enhanced geothermal systems test site (Fig. 1; Peacock et al. 2012, 2013; Thiel 2017). Site spacing was 1 km across most of the profiles, extending up to 3 km at the far east of the northern profile. The data, which is overall of good quality, is predominantly 2D based on phase tensor analysis (Caldwell, Bibby and Brown 2004) with subsequent rotation to a geoelectric strike of 5° to align with surficial faults. Data from the northern and southern lines was inverted separately using the OCCAM 2D inversion code (deGroot-Hedlin and Constable 1990; Fig. 2).


The resistivity structure revealed a 1–2 km thick layer of low resistivity (<10 Ωm) sediments above a highly resistive (>1,000 Ωm) ∼15–25 km thick layer with a predominantly 2D structure (R1), which extends along the azimuth of the geoelectric strike to a depth of ∼30 km. A subvertical mid-crustal conductor extends to about 10 km beneath the surface (C1), interpreted as a signature of fluid ponding beneath the brittle–ductile transition, which tends to resist the passage of fluids (Connolly and Podladchikov 2004). Narrow, 45°-angle conductive faults extend to the surface. The eastern mid-crust conductor beneath the Paralana enhanced geothermal system (C1) has a maximum width of about 10 km. The western region of the profiles contains very resistive (10,000 Ωm) unfractured basement rock. Along the western edge of both profiles in the mid-crust is a conductor (C4) that aligns with the range-bounding Paralana Fault (Fig. 1). The NNE–SSW profile ties in with the results of the main profiles with a lower crustal conductor and a subvertical low-conductivity extension into the mid-crust below the brittle–ductile transition (Fig. 2).

Conductivity anomalies within a few hundred kilometres of the survey area are associated with Mesoproterozoic magmatic events along the margin of the Archean–Proterozoic Gawler Craton to the west (Heinson, Direen and Gill 2006; Thiel and Heinson 2013). To the south, Delamerian-subduction associated serpentinisation and resultant enhanced magnetite within mafic and ultra-mafic rocks is interpreted to cause higher crustal conductivity (Robertson et al. 2015; Robertson et al. 2017). Low resistivity in the lower crust and upper mantle beneath the Newer Volcanic Province to the south is a result of recent decompressional melting within a lithospheric step (Aivazpourporgou et al. 2015). Recently, long-period MT AusLAMP data revealed a pervasive lower crustal conductor of the Curnamona Province and isolated conductors within the Nackara Arc province of the Ikara-Flinders Ranges (Robertson, Heinson and Thiel 2016; Robertson, Heinson and Thiel 2017). Both present and fossil crustal fluids may result in enhanced conductivity through precipitation of sulfides, graphite or magnetite on interconnected grain boundaries under the right conditions (Nover 2005; Thiel, Heinson and White 2005).

The coincident position of enhanced conductivity in the lower crust (C2; Fig. 2) and enhanced strain derived from geodynamic modelling (Célérier et al. 2005) beneath the area of maximum topographic flexure is suggestive of a connection between neotectonic uplift of the Flinders Ranges and reactivation of the lower crust beneath the survey area. The horizontal geometry of the lower crustal conductor (C2; Fig. 2) does not appear to have mantle connections; however, mantle resolution is limited due to the small profile aperture (~40 km). It is unlikely that the fluids are entirely meteoric due to the fluid flow resistant brittle–ductile transition (Connolly and Podladchikov 2004; Raimondo et al. 2013); a potential additional fluid source is geologically recent fluid release from metamorphism of lower crustal rocks.

The westernmost upper crustal pathway that branches up from C1 reaches the surface beneath the Wooltana–Poontana Fault (which has most likely been active since the Pliocene, Wülser et al. 2011) and the Beverley uranium mine (Fig. 2). The crustal architecture may therefore have some control on the location of uranium in near-surface sedimentary layers. Mobile reducing fluids may provide the required redox conditions near the neotectonic displacement along the Wooltana–Poontana Fault for uranium emplacement. The results indicate the significance of understanding the full crustal architecture when exploring for near-surface economic mineral deposits.


This study is another example of how geophysical surveys contribute to the scale-reduction process that is critical for mineral exploration success (e.g. McCuaig, Beresford and Hronsky 2010). Starting from the regional scale of the AusLAMP surveys, the follow-up high density MT surveys, such as that over the Beverley uranium mine, resolve crustal-scale faults that can control the location of mineral deposits.

As an outcome of the insights gained from this study and the AusLAMP survey, $500,000 will be invested through the PACE Copper initiative to acquire over 300 broadband MT stations across the Olympic Domain. The site layout will cover a 100 x 100 km area to map low conductivity zones to the surface in 3D. The data collection will likely commence in October/November 2017. Further in-fill surveys are planned to maximise the impact of MT methods in the process of unlocking the mineral potential of South Australia.


Petratherm provided MT data for the southern profile. Heathgate provided logistical support during the data acquisition of the northern profile. Anthony Reid (Geological Survey of South Australia) reviewed the article.


Aivazpourporgou S, Thiel S, Hayman PC, Moresi LN and Heinson G 2015. Decompression melting driving intraplate volcanism in Australia: evidence from magnetotelluric sounding. Geophysical Research Letters 42(2):346–354.

Caldwell TG, Bibby HM and Brown C 2004. The magnetotelluric phase tensor. Geophysical Journal International 158:457–469.

Célérier J, Sandiford M, Hansen DL and Quigley M 2005. Modes of active intraplate deformation, Flinders Ranges, Australia. Tectonics 24(6):TC6006–TC6006.

Connolly JAD and Podladchikov YY 2004. Fluid flow in compressive tectonic settings: implications for midcrustal seismic reflectors and downward fluid migration. Journal of Geophysical Research: Solid Earth 109(B4):B04201–B04201.

deGroot-Hedlin C and Constable S 1990. OCCAM's inversion to generate smooth, two-dimensional models from magnetotelluric data. Geophysics 55:1,613–1,624.

Heinson G, Direen N and Gill R 2006. Magnetotelluric evidence for a deep-crustal mineralizing system beneath the Olympic Dam iron oxide copper-gold deposit, southern Australia. Geology 34:573–576.

McCuaig TC, Beresford S and Hronsky J 2010. Translating the mineral systems approach into an effective exploration targeting system. Ore Geology Reviews 38:128–138.

Nover G 2005. Electrical properties of crustal and mantle rocks: a review of laboratory measurements and their explanation. Surveys in Geophysics 26:593–651.

Peacock JR, Thiel S, Reid P and Heinson G 2012. Magnetotelluric monitoring of a fluid injection: example from an enhanced geothermal system. Geophysical Research Letters 39(18):L18403–L18403.

Peacock JR, Thiel S, Heinson G and Reid P 2013. Time-lapse magnetotelluric monitoring of an enhanced geothermal system. Geophysics 78(3):B121–B130.

Raimondo T, Clark C, Hand M, Cliff J and Anczkiewicz R 2013. A simple mechanism for mid-crustal shear zones to record surface-derived fluid signatures. Geology 41(6):711–714.

Robertson K, Taylor D, Thiel S and Heinson G 2015. Magnetotelluric evidence for serpentinisation in a Cambrian subduction zone beneath the Delamerian Orogen, southeast Australia. Gondwana Research 28(2):601–611.

Robertson KE, Heinson GS, Taylor DH and Thiel S 2017. The lithospheric transition between the Delamerian and Lachlan orogens in western Victoria: new insights from 3D magnetotelluric imaging. Australian Journal of Earth Sciences 64(3):385–399.

Robertson K, Heinson G and Thiel S 2016. Lithospheric reworking at the Proterozoic-Phanerozoic transition of Australia imaged using AusLAMP magnetotelluric data. Earth and Planetary Science Letters 452:27–35.

Robertson K, Heinson G and Thiel S 2017. Mapping lithospheric alteration using AusLAMP MT data in the Ikara-Flinders Ranges and Curnamona Province. MESA Journal 83:4–7. Department of the Premier and Cabinet, Adelaide.

Thiel S 2017. Electromagnetic monitoring of hydraulic fracturing: relationship to permeability, seismicity and stress. Surveys in Geophysics.

Thiel S and Heinson G 2013. Electrical conductors in Archean mantle—result of plume interaction? Geophysical Research Letters 40:2,947–2,952.

Thiel S, Heinson G and White A 2005. Tectonic evolution of the southern Gawler Craton, South Australia, from electromagnetic sounding. Australian Journal of Earth Sciences 52:887–896.

Thiel S, Reid A, Robertson K and Heinson G 2017. Magnetotelluric characterization of cratonic lithosphere and controls on mineral deposits: examples from South Australia. IAPSO-IAMAS-IAGA Assembly, Cape Town, South Africa.

Thiel S and Robertson K 2017. Status update on AusLAMP deployment. In Department of the Premier and Cabinet, MESA Journal – News 2017, MESA Journal News 2017:42. Department of the Premier and Cabinet, South Australia, Adelaide.

Thiel S, Soeffky P, Krieger L, Regenauer-Lieb K, Peacock J and Heinson G 2016. Conductivity response to intraplate deformation: evidence for metamorphic devolatilization and crustal-scale fluid focusing. Geophysical Research Letters 43(21):11,236–11,244.

Wülser P-A, Brugger J, Foden J and Pfeifer H-R 2011. The sandstone-hosted Beverley uranium deposit, Lake Frome Basin, South Australia: mineralogy, geochemistry, and a time-constrained model for its genesis. Economic Geology 106(5):835–867.

Back to top

For more information, contact:

Stephan Thiel