The Upper Burdekin Basalt extents web service delivers province extents, detailed geology, spring locations and inferred regional groundwater contours for the formations of the Nulla and McBride Basalts. This work has been carried out as part of Geoscience Australia's Exploring for the Future program.

This web service provides access to groundwater raster products for the Upper Burdekin region, including: inferred relative groundwater recharge potential derived from weightings assigned to qualitative estimates of relative permeability based on mapped soil type and surface geology; Normalised Difference Vegetation Index (NDVI) used to map vegetation with potential access to groundwater in the basalt provinces, and; base surfaces of basalt inferred from sparse available data.

This report presents key results of groundwater level interpretations 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. This report interprets groundwater levels measured in both provinces by Geoscience Australia and the Queensland Government to provide recommendations for resource management. The NBP and MBP basalt aquifers are heterogeneous, fractured, vesicular systems. Several lava flows are mapped at surface in both provinces, and the degree of hydraulic connectivity between these flows is unclear. Although there was some uncertainty due to monitoring well construction issues, barometric efficiency analyses from supporting project documents suggest that the basalts of the NBP and MBP were unconfined where monitored during the EFTF project. That finding generally matches observations presented here. Longer term groundwater hydrographs suggest that groundwater levels have been declining in the NBP and MBP following major flooding in 2010-2011 related to one of the strongest La Niña events on record. Groundwater levels are yet to decline to pre-flood elevations in places. Importantly, these longer term hydrographs set the project in context: the EFTF monitoring period is only a small fraction of a much longer-functioning, dynamic groundwater system. Nulla Basalt Province The NBP is elongated east-west, and is situated entirely within the Burdekin River catchment. Volcanic vents in the west identify that area as the main extrusive centre. Regionally, groundwater migrates through the basalts of the NBP from the western high ground towards the Burdekin River in the east. Although lava flows of the NBP reach the Burdekin River, direct groundwater discharge in this area has not yet been proven. However, groundwater does discharge to various springs and surface watercourses in the NBP that are known tributaries of the Burdekin River. Despite the presence of many registered extraction bores, no clear signs of pumping were observed in groundwater hydrographs from the NBP during the EFTF monitoring period. Water levels in many bores responded to major rainfall events, ranging from a simple change in declining hydrograph slope to a water level increase of ~6.8 m in the central west. While some responses could have been induced by loading, electrical conductivity loggers and the extent of water level rise showed that many were clearly caused by recharge. At nested monitoring locations, groundwater levels remained commensurate with downward flow potentials throughout the EFTF monitoring period. McBride Basalt Province The MBP is approximately circular in plan, with volcanic vents present in a north-northeast trending band through the province centre. Lava flows extend away from the high ground of the province centre towards lower ground near the edges. In part due to its geometry, the MBP is situated within four river catchments; only surface water landing in the east flows into the Burdekin River. Regionally, groundwater migrates through the basalts of the MBP from the central high ground radially towards the edges. Direct groundwater discharge from the MBP basalts into the Burdekin River has been shown in this project. Similarly to the NBP, groundwater is also known to discharge to numerous springs and surface watercourses in the MBP. Water levels in many bores responded to major rainfall events. Responses ranged from a change in declining hydrograph slope to a water level increase of ~6.8 m in the southeast. While some responses could have been induced by loading, the extent of water level rise showed that others were clearly caused by recharge. No nested monitoring locations were installed for the EFTF project, so vertical head gradients are currently unknown. Although there are numerous registered extraction bores in the MBP, groundwater level response to pumping was only definitively identified in the east in bore RN12010016. However, several registered bores with high estimated yields have been installed in the northeast since EFTF fieldwork completion. It is possible that these higher yielding extraction bores may induce visible drawdown in monitoring bores in the future. Their high estimated yields may be associated with lava tubes; features not reported in the literature reviewed for this project for the NBP, but identified at surface and potentially in several Queensland Government bores drilled in the MBP. Conclusions and recommendations The Upper Burdekin Groundwater Project has provided abundant information on various aspects of the hydrogeology of the Nulla and McBride basalt provinces. General groundwater flow processes are understood at a regional scale for the EFTF monitoring period, but more detailed investigations and longer term monitoring are required to fully evaluate local conditions. One of the main observations of this study are the long term groundwater level declines in both the NBP and MBP following the 2010-2011 La Niña-associated floods. Groundwater levels are yet to reduce to pre-flood elevations in places, showing that the EFTF monitoring period represents only a small fraction of a much longer-functioning, dynamic groundwater system. It is unclear what, if any, contribution groundwater extraction has made to regional water level declines. Numerous correlations were assessed between groundwater hydrograph characteristics and potentially influencing factors, but the results were mostly inconclusive. There is uncertainty in hydraulic connectivity across lava flow boundaries and between intra-lava flow aquifers. Although interesting groundwater processes were identified at many bores, at the current bore spacing it is not generally possible to interpolate between locations with any certainty. Knowledge gaps and suggestions for further investigation are recorded in Section 5 of the report. The gaps identified should assist planning of future work to inform: - Further characterisation of groundwater resources. - Protection of groundwater dependent ecosystems. - Appropriate groundwater resource management.

