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.

Fluid-rock interaction during late stages of the Alpine exhumation

Résumé de la thèseLa fracturation des roches au cours de phases compressives ou extensives est un souvent évoquée pour expliquer la circulation de fluide au sein des roches cristallines. Dans le cadre de cette thèse, la circulation des fluides lors de l'exhumation tardive des Alpes a été étudiée en utilisant deux approches différentes: analyses structurales de la déformation fragile d'une part et analyses géochimiques des roches et des minéraux (isotopes stables, datations U/Pb, thermochronologie (U-Th)/He) d'autre part. Cette approche combinée a permis de mieux comprendre l'interaction existante entre les fluides métamorphiques et les fluides météoriques, ainsi que leur interaction avec les roches encaissantes. Le travail a été effectué dans la zone Pennique du Valais suisse.La première partie était focalisée sur la déformation fragile, le but étant de définir les différents types de déformations existantes et de déterminer l'âge relatif des différentes familles de failles. Dans la région d'étude, quatre domaines ont été distingués. Chacun d'eux comportent deux types de structures fragiles, certaines sont minéralisées alors que d'autre non. Au sein de chaque domaine, la direction principale des structures minéralisées correspond à l'orientation des accidents tectoniques majeurs de la région (Aosta- Ranzola Line au Sud, Rhône Line au Nord et Simplon Fault Zone à l'Est), alors que les structures non- minéralisées montrent des orientations plus variables. Ainsi, le premier type de structure est interprété comme résultant d'une dislocation tectonique alors que le deuxième type de structure résulterait d'une dislocation gravitaire locale. Il n'est néanmoins pas possible de classer chronologiquement la formation de ces deux types de structure ni d'attribuer un âge relatif aux changements d'orientation des contraintes majeures.La deuxième étude a été effectuée dans la région de la zone de faille du Simplon. Dans cette zone, la composition isotopique des minéraux ayant cristallisé à l'intérieur des fractures tardives permet de distinguer différents types de circulation de fluide. Les valeurs δ180 du quartz de la roche encaissante ainsi que ceux des veines tardives du bloque inférieur de la faille sont comparables. Ces valeurs indiquent un rééquilibrage et un tamponnage isotopique des fluides tardifs au contact de la roche encaissante lors de la fracturation de cette dernière et de la cristallisation des veines tardives. La même situation est observée dans la partie nord du bloque supérieur ainsi que dans sa partie sud. Ceci n'est néanmoins pas le cas pour la partie centrale du bloque supérieur où les valeurs isotopiques des minéraux dans les veines tardives sont approximativement 3 %o plus basses (avec des valeurs extrêmes négatifs), indiquant une contribution d'eau météorique aux fluides circulant dans les veines. Ces données suggèrent qu'une infiltration d'eau météorique a pu avoir lieu dans le bloque supérieur, où la fracturation des roches est plus intensive car le déplacement relatif le long de la faille y fut plus important, et la température maximale du métamorphisme plus basse. La troisième contribution traite de la géo-thermochronologie de la zone de contact entre la klippe de la Dent Blanche et la nappe de Tsaté. De petits zircons euhédraux ont été trouvés dans un plan de faille minéralisé (parallèle à la Faille du Rhône, voir première partie de l'étude), riche en hématite et quartz, de la zone d'étude. Les analyses U/Pb donnent des âges radiométriques autour de 270 - 280 Ma aux zircons extraits de la minéralisation ainsi que ceux extraits de la roche encaissante, ce qui correspond à l'âge de la nappe de la Dent Blanche et non celui de la nappe du Tsaté qui est elle-même classiquement interprétée comme une ophiolite Jurassique de l'Océan Liguro-Piémontais. Ces données suggèrent que les zircons contenus dans la veine ont été hérités de la roche encaissante. Les résultats (U-Th)/He indiquent un âge de refroidissement différent pour la roche encaissante (25.5 ± 2.0 Ma) que celui de la minéralisation (17.7 ±1.4 Ma). Le thermomètre isotopique quartz-hématite indique une température d'équilibre, et donc de mise en place de la minéralisation, d'environ 170 °C, température très proche de la température de -180 °C de fermeture du zircon pour le système (U-Th)/He. Ceci suggère que l'âge de refroidissement des zircons de la minéralisation correspond aussi à l'âge de formation de la faille.Thesis abstractFluid circulation in fractured rocks is a common process in geology, and it is generally the consequence of faulting and fracturing during both tectonic compression and extension. This thesis is focused on fluid circulation during late stages of the Alpine exhumation. After a structural analysis of the late brittle deformation of the studied samples, several analytical methods (stable isotope investigations, U/Pb radiometric dating, (U-Th)/He thermochronology) have been applied to understand the interaction of metamorphic and meteoric fluids with one another as well as with the host rock. This thesis is articulated around three study directions. All studies were conducted in the Penninic Zone of the Valais, Switzerland. The first study deals with late, brittle deformation and focuses on the different deformation styles and on the relative age of the different families of fractures. In order to do this, late brittle structures observed in four different domains have been subdivided as a function of the existence (or not) and type of mineralization. Comparisons between mineralized and non-mineralized strike directions for all four domains show that mineralized structures follow the strike orientation of major tectonic movements indicated in the Penninic Zone of the Valais (Aosta-Ranzola Line to the S, Rhône Line to the Ν and Simplon Fault Zone to the E), whereas non-mineralized fractures have a more variable strike orientation. This difference could be interpreted as indicative of tectonic-related faulting (mineralized structures) vs. local, collapse-related faulting (non-mineralized fractures), but it is not strong enough to indicate a relative age of the late brittle structures, and/or a change in the orientation of the strain field in post-Miocene times. The second studied area is focused on the Simplon Fault Zone (SFZ). Stable isotope analyses of minerals filling these late fractures indicate that there are two different fluid circulation systems in the footwall and hanging wall of the SFZ. In the footwall, δ180 values of quartz from both the host rock and the late veins range from +10 %o to +12 %o. This is consistent with buffering of circulating fluids by the host rock during fracturing and vein precipitation. In the hanging wall, δΙ80 values for quartz crystals from the host rock and the late veins are similar in both the northern and southern parts of the detachment that are both affected by the same degree of metamorphism (greenschist to the Ν and amphibolite to the S). This is not the case in the central part of the SFZ, where there is a jump from amphibolite facies in the footwall to greenschist facies in the hanging wall. δ,80 values for quartz from the hanging wall late veins are approximately 3.0 %o lower (down to negative values in some cases) than the values observed in the footwall These data suggest that infiltration of meteoric water may have occurred in the most fractured parts of the hanging wall, where relative displacement on the SFZ was the greatest and the peak temperature lower. In the less fractured footwall the δ180 values reflect a host rock-buffered system.The third study is focused on geo-thermochronology at the contact between the Dent Blanche nappe and the Tsaté nappe where small, euhedral zircons were found in a hematite- and quartz-rich mineralization on a late normal fault plane parallel to the Rhône Line (see first part of the study). U/Pb analysis indicates that the zircons - both in the late mineralization and in the host rock - have absolute radiometric ages clustering around 270 - 280 Ma, which is the accepted age for intrusive rocks from the Austroalpine Dent Blanche units but not for the Tsaté nappe. The latter is classically interpreted as an ophiolitic remnant of the Jurassic Liguro-Piemontais Ocean. U/Pb analyses suggest that zircons in late mineralization are all inherited from the host rock; however, results of (U-Th)/He analyses indicate that cooling ages for the host rocks are different to the cooling ages for the zircons in late mineralization. Indeed, the calculated cooling age for the Arolla gneiss is 25.5 ± 2.0 Ma, whilst the cooling age for the associated mineralized fault plane is 17.7 ±1.4 Ma. Oxygen stable isotope fractionation between quartz and hematite in the same late mineralization corresponds to temperatures of about 170 °C. The proximity of the calculated emplacement temperature for the mineralization and the lower accepted closure temperature for zircon in the (U-Th)/He system (-180 °C) imply that the age of 17.7 ± 1.4 Ma can also be interpreted as the formation age of this late brittle fault.Résumé grand publicLa circulation des fluides dans les roches fracturées est typique de nombreux processus géologiques, et très souvent est la conséquence de la fracturation des roches. Cette thèse aborde la question de la circulation des fluides pendant les dernières phases du soulèvement des Alpes. Après une analyse structurale de la fracturation directement sur le terrain, plusieurs méthodes géochimiques ont été appliquées pour comprendre l'interaction entre les différents fluides circulants, et avec leur propre roche mère. L'étude, concentrée sur trois directions principales, a été conduite dans la zone Pennique du Valais suisse. La première partie traite de la déformation cassante dans le secteur cité. L'analyse détaillée des fractures a permis de les subdiviser en structures minéralisées et non-minéralisées, sur quatre domaines différents. La comparaison entre les directions des structures minéralisées et non-minéralisées a permis de montrer que les premières suivent l'orientation des accidents tectoniques majeurs de la région, alors que les structures non- minéralisées ont une orientation plus variable. Cette différence pourrait être interprétée comme indication d'une dislocation tectonique (structures minéralisées) contre une dislocation gravitaire locale (structures non-minéralisées), mais elle n'est pas assez forte pour indiquer un âge relatif des structures tardives et/ou un changement de l'orientation des contraintes après -20 Ma vers le présent.A partir de ces observations, la deuxième étude est concentrée dans la région de la faille du Simplon. Les analyses géochimiques sur les minéraux remplissant les structures tardives indiquent qu'il y a deux différents systèmes de circulation des fluides dans les deux parties (toit et mur) de la faille. Dans le mur, les valeurs isotopiques des minéraux cristallisés à partir d'un fluide tardif sont les mêmes de ceux de la roche mère, donc il y a eu rééquilibration chimique entre fluide et roche pendant la fracturation de cette dernière et la précipitation des minéraux. Dans le toit, les valeurs isotopiques dans la roche mère et dans les minéraux des veines tardives sont comparables dans les parties Ν et S de la faille, où les roches du toit et du mur ont atteint une température maximale - pendant phase prograde de la formation des Alpes - comparable. Au contraire, dans la partie centrale, où le mur a atteint des températures maximales plus élevées par rapport au toit, les valeurs géochimiques des minéralisations tardives du toit sont parfois plus basses que les valeurs observées dans le mur. Ces données suggèrent que l'infiltration de l'eau de surface aurait pu se produire dans la partie plus fracturée du toit, où le déplacement relatif le long de la faille était majeur et les températures maximales mineures. Au contraire, les données géochimiques du mur de la partie centrale indiquent un système isotopique équilibré par la roche mère.La troisième partie de ce travail se base sur l'étude géochimique intégrée des isotopes stables d'Oxygène et radioactifs du Plomb, Uranium, Thorium et Hélium, auprès d'une faille normale minéralisée et des roches de la région à cheval entre deux nappes, la nappe de la Dent Blanche et la nappe de Tsaté. Ici, des petits zircons ont été trouvés dans la minéralisation citée, riche en hématite et quartz. L'analyse radiométrique Uranium/Plomb a montré que les zircons dans la minéralisation et dans les roches autour ont des âges comparables (autour 280 Ma). Cela signifie que les zircons dans la minéralisation tardive ont été hérités de la roche mère pendant la fracturation et la circulation des fluides tardives. De l'autre coté, les résultats des analyses Uranium-Thorium/Hélium indiquent que les âges de refroidissement pour les roches mères sont différents comparés aux âges de refroidissement pour les zircons dans la minéralisation tardive: ces derniers sont plus jeunes d'environ 8 Ma (autour 25 Ma et autour 17 Ma respectivement). Les analyses des isotopes de l'oxygène sur quartz et hématite dans la même minéralisation donnent une température de mise en place de cette dernière d'environ 170° C. La température de fermeture du système chimique des zircons dans le système (Uranium-Thorium)/Hélium est d'environ 180 °C: la proximité de ces deux températures implique que l'âge de refroidissement de la minéralisation tardive peut également être interprété comme âge de formation de la faille.

Dating brittle tectonic movements with cleft monazite: Fluid-rock interaction and formation of REE minerals

[1] Two millimeter-sized hydrothermal monazites from an open fissure (cleft) that developed late during a dextral transpressional deformation event in the Aar Massif, Switzerland, have been investigated using electron microprobe and ion probe. The monazites are characterized by high Th/U ratios typical of other hydrothermal monazites. Deformation events in the area have been subdivided into three phases: (D1) main thrusting including formation of a new schistosity, (D2) dextral transpression, and (D3) local crenulation including development of a new schistosity. The two younger deformational structures are related to a subvertically oriented intermediate stress axis, which is characteristic for strike slip deformation. The inferred stress environment is consistent with observed kinematics and the opening of such clefts. Therefore, the investigated monazite-bearing cleft formed at the end of D2 and/or D3, and during dextral movements along NNW dipping planes. Interaction of cleft-filling hydrothermal fluid with wall rock results in rare earth element (REE) mineral formation and alteration of the wall rock. The main newly formed REE minerals are Y-Si, Y-Nb-Ti minerals, and monazite. Despite these mineralogical changes, the bulk chemistry of the system remains constant and thus these mineralogical changes require redistribution of elements via a fluid over short distances (centimeter). Low-grade alteration enables local redistribution of REE, related to the stability of the accessory phases. This allows high precision isotope dating of cleft monazite. 232Th/208Pb ages are not affected by excess Pb and yield growth domain ages between 8.03 ± 0.22 and 6.25 ± 0.60 Ma. Monazite crystallization in brittle structures is coeval or younger than 8 Ma zircon fission track data and hence occurred below 280°C.

Fluid-rock-strain feedbacks in extensional shear zones

La présence de fluide météorique synchrone à l'activité du détachement (Farmin, 2003 ; Mulch et al., 2007 ; Gébelin et al., 2011), implique que les zones de cisaillement sont des systèmes ouverts avec des cellules de convections à l'échelle crustale et un intense gradient géothermique au sein du détachement (Morrison et Anderson, 1998, Gottardi et al., 2011). De plus, les réactions métamorphiques liées à des infiltrations fluides dans les zones de cisaillement extensionnel peuvent influencer les paramètres rhéologiques du système (White and Knipe, 1978), et impliquer la localisation de la déformation dans la croûte. Dans ce manuscrit, deux zones de cisaillement infiltrées par des fluides météoriques sont étudiées, l'une étant largement quartzitique, et l'autre de nature granitique ; les relations entre déformation, fluides, et roches s'appuient sur des approches structurales, microstructurales, chimiques et isotopiques. L'étude du détachement du Columbia river (WA, USA) met en évidence que la déformation mylonitique se développe en un million d'années. La phase de cisaillement principal s'effectue à 365± 30°C d'après les compositions isotopiques en oxygène du quartz et de la muscovite. Ces minéraux atteignent l'équilibre isotopique lors de leur recristallisation dynamique contemporaine à la déformation. La zone de cisaillement enregistre une baisse de température, remplaçant le mécanisme de glissement par dislocation par celui de dissolution- précipitation dans les derniers stades de l'activité du détachement. La dynamique de circulation fluide bascule d'une circulation pervasive à chenalisée, ce qui engendre localement la rupture des équilibres d'échange isotopiques. La zone de cisaillement de Bitterroot (MT, USA) présente une zone mylonitique de 600m d'épaisseur, progressant des protomylonites aux ultramylonites. L'intensité de la localisation de la déformation se reflète directement sur l'hydratation des feldspaths, réaction métamorphique majeure dite de « rock softening ». Une étude sur roche totale indique des transferts de masse latéraux au sein des mylonites, et d'importantes pertes de volume dans les ultramylonites. La composition isotopique en hydrogène des phyllosilicates met en évidence la présence (1) d'une source magmatique/métamorphique originelle, caractérisée par les granodiorites ayant conservé leur foliation magmatique, jusqu'aux protomylonites, et (2) une source météorique qui tamponne les valeurs des phyllosilicates des fabriques mylonitiques jusqu'aux veines de quartz non-déformées. Les compositions isotopiques en oxygène des minéraux illustrent le tamponnement de la composition du fluide météorique par l'encaissant. Ce phénomène cesse lors du processus de chloritisation de la biotite, puisque les valeurs des chlorites sont extrêmement négatives (-10 per mil). La thermométrie isotopique indique une température d'équilibre isotopique de la granodiorite entre 600-500°C, entre 500-300°C dans les mylonites, et entre 300 et 200°C dans les fabriques cassantes (cataclasites et veines de quartz). Basé sur les résultats issus de ce travail, nous proposons un modèle général d'interactions fluide-roches-déformation dans les zones de détachements infiltrées par des fluides météoriques. Les zones de détachements évoluent rapidement (en quelques millions d'années) au travers de la transition fragile-ductile ; celle-ci étant partiellement contrôlée par l'effet thermique des circulations de fluide météoriques. Les systèmes de détachements sont des lieux où la déformation et les circulations fluides sont couplées ; évoluant rapidement vers une localisation de la déformation, et de ce fait, une exhumation efficace. - The presence of meteoric fluids synchronous with the activity of extensional detachment zones (Famin, 2004; Mulch et al., 2007; Gébelin et al., 2011) implies that extensional systems involve fluid convection at a crustal scale, which results in high geothermal gradients within active detachment zones (Morrison and Anderson, 1998, Gottardi et al., 2011). In addition, the metamorphic reactions related to fluid infiltration in extensional shear zones can influence the rheology of the system (White and Knipe, 1978) and ultimately how strain localizes in the crust. In this thesis, two shear zones that were permeated by meteoric fluids are studied, one quartzite-dominated, and the other of granitic composition; the relations between strain, fluid, and evolving rock composition are addressed using structural, microstructural, and chemical/isotopic measurements. The study of the Columbia River detachment that bounds the Kettle core complex (Washington, USA) demonstrates that the mylonitic fabrics in the 100 m thick quartzite- dominated detachment footwall developed within one million years. The main shearing stage occurred at 365 ± 30°C when oxygen isotopes of quartz and muscovite equilibrated owing to coeval deformation and dynamic recrystallization of these minerals. The detachment shear zone records a decrease in temperature, and dislocation creep during detachment shearing gave way to dissolution-precipitation and fracturing in the later stages of detachment activity. Fluid flow switched from pervasive to channelized, leading to isotopic disequilibrium between different minerals. The Bitterroot shear zone detachment (Montana, USA) developed a 600 m thick mylonite zone, with well-developed transitions from protomylonite to ultramylonite. The localization of deformation relates directly to the intensity of feldspar hydration, a major rock- softening metamorphic reaction. Bulk-rock analyses of the mylonitic series indicate lateral mass transfer in the mylonite (no volume change), and significant volume loss in ultramylonite. The hydrogen isotope composition of phyllosilicates shows (1) the presence of an initial magmatic/metamorphic source characterized by the granodiorite in which a magmatic, and gneissic (protomylonite) foliation developed, and (2) a meteoric source that buffers the values of phyllosilicates in mylonite, ultramylonite, cataclasite, and deformed and undeformed quartz veins. The mineral oxygen isotope compositions were buffered by the host-rock compositions until chloritization of biotite started; the chlorite oxygen isotope values are negative (-10 per mil). Isotope thermometry indicates a temperature of isotopic equilibrium of the granodiorite between 600-500°C, between 500-300°C in the mylonite, and between 300 and 200°C for brittle fabrics (cataclasite and quartz veins). Results from this work suggest a general model for fluid-rock-strain feedbacks in detachment systems that are permeated by meteoric fluids. Phyllosilicates have preserved in their hydrogen isotope values evidence for the interaction between rock and meteoric fluids during mylonite development. Fluid flow generates mass transfer along the tectonic anisotropy, and mylonites do not undergo significant volume change, except locally in ultramylonite zones. Hydration of detachment shear zones attends mechanical grain size reduction and enhances strain softening and localization. Self-exhuming detachment shear zones evolve rapidly (a few million years) through the transition from ductile to brittle, which is partly controlled by the thermal effect of circulating surface fluids. Detachment systems are zones in the crust where strain and fluid flow are coupled; these systems. evolve rapidly toward strain localization and therefore efficient exhumation.

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.

