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.

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.

Assigning dates to thin gneissic veins in high-grade metamorphic terranes: A cautionary tale from Akilia, southwest Greenland

A granodiorite from Akilia, southwest Greenland, previously suggested to date putative life-bearing rocks to greater than or equal to3.84 Ga, is re-investigated using whole-rock major and trace-element geochemistry, and detailed cathodoluminescence image-guided secondary ion mass spectrometer analyses of zircon U-Th-Pb and rare earth elements. Complex zircon internal structure reveals three episodes of zircon growth and/or recrystallization dated to c. 3.84 Ga, 3.62 Ga and 2.71 Ga. Rare earth element abundances imply a significant role for garnet in zircon generation at 3.62 Ga and 2.71 Ga. The 3.62 Ga event is interpreted as partial melting of a c. 3.84 Ga grey gneiss precursor at granulite facies with residual garnet. Migration of this 3.62 Ga magma (or melt-crystal mush) away from the melt source places a maximum age limit on any intrusive relationship. These early Archaean relationships have been complicated further by isotopic reworking in the 2.71 Ga event, which could have included a further episode of partial melting. This study highlights a general problem associated with dating thin gneissic veins in polyphase metamorphic terranes, where field relationships may be ambiguous and zircon inheritance can be expected.

Chemical characterization of earth's most ancient clastic metasediments from the Isua Greenstone Belt, southern West Greenland

New and published major and trace element abundances of elastic metasediments (mainly garnet-biotite-plagioclase schists) from the similar to 3.8 Ga Isua Greenstone Belt (IGB), southern West Greenland, are used in an attempt to identify the compositional characteristics of the protoliths of these sediments. Compositionally, the metasediments are heterogeneous with enrichment of LREE (La/Sm-chord = 1.1-3.9) and variable enrichment and depletion of HREE (Gd/Yb-chord = 0.8-4.3). Chondrite-normalized Eu is also variable, spanning a range from relative Eu depletion to enrichment (Eu/Eu* = 0.6-1.3). A series of geochemical and geological criteria provides conclusive evidence for a sedimentary origin, in disagreement with some previous studies that questioned the presence of genuine elastic metasediments. In particular, trace element systematics of IGB metasediments show strong resemblance to other well-documented Archaean clastic sediments, and are consistent with a provenance consisting of ultramafic, malic and felsic igneous rocks. Two schists, identified as metasomatized mafic igneous rocks from petrographic and field evidence, show distinct compositional differences to the metasediments. Major element systematics document incipient-to-moderate source weathering in the majority of metasediments, while signs of secondary K-addition are rare. Detailed inspection of Eu/Eu*, Fe2O3 and CIW (chemical index of weathering) relationships reveals that elevated iron contents (when compared to averages for continental crust) and strong relative enrichment in Eu may be due to precipitation of marine Fe-oxyhydroxides during deposition or diagenesis on the seafloor. Some of the IGB metasediments have yielded anomalous Nd-142 and W-182 isotopic compositions that were respectively interpreted in terms of early mantle differentiation processes and the presence of a meteorite component. Alternatively, W and possibly Nd isotopes could have been affected by thermal neutron capture on the Hadean surface. The latter process was tested in this study by analysis of Sm isotope compositions, which serve as an effective monitor for neutron capture effects. As no anomalous variation from terrestrial values was detected, we infer that isotope systematics (including W-182 and Nd-142) of IGB metasediments were not affected by neutron capture, but reflect decay of radioactive parent isotopes. Copyright (c) 2005 Elsevier Ltd.

Finite element analysis of steady-state natural convection problems in fluid-saturated porous media heated from below

In this paper, a progressive asymptotic approach procedure is presented for solving the steady-state Horton-Rogers-Lapwood problem in a fluid-saturated porous medium. The Horton-Rogers-Lapwood problem possesses a bifurcation and, therefore, makes the direct use of conventional finite element methods difficult. Even if the Rayleigh number is high enough to drive the occurrence of natural convection in a fluid-saturated porous medium, the conventional methods will often produce a trivial non-convective solution. This difficulty can be overcome using the progressive asymptotic approach procedure associated with the finite element method. The method considers a series of modified Horton-Rogers-Lapwood problems in which gravity is assumed to tilt a small angle away from vertical. The main idea behind the progressive asymptotic approach procedure is that through solving a sequence of such modified problems with decreasing tilt, an accurate non-zero velocity solution to the Horton-Rogers-Lapwood problem can be obtained. This solution provides a very good initial prediction for the solution to the original Horton-Rogers-Lapwood problem so that the non-zero velocity solution can be successfully obtained when the tilted angle is set to zero. Comparison of numerical solutions with analytical ones to a benchmark problem of any rectangular geometry has demonstrated the usefulness of the present progressive asymptotic approach procedure. Finally, the procedure has been used to investigate the effect of basin shapes on natural convection of pore-fluid in a porous medium. (C) 1997 by John Wiley & Sons, Ltd.

