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.

Timing and duration of the last interglacial inferred from high resolution U-series chronology of stalagmite growth in Southern Hemisphere

High-precision Th-230-U-238 ages for a stalagmite from Newdegate Cave in southern Tasmania, Australia define a rare record of precipitation between 100 and 155 ka before the present. The fastest stalagmite growth occurred between 129.2 +/- 1.6 and 122.1 +/- 2.0 ka (similar to 61.5 mm/ka), coinciding with a time of prolific coral growth from Western Australia (128-122 ka). This is the first high-resolution continental record in the Southern Hemisphere that can be compared and correlated with the marine record. Such correlation shows that in southern Australia the onset of full interglacial sea level and the initiation of highest precipitation on land were synchronous. The stalagmite growth rate between 129.2 and 142.2 ka (similar to 5.9 mm/ka) was lower than that between 142.2 and 154.5 ka (similar to 18.7 mm/ka), implying drier conditions during the Penultimate Deglaciation, despite rising temperature and sea level. This asymmetrical precipitation pattern is caused by latitudinal movement of subtropical highs and an associated Westerly circulation, in response to a changing Equator-to-Pole temperature gradient. Both marine and continental records in Australia strongly suggest that the insolation maximum between 126 and 128 ka at 65 degreesN was directly responsible for the maintenance of full Last Interglacial conditions, although the triggers that initiated Penultimate Deglaciation (at similar to 142 ka) remain unsolved. (C) 2001 Elsevier Science B.V. All rights reserved.

The geochemistry of late Archaean microbial carbonate: Implications for ocean chemistry and continental erosion history

Trace element concentrations and combined Sr- and Nd-isotope compositions were determined on stromatolitic carbonates (microbialites) from the 2.52 Ga Campbellrand carbonate platform (South Africa). Shale-normalised rare earth element and yttrium patterns of the ancient samples are similar to those of modern seawater in having positive La and Y anomalies and in being depleted in light rare earth elements. In contrast to modem seawater (and microbialite proxies), the 2.52 Ga samples lack a negative Ce anomaly but possess a positive Eu anomaly. These latter trace element characteristics are interpreted to reflect anoxic deep ocean waters where, unlike today, hydrothermal Fe input was not oxidised, and scavenged and rare earth elements were not coprecipitated with Fe-oxyhydroxides. The persistence of a positive Eu anomaly in relatively shallow Campbellrand platform waters indicates a dramatic reversal from hydrothermally dominated (Archaean) to continental erosion-dominated (Phanerozoic) rare earth element flux ratio. The dominant hydrothermal input is also expressed in the initial Sr- and Nd-isotope ratios. There is collinear variation in Sr-Nd systematics, which range from primitive values (Sr-87/Sr-86 of 0.702386 and epsilon (Nd) of +2.1) to more evolved crustal ratios. Mixing calculations show that the range in trace element ratios (e.g., Y/Ho) and initial isotope ratios is not a result of contamination by trapped sediment, but that the chemical band isotopic variation reflects carbonate deposition in an environment where different water masses mixed. Calculated Nd flux ratios yield a hydrothermal input into the 2.52 Ga oceans one order of magnitude larger than continental input. Such a change in flux ratio most likely required substantially reduced continental inputs, which could, in turn, reflect a plate tectonic causation (e.g., reduced topography or expansion of epicontinental seas). Copyright (C) 2001 Elsevier Science Ltd.

Petrology and geochemistry of the Bandas del Sur Formation, Las Canadas Edifice, Tenerife (Canary Islands)

The Bandas del Sur Formation preserves a Quaternary extra-caldera record of central phonolitic explosive volcanism of the Las Canadas volcano at Tenerife. Volcanic rocks are bimodal in composition, being predominantly phonolitic pyroclastic deposits, several eruptions of which resulted in summit caldera collapse, alkali basaltic lavas erupted from many fissures around the flanks. For the pyroclastic deposits, there is a broad range of pumice glass compositions from phonotephrite to phonolite. The phonolite pyroclastic deposits are also characterized by a diverse, 7-8-phase phenocryst assemblage (alkali feldspar + biotite + sodian diopside + titanomagnetite + ilmenite + nosean-hauyne + titanite + apatite) with alkali feldspar dominant, in contrast to interbedded phonolite lavas that typically have lower phenocryst contents and lack hydrous phases. Petrological and geochemical data are consistent with fractional crystallization (involving the observed phenocryst assemblages) as the dominant process in the development of phonolite magmas. New stratigraphically constrained data indicate that petrological and geochemical differences exist between pyroclastic deposits of the last two explosive cycles of phonolitic volcanism. Cycle 2 (0.85-0.57 Ma) pyroclastic fall deposits commonly show a cryptic compositional zonation indicating that several eruptions tapped chemically, and probably thermally stratified magma systems. Evidence for magma mixing is most widespread in the pyroclastic deposits of Cycle 3 (0.37-0.17 Ma), which includes the presence of reversely and normally zoned phenocrysts, quenched mafic glass blebs in pumice, banded pumice, and bimodal to polymodal phenocryst compositional populations. Syn-eruptive mixing events involved mostly phonolite and tephriphonolite magmas, whereas a pre-eruptive mixing event involving basaltic magma is recorded in several banded pumice-bearing ignimbrites of Cycle 3. The periodic addition and mixing of basaltic magma ultimately may have triggered several eruptions. Recharge and underplating by basaltic magma is interpreted to have elevated sulphur contents (occurring as an exsolved gas phase) in the capping phonolitic magma reservoir. This promoted nosean-hauyne crystallization over nepheline, elevated SO3 contents in apatite, and possibly resulted in large, climatologically important SO2 emissions.

