The Cooper Basin is Australia's premier onshore hydrocarbon producing province and hosts a range of conventional and unconventional gas play types. This study investigates the petroleum generation potential of the basin's major Permian source rocks, to improve regional understanding of the basin's hydrocarbon prospectivity. Source rock distribution, thickness, present-day amount of total organic carbon (TOC), quality (Hydrogen Index) and maturity were mapped across the basin, together with original source quality maps prior to the on-set of generation. Results of the source rock property mapping and basin-specific kinetics were integrated with 1D burial and thermal history models and a 3D basin model to create a regional pseudo-3D petroleum system model for the basin. The modelling outputs quantify the spatial distribution of both the maximum possible hydrocarbon yield, as well as the oil/ gas expelled and retained, for ten Permian source rocks. Monte Carlo simulations were used to quantify the uncertainty associated with hydrocarbon yields and to highlight the sensitivity of results to each input parameter. The principal source rocks are the Permian coal and coaly shales of the Gidgealpa Group, with highest potential yields from the Patchawarra Formation coals and coaly shales. The broad extent of the Cooper Basin's Permian source kitchen and its large total generation potential (P50 scenario >2000 bboe) highlights the basin¿s significance as a world-class hydrocarbon province. The difference between the P90 (~800 bboe) and P10 (>4000 bboe) scenarios demonstrate the range of uncertainties inherent in this modelling.

In 2017, Queensland Fire and Emergency Services (QFES) completed the State Natural Hazard Risk Assessment which evaluated the risks presented to Queensland by seven in-scope natural hazards. This publication can be found at www.disaster.qld.gov.au. The risks presented by tsunami were not evaluated as part of this assessment as there were State and Commonwealth projects underway at the time that would better inform the understanding of the hazard. These have since been completed and now underpin this guide. Following the release of the State Natural Hazard Risk Assessment and through consultation with stakeholders at all levels of Queensland’s Disaster Management Arrangements, the need for consistent information regarding Queensland’s risk from tsunami impact and inundation was identified. Accordingly, this Tsunami Guide for Queensland was developed, with support from Geoscience Australia and the Department of Environment and Science’s Coastal Impacts Unit (CIU), through a consultative process which also helped contextualise the findings of Geoscience Australia’s Probabilistic Tsunami Hazard Assessment 2018 (PTHA18) for Queensland.

