Ductile deformation during carbonation of serpentinized peridotite.

Carbonated serpentinites (listvenites) in the Samail Ophiolite, Oman, record mineralization of 1-2 Gt of CO2, but the mechanisms providing permeability for continued reactive fluid flow are unclear. Based on samples of the Oman Drilling Project, here we show that listvenites with a penetrative foliation have abundant microstructures indicating that the carbonation reaction occurred during deformation. Folded magnesite veins mark the onset of carbonation, followed by deformation during carbonate growth. Undeformed magnesite and quartz overgrowths indicate that deformation stopped before the reaction was completed. We propose deformation by dilatant granular flow and dissolution-precipitation assisted the reaction, while deformation in turn was localized in the weak reacting mass. Lithostatic pore pressures promoted this process, creating dilatant porosity for CO2 transport and solid volume increase. This feedback mechanism may be common in serpentinite-bearing fault zones and the mantle wedge overlying subduction zones, allowing massive carbonation of mantle rocks.

Nested multiresolution hierarchical simulated annealing algorithm for porous media reconstruction.

Microstructure strongly influences flow and transport properties of porous media. Flow and transport simulations within porous media, therefore, requires accurate three-dimensional (3D) models of the pore and solid phase structure. To date, no imaging method can resolve all relevant heterogeneities from the nano- to the centimeter scale within complex heterogeneous materials such as clay, reservoir rocks (e.g., travertine, chalk, ...), hardened cement paste, and concrete. To reconstruct these porous materials it is thus necessary to merge information from different 2D and potentially 3D imaging methods. One porous media reconstruction methodology that has been around for at least two decades is simulated annealing (SA). However, realizations with SA typically suffer an artificially reduced long-range connectivity, while multiphase reconstructions are not feasible in most cases because of a prohibitive computational burden. To solve these problems we propose a hierarchical multiresolution and multiphase simulated annealing algorithm. To decrease the computational cost of multiphase simulation, our algorithm sequentially simulates one phase after another, in a hierarchical way, which enables handling multimodal distributions and topological relations. Building upon recent work, our algorithm improves long-range connectivity and CPU efficiency by simulating larger particles using a coarser resolution that is subsequently refined compared to standard SA; our proposed extension not only offers the possibility to perform multiphase reconstruction but also allows us (i) to improve binary reconstruction quality, as quantified, e.g., by multiple-point histograms by up to one order of magnitude and (ii) to achieve an overall speed-up. The proposed algorithm is also shown to outperform the direct sampling multiple-point statistics method for the generation of cement paste microstructure with respect to both generation time and quality.

Structure and evolution of a rocksalt-mudrock-tectonite: The haselgebirge in the Northern Calcareous Alps.

The Northern Calcareous Alps are part of the Eastern Alps in Austria and Germany. The Mesozoic units of this fold-and-thrust belt were detached, thrusted and stacked along the evaporitic Haselgebirge Formation. Exposed in salt mines, rocksalt and mudrock form a two component tectonite: The rock type "haselgebirge" consists of 10-70 wt % halite with silt- to gravel- or block-sized components within a halite matrix, and the "kerngebirge" with >70 wt % halite. All rock types studied are fault rocks. By use of a temperature-independent subgrain size piezometer, the paleo-differential stress of halite was calculated at ca. 2.5 MPa in Altaussee and ca. 4.5 MPa in Berchtesgaden. Including data from a grain-size piezometer, temperatures were estimated at ca. 150 ± 20 °C and 110 ± 10 °C. This implies very high strain rates, which are about 10-10-10-9 s-1. During the tectonic movement, the halite deformed, recrystallized, and crystallized as veins in mudrock fractures. We interpret high overpressure of the pore fluid to have significantly contributed to fracturing of the mudrock.

Weakening of rock salt by water during long-term creep.

The rheological properties of rock salt are of fundamental importance in predicting the long-term evolution of salt-based radioactive waste repositories and strategic storage caverns, and in modelling the formation of salt diapirs and associated oil traps1,2. The short-term, high-stress rheology of rock salt is well known from laboratory experiments; however, extrapolation to appropriately low stresses fails to predict the rapid flow seen in certain natural structures. Furthermore, experiments have failed to reproduce the recrystallized microstructure of naturally deformed salt. Here we report experiments indicating that the above discrepancies can be explained by taking into account the influence of trace amounts of brine. Trace brine is always present in natural salt but sometimes escapes during experiments. Our tests on dry dilated salt show more or less conventional dislocation creep behaviour, but brine-bearing samples show marked weakening at low strain rates. This is associated with dynamic recrystallization and a change of deformation mechanism to solution transfer creep. Because natural rock salt always contains some brine, these results cast substantial doubt on the validity of presently accepted dislocation creep laws for predicting the long-term rheological behaviour of salt in nature.