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.

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

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

Controls on mid-Holocene fringing reef growth and termination in a high latitude, estuarine setting, Wellington Point, Southeast Queensland

Several fringing coral reefs in Moreton Bay, Southeast Queensland, some 300 km south of the Great Barrier Reef (GBR), are set in a relatively high latitude, estuarine environment that is considered marginal for coral growth. Previous work indicated that these marginal reefs, as with many fringing reefs of the inner GBR, ceased accreting in the mid-Holocene. This research presents for the first time data from the subsurface profile of the mid-Holocene fossil reef at Wellington Point comprising U/Th dates of in situ and framework corals, and trace element analysis from the age constrained carbonate fragments. Based on trace element proxies the palaeo-water quality during reef accretion was reconstructed. Results demonstrate that the reef initiated more than 7,000 yr BP during the post glacial transgression, and the initiation progressed to the west as sea level rose. In situ micro-atolls indicate that sea level was at least 1 m above present mean sea level by 6,680 years ago. The reef remained in "catch-up" mode, with a seaward sloping upper surface, until it stopped aggrading abruptly at ca 6,000 yr BP; no lateral progradation occurred. Changes in sediment composition encountered in the cores suggest that after the laterite substrate was covered by the reef, most of the sediment was produced by the carbonate factory with minimal terrigenous influence. Rare earth element, Y and Ba proxies indicate that water quality during reef accretion was similar to oceanic waters, considered suitable for coral growth. A slight decline in water quality on the basis of increased Ba in the later stages of growth may be related to increased riverine input and partial closing up of the bay due to either tidal delta progradation, climatic change and/or slight sea level fall. The age data suggest that termination of reef growth coincided with a slight lowering of sea level, activation of ENSO and consequent increase in seasonality, lowering of temperatures and the constrictions to oceanic flushing. At the cessation of reef accretion the environmental conditions in the western Moreton Bay were changing from open marine to estuarine. The living coral community appears to be similar to the fossil community, but without the branching Acropora spp. that were more common in the fossil reef. In this marginal setting coral growth periods do not always correspond to periods of reef accretion due to insufficient coral abundance. Due to several environmental constraints modern coral growth is insufficient for reef growth. Based on these findings Moreton Bay may be unsuitable as a long term coral refuge for most species currently living in the GBR.

Mid-Holocene sea-level and coral reef demise : U-Th dating of subfossil corals in Moreton Bay, Australia

It is increasingly apparent that sea-level data (e.g. microfossil transfer functions, dated coral microatolls and direct observations from satellite and tidal gauges) vary temporally and spatially at regional to local scales, thus limiting our ability to model future sea-level rise for many regions. Understanding sealevel response at ‘far-field’ locations at regional scales is fundamental for formulating more relevant sea-level rise susceptibility models within these regions under future global change projections. Fossil corals and reefs in particular are valuable tools for reconstructing past sea levels and possible environmental phase shifts beyond the temporal constraints of instrumental records. This study used abundant surface geochronological data based on in situ subfossil corals and precise elevation surveys to determine previous sea level in Moreton Bay, eastern Australia, a far-field site. A total of 64 U-Th dates show that relative sea level was at least 1.1 m above modern lowest astronomical tide (LAT) from at least ˜6600 cal. yr BP. Furthermore, a rapid synchronous demise in coral reef growth occurred in Moreton Bay ˜5800 cal. yr BP, coinciding with reported reef hiatus periods in other areas around the Indo-Pacific region. Evaluating past reef growth patterns and phases allows for a better interpretation of anthropogenic forcing versus natural environmental/climatic cycles that effect reef formation and demise at all scales and may allow better prediction of reef response to future global change.

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.

An analytical framework for quantifying aquifer response time scales associated with transient boundary conditions

A major challenge in studying coupled groundwater and surface-water interactions arises from the considerable difference in the response time scales of groundwater and surface-water systems affected by external forcings. Although coupled models representing the interaction of groundwater and surface-water systems have been studied for over a century, most have focused on groundwater quantity or quality issues rather than response time. In this study, we present an analytical framework, based on the concept of mean action time (MAT), to estimate the time scale required for groundwater systems to respond to changes in surface-water conditions. MAT can be used to estimate the transient response time scale by analyzing the governing mathematical model. This framework does not require any form of transient solution (either numerical or analytical) to the governing equation, yet it provides a closed form mathematical relationship for the response time as a function of the aquifer geometry, boundary conditions, and flow parameters. Our analysis indicates that aquifer systems have three fundamental time scales: (i) a time scale that depends on the intrinsic properties of the aquifer; (ii) a time scale that depends on the intrinsic properties of the boundary condition, and; (iii) a time scale that depends on the properties of the entire system. We discuss two practical scenarios where MAT estimates provide useful insights and we test the MAT predictions using new laboratory-scale experimental data sets.

