Iron biogeochemistry and associated greenhouse gas evolution in a forested subtropical Australian coastal catchment : Poona Creek, Southeast Queensland

Iron (Fe) is the fourth most abundant element in the Earth’s crust. Excess Fe mobilization from terrestrial into aquatic systems is of concern for deterioration of water quality via biofouling and nuisance algal blooms in coastal and marine systems. Substantial Fe dissolution and transport involve alternate Fe(II) oxidation followed by Fe(III) reduction, with a diversity of Bacteria and Archaea acting as the key catalyst. Microbially-mediated Fe cycling is of global significance with regard to cycles of carbon (C), sulfur (S) and manganese (Mn). However, knowledge regarding microbial Fe cycling in circumneutral-pH habitats that prevail on Earth has been lacking until recently. In particular, little is known regarding microbial function in Fe cycling and associated Fe mobilization and greenhouse (CO2 and CH4, GHG) evolution in subtropical Australian coastal systems where microbial response to ambient variations such as seasonal flooding and land use changes is of concern. Using the plantation-forested Poona Creek catchment on the Fraser Coast of Southeast Queensland (SEQ), this research aimed to 1) study Fe cycling-associated bacterial populations in diverse terrestrial and aquatic habitats of a representative subtropical coastal circumneutral-pH (4–7) ecosystem; and 2) assess potential impacts of Pinus plantation forestry practices on microbially-mediated Fe mobilization, organic C mineralization and associated GHG evolution in coastal SEQ. A combination of wet-chemical extraction, undisturbed core microcosm, laboratory bacterial cultivation, microscopy and 16S rRNA-based molecular phylogenetic techniques were employed. The study area consisted primarily of loamy sands, with low organic C and dissolved nutrients. Total reactive Fe was abundant and evenly distributed within soil 0–30 cm profiles. Organic complexation primarily controlled Fe bioavailability and forms in well-drained plantation soils and water-logged, native riparian soils, whereas tidal flushing exerted a strong “seawater effect” in estuarine locations and formed a large proportion of inorganic Fe(III) complexes. There was a lack of Fe(II) sources across the catchment terrestrial system. Mature, first-rotation plantation clear-felling and second-rotation replanting significantly decreased organic matter and poorly crystalline Fe in well-drained soils, although variations in labile soil organic C fractions (dissolved organic C, DOC; and microbial biomass C, MBC) were minor. Both well-drained plantation soils and water-logged, native-vegetation soils were inhabited by a variety of cultivable, chemotrophic bacterial populations capable of C, Fe, S and Mn metabolism via lithotrophic or heterotrophic, (micro)aerobic or anaerobic pathways. Neutrophilic Fe(III)-reducing bacteria (FeRB) were most abundant, followed by aerobic, heterotrophic bacteria (heterotrophic plate count, HPC). Despite an abundance of FeRB, cultivable Fe(II)-oxidizing bacteria (FeOB) were absent in associated soils. A lack of links between cultivable Fe, S or Mn bacterial densities and relevant chemical measurements (except for HPC correlated with DOC) was likely due to complex biogeochemical interactions. Neither did variations in cultivable bacterial densities correlate with plantation forestry practices, despite total cultivable bacterial densities being significantly lower in estuarine soils when compared with well-drained plantation soils and water-logged, riparian native-vegetation soils. Given that bacterial Fe(III) reduction is the primary mechanism of Fe oxide dissolution in soils upon saturation, associated Fe mobilization involved several abiotic and biological processes. Abiotic oxidation of dissolved Fe(II) by Mn appeared to control Fe transport and inhibit Fe dissolution from mature, first-rotation plantation soils post-saturation. Such an effect was not observed in clear-felled and replanted soils associated with low SOM and potentially low Mn reactivity. Associated GHG evolution post-saturation mainly involved variable CO2 emissions, with low, but consistently increasing CH4 effluxes in mature, first-rotation plantation soil only. In comparison, water-logged soils in the riparian native-vegetation buffer zone functioned as an important GHG source, with high potentials for Fe mobilization and GHG, particularly CH4 emissions in riparian loam soils associated with high clay and crystalline Fe fractions. Active Fe–C cycling was unlikely to occur in lower-catchment estuarine soils associated with low cultivable bacterial densities and GHG effluxes. As a key component of bacterial Fe cycling, neutrophilic FeOB widely occurred in diverse aquatic, but not terrestrial, habitats of the catchment study area. Stalked and sheathed FeOB resembling Gallionella and Leptothrix were limited to microbial mat material deposited in surface fresh waters associated with a circumneutral-pH seep, and clay-rich soil within riparian buffer zones. Unicellular, Sideroxydans-related FeOB (96% sequence identity) were ubiquitous in surface and subsurface freshwater environments, with highest abundance in estuary-adjacent shallow coastal groundwater water associated with redox transition. The abundance of dissolved C and Fe in the groundwater-dependent system was associated with high numbers of cultivable anaerobic, heterotrophic FeRB, microaerophilic, putatively lithotrophic FeOB and aerobic, heterotrophic bacteria. This research represents the first study of microbial Fe cycling in diverse circumneutral-pH environments (terrestrial–aquatic, freshwater–estuarine, surface–subsurface) of a subtropical coastal ecosystem. It also represents the first study of its kind in the southern hemisphere. This work highlights the significance of bacterial Fe(III) reduction in terrestrial, and bacterial Fe(II) oxidation in aquatic catchment Fe cycling. Results indicate the risk of promotion of Fe mobilization due to plantation clear-felling and replanting, and GHG emissions associated with seasonal water-logging. Additional significant outcomes were also achieved. The first direct evidence for multiple biomineralization patterns of neutrophilic, microaerophilic, unicellular FeOB was presented. A putatively pure culture, which represents the first cultivable neutrophilic FeOB from the southern hemisphere, was obtained as representative FeOB ubiquitous in diverse catchment aquatic habitats.