Stability and asymptotic analysis of a fluid-particle interaction model

We are interested in coupled microscopic/macroscopic models describing the evolution of particles dispersed in a fluid. The system consists in a Vlasov-Fokker-Planck equation to describe the microscopic motion of the particles coupled to the Euler equations for a compressible fluid. We investigate dissipative quantities, equilibria and their stability properties and the role of external forces. We also study some asymptotic problems, their equilibria and stability and the derivation of macroscopic two-phase models.

Remediation of a marine shore tailings deposit and the importance of water-rock interaction on element cycling in the coastal aquifer

We present the study of the geochemical processes associated with the first successful remediation of a marine shore tailings deposit in a coastal desert environment (Bahia de Ite, in the Atacama Desert of Peru). The remediation approach implemented a wetland on top of the oxidized tailings. The site is characterized by a high hydrauliz gradient produced by agricultural irrigation on upstream gravel terraces that pushed river water (similar to 500 mg/L SO(4)) toward the sea and through the tailings deposit. The geochemical and isotopic (delta(2)H(water) and delta(18)O(water), delta(34)S(sulfate) , delta(18)O(sulfate)) approach applied here revealed that evaporite horizons (anhydrite and halite) in the gravel terraces are the source of increased concentrations of SO(4), Cl, and Na up to similar to 1500 mg/L in the springs at the base of the gravel terraces. Deeper groundwater interacting with underlying marine sequences increased the concentrations of SO(4), Cl, and Na up to 6000 mg/L and increased the alkalinity up to 923 mg/L CaCO(3) eq. in the coastal aquifer. These waters infiltrated into the tailings deposit at the shelf-tailings interface. Nonremediated tailings had a low-pH oxidation zone (pH 1-4) with significant accumulations of efflorescent salts (10-20 cm thick) at the surface because of upward capillary transport of metal cations in the arid climate. Remediated tailings were characterized by neutral pH and reducing conditions (pH similar to 7, Eh similar to 100 mV). As a result, most bivalent metals such as Cu, Zn, and Ni had very low concentrations (around 0.01 mg/L or below detection limit) because of reduction and sorption processes. In contrast, these reducing conditions increased the mobility of iron from two sources in this system: (1) The originally Fe(III)-rich oxidation zone, where Fe(II) was reduced during the remediation process and formed an Fe(II) plume, and (2) reductive dissolution of Fe(III) oxides present in the original shelf lithology formed an Fe-Mn plume at 10-m depth. These two Fe-rich plumes were pushed toward the shoreline where more oxidizing and higher pH conditions triggered the precipitation of Fe(HI)hydroxide coatings on silicates. These coatings acted as a filter for the arsenic, which naturally infiltrated with the river water (similar to 500 mu g/L As natural background) into the tailings deposit.

