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.

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.

Partitioning of deformation within a subduction channel during exhumation of high-pressure rocks: a case study from the Western Alps

The metamorphic belt of the Western Alps was subjected to widespread extensional tectonism at the end of the Eocene (ca. 45-35 Ma). Extension was accommodated by hinterland-directed movements along gently inclined extensional shear zones, which facilitated rapid exhumation of high-pressure and ultra-high-pressure rocks. This deformation resulted in a normal metamorphic sequence. Extension in the inner parts of the Western Alps was coeval with shortening at the front of the belt (foreland-directed thrusts), which took place during decompression, and emplaced higher grade metamorphic units over lower grade metamorphic rocks, thus forming an inverse metamorphic sequence. Two mechanisms for this extensional episode are discussed: (1) collapse of an overthickened lithosphere, and (2) internal readjustments within the orogenic wedge due to subduction channel dynamics. We favour the latter mechanism because it can account for the development of the observed inverse and normal metamorphic sequences along foreland-directed thrusts and hinterland-directed detachments, respectively. This hypothesis is supported by published structural, metamorphic and geochronological data from four geological transects through the Western Alps. This study also emphasizes the importance of post-shearing deformation (e.g. horizontal buckling versus vertical flattening), which can modify the distribution of hinterland- and foreland-directed shear zones in orogenic belts. (c) 2006 Elsevier Ltd. All rights reserved.

Towards an understanding of the age and origin of mesothermal gold mineralisation in the Etheridge Goldfield, Georgetown region, north Queensland

Small mesothermal vein quam-gold-base-metal sulfide deposits from which some 20 t of Au-Ag bullion have been extracted, are the most common gold deposits in the Georgetown region of north Queensland-several hundred were mined or prospected between 1870 and 1950. These deposits are mostly hosted by Proterozoic granitic and metamorphic rocks and are similar to the much larger Charters Towers deposits such as Day Dawn and Brilliant, and in some respects to the Motherlode deposits of California. The largest deposit in the region-Kidston (> 138 t of Au and Ag since 1985)- is substantially different. It is hosted by sheeted quartz veins and cavities in brecciated Silurian granite and Proterozoic metamorphics above nested high-level Carboniferous intrusives associated with a nearby cauldron subsidence structure. This paper provides new information (K-Ar and Rb-Sr isotopic ages, preliminary oxygen isotope and fluid-inclusion data) from some of the mesothermal deposits and compares it with the Kidston deposit. All six dated mesothermal deposits have Siluro-Devonian (about 425 to 400 Ma) ages. All nine of such deposits analysed have delta(18)O quartz values in the range 8.4 to 15.7 parts per thousand, Fluid-inclusion data indicate homogenisation temperatures in the range 230-350 degrees C. This information, and a re-interpretation of the spatial relationships of the deposits with various elements of the updated regional geology, is used to develop a preliminary metallogenic model of the mesothermal Etheridge Goldfield. The model indicates how the majority of deposits may have formed from hydrothermal systems initiated during the emplacement of granitic batholiths that were possibly, but not clearly, associated with Early Palaeozoic subduction, and that these fluid systems were dominated by substantially modified meteoric and/or magmatic fluids. The large Kidston deposit and a few small relatives are of Carboniferous age and formed more directly from magmatic systems much closer to the surface.

Coaxial flattening at deep levels of orogenic belts: evidence from blueschists and eclogites on Syros and Sifnos (Cyclades, Greece)

