Seismological data are used for a variety of purposes, from earthquake hazard zonation to mapping Earth structure and mineral resource exploration. The immense volumes of seismic data now available challenge the application of routine seismic analysis techniques using existing tools. These tools fail to take advantage of recent advances in computing hardware and data formats. Given the scale of data to process and the computational complexity of algorithms involved, a more efficient approach that scales on high-performance computing and data (HPC-HPD) platforms is needed. In addition, different agencies have tended to use bespoke and ad hoc data formats, data curation processes and quality standards, hindering large-scale analyses and modelling. High-performance seismological tools (HiPerSeis) facilitate the transformation of source seismological data into formats geared towards HPC-HPD platforms. HiPerSeis also implements optimised seismological workflows that can be run at large scale on HPC-HPD platforms. <b>Citation:</b> Hassan, R., Hejrani, B., Medlin, A., Gorbatov, A. and Zhang, F., 2020. High-performance seismological tools (HiPerSeis). In: Czarnota, K., Roach, I., Abbott, S., Haynes, M., Kositcin, N., Ray, A. and Slatter, E. (eds.) Exploring for the Future: Extended Abstracts, Geoscience Australia, Canberra, 1–4.

This paper presents the application of a novel trans-dimensional sampling approach to a time domain airborne electromagnetic (AEM) inverse problem to solve for plausible conductivities of the subsurface. Geophysical inverse problems, such as time domain AEM, are well known to have a large degree of non-uniqueness. Common least-squares optimization approaches fail to take this into account and provide a single solution with linearized estimates of uncertainty which can result in overly optimistic appraisal of the conductivity of the subsurface. In this new non-linear approach, the spatial complexity of a 2D profile is controlled directly by the data. By examining an ensemble of proposed conductivity profiles it accommodates non-uniqueness and provides more robust estimates of uncertainties. <b>Citation:</b> Hawkins, R., Brodie, R. C., & Sambridge, M. (2018). Trans-dimensional Bayesian inversion of airborne electromagnetic data for 2D conductivity profiles. <i>Exploration Geophysics</i>, 49(2), 134–147. https://doi.org/10.1071/EG16139

There is a growing recognition that lithospheric structure places first-order controls on the distribution of resources within the upper crust. While this structure is increasingly imaged using geophysical techniques, there is a paucity of geological constraints on its morphology and temporal evolution. Cenozoic intraplate volcanic rocks along Australia’s eastern seaboard provide a significant opportunity to constrain mantle conditions at the time of their emplacement and thereby benchmark geophysical constraints. This volcanic activity is subdivided into two types: age-progressive provinces generated by the passage of mantle plumes beneath the plate; and age-independent provinces, which may arise from edge-driven convection at a lithospheric step. In this study, we collected and analysed 78 igneous rock samples from both types of volcanoes across Queensland. We combined these analyses with previous studies to create and augment a comprehensive database of Australian Cenozoic volcanism. Geochemical modelling techniques were used to estimate mantle temperatures and lithospheric thicknesses beneath each province. Our results show that melting occurred at depths of 45–70 km across eastern Australia. Mantle temperatures are inferred to be ~50–100 °C higher beneath age-progressive provinces than beneath age-independent provinces. These results agree with geophysical observations used to aid resource assessments and indicate that upper mantle temperatures have varied over Cenozoic times. <b>Citation:</b> Ball, P.W., Czarnota, K., White, N.J. and Champion, D.C. 2020. Exploiting Cenozoic volcanic activity to quantify upper mantle structure beneath eastern Australia. In: Czarnota, K., Roach, I., Abbott, S., Haynes, M., Kositcin, N., Ray, A. and Slatter, E. (eds.) Exploring for the Future: Extended Abstracts, Geoscience Australia, Canberra, 1–4.

The magnetotelluric (MT) method is increasingly being applied to mineral exploration under cover with several case studies showing that mineral systems can be imaged from the lower crust to the near surface. Driven by this success, the Australian Lithospheric Architecture Magnetotelluric Project (AusLAMP) is delivering long-period data on a 0.5° grid across Australia, and derived continental scale resistivity models that are helping to drive investment in mineral exploration in frontier areas. Part of this investment includes higher-resolution broadband MT surveys to enhance resolution of features of interest and improve targeting. To help gain best value for this investment it is important to have an understanding of the ability and limitations of MT to resolve features on different scales. Here we present synthetic modelling of conductive, narrow, near-vertical faults 500 m to 1500 m wide, and show that a station spacing of around 14 km across strike is sufficient to resolve these into the upper crust. However, the vertical extent of these features is not well constrained, with near-vertical planar features commonly resolved as two separate features. This highlights the need for careful interpretation of anomalies in MT inversion. In particular, in an exploration scenario, it is important to consider that a lack of interconnectivity between a lower crustal/upper mantle conductor and conductors higher up in the crust and the surface might be apparent only, and may not reflect reduced mineral prospectivity. Appeared in Exploration Geophysics Journal 05 Dec 2022

