Mapping of shallow coastal groundwaters, their hydrology and environmental geochemistry : Pumicestone catchment, Southeast Queensland

Blooms of the toxic cyanobacterium majuscula Lyngbya in the coastal waters of southeast Queensland have caused adverse impacts on both environmental health and human health, and on local economies such as fishing and tourism. A number of studies have confirmed that the main limiting nutrients (“nutrients of concern”) that contribute to these blooms area Fe, DOC, N, P and also pH. This study is conducted to establish the distribution of these parameters in a typical southeast Queensland coastal setting. The study maps the geochemistry of shallow groundwater in the mainland Pumicestone catchment with an emphasis on the nutrients of concern to understand how these nutrients relate to aquifer materials, landuse and anthropogenic activities. The results of the study form a GIS information layer which will be incorporated into a larger GIS model being produced by Queensland Department of Environment and Resource Management (DERM) to support landuse management to avoid/minimize blooms of Lyngbya in Moreton Bay, southeast Queensland, and other similar settings. A total of 38 boreholes were established in the mainland Pumicestone region and four sampling rounds of groundwater carried out in both dry and wet conditions. These groundwater samples were measured in the field for physico-chemical parameters, and in the laboratory analyses for the nutrients of concern, and other major and minor ions. Aquifer materials were confirmed using the Geological Survey of Queensland digital geology map, and geomaterials were assigned to seven categories which are A (sands), B (silts, sandy silts), C (estuarine mud, silts), D (humid soils), E (alluvium), F (sandstone) and G (other bedrock). The results of the water chemistry were examined by use of the software package AquaChem/AqQA, and divided into six groundwater groups, based on groundwater chemical types and location of boreholes. The type of aquifer material and location, and proximity to waterways was found to be important because they affected physico-chemical properties and concentrations of nutrients of concern and dissolved ions. The analytical results showed that iron concentrations of shallow groundwaters were high due to acid sulfate soils, and also mud and silt, but were lower in sand materials. DOC concentrations of these shallow groundwaters in the sand material were high probably due to rapid infiltration. In addition, DOC concentrations in some boreholes were high because they were installed in organic rich wetlands. The pH values of boreholes were from acidic to near neutral; some boreholes with pH values were low (< 4), showing acid sulfate soils in these boreholes. Concentrations of total nitrogen and total phosphorus of groundwaters were generally low, and the main causes of elevated concentrations of total nitrogen and total phosphorus are largely due to animal and human wastes and tend to be found in localized source areas. Comparison of the relative percentage of nitrogen species (NH3/NH4< Org-N, NO3-N and NO2-N) demonstrated that they could be related to sources such as animal waste, residential and agricultural fertilizers, forest and vegetation, mixed residents and farms, and variable setting and vegetation covers. Total concentrations of dissolved ions in sampling round 3 (dry period) were higher than those in sampling round 2 (wet period) due to both evaporation of groundwater in the dry period and the dilution of rainfall in the wet period. This showed that the highest concentrations of nutrients of concern were due to acid sulfate soils, aquifer materials, landuse and anthropogenic activities and were typically in aquifer materials of E (alluvium) and C (estuarine muds) and locations of Burpengary, Caboolture, and Glass Mountain catchments.

Neutrophilic, microaerophilic Fe(II)-oxidizing bacteria are ubiquitous in aquatic habitats of a subtropical Australian coastal catchment (ubiquitous FeOB in catchment aquatic habitats)

We examined the abundance and distribution of neutrophilic, microaerophilic Fe(II)-oxidizing bacteria (FeOB) in aquatic habitats of a highly weathered, subtropical coastal catchment where Fe biogeochemistry is of environmental significance. Laboratory cultivation and microscopy indicated that stalked Gallionella and sheathed Leptothrix-like FeOB were present in microbial mats associated with a circumneutral-pH, groundwater seep and streambank surface sediment,whereas unicellular FeOB werewidespread in surface and subsurface waters, including a seep, shallow stream and estuary-adjacent groundwater. Direct Gallionella-specificPCR detected dominant bacterial members related to Sideroxydans paludicola (95% sequence identity, SI) and Gallionella capsiferriformans (96% SI) in the seep microbialmat. TGGE analysis indicated that themost common FeOB in water enrichment cultures were related to S. lithotrophicus (96% SI). The ubiquity of FeOB in Poona catchment aquatic habitats suggests bacterial Fe(II) oxidation is integral to catchment Fe biogeochemistry.

Review of 3D model development with application to groundwater studies in New Zealand 1996 – 2011

Three dimensional geological modelling techniques have been applied since 1996 with an aim to characterise the lithological and chronological units of New Zealand’s many diverse aquifers. Models of property-scattered data have also been applied to assess physical properties of aquifers and the distribution of groundwater chemistry, including groundwater age, to inform an understanding of groundwater systems. These models, fundamental to understanding groundwater recharge, flow and discharge have found many uses as outlined in this paper.

Modelling sea water intrusion in coastal aquifers using heterogeneous computing

