Interpretation of the Thomson Orogen and its context within the Tasmanides of eastern Australia is hampered by vast areas of deep sedimentary cover which also mask potential relationships between central and eastern Australia. Within covered areas, basement drill cores offer the only direct geological information. This study presents new detrital zircon isotopic data for these drill cores and poorly understood outcropping units to provide new age and provenance information on the Thomson Orogen. Two distinct detrital zircon signatures are revealed. One is dominated by Grenvillian-aged (1300900 Ma) zircons with a significant peak at ~1180 Ma and lesser peak at ~1070 Ma. These age peaks, along with Lu-Hf isotopic compositions (median Hf(t) = +1.5), dominantly mantle-like 18O values (median = 5.53) and model ages of ~1.89 Ga, support a Musgrave Province (central Australia) source. The dominance of Grenvillian-aged material additionally points to deposition during the Petermann Orogney (570530 Ma) when the Musgrave Province was uplifted shedding abundant material to the Centralian Superbasin. Comparable age spectra suggest that parts of the Thomson Orogen were connected to the Centralian Superbasin during this period. We use the term `Syn-Petermann to describe this signature which is observed in two drill cores adjacent to the North Australian Craton and scattered units in the outcropping Thomson Orogen. The second signature marks a significant provenance shift and is remarkably consistent throughout the Thomson Orogen. Age spectra exhibit dominant peaks at 600560 Ma, lesser 1300900 Ma populations and maximum depositional ages of ~496 Ma. This pattern is termed the `Pacific Gondwana detrital zircon signature and is recognised throughout eastern Australia, Antarctica and central Australia. LuHf isotope data for Thomson Orogen rocks with this signature is highly variable with Hf(t) values between -49 and +10 and dominantly supracrustal 18O values suggesting input from different and more diverse source regions.

