A multi-disciplinary, hydrogeological systems mapping approach has been developed to guide development of new geological and hydrogeological conceptual models, and provide a framework for understanding complex hydrogeological and hydrogeochemical processes. Integration of the 3D mapping with hydrochemical and hydrodynamic data provides critical new insights into surface-groundwater interactions and groundwater flow. Using this approach, it has been possible to develop a new understanding of recharge processes, and identify potential recharge and groundwater flow pathways. The new datasets, knowledge and hydrogeological conceptual models provide a reliable basis for the identification, characterisation and initial assessment of groundwater resources and MAR options. To meet the challenge of rapid identification and assessment of potential MAR targets and groundwater resources over the relatively large study area (7,541.5 sq km) within relatively short timeframes (18 months), the only cost-effective method with the ability to resolve key features of the hydrogeological system in the 0-150m depth range was airborne electromagnetics (AEM). The SkyTEM system is a high-resolution helicopter-borne time-domain electromagnetic system, and was developed specifically for high-resolution groundwater and environmental investigations. The SkyTEM survey, validated by borehole and ground geophysics and drilling, successfully delineated the key functional elements of the Darling Floodplain hydrogeological system, and identified potential groundwater resources, zones of river leakage, and a large number of potential MAR targets. The survey revealed significant heterogeneity in the sub-surface electrical conductivity structure, reflecting a complex geology. The survey mapped heterogeneity (and 'holes') within the near-surface aquifers and confining aquitards, while conductivity variations validated by drilling enabled five hydraulic classes (based on grain size) to be mapped within the main aquifers, as well as groundwater salinities. Locally, pump and slug tests, and NMR data were integrated with the AEM data to produce maps of interpreted hydraulic conductivity and aquifer transmissivity. Previously unrecognised faults, and landscape warping and tilting are observed to disrupt hydrostratigraphic units. These data necessitated development of a completely new hydrogeological conceptual model for the study area. This model shows the importance of faulting and erosional 'holes' in aquitards for recharge models. Discrete vertical fault offsets up to 20m produce localised inter-aquifer leakage. Sampling of rainfall, river, lake, groundwater and pore fluids has provided a comprehensive hydrochemical dataset for the alluvial aquifers of the Darling River floodplain. Major ion chemistry highlighted a mixing signature between river waters, the shallow unconfined aquifer and the underlying semi-confined target Calivil aquifer. Hydrochemical analysis including fuzzy-k means (FCM) cluster analysis, integrated with conventional hydrochemical and hydrodynamic analysis also provides invaluable new insights into groundwater processes. Recharge is dominated by river leakage during high flows, when scouring of riverbank mud veneers allows infiltration. In summary, the new hydrogeological conceptual model of the study area has enabled a number of MAR options to be identified and assessed. The integrated, multi-disciplinary approach provides critical insights for developing appropriate conceptual models for groundwater processes and dynamics. This approach provides an invaluable tool for the rapid identification and assessment of MAR options, particularly in shallow sedimentary systems. *Note: corresponding author is Ken Lawrie, as Ross S. Brodie is currently on leave until February.

In the southern half of Australia, recent droughts and predictions of a drier future under a number of climate change scenarios have led to the search for innovative strategies to identify more secure water supplies for regional communities and industries. This study was commissioned to investigate groundwater options for increasing the drought security for the city of Broken Hill. Investigations involved the assessment of a number of aquifers across a broad region, followed by the rapid mapping and assessment of potential MAR and/or groundwater extraction sites over a large data-poor area (>7,500 km2), of the Darling floodplain. Data acquisition included an airborne electromagnetics (AEM) survey (31,834 line km), a 7.5 km drilling program (100 sonic and rotary mud holes), and complementary field and laboratory measurements. Integrated studies found that surface infiltration approaches were not viable MAR options in this area due to the ubiquitous presence of thick near-surface aquitards. However, 3D mapping validated by drilling and complementary hydrogeological investigations identified >30 potential Aquifer Storage and Recovery (ASR) targets where Pliocene aquifers contain significant volumes of fresh groundwater and are sandwiched between confining aquitards. A pre-commissioning semi-quantitative residual risk assessment was carried out for a priority site (Jimargil), located within 20 km of existing power and surface water infrastructure at Menindee. Using national MAR guidelines, assessment of 12 hazard types included hydrogeological modelling, laboratory column clogging studies and geochemical assessment to identify source water treatment requirements. The study found that the residual scientific/technical risks for ASR at Jimargil are low. Key to project success was the development of new rapid mapping and assessment methodologies and workflows. It is our understanding that this is the first use of AEM as part of multi-disciplinary mapping and assessment of MAR targets. The investigations in this study also completely revised our understanding of the age, stratigraphy, structure and mode of deposition of the Darling floodplain sediments, with practical implications for the hydrogeological conceptual model underpinning the assessment of groundwater resources and MAR options.

