Groundwater flow model of the Logan river alluvial aquifer system Josephville, South East Queensland

The study focuses on an alluvial plain situated within a large meander of the Logan River at Josephville near Beaudesert which supports a factory that processes gelatine. The plant draws water from on site bores, as well as the Logan River, for its production processes and produces approximately 1.5 ML per day (Douglas Partners, 2004) of waste water containing high levels of dissolved ions. At present a series of treatment ponds are used to aerate the waste water reducing the level of organic matter; the water is then used to irrigate grazing land around the site. Within the study the hydrogeology is investigated, a conceptual groundwater model is produced and a numerical groundwater flow model is developed from this. On the site are several bores that access groundwater, plus a network of monitoring bores. Assessment of drilling logs shows the area is formed from a mixture of poorly sorted Quaternary alluvial sediments with a laterally continuous aquifer comprised of coarse sands and fine gravels that is in contact with the river. This aquifer occurs at a depth of between 11 and 15 metres and is overlain by a heterogeneous mixture of silts, sands and clays. The study investigates the degree of interaction between the river and the groundwater within the fluvially derived sediments for reasons of both environmental monitoring and sustainability of the potential local groundwater resource. A conceptual hydrogeological model of the site proposes two hydrostratigraphic units, a basal aquifer of coarse-grained materials overlain by a thick semi-confining unit of finer materials. From this, a two-layer groundwater flow model and hydraulic conductivity distribution was developed based on bore monitoring and rainfall data using MODFLOW (McDonald and Harbaugh, 1988) and PEST (Doherty, 2004) based on GMS 6.5 software (EMSI, 2008). A second model was also considered with the alluvium represented as a single hydrogeological unit. Both models were calibrated to steady state conditions and sensitivity analyses of the parameters has demonstrated that both models are very stable for changes in the range of ± 10% for all parameters and still reasonably stable for changes up to ± 20% with RMS errors in the model always less that 10%. The preferred two-layer model was found to give the more realistic representation of the site, where water level variations and the numerical modeling showed that the basal layer of coarse sands and fine gravels is hydraulically connected to the river and the upper layer comprising a poorly sorted mixture of silt-rich clays and sands of very low permeability limits infiltration from the surface to the lower layer. The paucity of historical data has limited the numerical modelling to a steady state one based on groundwater levels during a drought period and forecasts for varying hydrological conditions (e.g. short term as well as prolonged dry and wet conditions) cannot reasonably be made from such a model. If future modelling is to be undertaken it is necessary to establish a regular program of groundwater monitoring and maintain a long term database of water levels to enable a transient model to be developed at a later stage. This will require a valid monitoring network to be designed with additional bores required for adequate coverage of the hydrogeological conditions at the Josephville site. Further investigations would also be enhanced by undertaking pump testing to investigate hydrogeological properties in the aquifer.

Visualisation modelling in 3D of an alluvial aquifer system at a valley-wide scale to understand groundwater status : Lockyer, QLD

Visualisation provides a method to efficiently convey and understand the complex nature and processes of groundwater systems. This technique has been applied to the Lockyer Valley to aid in comprehending the current condition of the system. The Lockyer Valley in southeast Queensland hosts intensive irrigated agriculture sourcing groundwater from alluvial aquifers. The valley is around 3000 km2 in area and the alluvial deposits are typically 1-3 km wide and to 20-35 m deep in the main channels, reducing in size in subcatchments. The configuration of the alluvium is of a series of elongate “fingers”. In this roughly circular valley recharge to the alluvial aquifers is largely from seasonal storm events, on the surrounding ranges. The ranges are overlain by basaltic aquifers of Tertiary age, which overall are quite transmissive. Both runoff from these ranges and infiltration into the basalts provided ephemeral flow to the streams of the valley. Throughout the valley there are over 5,000 bores extracting alluvial groundwater, plus lesser numbers extracting from underlying sandstone bedrock. Although there are approximately 2500 monitoring bores, the only regularly monitored area is the formally declared management zone in the lower one third. This zone has a calibrated Modflow model (Durick and Bleakly, 2000); a broader valley Modflow model was developed in 2002 (KBR), but did not have extensive extraction data for detailed calibration. Another Modflow model focused on a central area river confluence (Wilson, 2005) with some local production data and pumping test results. A recent subcatchment simulation model incorporates a network of bores with short-period automated hydrographic measurements (Dvoracek and Cox, 2008). The above simulation models were all based on conceptual hydrogeological models of differing scale and detail.