This work presents new Structural data from a high-pressure/low-temperature (HP/LT) metamorphic terrane exposed on the islands of Syros and Sifnos (Cyclades, Greece). The structure and the metamorphism of a relatively coherent HP/LT rock section were studied in order to elucidate how strain was accommodated at deep crustal levels during the formation and exhumation of HP/LT rocks. At least three deformation phases associated with eclogite- and blueschist-facies conditions (P = 8-15 kbar; T = 400-550 degreesC) were recognised. The earliest deformation fabric (S1), preserved as inclusion trails within garnet porphyroblasts, is aligned to define a sub-vertical schistosity (at present orientation), which is frequently orthogonal to the flat matrix schistosity (S2), and may indicate that deep crustal thickening involved upright folding. The currently dominant fabric in the HP rock section, S2, is Usually moderately dipping and locally contains NW-trending glaucophane lineations, symmetric pressure-shadows and eclogitic boudins. The symmetric structures associated with this fabric seem to indicate coaxial vertical thinning, although the existence of non-coaxial structures out of the study area cannot be excluded. Glaucophane-bearing shear bands (S3), with top-to-NW sense of shearing, locally crosscut the earlier structures. The latest recognised fabric (D4) is scarce and often absent within the HP rocks. It is associated with top-to-NE kinematic criteria that formed at greenschist-facies conditions (P = 4-7 kbar; T = 400-450 degreesC). Based on these observations, it is suggested that partitioning of strain occurred at different crustal levels and at different times. Deep crustal deformation was governed by thickening via upright folding followed by coaxial vertical thinning, whereas non-coaxial shearing occurred when the rocks were already exhumed to relatively shallow crustal levels. The earliest fabrics (D1 to D3) pertain to Alpine orogenesis and possibly to syn-orogenic extension, whereas the latest correspond to whole-crust back-are extension. (C) 2002 Elsevier Science Ltd. All rights reserved.

Effects of hot intrusions on pore-fluid flow and heat transfer in fluid-saturated rocks

We use the finite element method to solve coupled problems between pore-fluid flow and heat transfer in fluid-saturated porous rocks. In particular, we investigate the effects of both the hot pluton intrusion and topographically driven horizontal flow on the distributions of the pore-flow velocity and temperature in large-scale hydrothermal systems. Since general mineralization patterns are strongly dependent on distributions of both the pore-fluid velocity and temperature fields, the modern mineralization theory has been used to predict the general mineralization patterns in several realistic hydrothermal systems. The related numerical results have demonstrated that: (1) The existence of a hot intrusion can cause an increase in the maximum value of the pore-fluid velocity in the hydrothermal system. (2) The permeability of an intruded pluton is one of the sensitive parameters to control the pore-fluid flow, heat transfer and ore body formation in hydrothermal systems. (3) The maximum value of the pore-fluid velocity increases when the bottom temperature of the hydrothermal system is increased. (4) The topographically driven flow has significant effects on the pore-fluid flow, temperature distribution and precipitation pattern of minerals in hydrothermal systems. (5) The size of the computational domain may have some effects on the pore-fluid flow and heat transfer, indicating that the size of a hydrothermal system may affect the pore-fluid flow and heat transfer within the system. (C) 2003 Elsevier Science B.V. All rights reserved.

An equivalent algorithm for simulating thermal effects of magma intrusion problems in porous rocks

An equivalent algorithm is proposed to simulate thermal effects of the magma intrusion in geological systems, which are composed of porous rocks. Based on the physical and mathematical equivalence, the original magma solidification problem with a moving boundary between the rock and intruded magma is transformed into a new problem without the moving boundary but with a physically equivalent heat source. From the analysis of an ideal solidification model, the physically equivalent heat source has been determined in this paper. The major advantage in using the proposed equivalent algorithm is that the fixed finite element mesh with a variable integration time step can be employed to simulate the thermal effect of the intruded magma solidification using the conventional finite element method. The related numerical results have demonstrated the correctness and usefulness of the proposed equivalent algorithm for simulating the thermal effect of the intruded magma solidification in geological systems. (C) 2003 Elsevier B.V. All rights reserved.

Offspring size affects the post-metamorphic performance of a colonial marine invertebrate

The positive relationship between offspring size and offspring fitness is a fundamental assumption of life-history theory, but it has received relatively little attention in the marine environment. This is surprising given that substantial intraspecific variation in offspring size is common in marine organisms and there are clear links between larval experience and adult performance. The metamorphosis of most marine invertebrates does not represent a newbeginning, and larval experiences can have effects that carry over to juvenile survival and growth. We show that larval size can have equally important carryover effects in a colonial marine invertebrate. In the bryozoan Bugula neritina, the size of the non-feeding larvae has a prolonged effect on colony performance after metamorphosis. Colonies that came from larger larvae survived better, grew faster, and reproduced sooner or produced more embryos than colonies that came from smaller larvae. These effects crossed generations, with colonies from larger larvae themselves producing larger larvae. These effects were found in two populations (in Australia and in the United States) in contrasting habitats.

