Investigating Microbial Trace-fossils and Abiotic Alteration in Hydrovolcanic Tuffs of the Fort Rock Volcanic Field, Oregon

Microbial ichnofossils in volcanic rocks provide a significant record of subsurface microbes and potentially extraterrestrial biosignatures. Here, the textures, mineralogy, and geochemistry of two continental basaltic hydrovolcanic deposits - Reed Rocks and Black Hills - in the Fort Rock Volcanic Field (FRVF) are investigated. Methods include petrographic microscopy, micro and powder X-ray diffraction, SEM/BSE/EDF imaging, energy dispersive spectroscopy, stable isotopes, and X-ray fluorescence. Petrographic analysis revealed granular and tubular textures with biogenic morphologies that include terminal enlargements, septate divisions, branching forms, spiral filaments, and ovoid bodies resembling endolithic microborings described in ocean basalts. They display evidence of behaviour and a geologic context expressing their relative age and syngenicity. Differences in abiotic alteration and the abundance/morphotype assemblage of putative microborings between the sites indicate that water/rock ratio, fluid composition and flux, temperature and secondary phase formation are influences on microboring formation. This study is the first report of reputed endolithic microborings in basalts erupted in a continental lacustrine setting.

Excited state electronic polarization and reappraisal of the n <- pi* emission of acetone in water

Electronic polarization of the acetone molecule in the excited n -> pi* state is considered and its influence on the solvent shift in the emission spectrum is analyzed. Using an iterative procedure the electronic polarizations of both the ground and the excited states are included and compared with previous results obtained with Car-Parrinello dynamics. Analysis of the emission transition obtained using CIS(D)/aug-cc-pVDZ on statistically uncorrelated solute-solvent structures, composed of acetone and twelve explicit water molecules embedded in the electrostatic field of remaining 263 water molecules, corroborates that the solvent effect is mild, calculated here between 80 and 380 cm (1). (c) 2010 Published by Elsevier B.V.

Polarization and Spectral Shift of Benzophenone in Supercritical Water

Monte Carlo simulation and quantum mechanics calculations based on the INDO/CIS and TD-DFT methods were utilized to study the solvatochromic shift of benzophenone when changing the environment from normal water to supercritical (P = 340.2 atm and T = 673 K) condition. Solute polarization increases the dipole moment of benzophenone, compared to gas phase, by 88 and 35% in normal and supercritical conditions, giving the in-solvent dipole value of 5.8 and 4.2 D, respectively. The average number of solute-solvent hydrogen bonds was analyzed, and a large decrease of 2.3 in normal water to only 0.8 in the supercritical environment was found. By using these polarized models of benzophenone in the two different conditions of water, we performed MC simulations to generate statistically uncorrelated configurations of the solute surrounded by the solvent molecules and subsequent quantum mechanics calculations on these configurations. When changing from normal to supercritical water environment, INDO/CIS calculations explicitly considering all valence electrons of the 235 solvent water molecules resulted in a solvatochromic shift of 1425 cm(-1) for the most intense transition of benzophenone, that is, slightly underestimated in comparison with the experimentally inferred result of 1700 cm(-1). TD-B3LYP/6-311+G(2d,p) calculations on the same configurations but with benzophenone electrostatically embedded in the 320 water molecules resulted in a solvatochromic shift of 1715 cm(-1) for this transition, in very good agreement with the experimental result. When using the unpolarized model of the benzophenone, this calculated solvatochromic shift was only 640 cm(-1). Additional calculations were also made by using BHandHLYP/6-311+G(2d,p) to analyze the effect of the asymptotic decay of the exchange functional. This study indicates that, contrary to the general expectation, there is a sizable solute polarization even in the low-density regime of supercritical condition and that the inclusion of this polarization is important for a reliable description of the spectral shifts considered here.

Understanding water and ion transport behaviour and permeability through poly(amide) thin film composite membrane