Determination of pore size distributions by regularization and finite element collocation

The problem of extracting pore size distributions from characterization data is solved here with particular reference to adsorption. The technique developed is based on a finite element collocation discretization of the adsorption integral, with fitting of the isotherm data by least squares using regularization. A rapid and simple technique for ensuring non-negativity of the solutions is also developed which modifies the original solution having some negativity. The technique yields stable and converged solutions, and is implemented in a package RIDFEC. The package is demonstrated to be robust, yielding results which are less sensitive to experimental error than conventional methods, with fitting errors matching the known data error. It is shown that the choice of relative or absolute error norm in the least-squares analysis is best based on the kind of error in the data. (C) 1998 Elsevier Science Ltd. All rights reserved.

Finite element modelling of temperature gradient driven rock alteration and mineralization in porous rock masses

We present finite element simulations of temperature gradient driven rock alteration and mineralization in fluid saturated porous rock masses. In particular, we explore the significance of production/annihilation terms in the mass balance equations and the dependence of the spatial patterns of rock alteration upon the ratio of the roll over time of large scale convection cells to the relaxation time of the chemical reactions. Special concepts such as the gradient reaction criterion or rock alteration index (RAI) are discussed in light of the present, more general theory. In order to validate the finite element simulation, we derive an analytical solution for the rock alteration index of a benchmark problem on a two-dimensional rectangular domain. Since the geometry and boundary conditions of the benchmark problem can be easily and exactly modelled, the analytical solution is also useful for validating other numerical methods, such as the finite difference method and the boundary element method, when they are used to dear with this kind of problem. Finally, the potential of the theory is illustrated by means of finite element studies related to coupled flow problems in materially homogeneous and inhomogeneous porous rock masses. (C) 1998 Elsevier Science S.A. All rights reserved.

Binding of YY1 and Oct1 to a novel element that downregulates expression of IL-5 in human T cells

Background: IL-5 controls development of eosinophilia and has been shown to be involved in the pathogenesis of allergic diseases. In both atopic and nonatopic asthma, elevated IL-5 has been detected in peripheral blood and the airways. IL-5 is produced mainly by activated T cells, and its expression is regulated at the transcriptional level. Objective: This study focuses on the functional analysis of the human IL-5 (hIL-5) promoter and characterization of eis-regulatory elements and transcription factors involved in the suppression of IL-5 transcription in T cells. Methods: Methods used in this study include DNase I footprint assays, electrophoretic mobility shift assays, and functional analysis by mammalian cell transfection involving deletion analysis and site-directed mutagenesis. Results: We identified 5 protein binding regions (BRs) located within the proximal hIL-5 promoter. Functional analysis indicates that the BRs are involved in control of hIL-5 promoter activity. Two of these regions, BR3 and BR4 located at positions -102 to -73, have not previously been described as regulators of IL-5 expression in T cells. We show that the BR3 sequence contains a novel negative regulatory element located at positions -90 to -79 of the hIL-5 promoter, which binds Oct1, octamer-like, and YY1 nuclear factors. Substitution mutations, which abolished binding of these proteins to the BR3 sequence, significantly increased hIL-5 promoter activity in activated T cells. Conclusion: We suggest that Oct1, YY1, and octamer-like factors binding to the -90/-79 sequence within the proximal IL-5 promoter are involved in suppression of IL-5 transcription in T cells.

Finite element analysis of flow patterns near geological lenses in hydrodynamic and hydrothermal systems

