This study demonstrates the importance of the conceptual hydrogeological model for the estimation of groundwater recharge rates in an alluvial system interconnected with an ephemeral or intermittent stream in south-east Queensland, Australia. The losing/gaining condition of these streams is typically subject to temporal and spatial variability, and knowledge of these hydrological processes is critical for the interpretation of recharge estimates. Recharge rate estimates of 76–182 mm/year were determined using the water budget method. The water budget method provides useful broad approximations of recharge and discharge fluxes. The chloride mass balance (CMB) method and the tritium method were used on 17 and 13 sites respectively, yielding recharge rates of 1–43 mm/year (CMB) and 4–553 mm/year (tritium method). However, the conceptual hydrogeological model confirms that the results from the CMB method at some sites are not applicable in this setting because of overland flow and channel leakage. The tritium method was appropriate here and could be applied to other alluvial systems, provided that channel leakage and diffuse infiltration of rainfall can be accurately estimated. The water-table fluctuation (WTF) method was also applied to data from 16 bores; recharge estimates ranged from 0 to 721 mm/year. The WTF method was not suitable where bank storage processes occurred.

The heterogeneous hydraulic conductivity (K) in water-bearing formations controls subsurface flow and solute transport processes. Geostatistical techniques are often employed to characterize the K distribution in space based on the correlation between K measurements. However, at the basin scale, there are often insufficient measurements for inferring the spatial correlation. This is a widespread problem that we address in this study using the example of the Betts Creek Beds (BCB) in the Galilee Basin, Australia. To address the lack of data, we use a 1D stochastic fluvial process-based model (SFPM) to quantify the total sediment thickness, Z(x), and the sandstone proportion over the total thickness, P-s(x), in the BCB. The semivariograms of Z(x) and P-s(x) are then extracted and used in sequential Gaussian simulation to construct the 2D spatial distribution of Z(x) and P-s(x). P-s(x) can be converted to a K distribution based on classical averaging methods. The results demonstrate that the combination of SFPM and geostatistical simulation allows for the evaluation of upscaled K distribution with a limited number of K measurements. Copyright (C) 2015 John Wiley & Sons, Ltd.

Coal-seam gas production requires groundwater extraction from coal-bearing formations to reduce the hydraulic pressure and improve gas recovery. In layered sedimentary basins, the coalbeds are often separated from freshwater aquifers by low-permeability aquitards. However, hydraulic connection between the coalbed and aquifers is possible due to the heterogeneity in the aquitard such as the existence of conductive faults or sandy channel deposits. For coal-seam gas extraction operations, it is desirable to identify areas in a basin where the probability of hydraulic connection between the coalbed and aquifers is low in order to avoid unnecessary loss of groundwater from aquifers and gas production problems. A connection indicator, the groundwater age indictor (GAI), is proposed, to quantify the degree of hydraulic connection. The spatial distribution of GAI can indicate the optimum positions for gas/water extraction in the coalbed. Depressurizing the coalbed at locations with a low GAI would result in little or no interaction with the aquifer when compared to the other positions. The concept of GAI is validated on synthetic cases and is then applied to the north Galilee Basin, Australia, to assess the degree of hydraulic connection between the Aramac Coal Measure and the water-bearing formations in the Great Artesian Basin, which are separated by an aquitard, the Betts Creek Beds. It is found that the GAI is higher in the western part of the basin, indicating a higher risk to depressurization of the coalbed in this region due to the strong hydraulic connection between the coalbed and the overlying aquifer.