The Tasselled Cap Wetness (TCW) percentage exceedance composite represents the behaviour of water in the landscape, as defined by the presence of water, moist soil or wet vegetation at each pixel through time. The summary shows the percentage of observed scenes where the Wetness layer of the Tasselled Cap transform is above the threshold, i.e. where each pixel has been observed as ‘wet’. Areas that retain surface water or wetness in the landscape during the dry season are potential areas of groundwater discharge and associated GDEs. The TCW exceedance composite was classified into percentage intervals to distinguish areas that were wet for different proportions of time during the 2013 dry season. Areas depicted in the dataset have been exaggerated to enable visibility.

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 technical report details the methods and results the drilling programs of the Upper Burdekin Groundwater Project 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. This report was written by Queensland Government collaborators in the Department of Environment and Science, and is published here as supplied to Geoscience Australia at the conclusion of the project. The drilling program itself was conducted by the Department of Environment and Science as part of the Upper Burdekin Groundwater Project. A total of 17 holes were drilled in 2017-18 at 13 sites with a total combined depth of 943.2 metres. These comprise selected locations across both the Nulla Basalt Province and McBride Basalt Province. A network of 15 monitoring bores were constructed with two test holes backfilled and decommissioned.

This is a raster representing the base surface of the McBride Basalt Province, inferred from sparse data available, dominated by private water bore records. This interpretation was conducted by a hydrogeologist from Geoscience Australia. Caveats • This is just one model, based on sparse data and considerable palaeotopographic interpretation • This model relies on the input datasets being accurate. However it is noted that substantial uncertainty exists both in the location of private bores and the use of drillers’ logs for identifying stratigraphic contacts. • The location of palaeothalwegs is imprecise, and often it is only indicative of the presence of a palaeovalley. • The purpose of this model is for visualisation purposes, so should not be considered a definitive depth prediction dataset.