Cenozoic continental weathering and its implications for the palaeoclimate: evidence from Ar-40/Ar-39 geochronology of supergene K-Mn oxides in Mt Tabor, central Queensland, Australia

One hundred and twenty-five mineral grains from 45 visually pure K-bearing Mn oxide (hollandite group) samples collected from weathering profiles in the Mt Tabor region of central Queensland, Australia, were analysed by the Ar-40/Ar-39 laser probe technique. These K-Mn oxides precipitated mainly through a process of cavity filling (direct precipitation from weathering solution), with botryoidal texture formed by micrometric mineral bands. Well-defined and reproducible plateau ages have been obtained for most samples, ranging from 27.2 +/- 0.8 to 6.8 +/- 0.5 Ma (2 sigma). Statistical analysis of the geochronological results by mixture modelling suggests an episodic mineral precipitation history, with two major peaks at 20.2 +/- 0.22 Ma and 16.5 +/- 0.17 Ma. The geochronological results, when combined with information on paragenetic relationships and mineralogical textures obtained from petrographic, scanning electron microscopy, and electron microprobe investigations, indicate that warm and humid palaeoclimatic conditions favourable to intense chemical weathering prevailed in central Queensland from late Oligocene to middle Miocene, particularly in the early Miocene. These results, in conjunction with previous and ongoing investigations in NW and eastern Queensland, suggest that most of Queensland was dominated by humid climates during the Miocene. (C) 2002 Elsevier Science BN. All rights reserved.

New W-isotope evidence for rapid terrestrial accretion and very early core formation

The short-lived Hf-182-W-182-isotope system is an ideal clock to trace core formation and accretion processes of planets. Planetary accretion and metal/silicate fractionation chronologies are calculated relative to the chondritic Hf-182-W-182-isotope evolution. Here, we report new high-precision W-isotope data for the carbonaceous chondrite Allende that are much less radiogenic than previously reported and are in good agreement with published internal Hf-W chronometry of enstatite chondrites. If the W-isotope composition of terrestrial rocks, representing the bulk silicate Earth, is homogeneous and 2.24 epsilon(182W) units more radiogenic than that of the bulk Earth, metal/silicate differentiation of the Earth occurred very early. The new W-isotope data constrain the mean time of terrestrial core formation to 34 million years after the start of solar system accretion. Early terrestrial core formation implies rapid terrestrial accretion, thus permitting formation of the Moon by giant impact while Hf-182 was still alive. This could explain why lunar W-isotopes are more radiogenic than the terrestrial value. Copyright (C) 2002 Elsevier Science Ltd.

On the overabundance of light rare earth elements in terrestrial zircons and its implication for Earth's earliest magmatic differentiation

We present whole-rock and zircon rare earth element (REE) data from two early Archaean gneisses (3.81 Ga and 3.64 Ga) from the Itsaq gneiss complex, south-west Greenland. Both gneisses represent extremely rare examples of unaltered, fresh and relatively undeformed igneous rocks of such antiquity. Cathodoluminescence imaging of their zircons indicates a single crystallisation episode with no evidence for either later metamorphic and/or anatectic reworking or inheritance of earlier grains. Uniform, single-population U/Pb age data confirm the structural simplicity of these zircons. One sample, a 3.64 Ga granodioritic gneiss from the Gothabsfjord, yields a chondrite-normalised REE pattern with a positive slope from La to Lu as well as substantial positive Ce and slight negative Eu anomalies, features generally considered to be typical of igneous zircon. In contrast, the second sample, a 3.81 Ga tonalite from south of the Isua Greenstone Belt, has variable but generally much higher light REE abundances, with similar middle to heavy REE. Calculation of zircon/melt distribution coefficients (D-REE(zircon/melt)) from each sample yields markedly different values for the trivalent REE (i.e. Ce and Eu omitted) and simple application of one set of D-REE(zircon/melt) to model the melt composition for the other sample yields concentrations that are in error by up to two orders of magnitude for the light REE (La-Nd). The observed light REE overabundance in the 3.81 Ga tonalite is a commonly observed feature in terrestrial zircons for which a number of explanations ranging from lattice strain to disequilibrium crystallisation have been proposed and are further investigated herein. Regardless of the cause of light REE overabundance, our study shows that simple application of zircon/melt distribution coefficients is not an unambiguous method for ascertaining original melt composition. In this context, recent studies that use REE data to claim that > 4.3 Ga Hadean detrital zircons originally crystallised from an evolved magma, in turn suggesting the operation of geological processes in the early Earth analogous to those of the present day (e.g. subduction and melting of hydrated oceanic crust), must be regarded with caution. Indeed, comparison of terrestrial Hadean and > 3.9 Ga lunar highland zircons shows remarkable similarities in the light REE, even though subduction processes that have been used to explain the terrestrial zircons have never operated on the Moon. (C) 2002 Elsevier Science B.V. All rights reserved.