Molecular dynamics (MD) together with the adaptive biasing force (ABF) and metadynamics free energy calculation methods was used to investigate the permeation properties of salt water through poly(amide) thin film composite reverse osmosis membranes. The thin films were generated by annealing an amorphous cell of poly(amide) chains through an MD method. The MD results showed they have typical structural properties of the active layer of thin film composite membranes and comparable water diffusivity (2.13×10^-5cm²/s for the film with a density of 1.06g/cm³) and permeability (9.27×10^-15cm³cm/cm²sPa) to experimental data. The simulations of water permeation through the films under different transmembrane pressures revealed the behaviours of water molecules in the thin films and the dynamic regimes of water permeation, including Brownian diffusion, flush and jump diffusion regimes. The intermolecular interactions of water and ions with poly(amide) chains showed a strong dependence on the local structure of films. The attraction between water and ploy(amide) molecules can be up to 8.5kcal/mol in dense polymer regions and 5kcal/mol in the pores of about 3nm. The ABF and metadynamics simulations produced the profiles of free energy potential of water and ions along the depth of the thin films, which provided important information for quantitatively determining the barrier energy required for water permeation and rejection of ions. The thin film with a density of 1.06g/cm³ and a thickness of 6nm offers a rejection to Na⁺ but a slight absorption of Cl^- (0.25kcal/mol) at 0.3-0.4nm distance to its surface. Water molecules must overcome 63kcal/mol energy to move to the centre of the film. The dependences of the barrier energy and the water-polymer interaction energy on the local free volume size in the thin film were analysed. The simulations of water permeation under high transmembrane pressures showed a nonlinear response of the concentration and distribution of water molecules in the film to the imposed pressure. Compaction of the film segments close to the porous substrate and water congestion in dense regions significantly influenced the water permeation when the membrane was operated under pressures of more than 3.0MPa.

Theoretical studies of liquids by computer simulations: the water-acetone mixture.

Monte Carlo simulation results for pure liquid acetone and water-acetone mixtures calculated in the isothermal and isobaric (NPT) ensemble at T=298K and p=1.0atm are presented. The TIP4P model was used for water and optimized potential for liquid simulation (OPLS) force field parameters used for acetone. The results obtained for the average configurational energy as a function of the mole fraction are in good accord with experimental data. Energy partitioning and co-ordination numbers results calculated for equimolar water-acetone solution are compared to similar data obtained for other water-organic liquid mixtures. These results show an increase in water-water interaction energy and co-ordination numbers when the interaction between water and organic liquid molecules decrease. Distribution functions for pure liquid acetone and water-acetone mixtures are presented. Dipole-dipole angular correlation functions obtained for pure liquid acetone show a predominance of dimers with parallel alignment of dipole moments. Radial distribution functions from water-acetone interaction show characteristic features of hydrogen bonded liquids. Radial and angular distribution functions for water-water correlation calculated in pure water and in equimolar water-acetone mixture are compared, showing very similar features in both systems. (C) 1999 Elsevier B.V. B.V. All rights reserved.

ONIOM study of dissociated hydrogen and water on ZnO surface

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

Interaction between wave and coastal structure: Validation of two lagrangian numerical models with experimental results marine 2011

Numerical modeling of the interaction among waves and coastal structures is a challenge due to the many nonlinear phenomena involved, such as, wave propagation, wave transformation with water depth, interaction among incident and reflected waves, run-up / run-down and wave overtopping. Numerical models based on Lagrangian formulation, like SPH (Smoothed Particle Hydrodynamics), allow simulating complex free surface flows. The validation of these numerical models is essential, but comparing numerical results with experimental data is not an easy task. In the present paper, two SPH numerical models, SPHysics LNEC and SPH UNESP, are validated comparing the numerical results of waves interacting with a vertical breakwater, with data obtained in physical model tests made in one of the LNEC's flume. To achieve this validation, the experimental set-up is determined to be compatible with the Characteristics of the numerical models. Therefore, the flume dimensions are exactly the same for numerical and physical model and incident wave characteristics are identical, which allows determining the accuracy of the numerical models, particularly regarding two complex phenomena: wave-breaking and impact loads on the breakwater. It is shown that partial renormalization, i.e. renormalization applied only for particles near the structure, seems to be a promising compromise and an original method that allows simultaneously propagating waves, without diffusion, and modeling accurately the pressure field near the structure.

