Finite element modeling of fluid-rock interaction problems in pore-fluid saturated hydrothermal/sedimentary basins

In order to use the finite element method for solving fluid-rock interaction problems in pore-fluid saturated hydrothermal/sedimentary basins effectively and efficiently, we have presented, in this paper, the new concept and numerical algorithms to deal with the fundamental issues associated with the fluid-rock interaction problems. These fundamental issues are often overlooked by some purely numerical modelers. (1) Since the fluid-rock interaction problem involves heterogeneous chemical reactions between reactive aqueous chemical species in the pore-fluid and solid minerals in the rock masses, it is necessary to develop the new concept of the generalized concentration of a solid mineral, so that two types of reactive mass transport equations, namely, the conventional mass transport equation for the aqueous chemical species in the pore-fluid and the degenerated mass transport equation for the solid minerals in the rock mass, can be solved simultaneously in computation. (2) Since the reaction area between the pore-fluid and mineral surfaces is basically a function of the generalized concentration of the solid mineral, there is a definite need to appropriately consider the dependence of the dissolution rate of a dissolving mineral on its generalized concentration in the numerical analysis. (3) Considering the direct consequence of the porosity evolution with time in the transient analysis of fluid-rock interaction problems; we have proposed the term splitting algorithm and the concept of the equivalent source/sink terms in mass transport equations so that the problem of variable mesh Peclet number and Courant number has been successfully converted into the problem of constant mesh Peclet and Courant numbers. The numerical results from an application example have demonstrated the usefulness of the proposed concepts and the robustness of the proposed numerical algorithms in dealing with fluid-rock interaction problems in pore-fluid saturated hydrothermal/sedimentary basins. (C) 2001 Elsevier Science B.V. All rights reserved.

Finite element modelling of reactive fluids mixing and mineralization in pore-fluid saturated hydrothermal/sedimentary basins

We present the finite element simulations of reactive mineral carrying fluids mixing and mineralization in pore-fluid saturated hydrothermal/sedimentary basins. In particular we explore the mixing of reactive sulfide and sulfate fluids and the relevant patterns of mineralization for Load, zinc and iron minerals in the regime of temperature-gradient-driven convective flow. Since the mineralization and ore body formation may last quite a long period of time in a hydrothermal basin, it is commonly assumed that, in the geochemistry, the solutions of minerals are in an equilibrium state or near an equilibrium state. Therefore, the mineralization rate of a particular kind of mineral can be expressed as the product of the pore-fluid velocity and the equilibrium concentration of this particular kind of mineral Using the present mineralization rate of a mineral, the potential of the modern mineralization theory is illustrated by means of finite element studies related to reactive mineral-carrying fluids mixing problems in materially homogeneous and inhomogeneous porous rock basins.

The application of electron microscopy and microanalysis in conjunction with organic petrology to further the understanding of organic-mineral association: examples from Mount Isa and McArthur Basins, Australia

Organic petrology supported by electron microscopical and micro-analytical techniques was applied to organic matter in Proterozoic sediments to better understand hydrothermal processes responsible for ore-grade mineralisation. It was shown that organic maturation was not only closely linked to the geological history of the sediments, but also highlighted heat transfer by convection as differentiated from conduction solely through sediment burial and step-wise subsidence. Water-rock ratios effect organic maturation in hydrothermal systems, and erratic reflectance profiles are indicators of convective heat transfer. Identification and characterisation of organic materials in terms of source rock and migrated hydrocarbons was shown to be a powerful tool in reconstructing the thermal history of sediments, identifying hydrothermal episodes, fluid pathways and heat source in the northern Australian Proterozoic basins. Higher reflectance of organic matter towards the central parts of the Mount Isa Basin and some of the most northerly parts point to proximity to higher heat flow at times, in contrast to relatively low temperatures prevailing in the western parts of the basin, next to the Murphy Inlier. A close correlation shown between peak organic reflectance values and super-sequence boundaries farther highlighted the valuable information to be gained from organic petrology, by allowing the separation of processes responsible for metal dissolution and transportation from those inducing precipitation. (C) 2001 Elsevier Science B.V All rights reserved.

