Logan Basin Water Resource Plan : environmental assessment of the southern Moreton Bay islands. Technical Advisory Panel Report: Stage 1.

This report is a technical assessment of the hydrological environment of the southern Moreton Bay islands and follows the terms of reference supplied by the then Queensland Department of Natural Resources and Water. The terms of reference describe stage 1 as a condition assessment and stage 2 as an assessment of the implications of water planning scenarios on future condition. This report is the first stage of a two-stage investigation whose primary purpose is to identify and assess groundwater dependent ecosystems (GDEs) and the groundwater flow regimes necessary to support them. Within this context, the groundwaters themselves are also considered and comment made on their condition. Information provided in this report will inform an amendment to the Logan Basin Water Resource Plan to incorporate the southern Moreton Bay islands. The study area is the water resource plan amendment area, which includes North and South Stradbroke islands and the smaller islands between these and the mainland, including the inhabited smaller rocky islands—namely, Macleay, Russell, Karragarra, Lamb and Coochiemudlo islands. This assessment is largely a desktop study based on existing information, but incorporates some field observations, input from experts in specific areas and community representatives, and the professional experience and knowledge of the authors. This report reviews existing research and information on the southern Moreton Bay area with an emphasis on North Stradbroke Island, as it represents the largest and most regionally significant groundwater resource in southern Moreton Bay. The report provides an assessment of key waterrelated environmental features, their condition and their degree of dependence on groundwater. This report also assesses the condition and status of ecosystems within this region. In addition, the report identifies information gaps, uncertainties and potential impacts; reviews groundwater models that have been developed for North Stradbroke Island; and makes recommendations on monitoring and research needs.

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.

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.

A case study of interactions between flood flow and buildings in an urban environment : Gardens Point during the 12-13 January 2011 flood of the Brisbane River (Australia)

Flood flows in inundated urban environment constitute a natural hazard. During the 12- 13 January 2011 flood of the Brisbane River, detailed water elevation, velocity and suspended sediment data were recorded in an inundated street at the peak of the flood. The field observations highlighted a number of unusual flow interactions with the urban surroundings. These included some slow fluctuations in water elevations and velocity with distinctive periods between 50 and 100 s caused by some local topographic effect (choking), superposed with some fast turbulent fluctuations. The suspended sediment data highlighted some significant suspended sediment loads in the inundated zone.

IODP Expedition 339 in the Gulf of Cadiz and off West Iberia : decoding the environmental significance of the Mediterranean outflow water and its global influence

IODP Expedition 339 drilled five sites in the Gulf of Cadiz and two off the west Iberian margin (November 2011 to January 2012), and recovered 5.5 km of sediment cores with an average recovery of 86.4%. The Gulf of Cadiz was targeted for drilling as a key location for the investigation of Mediterranean outflow water (MOW) through the Gibraltar Gateway and its influence on global circulation and climate. It is also a prime area for understanding the effects of tectonic activity on evolution of the Gibraltar Gateway and on margin sedimentation. We penetrated into the Miocene at two different sites and established a strong signal of MOW in the sedimentary record of the Gulf of Cadiz, following the opening of the Gibraltar Gateway. Preliminary results show the initiation of contourite deposition at 4.2–4.5 Ma, although subsequent research will establish whether this dates the onset of MOW. The Pliocene succession, penetrated at four sites, shows low bottom current activity linked with a weak MOW. Significant widespread unconformities, present in all sites but with hiatuses of variable duration, are interpreted as a signal of intensified MOW, coupled with flow confinement. The Quaternary succession shows a much more pronounced phase of contourite drift development, with two periods of MOW intensification separated by a widespread unconformity. Following this, the final phase of drift evolution established the contourite depositional system (CDS) architecture we see today. There is a significant climate control on this evolution of MOW and bottom-current activity. However, from the closure of the Atlantic–Mediterranean gateways in Spain and Morocco just over 6 Ma and the opening of the Gibraltar Gateway at 5.3 Ma, there has been an even stronger tectonic control on margin development, downslope sediment transport and contourite drift evolution. The Gulf of Cadiz is the world's premier contourite laboratory and thus presents an ideal testing ground for the contourite paradigm. Further study of these contourites will allow us to resolve outstanding issues related to depositional processes, drift budgets, and recognition of fossil contourites in the ancient record on shore. The expedition also verified an enormous quantity and extensive distribution of contourite sands that are clean and well sorted. These represent a relatively untapped and important exploration target for potential oil and gas reservoirs.

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.

Characterization of a DC-driven microplasma between a capillary tube and water surface

