Hydrogen isotope fractionations between amphiboles, micas, and fluids in alkaline igneous intrusions

RÉSUMÉ DE LA THÈSE Les teneurs des amphiboles en éléments majeurs et en isotopes stables ont été analysées dans plusieurs complexes ignés alcalins et hyperalcalins, dans le but de déterminer l'importance des variations de composition des minéraux pour le fractionnement isotopique de l'hydrogène dans un système naturel minéral-magma-fluide. Cette étude se concentre principalement sur les syénites néphéliniques de complexes intrusifs alcalins bien connus mais à chimie variable, dont les amphiboles, ainsi que d'autres silicates hydratés tels que micas et eudialytes, lorsque cela était possible, ont été séparés. L'intérêt principal s'est porté sur le complexe alcalin d'Ilímaussaq de la Province du Gardar, au Sud du Groenland. Dans une optique de comparaison, nous avons collecté et analysé d'autres échantillons provenant du complexe de Tugtutôq (Sud Groenland), des complexes de Khibina et Lovozero (Péninsule de Kola, Russie), du Mont St-Hilaire et du Mont Royal (Canada) et de 6 autres du nord-ouest de la Namibie (Cape Cross, Okenyenya, Messum, Etaneno, Kalkfeld,et Okorusu). Les compositions isotopiques de l'hydrogène des amphiboles des ces différentes zones présentent de grandes variations (-227 à -700/00), ce qui est atypique pour des magmas d'origine mantellique. Les valeurs comprises entre -80 et -400/00 indiquent une provenance du manteau. Ces larges variations de compositions ainsi que l'extrême appauvrissement en isotope lourd de l'hydrogène (D), en comparaison avec d'autres roches ignées, semblent être propres.aux roches alcalines et hyperalcalines de ce type, ce qui indiquerait un processus commun. Les différents complexes alcalins choisis présentent un large intervalle de composition chimique des amphiboles. La caractérisation des amphiboles par microscopie électronique et par spectroscopie Mössbauer contribuent à observer le contrôle du Fe sur le fractionnement des isotopes de l'hydrogène. En effet, cela a mis en évidence un contrôle du Fe sur le fractionnement et même, dans le cas du complexe hyperalcalin d'Ilímaussaq, une relation entre le rapport Fei3+/FeT et les variations du rapport D/H. Les complexes étudiés diffèrent de par leur index agpaïtique (Na+K/Al) et également de par leur contenu en fer. Les plus hautes valeurs en Fe (27-35 wt%) et en éléments alcalins dans les amphiboles, ainsi que les teneurs de D/H les plus basses et leur grande variation, sont celles du complexe d'Ilímaussaq. Les amphiboles de la Péninsule de Kola et du Canada sont similaires, mais toutefois moins appauvries en D. En ce qui concerne les amphiboles des complexes du NO de la Namibie, elles présentent des compositions isotopiques de l'hydrogène magmatiques normales (-73 à -100 0/00), contiennent moins de Fe (15-17 wt%) et sont fortement enrichies en Ca et moins en Na. Dans ce cas, l'alcalinité est moins importante en comparaison des autres complexes étudiés. En dehors des teneurs en éléments alcalins des amphiboles, l'alcalinité des fluides s'avère également un facteur important, ce qui est cohérent avec certaines suggestions à partir de systèmes expérimentaux. Afin de mieux contraindre ce facteur, des expériences d'échanges hydrothermaux entre les amphiboles et les fluides de salinité différente ont été effectuées en simulant des conditions naturelles. L'approximation d'amphiboles naturelles de complexes ignés alcalins, couplée aux expériences d'échange, aide à préciser les facteurs contrôlant le fractionnement des isotopes de l'hydrogène dans les roches alcalines. Les valeurs extrêmement basses de 3D des amphiboles de ces complexes alcalins peuvent être dues à une combinaison de différents facteurs, telles qu'une haute alcalinité, une haute teneur en Fe et une faible profondeur d'intrusion. Les grandes variations ainsi que les faibles valeurs de SD des amphiboles étudiées peuvent résulter d'un processus magmatique interne et il est peu probable que de l'eau météorique soit impliquée et/ou que le dégazage magmatique ait joué un rôle. THESIS ABSTRACT Major element and stable isotope compositions of amphiboles were analyzed from a number of alkaline and peralkaline igneous complexes in order to determine the importance of compositional variations in minerals to hydrogen isotope fractionations in natural mineral-melt-fluid systems. The thesis mainly focuses on nepheline syenites of well-studied, but chemically variable alkaline intrusive rocks, from which amphiboles and, if possible, other hydrous silicates such as micas and eudialytes were separated. The system of primary interest was the alkaline Ilímaussaq Complex of the Gardar Province of South Greenland. For the purpose of comparison additional samples were collected and examined from the Tugtutôq Complex (South Greenland), the Khibina and Lovozero Complexes (Kola Peninsular, Russia), Mount St-Hilaire and Mount Royal (Canada) and six further complexes from NW Namibia (Cape Cross, Okenyenya, Messum, Etaneno, Kalkfeld, and Okorusu). The hydrogen isotope compositions of amphiboles from the localities studied differ greatly, which is atypical for amphiboles from mantle, range between - 227 and - 700/00 (latter compatible with a simple mantle origin). As this wide range in compositions and the extreme depletion in the heavy hydrogen isotope (D) content relative to other igneous rocks appear to be unique to alkaline to peralkaline rocks of this type, a common process is indicated. The different alkaline complexes chosen cover a wide range of amphibole chemical compositions. Detailed chemical characterization of amphiboles by electron microprobe and Mössbauer spectroscopy analyses helped to constrain the control of Fe on the H-isotope fractionations. Complete characterization of the chemical compositions of the amphiboles support Fe-control on fractionations and at least for the peralkaline Ilímaussaq complex a relationship between Fe3+/FeT ratios and variations in D/H. The studied complexes differ in their agpaitic index (Na+K/Al) and also in their Fe-content. The most iron (27-35 wt. %) and alkaline element rich amphiboles, with the lowermost D/H ratio, as well with very wide range, are the ones from Ilímaussaq complex. Similar, but less D depleted amphiboles are from the Kola Peninsula and the Canadian localities. The complexes described from NW Namibia have amphiboles with normal magmatic hydrogen isotope composition (-730/00 to -1000/00), and have less Fe-content (15-17 wt. %), and are more Ca-and less Na-rich. In this case alkalinity is not that important in comparison to the other studied complexes. Beside the alkaline element contents in the amphiboles, the alkalinity of the fluids has been found to be an important factor, in conjunction with earlier suggestions from experimental systems. To further constrain this factor, hydrothermal exchange experiments between amphiboles and fluids of different salinity simulating natural conditions were performed. The approach of examining natural amphiboles from alkaline igneous complexes in parallel to performing exchange experiments - helped to further constrain the factors controlling the H-isotope fractionations in alkaline rocks. The observed changes between the hydrogen and oxygen isotope compositions of amphiboles and fluids before and after the experiments suggest that another phase was produced during the experiments, which influenced the final hydrogen isotope composition of the system. This presumably hydrous phase has also influenced the Fe3 +/Fe2+ ratio of the amphiboles, which became more oxidized. The extremely low SD values of amphiboles in these alkaline complexes may be due to a combination of different factors such as high alkalinity, high Fe-content, and shallow intrusion depths. This wide range and the low SD values of the amphiboles studied might be a result of internal, magmatic processes and it is unlikely that meteoric water was involved and/or magmatic degassing played an important role. RÉSUMÉ DE LA THÈSE (pour le grand public) Fractionnement isotopique de l'hydrogène entre amphiboles, micas et fluides dans des intrusions alcalines Zsófia Wáczek Directeur de thèse, Prof. Torsten W. Vennemann Institut de Minéralogie et Géochimie, Université de Lausanne Les roches alcalines et celles qui leurs sont associées sont des sources importantes de nombreux minéraux et minerais, tels l'apatite, le niobium, le diamant et autres pierres précieuses. Cette étude se concentre sur des complexes alcalins localisés dans le sud du Groenland, au Canada, dans la péninsule de Kola en Russie et au nord-ouest de la Namibie. Ces complexes sont composés de roches ayant cristallisé à partir de magmas et de fluides très enrichis en alcalins. Cet enrichissement permet la précipitation de minéraux inhabituels riches en potassium et/ou sodium, telles les amphiboles sodiques, également enrichies en fer. Les amphiboles étudiées ont des compositions calciques, sodi-calciques et sodiques, qui reflètent leurs différents environnements de formation. Des études précédentes ont révélé une large gamme de rapports isotopiques de l'hydrogène dans les amphiboles de roches hyperalcalines, dont certains extrêmement bas. Cette variation importante est très intrigante, sachant que des valeurs entre -40 et -800/00 correspondent à des silicates ignés hydratés et non altérés, alors que des valeurs descendant jusqu'a -1500/00 nécessiteraient une altération par de l'eau météorique et/ou une contamination par les roches environnantes ou des sédiments riches en matière organique. Dans lé cas précis du complexe d'Ilímaussaq (sud du Groenland), aucune de ces explications n'a pu être démontrée et des valeurs encore plus faibles ont été trouvées. Le complexe d'Ilímaussaq présente des valeurs de rapport isotopique de l'hydrogène entre -227 et -500/00 dans les amphiboles. Une origine mantellique permet d'expliquer les valeurs élevées, mais d'autres processus doivent entrer en jeu pour engendrer les valeurs les plus négatives. C'est à l'identification de ces processus que nous nous sommes attachés dans ce travail. Les grandes variations observées dans les teneurs en fer et dans le rapport Fe3+/FeT des roches et des minéraux de ces complexes sont corrélées avec d'autres paramètres chimiques, tels que la composition isotopique de l'hydrogène dans les amphiboles. Nous avons dès lors abordé les questions suivantes: quelle est la relation entre la teneur en fer des amphiboles et leur composition isotopique? Que nous apprennent les changements de la teneur en fer et les changements dans le rapport Fe3+/FeT sur les processus pétrologiques dans ces roches? Pour répondre à ces questions, nous avons analysé les compositions isotopiques de l'oxygène et de l'hydrogène dans les amphiboles et d'autres silicates hydratés. La composition chimique et le rapport Fe3+/FeT des amphiboles ont également été déterminés. Des expériences hydrothermales simulant des conditions naturelles ont été entreprises afin de mieux comprendre les processus de fractionnement isotopiques dans ces systèmes très alcalins. Nos conclusions sont les suivantes: (1) Les valeurs extrêmement faibles ainsi que les larges variations des rapports isotopiques de l'hydrogène des amphiboles de ces complexes alcalins sont dues à une combinaison de facteurs tels que la forte alcalinité, la haute teneur en fer et la profondeur très faible de l'intrusion. (2) Ces valeurs sont probablement le résultat de processus magmatiques internes. (3) Il est peu probable que les eaux météoriques et/ou le dégazage magmatique aient joué un rôle lors de la formation de ces amphiboles. (4) Certaines corrélations, en accord avec les études précédentes, ont pu être trouvées au niveau des concentrations en fer. (5) Dans le cas du complexe d'Ilímaussaq exclusivement, une relation a été trouvée entre le rapport Fe3+/FeT et la composition isotopique de l'hydrogène des amphiboles.

