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.

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.

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.

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.

Use of 3D visualisation to enhance groundwater stable isotope data and its interpretation, Lockyer Valley, southeast Queensland

The use of stable isotope ratios δ18O and δ2H are well established in assessment of groundwater systems and their hydrology. The conventional approach is based on x/y plots and relation to various MWL’s, and plots of either ratio against parameters such as Clor EC. An extension of interpretation is the use of 2D maps and contour plots, and 2D hydrogeological vertical sections. An enhancement of presentation and interpretation is the production of “isoscapes”, usually as 2.5D surface projections. We have applied groundwater isotopic data to a 3D visualisation, using the alluvial aquifer system of the Lockyer Valley. The 3D framework is produced in GVS (Groundwater Visualisation System). This format enables enhanced presentation by displaying the spatial relationships and allowing interpolation between “data points” i.e. borehole screened zones where groundwater enters. The relative variations in the δ18O and δ2H values are similar in these ambient temperature systems. However, δ2H better reflects hydrological processes, whereas δ18O also reflects aquifer/groundwater exchange reactions. The 3D model has the advantage that it displays borehole relations to spatial features, enabling isotopic ratios and their values to be associated with, for example, bedrock groundwater mixing, interaction between aquifers, relation to stream recharge, and to near-surface and return irrigation water evaporation. Some specific features are also shown, such as zones of leakage of deeper groundwater (in this case with a GAB signature). Variations in source of recharging water at a catchment scale can be displayed. Interpolation between bores is not always possible depending on numbers and spacing, and by elongate configuration of the alluvium. In these cases, the visualisation uses discs around the screens that can be manually expanded to test extent or intersections. Separate displays are used for each of δ18O and δ2H and colour coding for isotope values.

Hydrogeological framework, conceptual and numerical groundwater flow model of Laidley Creek catchment, Queensland, Australia

This thesis studies the water resources of Laidley Creek catchment within the Lockyer Valley where groundwater is used for intensive irrigation of crops. A holistic approach was used to consider groundwater within the total water cycle. The project mapped the geology, measured stream flows and groundwater levels, and analysed the chemistry of the waters. These data were integrated within a catchment-wide conceptual model, including historic and rainfall records. From this a numerical simulation was produced to test data validity and develop predictions of behaviour, which can support management decisions, particularly in times of variable climate.

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.

A general geomorphological recession flow model for river basins

Recession flows in a basin are controlled by the temporal evolution of its active drainage network (ADN). The geomorphological recession flow model (GRFM) assumes that both the rate of flow generation per unit ADN length (q) and the speed at which ADN heads move downstream (c) remain constant during a recession event. Thereby, it connects the power law exponent of -dQ/dt versus Q (discharge at the outlet at time t) curve, , with the structure of the drainage network, a fixed entity. In this study, we first reformulate the GRFM for Horton-Strahler networks and show that the geomorphic ((g)) is equal to D/(D-1), where D is the fractal dimension of the drainage network. We then propose a more general recession flow model by expressing both q and c as functions of Horton-Strahler stream order. We show that it is possible to have = (g) for a recession event even when q and c do not remain constant. The modified GRFM suggests that is controlled by the spatial distribution of subsurface storage within the basin. By analyzing streamflow data from 39 U.S. Geological Survey basins, we show that is having a power law relationship with recession curve peak, which indicates that the spatial distribution of subsurface storage varies across recession events. Key Points The GRFM is reformulated for Horton-Strahler networks. The GRFM is modified by allowing its parameters to vary along streams. Sub-surface storage distribution controls recession flow characteristics.

Water–rock interactions : an investigation of the relationships between mineralogy and groundwater composition and flow in a subtropical basalt aquifer

A holistic study of the composition of the basalt groundwaters of the Atherton Tablelands region in Queensland, Australia was undertaken to elucidate possible mechanisms for the evolution of these very low salinity, silica- and bicarbonate-rich groundwaters. It is proposed that aluminosilicate mineral weathering is the major contributing process to the overall composition of the basalt groundwaters. The groundwaters approach equilibrium with respect to the primary minerals with increasing pH and are mostly in equilibrium with the major secondary minerals (kaolinite and smectite), and other secondary phases such as goethite, hematite, and gibbsite, which are common accessory minerals in the Atherton basalts. The mineralogy of the basalt rocks, which has been examined using X-ray diffraction and whole rock geochemistry methods, supports the proposed model for the hydrogeochemical evolution of these groundwaters: precipitation + CO 2 (atmospheric + soil) + pyroxene + feldspars + olivine yields H 4SiO 4, HCO 3 -, Mg 2+, Na +, Ca 2+ + kaolinite and smectite clays + amorphous or crystalline silica + accessory minerals (hematite, goethite, gibbsite, carbonates, zeolites, and pyrite). The variations in the mineralogical content of these basalts also provide insights into the controls on groundwater storage and movement in this aquifer system. The fresh and weathered vesicular basalts are considered to be important in terms of zones of groundwater occurrence, while the fractures in the massive basalt are important pathways for groundwater movement.