Mineralogical and isotopic record of diagenesis from the Opalinus Clay formation at Benken, Switzerland: Implications for the modeling of pore-water chemistry in a clay formation

Argillaceous rocks are considered to be a suitable geological barrier for the long-term containment of wastes. Their efficiency at retarding contaminant migration is assessed using reactive-transport experiments and modeling, the latter requiring a sound understanding of pore-water chemistry. The building of a pore-water model, which is mandatory for laboratory experiments mimicking in situ conditions, requires a detailed knowledge of the rock mineralogy and of minerals at equilibrium with present-day pore waters. Using a combination of petrological, mineralogical, and isotopic studies, the present study focused on the reduced Opalinus Clay formation (Fm) of the Benken borehole (30 km north of Zurich) which is intended for nuclear-waste disposal in Switzerland. A diagenetic sequence is proposed, which serves as a basis for determining the minerals stable in the formation and their textural relationships. Early cementation of dominant calcite, rare dolomite, and pyrite formed by bacterial sulfate reduction, was followed by formation of iron-rich calcite, ankerite, siderite, glauconite, (Ba, Sr) sulfates, and traces of sphalerite and galena. The distribution and abundance of siderite depends heavily on the depositional environment (and consequently on the water column). Benken sediment deposition during Aalenian times corresponds to an offshore environment with the early formation of siderite concretions at the water/sediment interface at the fluctuating boundary between the suboxic iron reduction and the sulfate reduction zones. Diagenetic minerals (carbonates except dolomite, sulfates, silicates) remained stable from their formation to the present. Based on these mineralogical and geochemical data, the mineral assemblage previously used for the geochemical model of the pore waters at Mont Terri may be applied to Benken without significant changes. These further investigations demonstrate the need for detailed mineralogical and geochemical study to refine the model of pore-water chemistry in a clay formation.

Benchmark problems for reactive transport modeling of the generation and attenuation of acid rock drainage

Acid rock drainage (ARD) is a problem of international relevance with substantial environmental and economic implications. Reactive transport modeling has proven a powerful tool for the process-based assessment of metal release and attenuation at ARD sites. Although a variety of models has been used to investigate ARD, a systematic model intercomparison has not been conducted to date. This contribution presents such a model intercomparison involving three synthetic benchmark problems designed to evaluate model results for the most relevant processes at ARD sites. The first benchmark (ARD-B1) focuses on the oxidation of sulfide minerals in an unsaturated tailing impoundment, affected by the ingress of atmospheric oxygen. ARD-B2 extends the first problem to include pH buffering by primary mineral dissolution and secondary mineral precipitation. The third problem (ARD-B3) in addition considers the kinetic and pH-dependent dissolution of silicate minerals under low pH conditions. The set of benchmarks was solved by four reactive transport codes, namely CrunchFlow, Flotran, HP1, and MIN3P. The results comparison focused on spatial profiles of dissolved concentrations, pH and pE, pore gas composition, and mineral assemblages. In addition, results of transient profiles for selected elements and cumulative mass loadings were considered in the intercomparison. Despite substantial differences in model formulations, very good agreement was obtained between the various codes. Residual deviations between the results are analyzed and discussed in terms of their implications for capturing system evolution and long-term mass loading predictions.

Automatic multi-resolution shape modeling of multi-organ structures.