Groundwater recharge processes, water-rock interaction and the hydraulic connectivity between aquifers of the Galilee and Eromanga Basins in central Queensland, Australia, were investigated using stable (delta H-2, delta C-18 and (87)sr/Sr-86) and radiogenic ((CI)-C-36) isotopes and dissolved methane concentrations, complemented by major ion chemistry. The central Eromanga and the upper sequence of the Galilee basins are both sub-basins of the Great Artesian Basin (GAB), and the coal seams of the Galilee Basin are currently explored for their potential as commercial coal seam gas deposits. In order to understand the potential influence of depressurisation of coal seams required to release the gas on adjacent aquifers, a detailed understanding of recharge processes and groundwater hydraulics of these basins prior to any development is required. Each of the different isotope systems were used in this study to provide different information on specific processes. For example, the assessment of delta C-13 and Sr-87/Sr-86 ratios suggested that carbonate dissolution is one of the major processes controlling the water chemistry within some aquifers. In addition, the combined assessment of delta H-2, delta O-18 and major ion chemistry indicates that transpiration is the primary process controlling the solute concentration in the GAB recharge area, whereas evaporation appears to be less significant. Groundwaters in the Galilee Basin recharge area (outside the limits of the GAB) are different to any groundwater within the GAB units. This difference is attributed to the dissolution of potassium-bearing micas, which are absent in the GAB. Groundwater age estimates based on Cl-36/Cl ratios suggest that there is a steady increase along the flow paths, and this lack of anomalous age estimates from the recharge areas to the deeper parts of the basin indicates that there is no evidence for regional inter-aquifer mixing based in isotopes only. However, dissolved methane concentrations and groundwater chemistry near faults indicates the potential mix of gas and/or water from the coal seams into the GAB groundwaters, suggesting that there may be local influence of faults as gas/water conduits. (C) 2016 Elsevier B.V. All rights reserved.

Modelling fluvial processes is an effective way to reproduce basin evolution and to recreate riverbed morphology. However, due to the complexity of alluvial environments, deterministic modelling of fluvial processes is often impossible. To address the related uncertainties, we derive a stochastic fluvial process model on the basis of the convective Exner equation that uses the statistics (mean and variance) of river velocity as input parameters. These statistics allow for quantifying the uncertainty in riverbed topography, river discharge and position of the river channel. In order to couple the velocity statistics and the fluvial process model, the perturbation method is employed with a non-stationary spectral approach to develop the Exner equation as two separate equations: the first one is the mean equation, which yields the mean sediment thickness, and the second one is the perturbation equation, which yields the variance of sediment thickness. The resulting solutions offer an effective tool to characterize alluvial aquifers resulting from fluvial processes, which allows incorporating the stochasticity of the paleoflow velocity.

A spatial analysis of hydrochemical data of groundwater and surface water was undertaken to identify groundwater-surface water connectivity in the headwaters of the Condamine River catchment, Southeast Queensland, Australia. An assessment of long-term hydrochemical and water level data supplemented by stable- and radioisotope measurements following a prolonged dry period dominated by baseflow, helped in determining patterns of interaction in different tributaries of the upper Condamine catchment. A conceptual hydrological model representing the major hydrochemical processes and their implications for stream-aquifer connectivity was developed and tested using multiple lines of evidence. The results of a multivariate statistical analysis highlight that there are two main regions with distinct hydrochemical facies (salinity, alkalinity, and predominant ions) in surface water. Geomorphology, geology, anthropogenic and climate influence were identified as the most relevant controlling factors of the spatial variability in water quality. Stable isotope data confirmed a clear evaporation trend in almost all surface water samples during baseflow conditions. Two water types can be identified and separated by the degree of evaporation and the proximity of one group to the local meteoric water line. The results confirm the discharge of groundwater from aquifers recharged by rainfall and located upstream of the surface water sampling sites. Overall, 222Rn data show a trend of increased activity in surface water towards the upstream portions of these tributaries, validating the use of this tracer to estimate groundwater input to the local creeks. The proportion of groundwater contribution to stream flow calculated by 222Rn and chloride mass balance is in agreement, and ranges between 20-70% in tributaries in the northern areas, and between 8-50% in the upper reaches of the main river channel. This study shows the efficacy of an integrated approach combining long-term hydrochemical data interpreted via multivariate statistics, hydraulic water level data and stable and radiogenic isotope hydrology for the determination of groundwater-surface interactions in headwater catchments.