This data package, completed as part of Geoscience Australia’s National Groundwater Systems (NGS) Project, presents results of the second iteration of the 3D Great Artesian Basin (GAB) and Lake Eyre Basin (LEB) (Figure 1) geological and hydrogeological models (Vizy & Rollet, 2023) populated with volume of shale (Vshale) values calculated on 2,310 wells in the Surat, Eromanga, Carpentaria and Lake Eyre basins (Norton & Rollet, 2023). This provides a refined architecture of aquifer and aquitard geometry that can be used as a proxy for internal, lateral, and vertical, variability of rock properties within each of the 18 GAB-LEB hydrogeological units (Figure 2). These data compilations and information are brought to a common national standard to help improve hydrogeological conceptualisation of groundwater systems across multiple jurisdictions. This information will assist water managers to support responsible groundwater management and secure groundwater into the future. This 3D Vshale model of the GAB provides a common framework for further data integration with other disciplines, industry, academics and the public and helps assess the impact of water use and climate change. It aids in mapping current groundwater knowledge at a GAB-wide scale and identifying critical groundwater areas for long-term monitoring. The NGS project is part of the Exploring for the Future (EFTF) program—an eight-year, $225 million Australian Government funded geoscience data and precompetitive information acquisition program. The program seeks to inform decision-making by government, community, and industry on the sustainable development of Australia's mineral, energy, and groundwater resources, including those to support the effective long-term management of GAB water resources. This work builds on the first iteration completed as part of the Great Artesian Basin Groundwater project (Vizy & Rollet, 2022; Rollet et al., 2022), and infills previous data and knowledge gaps in the GAB and LEB with additional borehole, airborne electromagnetic and seismic interpretation. The Vshale values calculated on additional wells in the southern Surat and southern Eromanga basins and in the whole of Carpentaria and Lake Eyre basins provide higher resolution facies variability estimates from the distribution of generalised sand-shale ratio across the 18 GAB-LEB hydrogeological units. The data reveals a complex mixture of sedimentary environments in the GAB, and highlights sand body development and hydraulic characteristics within aquifers and aquitards. Understanding the regional extents of these sand-rich areas provides insights into potential preferential flow paths, within and between the GAB and LEB, and aquifer compartmentalisation. However, there are limitations that require further study, including data gaps and the need to integrate petrophysics and hydrogeological data. Incorporating major faults and other structures would also enhance our understanding of fluid flow pathways. The revised Vshale model, incorporating additional boreholes to a total of 2,310 boreholes, contributes to our understanding of groundwater flow and connectivity in the region, from the recharge beds to discharge at springs, and Groundwater Dependant Ecosystems (GDEs). It also facilitates interbasinal connectivity analysis. This 3D Vshale model offers a consistent framework for integrating data from various sources, allowing for the assessment of water use impacts and climate change at different scales. It can be used to map groundwater knowledge across the GAB and identify areas that require long-term monitoring. Additionally, the distribution of boreholes with gamma ray logs used for the Vshale work in each GAB and LEB units (Norton & Rollet, 2022; 2023) is used to highlight areas where additional data acquisition or interpretation is needed in data-poor areas within the GAB and LEB units. The second iteration of surfaces with additional Vshale calculation data points provides more confidence in the distribution of sand bodies at the whole GAB scale. The current model highlights that the main Precipice, Hutton, Adori-Springbok and Cadna-owie‒Hooray aquifers are relatively well connected within their respective extents, particularly the Precipice and Hutton Sandstone aquifers and equivalents. The Bungil Formation, the Mooga Sandstone and the Gubberamunda Sandstone are partial and regional aquifers, which are restricted to the Surat Basin. These are time equivalents to the Cadna-owie–Hooray major aquifer system that extends across the Eromanga Basin, as well as the Gilbert River Formation and Eulo Queen Group which are important aquifers onshore in the Carpentaria Basin. The current iteration of the Vshale model confirms that the Cadna-owie–Hooray and time equivalent units form a major aquifer system that spreads across the whole GAB. It consists of sand bodies within multiple channel belts that have varying degrees of connectivity' i.e. being a channelised system some of the sands will be encased within overbank deposits and isolated, while others will be stacked, cross-cutting systems that provide vertical connectivity. The channelised systemtransitions vertically and laterally into a shallow marine environment (Rollet et al., 2022). Sand-rich areas are also mapped within the main Poolowanna, Brikhead-Walloon and Westbourne interbasinal aquitards, as well as the regional Rolling Downs aquitard that may provide some potential pathways for upward leakage of groundwater to the shallow Winton-Mackunda aquifer and overlying Lake Eyre Basin. Further integration with hydrochemical data may help groundtruth some of these observations. This metadata document is associated with a data package including: • Seventeen surfaces with Vshale property (Table 1), • Seventeen surfaces with less than 40% Vshale property (Table 2), • Twenty isochore with average Vshale property (Table 3), • Twenty isochore with less than 40% Vshale property (Table 4), • Sixteen Average Vshale intersections of less than 40% Vshale property delineating potential connectivity between isochore (Table 5), • Sixteen Average Vshale intersections of less than 40% Vshale property delineating potential connectivity with isochore above and below (Table 6), • Seventeen upscaled Vshale log intersection locations (Table 7), • Six regional sections showing geology and Vshale property (Table 8), • Three datasets with location of boreholes, sections, and area of interest (Table 9).

Recent national and state assessments have concluded that sedimentary formations that underlie or are within the Great Artesian Basin (GAB) may be suitable for the storage of greenhouse gases. These same formations contain methane and naturally generated carbon dioxide that has been trapped for millions of years. The Queensland government has released exploration permits for Greenhouse Gas Storage in the Bowen and Surat basins. An important consideration in assessing the potential economic, environmental, health and safety risks of such projects is the potential impact CO2 migrating out of storage reservoirs could have on overlying groundwater resources. The risk and impact of CO2 migrating from a greenhouse gas storage reservoir into groundwater cannot be objectively assessed without knowledge of the natural baseline characteristics of the groundwater within these systems. Due to the phase behaviour of CO2, geological storage of carbon dioxide in the supercritical state requires depths greater than 800m, but there are no hydrochemical studies of such deeper aquifers in the prospective storage areas. Geoscience Australia (GA) and the Geological Survey of Queensland (GSQ), Queensland Department of Mines and Energy, worked collaboratively under the National Geoscience Agreement (NGA) to characterise the regional hydrochemistry of the Denison Trough and Surat Basin and trialled different groundwater monitoring strategies. The output from this Project constitutes part of a regional baseline reference set for future site-specific and semi-regional monitoring and verification programmes conducted by geological storage proponents. The dataset provides a reference of hydrochemistry for future competing resource users.

