The Great Artesian Basin (GAB), a hydrogeological entity that contains predominantly the Jurassic-Cretaceous Eromanga, Surat and Carpentaria geological basins, is the largest groundwater basin in Australia. It underlies one fifth of the continent, including parts of Queensland, New South Wales, South Australia and the Northern Territory. Groundwater from the GAB is a vital resource for agricultural and extractive industries, as well as for community water supply. It supports cultural values and sustains a range of groundwater-dependent ecosystems. Water managers from each jurisdiction regulate GAB resources using hydrogeological conceptualisations based on a diverse historical geoscientific nomenclature that is often unique to a jurisdiction. However, the basin and its resources are continuous across borders, and recent studies have shown high spatial variability in the hydrostratigraphic units across the basin. There is, therefore, a clear need to map the geological complexity consistently at a basin-wide scale in order to provide a hydrogeological framework to underpin effective long-term management of GAB water resources. The present study is part of the Australian Government funded project ‘Assessing the Status of Groundwater in the Great Artesian Basin’ to refine the basin conceptualisation and water balance estimates (Figure 1.1). This study focuses on an updated GAB hydrogeological architecture by compiling and standardising existing and newly interpreted biostratigraphic and well formation picks from geological logs, 2D seismic and airborne electromagnetic data in a consistent chronostratigraphic framework. This framework is used to correlate geological units across the GAB. The basin-wide correlation identifies age-equivalent sediments in different depositional settings encompassing transgressive and regressive phases. Biostratigraphic control using a common unified zonation scheme is used to identify lithological correlations. Rock properties are attributed based on sediment facies deposited during similar geological events. The approach provides a consistent way of mapping the distribution and properties of aquifers and aquitards across the GAB. In particular, the refined correlation of Jurassic and Cretaceous units between the Surat and Eromanga basins improves the resolution of hydrogeological unit geometry and lithological variation that may influence groundwater flow within and between aquifers. The 3D hydrogeological architecture developed provides a model for refining hydrogeological conceptualisations and assists in revising GAB water balance estimates. Key findings are: • The new 3D model of the GAB extends the connectivity of aquifers across the entire GAB, with potential implications for jurisdictional groundwater management. For example, the Adori Sandstone, which was previously mapped largely in the central and eastern Eromanga Basin, potentially has connectivity with the time-equivalent Springbok Sandstone in the Surat Basin across the boundary between the two basins (the Nebine Ridge). Coincident with the Nebine Ridge is a groundwater divide that tends to segregate groundwater flow between the two basins. However, cumulative impacts from excessive pumping could cause the groundwater divide to migrate due to the continuation of sandstone unit (and connectivity) across the Nebine Ridge. In addition, the Adori Sandstone is connected with the time-equivalent Algebuckina Sandstone found towards the western margin of the Eromanga Basin, which suggests there is potential for connectivity from basin margin to basin centre. This key finding improves estimates of volume and distribution of sandstone of this aquifer across all GAB jurisdictions. • The extent of other hydrogeological units have also been refined. For instance the Cadna-owie-Hooray aquifer of Ransley et al. (2015) is now separated into two units 1. Murta Formation/Hooray–Namur–Mooga sandstones aquifer and the 2. Cadna-owie–Bungil formation and equivalents aquifer. The updated mapping highlights that the upper Cadna-owie‒Bungil‒Wyandra aquifer extends across the whole GAB, and is potentially confined by the underlying Murta and lower Cadna-owie‒Bungil aquitards and overlying Rolling Downs aquitard. Higher resolution mapping of sub-units within the Cadna-owie–Bungil–Hooray and equivalents aquifer provides an improved understanding of lithological variability and the potential compartmentalisation of groundwater that may be isolated from from regional flow paths (i.e. ‘dead ends’). The lithological variability mapping within hydrogeological units highlights zones of potential connectivity where leakage may occur between the deeper and shallower aquifers, affecting upward loss of groundwater from GAB aquifers in areas distal to the outcropping recharge beds. • The new lithology mapping also highlights that the Birkhead and Westbourne formations, classified as interbedded aquitard and tight aquitard, respectively, in the Eromanga Basin, correlate laterally with time-equivalent intervals within the Algebuckina Sandstone aquifer, suggesting connection between the Hutton, Adori and Namur‒Hooray aquifers across the central and western Eromanga Basin. • The new 3D model updates hydrogeological conceptualisations in the GAB and improves groundwater balance estimates for the GAB (Ransley et al., 2022.). It is also used to constrain a regional-scale groundwater flow dynamics model for the region, including uncertainty analysis within a Bayesian framework (Knight et al., 2022). This aspect of the study is assessing a powerful approach for solving non-unique inverse problems in terms of quantifying model uncertainty. This is crucial in providing a context for, and awareness of, uncertainties in system conceptualisation that need to be accounted for, or at least acknowledged up front. • This study compiles, collates and integrates existing and newly acquired geoscientific data characterising Jurassic Cretaceous geological units that represent the hydrostratigraphy of the GAB. The updated stratigraphy improved correlations between the Eromanga, Surat and Carpentaria basins leading to better hydrogeological interpretations at the whole of GAB scale. The work draws upon the results of other recent studies to gain new insights into the geological architecture and depositional history, which have implications for groundwater occurrence and flow within and between key GAB aquifers. This updated understanding has basin-wide implications for water management, and plays a key role in revising water balance estimates for the whole GAB. The chronostratigraphic approach used here can be applied at a national scale to correlate consistently hydrostratigraphic units, providing a broader context for groundwater systems assessments.

