Hydrogeological Features of the Lockyer Valley

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.

Brave new green world : consequences of a carbon economy for the conservation of Australian biodiversity

Pricing greenhouse gas emissions is a burgeoning and possibly lucrative financial means for climate change mitigation. Emissions pricing is being used to fund emissions-abatement technologies and to modify land management to improve carbon sequestration and retention. Here we discuss the principal land-management options under existing and realistic future emissions-price legislation in Australia, and examine them with respect to their anticipated direct and indirect effects on biodiversity. The main ways in which emissions price-driven changes to land management can affect biodiversity are through policies and practices for (1) environmental plantings for carbon sequestration, (2) native regrowth, (3) fire management, (4) forestry, (5) agricultural practices (including cropping and grazing), and (6) feral animal control. While most land-management options available to reduce net greenhouse gas emissions offer clear advantages to increase the viability of native biodiversity, we describe several caveats regarding potentially negative outcomes, and outline components that need to be considered if biodiversity is also to benefit from the new carbon economy. Carbon plantings will only have real biodiversity value if they comprise appropriate native tree species and provide suitable habitats and resources for valued fauna. Such plantings also risk severely altering local hydrology and reducing water availability. Management of regrowth post-agricultural abandonment requires setting appropriate baselines and allowing for thinning in certain circumstances, and improvements to forestry rotation lengths would likely increase carbon-retention capacity and biodiversity value. Prescribed burning to reduce the frequency of high-intensity wildfires in northern Australia is being used as a tool to increase carbon retention. Fire management in southern Australia is not readily amenable for maximising carbon storage potential, but will become increasingly important for biodiversity conservation as the climate warms. Carbon price-based modifications to agriculture that would benefit biodiversity include reductions in tillage frequency and livestock densities, reductions in fertiliser use, and retention and regeneration of native shrubs; however, anticipated shifts to exotic perennial grass species such as buffel grass and kikuyu could have net negative implications for native biodiversity. Finally, it is unlikely that major reductions in greenhouse gas emissions arising from feral animal control are possible, even though reduced densities of feral herbivores will benefit Australian biodiversity greatly.

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.

Role of land use and seasonal factors in water quality degradation

Surface water and groundwater are the most important water sources in the natural environment. Land use and seasonal factors play an important role in influencing the quality of these water sources. An in-depth understanding of the role of these two influential factors can help to implement an effective catchment management strategy for the protection of these water sources. This paper discusses the outcomes of an extensive research study which investigated the role of land use and seasonal factors on surface water and groundwater pollution in a mixed land use coastal catchment. The study confirmed that the influence exerted on the water environment by seasonal factors is secondary to that of land use. Furthermore, the influence of land use and seasonal factors on surface water and groundwater quality varies with the pollutant species. This highlights the need to specifically take into consideration the targeted pollutants and the key influential factors for the effective protection of vulnerable receiving water environments.

Assessing resilience of water resources systems under the impacts of climate change

This paper presents an approach to assess the resilience of a water supply system under the impacts of climate change. Changes to climate characteristics such as rainfall, evapotranspiration and temperature can result in changes to the global hydrological cycle and thereby adversely impact on the ability of water supply systems to meet service standards in the future. Changes to the frequency and characteristics of floods and droughts as well as the quality of water provided by groundwater and surface water resources are the other consequences of climate change that will affect water supply system functionality. The extent and significance of these changes underline the necessity for assessing the future functionality of water supply systems under the impacts of climate change. Resilience can be a tool for assessing the ability of a water supply system to meet service standards under the future climate conditions. The study approach is based on defining resilience as the ability of a system to absorb pressure without going into failure state as well as its ability to achieve an acceptable level of function quickly after failure. In order to present this definition in the form of a mathematical function, a surrogate measure of resilience has been proposed in this paper. In addition, a step-by-step approach to estimate resilience of water storage reservoirs is presented. This approach will enable a comprehensive understanding of the functioning of a water storage reservoir under future climate scenarios and can also be a robust tool to predict future challenges faced by water supply systems under the consequence of climate change.

Determination of vertical hydraulic conductivity of aquitards in a multilayered leaky system using water-level signals in adjacent aquifers

This paper presents a methodology for determining the vertical hydraulic conductivity (Kv) of an aquitard, in a multilayered leaky system, based on the harmonic analysis of arbitrary water-level fluctuations in aquifers. As a result, Kv of the aquitard is expressed as a function of the phase-shift of water-level signals measured in the two adjacent aquifers. Based on this expression, we propose a robust method to calculate Kv by employing linear regression analysis of logarithm transformed frequencies and phases. The frequencies, where the Kv are calculated, are identified by coherence analysis. The proposed methods are validated by a synthetic case study and are then applied to the Westbourne and Birkhead aquitards, which form part of a five-layered leaky system in the Eromanga Basin, Australia.