Application of reactive fluid transport modeling to hydrothermal systems

A one-dimensional multi-component reactive fluid transport algorithm, 1DREACT (Steefel, 1993) was used to investigate different fluid-rock interaction systems. A major short coming of mass transport calculations which include mineral reactions is that solid solutions occurring in many minerals are not treated adequately. Since many thermodynamic models of solid solutions are highly non-linear, this can seriously impact on the stability and efficiency of the solution algorithms used. Phase petrology community saw itself faced with a similar predicament 10 years ago. To improve performance and reliability, phase equilibrium calculations have been using pseudo compounds. The same approach is used here in the first, using the complex plagioclase solid solution as an example. Thermodynamic properties of a varying number of intermediate plagioclase phases were calculated using ideal molecular, Al-avoidance, and non-ideal mixing models. These different mixing models can easily be incorporated into the simulations without modification of the transport code. Simulation results show that as few as nine intermediate compositions are sufficient to characterize the diffusional profile between albite and anorthite. Hence this approach is very efficient, and can be used with little effort. A subsequent chapter reports the results of reactive fluid transport modeling designed to constrain the hydrothermal alteration of Paleoproterozoic sediments of the Southern Lake Superior region. Field observations reveal that quartz-pyrophyllite (or kaolinite) bearing assemblages have been transformed into muscovite-pyrophyllite-diaspore bearing assemblages due to action of fluids migrating along permeable flow channels. Fluid-rock interaction modeling with an initial qtz-prl assemblage and a K-rich fluid simulates the formation of observed mineralogical transformation. The bulk composition of the system evolves from an SiO2-rich one to an Al2O3+K2O-rich one. Simulations show that the fluid flow was up-temperature (e.g. recharge) and that fluid was K-rich. Pseudo compound approach to include solid solutions in reactive transport models was tested in modeling hydrothermal alteration of Icelandic basalts. Solid solutions of chlorites, amphiboles and plagioclase were included as the secondary mineral phases. Saline and fresh water compositions of geothermal fluids were used to investigate the effect of salinity on alteration. Fluid-rock interaction simulations produce the observed mineral transformations. They show that roughly the same alteration minerals are formed due to reactions with both types of fluid which is in agreement with the field observations. A final application is directed towards the remediation of nitrate rich groundwaters. Removal of excess nitrate from groundwater by pyrite oxidation was modeled using the reactive fluid transport algorithm. Model results show that, when a pyrite-bearing, permeable zone is placed in the flow path, nitrate concentration in infiltrating water can be significantly lowered, in agreement with proposals from the literature. This is due to nitrogen reduction. Several simulations investigate the efficiency of systems with different mineral reactive surface areas, reactive barrier zone widths, and flow rates to identify the optimum setup.

Fluid-related inclusions in Alpine high-pressure peridotite reveal trace element recycling during subduction-zone dehydration of serpentinized mantle (Cima di Gagnone, Swiss Alps).