Cratonic margins host many of the natural resources upon which our society depends. Despite this, little is known about the dynamic evolution of these regions and the stability of substantial steps in plate thickness that delineate their boundaries with adjacent mantle. Here, we investigate the spatio-temporal evolution of Australian cratonic lithosphere and underlying asthenospheric mantle by using the geochemical composition of mafic volcanic or shallow intrusive rocks preserved throughout the continent’s history. We have collated a large database of mafic samples that were screened to remove data affected by crystal fractionation or assimilation of cumulate material. We use forward and inverse modelling of igneous trace element compositions to calculate the depth and extent of melting for 28 distinct igneous provinces in the North Australian Craton. These results are used to infer mantle potential temperature and lithospheric thickness at the time of eruption. The majority of Paleoproterozoic magmatic events record high mantle potential temperatures of 1350–1450 °C and relatively low lithospheric thicknesses of ≤50 km. In contrast, younger igneous provinces show a gradual decrease in potential temperature and an increase in lithospheric thickness with time. These constraints on the mantle lay the foundation for the development of a quantitative geodynamic understanding of the evolution of the Australian lithosphere and its resources. <b>Citation:</b> Klöcking, M., Czarnota, K., Champion, D.C., Jaques, A.L. and Davies, D. R., 2020. Mapping the cover in northern Australia: towards a unified national 3D geological model. In: Czarnota, K., Roach, I., Abbott, S., Haynes, M., Kositcin, N., Ray, A. and Slatter, E. (eds.) Exploring for the Future: Extended Abstracts, Geoscience Australia, Canberra, 1–4.

The product consists of 8,800 line kilometres of time‐domain airborne electromagnetic (AEM) geophysical data acquired over the far north part of South Australia known as the Musgrave Province. This product release includes: a) the measured AEM point located data, b) electrical conductivity depth images derived from the dataset, and c) the acquisition and processing report. The data were acquired using the airborne SkyTEM312 Dual Moment 275Hz/25Hz electromagnetic and magnetic system, which covered a survey area of ~14,000 km2, which includes the standard 1:250 000 map sheets of SG52-12 (Woodroffe), SG52-16 (Lindsay), SG53-09 (Alberga) and SG53-13 (Everard). The survey lines where oriented N-S and flown at 2km, 500m and 250m line spacing. A locality diagram for the survey is shown in Figure 1. This survey was funded by the Government of South Australia, as part of the Plan for Accelerating Exploration (PACE) Copper Initiative, through the Department of the Premier and Cabinet, (DPC) and the Goyder Institute of Water Research. Geoscience Australia managed the survey as part of a National Collaborative Framework project agreement with SA. The principal objective of this project was to capture a baseline geoscientific dataset to provide further information on the geological context and setting of the area for mineral systems as well as potential for groundwater resources, of the central part of the South Australian Musgrave Province. Geoscience Australia contracted SkyTEM (Australia) Pty. Ltd. to acquire SkyTEM312 electromagnetic data, between September and October 2016. The data were processed and inverted by SkyTEM using the AarhusInv inversion program (Auken et al., 2015) and the Aarhus Workbench Laterally Constrained Inversion (LCI) algorithm (Auken et al. 2005; Auken et al. 2002). The LCI code was run in multi-layer, smooth-model mode. In this mode the layer thicknesses are kept fixed and the data are inverted only for the resistivity of each layer. For this survey a 30 layer model was used. The thickness of the topmost layer was set to 2 m and the depth to the top of the bottommost (half-space) layer was set to 600 m. The layer thicknesses increase logarithmically with depth. The thicknesses and depths to the top of each layer are given in Table 1. The regional AEM survey data can be used to inform the distribution of cover sequences, and at a reconnaissance scale, trends in regolith thickness and variability, variations in bedrock conductivity, and conductivity values of key bedrock (lithology related) conductive units under cover. The data will also assist in assessing groundwater resource potential and the extent of palaeovalley systems known to exist in the Musgrave Province. A considerable area of the survey data has a small amplitude response due to resistive ground. It very soon becomes evident that lack of signal translates to erratic non-monotonic decays, quite opposite to the smooth transitional exponential decays that occur in conductive ground. Some sections of the data have been flown over what appears to be chargeable ground, hence contain what potentially can be identified as an Induced Polarization effect (airborne IP—AIP). For decades these decay sign changes, which characterize AIP, have not been accounted for in conventional AEM data processing and modelling (Viezzoli et al., 2017). Instead they have mostly been regarded as noise, calibration or levelling issues and are dealt with by smoothing, culling or applying DC shifts to the data. Not accounting for these effects is notable on the contractor’s conductivity-depth sections, where data can’t be modelled to fit the data hence large areas of blank-space have been used to substitute the conductivity structure. The selection of the survey area was undertaken through a consultative process involving the CSIRO, GOYDER Institute, Geological Survey of South Australia and the exploration companies currently active in the region (including industry survey partner PepinNini Minerals Ltd). The data will be available from Geoscience Australia’s web site free of charge. It will also be available through the South Australian Government’s SARIG website at https://map.sarig.sa.gov.au. The data will feed into the precompetitive exploration workflow developed and executed by the Geological Survey of South Australia (GSSA) and inform a new suite of value-added products directed at the exploration community.