The Cressbrook Creek alluvial aquifer system, Southeast Queensland : hydrochemistry and isotopes to determine hydrological processes and response to floods

This study developed an understanding of hydrological processes within the Cressbrook Creek catchment of the upper Brisbane River, in particular for the alluvial aquifers. Those aquifers within the lower catchment are used for intensive irrigation, and have been impacted by long-term drought followed by flooding. The study utilised water chemistry, isotopic characters and hydraulic measurements to determine factors such as recharge, links between creeks and groundwater, and variations in water quality. The catchment-wide study will enable improved management of the local water resources.

Simulation of Ground-Water Flow in the Cedar River Alluvial Aquifer Flow System, Cedar Rapids, Iowa Scientific

The Cedar River alluvial aquifer is the primary source of municipal water in the Cedar Rapids, Iowa, area. Since 1992, the U.S. Geological Survey, in cooperation with the City of Cedar Rapids, has investigated the hydrogeology and water quality of the Cedar River alluvial aquifer. This report describes a detailed analysis of the ground-water flow system in the alluvial aquifer, particularly near well field areas. The ground-water flow system in the Cedar Rapids area consists of two main components, the unconsolidated Quaternary deposits and the underlying carbonate bedrock that has a variable fracture density. Quaternary deposits consist of eolian sand, loess, alluvium, and glacial till. Devonian and Silurian bedrock aquifers overlie the Maquoketa Shale (Formation) of Ordovician age, a regional confining unit. Ground-water and surface-water data were collected during the study to better define the hydrogeology of the Cedar River alluvial aquifer and Devonian and Silurian aquifers. Stream stage and discharge, ground-water levels, and estimates of aquifer hydraulic properties were used to develop a conceptual ground-water flow model and to construct and calibrate a model of the flow system. This model was used to quantify the movement of water between the various components of the alluvial aquifer flow system and provide an improved understanding of the hydrology of the alluvial aquifer.

The contribution of the Finite Volume Point Dilution Method to the determination of solute mass flux in an alluvial aquifer.

Groundwater represents one of the most important resources of the world and it is essential to prevent its pollution and to consider remediation intervention in case of contamination. According to the scientific community the characterization and the management of the contaminated sites have to be performed in terms of contaminant fluxes and considering their spatial and temporal evolution. One of the most suitable approach to determine the spatial distribution of pollutant and to quantify contaminant fluxes in groundwater is using control panels. The determination of contaminant mass flux, requires measurement of contaminant concentration in the moving phase (water) and velocity/flux of the groundwater. In this Master Thesis a new solute flux mass measurement approach, based on an integrated control panel type methodology combined with the Finite Volume Point Dilution Method (FVPDM), for the monitoring of transient groundwater fluxes, is proposed. Moreover a new adsorption passive sampler, which allow to capture the variation of solute concentration with time, is designed. The present work contributes to the development of this approach on three key points. First, the ability of the FVPDM to monitor transient groundwater fluxes was verified during a step drawdown test at the experimental site of Hermalle Sous Argentau (Belgium). The results showed that this method can be used, with optimal results, to follow transient groundwater fluxes. Moreover, it resulted that performing FVPDM, in several piezometers, during a pumping test allows to determine the different flow rates and flow regimes that can occurs in the various parts of an aquifer. The second field test aiming to determine the representativity of a control panel for measuring mass flus in groundwater underlined that wrong evaluations of Darcy fluxes and discharge surfaces can determine an incorrect estimation of mass fluxes and that this technique has to be used with precaution. Thus, a detailed geological and hydrogeological characterization must be conducted, before applying this technique. Finally, the third outcome of this work concerned laboratory experiments. The test conducted on several type of adsorption material (Oasis HLB cartridge, TDS-ORGANOSORB 10 and TDS-ORGANOSORB 10-AA), in order to determine the optimum medium to dimension the passive sampler, highlighted the necessity to find a material with a reversible adsorption tendency to completely satisfy the request of the new passive sampling technique.

Recharge Areas and Geochemical Evolution of Groundwater in an Alluvial Aquifer System in the Sultanate of Oman

