Chalcophile element geochemistry of the Boggy Plain Zoned Pluton, southeastern Australia : a S-saturated barren compositionally diverse magma

The behavior of the platinum group elements (PGE) and Re in felsic magmas is poorly understood due to scarcity of data. We report the concentrations of Ni, Cu, Re, and PGE in the compositionally diverse Boggy Plain zoned pluton (BPZP), which shows a variation of rock type from gabbro through granodiorite and granite to aplite with a SiO2 range from 52 to 74 wt %. In addition, major silicate and oxide minerals were analyzed for Ni, Cu, and Re, and a systematic sulfide study was carried out to investigate the role of silicate, oxide, and sulfide minerals on chalcophile element geochemistry of the BPZP. Mass balance calculation shows that the whole rock Cu budget hosted by silicate and oxide minerals is <13 wt % and that Cu is dominantly located in sulfide phases, whereas most of the whole rock Ni budget (>70 wt %) is held in major silicate and oxide minerals. Rhenium is dominantly hosted by magnetite and ilmenite. Ovoid-shaped sulfide blebs occur at the boundary between pyroxene phenocrysts and neighboring interstitial phases or within interstitial minerals in the gabbro and the granodiorite. The blebs are composed of pyrrhotite, pyrite, chalcopyrite, and S-bearing Fe-oxide, which contain total trace metals (Co, Ni, Cu, Ag, Pb) up to ~16 wt %. The mineral assemblage, occurrence, shape, and composition of the sulfide blebs are a typical of magmatic sulfide. PGE concentrations in the BPZP vary by more than two orders of magnitude from gabbro (2.7–7.8 ppb Pd, 0.025–0.116 ppb Ir) to aplite (0.05 ppb Pd, 0.001 ppb Ir). Nickel, Cu, Re, and PGE concentrations are positively correlated with MgO in all the rock types although there is a clear discontinuity between the granodiorite and the granite in the trends for Ni, Rh, and Ir when plotted against MgO. Cu/Pd values gradually increase from 6,100 to 52,600 as the MgO content decreases. The sulfide petrology and chalcophile element geochemistry of the BPZP show that sulfide saturation occurred in the late gabbroic stage of magma differentiation. Segregation and distribution of these sulfide blebs controlled Cu and PGE variations within the BPZP rocks although the magma of each rock type may have experienced a different magma evolution history in terms of crustal assimilation and crystal fractionation. The sulfide melt locked in the cumulate rocks must have sequestered a significant portion of the chalcophile elements, which restricted the availability of these metals to magmatic-hydrothermal ore fluids. Therefore, we suggest that the roof rocks that overlay the BPZP were not prospective for magmatic-hydrothermal Cu, Au, or Cu–Au deposits.

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.

Assessment of flood-related recharge to alluvial aquifers of an irrigated catchment following extended drought conditions using 3D visualisation and environmental tracers, Lockyer Valley, southeast Queensland, Australia.

The Lockyer Valley in southeast Queensland, Australia, hosts an economically significant alluvial aquifer system which has been impacted by prolonged drought conditions (~1997 to ~ 2009). Throughout this time, the system was under continued groundwater extraction, resulting in severe aquifer depletion. By 2008, much of the aquifer was at <30% of storage but some relief occurred with rains in early 2009. However, between December 2010 and January 2011, most of southeast Queensland experienced unprecedented flooding, which generated significant aquifer recharge. In order to understand the spatial and temporal controls of groundwater recharge in the alluvium, a detailed 3D lithological property model of gravels, sands and clays was developed using GOCAD software. The spatial distribution of recharge throughout the catchment was assessed using hydrograph data from about 400 groundwater observation wells screened at the base of the alluvium. Water levels from these bores were integrated into a catchment-wide 3D geological model using the 3D geological modelling software GOCAD; the model highlights the complexity of recharge mechanisms. To support this analysis, groundwater tracers (e.g. major and minor ions, stable isotopes, 3H and 14C) were used as independent verification. The use of these complementary methods has allowed the identification of zones where alluvial recharge primarily occurs from stream water during episodic flood events. However, the study also demonstrates that in some sections of the alluvium, rainfall recharge and discharge from the underlying basement into the alluvium are the primary recharge mechanisms of the alluvium. This is indicated by the absence of any response to the flood, as well as the observed old radiocarbon ages and distinct basement water chemistry signatures at these locations. Within the 3D geological model, integration of water chemistry and time-series displays of water level surfaces before and after the flood suggests that the spatial variations of the flood response in the alluvium are primarily controlled by the valley morphology and lithological variations within the alluvium. The integration of time-series of groundwater level surfaces in the 3D geological model also enables the quantification of the volumetric change of groundwater stored in the unconfined sections of this alluvial aquifer during drought and following flood events. The 3D representation and analysis of hydraulic and recharge information has considerable advantages over the traditional 2D approach. For example, while many studies focus on singular aspects of catchment dynamics and groundwater-surface water interactions, the 3D approach is capable of integrating multiple types of information (topography, geological, hydraulic, water chemistry and spatial) into a single representation which provides valuable insights into the major factors controlling aquifer processes.

Understanding hydrochemical processes using 3D models and hydrochemical facies in Wairau Plain groundwater, Marlborough

Three dimensional models and groundwater quality are combined to better understand and conceptualise groundwater systems in complex geological settings in the Wairau Plain, Marlborough. Hydrochemical facies, which are characteristic of distinct evolutionary pathways and a common hydrologic history of groundwaters, are identified within geological formations to assess natural water-rock interactions, redox potential and human agricultural impact on groundwater quality in the Wairau Plain.

Sources of hydrochemical variability and implications associated with CSG water extraction

Coal Seam Gas (CSG) production is achieved by extracting groundwater to depressurize coal seam aquifers in order to promote methane gas desorption from coal micropores. CSG waters are characteristically alkaline, have a neutral pH (~7), are of the Na-HCO3-Cl type, and exhibit brackish salinity. In 2004, a CSG exploration company carried out a gas flow test in an exploration well located in Maramarua (Waikato Region, New Zealand). This resulted in 33 water samples exhibiting noteworthy chemical variations induced by pumping. This research identifies the main causes of hydrochemical variations in CSG water, makes recommendations to manage this effect, and discusses potential environmental implications. Hydrochemical variations were studied using Factor Analysis and this was supported with hydrochemical modelling and a laboratory experiment. This reveals carbon dioxide (CO2) degassing as the principal source of hydrochemical variability (about 33%). Factor Analysis also shows that major ion variations could also reflect changes in hydrochemical composition induced by different pumping regimes. Subsequent chloride, calcium, and TDS variations could be a consequence of analytical errors potentially committed during laboratory determinations. CSG water chemical variations due to degassing during pumping can be minimized with good completion and production techniques; variations due to sample degassing can be controlled by taking precautions during sampling, transit, storage and analysis. In addition, the degassing effect observed in CSG waters can lead to an underestimation of their potential environmental effect. Calcium precipitation due to exposure to normal atmospheric pressure results in a 23% increase in SAR values from Maramarua CSG water samples.

GWCarb v1.0: carbonate speciation tool

This spreadsheet calculates carbonate speciation using carbonate equilibrium equations at standard conditions (T=25°C) with ionic strength corrections. The user will typically be able to calculate the different carbonate species by entering total alkalinity and pH. This spreadsheet contains additional tools to calculate the Langelier Index for calcium and the SAR of the water. Note that in this last calculation the potential for calcium precipitation is not taken into account. The last tool presented here is a carbonate speciation tool in open systems (e.g. open to the atmosphere) which takes into account atmospheric pressure.