Geochemical mass-balance and oxygen-isotope constrains on silcrete formation and its paleoclimatic implications in Southern Australia

Two geographically distinct silcrete associations are present in southern Australia, inland and eastern; these were sampled in central South Australia and central Victoria, respectively, At each site, both silicified and immediately adjacent unsilicified parent material were collected. Analytical data from these pairs were used to construct isocons, assuming Zr immobility, and to calculate the volume change and amount of silica introduced during silicification, These results, together with whole-rock oxygen isotope compositions, were used to determine the delta(18)O of th, introduced silica, The results show that the eastern silcretes in central Victoria are probably linked genetically to the associated basalts, weathering of which supplied the introduced silica, This conclusion is based on the close spatial connection between the two, as well as the substantial amount of introduced silica in the silcretes (greater than in the inland silcretes), resulting in volume increases in some eastern silcretes, Oxygen isotopic calculations for the silcretes indicate that the silica precipitated from groundwaters at temperatures slightly higher than present conditions. Silcrete formation apparently occurred during the Miocene and Pliocene (basalts in Victoria younger than Pliocene lack associated silcrete) and may reflect the much wetter climate in southeastern Australia at that time. The inland silcretes of central South Australia can be divided into pedogenic (the most common) and groundwater varieties. The pedogenic silcretes, which show typical soil features like columnar and nodular textures, contain moderate amounts of introduced silica that precipitated by evaporation from saline groundwaters, For the groundwater silcretes, which have massive textures and formed at or close to the water table, insufficient data are available to determine the mode of formation. The inland pedogenic silcretes have probably been farming from the Eocene-Miocene to the present, implying that conditions of seasonally high evaporation have occurred in central Australia during this time period. Thus silcrete formation depends on a complex interplay between climate and silica supply, and it is impossible to generalize that the presence of silcrete is indicative of a particular climate. Likewise, the elemental composition of silcretes, particularly Ti content, is not necessarily of climatic significance, Nevertheless, detailed geochemical and oxygen isotopic studies of a silcrete and its parent material can elucidate the mechanisms of silcrete formation, and if evaporation is indicated as a major factor in silcrete formation, then the climate at the time was likely to have been at least seasonally arid.

Computer simulations of coupled problems in geological and geochemical systems

The finite element method is used to simulate coupled problems, which describe the related physical and chemical processes of ore body formation and mineralization, in geological and geochemical systems. The main purpose of this paper is to illustrate some simulation results for different types of modelling problems in pore-fluid saturated rock masses. The aims of the simulation results presented in this paper are: (1) getting a better understanding of the processes and mechanisms of ore body formation and mineralization in the upper crust of the Earth; (2) demonstrating the usefulness and applicability of the finite element method in dealing with a wide range of coupled problems in geological and geochemical systems; (3) qualitatively establishing a set of showcase problems, against which any numerical method and computer package can be reasonably validated. (C) 2002 Published by Elsevier Science B.V.

Clay mineralogical, geochemical and isotopic tracing of the evolution of the Woodleigh impact structure, Southern Carnarvon Basin, Western Australia

Chaotically structured diamictite from the inner ring syncline surrounding the central uplift of the Woodleigh impact structure contains shocked metamorphic and impact melt-rock fragments, largely derived from Ordovician and Devonian target sandstones. Coarse illite fractions (< 2 mu m) from the sandstones containing no K-feldspar yield K-Ar ages of around 400 Ma, whereas the K-Ar ages of authigenic clays of > 0.2 mu m fractions from the diamictite without smectite and K-feldspar cluster around 360 Ma, consistent with Rb-Sr data. Crystallisation of newly formed illite in the impact melt rock clasts and recrystallisation of earlier formed illite in the sandstone clasts preserved in the diamictite, are attributed to impact-induced hydrothermal processes in the Late Devonian. The illitic clays from the diamictite and from the sandstones have very similar trace element compositions, with significantly enriched incompatible lithophile elements, which increase in concentrations correlatively with those of the compatible ferromagnesian elements. The unusual trace element associations in the clays may be due to the involvement of hot gravity-driven basinal fluids that interacted with rocks of the Precambrian craton to the east of the study area, or with such material transported and reworked in the studied sedimentary succession.