Serpentinites release at sub-arc depths volatiles and several fluid-mobile trace elements found in arc magmas. Constraining element uptake in these rocks and defining the trace element composition of fluids released upon serpentinite dehydration can improve our understanding of mass transfer across subduction zones and to volcanic arcs. The eclogite-facies garnet metaperidotite and chlorite harzburgite bodies embedded in paragneiss of the subduction melange from Cima di Gagnone derive from serpentinized peridotite protoliths and are unique examples of ultramafic rocks that experienced subduction metasomatism and devolatilization. In these rocks, metamorphic olivine and garnet trap polyphase inclusions representing the fluid released during high-pressure breakdown of antigorite and chlorite. Combining major element mapping and laser-ablation ICP-MS bulk inclusion analysis, we characterize the mineral content of polyphase inclusions and quantify the fluid composition. Silicates, Cl-bearing phases, sulphides, carbonates, and oxides document post-entrapment mineral growth in the inclusions starting immediately after fluid entrapment. Compositional data reveal the presence of two different fluid types. The first (type A) records a fluid prominently enriched in fluid-mobile elements, with Cl, Cs, Pb, As, Sb concentrations up to 10(3) PM (primitive mantle), similar to 10(2) PM Tit Ba, while Rb, B, Sr, Li, U concentrations are of the order of 10(1) PM, and alkalis are similar to 2 PM. The second fluid (type B) has considerably lower fluid-mobile element enrichments, but its enrichment patterns are comparable to type A fluid. Our data reveal multistage fluid uptake in these peridotite bodies, including selective element enrichment during seafloor alteration, followed by fluid-rock interaction along with subduction metamorphism in the plate interface melange. Here, infiltration of sediment-equilibrated fluid produced significant enrichment of the serpentinites in As, Sb, B, Pb, an enriched trace element pattern that was then transferred to the fluid released at greater depth upon serpentine dehydration (type A fluid). The type B fluid hosted by garnet may record the composition of the chlorite breakdown fluid released at even greater depth. The Gagnone study-case demonstrates that serpentinized peridotites acquire water and fluid-mobile elements during ocean floor hydration and through exchange with sediment-equilibrated fluids in the early subduction stages. Subsequent antigorite devolatilization at subarc depths delivers aqueous fluids to the mantle wedge that can be prominently enriched in sediment-derived components, potentially triggering arc magmatism without the need of concomitant dehydration/melting of metasediments or altered oceanic crust.

87Sr/86Sr ratios, Sr and Rb concentrations and stable oxygen isotope values for some basalts from ODP Hole 504B (Table 1)

Seventeen whole-rock samples, generally taken at 25-50 m intervals from 5 to 560 m sub-basement in Hole 504B, drilled in 6.2 m.y. old crust, were analysed for 87Sr/86Sr ratios, Sr and Rb concentrations, and 18O/16O ratios. Sr isotope ratios for 8 samples from the upper 260 m of the hole range from 0.70287 to 0.70377, with a mean of 0.70320. In the 330-560 m interval, 5 samples have a restricted range of 0.70255-0.70279, with a mean of 0.70266, the average value for fresh mid-ocean ridge basalts (MORB). In the 260-330 m interval, approximately intermediate Sr isotopic ratios are found. Delta18O values (?) range from 6.4 to 7.8 in the upper 260 m, 6.2-6.4 in the 270-320 m interval, and 5.8-6.2 in the 320-560 m interval. The values in the upper 260 m are typical for basalts which have undergone low-temperature seawater alteration, whereas the values for the 320-560 m interval correspond to MORB which have experienced essentially no oxygen isotopic alteration. The higher 87Sr/86Sr and 18O/16O ratios in the upper part of the hole can be interpreted as the result of a greater overall water/rock ratio in the upper part of the Hole 504B crust than in the lower part. Interaction of basalt with seawater (87Sr/86Sr = 0.7091) increased basalt 87Sr/86Sr ratios and produced smectitic alteration products which raised whole-rock delta18O values. Seawater circulation in the lower basalts may have been partly restricted by the greater number of relatively impermeable massive lava flows below about 230 m sub-basement. These flows may have helped to seal off lower basalts from through-flowing seawater.

Rb-Sr and K-Ar data for Upper Vendian clayey rocks of the Russian Platform