Groundwater from the Great Artesian Basin (GAB) is a vital resource for pastoral, agricultural and extractive industries, as well as for town water supplies, supporting at least $12.8 billion in economic activity annually (Frontier Economics, 2016). It is an essential resource that supports Indigenous cultural values and sustains a range of groundwater-dependent ecosystems. The complex nature and large size of the GAB, in conjunction with increasing and competing demands for water to support new or expanding industries, communities and the environment, present complex challenges for the long-term management of the basin’s groundwater resources. Although considerable research has previously been undertaken to improve our understanding of the GAB groundwater systems, in large part by the individual jurisdictions, current knowledge varies across the basin, and there remain major hydrogeological knowledge gaps that limit management of groundwater resources in the GAB. A key challenge is to manage the groundwater resource in a way that protects existing values, accounting for inflows and outflows for the basin and ensures long-term access to artesian groundwater. The project ‘Assessing the Status of Groundwater in the Great Artesian Basin’ (referred to as the ‘Project’), was funded by the Australian Government through offsets from the former National Water Infrastructure Development Fund. Work under the Project is informing the Science Program of the National Water Grid Authority. The Project assessed existing and new geoscientific data and technologies, including satellite data, to improve our understanding of the groundwater system and water balance in the GAB, with focus areas in the northern Surat Basin (Queensland) and western Eromanga Basin (South Australia). This Project has revised and updated fundamental aspects of the GAB hydrogeological system understanding to underpin ongoing groundwater assessments and to guide water policy and resource planning in the basin. - Hydrogeological framework An updated classification of GAB aquifers and aquitards was produced, linking the hydrostratigraphic classification used in Queensland (Surat Basin) with that used in South Australia (western Eromanga Basin). This updated hydrogeological framework was produced at the whole-of-GAB scale, through the development and application of an integrated basin analysis workflow, resulting in an updated whole-of-GAB stratigraphic interpretation, consistent across jurisdictional boundaries. A total of 19 updated geological and hydrogeological surfaces were generated and used to produce a new three-dimensional hydrogeological model of the GAB. - Groundwater recharge estimation Regional groundwater recharge volumes were revised for aquifer outcrop areas along the entire eastern GAB recharge area using an improved groundwater recharge rate mapping method, integrating chloride concentration in groundwater, rainfall, soil clay content, vegetation type and surficial geology. The modelled 50th percentile (‘median’) of new groundwater recharge rate estimates, for the eastern GAB intake beds, range spatially from 0.79 mm/yr to 458.64 mm/yr, with a mean of 15 mm/yr. The modelled groundwater recharge rate output maps are subject to uncertainty. Using 1000 model replicates, the 5th percentile groundwater recharge rate map ranges from 0.06 mm/yr to 349.73 mm/yr (mean of 9 mm/yr) and the 95th percentile groundwater recharge rate map ranges from 1.19 mm/yr to 678.48 mm/yr (mean of 26 mm/yr). Groundwater recharge flow pathways into the main GAB aquifers were assessed using groundwater sample hydrochemical and environmental tracer analyses. Significant revisions to the mapped geometry and heterogeneity of the groundwater recharge beds were made using regional scale airborne electromagnetic (AEM) geophysical data, which identified the geometry of and potential connectivity between aquifers, possible structural controls on groundwater flow paths and plausible groundwater sources of spring discharge. - Groundwater system conceptualisation Revised groundwater system conceptual models of groundwater recharge processes and groundwater flow were developed. These revised groundwater system conceptualisations illustrate aquifer architecture and potential stratigraphic and/or structural variability which has the potential to affect groundwater flow paths. - Water balance estimates The water balance presented here has been undertaken to test incorporating and communicating uncertainty estimates as part of a whole-of-GAB water balance. The water balance incorporates new work, where available, from the jurisdictions and from this project (specifically groundwater recharge estimates). This project has produced quantified uncertainty for a single component of the water balance - groundwater recharge, the largest input component of the water balance. The water balance, as presented in this report, is not intended to represent a comprehensive critical appraisal of the techniques used by previous workers to estimate each element of the water balance, nor was it intended to develop new techniques for the estimation of the water balance elements other than groundwater recharge. The incorporation of a component of the uncertainty in the water balance of the GAB is new. The uncertainty in the water balance has always been there, in all past iterations, though it has not been quantified. This estimation of uncertainty is important in the communication of the water balance, as it highlights key issues such as: 1) a whole-of-basin water balance for the GAB using current information is not sufficiently detailed to be of use to water managers; and 2) local monitoring of groundwater levels and pressures remains the primary management tool for monitoring groundwater resources in the GAB. The Project produced a point-in-time assessment of the water balance of the GAB, comparing inflows (including long-term average groundwater recharge) and outflows to the main regional aquifers, for the year 2019 (The year 2019 was the latest year for which data was available at the start of the Project). The whole-of-GAB water balance, calculated using the 5th, 50th and 95th percentiles of modelled groundwater recharge rates, estimates a range of storage change volumes of -859 GL, -29 GL and +1,212 GL respectively, in 2019. The large variation in estimated storage volumes, ranging from a decreasing to increasing groundwater storage change, highlights the large uncertainty associated with the water balance when considering the groundwater recharge rate uncertainty. GAB sub-basin water balances also show a range of groundwater storage volume changes based on modelled groundwater recharge rates for the Eromanga Basin (-229 GL 5th percentile; 51 GL 50th percentile and 424 GL 95th percentile), Carpentaria Basin (-413 GL 5th percentile; -72 GL 50th percentile and 511 GL 95th percentile) and Surat Basin (-217 GL 5th percentile; -9 GL 50th percentile and 277GL 95th percentile). Using 50th percentile modelled groundwater recharge rates for major aquifer groups across the basin, water balance estimates for the Cadna-owie Aquifer Group and Precipice Aquifer Group suggest negative change in storage volumes, while water balance estimates for the Hutton - Injune Creek Aquifer Group and the Rolling Downs Aquifer Group suggest an increasing change in storage volumes. While the whole-of-GAB, sub-basin and major aquifer water balances provide basin-wide perspectives of the groundwater resource components, they also highlight the uncertainties associated with estimating groundwater recharge at a regional scale. The large range in groundwater storage values calculated for the water balance presented here, are too great to confidently provide a whole-of-GAB scale assessment of groundwater resources. - Assessment of new techniques for whole-of-GAB groundwater evaluation Assessments of new techniques, including spatial and temporal satellite data, show promising results for remote monitoring of groundwater levels at a whole-of-GAB scale. The new monitoring techniques are not currently operationalised and require further work to allow them to be integrated into current monitoring programs. The new techniques rely on ongoing groundwater level and pressure data, making existing on-ground groundwater monitoring networks essential for managing groundwater resources in the GAB. Gravity Recovery and Climate Experiment (GRACE) Gravity Recovery and Climate Experiment (GRACE) satellite derived groundwater storage change estimates were found to be largely consistent with calculated water balance groundwater storage change estimates, for GAB aquifers at the whole-of-GAB and sub-basin scale. The accuracy of GAB GRACE groundwater storage estimates is largely dependent on the accuracy of supplementary datasets required to account for gravity signals not associated with the GAB groundwater (e.g. surface water, soil moisture and shallow groundwater). The assessments undertaken for the Project indicate GRACE satellite observations are a powerful tool to remotely monitor confined groundwater storage change trends, at the whole-of-GAB and sub-basin scale over monthly to decadal time-scales. Interferometric Synthetic Aperture Radar (InSAR) Based on Interferometric Synthetic Aperture Radar (InSAR) satellite data, downward ground surface motion (subsidence) was shown to be associated with decreases in groundwater levels (drawdown) in aquifers of the northern Surat Basin focus area. In the western Eromanga Basin focus area, InSAR derived ground surface measurements were stable, consistent with stable groundwater levels over time. The Project has delivered the largest consistent mosaic of InSAR derived ground surface movement in Australia, with InSAR data processed for the majority of the GAB (~90% coverage). This assessment indicates that, where appropriate datasets for local scale corrections and interpretation are available, InSAR is a useful tool to remotely sense groundwater level changes over time. Groundwater flow model An assessment of whole-of-GAB scale groundwater flow modelling was undertaken using newly developed open-source geodynamic modelling code ‘Underworld’. The ‘Underworld’ code utilises high-performance computing capabilities and has the potential to produce groundwater flow simulations at unprecedented scale and resolution. A Bayesian-approach was applied to model simulations to characterise uncertainties with model results. The model outputs fit the hydraulic head input data acceptably and the results were consistent with the current groundwater system conceptualisation of the basin. The model has been developed as a proof-of-concept steady state model that currently has limitations. However, with further development the technique has potential to reconstruct past GAB groundwater flow regimes, testing revisions to GAB-scale hydrogeological conceptualisations and could be a useful tool to simulate basin-scale groundwater movement to complement and provide a broader context for local-scale groundwater flow models within the basin. - Data and knowledge gaps and recommendations for future work Remaining data and knowledge gaps were identified through the Project. Recommendations for future work are listed by theme below and include data acquisition, data integration and data processing at local, regional and whole-of-GAB scales to better constrain key groundwater system processes in order to improve sustainable resource management. Hydrogeological framework • Further updates to the geological framework, in particular New South Wales and Northern Territory, may be necessary to reduce lithostratigraphic interpretation uncertainty in areas where scarce palynological data combined with the presence of sandy units makes it difficult to interpret and to maintain consistency across jurisdiction boundaries. • Expand mapping of aquifer sand/shale ratios to quantify variability and connectivity within and between aquifers in targeted areas, particularly across jurisdiction boundaries. • Refine three-dimensional hydrogeological model in areas identified as having high uncertainty, such as the Carpentaria Basin. Groundwater recharge evaluation • Quantify the effects that aquifer geometry, lithological heterogeneity and structural influences have on groundwater recharge rates and pathways within the main GAB aquifers through acquisition of targeted complimentary geophysical, geological and hydrogeological data. • Rainwater sample acquisition and chemical analysis to reduce uncertainty of calculated regional groundwater recharge rates. • Acquisition of groundwater hydrochemistry and environmental tracer data to quantify key groundwater processes, including groundwater recharge rates and connectivity within and between aquifers. Groundwater system conceptualisation • Improved conceptualisation of saprolite, including hydraulic properties, extent, thickness, erosional variation, and structural disruption, to quantify the effect of saprolite on groundwater recharge processes. • Measure vertical pressure gradients and groundwater pressure-elevation profiles between GAB aquifers to conceptualise groundwater leakage, and in which direction. Water balance estimates • Standardise the methods used to estimate water balance components across the GAB and report uncertainty to increase confidence in the water balance values. • Estimation of uncertainty for all of the inflow and outflow parameters. • More exploration of how to incorporate uncertainty in the water balance calculations. Assessment of new techniques for whole-of-GAB groundwater evaluation • Develop large scale whole-of-GAB scale datasets to effectively remove non-groundwater components of GRACE signal and independently assess the GRACE estimated groundwater storage change. • Integrate InSAR with complementary datasets to correct for non-groundwater effects and identify groundwater signals at local and regional scales. • Apply the large scale groundwater flow model ensemble to estimate groundwater residence times and fluxes between aquifers as well as assess the role of hydraulic conductivity variations on groundwater flow paths with GAB aquifers. Science knowledge sharing and communication activities • Dedicated workshops with each State and Territory to share findings and identify major changes related to the new work that may have implications for water bore aquifer attribution, water level and pressure monitoring and water management. Groundwater science data and knowledge supporting the Great Artesian Basin Strategic Management Plan Australian and GAB state/territory governments have developed the Great Artesian Basin Strategic Management Plan (GAB SMP) in 2020 (DAWE, 2020). This Strategic Management Plan takes a principles-based approach to guiding governments, industry and the community in managing this important resource together to achieve economic, environmental, cultural and social outcomes for the GAB and its users. The outputs of the current project will contribute to the Great Artesian Basin Strategic Management Plan Principle 6 - Information and knowledge generation ensures that accurate, timely and readily accessible information supports good management of the Great Artesian Basin. This summary report is one of a number of products released as part of the GAB groundwater project. In addition to the key findings and outcomes presented here, companion technical reports, datasets and associated value-added products for the project are available for public use.