Moreton Island and several other large siliceous sand dune islands and mainland barrier deposits in SE Queensland represent the distal, onshore component of an extensive Quaternary continental shelf sediment system. This sediment has been transported up to 1000 km along the coast and shelf of SE Australia over multiple glacioeustatic sea-level cycles. Stratigraphic relationships and a preliminary Optically Stimulated Luminance (OSL) chronology for Moreton Island indicate a middle Pleistocene age for the large majority of the deposit. Dune units exposed in the centre of the island and on the east coast have OSL ages that indicate deposition occurred between approximately 540 ka and 350 ka BP, and at around 96±10 ka BP. Much of the southern half of the island has a veneer of much younger sediment, with OSL ages of 0.90±0.11 ka, 1.28±0.16 ka, 5.75±0.53 ka and <0.45 ka BP. The younger deposits were partially derived from the reworking of the upper leached zone of the much older dunes. A large parabolic dune at the northern end of the island, OSL age of 9.90±1.0 ka BP, and palaeosol exposures that extend below present sea level suggest the Pleistocene dunes were sourced from shorelines positioned several to tens of metres lower than, and up to few kilometres seaward of the present shoreline. Given the lower gradient of the inner shelf a few km seaward of the island, it seems likely that periods of intermediate sea level (e.g. ~20 m below present) produced strongly positive onshore sediment budgets and the mobilisation of dunes inland to form much of what now comprises Moreton Island. The new OSL ages and comprehensive OSL chronology for the Cooloola deposit, 100 km north of Moreton Island, indicate that the bulk of the coastal dune deposits in SE Queensland were emplaced between approximately 540 ka BP and prior to the Last Interglacial. This chronostratigraphic information improves our fundamental understanding of long-term sediment transport and accumulation on large-scale continental shelf sediment systems.

The Galilee and Eromanga basins are sub-basins of the Great Artesian Basin (GAB). In this study, a multivariate statistical approach (hierarchical cluster analysis, principal component analysis and factor analysis) is carried out to identify hydrochemical patterns and assess the processes that control hydrochemical evolution within key aquifers of the GAB in these basins. The results of the hydrochemical assessment are integrated into a 3D geological model (previously developed) to support the analysis of spatial patterns of hydrochemistry, and to identify the hydrochemical and hydrological processes that control hydrochemical variability. In this area of the GAB, the hydrochemical evolution of groundwater is dominated by evapotranspiration near the recharge area resulting in a dominance of the Na-Cl water types. This is shown conceptually using two selected cross-sections which represent discrete groundwater flow paths from the recharge areas to the deeper parts of the basins. With increasing distance from the recharge area, a shift towards a dominance of carbonate (e.g. Na-HCO3 water type) has been observed. The assessment of hydrochemical changes along groundwater flow paths highlights how aquifers are separated in some areas, and how mixing between groundwater from different aquifers occurs elsewhere controlled by geological structures, including between GAB aquifers and coal bearing strata of the Galilee Basin. The results of this study suggest that distinct hydrochemical differences can be observed within the previously defined Early Cretaceous-Jurassic aquifer sequence of the GAB. A revision of the two previously recognised hydrochemical sequences is being proposed, resulting in three hydrochemical sequences based on systematic differences in hydrochemistry, salinity and dominant hydrochemical processes. The integrated approach presented in this study which combines different complementary multivariate statistical techniques with a detailed assessment of the geological framework of these sedimentary basins, can be adopted in other complex multi-aquifer systems to assess hydrochemical evolution and its geological controls.

A combination of multivariate statistical techniques, simple hydrochemical mixing models and inverse geochemical modelling was used to investigate the major hydrochemical evolutionary pathways of a large alluvial aquifer, the upper Condamine River alluvium, south-east Queensland, Australia. Hydrochemical similarities between alluvium and sedimentary bedrock groundwater imply some mixing between alluvial and sedimentary bedrock aquifers, but spatial assessment showed that this was localised around outcrops of sedimentary bedrock in upstream areas. Within the alluvium, a distinct shift towards a low salinity Na-HCO3 water type and a brackish Na-HCO3-Cl water type was obvious in two separate locations. Both of these water types are unique to the alluvium, and inverse modelling shows that they can evolve via a combination of in situ alluvial processes, including diffuse recharge of rainfall or river water or the evolution of basalt-derived groundwater via gypsum dissolution plagioclase weathering, cation exchange and some carbonate precipitation/dissolution. The evolution of these water types is potentially influenced by overlying sodic alkaline soils, and often is associated with a source of sulfate. Evapotranspiration is the dominant salinization process in the alluvium and increases in calcium cations during salinization indicate that brackish Na-HCO3-Cl groundwater in the underlying Walloon Coal Measures are unlikely to have a major influence on salinization in the alluvium. The most saline water types observed were endemic to shallow zones of the alluvium where evapotranspiration is likely. Results demonstrate that a combination of multivariate statistics and inverse geochemical modelling can be successfully used to delineate hydrochemical pathways in complex hydrogeological settings where a range of environmental and anthropogenic factors may be influencing the evolution of water types with similar hydrochemical compositions.