ADVANCED NUMERICAL METHODS FOR SIMULATING MICROBUBBLE IN AN AERATED MIXING TANK

Fluid mixing in mechanically agitated tanks is one of the major unit operations in many industries. Bubbly flows have been of interest among researchers in physics, medicine, chemistry and technology over the centuries. The aim of this thesis is to use advanced numerical methods for simulating microbubble in an aerated mixing tank. Main components of the mixing tank are a cylindrical vessel, a rotating Rushton turbine and the air nozzle. The objective of Computational Fluid Dynamics (CFD) is to predict fluid flow, heat transfer, mass transfer and chemical reactions. The CFD simulations of a turbulent bubbly flow are carried out in a cylindrical mixing tank using large eddy simulation (LES) and volume of fluid (VOF) method. The Rushton turbine induced flow is modeled by using a sliding mesh method. Numerical results are used to describe the bubbly flows in highly complex liquid flow. Some of the experimental works related to turbulent bubbly flow in a mixing tank are briefly reported. Numerical simulations are needed to complete and interpret the results of the experimental work. Information given by numerical simulations has a major role in designing and scaling-up mixing tanks. The results of this work have been reported in the following scientific articles: ·Honkanen M., Koohestany A., Hatunen T., Saarenrinne P., Zamankhan P., Large eddy simulations and PIV experiments of a two-phase air-water mixer, in Proceedings of ASME Fluids Engineering Summer Conference (2005). ·Honkanen M., Koohestany A., Hatunen T., Saarenrinne P., Zamankhan P., Dynamical States of Bubbling in an Aerated Stirring Tank, submitted to J. Computational Physics.