Relations between coal petrology and gas content in the Upper Newlands Seam, central Queensland, Australia

The Upper Newlands Seam in the northern Bowen Basin, Queensland, Australia consists of six benches (A-F) that have different petrographic assemblages. Benches C and E contain relatively abundant inertodetrinite and mineral matter, as well as anomalously high reflectance values; these characteristics support a largely allochthonous, detrital origin for the C and E benches. Fractures and cleats in the seam show a consistent orientation of northeast- southwest for face cleats, and a wide range of orientations for fractures. Cleat systems are well developed in bright bands, with poor continuity in the dull coal. Both maceral content and cleat character are suggested to influence gas drainage in the Upper Newlands Seam. A pronounced positive correlation between vitrinite abundance and gas desorption data suggests more efficient drainage from benches with abundant vitrinite. Conversely, inertinite-rich benches are suggested to have less efficient drainage, and possibly retain gas within pore spaces, which could increase the outburst potential of the coal. (C) 2001 Elsevier Science BN. All rights reserved.

Sedimentary cyclicity in CRP drillcore, Victoria Land Basin, Antarctica

The upper 1200 m of pre-Pliocene sediment recovered by Cape Roberts Project (CRP) drilling off the Victoria Land coast of Antarctica between 1997-1999 has been subdivided into 54 unconformity-bound stratigraphic sequences, spanning the period c. 32 to 17 Ma. The sequences are recognised on the basis of the cyclical vertical stacking of their constituent lithofacies, which are enclosed by erosion surfaces produced during the grounding of the advancing ice margin onto the sea floor. Each sequence represents deposition in a range of offshore shelf to coastal glacimarine sedimentary environments during oscillations in the ice margin across the Western Ross Sea shelf, and coeval fluctuations in water depth. This paper applies spectral analysis techniques to depth- and time-series of sediment grain size (500 samples) for intervals of the core with adequate chronological data. Time series analysis of 0.5-1.0m-spaced grainsize data spanning sequences 9-11 (CRP-2/2A) and sequences 1-7 (CRP-3) suggests that the length of individual sequences correspond to Milankovitch frequencies, probably 41 k.y., but possibly as low as 100 k.y. Higher frequency periodic components at 23 k.y. (orbital precession) and 15-10 k.y. (sub-orbital) are recognised at the intrasequence-scale, and may represent climatic cycles akin to the ice rafting episodes described in the North Atlantic Ocean during the Quaternary. The cyclicity recorded by glacimarine sequences in CRP core provides direct evidence from the periphery of Antarctica for orbital oscillations in the size of the Oligocene-Early Miocene East Antarctic Ice Sheet.

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.

Geological and trace element evidence for a marine sedimentary environment of deposition and biogenicity of 3.45 Ga stromatolitic carbonates in the Pilbara Craton, and support for a reducing Archaean ocean

Bedded carbonate rocks from the 3.45 Ga Warrawoona Group, Pilbara Craton, contain structures that have been regarded either as the oldest known stromatolites or as abiotic hydrothermal deposits. We present new field and petrological observations and high-precision REE + Y data from the carbonates in order to test the origin of the deposits. Trace element geochemistry from a number of laminated stromatolitic dolomite samples of the c. 3.40 Ga Strelley Pool Chert conclusively shows that they precipitated from anoxic seawater, probably in a very shallow environment consistent with previous sedimentological observations. Edge-wise conglomerates in troughs between stromatolites and widespread cross-stratification provide additional evidence of stromatolite construction, at least partly, from layers of particulate sediment, rather than solely from rigid crusts. Accumulation of particulate sediment on steep stromatolite sides in a high-energy environment suggests organic binding of the surface. Relative and absolute REE + Y contents are exactly comparable with Late Archaean microbial carbonates of widely agreed biological origin. Ankerite from a unit of bedded ankerite–chert couplets from near the top of the stratigraphically older (3.49 Ga) Dresser Formation, which immediately underlies wrinkly stromatolites with small, broad, low-amplitude domes, also precipitated from anoxic seawater. The REE + Y data of carbonates from the Strelley Pool Chert and Dresser Formation contrast strongly with those from siderite layers in a jasper–siderite–Fe-chlorite banded iron-formation from the base of the Panorama Formation (3.45 Ga), which is clearly hydrothermal in origin. The geochemical results, together with sedimentological data, strongly support: (1) deposition of Dresser Formation and Strelley Pool Chert carbonates from Archaean seawater, in part as particulate carbonate sediment; (2) biogenicity of the stromatolitic carbonates; (3) a reducing Archaean atmosphere; (4) ongoing extensive terrestrial erosion prior to ∼3.45 Ga.