A regional hydrogeochemical model was developed to evaluate the geochemical evolution of different groundwaters in an alluvial aquifer system in the Interior of Oman. In combination with environmental isotopes the model is able to extract qualitative and quantitative information about recharge, groundwater flow paths and hydraulic connections between different aquifers. The main source of water to the alluvial aquifer along the flow paths ofWadi Abyadh andWadi M’uaydin in the piedmont is groundwater from the high-altitude areas of the Jabal Akhdar and local infiltration along the wadi channels. In contrast, the piedmont alluvial aquifer alongWadi Halfayn is primarily replenished by lateral recharge from the ophiolite foothills to the east besides smaller contributions from the Jabal Akhdar and local infiltration. Further down gradient in the Southern Alluvial Plain aquifer a significant source of recharge is direct infiltration of rain and surface runoff, originating from a moisture source that approaches Oman from the south. The model shows that the main geochemical evolution of the alluvial groundwaters occurs along the flow path from the piedmont to the Southern Alluvial Plain, where dedolomitization is responsible for the observed changes in the chemical and carbon isotope composition in these waters.

Biogeochemistry and microbiology of groundwater during acetate and bicarbonate amendment to an alluvial aquifer

These data are from a field experiment conducted in a shallow alluvial aquifer along the Colorado River in Rifle, Colorado, USA. In this experiment, bicarbonate-promoted uranium desorption and acetate amendment were combined and compared to an acetate amendment-only experiment in the same experimental plot. Data include names and location data for boreholes, geochemical data for all the boreholes between June 1, 2010 and January 1, 2011, microarray data provided as signal to noise ratio (SNR) for individual microarray probes, microarray data provided as signal to noise ratio (SNR) by Genus.

Bacterially mediated iron cycling and associated biogeochemical processes in a subtropical shallow coastal aquifer: implications for groundwater quality

Bacterially mediated iron redox cycling exerts a strong influence on groundwater geochemistry, but few studies have investigated iron biogeochemical processes in coastal alluvial aquifers from a microbiological viewpoint. The shallow alluvial aquifer located adjacent to Poona estuary on the subtropical Southeast Queensland coast represents a redox-stratified system where iron biogeochemical cycling potentially affects water quality. Using a 300 m transect of monitoring wells perpendicular to the estuary, we examined groundwater physico-chemical conditions and the occurrence of cultivable bacterial populations involved in iron (and manganese, sulfur) redox reactions in this aquifer. Results showed slightly acidic and near-neutral pH, suboxic conditions and an abundance of dissolved iron consisting primarily of iron(II) in the majority of wells. The highest level of dissolved iron(III) was found in a well proximal to the estuary most likely a result of iron curtain effects due to tidal intrusion. A number of cultivable, (an)aerobic bacterial populations capable of diverse carbon, iron, or sulfur metabolism coexisted in groundwater redox transition zones. Our findings indicated aerobic, heterotrophic respiration and bacterially mediated iron/sulfur redox reactions were integral to carbon cycling in the aquifer. High abundances of dissolved iron and cultivable iron and sulfur bacterial populations in estuary-adjacent aquifers have implications for iron transport to marine waters. This study demonstrated bacterially mediated iron redox cycling and associated biogeochemical processes in subtropical coastal groundwaters using culture-based methods.

Assessment of flood-related recharge to alluvial aquifers of an irrigated catchment following extended drought conditions using 3D visualisation and environmental tracers, Lockyer Valley, southeast Queensland, Australia.