Variable effects of larval size on post-metamorphic performance in the field

Larval quality may be capable of explaining much of the variation in the recruitment and subsequent population dynamics of benthic marine invertebrates. Whilst the effects of larval nutritional condition on adult performance have received the most attention, recent work has shown that larval size may also be an important and ubiquitous source of variation in larval quality. We examined the effects of variation in larval size on the post-metamorphic survival and growth of Watersipora subtorquata in 2 very different habitats - experimental substrata and pier pilings. We found strong effects of larval size on colony performance, although these varied among experiments. For colonies on experimental substrata, larval size positively affected adult survival and, initially, growth. However, after 3 wk in the field, there was no relationship between larval size and colony size, possibly because colonies were completely surrounded by newly settled organisms. Larval size also positively affected post-metamorphic growth of colonies on pier pilings, but, surprisingly, colonies that came from larger larvae had lower survival than colonies from smaller larvae. Overall, variation in larval size will strongly affect the recruitment and subsequent performance of adults in this species, although this may vary among different habitats. This study highlights the importance of examining the effects of larval quality on adult performance in as realistic conditions as possible, because of the strong interaction between larval size effects and the environment.

Numerical modelling of chemical effects of magma solidification problems in porous rocks

The solidification of intruded magma in porous rocks can result in the following two consequences: (1) the heat release due to the solidification of the interface between the rock and intruded magma and (2) the mass release of the volatile fluids in the region where the intruded magma is solidified into the rock. Traditionally, the intruded magma solidification problem is treated as a moving interface (i.e. the solidification interface between the rock and intruded magma) problem to consider these consequences in conventional numerical methods. This paper presents an alternative new approach to simulate thermal and chemical consequences/effects of magma intrusion in geological systems, which are composed of porous rocks. In the proposed new approach and algorithm, the original magma solidification problem with a moving boundary between the rock and intruded magma is transformed into a new problem without the moving boundary but with the proposed mass source and physically equivalent heat source. The major advantage in using the proposed equivalent algorithm is that a fixed mesh of finite elements with a variable integration time-step can be employed to simulate the consequences and effects of the intruded magma solidification using the conventional finite element method. The correctness and usefulness of the proposed equivalent algorithm have been demonstrated by a benchmark magma solidification problem. Copyright (c) 2005 John Wiley & Sons, Ltd.

A sensitive vorticity gauge using rotated porphyroblasts, and its application to rocks adjacent to the Alpine Fault, New Zealand

Variable aspect ratio porphyroblasts deformed in non-coaxial flow. and internally containing rotated relicts of an external foliation, can be used to characterise plane strain flow regimes. The distribution obtained by plotting the orientation of the long axis of such grains, classified by aspect ratio, against the orientation of the internal foliation is potentially a sensitive gauge of both the bulk shear strain (as previously suggested) and kinematic vorticity number. We illustrate the method using rotated biotite porphyroblasts in the Alpine Schist: a sequence of mid-crustal rocks that have been ramped to the surface along the Alpine Fault. a major transpressional plate boundary. Results indicate that, at distances greater than or equal to similar to1 km from the fault, the rocks have undergone a combination of irrotational fattening and dextral-oblique, normal-sense shear, with a bulk shear strain of similar to0.6 and kinematic vorticity number of similar to0.2. The vorticity analysis is compatible with estimates of strongly oblate bulk strain of similar to 75% maximum shortening. Dextral-reverse transpressional flow characterises higher strain S-tectonite mylonite within similar to1 km of the Alpine Fault. These relationships provide insight into the kinematics of flow and distribution of strain in the hangingwall of the Alpine Fault and place constraints on numerical mechanical models for the exhumation of these mid-crustal rocks. (C) 2001 Elsevier Science Ltd. All rights reserved.

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.