The Mount Lofty and Flinders Ranges of South Australia are bound on the east and the west by reverse faults that thrust Proterozoic and/or Cambrian basement rocks over Quaternary sediment. These faults range from a few tens to almost one hundred kilometres in length and tend to be spaced significantly less than a fault length apart. In the few instances where the thickness of overthrust sediment can be estimated, total neotectonic throws are in the order of 100-200 m. Slip rates on individual faults range from 0.02-0.17 mm/a, with one unconfirmed estimate as high as 0.7 mm/a. Taking into account the intermittent nature of faulting in Australia, it has been suggested that 30-50% of the present-day elevation of the Flinders and Mount Lofty Ranges relative to adjacent piedmonts has developed in the last 5 Ma. Uplifted last interglacial shorelines (ca. 120 ka) along the southern coastline of the Mount Lofty Ranges indicate that deformation is ongoing. Palaeoseismological investigations provide important insight into the characteristics of the large earthquakes responsible for deformation events. Single event displacements of 1.8 m have been measured on the Williamstown-Meadows Fault and the Alma Fault, with the former relating to a surface rupture length of a least 25 km. Further to the south in Adelaide's eastern suburbs, a 5 km section of scarp, potentially relating to a single event slip on the Eden-Burnside Fault, is preserved in ca. 120 ka sediments. Where the Eden-Burnside Fault meets the coast at Port Stanvac 20 kilometres south, the last interglacial shoreline is uplifted by 2 m relative to its expected position. At Normanville, on the uplifted side of the Willunga Fault, the last interglacial shoreline is over 10 m above its expected position, implying perhaps five or more surface rupturing events in the last ca. 120 ka on this >50 km long fault. On the eastern range front, a very large single event displacement of 7 m is inferred on the 54 km long Milendella Fault, and the 79 km long Encounter Fault displaces last interglacial shorelines by up to 11 m. There is abundant evidence for large surface-breaking earthquakes on many faults within 100 km of the Adelaide CBD. Slip rates are low by plate margin standards, implying a low rate of recurrence for M7+ events on individual faults (perhaps 10,000 years or more). However, a proximal moderate-sized event or even a large event at distance has the potential to cause significant damage to Adelaide, particularly given its construction types and local site conditions.

Vertical geochemical profiling of the marine Toolebuc Formation, Eromanga Basin - implications for shale gas/oil potential The regionally extensive, marine, mid-Cretaceous (Albian) Toolebuc Formation, Eromanga Basin hosts one of Australia's most prolific potential source rocks. However, its general low thermal maturity precludes pervasive petroleum generation, although regions of high heat flow and/or deeper burial may make it attractive for unconventional (shale gas and shale oil) hydrocarbon exploration. Previous studies have provided a good understanding of the geographic distribution of the marine organic matter in the Toolebuc Formation where total organic carbon (TOC) contents range to over 20% with approx. half being of labile carbon and convertible to gas and oil. This study focuses on the vertical profiling, at the decimetre to metre scale, of the organic and inorganic geochemical fingerprints within the Toolebuc Formation with a view to quantify fluctuations in the depositional environment and mode of preservation of the organic matter and how these factors influence hydrocarbon generation thresholds. The Toolebuc Formation from three wells, Julia Creek-2 and Wallimbulla-2 and -3, was sampled over an interval from 172 to 360m depth. The total core length was 27m from which 60 samples were selected. Cores from the underlying Wallumbilla Formation (11 samples over 13m) and the overlying Allaru Mudstone (3 samples) completed the sample set. Bulk geochemical analyses included %TOC, %carbonate, %total S, -15N kerogen, -13C kerogen, -13C carbonate, -18O carbonate, and major, minor and tracer elements and quantitative mineralogy. More detailed organic geochemical analyses involved molecular fossils (saturated and aromatic hydrocarbons, and metalloporphyrins), compound specific carbon isotopes of n-alkanes, pyrolysis-gas chromatography and compositional kinetics. etc.

In this study, airborne electromagnetics (AEM), high resolution LiDAR, and drilling (100 bores) were acquired to map and assess groundwater resources and managed aquifer recharge options in the River Darling Floodplain. Neotectonic faulting and uplift has previously been described along the north-western margin of the Murray Basin along the adjacent Darling Lineament, however no evidence of neotectonics had previously been identified in the study area. Initial inversions of the AEM data revealed a multi-layered conductivity structure broadly consistent with the hydrostratigraphy identified in drilling. However, initial laterally and spatially constrained inversions showed only moderate correlations with ground data in the near-surface (~20m). As additional information from drilling and ground and borehole geophysical surveys became available, various horizontal and vertical constraints were trialled using a new Wave Number Domain Approximate Inversion procedure with a 1D multi-layer model and constraints in 3D. The resultant 3D conductivity model revealed that an important Pleistocene aquitard (Blanchetown Clay) confining the main aquifer of interest (Calivil Formation), has an undulating surface, which is locally sharply offset. An interpreted top surface suggests that it has been affected by significant warping and faulting, as well as regional tilting due to basin subsidence or margin uplift. Overall, the top surface of the Blanchetown Clay varies in elevation by 60m. Many of the sharp offsets in the conductivity layers are coincident with lineaments observed in the LiDAR data, and with underlying basement faults mapped from airborne magnetic data. The identification of neotectonics in this area was made possible through the acquisition of high resolution AEM data, and the selection of appropriate horizontal and vertical constraints in inversion procedures. Recognition of faulting in the unconsolidated sedimentary sequence helps explain the rapid recharge of underlying Pliocene aquifers, with neotectonics recognised as a key component of the hydrogeological conceptual model.