This report presents key results of groundwater barometric response function development and interpretation from the Upper Burdekin Groundwater Project in North Queensland, conducted as part of Exploring for the Future (EFTF)—an eight year, $225 million Australian Government funded geoscience data and information acquisition program focused on better understanding the potential mineral, energy and groundwater resources across Australia. The Upper Burdekin Groundwater Project is a collaborative study between Geoscience Australia and the Queensland Government. It focuses on basalt groundwater resources in two geographically separate areas: the Nulla Basalt Province (NBP) in the south and the McBride Basalt Province (MBP) in the north. The NBP and MBP basalt aquifers are heterogeneous, fractured, vesicular systems. This report assesses how water levels in monitoring bores in the NBP and MBP respond to barometric pressure changes to evaluate the degree of formation confinement. The main process used to evaluate water level response to barometric pressure in this study is based on barometric efficiency (BE). The BE of a formation is calculated by dividing the change in monitoring bore water level by the causative barometric pressure change. Both parameters are expressed in the same units, so BE will typically be some fraction between zero and one. BE is not necessarily constant over time; the way BE changes following a theoretical step change in barometric pressure can be described using a barometric response function (BRF). BRFs were calculated in the time domain and plotted as BE against time lag for interpretation. The BRF shape was used to assess the degree of formation confinement. Although there is some uncertainty due to monitoring bore construction issues (including long effective screens) and potentially air or gas trapped in the saturated zone, all BRFs in the current project are interpreted to indicate unconfined conditions. This finding is supported by the identification of recharge at many monitoring bores through hydrograph analysis in other EFTF project components. We conclude that formations are likely to be unconfined at many project monitoring bores assessed in this study.

Geoscience Australia’s Exploring for the Future (EFTF) 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. Further detail is available at http://www.ga.gov.au/eftf. The National Groundwater Systems (NGS) project, is part of the Australian Government’s Exploring for the Future (EFTF) program, led by Geoscience Australia (https://www.eftf.ga.gov.au/national-groundwater-systems), to improve understanding of Australia’s groundwater resources to better support responsible groundwater management and secure groundwater resources into the future. The project is developing new national data coverages to constrain groundwater systems, develop a new map of Australian groundwater systems and improve data standards and workflows of groundwater assessment to populate a consistent data discovery tool and web-based mapping portal to visualise, analyse and download hydrogeological information. While our hydrogeological conceptual understanding of Australian groundwater systems continues to grow in each State and Territory jurisdiction, in addition to legacy data and knowledge from the 1970s, new information provided by recent studies in various parts of Australia highlights the level of geological complexity and spatial variability in stratigraphic and hydrostratigraphic units across the continent. We recognise the need to standardise individual datasets, such as the location and elevation of boreholes recorded in different datasets from various sources, as well as the depth and nomenclature variations of stratigraphic picks interpreted across jurisdictions to map such geological complexity in a consistent, continent-wide stratigraphic framework that can support effective long-term management of water resources and integrated resource assessments. This stratigraphic units data compilation at a continental scale forms a single point of truth for basic borehole data including 47 data sources with 1 802 798 formation picks filtered to 1 001 851 unique preferred records from 171 367 boreholes. This data compilation provides a framework to interpret various borehole datasets consistently, and can then be used in a 3D domain as an input to improve the 3D aquifer geometry and the lateral variation and connectivity in hydrostratigraphic units across Australia. The reliability of each data source is weighted to use preferentially the most confident interpretation. Stratigraphic units are standardised to the Australian Stratigraphic Units Database (ASUD) nomenclature (https://asud.ga.gov.au/search-stratigraphic-units) and assigned the corresponding ASUD code to update the information more efficiently when needed. This dataset will need to be updated as information grows and is being revised over time. This dataset provides: 1. ABSUC_v1 Australian stratigraphic unit compilation dataset (ABSUC) 2. ABSUC_v1_TOP A subset of preferred top picks from the ABSUC_v1 dataset 3. ABSUC_v1_BASE A subset of preferred base picks from the ABSUC_v1 dataset 4. ABSUC_BOREHOLE_v1 ABSUC Borehole collar dataset 5. ASUD_2023 A subset of the Australia Stratigraphic Units Database (ASUD) This consistent stratigraphic units compilation has been used to refine the Great Artesian Basin geological and hydrogeological surfaces in this region and will support the mapping of other regional groundwater systems and other resources across the continent. It can also be used to map regional geology consistently for integrated resource assessments.

