The calculation of Afρ and mass loss rate for comets

Ab initio calculations of Afρ are presented using Mie scattering theory and a Direct Simulation Monte Carlo (DSMC) dust outflow model in support of the Rosetta mission and its target 67P/Churyumov-Gerasimenko (CG). These calculations are performed for particle sizes ranging from 0.010 μm to 1.0 cm. The present status of our knowledge of various differential particle size distributions is reviewed and a variety of particle size distributions is used to explore their effect on Afρ , and the dust mass production View the MathML sourcem˙. A new simple two parameter particle size distribution that curtails the effect of particles below 1 μm is developed. The contributions of all particle sizes are summed to get a resulting overall Afρ. The resultant Afρ could not easily be predicted a priori and turned out to be considerably more constraining regarding the mass loss rate than expected. It is found that a proper calculation of Afρ combined with a good Afρ measurement can constrain the dust/gas ratio in the coma of comets as well as other methods presently available. Phase curves of Afρ versus scattering angle are calculated and produce good agreement with observational data. The major conclusions of our calculations are: – The original definition of A in Afρ is problematical and Afρ should be: qsca(n,λ)×p(g)×f×ρqsca(n,λ)×p(g)×f×ρ. Nevertheless, we keep the present nomenclature of Afρ as a measured quantity for an ensemble of coma particles.– The ratio between Afρ and the dust mass loss rate View the MathML sourcem˙ is dominated by the particle size distribution. – For most particle size distributions presently in use, small particles in the range from 0.10 to 1.0 μm contribute a large fraction to Afρ. – Simplifying the calculation of Afρ by considering only large particles and approximating qsca does not represent a realistic model. Mie scattering theory or if necessary, more complex scattering calculations must be used. – For the commonly used particle size distribution, dn/da ∼ a−3.5 to a−4, there is a natural cut off in Afρ contribution for both small and large particles. – The scattering phase function must be taken into account for each particle size; otherwise the contribution of large particles can be over-estimated by a factor of 10. – Using an imaginary index of refraction of i = 0.10 does not produce sufficient backscattering to match observational data. – A mixture of dark particles with i ⩾ 0.10 and brighter silicate particles with i ⩽ 0.04 matches the observed phase curves quite well. – Using current observational constraints, we find the dust/gas mass-production ratio of CG at 1.3 AU is confined to a range of 0.03–0.5 with a reasonably likely value around 0.1.

Redox front geochemistry and weathering: theory with application to the Osamu Utsumi uranium mine, Poços de Caldas, Brazil

Numerical calculations describing weathering of the Poços de Caldas alkaline complex (Minas Gerais, Brazil) by infiltrating groundwater are carried out for time spans up to two million years in the absence of pyrite, and up to 500,000 years with pyrite present. Deposition of uranium resulting from infiltration of oxygenated, uranium bearing groundwater through the hydrothermally altered phonolitic host rock at the Osamu Utsumi uranium mine is also included in the latter calculation. The calculations are based on the quasi-stationary state approximation to mass conservation equations for pure advective transport. This approximation enables the prediction of solute concentrations, mineral abundances and porosity as functions of time and distance over geologic time spans. Mineral reactions are described by kinetic rate laws for both precipitation and dissolution. Homogeneous equilibrium is assumed to be maintained within the aqueous phase. No other constraints are imposed on the calculations other than the initial composition of the unaltered host rock and the composition of the inlet fluid, taken as rainwater modified by percolation through a soil zone. The results are in qualitative agreement with field observations at the Osamu Utsumi uranium mine. They predict a lateritic cover followed by a highly porous saprolitic zone, a zone of oxidized rock with pyrite replaced by iron-hydroxide, a sharp redox front at which uranium is deposited, and the reduced unweathered host rock. Uranium is deposited in a narrow zone located on the reduced side of the redox front in association with pyrite, in agreement with field observations. The calculations predict the formation of a broad dissolution front of primary kaolinite that penetrates deep into the host rock accompanied by the precipitation of secondary illite. Secondary kaolinite occurs in a saprolitic zone near the surface and in the vicinity of the redox front. Gibbsite forms a bi-modal distribution consisting of a maximum near the surface followed by a thin tongue extending downward into the weathered profile in agreement with field observations. The results are found to be insensitive to the kinetic rate constants used to describe mineral reactions.

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.