A micro-level indexing model for the assessment of sustainable urban ecosystems

During the last several decades, the quality of natural resources and their services have been exposed to significant degradation from increased urban populations combined with the sprawl of settlements, development of transportation networks and industrial activities (Dorsey, 2003; Pauleit et al., 2005). As a result of this environmental degradation, a sustainable framework for urban development is required to provide the resilience of natural resources and ecosystems. Sustainable urban development refers to the management of cities with adequate infrastructure to support the needs of its population for the present and future generations as well as maintain the sustainability of its ecosystems (UNEP/IETC, 2002; Yigitcanlar, 2010). One of the important strategic approaches for planning sustainable cities is „ecological planning‟. Ecological planning is a multi-dimensional concept that aims to preserve biodiversity richness and ecosystem productivity through the sustainable management of natural resources (Barnes et al., 2005). As stated by Baldwin (1985, p.4), ecological planning is the initiation and operation of activities to direct and control the acquisition, transformation, disruption and disposal of resources in a manner capable of sustaining human activities with a minimum disruption of ecosystem processes. Therefore, ecological planning is a powerful method for creating sustainable urban ecosystems. In order to explore the city as an ecosystem and investigate the interaction between the urban ecosystem and human activities, a holistic urban ecosystem sustainability assessment approach is required. Urban ecosystem sustainability assessment serves as a tool that helps policy and decision-makers in improving their actions towards sustainable urban development. There are several methods used in urban ecosystem sustainability assessment among which sustainability indicators and composite indices are the most commonly used tools for assessing the progress towards sustainable land use and urban management. Currently, a variety of composite indices are available to measure the sustainability at the local, national and international levels. However, the main conclusion drawn from the literature review is that they are too broad to be applied to assess local and micro level sustainability and no benchmark value for most of the indicators exists due to limited data availability and non-comparable data across countries. Mayer (2008, p. 280) advocates that by stating "as different as the indices may seem, many of them incorporate the same underlying data because of the small number of available sustainability datasets". Mori and Christodoulou (2011) also argue that this relative evaluation and comparison brings along biased assessments, as data only exists for some entities, which also means excluding many nations from evaluation and comparison. Thus, there is a need for developing an accurate and comprehensive micro-level urban ecosystem sustainability assessment method. In order to develop such a model, it is practical to adopt an approach that uses a method to utilise indicators for collecting data, designate certain threshold values or ranges, perform a comparative sustainability assessment via indices at the micro-level, and aggregate these assessment findings to the local level. Hereby, through this approach and model, it is possible to produce sufficient and reliable data to enable comparison at the local level, and provide useful results to inform the local planning, conservation and development decision-making process to secure sustainable ecosystems and urban futures. To advance research in this area, this study investigated the environmental impacts of an existing urban context by using a composite index with an aim to identify the interaction between urban ecosystems and human activities in the context of environmental sustainability. In this respect, this study developed a new comprehensive urban ecosystem sustainability assessment tool entitled the „Micro-level Urban-ecosystem Sustainability IndeX‟ (MUSIX). The MUSIX model is an indicator-based indexing model that investigates the factors affecting urban sustainability in a local context. The model outputs provide local and micro-level sustainability reporting guidance to help policy-making concerning environmental issues. A multi-method research approach, which is based on both quantitative analysis and qualitative analysis, was employed in the construction of the MUSIX model. First, a qualitative research was conducted through an interpretive and critical literature review in developing a theoretical framework and indicator selection. Afterwards, a quantitative research was conducted through statistical and spatial analyses in data collection, processing and model application. The MUSIX model was tested in four pilot study sites selected from the Gold Coast City, Queensland, Australia. The model results detected the sustainability performance of current urban settings referring to six main issues of urban development: (1) hydrology, (2) ecology, (3) pollution, (4) location, (5) design, and; (6) efficiency. For each category, a set of core indicators was assigned which are intended to: (1) benchmark the current situation, strengths and weaknesses, (2) evaluate the efficiency of implemented plans, and; (3) measure the progress towards sustainable development. While the indicator set of the model provided specific information about the environmental impacts in the area at the parcel scale, the composite index score provided general information about the sustainability of the area at the neighbourhood scale. Finally, in light of the model findings, integrated ecological planning strategies were developed to guide the preparation and assessment of development and local area plans in conjunction with the Gold Coast Planning Scheme, which establishes regulatory provisions to achieve ecological sustainability through the formulation of place codes, development codes, constraint codes and other assessment criteria that provide guidance for best practice development solutions. These relevant strategies can be summarised as follows: • Establishing hydrological conservation through sustainable stormwater management in order to preserve the Earth’s water cycle and aquatic ecosystems; • Providing ecological conservation through sustainable ecosystem management in order to protect biological diversity and maintain the integrity of natural ecosystems; • Improving environmental quality through developing pollution prevention regulations and policies in order to promote high quality water resources, clean air and enhanced ecosystem health; • Creating sustainable mobility and accessibility through designing better local services and walkable neighbourhoods in order to promote safe environments and healthy communities; • Sustainable design of urban environment through climate responsive design in order to increase the efficient use of solar energy to provide thermal comfort, and; • Use of renewable resources through creating efficient communities in order to provide long-term management of natural resources for the sustainability of future generations.