Numerical and experimental investigations of three-dimensional container filling with Newtonian viscous fluids

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

(Table 2) Composition of the deapest pore fluids from ODP Holes 168-1030B and 168-1031A

On the eastern flank of the Juan de Fuca Ridge, reaction between upwelling basement fluid and sediment alters hydrothermal fluxes of Ca, SiO2(aq), SO4, PO4, NH4, and alkalinity. We used the Global Implicit Multicomponent Reactive Transport (GIMRT) code to model the processes occurring in the sediment column (diagenesis, sediment burial, fluid advection, and multicomponent diffusion) and to estimate net seafloor fluxes of solutes. Within the sediment section, the reactions controlling the concentrations of the solutes listed above are organic matter degradation via SO4 reduction, dissolution of amorphous silica, reductive dissolution of amorphous Fe(III)-(hydr)oxide, and precipitation of calcite, carbonate fluorapatite, and amorphous Fe(II)-sulfide. Rates of specific discharge estimated from pore-water Mg profiles are 2 to 3 mm/yr. At this site the basement hydrothermal system is a source of NH4, SiO2(aq), and Ca, and a sink of SO4, PO4, and alkalinity. Reaction within the sediment column increases the hydrothermal sources of NH4 and SiO2(aq), increases the hydrothermal sinks of SO4 and PO4, and decreases the hydrothermal source of Ca. Reaction within the sediment column has a spatially variable effect on the hydrothermal flux of alkalinity. Because the model we used was capable of simulating the observed pore-water chemistry by using mechanistic descriptions of the biogeochemical processes occurring in the sediment column, it could be used to examine the physical controls on hydrothermal fluxes of solutes in this setting. Two series of simulations in which we varied fluid flow rate (1 to 100 mm/yr) and sediment thickness (10 to 100 m) predict that given the reactions modeled in this study, the sediment section will contribute most significantly to fluxes of SO4 and NH4 at slow flow rates and intermediate sediment thickness and to fluxes of SiO2(aq) at slow flow rates and large sediment thickness. Reaction within the sediment section could approximately double the hydrothermal sink of PO4 over a range of flow rates and sediment thickness, and could slightly decrease (by

Analysis Of The Fluorescence Of Body Fluids On Different Surfaces And Times.

The use of screening techniques, such as an alternative light source (ALS), is important for finding biological evidence at a crime scene. The objective of this study was to evaluate whether biological fluid (blood, semen, saliva, and urine) deposited on different surfaces changes as a function of the age of the sample. Stains were illuminated with a Megamaxx™ ALS System and photographed with a Canon EOS Utility™ camera. Adobe Photoshop™ was utilized to prepare photographs for analysis, and then ImageJ™ was used to record the brightness values of pixels in the images. Data were submitted to analysis of variance using a generalized linear mixed model with two fixed effects (surface and fluid). Time was treated as a random effect (through repeated measures) with a first-order autoregressive covariance structure. Means of significant effects were compared by the Tukey test. The fluorescence of the analyzed biological material varied depending on the age of the sample. Fluorescence was lower when the samples were moist. Fluorescence remained constant when the sample was dry, up to the maximum period analyzed (60 days), independent of the substrate on which the fluid was deposited, showing the novelty of this study. Therefore, the forensic expert can detect biological fluids at the crime scene using an ALS even several days after a crime has occurred.

Comparison Between Automated And Microscopic Analysis In Body Fluids Cytology.

37

Verification of a mixed high-order accurate DNS code for laminar turbulent transition by the method of manufactured solutions

This paper presents results on a verification test of a Direct Numerical Simulation code of mixed high-order of accuracy using the method of manufactured solutions (MMS). This test is based on the formulation of an analytical solution for the Navier-Stokes equations modified by the addition of a source term. The present numerical code was aimed at simulating the temporal evolution of instability waves in a plane Poiseuille flow. The governing equations were solved in a vorticity-velocity formulation for a two-dimensional incompressible flow. The code employed two different numerical schemes. One used mixed high-order compact and non-compact finite-differences from fourth-order to sixth-order of accuracy. The other scheme used spectral methods instead of finite-difference methods for the streamwise direction, which was periodic. In the present test, particular attention was paid to the boundary conditions of the physical problem of interest. Indeed, the verification procedure using MMS can be more demanding than the often used comparison with Linear Stability Theory. That is particularly because in the latter test no attention is paid to the nonlinear terms. For the present verification test, it was possible to manufacture an analytical solution that reproduced some aspects of an instability wave in a nonlinear stage. Although the results of the verification by MMS for this mixed-order numerical scheme had to be interpreted with care, the test was very useful as it gave confidence that the code was free of programming errors. Copyright (C) 2009 John Wiley & Sons, Ltd.

DEVELOPMENT OF STATIC MIXERS FOR MISCIBLE FLUIDS IN LAMINAR FLOW WITH THE USE OF COMPUTATIONAL FLUID DYNAMICS (CFD)

Static mixers with improved performance were developed from CFD simulations in a stepwise approach. The relevant geometric features of simple mixer designs and the corresponding mixing mechanisms-laminar shear, elongational flow, and distributive mixing-were identified first. This information was used to formulate guidelines for the development of new geometries. The solid elements of the static mixer should: (a) provide restrictions to the flow; (b) deflect the flow; (c) be sequentially rotated around the flow direction to provide symmetry; (d) extend from the center of the pipe to the vicinity of the walls to avoid short-circuiting; and (e) distribute and remix the flow. Based on these guidelines, two improved mixer designs were developed: the DS A-I mixer has a good mixing efficiency and an acceptable pressure drop; the Fins 35 degrees mixer is more efficient and compact, but requires a larger pressure drop. Their performance indicates that their use is possible on industrial applications.

Theory of modulational instability in Bragg gratings with quadratic nonlinearity

Modulational instability in optical Bragg gratings with a quadratic nonlinearity is studied. The electric field in such structures consists of forward and backward propagating components at the fundamental frequency and its second harmonic. Analytic continuous wave (CW) solutions are obtained, and the intricate complexity of their stability, due to the large number of equations and number of free parameters, is revealed. The stability boundaries are rich in structures and often cannot be described by a simple relationship. In most cases, the CW solutions are unstable. However, stable regions are found in the nonlinear Schrodinger equation limit, and also when the grating strength for the second harmonic is stronger than that of the first harmonic. Stable CW solutions usually require a low intensity. The analysis is confirmed by directly simulating the governing equations. The stable regions found have possible applications in second-harmonic generation and dark solitons, while the unstable regions maybe useful in the generation of ultrafast pulse trains at relatively low intensities. [S1063-651X(99)03005-6].

Which one of these is not like the others? An inter-hospital study of the viscosity of thickened fluids

Engineering This investigation examined the rheological (viscosity and yield stress) and material property (density) characteristics of the thickened meal-time and videofluorscopy fluids provided by 10 major metropolitan hospitals. Differences in the thickness of thickened fluids were considered as a source of variability and potential hazard for inter-hospital transfers of dysphagic patients. The results indicated considerable differences in the viscosity, density, and yield stress of both meal-time and videofluoroscopy fluids. In theory, the results suggest that dysphagic patients transferred between hospitals could be placed on inappropriate levels of fluid thickness because of inherent differences in the rheology and material property characteristics of the fluids provided by different hospitals. Slowed improvement or medical complications are potential worst-case scenarios for dysphagic patients if the difference between the thick fluids offered by 2 hospitals are extreme. The investigation outlines the most appropriate way to assess the rheological and material property characteristics of thickened fluids. In addition, it suggests a plan of quality improvement to reduce the variability of the thickness of fluids offered at different hospitals.

How thick is thick? Multicenter study of the rheological and material property characteristics of mealtime fluids and videofluoroscopy fluids

Objective rheological assessment of fluids given to dysphagic patients at mealtime and during videofluoroscopy was carried out using a multicenter format. Thin, quarter-thick, half-thick and full-thick fluids were examined for the degree of correlation between mealtime fluids and their allegedly matched videofluoroscopy counterparts. The study was carried out to determine whether perceived subjective differences between mealtime fluids and videofluoroscopy fluids could be quantified using the rheological parameters of viscosity, density, and yield stress. The results showed poor correlation between mealtime fluids and videofluoroscopy fluids over all parameters. In general, the videofluoroscopy fluids were more viscous, more dense, and showed higher yield stress values than their mealtime counterparts. Given these results, it is reasonable to assume that the fluids used during videofluoroscopy do not provide an accurate indication of swallowing ability at mealtime. Therefore, it is suggested that clinicians use objective methods to theologically match videofluoroscopy fluids to mealtime fluids.

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.

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.

A 2D numerical model for simulating the physics of fault systems

Simulations provide a powerful means to help gain the understanding of crustal fault system physics required to progress towards the goal of earthquake forecasting. Cellular Automata are efficient enough to probe system dynamics but their simplifications render interpretations questionable. In contrast, sophisticated elasto-dynamic models yield more convincing results but are too computationally demanding to explore phase space. To help bridge this gap, we develop a simple 2D elastodynamic model of parallel fault systems. The model is discretised onto a triangular lattice and faults are specified as split nodes along horizontal rows in the lattice. A simple numerical approach is presented for calculating the forces at medium and split nodes such that general nonlinear frictional constitutive relations can be modeled along faults. Single and multi-fault simulation examples are presented using a nonlinear frictional relation that is slip and slip-rate dependent in order to illustrate the model.