Stable isotope (C, O, S) systematics of the mercury mineralization at Idrija, Slovenia: constraints on fluid source and alteration processes

The world-class Idrija mercury deposit (western Slovenia) is hosted by highly deformed Permocarboniferous to Middle Triassic sedimentary rocks within a complex tectonic structure at the transition between the External Dinarides and the Southern Alps. Concordant and discordant mineralization formed concomitant with Middle Triassic bimodal volcanism in an aborted rift. A multiple isotopic (C, O, S) investigation of host rocks and ore minerals was performed to put constraints on the source and composition of the fluid, and the hydrothermal alteration. The distributions of the delta(13)C and delta(18)O values of host and gangue carbonates are indicative of a fracture-controlled hydrothermal system, with locally high fluid-rock ratios. Quantitative modeling of the delta(13)C and delta(18)O covariation for host carbonates during temperature dependent fluid-rock interaction, and concomitant precipitation of void-filling dolomites points to a slightly acidic hydrothermal fluid (delta(13)Capproximate to-4parts per thousand and delta(18)Oapproximate to+10parts per thousand), which most likely evolved during isotopic exchange with carbonates under low fluid/rock ratios. The delta(34)S values of hydrothermal and sedimentary sulfur minerals were used to re-evaluate the previously proposed magmatic and evaporitic sulfur sources for the mineralization, and to assess the importance of other possible sulfur sources such as the contemporaneous seawater sulfate, sedimentary pyrite, and organic sulfur compounds. The delta(34)S values of the sulfides show a large variation at deposit down to hand-specimen scale. They range for cinnabar and pyrite from -19.1 to +22.8parts per thousand, and from -22.4 to +59.6parts per thousand, respectively, suggesting mixing of sulfur from different sources. The peak of delta(34)S values of cinnabar and pyrite close to 0parts per thousand is compatible with ore sulfur derived dominantly from a magmatic fluid and/or from hydrothermal leaching of basement rocks. The similar stratigraphic trends of the delta(34)S values of both cinnabar and pyrite suggest a minor contribution of sedimentary sulfur (pyrite and organic sulfur) to the ore formation. Some of the positive delta(34)S values are probably derived from thermochemical reduction of evaporitic and contemporaneous seawater sulfates.