A parcel-scale assessment tool to measure sustainability through urban ecosystem components : the MUSIX model

In recent years a number of urban sustainability assessment frameworks are developed to better inform policy formulation and decision-making processes. This paper introduces one of these attempts in developing a comprehensive assessment tool—i.e., Micro-level Urban-ecosystem Sustainability IndeX (MUSIX). Being an indicator-based indexing model, MUSIX investigates the environmental impacts of land-uses on urban sustainability by measuring urban ecosystem components in local scale. The paper presents the methodology of MUSIX and demonstrates the performance of the model in a pilot test-bed—i.e., in Gold Coast, Australia. The model provides useful insights on the sustainability performance of the test-bed area. The parcel-scale findings of the indicators are used to identify local problems considering six main issues of urban development—i.e., hydrology; ecology; pollution; location; design, and; efficiency. The composite index score is used to propose betterment strategies to guide the development of local area plans in conjunction with the City's Planning Scheme. In overall, this study has shown that parcel-scale environmental data provides an overview of the local sustainability in urban areas as in the example of Gold Coast, which can also be used for setting environmental policy, objectives and targets.

Exact series solutions of reactive transport models with general initial conditions

Exact solutions of partial differential equation models describing the transport and decay of single and coupled multispecies problems can provide insight into the fate and transport of solutes in saturated aquifers. Most previous analytical solutions are based on integral transform techniques, meaning that the initial condition is restricted in the sense that the choice of initial condition has an important impact on whether or not the inverse transform can be calculated exactly. In this work we describe and implement a technique that produces exact solutions for single and multispecies reactive transport problems with more general, smooth initial conditions. We achieve this by using a different method to invert a Laplace transform which produces a power series solution. To demonstrate the utility of this technique, we apply it to two example problems with initial conditions that cannot be solved exactly using traditional transform techniques.

How much geological detail is needed for gasifier success?

The incipient Underground Coal Gasification (UCG) industry in Queensland, Australia, undertook three trial projects in two Mesozoic basins of southeast Queensland. The experiences of these three operations provide useful retrospective insight into gasifier productivity. This paper identifies key output measures of gasifier ‘success’ including output gas composition, presence of contaminants in groundwater and consistency of chamber operation. Likewise, a review of the geological and hydrogeological understanding of each site prior to gasifier commissioning was undertaken. Productivity parameters from gasification were then correlated against the level of baseline geological/hydrogeological understanding for each site. The aim of the study was to identify the optimum scope of geological and hydrogeological understanding required at the site assessment phase to ensure safe, maximum gasifier output during production phase. This approach allows identification of poor or unexpected performance that is attributable to pre-existing uncertainty. A historical review of gasifier conditions inferred from the three trial projects is presented. Hence from the Queensland experiences it is possible to identify what aspects of baseline geological understanding should be clearly understood at the site selection phase in order to limit anomalous gasifier performance and undesirable deviations, and maximise production output.

Electrical imaging and fluid modeling of convective fingering in a shallow water-table aquifer

Unstable density-driven flow can lead to enhanced solute transport in groundwater. Only recently has the complex fingering pattern associated with free convection been documented in field settings. Electrical resistivity (ER) tomography has been used to capture a snapshot of convective instabilities at a single point in time, but a thorough transient analysis is still lacking in the literature. We present the results of a 2 year experimental study at a shallow aquifer in the United Arab Emirates that was designed to specifically explore the transient nature of free convection. ER tomography data documented the presence of convective fingers following a significant rainfall event. We demonstrate that the complex fingering pattern had completely disappeared a year after the rainfall event. The observation is supported by an analysis of the aquifer halite budget and hydrodynamic modeling of the transient character of the fingering instabilities. Modeling results show that the transient dynamics of the gravitational instabilities (their initial development, infiltration into the underlying lower-density groundwater, and subsequent decay) are in agreement with the timing observed in the time-lapse ER measurements. All experimental observations and modeling results are consistent with the hypothesis that a dense brine that infiltrated into the aquifer from a surficial source was the cause of free convection at this site, and that the finite nature of the dense brine source and dispersive mixing led to the decay of instabilities with time. This study highlights the importance of the transience of free convection phenomena and suggests that these processes are more rapid than was previously understood.