<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).

Geoscience Australia’s regional assessments and basin inventories are investigating Australia’s groundwater systems to improve knowledge of the nation’s groundwater potential under the Exploring for the Future (EFTF) Program and Geoscience Australia’s Strategy 2028. Where applicable, integrated basin analysis workflows are being used to build geological architecture advancing our understanding of hydrostratigraphic units and tie them to a nationally consistent chronostratigraphic framework. Here we focus on the Great Artesian Basin (GAB) and overlying Lake Eyre Basin (LEB), where groundwater is vital for pastoral, agricultural and extractive industries, community water supplies, as well as supporting indigenous cultural values and sustaining a range of groundwater dependent ecosystems such as springs and vegetation communities. Geoscience Australia continued to revise the chronostratigraphic framework and hydrostratigraphy for the GAB infilling key data and knowledge gaps from previous compilations. In collaboration with Commonwealth and State government agencies, we compiled and standardised thousands of boreholes, stratigraphic picks, 2D seismic and airborne electromagnetic data across the GAB. We undertook a detailed stratigraphic review on hundreds of key boreholes with geophysical logs to construct consistent regional transects across the GAB and LEB, using geological time constraints from hundreds of boreholes with existing and newly interpreted biostratigraphic data. We infilled the stratigraphic correlations along key transects across Queensland, New South Wales, South Australia and Northern Territory borders to refine nomenclature and stratigraphic relationships between the Surat, Eromanga and Carpentaria basins, improving chronostratigraphic understanding within the Jurassic to Cretaceous units. We extended the GAB geological framework to the overlying LEB to better resolve the Cenozoic stratigraphy and potential hydrogeological connectivity. New data and information fill gaps and refine the previous 3D hydrogeological model of the entire GAB and LEB. The new 3D geological and hydrostratigraphic model provides a framework to integrate additional hydrogeological and rock property data. It assists in refining hydraulic relationships between aquifers within the GAB and provides a basis for developing more detailed hydrogeological system conceptualisations. This is a step towards the future goal of quantifying hydraulic linkages with underlying basins, and overlying Cenozoic aquifers to underpin more robust understanding of the hydrogeological systems within the GAB. This approach can be extended to other regional hydrogeological systems. This Abstract was submitted/presented at the 2023 Australasian Exploration Geoscience Conference (AEGC) 13-18 March (https://2023.aegc.com.au/)

<div>This report summarises information regarding groundwater processes considered to have direct influence on the water balance for the Great Artesian Basin (GAB). These processes are recharge, discharge, and connectivity within the GAB sequence, as well as connectivity with underlying basins and overlying cover. </div><div>The substantial body of literature available on the GAB gives the impression that there is a considerable degree of understanding of the GAB groundwater system. This is, however, misleading. The reality is that many reports and reviews have been cited or reworked from pre-existing studies without carrying over the original uncertainties. Over time, the scale of knowledge gaps has been reduced only incrementally, while there has been a growing appreciation of the complexities in the system. With so much conceptual and quantitative uncertainty, much additional investigation is still required.</div><div><br></div>

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).