Importance of water for Archaean granitoid petrology: a comparative study of TTG and potassic granitoids from Barberton Mountain Land, South Africa

A new model for Archaean granitoid magmatism is presented which reconciles the most important geochemical similarities and differences between tonalite-trondhjemite-granodiorite (TTG) and potassic granitoids. Trace element abundances reveal a strong arc magmatism signature in all studied granitoids from Barberton Mountain Land. Characteristic features include HFSE depletion as well as distinct enrichment peaks of fluid-sensitive trace elements such as Pb in N-MORB normalisation, clearly indicating that all studied granitoids are derived from refertilised mantle above subduction zones. We envisage hydrous basaltic (s.l.) melts as parental liquids, which underwent extensive fractional crystallisation. Distinctive residual cumulates evolved depending on initial water content. High-H2O melts crystallised garnet/amphibole together with pyroxenes and minor plagioclase, but no olivine. This fractionation path ultimately led to TTG-like melts. Less hydrous basaltic melts also crystallised garnet/amphibole, but the lower compatible element content indicates that olivine was also a liquidus phase. Pronounced negative Eu-anomalies of the granitic melts, correlating with Na, Ca and Al, indicate plagioclase to be of major importance. In the context of our model, the post-Archaean disappearance of TTG and concomitant preponderance of granites (s.l.), therefore, is explained with secular decrease of aqueous fluid transport into subduction zones and/or efficiency of deep fluid release from slabs.

Late Quaternary uplift and subsidence of the west coast of Tanna, south Vanuatu, southwest Pacific: U-Th ages of raised coral reefs in the Median Sedimentary Basin

Twelve Late Quaternary TIMS U-Th ages are reported here from 10 coral samples collected in situ from five transgressive coral/algal raised reefs (height: max. 113 m, min. 8 m) and two raised lagoonal deposits (height: max. 18 m, min. 8 m) along and near the west coast of Tanna, which lies in the Median Sedimentary Basin of South Vanuatu, southwest Pacific. These reefs and raised lagoonal deposits represent several age groups: (i) 215 ka (marine oxygen-isotope stage 7) penultimate interglacial (highest elevation and oldest); (ii) one lagoonal deposit of ca 127 ka (marine oxygen-isotope stage 5e); (iii) three last interglacial reefs with ages 102, 89 and 81 ka (representing marine oxygen-isotope stages 5c, 5b and 5a, respectively, of the latter part of the last interglacial); (iv) a lagoonal deposit with a 92 ka age (5b); and (v) a Holocene reef (age >5.7-5.0 ka) (lowest elevation and youngest). A ca 4.9 ka regressive reef (at elevation of 1.5 m above sea-level) is consistent with an island-wide 6.5 m uplift (probably largely coseismic), and a probable further island-wide uplift occurred in the late Holocene. The U-series ages taken together with the heights of transgressive reefs show that uplift since 215 ka was, on average, at similar to0.52 mm/y; although since 5 ka the uplift rate was, on average, similar to1.6 mm/y (the assumption being that a 1.5 m above sea-level reef has a coseismic origin). Elevation of transgressive reefs 5a, 5b and 5c and their ages indicates an island-wide subsidence during the period ?124-89 ka (i.e. Late Quaternary uplift/subsidence was jerky). Late Quaternary uplift/subsidence on the northwest coast of Tanna is considered to be due to irregular thicknesses of crust being subducted beneath Tanna.

The occurrence of hopanoids in planctomycetes: implications for the sedimentary biomarker record