Terrestrial invertebrates of dry river beds are not simply subsets of riparian assemblages

Dry river beds are common worldwide and are rapidly increasing in extent due to the effects of water management and prolonged drought periods due to climate change. While attention has been given to the responses of aquatic invertebrates to drying rivers, few studies exist on the terrestrial invertebrates colonizing dry river beds. Dry river beds are physically harsh and they often differ substantially in substrate, topography, microclimate and inundation frequency from adjacent riparian zones. Given these differences, we predicted that dry river beds provide a unique habitat for terrestrial invertebrates, and that their assemblage composition differs from that in adjacent riparian zones. Dry river beds and riparian zones in Australia and Italy were sampled for terrestrial invertebrates with pitfall traps. Sites differed in substrate type, climate and flow regime. Dry river beds contained diverse invertebrate assemblages and their composition was consistently different from adjacent riparian zones, irrespective of substrate, climate or hydrology. Although some taxa were shared between dry river beds and riparian zones, 66 of 320 taxa occurred only in dry river beds. Differences were due to species turnover, rather than shifts in abundance, indicating that dry river bed assemblages are not simply subsets of riparian assemblages. Some spatial patterns in invertebrate assemblages were associated with environmental variables (irrespective of habitat type), but these associations were statistically weak. We suggest that dry river beds are unique habitats in their own right. We discuss potential human stressors and management issues regarding dry river beds and provide recommendations for future research.

Analysis and modelling of the hydraulic conductivity in aquitards : application to the Galilee Basin and the Great Artesian Basin, Australia

In large sedimentary basins with layers of different rocks, the groundwater flow between aquifers depends on the hydraulic conductivity (K) of the separating low-permeable rocks, or aquitards. Three methods were developed to evaluate K in aquitards for areas with limited field data: • Coherence and harmonic analysis: estimates the regional-scale K based on water-level fluctuations in adjacent aquifers. • Cokriging and Bayes' rule: infers K from downhole geophysical logs. • Fluvial process model: reproduces the lithology architecture of sediment formations which can be converted to K. These proposed methods enable good estimates of K and better planning of further drillholes.

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.

Chemical and bioanalytical assessment of coal seam gas associated water

A comprehensive study was undertaken involving chemical (inorganic and organic) and bioanalytical (a suite of 14 in vitro bioassays) assessments of coal seam gas (coal bed methane) associated water (CSGW) in Queensland, Australia. CSGW is a by-product of the gas extraction process and is generally considered as water of poor quality. This was done to better understand what is known about the potential biological and environmental effects associated with the organic constituents of CSGW in Australia. In Queensland, large amounts of associated water must be withdrawn from coal seams to allow extraction of the gas. CSGW is disposed of via release to surface water, reinjected to groundwater or reused for irrigation of crops or pasture, supplied for power station cooling and or reinjected specifically to augment drinking water aquifers. Groundwater samples were collected from private wells tapping into the Walloon Coal Measures, the same coal aquifer exploited for coal seam gas production in the Surat Basin, Australia. The inorganic characteristics of these water samples were almost identical to the CSGW entering the nearby gas company operated Talinga-Condabri Water Treatment Facility. The water is brackish with a pH of 8 to 9, high sodium, bicarbonate and chloride concentrations but low calcium, magnesium and negligible sulphate concentrations. Only low levels of polyaromatic hydrocarbons (PAHs) were detected in the water samples, and neither phenols nor volatile organic compounds were found. Results from the bioassays showed no genotoxicity, protein damage, or activation of hormone receptors (with the exception of the estrogen receptor). However, five of the 14 bioassays gave positive responses: an arylhydrocarbon-receptor gene activation assay (AhR-CAFLUX), estrogenic endocrine activity (ERα-CALUX), oxidative stress response (AREc32), interference with cytokine production (THP1-CPA) and non-specific toxicity (Microtox). The observed effects were benchmarked against known water sources and were similar to secondary treated wastewater effluent, stormwater and surface water. As mixture toxicity modelling demonstrated, the detected PAHs explained less than 5% of the observed biological effects.