Three-dimensional geological modelling and multivariate statistical analysis of water chemistry data to analyse and visualise aquifer structure and groundwater composition in the Wairau Plain, Marlborough District, New Zealand

Concerns regarding groundwater contamination with nitrate and the long-term sustainability of groundwater resources have prompted the development of a multi-layered three dimensional (3D) geological model to characterise the aquifer geometry of the Wairau Plain, Marlborough District, New Zealand. The 3D geological model which consists of eight litho-stratigraphic units has been subsequently used to synthesise hydrogeological and hydrogeochemical data for different aquifers in an approach that aims to demonstrate how integration of water chemistry data within the physical framework of a 3D geological model can help to better understand and conceptualise groundwater systems in complex geological settings. Multivariate statistical techniques(e.g. Principal Component Analysis and Hierarchical Cluster Analysis) were applied to groundwater chemistry data to identify hydrochemical facies which are characteristic of distinct evolutionary pathways and a common hydrologic history of groundwaters. Principal Component Analysis on hydrochemical data demonstrated that natural water-rock interactions, redox potential and human agricultural impact are the key controls of groundwater quality in the Wairau Plain. Hierarchical Cluster Analysis revealed distinct hydrochemical water quality groups in the Wairau Plain groundwater system. Visualisation of the results of the multivariate statistical analyses and distribution of groundwater nitrate concentrations in the context of aquifer lithology highlighted the link between groundwater chemistry and the lithology of host aquifers. The methodology followed in this study can be applied in a variety of hydrogeological settings to synthesise geological, hydrogeological and hydrochemical data and present them in a format readily understood by a wide range of stakeholders. This enables a more efficient communication of the results of scientific studies to the wider community.

Limestone and speleothem trace element geochemistry as tools for palaeoclimatic reconstruction, Mount Etna region, central-coastal Queensland

