Building on newly acquired airborne electromagnetic and seismic reflection data during the Exploring for the Future (EFTF) program, Geoscience Australia (GA) generated a cover model across the Northern Territory and Queensland, in the Tennant Creek – Mount Isa (TISA) area (Figure 1; between 13.5 and 24.5⁰ S of latitude and 131.5 and 145⁰ E of longitude) (Bonnardot et al., 2020). The cover model provides depth estimates to chronostratigraphic layers, including: Base Cenozoic, Base Mesozoic, Base Paleozoic and Base Neoproterozoic. The depth estimates are based on the interpretation, compilation and integration of borehole, solid geology, reflection seismic, and airborne electromagnetic data, as well as depth to magnetic source estimates. These depth estimates in metres below the surface (relative to the Australian Height Datum) are consistently stored as points in the Estimates of Geophysical and Geological Surfaces (EGGS) database (Matthews et al., 2020). The data points compiled in this data package were extracted from the EGGS database. Preferred depth estimates were selected to ensure regional data consistency and aid the gridding. Two sets of cover depth surfaces (Bonnardot et al., 2020) were generated using different approaches to map megasequence boundaries associated with the Era unconformities: 1) Standard interpolation using a minimum-curvature gridding algorithm that provides minimum misfit where data points exist, and 2) Machine learning approach (Uncover-ML, Wilford et al., 2020) that allows to learn about relationships between datasets and therefore can provide better depth estimates in areas of sparse data points distribution and assess uncertainties. This data package includes the depth estimates data points compiled and used for gridding each surface, for the Base Cenozoic, Base Mesozoic, Base Paleozoic and Base Neoproterozoic (Figure 1). To provide indicative trends between the depth data points, regional interpolated depth surface grids are also provided for the Base Cenozoic, Base Mesozoic, Base Paleozoic and Base Neoproterozoic. The grids were generated with a standard interpolation algorithm, i.e. minimum-curvature interpolation method. Refined gridding method will be necessary to take into account uncertainties between the various datasets and variable distances between the points. These surfaces provide a framework to assess the depth and possible spatial extent of resources, including basin-hosted mineral resources, basement-hosted mineral resources, hydrocarbons and groundwater, as well as an input to economic models of the viability of potential resource development.

The Great Artesian Basin Research Priorities Workshop, organised by Geoscience Australia (GA), was held in Canberra on 27 and 28 April 2016. Workshop attendees represented a spectrum of stakeholders including government, policy, management, scientific and technical representatives interested in GAB-related water management. This workshop was aimed at identifying and documenting key science issues and strategies to fill hydrogeological knowledge gaps that will assist federal and state/territory governments in addressing groundwater management issues within the GAB, such as influencing the development of the next Strategic Management Plan for the GAB. This report summarises the findings out of the workshop.

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

The Upper Burdekin Chloride Mass Balance Recharge web service depicts the recharge rates have been estimated at borehole locations in the Nulla and McBride basalt provinces. Using rainfall rates, rainfall chemistry and groundwater chemistry, the recharge rates have been estimated through the Chloride Mass Balance approach.

The Severe Wind Hazard Assessment for Queensland - SWHA(Q), is a collaborative project between Geoscience Australia (GA) and Queensland Fire and Emergency Services (QFES) which aims to provide realistic and tangible information on the potential physical impacts of tropical cyclones on Queensland communities. This is intended to enable the emergency management sector and Local Governments to more effectively engage with the community on the current and future risks posed by cyclones and inform long term strategic risk management strategies.

From June 23rd to November 4th 2016 Geotech Ltd. carried out a helicopter-borne geophysical survey over part of East Isa in Queensland (figure 1). Operations were based at Cloncurry, Queensland. The traverse lines were flown in an east to west (N 90° E azimuth) direction with 2km and 2.5km traverse line spacings, with three Tie lines flown perpendicular to the traverse lines. During the survey the helicopter was maintained at a mean altitude of 76 metres above the ground with an average survey speed of 90 km/hour. This allowed for an actual average EM Transmitter-receiver loop terrain clearance of 38 metres and a magnetic sensor clearance of 68 metres. The principal geophysical sensors included a versatile time domain electromagnetic (VTEMTMPlus) full receiver-waveform system, and a caesium magnetometer. Ancillary equipment included a GPS navigation system, laser and radar altimeters, and inclinometer. A total of 15697 line-kilometres of geophysical data were acquired during the survey. The electromagnetic system is a Geotech Time Domain EM (VTEMplus) with full receiver-waveform streamed data recording at 192 kHz. The "full waveform VTEM system" uses the streamed half-cycle recording of transmitter current and receiver voltage waveforms to obtain a complete system response calibration throughout the entire survey flight. The VTEM transmitter loop and Z-component receiver coils are in a concentric-coplanar configuration and their axes are nominally vertical. An X-component receiver coil is also installed in the centre of the transmitter loop, with its axis nominally horizontal and in the flight line direction. The receiver coils measure the dB/dt response, and a B-Field response is calculated during the data processing. In-field data quality assurance and preliminary processing were carried out on a daily basis during the acquisition phase. Preliminary and final data processing, including generation of final digital data products were undertaken from the office of Geotech Ltd. in Aurora, Ontario. A set of Conductivity Depth Images (CDI) were generated using EM Flow version 3.3, developed by Encom Technologies Pty Ltd. A total of forty-five (45) dB/dt Z component channels, starting from channel 4 (21 µsec) to channel 48 (10667 µsec), were used for the CDI calculation. An averaged waveform at the receiver was used for the calculation since it was consistent for the majority of the flights with minor deviation from the average. Digital data includes all electromagnetic and magnetic data, conductivity imaging products, mulitplots plus ancillary data including the waveform.

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.

Geoscience Australia, in collaboration with state government agencies, has been collecting magnetotelluric (MT) data as part of the Australian Lithospheric Architecture Magnetotelluric Project (AusLAMP) for several years. This program aims to map the electrical resistivity of the rock layers, at depths from ten kilometres to hundreds of kilometres, across the entire continent. AusLAMP sites are each about 55 km apart from each other. Locations are chosen in consultation with landholders and other stakeholders to minimise impacts and avoid disturbance.MT data is collected using sensors that record naturally occurring variations of the Earth’s magnetic and electric fields. The equipment does not produce or transmit and signals. After four to six weeks the equipment is retrieved and the site restored to its original condition.

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