Effect of oxygen fugacity on the H2O storage capacity of forsterite in the carbon-saturated systems

High pressure experiments have been performed in the systems Mg₂SiO₄-C-O-H and Mg₂SiO₄-K₂CO₃-C at 6.3 GPa and 1200 to 1600 °C using a split-sphere multi-anvil apparatus. In the Mg₂SiO₄-C-O-H system the composition of fluid was modeled by adding different amounts of water and stearic acid. The fO₂ was controlled by the Mo-MoO₂ or Fe-FeO oxygen buffers. Several experiments in the Mg₂SiO₄-C-O-H system and all experiments in the Mg₂SiO₄-K₂CO₃-C system have been conducted without buffering the fO₂. Forsterite in the system Mg₂SiO₄-K₂CO₃-C does not reveal OH absorption bands in the IR spectra, while forsterite coexisting with carbon-bearing fluid and silicate melt at logfO₂ from FMQ-2 to FMQ-5 (from 2 to 5 log units below fayalite-magnetite-quartz oxygen buffer) contains 800-1850 wt. ppm H₂O. The maximum concentrations were detected at 1400 °C and FMQ-3.5. We observed an increase in the solidus temperature in the system Mg₂SiO₄-C-O-H from 1200 to above 1600 °C with log fO₂ decreasing from FMQ-2 to FMQ-5. The increase of the solidus temperature and the broadening of the stability field of the H₂O-H₂-CH₄ subsolidus fluid phase at 1400-1600 °C explain the high H₂O storage capacity of forsterite relative to that crystallized from carbon-free, oxidized, hydrous, silicic melt. At temperatures above 1400 °C liquidus forsterite precipitated along with diamond from oxidized (FMQ-1) carbonate-silicate melt and from silicate melt dissolving the moderately reduced C-O-H fluid (from FMQ-2 to FMQ-3.5). Formation of diamond was not detected under ultra-reduced conditions (FMQ-5) at 1200-1600 °C. Olivine co-precipitating with diamond from dry carbonate-silicate or hydrous-silicic fluid/melt can provide information on the H₂O contents and speciation of the diamond-forming media in the mantle. The conditions for minimum post-crystallization alteration of olivine and its hydrogen content are discussed.     

Large‐scale magnetic flux concentrations from turbulent stresses