Geochemical behavior of Rb-Sr and K-Ar systems in Upper Vendian clayey rocks of the Russian Platform is under consideredation. The use of additional data on grain size fractions of sedimentary rocks recovered from boreholes drilled in the Gavrilov Yam area made it possible to confirm the previous conclusion on two stages of epigenetic matter transformation (approximately 600 and 400 Ma ago). Distortions are related to transformation of sediments due to interaction in the water-rock system. Interaction degree was more intense in the upper part of the sedimentary section relative to its lower strata. These conclusions are substantiated by materials from boreholes that characterize different types of Vendian sections and different tectonic zones.

Geochemistry of lavas and dikes from the upper oceanic crust of Hess Deep Rift, eastern Pacific Ocean

Strontium isotopes are useful tracers of fluid-rock interaction in marine hydrothermal systems and provide a potential way to quantify the amount of seawater that passes through these systems. We have determined the whole-rock Sr-isotopic compositions of a section of upper oceanic crust that formed at the fast-spreading East Pacific Rise, now exposed at Hess Deep. This dataset provides the first detailed comparison for the much-studied Ocean Drilling Program (ODP) drill core from Site 504B. Whole-rock and mineral Sr concentrations indicate that Sr-exchange between hydrothermal fluids and the oceanic crust is complex, being dependent on the mineralogical reactions occurring; in particular, epidote formation takes up Sr from the fluid increasing the 87Sr/86Sr of the bulk-rock. Calculating the fluid-flux required to shift the Sr-isotopic composition of the Hess Deep sheeted-dike complex, using the approach of Bickle and Teagle (1992, doi:10.1016/0012-821X(92)90221-G) gives a fluid-flux similar to that determined for ODP Hole 504B. This suggests that the level of isotopic exchange observed in these two regions is probably typical for modern oceanic crust. Unfortunately, uncertainties in the modeling approach do not allow us to determine a fluid-flux that is directly comparable to fluxes calculated by other methods.

Mineralogical composition, copper and d18O at DSDP Site 59-448

Alteration in a submarine remnant volcanic arc should leave an important record of (1) the mineralogy of sea water-volcanic arc rock interaction; (2) the chemistry of solid reaction products; (3) the isotopic characteristics of such reactions (Muehlenbachs and Clayton, 1972; Spooner, Beckinsale, et al , 1977; Spooner, Chapman, et al., 1977); (4) the metallogenesis within such a sequence (Mitchell and Bell, 1973); and (5) the geothermal gradient during the alteration. The volcaniclastic breccias, tuffs, and igneous units of Sites 448 (993 m) and 451 (930.5 m) on the Palau-Kyushu and West Mariana ridges, respectively, are particularly suited for such studies because the thick sequences have remained submarine throughout their history, seemingly unaffected by magmatic or hydrothermal events after cessation of volcanic activity. Also, shipboard observations indicated a change in alteration products with depth. At both sites the igneous units and volcaniclastic rocks were altered to brownish clays and zeolites near the top of the volcanic sequence; to bright blue green clays and zeolites at moderate depths; and to very dark, nearly opaque, forest green clays and zeolites at still greater depths. Native copper occurs both as disseminated pockets in the volcaniclastic breccias and vesicular basalts and as veins in the breccias; native copper is restricted to stratigraphic levels characterized by the absence of sulfides or oxides of copper and iron. Although some native copper is found in vesicles of basalts and may be orthomagmatic, most of it is clearly secondary. Near dikes and sills, higher sulfur fugacity conditions caused the precipitation of iron and copper sulfides with an absence of native copper (Garrels and Christ, 1965). The occurrence of native copper may be an initial stage of Cu metallogenesis that forms porphyry coppers in island arcs (Mitchell and Bell, 1973). This study will address primarily the possibility that hydrothermal sea water interaction with volcanic arc rocks has created the mineralogical and isotopic zonation in Leg 59 cores. Hydrothermal activity can be expected in a rapidly growing island arc and is probably the result of a high geothermal gradient prevalent during arc magmatic activity. The chemical character of the alteration is further discussed by Hajash (1981).

