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.

A method of lines approach for modelling saltwater intrusion in coastal aquifiers

The equations governing saltwater intrusion in coastal aquifers are complex. Backward Euler time stepping approaches are often used to advance the solution to these equations in time, which typically requires that small time steps be taken in order to ensure that an accurate solution is obtained. We show that a method of lines approach incorporating variable order backward differentiation formulas can greatly improve the efficiency of the time stepping process.

Visualisation of groundwater flow using texture based visualisation techniques

Texture based techniques for visualisation of unsteady vector fields have been applied for the visualisation of a Finite volume model for variably saturated groundwater flow through porous media. This model has been developed by staff in the School of Mathematical Sciences QUT for the study of salt water intrusion into coastal aquifers. This presentation discusses the implementation and effectiveness of the IBFV algorithm in the context of visualisation of the groundwater simulation outputs.

Does sea-level rise have an impact on saltwater intrusion?

Climate change effects are expected to substantially raise the average sea level. It is widely assumed that this raise will have a severe adverse impact on saltwater intrusion processes in coastal aquifers. In this study we hypothesize that a natural mechanism, identified as the “lifting process” has the potential to mitigate or in some cases completely reverse the adverse intrusion effects induced by sea-level rise. A detailed numerical study using the MODFLOW-family computer code SEAWAT, was completed to test this hypothesis and to understand the effects of this lifting process in both confined and unconfined systems. Our conceptual simulation results show that if the ambient recharge remains constant, the sea-level rise will have no long-term impact (i.e., it will not affect the steady-state salt wedge) on confined aquifers. Our transient confined flow simulations show a self-reversal mechanism where the wedge which will initially intrude into the formation due to the sea-level rise would be naturally driven back to the original position. In unconfined systems, the lifting process would have a lesser influence due to changes in the value of effective transmissivity. A detailed sensitivity analysis was also completed to understand the sensitivity of this self-reversal effect to various aquifer parameters.

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.

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.

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.

New particle formation and growth at a remote, sub-tropical coastal location

A month-long intensive measurement campaign was conducted in March/April 2007 at Agnes Water, a remote coastal site just south of the Great Barrier Reef on the east coast of Australia. Particle and ion size distributions were continuously measured during the campaign. Coastal nucleation events were observed in clean, marine air masses coming from the south-east on 65% of the days. The events usually began at ~10:00 local time and lasted for 1-4 hrs. They were characterised by the appearance of a nucleation mode with a peak diameter of ~10 nm. The freshly nucleated particles grew within 1-4 hrs up to sizes of 20-50 nm. The events occurred when solar intensity was high (~1000 W m-2) and RH was low (~60%). Interestingly, the events were not related to tide height. The volatile and hygroscopic properties of freshly nucleated particles (17-22.5 nm), simultaneously measured with a volatility-hygroscopicity-tandem differential mobility analyser (VH-TDMA), were used to infer chemical composition. The majority of the volume of these particles was attributed to internally mixed sulphate and organic components. After ruling out coagulation as a source of significant particle growth, we conclude that the condensation of sulphate and/or organic vapours was most likely responsible for driving particle growth during the nucleation events. We cannot make any direct conclusions regarding the chemical species that participated in the initial particle nucleation. However, we suggest that nucleation may have resulted from the photo-oxidation products of unknown sulphur or organic vapours emitted from the waters of Hervey Bay, or from the formation of DMS-derived sulphate clusters over the open ocean that were activated to observable particles by condensable vapours emitted from the nutrient rich waters around Fraser Island or Hervey Bay. Furthermore, a unique and particularly strong nucleation event was observed during northerly wind. The event began early one morning (08:00) and lasted almost the entire day resulting in the production of a large number of ~80 nm particles (average modal concentration during the event was 3200 cm-3). The Great Barrier Reef was the most likely source of precursor vapours responsible for this event.

Enriching relationships : community engagment toward sustainable coastal futures

Coastal communities face the social, cultural and environmental challenges of managing rapid urban and industrial development, expanding tourism, and sensitive ecological environments. Enriching relationships between communities and universities through a structured engagement process can deliver integrated options towards sustainable coastal futures. This process draws on the embedded knowledge and values of all participants in the relationship, and offers a wide and affordable range of options for the future. This paper reviews lessons learnt from two projects with coastal communities, and discusses their application in a third. Queensland University of Technology has formed collaborative partnerships with industry in Queensland's Wide Bay-Burnett region to undertake a series of planning and design projects with community engagement as a central process. Senior students worked with community and produced design and planning drawings and reports outlining future options for project areas. A reflective approach has been adopted by the authors to assess the engagement process and outcomes of each project to learn lessons to apply in the next. Methods include surveying community and student participants regarding the value they place on process and outcomes respectively in planning for a sustainable future. All project participants surveyed have placed high importance on the process of engagement, emphasising the value of developing relationships between all project partners. The quality of these relationships is central to planning for sustainable futures, and while the outcomes the students deliver are valued, it is as much for their catalytic role as for their contents. Design and planning projects through community engagement have been found to develop innovative responses to the challenges faced by coastal communities seeking direction toward sustainable futures. The enrichment of engagement relationships and processes has an important influence on the quality of these design and planning responses.