This study investigated potential palaeoclimate proxies provided by rare earth element (REE) geochemistry in speleothems and in clay mineralogy of cave sediments. Speleothem and sediment samples were collected from a series of cave fill deposits that occurred with rich vertebrate fossil assemblages in and around Mount Etna National Park, Rockhampton (central coastal Queensland). The fossil deposits range from Plio- Pleistocene to Holocene in age (based on uranium/thorium dating) and appear to represent depositional environments ranging from enclosed rainforest to semi-arid grasslands. Therefore, the Mount Etna cave deposits offer the perfect opportunity to test new palaeoclimate tools as they include deposits that span a known significant climate shift on the basis of independent faunal data. The first section of this study investigates the REE distribution of the host limestone to provide baseline geochemistry for subsequent speleothem investigations. The Devonian Mount Etna Beds were found to be more complex than previous literature had documented. The studied limestone massif is overturned, highly recrystallised in parts and consists of numerous allochthonous blocks with different spatial orientations. Despite the complex geologic history of the Mount Etna Beds, Devonian seawater-like REE patterns were recovered in some parts of the limestone and baseline geochemistry was determined for the bulk limestone for comparison with speleothem REE patterns. The second part of the study focused on REE distribution in the karst system and the palaeoclimatic implications of such records. It was found that REEs have a high affinity for calcite surfaces and that REE distributions in speleothems vary between growth bands much more than along growth bands, thus providing a temporal record that may relate to environmental changes. The morphology of different speleothems (i.e., stalactites, stalagmites, and flowstones) has little bearing on REE distributions provided they are not contaminated with particulate fines. Thus, baseline knowledge developed in the study suggested that speleothems were basically comparable for assessing palaeoclimatically controlled variations in REE distributions. Speleothems from rainforest and semi-arid phases were compared and it was found that there are definable differences in REE distribution that can be attributed to climate. In particular during semiarid phases, total REE concentration decreased, LREE became more depleted, Y/Ho increased, La anomalies were more positive and Ce anomalies were more negative. This may reflect more soil development during rainforest phases and more organic particles and colloids, which are known to transport REEs, in karst waters. However, on a finer temporal scale (i.e. growth bands) within speleothems from the same climate regime, no difference was seen. It is suggested that this may be due to inadequate time for soil development changes on the time frames represented by differences in growth band density. The third part of the study was a reconnaissance investigation focused on mineralogy of clay cave sediments, illite/kaolinite ratios in particular, and the potential palaeoclimatic implications of such records. Although the sample distribution was not optimal, the preliminary results suggest that the illite/kaolinite ratio increased during cold and dry intervals, consistent with decreased chemical weathering during those times. The study provides a basic framework for future studies at differing latitudes to further constrain the parameters of the proxy. The identification of such a proxy recorded in cave sediment has broad implications as clay ratios could potentially provide a basic local climate proxy in the absence of fossil faunas and speleothem material. This study suggests that REEs distributed in speleothems may provide information about water throughput and soil formation, thus providing a potential palaeoclimate proxy. It highlights the importance of understanding the host limestone geochemistry and broadens the distribution and potential number of cave field sites as palaeoclimate information no longer relies solely on the presence of fossil faunas and or speleothems. However, additional research is required to better understand the temporal scales required for the proxies to be recognised.

Microcosm study of ion mobilization and greenhouse gas evolution in soils of a plantation-forested subtropical coastal catchment

This study examined the potential for Fe mobilization and greenhouse gas (GHG, e.g. CO2, and CH4) evolution in SEQ soils associated with a range of plantation forestry practices and water-logged conditions. Intact, 30-cm-deep soil cores collected from representative sites were saturated and incubated for 35 days in the laboratory, with leachate and headspace gas samples periodically collected. Minimal Fe dissolution was observed in well-drained sand soils associated with mature, first-rotation Pinus and organic Fe complexation, whereas progressive Fe dissolution occurred over 14 days in clear-felled and replanted Pinus soils with low organic matter and non-crystalline Fe fractions. Both CO2 and CH4 effluxes were relatively lower in clear-felled and replanted soils compared with mature, first-rotation Pinus soils, despite the lack of statistically significant variations in total GHG effluxes associated with different forestry practices. Fe dissolution and GHG evolution in low-lying, water-logged soils adjacent to riparian and estuarine, native-vegetation buffer zones were impacted by mineral and physical soil properties. Highest levels of dissolved Fe and GHG effluxes resulted from saturation of riparian loam soils with high Fe and clay content, as well as abundant organic material and Fe-metabolizing bacteria. Results indicate Pinus forestry practices such as clear-felling and replanting may elevate Fe mobilization while decreasing CO2 and CH4 emissions from well-drained, SEQ plantation soils upon heavy flooding. Prolonged water-logging accelerates bacterially mediated Fe cycling in low-lying, clay-rich soils, leading to substantial Fe dissolution, organic matter mineralization, and CH4 production in riparian native-vegetation buffer zones.

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.

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.

Preparing for coastal change

Coastal areas are dynamic environments that are home to billions of people worldwide and provide areas of unique natural importance. As such, coastal change is of considerable local and global interest, not only within the geological realm, but also in terms of socioeconomic and biodiversity impacts. An accurate understanding of how changes in relative sea level, geological processes and extreme events, such as storms and tsunamis, have interacted to shape and change the Earth’s coastlines over millennia is fundamental to future projections of coastal change. On the basis of this, researchers in these, and various other aspects of coastal change were brought together in late 2010 at the University of Hong Kong for the first meeting of International Geoscience Program Project 588 (IGCP588) e Preparing for Coastal Change. This special issue showcases some of the results presented at this meeting.

Increasing disaster resilience of people and property through fostering collaboration in coastal high hazard zones