Using a combination of multivariate statistical techniques and the graphical assessment of major ion ratios, the influences on hydrochemical variability of coal seam gas (or coal bed methane) groundwaters from several sites in the Surat and Clarence-Moreton basins in Queensland, Australia, were investigated. Several characteristic relationships between major ions were observed: 1) strong positive linear correlation between the Na/Cl and alkalinity/Cl ratios; 2) an exponentially decaying trend between the Na/Cl and Na/alkalinity ratios; 3) inverse linear relationships between increasing chloride concentrations and decreasing pH for high salinity groundwaters; and 4) high residual alkalinity for lower salinity waters, and an inverse relationship between decreasing residual alkalinity and increasing chloride concentrations for more saline waters. The interpretation of the hydrochemical data provides invaluable insights into the hydrochemical evolution of coal seam gas (CSG) groundwaters that considers both the source of major ions in coals and the influence of microbial activity. Elevated chloride and sodium concentrations in more saline groundwaters appear to be influenced by organic-bound chlorine held in the coal matrix; a sodium and chloride ion source that has largely been neglected in previous CSG groundwater studies. However, contrastingly high concentrations of bicarbonate in low salinity waters could not be explained, and are possibly associated with a number of different factors such as coal degradation, methanogenic processes, the evolution of high-bicarbonate NaHCO3 water types earlier on in the evolutionary pathway, and variability in gas reservoir characteristics. Using recently published data for CSG groundwaters in different basins, the characteristic major ion relationships identified for new data presented in this study were also observed in other CSG groundwaters from Australia, as well as for those in the Illinois Basin in the USA. This observation suggests that where coal maceral content and the dominant methanogenic pathway are similar, and where organic-bound chlorine is relatively abundant, distinct hydrochemical responses may be observed. Comparisons with published data of other NaHCO3 water types in non-CSG environments suggest that these characteristic major ion relationships described here can: i) serve as an indicator of potential CSG groundwaters in certain coal-bearing aquifers that contain methane; and ii) help in the development of strategic sampling programmes for CSG exploration and to monitor potential impacts of CSG activities on groundwater resources.

In areas of potential coal seam gas (CSG) development, understanding interactions between coal-bearing strata and adjacent aquifers and streams is of highest importance, particularly where CSG formations occur at shallow depth. This study tests a combination of hydrochemical and isotopic tracers to investigate the transient nature of hydrochemical processes, inter-aquifer mixing and recharge in a catchment where the coal-bearing aquifer is in direct contact with the alluvial aquifer and surface drainage network. A strong connection was observed between the main stream and underlying alluvium, marked by a similar evolution from fresh Ca-Mg-HCO3 waters in the headwaters towards brackish Ca-Na-Cl composition near the outlet of the catchment, driven by evaporation and transpiration. In the coal-bearing aquifer, by contrast, considerable site-to-site variations were observed, although waters generally had a Na-HCO3-Cl facies and high residual alkalinity values. Increased salinity was controlled by several coexisting processes, including transpiration by plants, mineral weathering and possibly degradation of coal organic matter. Longer residence times and relatively enriched carbon isotopic signatures of the downstream alluvial waters were suggestive of potential interactions with the shallow coal-bearing aquifer. The examination of temporal variations in deuterium excess enabled detection of rapid recharge of the coal-bearing aquifer through highly fractured igneous rocks, particularly at the catchment margins. Most waters collected from the coal-bearing aquifer also showed an enhanced influence of weathering during the wet season, which was likely triggered by the water-rock interaction with fresh recharge waters. An increase in both residual alkalinity and carbon isotopic ratios at two locations indicated inter-aquifer mixing between alluvium and bedrock during the wet season. The results of this study emphasise the need for conducting baseline hydrochemical surveys prior to CSG development in order to describe the transient nature of recharge and inter-aquifer mixing processes. (C) 2015 Elsevier Ltd. All rights reserved.