Modeling of the Fluid-Structure Interaction in Inelastic Piping Systems

In many engineering applications, compliant piping systems conveying liquids are subjected to inelastic deformations due to severe pressure surges such as plastic tubes in modern water supply transmission lines and metallic pipings in nuclear power plants. In these cases the design of such systems may require an adequate modeling of the interactions between the fluid dynamics and the inelastic structural pipe motions. The reliability of the prediction of fluid-pipe behavior depends mainly on the adequacy of the constitutive equations employed in the analysis. In this paper it is proposed a systematic and general approach to consistently incorporate different kinds of inelastic behaviors of the pipe material in a fluid-structure interaction analysis. The main feature of the constitutive equations considered in this work is that a very simple numerical technique can be used for solving the coupled equations describing the dynamics of the fluid and pipe wall. Numerical examples concerning the analysis of polyethylene and stainless steel pipe networks are presented to illustrate the versatility of the proposed approach.

Genesis and Evolution of Brittle Structures In Southwestern Finland and Western South Africa. Insights into Fault Reactivation, Fluid Flow and Structural Maturity in Precambrian Cratons

The bedrock of old crystalline cratons is characteristically saturated with brittle structures formed during successive superimposed episodes of deformation and under varying stress regimes. As a result, the crust effectively deforms through the reactivation of pre-existing structures rather than by through the activation, or generation, of new ones, and is said to be in a state of 'structural maturity'. By combining data from Olkiluoto Island, southwestern Finland, which has been investigated as the potential site of a deep geological repository for high-level nuclear waste, with observations from southern Sweden, it can be concluded that the southern part of the Svecofennian shield had already attained structural maturity during the Mesoproterozoic era. This indicates that the phase of activation of the crust, i.e. the time interval during which new fractures were generated, was brief in comparison to the subsequent reactivation phase. Structural maturity of the bedrock was also attained relatively rapidly in Namaqualand, western South Africa, after the formation of first brittle structures during Neoproterozoic time. Subsequent brittle deformation in Namaqualand was controlled by the reactivation of pre-existing strike-slip faults.In such settings, seismic events are likely to occur through reactivation of pre-existing zones that are favourably oriented with respect to prevailing stresses. In Namaqualand, this is shown for present day seismicity by slip tendency analysis, and at Olkiluoto, for a Neoproterozoic earthquake reactivating a Mesoproterozoic fault. By combining detailed field observations with the results of paleostress inversions and relative and absolute time constraints, seven distinctm superimposed paleostress regimes have been recognized in the Olkiluoto region. From oldest to youngest these are: (1) NW-SE to NNW-SSE transpression, which prevailed soon after 1.75 Ga, when the crust had sufficiently cooled down to allow brittle deformation to occur. During this phase conjugate NNW-SSE and NE-SW striking strike-slip faults were active simultaneous with reactivation of SE-dipping low-angle shear zones and foliation planes. This was followed by (2) N-S to NE-SW transpression, which caused partial reactivation of structures formed in the first event; (3) NW-SE extension during the Gothian orogeny and at the time of rapakivi magmatism and intrusion of diabase dikes; (4) NE-SW transtension that occurred between 1.60 and 1.30 Ga and which also formed the NW-SE-trending Satakunta graben located some 20 km north of Olkiluoto. Greisen-type veins also formed during this phase. (5) NE-SW compression that postdates both the formation of the 1.56 Ga rapakivi granites and 1.27 Ga olivine diabases of the region; (6) E-W transpression during the early stages of the Mesoproterozoic Sveconorwegian orogeny and which also predated (7) almost coaxial E-W extension attributed to the collapse of the Sveconorwegian orogeny. The kinematic analysis of fracture systems in crystalline bedrock also provides a robust framework for evaluating fluid-rock interaction in the brittle regime; this is essential in assessment of bedrock integrity for numerous geo-engineering applications, including groundwater management, transient or permanent CO2 storage and site investigations for permanent waste disposal. Investigations at Olkiluoto revealed that fluid flow along fractures is coupled with low normal tractions due to in-situ stresses and thus deviates from the generally accepted critically stressed fracture concept, where fluid flow is concentrated on fractures on the verge of failure. The difference is linked to the shallow conditions of Olkiluoto - due to the low differential stresses inherent at shallow depths, fracture activation and fluid flow is controlled by dilation due to low normal tractions. At deeper settings, however, fluid flow is controlled by fracture criticality caused by large differential stress, which drives shear deformation instead of dilation.