Geochemical changes recorded in Lynch's Crater, Northeastern Australia, over the past 50 ka

There are many geochemical reconstructions of environmental change in the mid and high latitudes but relatively few in the tropical latitudes, despite their considerable potential for reconstructing environmental processes that cannot be identified using more traditional proxies. Here we present one reconstruction of environmental change for the tropics. This reconstruction covers the past 50 ka using a suite of geochemical data from the high-resolution sequence of Lynch's Crater in northeast Queensland, Australia, a region highly sensitive to El Nino-Southern Oscillation (ENSO) activity. The 23 major oxides and trace elements measured Could be summarised by extracting three axes using principal components analysis (accounting for 72% of the variability). The data indicate that the greatest variability in the geochemical data accounted for erosional activity within the catchment that was associated with past changes in the frequency of ENSO activity (though this was less sensitive during wetter periods, probably as a result of buffering by high vegetation cover). The remaining variability was largely explained by elements that form complexes with organic compounds (e.g., humic acids) and those that are important nutrients for specific vegetation types (and therefore a measure of vegetation distribution). For more detailed reconstructions, further work is required to disentangle the complex controls of clements within sedimentary sequences. (c) 2005 Elsevier B.V. All rights reserved.

Finite element modelling of rock alteration and metamorphic process in hydrothermal systems

We use the finite element method to simulate the rock alteration and metamorphic process in hydrothermal systems. In particular, we consider the fluid-rock interaction problems in pore-fluid saturated porous rocks. Since the fluid rock interaction takes place at the contact interface between the pore-fluid and solid minerals, it is governed by the chemical reaction which usually takes place very slowly at this contact interface, from the geochemical point of view. Due to the relative slowness of the rate of the chemical reaction to the velocity of the pore-fluid flow in the hydrothermal system to be considered, there exists a retardation zone, in which the conventional static theory in geochemistry does not hold true. Since this issue is often overlooked by some purely numerical modellers, it is emphasized in this paper. The related results from a typical rock alteration and metamorphic problem in a hydrothermal system have shown not only the detailed rock alteration and metamorphic process, but also the size of the retardation zone in the hydrothermal system. Copyright (C) 2001 John Wiley & Sons, Ltd.

Numerical simulation of double-diffusion driven convective flow and rock alteration in three-dimensional fluid-saturated geological fault zones

Numerical methods are used to simulate the double-diffusion driven convective pore-fluid flow and rock alteration in three-dimensional fluid-saturated geological fault zones. The double diffusion is caused by a combination of both the positive upward temperature gradient and the positive downward salinity concentration gradient within a three-dimensional fluid-saturated geological fault zone, which is assumed to be more permeable than its surrounding rocks. In order to ensure the physical meaningfulness of the obtained numerical solutions, the numerical method used in this study is validated by a benchmark problem, for which the analytical solution to the critical Rayleigh number of the system is available. The theoretical value of the critical Rayleigh number of a three-dimensional fluid-saturated geological fault zone system can be used to judge whether or not the double-diffusion driven convective pore-fluid flow can take place within the system. After the possibility of triggering the double-diffusion driven convective pore-fluid flow is theoretically validated for the numerical model of a three-dimensional fluid-saturated geological fault zone system, the corresponding numerical solutions for the convective flow and temperature are directly coupled with a geochemical system. Through the numerical simulation of the coupled system between the convective fluid flow, heat transfer, mass transport and chemical reactions, we have investigated the effect of the double-diffusion driven convective pore-fluid flow on the rock alteration, which is the direct consequence of mineral redistribution due to its dissolution, transportation and precipitation, within the three-dimensional fluid-saturated geological fault zone system. (c) 2005 Elsevier B.V. All rights reserved.

Three-dimensional model for multi-component reactive transport with variable density groundwater flow