With the aim of elucidating the seasonal behaviour of rare earth elements (REEs), surface and groundwaters were collected under dry and wet conditions in different hydrological units of the Teviot Brook catchment (Southeast Queensland, Australia). Sampled waters showed a large degree of variability in both REE abundance and normalised patterns. Overall REE abundance ranged over nearly three orders of magnitude, and was consistently lower in the sedimentary bedrock aquifer (18ppt<∑REE<477ppt) than in the other hydrological systems studied. Abundance was greater in springs draining rhyolitic rocks (∑REE=300 and 2054ppt) than in springs draining basalt ranges (∑REE=25 and 83ppt), yet was highly variable in the shallow alluvial groundwater (16ppt<∑REE<5294ppt) and, to a lesser extent, in streamwater (85ppt<∑REE<2198ppt). Generally, waters that interacted with different rock types had different REE patterns. In order to obtain an unbiased characterisation of REE patterns, the ratios between light and middle REEs (R(M/L)) and the ratios between middle and heavy REEs (R(H/M)) were calculated for each sample. The sedimentary bedrock aquifer waters had highly evolved patterns depleted in light REEs and enriched in middle and heavy REEs (0.17<R(M/L)<1.00 and −0.16<R(H/M)<0.93), whereas the springs draining intrusive and extrusive rocks had relatively flat patterns (0.20<R(M/L)<0.38 and −0.16<R(H/M)<0.09). Surface waters were generally enriched in middle REEs (median R(M/L)=0.35 and median R(H/M)=−0.04), and waters from the shallow alluvial aquifer had very diverse patterns with important spatial variations. Samples collected from the alluvium exhibited an increasing influence of the sedimentary bedrock from upgradient to downgradient; typically they showed flat patterns in the upstream section of the alluvium (median R(M/L)=0.21 and median R(H/M)=-0.06) gradually evolving towards patterns depleted in light REEs and enriched in middle and heavy REEs downgradient (median R(M/L)=0.48 and median R(H/M)=0.38). To document the seasonal variations in REE patterns, the difference in ratios between dry and wet sampling campaigns was determined for each repeated sampling location. Contributions from the sedimentary bedrock water to the alluvium during the wet season were identified at two locations (increase from R(H/M)=0.03 and 0.35 to R(H/M)=0.62 and 0.89). The effect of recharge through fractured igneous rocks was also observed in two boreholes intercepting the sedimentary bedrock, where the freshly recharged waters likely contributed to the deeper groundwater flow during the wet season (decrease from R(M/L)=0.81 and 0.56 to R(M/L)=0.46 and 0.17). Results from this study suggest that REEs may be usefully applied as indicators of recharge processes and inter-aquifer mixing. They also underline the importance of conducting seasonal sampling campaigns to capture possible short-term variations in REE patterns and abundance, which is essential to enable a better understanding of hydrological and hydrochemical processes in complex geological settings.

With the predicted acceleration of drought and flood frequency worldwide, it is critical to build knowledge on how shallow groundwater systems respond to rainfall variability. In this study, we characterise the hydrodynamic response of an agricultural catchment located in southeast Queensland, Australia, to the low and high frequency fluctuations in precipitation that occurred in the past 25 years. Strong inter-annual variations in the precipitation input affected surface water flow more substantially than groundwater levels (GWLs). There was also a likely influence of groundwater abstraction on stream baseflow and GWLs, although it was difficult to quantify. Several low frequency oscillations were apparent in the GWL records, particularly in the upstream section of the aquifer, which were neither detectable in rainfall nor in discharge records. Statistically significant episodes of coherence were found at the 2-4-year band between the Nino3.4 index and GWLs for those upstream bores, especially during the 1995-2000 interval, which may be related to a strong La Nina event. In the downstream section of the catchment, higher groundwater persistence probably led to higher filtering of these low frequency oscillations. This paper also proposes a methodology for assessing the dynamic response time of GWLs in shallow aquifer systems to the precipitation input. Response time in a downstream borehole was highly variable temporally, ranging from 11 to 121 days. Rainfall amount was found to significantly affect the short-term dynamics of response time, with elevated rainfall leading to a decreased response time. Annually averaged response time was correlated with the annual number of days with rainfall >30 mm, which was interpreted as the potential for diffuse recharge. Higher recharge potential induced longer response times, probably because of the larger amplitude in GWL variations. (C) 2015 Elsevier B.V. All rights reserved.

