Relating in situ gas measurements to the surface outgassing properties of cometary nuclei

The sensitivity of the gas flow field to changes in different initial conditions has been studied for the case of a highly simplified cometary nucleus model. The nucleus model simulated a homogeneously outgassing sphere with a more active ring around an axis of symmetry. The varied initial conditions were the number density of the homogeneous region, the surface temperature, and the composition of the flow (varying amounts of H2O and CO2) from the active ring. The sensitivity analysis was performed using the Polynomial Chaos Expansion (PCE) method. Direct Simulation Monte Carlo (DSMC) was used for the flow, thereby allowing strong deviations from local thermal equilibrium. The PCE approach can be used to produce a sensitivity analysis with only four runs per modified input parameter and allows one to study and quantify non-linear responses of measurable parameters to linear changes in the input over a wide range. Hence the PCE allows one to obtain a functional relationship between the flow field properties at every point in the inner coma and the input conditions. It is for example shown that the velocity and the temperature of the background gas are not simply linear functions of the initial number density at the source. As probably expected, the main influence on the resulting flow field parameter is the corresponding initial parameter (i.e. the initial number density determines the background number density, the temperature of the surface determines the flow field temperature, etc.). However, the velocity of the flow field is also influenced by the surface temperature while the number density is not sensitive to the surface temperature at all in our model set-up. Another example is the change in the composition of the flow over the active area. Such changes can be seen in the velocity but again not in the number density. Although this study uses only a simple test case, we suggest that the approach, when applied to a real case in 3D, should assist in identifying the sensitivity of gas parameters measured in situ by, for example, the Rosetta spacecraft to the surface boundary conditions and vice versa.

Comparison of (31)P saturation and inversion magnetization transfer in human liver and skeletal muscle using a clinical MR system and surface coils.

(31)P MRS magnetization transfer ((31)P-MT) experiments allow the estimation of exchange rates of biochemical reactions, such as the creatine kinase equilibrium and adenosine triphosphate (ATP) synthesis. Although various (31)P-MT methods have been successfully used on isolated organs or animals, their application on humans in clinical scanners poses specific challenges. This study compared two major (31)P-MT methods on a clinical MR system using heteronuclear surface coils. Although saturation transfer (ST) is the most commonly used (31)P-MT method, sequences such as inversion transfer (IT) with short pulses might be better suited for the specific hardware and software limitations of a clinical scanner. In addition, small NMR-undetectable metabolite pools can transfer MT to NMR-visible pools during long saturation pulses, which is prevented with short pulses. (31)P-MT sequences were adapted for limited pulse length, for heteronuclear transmit-receive surface coils with inhomogeneous B1 , for the need for volume selection and for the inherently low signal-to-noise ratio (SNR) on a clinical 3-T MR system. The ST and IT sequences were applied to skeletal muscle and liver in 10 healthy volunteers. Monte-Carlo simulations were used to evaluate the behavior of the IT measurements with increasing imperfections. In skeletal muscle of the thigh, ATP synthesis resulted in forward reaction constants (k) of 0.074 ± 0.022 s(-1) (ST) and 0.137 ± 0.042 s(-1) (IT), whereas the creatine kinase reaction yielded 0.459 ± 0.089 s(-1) (IT). In the liver, ATP synthesis resulted in k = 0.267 ± 0.106 s(-1) (ST), whereas the IT experiment yielded no consistent results. ST results were close to literature values; however, the IT results were either much larger than the corresponding ST values and/or were widely scattered. To summarize, ST and IT experiments can both be implemented on a clinical body scanner with heteronuclear transmit-receive surface coils; however, ST results are much more robust against experimental imperfections than the current implementation of IT.

Helium in solubility equilibrium with quartz and porefluids in rocks - a new approach in hydrology

Quartz crystals in sandstones at depths of 1200 m–1400 m below the surface appear to reach a solubility equilibrium with the 4He-concentration in the surrounding pore- or groundwater after some time. A rather high 4Heconcentration of 4.5x10E-3 cc STP 4He/cm3 of water measured in a groundwater sample would for instance maintain a He pressure of 0.47 atm in a related volume. This value is equal within analytical error to the pressure deduced from the measured helium content of the quartz and its internal helium-accessible volume. To determine this volume, quartz crystals of 0.1 to 1 mm were separated from sandstones and exposed to a helium gas pressure of 32 atm at a temperature of 290°C for up to 2 months. By crushing, melting or isothermal heating the helium was then extracted from the helium saturated samples. Avolume on the order of 0.1% of the crystal volume is only accessible to helium atoms but not to argon atoms or water molecules. By monitoring the diffusive loss of He from the crystals at 350°C an effective diffusion constant on the order of 10E-9 cm2/s is estimated. Extrapolation to the temperature of 70°C in the sediments at a depth of 1400 m gives a typical time of about 100 000 years to reach equilibrium between helium in porewaters and the internal He-accessible volume of quartz crystals. In a geologic situation with stagnant pore- or groundwaters in sediments it therefore appears to be possible with this new method to deduce a 4He depth profile for porewaters in impermeable rocks based on their mineral record.

Barium isotope fractionation during witherite (BaCO3) dissolution, precipitation and at equilibrium

This study examines the behavior of Ba isotope fractionation between witherite and fluid during mineral dissolution, precipitation and at chemical equilibrium. Experiments were performed in batch reactors at 25 oC in 10-2 M NaCl solution where the pH was adjusted by continuous bubbling of a water saturated gas phase of CO2 or atmospheric air. During witherite dissolution no Ba isotope fractionation was observed between solid and fluid. In contrast, during witherite precipitation, caused by a pH increase, a preferential uptake of the lighter 134Ba isotopomer in the solid phase was observed. In this case, the isotope fractionation factor αwitherite-fluid is calculated to be 0.99993 ± 0.00004 (or Δ137/134Bawitherite-fluid ≈ -0.07 ± 0.04 ‰, 2sd). The most interesting feature of this study, however, is that after the attainment of chemical equilibrium, the Ba isotope composition of the aqueous phase is progressively becoming lighter, indicating a continuous exchange of Ba2+ ions between witherite and fluid. Mass balance calculations indicate that the detachment of Ba from the solid is not only restricted to the outer surface layer of the solid, but affects several (~7 unit cells) subsurface layers of the crystal. This observation comes in excellent agreement with the concept of a dynamic system at chemical equilibrium in a mineral-fluid system, denoting that the time required for the achievement of isotopic equilibrium in the witherite-fluid system is longer compared to that observed for chemical equilibrium. Overall, these results indicate that the isotopic composition of Ba bearing carbonates in natural environments may be altered due to changes in fluid composition without a net dissolution/precipitation to be observed.