We use theoretical and numerical methods to investigate the general pore-fluid flow patterns near geological lenses in hydrodynamic and hydrothermal systems respectively. Analytical solutions have been rigorously derived for the pore-fluid velocity, stream function and excess pore-fluid pressure near a circular lens in a hydrodynamic system. These analytical solutions provide not only a better understanding of the physics behind the problem, but also a valuable benchmark solution for validating any numerical method. Since a geological lens is surrounded by a medium of large extent in nature and the finite element method is efficient at modelling only media of finite size, the determination of the size of the computational domain of a finite element model, which is often overlooked by numerical analysts, is very important in order to ensure both the efficiency of the method and the accuracy of the numerical solution obtained. To highlight this issue, we use the derived analytical solutions to deduce a rigorous mathematical formula for designing the computational domain size of a finite element model. The proposed mathematical formula has indicated that, no matter how fine the mesh or how high the order of elements, the desired accuracy of a finite element solution for pore-fluid flow near a geological lens cannot be achieved unless the size of the finite element model is determined appropriately. Once the finite element computational model has been appropriately designed and validated in a hydrodynamic system, it is used to examine general pore-fluid flow patterns near geological lenses in hydrothermal systems. Some interesting conclusions on the behaviour of geological lenses in hydrodynamic and hydrothermal systems have been reached through the analytical and numerical analyses carried out in this paper.

Finite element modelling of reactive mass transport problems in fluid-saturated porous media

We use the finite element method to solve reactive mass transport problems in fluid-saturated porous media. In particular, we discuss the mathematical expression of the chemical reaction terms involved in the mass transport equations for an isothermal, non-equilibrium chemical reaction. It has turned out that the Arrhenius law in chemistry is a good mathematical expression for such non-equilibrium chemical reactions especially from the computational point of view. Using the finite element method and the Arrhenius law, we investigate the distributions of PH (i.e. the concentration of H+) and the relevant reactive species in a groundwater system. Although the main focus of this study is on the contaminant transport problems in groundwater systems, the related numerical techniques and principles are equally applicable to the orebody formation problems in the geosciences. Copyright (C) 1999 John Wiley & Sons, Ltd.

Finite element modelling of dissipative structures for nonequilibrium chemical reactions in fluid-saturated porous media

We use the finite element method to model and predict the dissipative structures of chemical species for a nonequilibrium chemical reaction system in a fluid-saturated porous medium. In particular, we explore the conditions under which dissipative structures of the species may exist in the Brusselator type of nonequilibrium chemical reaction. Since this is the first time the finite element method and related strategies have been used to study the chemical instability problems in a fluid-saturated porous medium, it is essential to validate the method and strategies before they are put into application. For this purpose, we have rigorously derived the analytical solutions for dissipative structures of chemical species in a benchmark problem, which geometrically is a square. Comparison of the numerical solutions with the analytical ones demonstrates that the proposed numerical method and strategy are robust enough to solve chemical instability problems in a fluid-saturated porous medium. Finally, the related numerical results from two application examples indicate that both the regime and the magnitude of pore-fluid flow have significant effects on the nature of the dissipative structures that developed for a nonequilibrium chemical reaction system in a fluid-saturated porous medium. The motivation for this study is that self-organization under conditions of pore-fluid flow in a porous medium is a potential mechanism of the orebody formation and mineralization in the upper crust of the Earth. (C) 2000 Elsevier Science S.A. All rights reserved.

Finite element analysis of heat transfer and mineralization in layered hydrothermal systems with upward throughflow

We use the finite element method to solve the coupled problem between convective pore-fluid flow, heat transfer and mineralization in layered hydrothermal systems with upward throughflow. In particular, we present the improved rock alteration index (IRAI) concept for predicting the most probable precipitation and dissolution regions of gold (Au) minerals in the systems. To validate the numerical method used in the computation, analytical solutions to a benchmark problem have been derived. After the numerical method is validated, it is used to investigate the pattern of pore-fluid Aom, the distribution of temperature and the mineralization pattern of gold minerals in a layered hydrothermal system with upward throughflow. The related numerical results have demonstrated that the present concept of IRAI is useful and applicable for predicting the most probable precipitation and dissolution regions of gold (Au) minerals in hydrothermal systems. (C) 2000 Elsevier Science S.A. All rights reserved.

Finite element modelling of three-dimensional convection problems in fluid-saturated porous media heated from below

We use the finite element method to model three-dimensional convective pore-fluid flow in fluid-saturated porous media when they are heated from below. In particular, we employ the particle-tracking technique to mimic the trajectories of particles in three-dimensional fluid flow problems. The related numerical results demonstrated that: (1) The progressive asymptotic approach procedure, which was previously developed for the finite element modelling of two-dimensional convective pore-fluid flow problems, is equally applicable to the finite element modelling of three-dimensional convective pore-fluid flow in fluid-saturated porous media heated from below. (2) The perturbation of gravity at different planes has a significant effect on the pattern of three-dimensional convective pore-fluid flow and therefore, may influence the pattern of orebody formation and mineralization in three-dimensional hydrothermal systems. Copyright (C) 2001 John Wiley & Sons, Ltd.

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.