Groundwater-dependent ecosystems (GDEs) rely on access to groundwater on a permanent or intermittent basis for some or all of their water requirements (Queensland Government, 2018). Remotely sensed data from Digital Earth Australia (DEA) (Geoscience Australia, 2018) were used to map potential aquatic and other GDEs and enhance understanding of surface water – groundwater interactions in the Upper Burdekin region. Two Landsat TM satellite products (Water Observations from Space (WOfS; Mueller et al. 2016) summary statistic and Tasselled Cap Index (TCI) wetness summary)) were used to investigate the persistence of surface water and soil moisture in the landscape to identify perennial streams, springs and other parts of the landscape that may rely on groundwater discharge. The WOfS summary statistic represents, for each pixel, the percentage of time that water is detected at the surface relative to the total number of clear observations. Due to the 25-m by 25-m pixel size of Landsat data, only features at least 25 m wide are detected and only features covering multiple pixels are consistently detected. The WOfS summary statistic was produced over the McBride and Nulla Basalt provinces for the entire period of available data (1987 to 2018). Pixels were polygonised and classified in order to visually enhance key data in the imagery, such as the identification of standing water for at least 80% of the time. The TCI is a method of reducing six surface reflectance bands of satellite data to three bands (Brightness, Greenness, Wetness) using a Principal Components Analysis (PCA) and Procrustes' Rotation (Roberts et al., 2018). The published coefficients of Crist (1985) are applied to DEA's Landsat data to generate a TCI composite. The resulting Tasselled Cap bands are a linear combination of the original surface reflectance bands that correlate with the Brightness (bare earth), Greenness and Wetness of the landscape. The TCI wetness summary (or Tasselled Cap Wetness (TCW) percentage exceedance composite), derived from the Wetness band, represents the behaviour of water in the landscape, as defined by the presence of water, moist soil or wet vegetation at each pixel through time. The summary shows the percentage of observed scenes where the Wetness layer of the Tasselled Cap transform is above the threshold, i.e. where each pixel has been observed as ‘wet’ according to the TCI. Areas that retain surface water or wetness in the landscape during the dry season are potential areas of groundwater discharge and associated GDEs. The TCW threshold is set at -600 to calculate the percentage exceedance. This threshold is based on scientific judgment and is currently in the research/testing phase. It is based on Australian conditions and conservative in nature. The dry season, when surface runoff to streams and rainfall are minimal, is particularly useful for identifying and mapping groundwater-fed streams, springs and other ecosystems that rely on access to groundwater during periods of limited rainfall. The Upper Burdekin region was especially dry between May and October 2013, with low rainfall totals in the months preceding this dry season and overall below-average rainfall conditions (i.e. decline in rainfall residual mass). The TCW exceedance composite was classified into percentage intervals to distinguish areas that were wet for different proportions of time during the 2013 dry season. Field validation of the remote sensing data products would be required to confirm the preliminary identification of parts of the landscape where groundwater discharges to the surface and potentially supports GDEs. This release includes the classified WOfS summary statistic and classified TCW percentage exceedance composite (May-October 2013) data products for the McBride and Nulla basalt provinces in the Upper Burdekin region, North Queensland. <b>References: </b> Crist EP (1985) A TM Tasseled Cap equivalent transformation for reflectance factor data. Remote Sensing of Environment 17(3), 301–306. Doi: 10.1016/0034-4257(85)90102-6. Geoscience Australia (2018) Digital Earth Australia. Geoscience Australia, http://www.ga.gov.au/dea. Mueller, N., Lewis, A., Roberts, D., Ring, S., Melrose, R., Sixsmith, J., Lymburner, L., McIntyre, A., Tan, P., Curnow, S. and Ip, A. (2016) Water observations from space: Mapping surface water from 25 years of Landsat imagery across Australia. Remote Sensing of Environment 174, 341-352, ISSN 0034-4257. Queensland Government (2018) Groundwater dependent ecosystems, WetlandInfo 2014. Queensland Government, Brisbane, https://wetlandinfo.des.qld.gov.au/wetlands/ecology/aquatic-ecosystems-natural/groundwater-dependent/. Roberts D, Dunn B and Mueller N (2018) Open Data Cube Products Using High-Dimensional Statistics of Time Series. International Geoscience and Remote Sensing Symposium. Valencia, Spain: IEEE Geoscience and Remote Sensing Society.

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 Tasselled Cap Wetness (TCW) percentage exceedance composite represents the behaviour of water in the landscape, as defined by the presence of water, moist soil or wet vegetation at each pixel through time. The summary shows the percentage of observed scenes where the Wetness layer of the Tasselled Cap transform is above the threshold, i.e. where each pixel has been observed as ‘wet’. Areas that retain surface water or wetness in the landscape during the dry season are potential areas of groundwater discharge and associated GDEs. The TCW exceedance composite was classified into percentage intervals to distinguish areas that were wet for different proportions of time during the 2013 dry season. Areas depicted in the dataset have been exaggerated to enable visibility.