Hopanoids have generally been found in aerobic bacteria (i.e. methanotrophs, heterotrophs and cyanobacteria). Here we show that a variety of hopanoids (i.e. bacteriohopanetetrol, diplopterol, diploptene and a C-27 hopanoid ketone) occur in planctomycetes, including strictly anaerobic bacteria capable of anaerobic ammonium oxidation. Since planctomycetes have a widespread occurrence in Nature, this indicates that they may be an additional source for sedimentary hopanoids. (C) 2004 Elsevier Ltd. All rights reserved.

Petrology and geochemistry of early cretaceous bimodal continental flood volcanism of the NW Etendeka, Namibia. Part 2: Characteristics and petrogenesis of the high-Ti latite and high-Ti and low-Ti voluminous quartz latite eruptives

As a result of their relative concentration towards the respective Atlantic margins, the silicic eruptives of the Parana (Brazil)-Etendeka large igneous province are disproportionately abundant in the Etendeka of Namibia. The NW Etendeka silicic units, dated at similar to132 Ma, occupy the upper stratigraphic levels of the volcanic sequences, restricted to the coastal zone, and comprise three latites and five quartz latites (QL). The large-volume Fria QL is the only low-Ti type. Its trace element and isotopic signatures indicate massive crustal input. The remaining NW Etendeka silicic units are enigmatic high-Ti types, geochemically different from low-Ti types. They exhibit chemical affinities with the temporally overlapping Khumib high-Ti basalt (see Ewart et al. Part 1) and high crystallization temperatures (greater than or equal to980 to 1120degreesC) inferred from augite and pigeonite phenocrysts, both consistent with their evolution from a mafic source. Geochemically, the high-Ti units define three groups, thought genetically related. We test whether these represent independent liquid lines of descent from a common high-Ti mafic parent. Although the recognition of latites reduces the apparent silica gap, difficulty is encountered in fractional crystallization models by the large volumes of two QL units. Numerical modelling does, however, support large-scale open-system fractional crystallization, assimilation of silicic to basaltic materials, and magma mixing, but cannot entirely exclude partial melting processes within the temporally active extensional environment. The fractional crystallization and mixing signatures add to the complexity of these enigmatic and controversial silicic magmas. The existence, however, of temporally and spatially overlapping high-Ti basalts is, in our view, not coincidental and the high-Ti character of the silicic magmas ultimately reflects a mantle signature.

Petrology and geochemistry of early cretaceous bimodal continental flood volcanism of the NW Etendeka, Namibia. Part 1: Introduction, mafic lavas and re-evaluation of mantle source components

The bimodal NW Etendeka province is located at the continental end of the Tristan plume trace in coastal Namibia. It comprises a high-Ti (Khumib type) and three low-Ti basalt (Tafelberg, Kuidas and Esmeralda types) suites, with, at stratigraphically higher level, interstratified high-Ti latites (three units) and quartz latites (five units), and one low-Ti quartz latite. Khumib basalts are enriched in high field strength elements and light rare earth elements relative to low-Ti types and exhibit trace element affinities with Tristan da Cunha lavas. The unradiogenic Pb-206/Pb-204 ratios of Khumib basalts are distinctive, most plotting to the left of the 132 Ma Geochron, together with elevated Pb-207/Pb-204 ratios, and Sr-Nd isotopic compositions plotting in the lower Nd-143/Nd-144 part of mantle array (EM1-like). The low-Ti basalts have less coherent trace element patterns and variable, radiogenic initial Sr (similar to0.707-0.717) and Pb isotope compositions, implying crustal contamination. Four samples, however, have less radiogenic Pb and Sr that we suggest approximate their uncontaminated source. All basalt types, but particularly the low-Ti types, contain samples with trace element characteristics (e.g. Nb/Nb-*) suggesting metasediment input, considered source-related. Radiogenic isotope compositions of these samples require long-term isolation of the source in the mantle and depletions (relative to unmodified sediment) in certain elements (e.g. Cs, Pb, U), which are possibly subduction-related. A geodynamic model is proposed in which the emerging Tristan plume entrained subducted material in the Transition Zone region, and further entrained asthenosphere during plume head expansion. Mixing calculations suggest that the main features of the Etendeka basalt types can be explained without sub-continental lithospheric mantle input. Crustal contamination is evident in most low-Ti basalts, but is distinct from the incorporation of a metasedimentary source component at mantle depths.