The Cressbrook Creek alluvial aquifer system, Southeast Queensland : hydrochemistry and isotopes to determine hydrological processes and response to floods

This study developed an understanding of hydrological processes within the Cressbrook Creek catchment of the upper Brisbane River, in particular for the alluvial aquifers. Those aquifers within the lower catchment are used for intensive irrigation, and have been impacted by long-term drought followed by flooding. The study utilised water chemistry, isotopic characters and hydraulic measurements to determine factors such as recharge, links between creeks and groundwater, and variations in water quality. The catchment-wide study will enable improved management of the local water resources.

Three-dimensional geological modelling of the Galilee and central Eromanga basins, Australia: New insights into aquifer/aquitard geometry and potential influence of faults on inter-connectivity

Study region The Galilee and Eromanga basins are located in central Queensland, Australia. Both basins are components of the Great Artesian Basin which host some of the most significant groundwater resources in Australia. Study focus This study evaluates the influence of regional faults on groundwater flow in an aquifer/aquitard interbedded succession that form one of the largest Artesian Basins in the world. In order to assess the significance of regional faults as potential barriers or conduits to groundwater flow, vertical displacements of the major aquifers and aquitards were studied at each major fault and the general hydraulic relationship of units that are juxtaposed by the faults were considered. A three-dimensional (3D) geological model of the Galilee and Eromanga basins was developed based on integration of well log data, seismic surfaces, surface geology and elevation data. Geological structures were mapped in detail and major faults were characterised. New hydrological insights for the region Major faults that have been described in previous studies have been confirmed within the 3D geological model domain and a preliminary assessment of their hydraulic significance has been conducted. Previously unknown faults such as the Thomson River Fault (herein named) have also been identified in this study.

Mississippi Delta subsidence primarily caused by compaction of Holocene strata

Coastal subsidence causes sea-level rise, shoreline erosion and wetland loss, which poses a threat to coastal populations. This is especially evident in the Mississippi Delta in the southern United States, which was devastated by Hurricane Katrina in 2005. The loss of protective wetlands is considered a critical factor in the extensive flood damage. The causes of subsidence in coastal Louisiana, attributed to factors as diverse as shallow compaction and deep crustal processes, remain controversial. Current estimates of subsidence rates vary by several orders of magnitude. Here, we use a series of radiocarbon-dated sediment cores from the Mississippi Delta to analyse late Holocene deposits and assess compaction rates. We find that millennial-scale compaction rates primarily associated with peat can reach 5mm per year, values that exceed recent model predictions. Locally and on timescales of decades to centuries, rates are likely to be 10 mm or more per year. We conclude that compaction of Holocene strata contributes significantly to the exceptionally high rates of relative sea-level rise and coastal wetland loss in the Mississippi Delta, and is likely to cause subsidence in other organic-rich and often densely populated coastal plains.

Mid-Holocene Aboriginal occupation of offshore islands in northern Australia? A reassessment of Wurdukanhan, Mornington Island, southern Gulf of Carpentaria, Australia

Claims for mid-Holocene Aboriginal occupation at the shell matrix site of Wurdukanhan, Mornington Island, Gulf of Carpentaria, Australia, are reassessed through an analysis of the excavated assemblage coupled with new surveys and an extensive dating program. Memmott et al. (2006, pp. 38, 39) reported basal ages of c.5000–5500 years from Wurdukanhan as 'the oldest date yet obtained for any archaeological site on the coast of the southern Gulf of Carpentaria' and used these dates to argue for 'a relatively lengthy occupation since at least the mid-Holocene'. If substantiated, with the exception of western Torres Strait, these claims make Mornington Island the only offshore island used across northern Australia in the mid-Holocene where it is conventionally thought that Aboriginal people only (re)colonised islands after sea-level maximum was achieved after the mid-Holocene. Our analysis of Wurdukanhan demonstrates high shellfish taxa diversity, high rates of natural shell predation and high densities of foraminifera throughout the deposit demonstrating a natural origin for the assemblage. Results are considered in the context of other dated shell matrix sites in the area and a geomorphological model for landscape development of the Sandalwood River catchment.

Environmental context for late Holocene human occupation of the South Wellesley Archipelago, Gulf of Carpentaria, northern Australia

A 2400 year record of environmental change is reported from a wetland on Bentinck Island in the southern Gulf of Carpentaria, northern Australia. Three phases of wetland development are identified, with a protected coastal setting from ca. 2400 to 500 years ago, transitioning into an estuarine mangrove forest from ca. 500 years ago to the 1940s, and finally to a freshwater swamp over the past +60 years. This sequence reflects the influence of falling sea-levels, development of a coastal dune barrier system, prograding shorelines, and an extreme storm (cyclone) event. In addition, there is clear evidence of the impacts that human abandonment and resettlement have on the island's fire regimes and vegetation. A dramatic increase in burning and vegetation thickening was observed after the cessation of traditional Indigenous Kaiadilt fire management practices in the 1940s, and was then reversed when people returned to the island in the 1980s. In terms of the longer context for human occupation of the South Wellesley Archipelago, it is apparent that the mangrove phase provided a stable and productive environment that was conducive for human settlement of this region over the past 1000 years.