A microplasma generated between a stainless-steel capillary and water surface in ambient air with flowing argon as working gas appears as a bright spot at the tube orifice and expands to form a larger footprint on the water surface, and the dimensions of the bell-shaped microplasma are all below 1 mm. The electron density of the microplasma is estimated to be ranging from 5.32 × 109 cm−3 to 2.02 × 1014 cm−3 for the different operating conditions, which is desirable for generating abundant amounts of reactive species. A computational technique is adopted to fit the experimental emission from the N2 second positive system with simulation results. It is concluded that the vibrational temperature (more than 2000 K) is more than twice the gas temperature (more than 800 K), which indicates the non-equilibrium state of the microplasma. Both temperatures showed dependence on the discharge parameters (i.e., gas flow and discharge current). Such a plasma device could be arranged in arrays for applications utilizing plasmainduced liquid chemistry.

CFD convective flow simulation of the varying properties of CO2-H20 mixtures in geothermal systems

Numerical simulation of a geothermal reservoir, modelled as a bottom-heated square box, filled with water-CO2 mixture is presented in this work. Furthermore, results for two limiting cases of a reservoir filled with either pure water or CO2 are presented. Effects of different parameters including CO2 concentration as well as reservoir pressure and temperature on the overall performance of the system are investigated. It has been noted that, with a fixed reservoir pressure and temperature, any increase in CO2concentration leads to better performance, i.e. stronger convection and higher heat transfer rates. With a fixed CO2 concentration, however, the reservoir pressure and temperature can significantly affect the overall heat transfer and flow rate from the reservoir. Details of such variations are documented and discussed in the present paper.

A robust methodology for developing an operational systems view of mine water interactions

The Minerals Council of Australia’s (MCA) Water Accounting Framework (WAF) is an industry lead initiative to enable cross company communication and comparisons of water management performance. The WAF consists of two models, the Input-Output Model that represents water interactions between an operation and its surrounding environment and the Operational Model that represents water interactions within an operation. Recently, MCA member companies have agreed to use the Input-Output Model to report on their external water interactions in Australian operations, with some adopting it globally. The next step will be to adopt the Operational Model. This will expand the functionality of the WAF from corporate reporting to allowing widespread identification of inefficiencies and to connect internal and external interactions. Implementing the WAF, particularly the Operational Model, is non-trivial. It can be particularly difficult for operations that are unfamiliar with the WAF definitions and methodology, lack information pertaining to flow volumes or contain unusual configurations. Therefore, there is a need to help industry with its implementation. This work presents a step-by-step guide to producing the Operational Model. It begins by describing a methodology for implementing the Operational Model by describing the identification of pertinent objects (stores, tasks and treatments), quantification of flows, aggregation of objects and production of reports. It then discusses how the Operational Model can represent a series of challenging scenarios and how it can be connected with Input-Output Model to improve water management.

Onset of transition in mixed convection of a lid-driven trapezoidal enclosure filled with water-Al2O3 nanofluid

A numerical study is carried out to investigate the transition from laminar to chaos in mixed convection heat transfer inside a lid-driven trapezoidal enclosure. In this study, the top wall is considered as isothermal cold surface, which is moving in its own plane at a constant speed, and a constant high temperature is provided at the bottom surface. The enclosure is assumed to be filled with water-Al2O3 nanofluid. The governing Navier–Stokes and thermal energy equations are expressed in non-dimensional forms and are solved using Galerkin finite element method. Attention is paid in the present study on the pure mixed convection regime at Richandson number, Ri = 1. The numerical simulations are carried out over a wide range of Reynolds (0.1 ≤ Re ≤ 103) and Grashof (0.01 ≤ Gr ≤ 106) numbers. Effects of the presence of nanofluid on the characteristics of mixed convection heat transfer are also explored. The average Nusselt numbers of the heated wall are computed to demonstrate the influence of flow parameter variations on heat transfer. The corresponding change of flow and thermal fields is visualized from the streamline and the isotherm contour plots.

Hydrochemical evolution and groundwater flow processes in the Galilee and Eromanga basins, Great Artesian Basin, Australia: A multivariate statistical approach

The Galilee and Eromanga basins are sub-basins of the Great Artesian Basin (GAB). In this study, a multivariate statistical approach (hierarchical cluster analysis, principal component analysis and factor analysis) is carried out to identify hydrochemical patterns and assess the processes that control hydrochemical evolution within key aquifers of the GAB in these basins. The results of the hydrochemical assessment are integrated into a 3D geological model (previously developed) to support the analysis of spatial patterns of hydrochemistry, and to identify the hydrochemical and hydrological processes that control hydrochemical variability. In this area of the GAB, the hydrochemical evolution of groundwater is dominated by evapotranspiration near the recharge area resulting in a dominance of the Na–Cl water types. This is shown conceptually using two selected cross-sections which represent discrete groundwater flow paths from the recharge areas to the deeper parts of the basins. With increasing distance from the recharge area, a shift towards a dominance of carbonate (e.g. Na–HCO3 water type) has been observed. The assessment of hydrochemical changes along groundwater flow paths highlights how aquifers are separated in some areas, and how mixing between groundwater from different aquifers occurs elsewhere controlled by geological structures, including between GAB aquifers and coal bearing strata of the Galilee Basin. The results of this study suggest that distinct hydrochemical differences can be observed within the previously defined Early Cretaceous–Jurassic aquifer sequence of the GAB. A revision of the two previously recognised hydrochemical sequences is being proposed, resulting in three hydrochemical sequences based on systematic differences in hydrochemistry, salinity and dominant hydrochemical processes. The integrated approach presented in this study which combines different complementary multivariate statistical techniques with a detailed assessment of the geological framework of these sedimentary basins, can be adopted in other complex multi-aquifer systems to assess hydrochemical evolution and its geological controls.