Point Distribution Models (PDM) are among the most popular shape description techniques and their usefulness has been demonstrated in a wide variety of medical imaging applications. However, to adequately characterize the underlying modeled population it is essential to have a representative number of training samples, which is not always possible. This problem is especially relevant as the complexity of the modeled structure increases, being the modeling of ensembles of multiple 3D organs one of the most challenging cases. In this paper, we introduce a new GEneralized Multi-resolution PDM (GEM-PDM) in the context of multi-organ analysis able to efficiently characterize the different inter-object relations, as well as the particular locality of each object separately. Importantly, unlike previous approaches, the configuration of the algorithm is automated thanks to a new agglomerative landmark clustering method proposed here, which equally allows us to identify smaller anatomically significant regions within organs. The significant advantage of the GEM-PDM method over two previous approaches (PDM and hierarchical PDM) in terms of shape modeling accuracy and robustness to noise, has been successfully verified for two different databases of sets of multiple organs: six subcortical brain structures, and seven abdominal organs. Finally, we propose the integration of the new shape modeling framework into an active shape-model-based segmentation algorithm. The resulting algorithm, named GEMA, provides a better overall performance than the two classical approaches tested, ASM, and hierarchical ASM, when applied to the segmentation of 3D brain MRI.

Atomistic Modeling of Effect of Mg on Oxygen Vacancy Diffusion in α-Alumina

Oxygen diffusion plays an important role in grain growth and densification during the sintering of alumina ceramics and governs high-temperature processes such as creep. The atomistic mechanism for oxygen diffusion in alumina is, however, still debated; atomistic calculations not being able to match experimentally determined activation energies for oxygen vacancy diffusion. These calculations are, however, usually performed for perfectly pure crystals, whereas virtually every experimental alumina sample contains a significant fraction of impurity/dopants ions. In this study, we use atomistic defect cluster and nudged elastic band (NEB) calculations to model the effect of Mg impurities/dopants on defect binding energies and migration barriers. We find that oxygen vacancies can form energetically favorable clusters with Mg, which reduces the number of mobile species and leads to an additional 1.5 eV energy barrier for the detachment of a single vacancy from Mg. The migration barriers of diffusive jumps change such that an enhanced concentration of oxygen vacancies is expected around Mg ions. Mg impurities were also found to cause destabilization of certain vacancy configurations as well as enhanced vacancy–vacancy interaction.

3D-modeling of Mercury's solar wind sputtered surface-exosphere environment

The efficiency of sputtered refractory elements by H+ and He++ solar wind ions from Mercury's surface and their contribution to the exosphere are studied for various solar wind conditions. A 3D solar wind-planetary interaction hybrid model is used for the evaluation of precipitation maps of the sputter agents on Mercury's surface. By assuming a global mineralogical surface composition, the related sputter yields are calculated by means of the 2013 SRIM code and are coupled with a 3D exosphere model. Because of Mercury's magnetic field, for quiet and nominal solar wind conditions the plasma can only precipitate around the polar areas, while for extreme solar events (fast solar wind, coronal mass ejections, interplanetary magnetic clouds) the solar wind plasma has access to the entire dayside. In that case the release of particles form the planet's surface can result in an exosphere density increase of more than one order of magnitude. The corresponding escape rates are also about an order of magnitude higher. Moreover, the amount of He++ ions in the precipitating solar plasma flow enhances also the release of sputtered elements from the surface in the exosphere. A comparison of our model results with MESSENGER observations of sputtered Mg and Ca elements in the exosphere shows a reasonable quantitative agreement. (C) 2015 Elsevier Ltd. All rights reserved.

Collisional Formation and Modeling of Asteroid Families

In the last decade, thanks to the development of sophisticated numerical codes, major breakthroughs have been achieved in our understanding of the formation of asteroid families by catastrophic disruption of large parent bodies. In this review, we describe numerical simulations of asteroid collisions that reproduced the main properties of families, accounting for both the fragmentation of an asteroid at the time of impact and the subsequent gravitational interactions of the generated fragments. The simulations demonstrate that the catastrophic disruption of bodies larger than a few hundred meters in diameter leads to the formation of large aggregates due to gravitational reaccumulation of smaller fragments, which helps explain the presence of large members within asteroid families. Thus, for the first time, numerical simulations successfully reproduced the sizes and ejection velocities of members of representative families. Moreover, the simulations provide constraints on the family dynamical histories and on the possible internal structure of family members and their parent bodies.