Stable oxygen isotope ratios in zircons from oceanic crust

Lower ocean crust is primarily gabbroic, although 1-2% felsic igneous rocks that are referred to collectively as plagiogranites occur locally. Recent experimental evidence suggests that plagiogranite magmas can form by hydrous partial melting of gabbro triggered by seawater-derived fluids, and thus they may indicate early, high-temperature hydrothermal fluid circulation. To explore seawater-rock interaction prior to and during the genesis of plagiogranite and other late-stage magmas, oxygen-isotope ratios preserved in igneous zircon have been measured by ion microprobe. A total of 197 zircons from 43 plagiogranite, evolved gabbro, and hydrothermally altered fault rock samples have been analyzed. Samples originate primarily from drill core acquired during Ocean Drilling Program and Integrated Ocean Drilling Program operations near the Mid-Atlantic and Southwest Indian Ridges. With the exception of rare, distinctively luminescent rims, all zircons from ocean crust record remarkably uniform d18O with an average value of 5.2 ± 0.5 per mil (2SD). The average d18O(Zrc) would be in magmatic equilibrium with unaltered MORB [d18O(WR) ~5.6-5.7 per mil], and is consistent with the previously determined value for equilibrium with the mantle. The narrow range of measured d18O values is predicted for zircon crystallization from variable parent melt compositions and temperatures in a closed system, and provides no indication of any interactions between altered rocks or seawater and the evolved parent melts. If plagiogranite forms by hydrous partial melting, the uniform mantle-like d18O(Zrc) requires melting and zircon crystallization prior to significant amounts of water-rock interactions that alter the protolith d18O. Zircons from ocean crust have been proposed as a tectonic analog for >3.9 Ga detrital zircons from the earliest (Hadean) Earth by multiple workers. However, zircons from ocean crust are readily distinguished geochemically from zircons formed in continental crustal environments. Many of the >3.9 Ga zircons have mildly elevated d18O (6.0-7.5 per mil), but such values have not been identified in any zircons from the large sample suite examined here. The difference in d18O, in combination with newly acquired lithium concentrations and published trace element data, clearly shows that the >3.9 Ga detrital zircons did not originate by processes analogous to those in modern mid-ocean ridge settings.

Hydrogen isotope composition of pore waters and interlayer water in sediments of ODP Leg 129 holes

Hydrogen isotope compositions have been measured on pore waters from sediments of Leg 129 sites in the Pigafetta and East Mariana basins (central western Pacific). Total water (pore + sorbed waters) contents and their dD have been analyzed for three samples that contain smectite but no zeolite so that sorbed water can be attributed to interlayer water. The H budget for pore and total waters implies that interlayer water is 20 per mil to 30 per mil depleted in D compared to pore water. Because the interlayer/total water molar ratio (0.25 to 0.5) in smectitic sediments is very high, interlayer water represents an important reservoir of D-depleted water in sediments. dD depth profiles for pore water at Sites 800 and 801 show breaks related to chert and radiolarite layers and are relatively vertical below. Above these chert units, pore waters are similar to modern seawater but below, they are between -10 per mil and -5.5 per mil. These values could represent little modified pre-Miocene seawater values, which were D-depleted because of the absence of polar caps, and were preserved from diffusive exchange with modern seawater by the relatively impermeable overlying chert layers. At Site 802, dD values of the pore waters show a decrease in the Miocene tuffs from 0 per mil values at the top to -8 per mil at 250 mbsf. Below, dD values are relatively uniform at about -8ë. Miocene tuffs are undergoing low water/rock alteration. A positive covariation of dD and Cl content of pore water in the tuffs suggests that the increase of dD values could result from secondary smectite formation. Low diffusive exchange coupled with D enrichment due to alteration of preglacial waters could explain the observed profile.