Small back-barrier sand islands are poorly known in terms of hydrogeology and have been overlooked in more extensive studies of coastal groundwater systems that include larger barrier island complexes. This study employs a three-fold sequential approach to aquifer characterisation in a back-barrier sand island. A three-dimensional stratigraphic model forms the foundation framework, being derived from a multidisciplinary approach to sedimentary analysis and the construction of a depositional chronology. A conceptual hydrostratigraphic model is formulated based on the translation of sedimentary facies to hydrofacies, combined with density dependent flow calculations and tidal oscillation measurements. Groundwater hydrochemical data and mineral geochemistry are integrated with the resulting hydrogeological model to examine water-rock interaction and solute transport mechanisms. The study area is Toorbul Island, a small back-barrier sand mass of ~5 km2 with a maximum surface elevation of ~3.5 m AHD, located in the Pumicestone Passage of Southeast Queensland. The island hosts a dual aquifer system consisting of an unconfined island freshwater lens, underlain by a semi-confined palaeovalley-fill aquifer. Groundwater in the semi-confined aquifer is hyper-saline, carrying high concentrations of dissolved metals, with iron, in particular, ranging from 40 to < 200 mg l-1. This is of significant interest for both human health and environmental management, because iron is an important nutrient source for toxic algal bacteria such as Lyngbya majuscula. Conceptual modelling demonstrates that iron oxides and hydroxides are the main source of iron in the semi-confined aquifer, with a contribution from ferruginous chlorite dissolution. Aqueous manganese and a proportion of the aqueous iron are derived from the dissolution of manganoan ilmenite. Ferric iron minerals also contribute a significant proportion of dissolved iron in the deeper regions of the unconfined aquifer. Aqueous iron in the shallow unconfined groundwater is limited by iron sulphides, which also regulate acidity and indirectly limit dissolved aluminium concentrations. Groundwater redox state governed by seasonal climatic fluxes is the most significant control on iron-bearing mineral phase stability. Transport of dissolved metals to the surrounding estuary and the adjacent barrier island groundwater system is limited by the rate of ion diffusion across transition zone boundaries. The overall conclusions derived from this research show that back-barrier islands should be evaluated as discrete hydrogeological entities. The stratigraphic complexity that may be apparent within these island landforms should not be underestimated and the model domain should not necessarily be treated as a homogeneous system. This complexity is exemplified by the relationship between the upper and lower aquifers on Toorbul Island and the associated distribution of groundwater compositional heterogeneity. The complex stratigraphy within the sedimentary pile is derived from the presence of a sub-surface palaeovalley and the sedimentary response to changing sea-level over time. Considering the current widespread distribution of estuarine systems, complex hydrogeology as exhibited by Toobul Island, may be common in many small back-barrier island groundwater systems. The aquifer characteristics and their influence on solute transport and delivery can have significant ramifications for the exploitation of the adjacent coastal plain and barrier island aquifers. The potential influence on the latter is of particular concern due to the pressure imposed on potable groundwater supplies by increasing population densities in coastal areas.

Accessibility to consistent subsurface hydrostratigraphic information is crucial for the development of robust groundwater flow and contaminant transport models. However, full three-dimensional understanding of the subsurface geology is often the missing link. Construction of watershed-scale hydrostratigraphic models continues to be limited by the quality and density of borehole data which often lack detailed geologic information. This can become a serious problem where rapid sediment facies changes and intricate sediment architecture occur. This research is motivated by the idea that if we can understand more about the distribution of sediments and structures of complex deposits, we learn more about depositional processes and how they affect the internal geometry of a deposit and the distribution of hydraulic properties. One approach is to study surficial excavations (e.g. sand and gravel pits) that often punctuate shallow aquifers. The purpose of this study is to develop and test a method of integrating high-resolution georeferenced stratigraphic and sedimentologic information from sand and gravel pits as a means to better document sedimentologic data and improve understanding of the depositional environments. The study area is located within the Waterloo Moraine, in southwestern Ontario, and is an unconsolidated shallow aquifer system with a complex internal architecture and sediment heterogeneity. The method involves the integration of high-resolution field data with borehole and geophysical information in a computer-based 3D environment. A total of fourteen virtual sedimentary sections were constructed by georegistering digital photographs within a framework of georeferenced positions collected using a reflectorless total station and GPS. Fourteen sediment facies have been described in the field. These include crudely stratified gravel beds, planar and cross-laminated sandy strata (ripple and dune scales), along with laminated and massive silty and clayey beds. Calculated hydraulic conductivities span over seven orders of magnitude. The analysis of a single excavation has shown contrasting sediment assemblages from one end of the pit to the other, highlighting the complexity of the Waterloo Moraine. The heterogeneous and deformed layers of gravel, sand, and mud may be the product of an ice-contact to ice-proximal environment, whereas the extensive sandy assemblages may reflect an intermediate subaqueous fan region. The results also suggest that the borehole database overestimates the amount of fine-grained material in the study area. Finally, this research demonstrates that it is possible to build in a timely manner a 3D virtual sedimentologic database. New emerging technologies will lead to increased resolution and accuracy, and will help streamline the process even further. The possibility of expanding the 3D geodatabase to other excavations across the region in a timely manner is likely to lead to improved hydrostratigraphic models and, by extension, to more efficient strategies in water resources planning, management and protection.

Accessibility to consistent subsurface hydrostratigraphic information is crucial for the development of robust groundwater flow and contaminant transport models. However, full three-dimensional understanding of the subsurface geology is often the missing link. Construction of watershed-scale hydrostratigraphic models continues to be limited by the quality and density of borehole data which often lack detailed geologic information. This can become a serious problem where rapid sediment facies changes and intricate sediment architecture occur. This research is motivated by the idea that if we can understand more about the distribution of sediments and structures of complex deposits, we learn more about depositional processes and how they affect the internal geometry of a deposit and the distribution of hydraulic properties. One approach is to study surficial excavations (e.g. sand and gravel pits) that often punctuate shallow aquifers. The purpose of this study is to develop and test a method of integrating high-resolution georeferenced stratigraphic and sedimentologic information from sand and gravel pits as a means to better document sedimentologic data and improve understanding of the depositional environments. The study area is located within the Waterloo Moraine, in southwestern Ontario, and is an unconsolidated shallow aquifer system with a complex internal architecture and sediment heterogeneity. The method involves the integration of high-resolution field data with borehole and geophysical information in a computer-based 3D environment. A total of fourteen virtual sedimentary sections were constructed by georegistering digital photographs within a framework of georeferenced positions collected using a reflectorless total station and GPS. Fourteen sediment facies have been described in the field. These include crudely stratified gravel beds, planar and cross-laminated sandy strata (ripple and dune scales), along with laminated and massive silty and clayey beds. Calculated hydraulic conductivities span over seven orders of magnitude. The analysis of a single excavation has shown contrasting sediment assemblages from one end of the pit to the other, highlighting the complexity of the Waterloo Moraine. The heterogeneous and deformed layers of gravel, sand, and mud may be the product of an ice-contact to ice-proximal environment, whereas the extensive sandy assemblages may reflect an intermediate subaqueous fan region. The results also suggest that the borehole database overestimates the amount of fine-grained material in the study area. Finally, this research demonstrates that it is possible to build in a timely manner a 3D virtual sedimentologic database. New emerging technologies will lead to increased resolution and accuracy, and will help streamline the process even further. The possibility of expanding the 3D geodatabase to other excavations across the region in a timely manner is likely to lead to improved hydrostratigraphic models and, by extension, to more efficient strategies in water resources planning, management and protection.