Real advances in understanding geology for mineral, energy and groundwater resource potential of the Australian continent will come from unveiling what lies at depth, especially in the extensive under-explored regions that are obscured by cover. In this context, airborne electromagnetics (AEM) is a geophysical method at the forefront in addressing the challenge of exploration undercover. In collaboration with the state and territory geological surveys, Geoscience Australia has led a national initiative whose goal is to acquire AEM data at 20 km line spacing across Australia. This initiative, AusAEM, represents the world’s largest AEM survey flown to date; it has covered ~2.5 million km2, a substantial part of northern Australia, and is providing new insights in areas that have not been extensively explored. Regional models of subsurface electrical conductivity derived from AusAEM data support a range of applications that include geological mapping, mineral and petroleum exploration, watershed management and environmental studies. <b>Citation:</b> Ley Cooper, A.Y.. and Brodie, R.C.., 2020. AusAEM: imaging the near-surface from the world’s largest airborne electromagnetic survey. In: Czarnota, K., Roach, I., Abbott, S., Haynes, M., Kositcin, N., Ray, A. and Slatter, E. (eds.) Exploring for the Future: Extended Abstracts, Geoscience Australia, Canberra, 1–4.

Airborne electromagnetic data generated by the AusAEM Survey are shown to map mineral deposit host rocks and regional geological features within the AusAEM Survey area. We have developed new functionality in Geoscience Australia’s sample-by-sample layered earth inversion algorithm, allowing inversion of the magnitude of the combined vector sum of the X- and Z-components of TEMPEST AEM data. This functionality improves the clarity of inverted interpretation products by reducing the degree of along-line incoherency inherent to stitched 1D inversions. The new inversion approach improves the interpretability of sub-horizontal conductors, allowing better mapping of geological features under cover. Examples of geological mapping by the AusAEM survey highlight the utility of wide line spacing, regional AEM surveying to improve geological, mineral systems and groundwater resource understanding in the regions flanking outcropping mineral deposit host rocks in northern Australia. Presented at the 2019 Australasian Exploration Geoscience Conference

<div>Airborne electromagnetic (AEM) surveying provides a rapid means of imaging shallow subsurface geology as represented by changes in electrical conductivity within the earth. Aircraft-borne systems fly at different heights and with different speeds, and the exciting transients for time-domain AEM systems provide different spectral content to image the earth with. Geoscience Australia</div><div>operates a test range over one of the Menindee lakes, in New South Wales, Australia, where different AEM systems have been flown over nearly a decade. Due to well studied geology and downhole induction data available in the area, this test range provides a useful proving ground for new AEM technology. For every test survey, certain lines within the range are repeatedly flown, and high-altitude lines are also acquired, such that robust data noise statistics can be established for all overflying AEM systems.</div><div><br></div><div>Test-range data and noise for various systems naturally allows us to compare AEM derived subsurface images of the test line. This study presents the results of both deterministic as well as Bayesian stochastic inversion over the same 13 km stretch of land, with six different systems flown between 2014-2023. While a deterministic inversion provides a first-pass image for comparing AEM systems, far more information is provided by the full posterior distribution of inverted conductivities, and in particular, the marginal quantiles of median and extremal conductivities over the entire image section. </div><div><br></div><div>Our findings indicate that there is generally good agreement with borehole logs, and the posterior conductivities for all systems agree well at the regional scale. The uncertainty (or the lack thereof) around ambiguous features in deterministic inversions is revealed through the stochastic inversions. Finally, we note that examination of water volumes in Menindee lakes do not show a simple relationship with inferred conductivity, indicating that unentangling environmental factors and system differences is a non-trivial matter. Presented at the 2024 Australian Society of Exploration Geophysicists (ASEG) Discover Symposium