<div>Lithospheric structure and composition have direct relevance for our understanding of mineral prospectivity. Aspects of the lithosphere can be imaged using geophysical inversion or analysed from exhumed samples at the surface of the Earth, but it is a challenge to ensure consistency between competing models and datasets. The LitMod platform provides a probabilistic inversion framework that uses geology as the fabric to unify multiple geophysical techniques and incorporates a priori geochemical information. Here, we present results from the application of LitMod to the Australian continent. The rasters summarise the results and performance of a Markov-chain Monte Carlo sampling from the posterior model space. Release KY22 is developed using the primary-mode Rayleigh phase velocity grids of Yoshizawa (2014).</div><div><br></div><div>Geoscience Australia's Exploring for the Future program provides precompetitive information to inform decision-making by government, community and industry on the sustainable development of Australia's mineral, energy and groundwater resources. By gathering, analysing and interpreting new and existing precompetitive geoscience data and knowledge, we are building a national picture of Australia's geology and resource potential. This leads to a strong economy, resilient society and sustainable environment for the benefit of all Australians. This includes supporting Australia's transition to a low emissions economy, strong resources and agriculture sectors, and economic opportunities and social benefits for Australia's regional and remote communities. The Exploring for the Future program, which commenced in 2016, is an eight year, $225m investment by the Australian Government.</div>