This report, completed as part of Geoscience Australia’s Exploring for the Future Program National Groundwater Systems (NGS) Project, presents results of the second iteration of 3D geological and hydrogeological surfaces across eastern Australian basins. 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. The datasets incorporate infills of data and knowledge gaps in the Great Artesian Basin (GAB), Lake Eyre Basin (LEB), Upper Darling Floodplain (UDF) and existing data in additional basins in eastern Australia. The study area extends from the offshore Gulf of Carpentaria in the north to the offshore Bight, Otway, and Gippsland basins in the South and from the western edge of the GAB in the west to the eastern Australian coastline to the east. The revisions are an update to the surface extents and thicknesses for 18 region-wide hydrogeological units produced by Vizy & Rollet, 2022. The second iteration of the 3D model surfaces further unifies geology across borders and provides the basis for a consistent hydrogeological framework at a basin-wide, and towards a national-wide, scale. The stratigraphic nomenclature used follows geological unit subdivisions applied: (1) in the Surat Cumulative Management Area (OGIA - Office of Groundwater Impact Assessment, 2019) to correlate time equivalent regional hydrogeological units in the GAB and other Jurassic and Cretaceous time equivalent basins in the study area and (2) in the LEB to correlate Cenozoic time equivalents in the study area. Triassic to Permian and older basins distribution and thicknesses are provided without any geological and hydrogeological unit sub-division. Such work helps to (1) reconcile legacy and contemporary regional studies under a common stratigraphic framework, (2) support the effective management of groundwater resources, and (3) provide a regional geological context for integrated resource assessments. The 18 hydrogeological units were constructed using legacy borehole data, 2D seismic and airborne electromagnetic (AEM) data that were compiled for the first iteration of the geological and hydrogeological surfaces under the GAB groundwater project (Vizy & Rollet, 2022a) with the addition of: • New data collected and QC’d from boreholes (including petroleum, CSG [Coal Seam Gas], stratigraphic, mineral and water boreholes) across Australia (Vizy & Rollet, 2023a) since the first iteration, including revised stratigraphic correlations filling data and knowledge gaps in the GAB, LEB, UDF region (Norton & Rollet, 2023) with revised palynological constraints (Hannaford & Rollet 2023), • Additional AEM interpretation since the first iteration in the GAB, particularly in the northern Surat (McPherson et al., 2022b), as well as in the LEB (Evans et al., in prep), in the southern Eromanga Basin (Wong et al., 2023) and in the UDF region (McPherson et al., 2022c), and • Additional 2D seismic interpretation in the Gulf of Carpentaria (Vizy & Rollet, 2023b) and in the western and central Eromanga Basin (Szczepaniak et al., 2023). These datasets were then analysed and interpreted in a common 3D domain using a consistent chronostratigraphic framework tied to the geological timescale of 2020, as defined by Hannaford et al. (2022). Confidence maps were also produced to highlight areas that need further investigation due to data gaps, in areas where better seismic depth conversion or improved well formation picks are required. New interpretations from the second iteration of the 18 surfaces include (1) new consistent and regionally continuous surfaces of Cenozoic down to Permian and older sediments beyond the extent of the GAB across eastern Australia, (2) revised extents and thicknesses of Jurassic and Cretaceous units in the GAB, including those based on distributed thickness, (3) revised extents and thicknesses of Cenozoic LEB units constrained by the underlying GAB 3D model surfaces geometry. These data constraints were not used in the model surfaces generated for the LEB detailed inventory (Evans et al., 2023), and (4) refinements of surfaces due to additional seismic and AEM interpretation used to infill data and knowledge gaps. Significant revisions include: • The use of additional seismic data to better constrain the base of the Poolowanna-Evergreen formations and equivalents and the top of Cadna-owie Formation and equivalents in the western and central Eromanga Basin, and the extent and thicknesses of the GAB units and Cenozoic Karumba Basin in the Gulf of Carpentaria, • The use of AEM interpretations to refine the geometry of outcropping units in the northern Surat Basin and the basement surface underneath the UDF region, and • A continuous 3D geological surface of base Cenozoic sediments across eastern Australia including additional constraints for the Lake Eyre Basin (borehole stratigraphy review), Murray Basin (AEM interpretation) and Karumba Basin (seismic interpretation). These revisions to the 18 geological and hydrogeological surfaces will help improve our understanding on the 3D spatial distribution of aquifers and aquitards across eastern Australia, from the groundwater recharge areas to the deep confined aquifers. These data compilations and information brought to a common national standard help improve hydrogeological conceptualisation of groundwater systems across multiple jurisdictions to assist water managers to support responsible groundwater management and secure groundwater into the future. These 3D geological and hydrogeological modelled surfaces also provide a tool for consistent data integration from multiple datasets. These modelled surfaces bring together variable data quality and coverage from different databases across state and territory jurisdictions. Data integration at various scale is important to assess potential impact of different water users and climate change. The 3D modelled surfaces can be used as a consistent framework to map current groundwater knowledge at a national scale and help highlight critical groundwater areas for long-term monitoring of potential impacts on local communities and Groundwater Dependant Ecosystems. The distribution and confidence on data points used in the current iteration of the modelled surfaces highlight where data poor areas may need further data acquisition or additional interpretation to increase confidence in the aquifers and aquitards geometry. The second iteration of surfaces highlights where further improvements can be made, notably for areas in the offshore Gulf of Carpentaria with further seismic interpretation to better constrain the base of the Aptian marine incursion (to better constrain the shape and offshore extent of the main aquifers). Inclusion of more recent studies in the offshore southern and eastern margins of Australia will improve the resolution and confidence of the surfaces, up to the edge of the Australian continental shelf. Revision of the borehole stratigraphy will need to continue where more recent data and understanding exist to improve confidence in the aquifer and aquitard geometry and provide better constraints for AEM and seismic interpretation, such as in the onshore Carpentaria, Clarence-Moreton, Sydney, Murray-Darling basins. Similarly adding new seismic and AEM interpretation recently acquired and reprocessed, such as in the eastern Eromanga Basin over the Galilee Basin, would improve confidence in the surfaces in this area. Also, additional age constraints in formations that span large periods of time would help provide greater confidence to formation sub-divisions that are time equivalent to known geological units that correlate to major aquifers and aquitards in adjacent basins, such as within the Late Jurassic‒Early Cretaceous in the Eromanga and Carpentaria basins. Finally, incorporating major faults and structures would provide greater definition of the geological and hydrogeological surfaces to inform with greater confidence fluid flow pathways in the study area. 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 1.1), • Twenty-one geological and hydrogeological unit thickness maps from the top crystalline basement to the surface (Figure 3.1 to Figure 3.21), • The formation picks and constraining data points (i.e., from boreholes, seismic, AEM and outcrops) compiled and used for gridding each surface (Table 2.7). Detailed explanation of methodology and processing is described in the associated report (Vizy & Rollet, 2023).