<p>Geoscience Australia (GA) generated a series of gravity and magnetic grids and enhancements covering Northern Australia. Several derivative gravity datasets have been generated 1) for the North-West Shield Western Australia region (approximately between latitudes 7‒26⁰ S and longitudes 110‒130⁰ E), 2) for the Northern Territory (approximately between latitudes 7‒26⁰ S and longitudes 125.5‒141⁰ E) and for Queensland (approximately between latitudes 7‒30⁰ S and longitudes 135‒160⁰ E). The magnetic dataset has been generated only for the North-West Shield Western Australia region (approximately between latitudes 7‒26⁰ S and longitudes 110‒130⁰ E). The magnetic and gravity data were downloaded from the Geophysical Archive Data Delivery System (GADDS), website (http://www.geoscience.gov.au/cgi-bin/mapserv?map=/nas/web/ops/prod/apps/mapserver/gadds/wms_map/gadds.map&mode=browse). Satellite Free-air (FA) gravity v27.1 (released March 11, 2019) and Satellite Topography v19.1 (released January 14, 2019) data were sourced from Sandwell et al. (2014) and downloaded from the Scripps Institution of Oceanography (SIO), National Oceanic and Atmospheric Administration (NOAA), U.S. Navy and National Geospatial-Intelligence Agency (NGA) (SIO Satellite Geodesy, website, http://topex.ucsd.edu/WWW_html/mar_grav.html). The Satellite Bouguer gravity grid with onshore correction density of 2.67 gcm-3 and offshore correction density of 2.20 gcm-3 was derived from the Free-air gravity v27.1 and Topography data V19.1. This Bouguer gravity grid was used for filling areas of data gaps in the offshore region. <p>Data evaluation and processing of gravity and magnetic data available in the area of interest resulted in the production of stitched onshore-offshore Bouguer gravity grid derived from offshore satellite Bouguer gravity grid and GA’s onshore ground and airborne gravity survey data and a stitched Total Magnetic Intensity (TMI) grid derived from airborne and shipborne surveys (Tables 1 and 5). A Reduction to the Pole (RTP) grid was derived from the stitched TMI grid. The TMI, RTP, FA and terrain corrected Bouguer gravity anomalies are standard datasets for geological analysis. The free-air gravity anomaly provides the raw and basic gravity information. Images of free-air gravity are useful for first-pass interpretation and the data is used for gravity modelling. Magnetic anomalies provide information on numerous magnetic sources, including deep sources as arising from the structure and composition of magnetic basement and shallow sources such as intra-sedimentary magnetic units (e.g. volcanics, intrusions, and magnetic sedimentary layers). A standard TMI image will contain information from all these sources. Geosoft Oasis montaj software was used throughout the data processing and enhancement procedure and the montaj GridKnit module was used to generate the stitched gravity and magnetic grids. <p>Enhancement techniques have been applied to the final processed Bouguer gravity and RTP magnetic grids to highlight subtle features from various sources and to separate anomalies from different source depths. These enhancement techniques are described in the next section. <p>Enhancement processing techniques and results <p>A summary of image processing techniques used to achieve various outcomes is described in Table 1. <p>Data type Filter applied Enhancement/outcome <p>Gravity/Magnetic First vertical derivative (1VD) Near surface features (e.g. intrabasinal) <p>Gravity/Magnetic Upward continuation Noise reduction in data <p>Gravity/Magnetic Low pass filter, or large distance upward continuation Enhancement of deep features (e.g. basement) <p>Gravity/Magnetic High pass filter Enhancement of shallow features (e.g. surface anomalies) <p>Gravity/Magnetic Tilt filter and 1VD Enhancement of structure (e.g. in basement) <p>Gravity/Magnetic ZS-Edgezone and ZS-Edge filters Enhancement of edges <p>Gravity/Magnetic horizontal modulus / horizontal gradient Enhancement of boundaries <p>Magnetic RTP (reduction to the pole), Compound Anomaly, and Analytic Signal filter Accurate location of sources

A SkyTEM airborne electromagnetic (AEM) survey was flown over Surat and Galilee regions of Queensland, Australia during July 2017. The Surat-Galilee area was surveyed, with multiple survey sites in Orana, Injune and Galilee. The area is comprised of 1293 line kilometres in Orana, 552 line kilometres in Injune and 2922 line kilometres in Galilee. A total of 4767 line kilometres were flown for this survey. The projected grid coordinates have been supplied in GDA94 MGA Zone 55 for Galilee and Injune and in MGA zone 56 for Orana. The aim of the survey is to provide at a reconnaissance scale: a) trends in regolith thickness and variability b) variations in bedrock conductivity c) conductivity of key bedrock (lithology related) conductive units under cover d) the groundwater resource potential of the region e) palaeovalley systems known to exist in the region.

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.