The Lockyer Valley in southeast Queensland, Australia, hosts an economically significant alluvial aquifer system which has been impacted by prolonged drought conditions (~1997 to ~ 2009). Throughout this time, the system was under continued groundwater extraction, resulting in severe aquifer depletion. By 2008, much of the aquifer was at <30% of storage but some relief occurred with rains in early 2009. However, between December 2010 and January 2011, most of southeast Queensland experienced unprecedented flooding, which generated significant aquifer recharge. In order to understand the spatial and temporal controls of groundwater recharge in the alluvium, a detailed 3D lithological property model of gravels, sands and clays was developed using GOCAD software. The spatial distribution of recharge throughout the catchment was assessed using hydrograph data from about 400 groundwater observation wells screened at the base of the alluvium. Water levels from these bores were integrated into a catchment-wide 3D geological model using the 3D geological modelling software GOCAD; the model highlights the complexity of recharge mechanisms. To support this analysis, groundwater tracers (e.g. major and minor ions, stable isotopes, 3H and 14C) were used as independent verification. The use of these complementary methods has allowed the identification of zones where alluvial recharge primarily occurs from stream water during episodic flood events. However, the study also demonstrates that in some sections of the alluvium, rainfall recharge and discharge from the underlying basement into the alluvium are the primary recharge mechanisms of the alluvium. This is indicated by the absence of any response to the flood, as well as the observed old radiocarbon ages and distinct basement water chemistry signatures at these locations. Within the 3D geological model, integration of water chemistry and time-series displays of water level surfaces before and after the flood suggests that the spatial variations of the flood response in the alluvium are primarily controlled by the valley morphology and lithological variations within the alluvium. The integration of time-series of groundwater level surfaces in the 3D geological model also enables the quantification of the volumetric change of groundwater stored in the unconfined sections of this alluvial aquifer during drought and following flood events. The 3D representation and analysis of hydraulic and recharge information has considerable advantages over the traditional 2D approach. For example, while many studies focus on singular aspects of catchment dynamics and groundwater-surface water interactions, the 3D approach is capable of integrating multiple types of information (topography, geological, hydraulic, water chemistry and spatial) into a single representation which provides valuable insights into the major factors controlling aquifer processes.

Assessment of interactions between alluvial, Great Artesian Basin (GAB) and volcanic aquifers using 3D visualisation and environmental tracers, Lockyer Valley, southeast Queensland, Australia

A detailed 3D lithological model framework was developed using GOCAD software to understand interactions between alluvial, volcanic and GAB aquifers and the spatial and temporal distribution of groundwater recharge to the alluvium of the Lockyer Valley. Groundwater chemistry, isotope data (H20-δ2H and δ18O , 87Sr/86Sr, 3H and 14C) and groundwater level time-series data from approximately 550 observation wells were integrated into the catchment-wide 3D model to assess the recharge processes involved. This approach enabled the identification of zones where recharge to the alluvium primarily occurs from stream water during episodic flood events. Importantly, the study also demonstrates that in some sections of the alluvium recharge is also from storm rainfall and seepage discharge from the underlying GAB aquifers. These other sources of recharge are indicated by (a) the absence of a response of groundwater levels to flooding in some areas, (b) old radiocarbon ages, and (c) distinct bedrock water chemistry and δ2H and δ18O signatures in alluvial groundwater at these locations. Integration of isotopes, water chemistry and time-series displays of groundwater levels before and after the 2010/2011 flood into the 3D model suggest that the spatial variations in the alluvial groundwater response are mostly controlled by valley morphology and lithological (i.e. permeability) variations within the alluvium. Examination of the groundwater level variations in the 3D model also enabled quantification of the volumetric change of groundwater stored in the unconfined alluvial aquifer prior to and post-flood events.

Multivariate statistical analysis of hydrochemical data to assess alluvial aquifer–stream connectivity during drought and flood : Cressbrook Creek, southeast Queensland, Australia

A catchment-scale multivariate statistical analysis of hydrochemistry enabled assessment of interactions between alluvial groundwater and Cressbrook Creek, an intermittent drainage system in southeast Queensland, Australia. Hierarchical cluster analyses and principal component analysis were applied to time-series data to evaluate the hydrochemical evolution of groundwater during periods of extreme drought and severe flooding. A simple three-dimensional geological model was developed to conceptualise the catchment morphology and the stratigraphic framework of the alluvium. The alluvium forms a two-layer system with a basal coarse-grained layer overlain by a clay-rich low-permeability unit. In the upper and middle catchment, alluvial groundwater is chemically similar to streamwater, particularly near the creek (reflected by high HCO3/Cl and K/Na ratios and low salinities), indicating a high degree of connectivity. In the lower catchment, groundwater is more saline with lower HCO3/Cl and K/Na ratios, notably during dry periods. Groundwater salinity substantially decreased following severe flooding in 2011, notably in the lower catchment, confirming that flooding is an important mechanism for both recharge and maintaining groundwater quality. The integrated approach used in this study enabled effective interpretation of hydrological processes and can be applied to a variety of hydrological settings to synthesise and evaluate large hydrochemical datasets.

Seasonal and spatial variations in rare earth elements to identify inter-aquifer linkages and recharge processes in an Australian catchment

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