<b>This data package is superseded by a second iteration presenting updates on 3D geological and hydrogeological surfaces across eastern Australia that can be accessed through </b><a href="https://dx.doi.org/10.26186/148552">https://dx.doi.org/10.26186/148552</a> The Australian Government, through the National Water Infrastructure Fund – Expansion, commissioned Geoscience Australia to undertake the project ‘Assessing the Status of Groundwater in the Great Artesian Basin’ (GAB). The project commenced in July 2019 and will finish in June 2022, with an aim to develop and evaluate new tools and techniques to assess the status of GAB groundwater systems in support of responsible management of basin water resources. While our hydrogeological conceptual understanding of the GAB continues to grow, in many places we are still reliant on legacy data and knowledge from the 1970s. Additional information provided by recent studies in various parts of the GAB highlights the level of complexity and spatial variability in hydrostratigraphic units across the basin. We now recognise the need to link these regional studies to map such geological complexity in a consistent, basin-wide hydrostratigraphic framework that can support effective long-term management of GAB water resources. Geological unit markers have been compiled and geological surfaces associated with lithostratigraphic units have been correlated across the GAB to update and refine the associated hydrogeological surfaces. Recent studies in the Surat Basin in Queensland and the Eromanga Basin in South Australia are integrated with investigations from other regions within the GAB. These bodies of work present an opportunity to link regional studies and develop a revised, internally consistent geological framework to map geological complexity across the GAB. Legacy borehole data from various sources, seismic and airborne electromagnetic (AEM) data were compiled, then combined and analysed in a common 3D domain. Correlation of interpreted geological units and stratigraphic markers from these various data sets are classified using a consistent nomenclature. This nomenclature uses geological unit subdivisions applied in the Surat Cumulative Management Area (OGIA (Office of Groundwater Impact Assessment), 2019) to correlate time equivalent regional hydrogeological units. Herein we provide an update of the surface extents and thicknesses for key hydrogeological units, reconciling geology across borders and providing the basis for a consistent hydrogeological framework at a basin-wide scale. The new surfaces can be used for facilitating an integrated basin systems assessment to improve our understanding of potential impacts from exploitation of sub-surface resources (e.g., extractive industries, agriculture and injection of large volumes of CO2 into the sub-surface) in the GAB and providing a basis for more robust water balance estimates. This report is associated with a data package including (Appendix A – Supplementary material): • Nineteen geological and hydrogeological surfaces from the Base Permo-Carboniferous, Top Permian, Base Jurassic, Base Cenozoic to the surface (Table 2.1), • Twenty-one geological and hydrogeological unit thickness maps from the top crystalline basement to the surface (Figure 3.7 to Figure 3.27), • The formation picks and constraining data points (i.e., from boreholes, seismic, AEM and outcrops) compiled and used for gridding each surface (Table 3.8).

This data release presents regional scale groundwater contours developed for the Upper Burdekin Groundwater Project in North Queensland, conducted as part of Exploring for the Future (EFTF), an Australian Government funded geoscience data and information acquisition program. The four-year (2016-20) program focused on better understanding the potential mineral, energy and groundwater resources in northern Australia. The Upper Burdekin Groundwater Project is a collaborative study between Geoscience Australia and the Queensland Government. It focuses on basalt groundwater resources in two geographically separate areas: the Nulla Basalt Province (NBP) in the south and the McBride Basalt Province (MBP) in the north. This data release includes separate, regional-scale groundwater contour datasets for the Nulla and McBride basalt provinces developed by Geoscience Australia in: Cook, S. B. & Ransley, T. R., 2020. Exploring for the Future—Groundwater level interpretations for the McBride and Nulla basalt provinces: Upper Burdekin region, North Queensland. Geoscience Australia, Canberra, https://pid.geoscience.gov.au/dataset/ga/135439. As detailed in that document, the groundwater contours were drawn by hand based on: - Groundwater levels from monitoring bores measured mostly on 17 February 2019 following extensive rainfall. - Surface topography. - Surface water features (rivers and springs). - Remote sensing data. The inferred groundwater contours were used in various Upper Burdekin Groundwater Project components to frame hydrogeological discussions. It is important to note that they were drawn following a wet period; groundwater contours are temporally variable and those presented in this data release therefore only represent part of the regional groundwater flow system.