The geology and mineral prospectivity of the southern Thomson Orogen is poorly understood because the vast majority of its extent is buried beneath younger regolith and/or sedimentary rocks. To address this issue a collaborative program to drill 16 stratigraphic boreholes was proposed to collect samples of the basement geology that can be comprehensively analysed to improve the understanding of the geological evolution of this region. To reduce the uncertainty associated with intersecting the target stratigraphy at each of the borehole sites, estimates of the cover thickness were obtained by applying the geophysical techniques of refraction seismic, audio-magnetotellurics (AMT) and targeted magnetic inversion modelling (TMIM) prior to drilling. Refraction seismic was acquired at all 16 proposed borehole sites using a system with 48 single-component geophones and a propelled weight drop primary-wave source. At 14 of the sites clear basement refractors were observed in the data. At the two other sites, Nantilla 1 and Barrygowan 1, loss of signal due to seismic attenuation at far offsets meant that a clear basement refractor was not observed. With the exception of these two sites, three distinct refractors are generally observed in the data. Those with velocities ranging from 0.4 km/s to 1.5 km/s are interpreted as regolith, those ranging from 1.8 km/s to 2.4 km/s are interpreted as Eromanga Basin sediments, and those ranging from 3.9 km/s to 5.7 km/s are interpreted as metamorphic/igneous basement. Two-dimensional velocity models of the subsurface geology were then generated using the time-term inversion method, which allowed for the thickness of each layer to be estimated. Cover thickness estimates using refraction data vary widely from site to site, with the shallowest estimate being Overshot 1 (49 m - 55 m) and the deepest Adventure Way 1 (295 m - 317 m). These variations in cover thickness estimates from site to site are indicative of basement topography variations and are not error margins. Audio-magnetotelluric data was collected at ten sites by simultaneously deploying four porous pot electrodes, to collect the two orthogonal components of telluric data (Ex and Ey), and three magnetic induction coils, to collect the three components of magnetic data (Hx, Hy and Hz). For each dataset, a one-dimensional inversion model was produced, from which resistivity contrasts were identified and used to describe electrical conductivity discontinuities in the subsurface geology. In general, the models show a near-surface conductive layer with resistivity values ≤10 Ω·m overlying layers with continuously increasing resistivities with depth (up to 102-103 Ω·m). Those layers which were >10 Ω·m were interpreted as metamorphic/igneous basement rocks and were observed to occur at depths of ~100 m to ~300 m across the survey sites, except at Overshot 1 (38 m ±10%) and Barrygowan 1 (480 m ±10%). Targeted magnetic inversion modelling (TMIM) was applied to freely available, good quality, regional airborne magnetic survey data. Depth to magnetic source estimates were generated for 53 targets, with confidence ratings, using a dipping tabular source body to model targeted magnetic anomalies in the vicinity of the borehole sites. A combined depth estimate was generated using a distance and confidence weighted average from multiple depth estimates at all but two borehole sites. Only a single depth estimate was available at Adventure Way 1 while no depth estimates were generated at Eulo 1. These combined depth estimates provide cover thickness estimates at the sites as they are likely sourced from, or near, the top of basement. Of the ten proposed borehole sites with coincident AMT and refraction seismic data, five sites have overlapping cover thickness estimates. Cover thickness estimates from the TMIM overlap both the AMT and refraction data at four sites and at two sites where only the refraction depth estimates were available. 2 Estimating Cover Thickness in the Southern Thomson Orogen The cover thickness estimates presented in this report lower the risks associated with the proposed southern Thomson Orogen stratigraphic drilling program by reducing the uncertainty in intersecting the target stratigraphy at each of the borehole sites as well as allowing for better project and program planning. Successful completion of the stratigraphic drilling program in the southern Thomson Orogen will allow for each of these geophysical methods for estimating cover thickness to be benchmarked using actual cover thicknesses measured in the boreholes.

This GA Record report is one of a series of 4 reports undertaken by the GA Groundwater Group under the National Collaboration Framework Project Agreement with the Office of Water Science (in the Department of the Environment). The report was originally submitted to the OWS in July 2013, and subsequently reviewed by Queensland government. The Maryborough Basin in Queensland is a priority coal-bearing sedimentary basin that is not currently slated for Bioregional Assessment.

This dataset includes point estimates of groundwater recharge in mm/year. Recharge rates have been estimated at monitoring bore locations in the basaltic aquifers of the Nulla and McBride basalt provinces. Recharge estimates have been calculated using the “chloride mass balance” method. The chloride mass balance process assumes that the chloride ion is a conservative tracer in precipitation, evapotranspiration, recharge and runoff; and that all the chloride is from rainfall, instead of for example halite saturation or dissolution processes. So the volumetric water balance and the flux of chloride balance must both be true. Assuming that runoff and evapotranspiration are negligible (so approximated by zero), the equation is simplified: Water balance P=ET+R+Q Water balance multiplied by chloride concentrations (chloridefluxbalance) P∙Cl_ppt=ET∙Cl_ET+R∙Cl_gw+Q∙Cl_riv | ΔCl_reac≈0 Assumptions to simplify equation P∙Cl_ppt=R∙Cl_gw | Q≈0 & ET≈0 Rearranging for recharge rate (unknown) R=P∙(Cl_ppt)/(Cl_gw ) | Q≈0 & ET≈0 Where P = precipitation rate; ET = evapotranspiration rate; R = recharge rate; Q = runoff to streams; Clppt = concentration of Cl in precipitation; ClET = concentration of chloride in evapotranspiration; Clgw = concentration of Cl in groundwater; Clriv = concentration of chloride in river runoff; ΔClreac = change in chloride concentrations from reactions.

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

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.