Geoscience Australia has undertaken a regional seismic mapping study of the offshore Otway Basin extending across the explored inner basin to the frontier deep-water region. Seismic interpretation covers over 18000 line-km of new and reprocessed data acquired in the 2020 Otway Basin seismic program and over 40000 line-km of legacy 2D seismic data. We present new basin-scale isochore maps that show the distribution of the Cretaceous depocentres. Maps for the Lower Cretaceous Crayfish and Eumeralla supersequences, together with those recently published for the Upper Cretaceous Shipwreck and Sherbrook Supersequences, completes the set of isochore maps for the main tectonostratigraphic basin intervals. Mapping of basement involved faults has revealed structural fabrics that have influenced depocentre development. The tectonostratigraphic development of depocentres and maps of deep crustal units delineate crustal thinning trends related to late Cretaceous extension phases. This work highlights the need to review and update structural elements. For example, the boundary between the Otway and Sorell basins is now geologically constrained. The refinements to the tectonostratigraphic evolution of the Otway Basin presented here have important implications for the distribution and potential maturity of petroleum systems, especially with regard to heat flow associated with crustal extension. Presented at the 2024 Australian Energy Producers Conference & Exhibition (AEP) (https://energyproducersconference.au/conference/)

<div>The Australian Government’s Data Driven Discoveries program has reprocessed 60 selected multi-era legacy seismic lines, covering approximately 2,520 km across the Adavale Basin, south-central Queensland. Reprocessing of legacy seismic data from the Adavale Basin aims to create a modern, consistent and integrated seismic dataset that provides new insights into the geological structure of the basin and deepens our understanding of the basin’s minerals, energy, underground storage and groundwater potential.</div><div><br></div><div>The reprocessed lines were chosen to tie into 5 wells that were previously sampled for chemostratigraphic analysis through the Data Driven Discoveries program (Riley et al., 2023, eCat 147773), including Allendale 1, Boree 1, Gilmore 1, Quilberry 1 and Stafford 1. The Adavale Basin 2D Reprocessed Seismic Data Package also complements new deep crustal seismic data being acquired in the Adavale Basin by the program.</div><div><br></div><div>The reprocessing workflow prioritised enhancing the image quality of the selected legacy seismic lines, reducing noise, and fine-tuning frequency content for specific target depths. Techniques employed included creating a 3D static model, applying noise attenuation methods, surface-consistent deconvolution, and constructing an accurate velocity model to optimise pre-stack time and depth migration. </div><div><br></div><div>Both stacks and gather data are provided in SEG-Y format, along with navigation data, velocity, and statics.</div><div><strong>&nbsp;</strong></div><div><strong>Processed gather data for this survey are available on request from clientservices@ga.gov.au - Quote eCat# 149018</strong></div>

<div>Raster datasets of inferred magnesium number for the bulk lithospheric mantle across the Australian continent. The magnesium number is an elemental ratio, defined by Mg / (Mg + Fe), which correlates to the relative enrichment or depletion in incompatible elements. Oxide concentrations are inferred in from thermo-chemical inverse modelling of Rayleigh phase velocities, surface heat flow, geoid anomalies, and topography. The magnesium number rasters summarise the results of a Markov-chain Monte Carlo sampling of the posterior model space from an ensemble of plausible candidate models. Model release 'FR23' is developed using primary-mode Rayleigh phase velocity grids adapted from Fishwick & Rawlinson (2012; "3-D structure of the Australian lithosphere from evolving seismic datasets"). Model release 'KY22' is developed using the primary-mode Rayleigh phase velocity grids of Yoshizawa (2014; "Radially anisotropic 3-D shear wave structure of the Australian lithosphere and asthenosphere from multi-mode surface waves"). All models are products of the Exploring for the Future program.</div>

The L210 South Nicholson 1096 km-long deep seismic reflection lines were acquired from 6 June to 14 August, 2017. The survey involved the acquisition of seismic reflection and gravity data along five traverses, 17GA-SN1 (375 km), 17GA-SN2 (213 km), 17GA-SN3 (58 km), 17GA-SN3 (98 km), and 17GA-SN5 (352 km). The South Nicholson seismic survey was undertaken in collaboration with and funded by: The energy theme in Geoscience Australia - Exploring for the Future; Northern Territory Geological Survey (NTGS); Department of Natural Resources and Mines - through the Geological Survey of Queensland (GSQ); and AuScope. <b>Raw data for this survey are available on request from clientservices@ga.gov.au - Quote eCat# 116881</b>

Following deep seismic reflection surveys on the Yilgarn and Pilbara cratons by Geoscience Australia with the Geological Survey of Western Australia and on the Superior Craton by the Canadian Lithoprobe program, these cratons are now some of the best surveyed Archean regions on Earth. We present seismic images that highlight how variations in crustal architecture relate to differences in Archean tectonic processes between cratons. All cratons are characterized by a mostly non-reflective 4–12 km-thick uppermost crust due to the presence of large granitoid plutons and gneissic domains. Localized regions of upper crustal seismic reflectivity are typically interpreted as supracrustal rocks and mafic sills or faults and shear zones. The middle and lower Archean crust contains variably complex geometries of relatively high amplitude reflections, though in some regions, such as the Eastern Goldfields Superterrane and the Abitibi Greenstone Belt, the lower crust appears less reflective than the middle crust. Crustal thicknesses vary from 30 km in the eastern Pilbara to 35–40 km across much of the Yilgarn and Superior, though thicknesses as great as 45–52 km occur locally in the latter two cratons. The characteristics of the Archean crust-mantle boundary, or Moho, which is commonly well-defined, differs between cratons, indicating significant variations in the tectonic processes that have driven the final stages of crustal evolution. Dipping reflections in the uppermost mantle linked to convergent crustal structures are interpreted as relict subduction scars. In the southern Superior Craton, Moho offsets and northdipping reflections in the middle and lower crust arose through successive underthrusting of Meso-Neoarchean island arcs, oceanic plateaux and microcontinental fragments, as they accreted against a pre-existing northern nucleus (e.g. North Caribou and Opatica terranes). Seismic reflection lines reveal a doubly vergent orogen above north-dipping mantle reflections that indicate subduction drive accretion. Post-orogenic crustal extension, which is inferred from crustal-scale normal shear zones and dropped greenstone belts, has not erased the original accretionary crustal architecture. In contrast, in the Yilgarn Craton interior, accretionary structures are less clear and there are no prominent offsets in the Moho. In the Youanmi Terrane, which represents the cratonic nucleus, a pervasive fabric of listric east-dipping mid-crustal reflections soles out into the upper part of subhorizontal lower crustal reflections. We interpret this reflective fabric to be the result of widespread crustal collapse during the late stage of craton evolution at c. 2.65–2.6 Ga that also produced subsidence of the upper crust. Though terrane boundaries can be identified in seismic data across the Eastern Goldfields Superterrane, these boundaries have commonly been modified by extension, which also overprinted any accretionary lower crustal structures, perhaps simultaneous with widespread intrusion of post-tectonic melts. Exhumation of moderately reflective, amphibolite to granulite facies crust in the Narryer Terrane above dipping mantle reflectors indicates that shortening along the northwestern edge of the Yilgarn Craton was subduction driven. In the eastern Pilbara Craton, shallowly dipping to subhorizontal reflections in the middle and lower crust preclude crustal-scale vertical tectonic movements and imply that the vertical displacements inferred from surface mapping were largely confined to the upper crust. <div>The abstract accompanies a talk the describes the architecture and and related tectonic processes of several Archean cratons based on reflection seismic interpretations. </div> This Abstract was submitted to & presented at the 2023 6th International Archean Symposium (6IAS) 25 - 27 July (https://6ias.org/)

<div><strong>Output Type:</strong> Exploring for the Future Extended Abstract</div><div><br></div><div><strong>Short Abstract: </strong>Knowledge of lithospheric structure is crucial information for resources exploration and deepening understanding of natural hazards. Available tomographic models of the Australian lithosphere often agree on large scale features, but in detail significant differences remain. Consequently, there is a growing need for a fully verifiable lithospheric model of Australia. Geoscience Australia has committed to develop such a model and share all results and datasets involved in model building. Here we present the first results of a full waveform inversion tomography model of Australia lithosphere down to a period of 70 s potentially able to resolve half wavelengths across continental Australia. Our model is based on seismic records from the National Seismic Network and legacy datasets with the addition of data from the currently deployed continental-scale 2° AusArray survey, which includes stations installed in previously inaccessible areas. We start with 193 earthquakes (moment magnitude (Mw) 6.2–7.5) and add 165 more earthquakes (Mw &gt;5.0) once the model progressed to a period of 70 s. Model resolution will improve over time as more data become available and more time is allowed for computation and quality control. As further iterations continue, and the inversion frequency range expands to higher frequencies, body waves can be exploited in full to constrain the model in detail and provide enough information for all components of the wavefield, building high-resolution tomographic models at a period of 40 s and below. Our model reveals previously observed first order features while revealing finer detail across much of continental Australia.</div><div><br></div><div><strong>Citation: </strong>Holzschuh, J., Gorbatov, A., Hejrani, B., Boehm, C. &amp; Hassan, R., 2024. Tomographic model of the Australian region from seismic full waveform inversion. In: Czarnota, K. (ed.) Exploring for the Future: Extended Abstracts, Geoscience Australia, Canberra, https://doi.org/10.26186/149404</div>

<div>The ‘Major crustal boundaries of Australia’ data set synthesises more than 40 years of acquisition of deep seismic reflection data across Australia, where major crustal-scale breaks, often inferred to be relict sutures between different crustal blocks, have been interpreted in the seismic reflection profiles. The widespread coverage of the seismic profiles now provides the opportunity to construct a map of major crustal boundaries across Australia. Starting with the locations of the crustal breaks identified in the seismic profiles, geological (e.g. outcrop mapping, drill hole, geochronology, isotope) and geophysical (e.g. gravity, aeromagnetic, magnetotelluric, passive seismic) data are used to map the crustal boundaries, in map view, away from the seismic profiles. For some of these boundaries, a high level of confidence can be placed on the location, whereas the location of other boundaries can only be considered to have medium or low confidence. In other areas, especially in regions covered by thick sedimentary successions, the locations of some crustal boundaries are essentially unconstrained. </div><div>The ‘Major crustal boundaries of Australia’ map shows the locations of inferred ancient plate boundaries, and will provide constraints on the three dimensional architecture of Australia. It allows a better understanding of how the Australian continent was constructed from the Mesoarchean through to the Phanerozoic, and how this evolution and these boundaries have controlled metallogenesis. It is best viewed as a dynamic dataset, which will need to be refined and updated as new information, such as new seismic reflection data, becomes available.</div><div><br></div>

<div>The architecture of the lithosphere controls the distribution of thermal, compositional and rheological interfaces. It therefore plays a fundamental role in modulating key ore-forming processes including the generation, transport, fractionation, and contamination of melts.&nbsp;Recognition of its importance has led to renewed efforts in recent years to incorporate constraints on lithospheric structure into the targeting of prospective regions for mineral exploration. One example is a suggested relationship between the genesis of porphyry copper deposits – known to be associated with evolved, silica-rich magmas – and the thickness of the crust.&nbsp;Here, using a new compilation of spot measurements, we explore the utility of crustal thickness as an exploration tool for porphyry copper deposits.</div> This Abstract was submitted & presented at the 2022 American Geophysical Union (AGU) Fall Meeting 12-16 December (https://www.agu.org/Fall-Meeting-2022)

<div><strong>Output Type: </strong>Exploring for the Future Extended Abstract</div><div><br></div><div><strong>Short Abstract: </strong>The thickness and thermal structure of continental lithosphere influences the location of seismic and volcanic hazards and is important for predicting long-term evolution of landscapes, sedimentary basins, and the distribution of natural resources. In this project, we have developed new, continental-scale models of the thermomechanical structure of the Australian plate. We begin by compiling an inventory of >15,000 geochemical analyses of peridotitic xenoliths and xenocrysts from across the continent that have been carried up to the surface in volcanic eruptions. We apply thermobarometric techniques to constrain their pressure and temperature of equilibration and perform steady-state heat flow modelling to assess the paleogeotherm beneath these sites. We subsequently use the paleogeotherms as constraints in a Bayesian calibration of anelasticity at seismic frequencies to provide a mapping between seismic velocity and temperature as a function of pressure. We apply this method to several regional-scale seismic tomography models, allowing the temperature to be continuously mapped throughout the Australian lithospheric and asthenospheric mantle. Our models include assessment of uncertainties and can be used to query thermomechanical properties, such as lithospheric thickness, heat flow through the Moho, and the Curie depth.</div><div><br></div><div><strong>Citation: </strong>Hoggard, M.J., Hazzard, J., Sudholz, Z., Richards, F., Duvernay, T., Austermann, J., Jaques, A.L., Yaxley, G., Czarnota, K. & Haynes, M., 2024. Thermochemical models of the Australian plate. In: Czarnota, K. (ed.) Exploring for the Future: Extended Abstracts, Geoscience Australia, Canberra. https://doi.org/10.26186/149411</div>