The Fluid Structure Interaction Analysis of a Peristaltic Pump

The thesis work models the squeezing of the tube and computes the fluid motion of a peristaltic pump. The simulations have been conducted by using COMSOL Multiphysics FSI module. The model is setup in axis symmetric with several simulation cases to have a clear understanding of the results. The model captures total displacement of the tube, velocity magnitude, and average pressure fluctuation of the fluid motion. A clear understanding and review of many mathematical and physical concepts are also discussed with their applications in real field. In order to solve the problems and work around the resource constraints, a thorough understanding of mass balance and momentum equations, finite element concepts, arbitrary Lagrangian-Eulerian method, one-way coupling method, two-way coupling method, and COMSOL Multiphysics simulation setup are understood and briefly narrated.

The mineralizing fluid in the Piaba gold deposit, São Luís cratonic fragment (NW-Maranhão, Brazil) based on fluid inclusion studies in quartz veins

O depósito aurífero de Piaba tornou-se a primeira mina em operação no fragmento cratônico São Luís, noroeste do Maranhão. Seu ambiente geológico compreende rochas metavulcanossedimentares do Grupo Aurizona e granitoides da Suíte Tromaí, entre outras unidades menores, formadas em ambiente de arcos de ilhas entre 2240 e 2150 Ma, juntamente com outras unidades menores. A mineralização é caracterizada por uma trama stockwork de veios e vênulas de quartzo com seus halos de alteração (clorita + muscovita + carbonato + pirita + calcopirita e ouro) hospedada em um granodiorito granofírico fino (Granófiro Piaba) e em rocha subvulcânica andesítica do Grupo Aurizona. O corpo mineralizado é espacialmente limitado à zona de cisalhamento ENE-WSW rúptil-dúctil (Falha Piaba). Estudos petrográficos, microtermométricos e por espectroscopia microRaman no quartzo definiram inclusões aquo-carbônicas bifásicas e trifásicas, produzidas por aprisionamento heterogêneo durante separação de fases, e fluidos aquosos tardios. A solução mineralizadora corresponde a um fluido aquo-carbônico composto por CO₂ (5 - 24 mol%, densidade de 0,96 - 0,99 g/cm³), H₂O (74 - 93 mol%), N₂ (< 1 mol%), CH₄ (<1mol%) e 5,5 % em peso NaCl equivalente. O minério depositou a 267 - 302ºC e 1,25 - 2,08 kbar, correspondendo a profundidades de 4 a 7 km, em consonância com o regime estrutural. A composição e o intervalo de P-T do fluido mineralizador, combinados com o caráter redutor (log ƒ_O2 -31,3 a -34,3) e a sulfetação da rocha hospedeira, sugerem que o ouro foi transportado como um complexo sulfetado. O minério foi depositado em consequência da separação de fase, redução da atividade de enxofre e da ƒO₂ pela interação fluido-rocha.