Large variations in the Holocene marine radiocarbon reservoir effect reflect ocean circulation and climatic changes

Accurate radiocarbon dating of marine samples requires knowledge of the marine radiocarbon reservoir effect. This effect for a particular site/region is generally assumed constant through time when calibrating marine 14C ages. However, recent studies have shown large temporal variations of several hundred to a couple of thousand years in this effect for a number of regions during the late Quaternary and Holocene. Here we report marine radiocarbon reservoir correction (ΔRΔR) for Heron Reef and Moreton Bay in southwestern (SW) Pacific for the last 8 ka derived from 14C analysis of 230Th-dated corals. Most of our ΔRΔR for the last ∼5.4 ka agree well with their modern value, but large ΔRΔR variability of ∼410 yr (from trough to peak) with possible decadal/centennial fluctuations is evident for the period ∼5.4–8 ka. The latter time interval also has significant variations with similar features in previously published ΔRΔR values for other sites in the Pacific, including southern Peru–northern Chile in southeastern (SE) Pacific, the South China Sea, Vanuatu and Papua New Guinea, with the largest magnitude of ∼920 yr from SE Pacific. The mechanisms for these large ΔRΔR variations across the Pacific during the mid-Holocene are complex processes involving (1) changes in the quantity and 14C content of upwelled waters in tropical east Pacific (TEP) (frequency and intensity of ocean upwelling in the TEP, and contribution of Subantarctic Mode Water to the upwelled waters, which is influenced by the intensity and position of southern westerly winds), and (2) variations in ocean circulation associated with climate change (La Niña/El Niño conditions, intensity of easterly trade winds, positions of the Intertropical Convergence Zone and the South Pacific Convergence Zone), which control the spreading of the older upwelled surface waters in the TEP to the western sites. Our results imply the need for employing temporal changes in ΔRΔR values, instead of constant (modern) values, for age calibration of Holocene marine samples not only for the SW Pacific sites but also for other tropical and subtropical sites in the Pacific.

Seasonal and spatial variations in rare earth elements to identify inter-aquifer linkages and recharge processes in an Australian catchment

With the aim of elucidating the seasonal behaviour of rare earth elements (REEs), surface and groundwaters were collected under dry and wet conditions in different hydrological units of the Teviot Brook catchment (Southeast Queensland, Australia). Sampled waters showed a large degree of variability in both REE abundance and normalised patterns. Overall REE abundance ranged over nearly three orders of magnitude, and was consistently lower in the sedimentary bedrock aquifer (18ppt<∑REE<477ppt) than in the other hydrological systems studied. Abundance was greater in springs draining rhyolitic rocks (∑REE=300 and 2054ppt) than in springs draining basalt ranges (∑REE=25 and 83ppt), yet was highly variable in the shallow alluvial groundwater (16ppt<∑REE<5294ppt) and, to a lesser extent, in streamwater (85ppt<∑REE<2198ppt). Generally, waters that interacted with different rock types had different REE patterns. In order to obtain an unbiased characterisation of REE patterns, the ratios between light and middle REEs (R(M/L)) and the ratios between middle and heavy REEs (R(H/M)) were calculated for each sample. The sedimentary bedrock aquifer waters had highly evolved patterns depleted in light REEs and enriched in middle and heavy REEs (0.17aquifer had very diverse patterns with important spatial variations. Samples collected from the alluvium exhibited an increasing influence of the sedimentary bedrock from upgradient to downgradient; typically they showed flat patterns in the upstream section of the alluvium (median R(M/L)=0.21 and median R(H/M)=-0.06) gradually evolving towards patterns depleted in light REEs and enriched in middle and heavy REEs downgradient (median R(M/L)=0.48 and median R(H/M)=0.38). To document the seasonal variations in REE patterns, the difference in ratios between dry and wet sampling campaigns was determined for each repeated sampling location. Contributions from the sedimentary bedrock water to the alluvium during the wet season were identified at two locations (increase from R(H/M)=0.03 and 0.35 to R(H/M)=0.62 and 0.89). The effect of recharge through fractured igneous rocks was also observed in two boreholes intercepting the sedimentary bedrock, where the freshly recharged waters likely contributed to the deeper groundwater flow during the wet season (decrease from R(M/L)=0.81 and 0.56 to R(M/L)=0.46 and 0.17). Results from this study suggest that REEs may be usefully applied as indicators of recharge processes and inter-aquifer mixing. They also underline the importance of conducting seasonal sampling campaigns to capture possible short-term variations in REE patterns and abundance, which is essential to enable a better understanding of hydrological and hydrochemical processes in complex geological settings