PHWAT is a new model that couples a geochemical reaction model (PHREEQC-2) with a density-dependent groundwater flow and solute transport model (SEAWAT) using the split-operator approach. PHWAT was developed to simulate multi-component reactive transport in variable density groundwater flow. Fluid density in PHWAT depends not on only the concentration of a single species as in SEAWAT, but also the concentrations of other dissolved chemicals that can be subject to reactive processes. Simulation results of PHWAT and PHREEQC-2 were compared in their predictions of effluent concentration from a column experiment. Both models produced identical results, showing that PHWAT has correctly coupled the sub-packages. PHWAT was then applied to the simulation of a tank experiment in which seawater intrusion was accompanied by cation exchange. The density dependence of the intrusion and the snow-plough effect in the breakthrough curves were reflected in the model simulations, which were in good agreement with the measured breakthrough data. Comparison simulations that, in turn, excluded density effects and reactions allowed us to quantify the marked effect of ignoring these processes. Next, we explored numerical issues involved in the practical application of PHWAT using the example of a dense plume flowing into a tank containing fresh water. It was shown that PHWAT could model physically unstable flow and that numerical instabilities were suppressed. Physical instability developed in the model in accordance with the increase of the modified Rayleigh number for density-dependent flow, in agreement with previous research. (c) 2004 Elsevier Ltd. All rights reserved.

Yarrol terrane of the northern New England Fold Belt: forearc or backarc?

The Upper Devonian to Lower Carboniferous volcanosedimentary rocks of the Yarrol terrane of the northern New England Fold Belt have previously been ascribed to a forearc basin setting. New data presented here, however, suggest that the Yarrol terrane developed as a backarc basin during the Middle to early Late Devonian. Based on field studies, we recognise four regionally applicable strati graphic units: (i) a basal, ?Middle to Upper Devonian submarine mafic volcanic suite (Monal volcanic facies association); (ii) the lower Frasnian Lochenbar beds that locally unconformably overlie the Monal volcanic facies association: (iii) the Three Moon Conglomerate (Upper Devonian - Lower Carboniferous): and (iv) the Lower Carboniferous Rockhampton Group characterised by the presence of oolitic limestone. Stratigraphic and compositional differences suggest the Monal volcanic facies association post-dates Middle Devonian silicic-dominated magmatism that was coeval with gold-copper mineralisation at Mt Morgan. The Lochenbar beds, Three Moon Conglomerate and Rockhampton Group represent a near-continuous sedimentary record of volcanism that changed in composition and style from mafic effusive (Late Devonian) to silicic explosive volcanism (Early Carboniferous). Palaeocurrent data from the Three Moon Conglomerate and Rockhampton Group indicate dispersal of sediment to the west and northwest, and are inconsistent with derivation from a volcanic-are source situated to the west (Connors-Auburn Arch). Geochemical data show that the Monal volcanic facies association ranges from tholeiitic subalkaline basalts to calc-alkaline basaltic andesite. Trace and rare-earth element abundances are distinctly MORE-like (e.g, light rare earth element depletion), with only moderate enrichment of the large-ion lithophile elements in some units, and negative Nb anomalies, suggesting a subduction-related signature. Basalts of the Monal volcanic facies association are best described as transitional between calc-alkali basalts and N-MORB. The elevated high field strength element contents (e.g. Zr, Y, Ti) are higher than modern island-are basalts, but comparable to basalts that floor modern backarc basins. This geochemical study, coupled with stratigraphic relationships, suggest that the eruption of backarc basin basalts followed widespread Middle Devonian, extension-related silicic magmatism (e.g. Retreat Batholith, Mt Morgan), and floored the Yarrol terrane. The Monal volcanic facies association thus shows similarities in its tectonic environment to the Lower Permian successions (e.g. Rookwood Volcanics) of the northern New England Fold Belt. These mafic volcanic sequences are interpreted to record two backarc basin-forming periods (Middle - Late Devonian and Late Carboniferous - Early Permian) during the Late Palaeozoic history of the New England Orogen. Silicic-dominated explosive volcanism, occurring extensively across the northern New England Fold Belt in the Early Carboniferous (Varrol terrane, Campwyn Volcanics, Drummond and Burdekin Basins), reflects another period of crustal melting and extension, most likely related to the opening of the Drummond Basin.