Study region The Galilee and Eromanga basins are located in central Queensland, Australia. Both basins are components of the Great Artesian Basin which host some of the most significant groundwater resources in Australia. Study focus This study evaluates the influence of regional faults on groundwater flow in an aquifer/aquitard interbedded succession that form one of the largest Artesian Basins in the world. In order to assess the significance of regional faults as potential barriers or conduits to groundwater flow, vertical displacements of the major aquifers and aquitards were studied at each major fault and the general hydraulic relationship of units that are juxtaposed by the faults were considered. A three-dimensional (3D) geological model of the Galilee and Eromanga basins was developed based on integration of well log data, seismic surfaces, surface geology and elevation data. Geological structures were mapped in detail and major faults were characterised. New hydrological insights for the region Major faults that have been described in previous studies have been confirmed within the 3D geological model domain and a preliminary assessment of their hydraulic significance has been conducted. Previously unknown faults such as the Thomson River Fault (herein named) have also been identified in this study.

In this paper, semi-analytical expressions of the effective hydraulic conductivity (\(K^{E})\) and macrodispersivity (\(\alpha ^{E})\) for 3D steady-state density-dependent groundwater flow are derived using a stationary spectral method. Based on the derived expressions, we present the dependence of \(K^{E}\) and \(\alpha ^{E}\) on the density of fluid under different dispersivity and spatial correlation scale of hydraulic conductivity. The results show that the horizontal \(K^{E}\) and \(\alpha ^{E}\) are not affected by density-induced flow. However, due to gravitational instability of the fluid induced by density contrasts, both vertical \(K^{E}\) and \(\alpha ^{E}\) are found to be reduced slightly when the density factor (\(\gamma \)) is less than 0.01, whereas significant decreases occur when \(\gamma \) exceeds 0.01. Of note, the variation of \(K^{E}\) and \(\alpha ^{E}\) is more significant when local dispersivity is small and the correlation scale of hydraulic conductivity is large.

This study uses borehole geophysical log data of sonic velocity and electrical resistivity to estimate permeability in sandstones in the northern Galilee Basin, Queensland. The prior estimates of permeability are calculated according to the deterministic log–log linear empirical correlations between electrical resistivity and measured permeability. Both negative and positive relationships are influenced by the clay content. The prior estimates of permeability are updated in a Bayesian framework for three boreholes using both the cokriging (CK) method and a normal linear regression (NLR) approach to infer the likelihood function. The results show that the mean permeability estimated from the CK-based Bayesian method is in better agreement with the measured permeability when a fairly apparent linear relationship exists between the logarithm of permeability and sonic velocity. In contrast, the NLR-based Bayesian approach gives better estimates of permeability for boreholes where no linear relationship exists between logarithm permeability and sonic velocity.

We applied a texture-based flow visualisation technique to a numerical hydrodynamic model of the Pumicestone Passage in southeast Queensland, Australia. The quality of the visualisations using our flow visualisation tool, are compared with animations generated using more traditional drogue release plot and velocity contour and vector techniques. The texture-based method is found to be far more effective in visualising advective flow within the model domain. In some instances, it also makes it easier for the researcher to identify specific hydrodynamic features within the complex flow regimes of this shallow tidal barrier estuary as compared with the direct and geometric based methods.

The Lockyer Valley is situated 80 km west of Brisbane and is bounded on the sou th and west by the Great Dividing Range. The valley is a major western sub - catchment of the larger Brisbane River drainage system and is drained by the Lockyer Creek. The Lockyer catchment forms approximately 20% of the total Brisbane River catchment and has an area of around 2900 km2. The Lockyer Creek is an ephemeral drainage system, and the stream and associated alluvium are the main source for irrigation water supply in the Lockyer Valley. The catchment is comprised of a number of well -defined, elongate tributaries in the south, and others in the north, which are more meandering in nature.