The objective of this PhD research program is to investigate numerical methods for simulating variably-saturated flow and sea water intrusion in coastal aquifers in a high-performance computing environment. The work is divided into three overlapping tasks: to develop an accurate and stable finite volume discretisation and numerical solution strategy for the variably-saturated flow and salt transport equations; to implement the chosen approach in a high performance computing environment that may have multiple GPUs or CPU cores; and to verify and test the implementation. The geological description of aquifers is often complex, with porous materials possessing highly variable properties, that are best described using unstructured meshes. The finite volume method is a popular method for the solution of the conservation laws that describe sea water intrusion, and is well-suited to unstructured meshes. In this work we apply a control volume-finite element (CV-FE) method to an extension of a recently proposed formulation (Kees and Miller, 2002) for variably saturated groundwater flow. The CV-FE method evaluates fluxes at points where material properties and gradients in pressure and concentration are consistently defined, making it both suitable for heterogeneous media and mass conservative. Using the method of lines, the CV-FE discretisation gives a set of differential algebraic equations (DAEs) amenable to solution using higher-order implicit solvers. Heterogeneous computer systems that use a combination of computational hardware such as CPUs and GPUs, are attractive for scientific computing due to the potential advantages offered by GPUs for accelerating data-parallel operations. We present a C++ library that implements data-parallel methods on both CPU and GPUs. The finite volume discretisation is expressed in terms of these data-parallel operations, which gives an efficient implementation of the nonlinear residual function. This makes the implicit solution of the DAE system possible on the GPU, because the inexact Newton-Krylov method used by the implicit time stepping scheme can approximate the action of a matrix on a vector using residual evaluations. We also propose preconditioning strategies that are amenable to GPU implementation, so that all computationally-intensive aspects of the implicit time stepping scheme are implemented on the GPU. Results are presented that demonstrate the efficiency and accuracy of the proposed numeric methods and formulation. The formulation offers excellent conservation of mass, and higher-order temporal integration increases both numeric efficiency and accuracy of the solutions. Flux limiting produces accurate, oscillation-free solutions on coarse meshes, where much finer meshes are required to obtain solutions with equivalent accuracy using upstream weighting. The computational efficiency of the software is investigated using CPUs and GPUs on a high-performance workstation. The GPU version offers considerable speedup over the CPU version, with one GPU giving speedup factor of 3 over the eight-core CPU implementation.

How long does it take for aquifer recharge or aquifer discharge processes to reach steady state?

Groundwater flow models are usually characterized as being either transient flow models or steady state flow models. Given that steady state groundwater flow conditions arise as a long time asymptotic limit of a particular transient response, it is natural for us to seek a finite estimate of the amount of time required for a particular transient flow problem to effectively reach steady state. Here, we introduce the concept of mean action time (MAT) to address a fundamental question: How long does it take for a groundwater recharge process or discharge processes to effectively reach steady state? This concept relies on identifying a cumulative distribution function, $F(t;x)$, which varies from $F(0;x)=0$ to $F(t;x) \to \infty$ as $t\to \infty$, thereby providing us with a measurement of the progress of the system towards steady state. The MAT corresponds to the mean of the associated probability density function $f(t;x) = \dfrac{dF}{dt}$, and we demonstrate that this framework provides useful analytical insight by explicitly showing how the MAT depends on the parameters in the model and the geometry of the problem. Additional theoretical results relating to the variance of $f(t;x)$, known as the variance of action time (VAT), are also presented. To test our theoretical predictions we include measurements from a laboratory–scale experiment describing flow through a homogeneous porous medium. The laboratory data confirms that the theoretical MAT predictions are in good agreement with measurements from the physical model.

Groundwater table response to sea level rise and its impact on pavement structure

It has been predicted that sea level will rise about 0.8 m by 2100. Consequently, seawater can intrude into the coastal aquifers and change the level of groundwater table. A raise in groundwater table due to seawater intrusion threats the coastal infrastructure such as road pavements. The mechanical properties of subgrade materials will change due to elevated rise of groundwater table, leading to pavement weakening and decreasing the subgrade strength and stiffness. This paper presents an assessment of the vulnerability of subgrade in coastal areas to change in groundwater table due to sea-level rise. A simple bathtub approach is applied for estimating the groundwater level changes according to sea-level rise. Then the effect of groundwater level changes on the soil water content (SWC) of a single column of fine-sand soil is simulated using MIKE SHE. The impact of an increase in moisture content on subgrade strength/stiffness is assessed for a number of scenarios.

Road infrastructure vulnerability to groundwater table variation due to sea level rise

This study was conducted to assess the vulnerability of coastal road infrastructures due to climate change induced sea level rise and extreme weather conditions through the estimation of road subgrade strength reduction as a result of changes in soil moisture content. The study area located in the Gold Coast, Australia highlighted that the risk is significant. In wet seasons or areas with wet condition, the groundwater table is already high, so even a small change in the groundwater table can raise the risk of inundation; particularly, in areas with existing shallow groundwater. The predicted risk of a high groundwater table on road infrastructure is a long-term hazard. Therefore, there is time to undertake some management plans to decrease the possible risks, for instance, some deep root plants could be planted along the roads with a high level of risk, to decrease the groundwater table elevation.

An investigation of nucleation events in a coastal urban environment in the southern hemisphere

The occurrence of and conditions favourable to nucleation were investigated at an industrial and commercial coastal location in Brisbane, Australia during five different campaigns covering a total period of 13 months. To identify potential nucleation events, the difference in number concentration in the size range 14-30 nm (N14-30) between consecutive observations was calculated using first-order differencing. The data showed that nucleation events were a rare occurrence, and that in the absence of nucleation the particle number was dominated by particles in the range 30-300 nm. In many instances, total particle concentration declined during nucleation. There was no clear pattern in change in NO and NO2 concentrations during the events. SO2 concentration, in the majority of cases, declined during nucleation but there were exceptions. Most events took place in summer, followed by winter and then spring, and no events were observed for the autumn campaigns. The events were associated with sea breeze and long-range transport. Roadside emissions, in contrast, did not contribute to nucleation, probably due to the predominance of particles in the range 50-100 nm associated with these emissions.