The GAB covers parts of New South Wales, Northern Territory, Queensland and South Australia, each with their own water management regimes. Regular assessment of the groundwater resources of the GAB at a whole-of-basin scale is important for reviewing the effectiveness and interaction of the different management approaches. However, historically this type of assessment has been undertaken on an ad hoc basis, by different agencies, for different reporting units, using different methodologies and input datasets, meaning it is difficult to compare assessments. The water balance presented here has been undertaken to test incorporating and communicating uncertainty estimates as part of a whole-of-GAB water balance. The water balance incorporates new work, where available, from the jurisdictions and from this project (specifically groundwater recharge estimates). This project has produced quantified uncertainty for a single component of the water balance - groundwater recharge, the largest input component of the water balance. The water balance, as presented in this report, is not intended to represent a comprehensive critical appraisal of the techniques used by previous workers to estimate each element of the water balance, nor was it intended to develop new techniques for the estimation of the water balance elements other than groundwater recharge. The incorporation of a component of the uncertainty in the water balance of the GAB is new. The uncertainty in the water balance has always been there, in all past iterations, though it has not been quantified. This estimation of uncertainty is important in the communication of the water balance, as it highlights several key issues: • using current information, a whole-of-basin water balance for the GAB is not sufficiently detailed to be of use to water managers • local monitoring of groundwater levels and pressures remains the primary management tool for monitoring groundwater resources in the GAB. The Project produced a point-in-time assessment of the water balance of the GAB, comparing inflows (including long-term average groundwater recharge) and outflows to the main regional aquifers, for the year 2019 (The year 2019 was the latest year for which data was available at the start of the Project). Most components of the water balance cannot be directly measured and in some cases reported values are long-term averages (e.g. groundwater recharge) and estimated values for the reporting time period. As such, the water balance relies on indirect measurements, long term averages and approximations, which have a level of inherent uncertainty resulting from the underlying assumptions used in their estimation. To calculate a whole-of GAB water balance, the project divided all the formations in the GAB into four different ‘aquifer groups’: Rolling Downs Aquifer Group; Cadna-owie Aquifer Group; Hutton – Injune Creek Aquifer Group; and Precipice Aquifer Group. This approach was developed by KCB (KCB 2016b,c,d) in Queensland and has been adopted for this Project, as it ensures all inflows and outflows from the GAB were included in the water balance. The whole-of-GAB water balance has been calculated, based on water balances estimates calculated for each constituent sub-basin – the Eromanga, Carpentaria and Surat basins. This was done to provide a picture of variations in the water balance across different parts of the GAB. The whole-of-GAB water balance, calculated using the 5th, 50th and 95th percentiles of modelled long term average groundwater recharge rates, estimates a range of storage change volumes of -859 GL, -29 GL and +1,212 GL respectively, in 2019. The large variation in estimated storage volumes, ranging from a negative to positive groundwater storage change, highlights the large uncertainty associated with the water balance. Using 50th percentile modelled groundwater recharge rates, sub-basin water balances shows an increase in groundwater storage for the Eromanga Basin (-229 GL 5th percentile; 51 GL 50th percentile and 424 GL 95th percentile) and a decrease in groundwater storage of the Carpentaria Basin (-413 GL 5th percentile; -72 GL 50th percentile and 511 GL 95th percentile) and Surat Basin (-217 GL 5th percentile; -9 GL 50th percentile and 277GL 95th percentile). Using 50th percentile modelled groundwater recharge rates for major aquifer groups across the basin, water balance estimates for the Cadna-owie Aquifer Group and Precipice Aquifer Group suggest a decrease in storage volumes, while water balance estimates for the Hutton - Injune Creek Aquifer Group and the Rolling Downs Aquifer Group suggest an increase in storage volumes. While the whole-of-GAB, sub-basin and major aquifer water balances provide basin-wide perspectives of the groundwater resource components, they also highlight the high uncertainties associated with estimating groundwater recharge at a regional scale. The large range in groundwater storage values calculated for the water balance presented here, are too great to confidently provide a whole-of-GAB scale assessment of groundwater resources. The techniques used to estimate some water balance components have improved, for example Queensland has developed a repeatable methodology for estimating unmetered groundwater extraction. This project shows that our ability to confidently model groundwater recharge to the GAB is still evolving, and will continue to improve as further investigations are undertaken. Limited hydrograph analysis showed that areas where formerly free-flowing artesian bores have been rehabilitated, through the Great Artesian Basin Sustainability Initiative and predecessor programs, are seeing stabilisation or even recovery of water levels. Monitoring bores and hydrograph analysis continue to be important elements of any water resource assessment. Undertaking the water balance for the GAB and sub-basins has been useful for highlighting components of the water balance that should be considered carefully by groundwater resource managers and where necessary, targeted for future research eg. groundwater recharge, consistent GAB wide bore discharge estimates and evapotranspiration.