Aquifer characterisation and groundwater resources of Dili, East Timor

Sub-glacial geothermal heat flow is acknowledged to be a critical, yet poorly constrained, boundary parameter influencing ice sheet behaviour (Winsborrow et al 2010). Geothermal heat flow is the sum of residual heat from the formation of the Earth and the natural heat generated within the Earth from the radiogenic decay of the major heat producing elements (HPEs), U, Th and K. Estimates of the sub-glacial geothermal heat flow in Antarctica are largely deduced from remotely-sensed low-resolution datasets such as seismic tomography or satellite-based geomagnetics. These methods provide broad regional estimates of geothermal heat flow reflecting variations in the mantle contribution as a function of thickness of a thermally homogeneous crust. These estimates of sub-glacial geothermal heat flow, although widely utilised in ice sheet modelling studies, fail to account for lateral and vertical heterogeneity of heat production within the crust where HPEs are concentrated and that are known to significantly impact regional geothermal heat flow values. Significant variations in regional geothermal heat flow due to heterogeneous crustal distribution of HPEs have been recognised within southern Australia (e.g. McLaren et al., 2006), a region that was connected to east Antarctica along the George V, Adélie and Wilkes Lands coastline prior to breakup of Gondwana. The South Australian Heat Flow Anomaly (SAHFA; e.g. Neumann et al., 2000) is characterized by surface heat flows as high as 126 mWm-2, some '2-3 times' that of typical continental values, due to local enrichment of HPEs. The SAHFA forms part of a once contiguous continental block called the Mawson Continent, a now dismembered crustal block that is known, from geological and geophysical evidence, to extend deep into the sub-glacial interior of the Antarctic. It is highly probable that the high geothermal heat flow characteristics of the SAHFA also extend into the sub-glacial hinterland of Terra Adélie and George V lands, a possibility that has not been previously considered in ice sheet studies. In order to account for the occurrence of several sub-glacial lakes in Adélie Land, Siegert & Dowdeswell (1996) concluded that 'a further 25-50 mWm-2 of equivalent geothermal heat' was required over the assumed local geothermal heat flow of ca. 54 mWm-2. Although that study concluded that the additional heat required for basal melting was derived from internal ice deformation, they also acknowledged the possible role of variations in geothermal heat flow, and now that the SAHFA is well characterised, this is a possibility that appears very likely.

In 1994, the United Nations Regional Cartographic Conference for Asia and the Pacific resolved to establish a Permanent Committee comprising of national surveying and mapping agencies to address the concept of establishing a common geographic information infrastructure for the region. This resolution subsequently led to the establishment of the Permanent Committee for GIS Infrastructure for the Asia and Pacific (PCGIAP). One of the goals of the PCGIAP was to establish and maintain a precise understanding of the relationship between permanent geodetic stations across the region. To this end, campaign-style geodetic-GPS observations, coordinated by Geoscience Australia, have been undertaken throughout the region since 1997. In this presentation, we discuss the development of an Asia Pacific regional reference frame based on the PCGIAP GPS campaign data, which now includes data from 417 non-IGS GPS stations and provides long term crustal deformation estimates for over 200 GPS stations throughout the region. We overview and evaluate: our combination strategy with particular emphasis on the alignment of the solution onto the International Terrestrial Reference Frame (ITRF); the sensitivity of the solution to reference frame site selection; the treatment of regional co-seismic and post-seismic deformation; and the Asia-Pacific contribution to the International Association of Geodesy (IAG) Working Group on "Regional Dense Velocity Fields". The level of consistency of the coordinate estimates with respect to ITRF2005 is 6, 5, 15 mm, in the east, north and up components, respectively, while the velocity estimates are consistent at 2, 2, 6 mm/yr in the east, north and up components, respectively.

North Queensland Geodynamic and Mineral System Synthesis

A depth to magnetic basement map has been produced for the Gawler-Curnamona region of South Australia. The map combines depth to magnetic source estimates with outcrop, drill hole and seismic data. The spectral domain method of analysing the slope of straight line segments in the power spectrum was used to produce the majority of the magnetic source depth estimates. The spectral domain method was incorporated into a semi-automated in-house software package to rapidly produce the regional scale map. The reliability of the depth to magnetic basement map is heavily dependent on the reliability of the depth to magnetic source estimation methods. There are a number of factors that can lead to errors, such as data quality and wrongly assigning magnetic sources to the cover or basement. The spectral domain method tends to slightly over estimate depths, however the average absolute errors are less than %30 when compared to known depths which is considered reasonable for the production of this type of regional scale map. The map delineates large areas of prospective Gawler Craton and Curnamona Province basement beneath less than 300 m of cover material, providing a useful tool for the mineral explorer. The map also delineates large areas under thick sequences of sediments, greater than 1000 m, which may prove of interest for the hydrocarbon explorer or act as a thermal blanket for the geothermal explorer.