In this study we provide the first numerical demonstration of the effects of turbulence on the mean Lorentz force and the resulting formation of large‐scale magnetic structures. Using three‐dimensional direct numerical simulations (DNS) of forced turbulence we show that an imposed mean magnetic field leads to a decrease of the turbulent hydromagnetic pressure and tension. This phenomenon is quantified by determining the relevant functions that relate the sum of the turbulent Reynolds and Maxwell stresses with the Maxwell stress of the mean magnetic field. Using such a parameterization, we show by means of two‐dimensional and three‐dimensional mean‐field numerical modelling that an isentropic density stratified layer becomes unstable in the presence of a uniform imposed magnetic field. This large‐scale instability results in the formation of loop‐like magnetic structures which are concentrated at the top of the stratified layer. In three dimensions these structures resemble the appearance of bipolar magnetic regions in the Sun. The results of DNS and mean‐field numerical modelling are in good agreement with theoretical predictions. We discuss our model in the context of a distributed solar dynamo where active regions and sunspots might be rather shallow phenomena (© 2010 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Pore to pore validation of pore network modelling against micromodel experiment results

Pore network modelling (PNM) has been widely used to study the multiphase flow and transport in porous media. Although a number of recent papers discussed the PNM validation on core-scale parameters such as permeability, relative permeability and capillary pressure; quantitative predictive potential of PNM on pore by pore basis has rarely been studied. The aim of this paper is to present a direct comparison between PNM simulations and corresponding micro-model experiments at the same scale and the same geometry. A number of well-defined and constrained two-phase flow in porous medium experimental scenarios were utilized to validate the physics solving part in PNM (filling rules, capillary and viscous pressure). This work validates that a dynamic pore network flow solver can predict two-phase flow displacements for these experiments for drainage situations at both pore and plug scales. A glass-etched micro-model is used to quantify the accuracy of a dynamic PNM solver on pore and core levels. Two-phase drainage micro fluidic experiments at different flow conditions are performed on micro-models. PNM simulations are performed on the same pattern and flow conditions as used in micro-model experiments. The two-phase distribution extracted from experiment images is registered onto rsults of PNM simulations for direct pore to pore comparison. Pore-scale matching level is found at around 75 % for all three test cases. The matching level of core-scale parameters such as S _{w c} and oil-phase permeability varies from case to case; the relative error to micro-model experiment measurements varies from 15 to 60 %. Possible reasons leading to discrepancies on core-scale parameters are discussed: missing considerations during validation of the combination of uncertainty in both simulator input parameters and experiments are seen as the principal factors.     

A re-examination of petrogenesis and 40Ar/39Ar systematics in the Chain of Ponds K-feldspar: “diffusion domain” archetype versus polyphase hygrochronology

K-feldspar (Kfs) from the Chain of Ponds Pluton (CPP) is the archetypal reference material, upon which thermochronological modeling of Ar diffusion in discrete “domains” was founded. We re-examine the CPP Kfs using cathodoluminescence and back-scattered electron imaging, transmission electron microscopy, and electron probe microanalysis. ⁴⁰Ar/³⁹Ar stepwise heating experiments on different sieve fractions, and on handpicked and unpicked aliquots, are compared. Our results reproduce the staircase-shaped age spectrum and the Arrhenius trajectory of the literature sample, confirming that samples collected from the same locality have an identical Ar isotope record. Even the most pristine-looking Kfs from the CPP contains successive generations of secondary, metasomatic/retrograde mineral replacements that post-date magmatic crystallization. These chemically and chronologically distinct phases are responsible for its staircase-shaped age spectra, which are modified by handpicking. While genuine within-grain diffusion gradients are not ruled out by these data, this study demonstrates that the most important control on staircase-shaped age spectra is the simultaneous presence of heterochemical, diachronous post-magmatic mineral growth. At least five distinct mineral species were identified in the Kfs separate, three of which can be traced to external fluids interacting with the CPP in a chemically open system. Sieve fractions have size-shifted Arrhenius trajectories, negating the existence of the smallest “diffusion domains.” Heterochemical phases also play an important role in producing nonlinear trajectories. In vacuo degassing rates recovered from Arrhenius plots are neither related to true Fick’s Law diffusion nor to the staircase shape of the age spectra. The CPP Kfs used to define the “diffusion domain” model demonstrates the predominance of metasomatic alteration by hydrothermal fluids and recrystallization in establishing the natural Ar distribution among different coexisting phases that gives rise to the staircase-shaped age spectrum. Microbeam imaging of textures is as essential for ⁴⁰Ar/³⁹Ar hygrochronology as it is for U–Pb geochronology.     

Mercury’s magnetospheric magnetic field after the first two MESSENGER flybys

The “paraboloid” model of Mercury’s magnetospheric magnetic field is used to determine the best-fit magnetospheric current system and internal dipole parameters from magnetic field measurements taken during the first and second MESSENGER flybys of Mercury on 14 January and 6 October 2008. Together with magnetic field measurements taken during the Mariner 10 flybys on 29 March 1974 and 16 March 1975, there exist three low-latitude traversals separated in longitude and one high-latitude encounter. From our model formulation and fitting procedure a Mercury dipole moment of 196 nT ·  (where R_M is Mercury’s radius) was determined. The dipole is offset from Mercury’s center by 405 km in the northward direction. The dipole inclination to Mercury’s rotation axis is relatively small, ∼4°, with an eastern longitude of 193° for the dipole northern pole. Our model is based on the a priori assumption that the dipole position and the moment orientation and strength do not change in time. The root mean square (rms) deviation between the Mariner 10 and MESSENGER magnetic field measurements and the predictions of our model for all four flybys is 10.7 nT. For each magnetic field component the rms residual is ∼6 nT or about 1.5% of the maximum measured magnetic field, ∼400 nT. This level of agreement is possible only because the magnetospheric current system parameters have been determined separately for each flyby. The magnetospheric stand-off distance, the distance from the planet’s center to the inner edge of the tail current sheet, the tail lobe magnetic flux, and the displacement of the tail current sheet relative to the Mercury solar-magnetospheric equatorial plane have been determined independently for each flyby. The magnetic flux in the tail lobes varied from 3.8 to 5.9 MWb; the subsolar magnetopause stand-off distance from 1.28 to 1.43 R_M; and the distance to the inner edge of the current sheet from 1.23 to 1.32 R_M. The differences in the current systems between the first and second MESSENGER flybys are attributed to the effects of strong magnetic reconnection driven by southward interplanetary magnetic field during the latter flyby.     

Nonhelical mean‐field dynamos in a sheared turbulence

Mechanisms of nonhelical large‐scale dynamos (shear‐current dynamo and effect of homogeneous kinetic helicity fluctuations with zero mean) in a homogeneous turbulence with large‐scale shear are discussed. We have found that the shearcurrent dynamo can act even in random flows with small Reynolds numbers. However, in this case mean‐field dynamo requires small magnetic Prandtl numbers (i.e., when Pm < Pm^cr < 1). The threshold in the magnetic Prandtl number, Pm^cr = 0.24, is determined using second order correlation approximation (or first‐order smoothing approximation) for a background random flow with a scale‐dependent viscous correlation time τ_c = (νk 2)^–1 (where ν is the kinematic viscosity of the fluid and k is the wave number). For turbulent flows with large Reynolds numbers shear‐current dynamo occurs for arbitrary magnetic Prandtl numbers. This dynamo effect represents a very generic mechanism for generating large‐scale magnetic fields in a broad class of astrophysical turbulent systems with large‐scale shear. On the other hand, mean‐field dynamo due to homogeneous kinetic helicity fluctuations alone in a sheared turbulence is not realistic for a broad class of astrophysical systems because it requires a very specific random forcing of kinetic helicity fluctuations that contains, e.g., low‐frequency oscillations. (© 2008 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Modeling Arctic Ocean heat transport and warming episodes in the 20th century caused by the intruding Atlantic Water

This study investigates the Arctic Ocean warming episodes in the 20th century using both a high-resolution coupled global climate model and historical observations. The model, with no flux adjustment, reproduces well the Atlantic Water core temperature （AWCT） in the Arctic Ocean and shows that four largest decadalscale warming episodes occurred in the 1930s, 70s, 80s, and 90s, in agreement with the hydrographic observational data. The difference is that there was no pre-warming prior to the 1930s episode, while there were two pre-warming episodes in the 1970s and 80s prior to the 1990s, leading the 1990s into the largest and prolonged warming in the 20th century. Over the last century, the simulated heat transport via Fram Strait and the Barents Sea was estimated to be, on average, 31.32 TW and 14.82 TW, respectively, while the Bering Strait also provides 15.94 TW heat into the west- ern Arctic Ocean. Heat transport into the Arctic Ocean by the Atlantic Water via Fram Strait and the Barents Sea correlates significantly with AWCT （ C = 0.75 ） at 0- lag. The modeled North Atlantic Oscillation （NAO） index has a significant correlation with the heat transport （ C = 0.37 ）. The observed AWCT has a significant correlation with both the modeled AWCT （ C =0.49） and the heat transport （ C =0.41 ）. However, the modeled NAO index does not significantly correlate with either the observed AWCT （ C = 0.03 ） or modeled AWCT （ C = 0.16 ） at a zero-lag, indicating that the Arctic climate system is far more complex than expected.     

V475 Sct (Nova Scuti 2003) – Binary or Triple System?

UBV RI photometry and 4600–9000 Å spectroscopy of nova V475 Sct taken in the first 3 month after discovery is presented. The object can be classified as a Fe II type slow nova with t _2,V = 48 days, t _3,V = 53 days. The absolute magnitude of the nova at maximum, its colour excess and distance were determined. The observed 13.4 day periodicity of flares can be explained by the mass transfer bursts from the red to the white dwarf, probably caused by the periastron passage of a third body. Two sets of absorptions are seen in the P Cyg-type Hα line profile. They arise in the expanding shell of the nova.     

Mo isotope composition in Mo-rich high- and low-T hydrothermal systems from the Swiss Alps

We analyzed the molybdenum (Mo) isotope compositions (IC) of 59 samples from two molybdenite mineralizations (Alpjahorn and Grimsel) and from a Mo-rich hydrothermal breccia (Grimsel) from the Aar Massif, Switzerland. The formation temperature of the Late Paleozoic Mo mineralizations (300-600 °C) is much higher than that of the Pliocene breccia (100-160 °C). The Mo IC of the molybdenites varies over 1.35‰. Even in a single hand specimen it spans 0.45‰, indicating that fractionation processes during molybdenite precipitation can vary on a cm scale. The Mo IC of most molybdenites analyzed here are significantly heavier than that of the host rock (δ^98/95Mo = (0.05 ± 0.1)‰) and show a bimodal distribution centered around δ^98/95Mo ≈ 1.1‰ and 0.2‰. This result rules out single stage Rayleigh fractionation as the relevant formation mechanism and instead, redox variations are suggested to be a main factor controlling the Mo IC of the studied high-temperature Mo deposits. The range of the Mo IC in one single deposit, the Alpjahorn, overlaps with the variation range of almost all other published values for Mo IC in Mo deposits. Compared to the molybdenites, the breccia shows an even wider variation of 3.0‰ (δ^98/95Mo between −1.6‰ and +1.4‰). In contrast to the high-T molybdenite deposits, here the Mo was transported via oxidized surface waters into the breccia system, where it was reduced and precipitated. This indicates that oxidation and reduction of Mo complexes may lead to highly variable Mo IC in hydrothermal systems.