4D GPR tracking of water infiltration in fractured high-porosity limestone

Three thousand liters of water were infiltrated from a 4 m diameter pond to track flow and transport inside fractured carbonates with 20-40 % porosity. Sixteen time-lapse 3D Ground Penetrating Radar (GPR) surveys with repetition intervals between 2 hrs and 5 days monitored the spreading of the water bulb in the subsurface. Based on local travel time shifts between repeated GPR survey pairs, localized changes of volumetric water content can be related to the processes of wetting, saturation and drainage. Deformation bands consisting of thin sub vertical sheets of crushed grains reduce the magnitude of water content changes but enhance flow in sheet parallel direction. This causes an earlier break through across a stratigraphic boundary compared to porous limestone without deformation bands. This experiment shows how time-lapse 3D GPR or 4D GPR can non-invasively track ongoing flow processes in rock-volumes of over 100 m3.

High-resolution GNSS-tracked drifter for studying surface dispersion in shallow water

The use of GNSS tracked Lagrangian drifters allows more realistic quantification of fluid motion and dispersion coefficients than Eulerian techniques because such drifters are analogues of particles that are relevant to flow field characterisation and pollutant dispersion. Using the fast growing Real Time Kinematic (RTK) positioning technique derived from Global Satellite Navigation Systems (GNSS), drifters are developed for high frequency (10 Hz) sampling with position estimates to centimetre accuracy. The drifters are designed with small size and less direct wind drag to follow the sub-surface flow which characterizes dispersion in shallow waters. An analysis of position error from stationary observation indicates that the drifter can efficiently resolve motion up to 1 Hz. The result of the field deployments of the drifter in conjunction with acoustic Eulerian devices shows higher estimate of the drifter streamwise velocities. Single particle statistical analysis of field deployments in a shallow estuarine zone yielded dispersion coefficients estimate comparable to those of dye tracer studies. The drifters capture the tidal elevation during field studies in a tidal estuary.

Estimate of Lagrangian integral scales in shallow tidal water using high resolution GPS-tracked drifters

In estuaries and natural water channels, the estimate of velocity and dispersion coefficients is critical to the knowledge of scalar transport and mixing. This estimate is rarely available experimentally at sub-tidal time scale in shallow water channels where high frequency is required to capture its spatio-temporal variation. This study estimates Lagrangian integral scales and autocorrelation curves, which are key parameters for obtaining velocity fluctuations and dispersion coefficients, and their spatio-temporal variability from deployments of Lagrangian drifters sampled at 10 Hz for a 4-hour period. The power spectral densities of the velocities between 0.0001 and 0.8 Hz were well fitted with a slope of 5/3 predicted by Kolmogorov’s similarity hypothesis within the inertial subrange, and were similar to the Eulerian power spectral previously observed within the estuary. The result showed that large velocity fluctuations determine the magnitude of the integral time scale, TL. Overlapping of short segments improved the stability of the estimate of TL by taking advantage of the redundant data included in the autocorrelation function. The integral time scales were about 20 s and varied by up to a factor of 8. These results are essential inputs for spatial binning of velocities, Lagrangian stochastic modelling and single particle analysis of the tidal estuary.

Effect of hydrodynamics and surface roughness on the electrochemical behaviour of carbon steel in CSG produced water

The influence of fluid flow, surface roughness and immersion time on the electrochemical behaviour of carbon steel in coal seam gas produced water under static and hydrodynamic conditions has been studied. The disc electrode surface morphology before and after the corrosion test was characterized using scanning electron microscopy (SEM). The corrosion product was examined using X-ray photoelectron spectroscopy (XPS) and X-ray diffractometry (XRD).The results show that the anodic current density increased with increasing surface roughness and consequently a decrease in corrosion surface resistance. Under dynamic flow conditions, the corrosion rate increased with increasing rotating speed due to the high mass transfer coefficient and formation of non-protective akaganeite β- FeO(OH) and goethite α- FeO(OH) corrosion scale at the electrode surface.The corrosion rate was lowest at 0 rpm.The corrosion rate decreased in both static and dynamic conditions with increasing immersion time. The decrease in corrosion rate is attributed to the deposition of corrosion products on the electrode surface. SEM results revealed that the rougher surface exhibited a great tendency toward pitting corrosion.