As evidenced with the 2011 floods the state of Queensland in Australia is quite vulnerable to this kind of disaster. Climate change will increase the frequency and magnitude of such events and will have a variety of other impacts. To deal with these governments at all levels need to be prepared and work together. Since most of the population of the state is located in the coastal areas and these areas are more vulnerable to the impacts of climate change this paper examines climate change adaptation efforts in coastal Queensland. The paper is part of a more comprehensive project which looks at the critical linkages between land use and transport planning in coastal Queensland, especially in light of increased frequencies of cyclonic activity and other impacts associated with climate change. The aim is improving coordination between local and state government in addressing land use and transport planning in coastal high hazard areas. By increasing the ability of local governments and state agencies to coordinate planning activities, we can help adapt to impacts of climate change. Towards that end, we will look at the ways that these groups currently interact, especially with regard to issues involving uncertainty related to climate change impacts. Through surveys and interviews of Queensland coastal local governments and state level planning agencies on how they coordinate their planning activities at different levels as well as how much they take into account the linkage of transportation and land use we aim to identify the weaknesses of the current planning system in responding to the challenges of climate change adaptation. The project will identify opportunities for improving the ways we plan and coordinate planning, and make recommendations to improve resilience in advance of disasters so as to help speed up recovery when they occur.

Linking three-dimensional geological modelling and multivariate statistical analysis to define groundwater chemistry baseline and identify inter-aquifer connectivity within the Clarence-Moreton Basin, Southeast Queensland, Australia

The Clarence-Moreton Basin (CMB) covers approximately 26000 km2 and is the only sub-basin of the Great Artesian Basin (GAB) in which there is flow to both the south-west and the east, although flow to the south-west is predominant. In many parts of the basin, including catchments of the Bremer, Logan and upper Condamine Rivers in southeast Queensland, the Walloon Coal Measures are under exploration for Coal Seam Gas (CSG). In order to assess spatial variations in groundwater flow and hydrochemistry at a basin-wide scale, a 3D hydrogeological model of the Queensland section of the CMB has been developed using GoCAD modelling software. Prior to any large-scale CSG extraction, it is essential to understand the existing hydrochemical character of the different aquifers and to establish any potential linkage. To effectively use the large amount of water chemistry data existing for assessment of hydrochemical evolution within the different lithostratigraphic units, multivariate statistical techniques were employed.

Disaster resilience and local adaptation capacity in coastal Queensland

The 2011 floods illustrated once again Queensland’s vulnerability to flooding and similar disasters. Climate change will increase the frequency and magnitude of such events and will have a variety of other impacts. To deal with these impacts governments at all levels need to be prepared and work together. Like the rest of the nation most of the population of the state is located in the coastal areas and these areas are more vulnerable to the impacts of climate change. This paper examines climate change adaptation efforts in coastal Queensland. The aim is increasing local disaster resilience of people and property through fostering coordination between local and state government planning activities in coastal high hazard areas. By increasing the ability of local governments and state agencies to coordinate planning activities, we can help adapt to impacts of climate change. Towards that end, we will look at the ways that these groups currently interact, especially with regard to issues involving uncertainty related to climate change impacts. Through an examination of climate change related activities by Queensland’s coastal local governments and state level planning agencies and how they coordinate their planning activities at different levels we aim to identify the weaknesses of the current planning system in responding to the challenges of climate change adaptation and opportunities for improving the ways we plan and coordinate planning, and make recommendations to improve resilience in advance of disasters so as to help speed up recovery when they occur.

Modelling surface and shallow groundwater interactions in an ungauged subtropical coastal catchment using the SWAT model, Elimbah Creek, southeast Queensland, Australia

The study presented here applies the highly parameterised semi-distributed U.S. Department of Agriculture Soil and Water Assessment Tool (SWAT) to an Australian subtropical catchment. SWAT has been applied to numerous catchments worldwide and is considered to be a useful tool that is under ongoing development with contributions coming from different research groups in different parts of the world. In a preliminary run the SWAT model application for the Elimbah Creek catchment has estimated water yield for the catchment and has quantified the different sources. For the modelling period of April 1999 to September 2009 the results show that the main sources of water in Elimbah Creek are total surface runoff and lateral flow (65%). Base-flow contributes 36% to the total runoff. On a seasonal basis modelling results show a shift in the source of water contributing to Elimbah Creek from surface runoff and lateral flow during intense summer storms to base-flow conditions during dry months. Further calibration and validation of these results will confirm that SWAT provides an alternative to Australian water balance models.

Integrated coastal zone management and sustainable development of coastal area : a short overview of international legal framework

This chapter contains sections titled: Introduction ICZM and sustainable development of coastal zone International legal framework for ICZM Implementation of international legal obligations in domestic arena Concluding remarks References

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.