This grid dataset is an estimation of the relative surface potential for recharge within the McBride Basalt Province. This process combined numerous factors together as to highlight the areas likely to have higher potential for recharge to occur. Soil permeability and surface geology are the primary inputs. Vegetation and slope were excluded from consideration, as these were considered to add too much complexity. Furthermore, this model does not include rainfall intensity – although this is known to vary spatially through average rainfall grids, this model is a depiction of the ground ability for recharge to occur should a significant rainfall event occur in each location. The relative surface potential recharge presented is estimated through a combination of soil and geological factors, weighting regions that are considered likely to have greater potential for recharge (e.g. younger basalts, vent-proximal facies, and highly permeable soils). Near-surface permeability of soil layers has been considered as a quantified input to the ability for water to infiltrate soil strata. It was hypothesised that locations proximal to volcanic vents would be preferential recharge sites, due to deeply penetrative columnar jointing. This suggestion is based on observations in South Iceland, where fully-penetrating columnar joint sets are more prevalent in proximal facies compared to distal facies in South Iceland (Bergh & Sigvaldson 1991). To incorporate this concept, preferential recharge sites are assumed to be within the polygons of vent-proximal facies as derived from detailed geological mapping datasets. Remaining geology has been categorised to provide higher potential recharge through younger lava flows. As such, a ranking between geological units has been used to provide the variation in potential recharge estimates. <b>References</b> Bergh, S. G., & Sigvaldason, G. E. (1991). Pleistocene mass-flow deposits of basaltic hyaloclastite on a shallow submarine shelf, South Iceland. Bulletin of Volcanology, 53(8), 597-611. doi:10.1007/bf00493688

The middle to lower Jurassic sequence in Australia's Surat Basin has been identified as a potential reservoir system for geological CO2 storage. The sequence comprises three major formations with distinctly different mineral compositions, and generally low salinity formation water (TDS<3000 mg/L). Differing geochemical responses between the formations are expected during geological CO2 storage. However, given the prevailing use of saline reservoirs in CCS projects elsewhere, limited data are available on CO2-water-rock dynamics during CO2 storage in such low-salinity formations. Here, a combined batch experiment and numerical modelling approach is used to characterise reaction pathways and to identify geochemical tracers of CO2 migration in the low-salinity Jurassic sandstone units. Reservoir system mineralogy was characterized for 66 core samples from stratigraphic well GSQ Chinchilla 4, and six representative samples were reacted with synthetic formation water and high-purity CO2 for up to 27 days at a range of pressures. Low formation water salinity, temperature, and mineralization yield high solubility trapping capacity (1.18 mol/L at 45°C, 100 bar), while the paucity of divalent cations in groundwater and the silicate reservoir matrix results in very low mineral trapping capacity under storage conditions. Formation water alkalinity buffers pH at elevated CO2 pressures and exerts control on mineral dissolution rates. Non-radiogenic, regional groundwater-like 87Sr/86Sr values (0.7048-0.7066) indicate carbonate and authigenic clay dissolution as the primary reaction pathways regulating solution composition, with limited dissolution of the clastic matrix during the incubations. Several geochemical tracers are mobilised in concentrations greater than found in regional groundwater, most notably cobalt, concentrations of which are significantly elevated regardless of CO2 pressure or sample mineralogy.

To meet the increasing demand for natural resources globally, industry faces the challenge of exploring new frontier areas that lie deeper undercover. Here, we present an approach to, and initial results of, modelling the depth of four key chronostratigraphic packages that obscure or host mineral, energy and groundwater resources. Our models are underpinned by the compilation and integration of ~200 000 estimates of the depth of these interfaces. Estimates are derived from interpretations of newly acquired airborne electromagnetic and seismic reflection data, along with boreholes, surface and solid geology, and depth to magnetic source investigations. Our curated estimates are stored in a consistent subsurface data repository. We use interpolation and machine learning algorithms to predict the distribution of these four packages away from the control points. Specifically, we focus on modelling the distribution of the base of Cenozoic-, Mesozoic-, Paleozoic- and Neoproterozoic-age stratigraphic units across an area of ~1.5 million km2 spanning the Queensland and Northern Territory border. Our repeatable and updatable approach to mapping these surfaces, together with the underlying datasets and resulting models, provides a semi-national geometric framework for resource assessment and exploration. <b>Citation:</b> Bonnardot, M.-A., Wilford, J., Rollet, N., Moushall, B., Czarnota, K., Wong, S.C.T. and Nicoll, M.G., 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.