Albert Griera

Strain localization and porphyroclast rotation

It has been debated for decades whether rigid inclusions, such as porphyroclasts and porphyroblasts, do or do not rotate in a softer matrix during deformation. Experiments and numerical simulations with viscous matrix rheologies show ongoing rotation of circular inclusions, whereas using Mohr-Coulomb plasticity results in nonrotation. Because the rocks in which inclusions are found normally undergo deformation by dislocation creep, we applied a full-field crystal plasticity approach to investigate the rotation behavior of rigid circular inclusions. We show that the inclusion9s rotation strongly depends on the anisotropy of the matrix minerals. Strongly anisotropic minerals will develop shear bands that reduce the rotation of inclusions. Inhibition of rotation can only occur after a significant amount of strain. Our results may help to explain why geologic rigid objects often show evidence for rotation, but not necessarily in accordance with the viscous theory that is usually applied to these systems.

When do folds unfold during progressive shear

Folded layers in rocks can be stretched again, potentially unfolding the folds back to straight layers. Little is known, however, about how to recognize partly or even entirely unfolded layers. When folded layers can unfold, what determines their mechanical behavior, and how can we recognize them in the field? In order to address these questions, we present a series of numerical simulations of the stretching of previously folded single layers and multilayers in simple shear. Layers do not completely unfold when they undergo softening before or during the stretching process, or when adjacent competent layers prevent them from unfolding. Intrafolial folds and cusp-like folds adjacent to straight layers as well as variations in fold amplitudes and limb lengths of irregular folds are indicative of stretching of a fold train.

Competition between grain growth and grain-size reduction in polar ice

Static (or ‘normal’) grain growth, i.e. grain boundary migration driven solely by grain boundary energy, is considered to be an important process in polar ice. Many ice-core studies report a continual increase in average grain size with depth in the upper hundreds of metres of ice sheets, while at deeper levels grain size appears to reach a steady state as a consequence of a balance between grain growth and grain-size reduction by dynamic recrystallization. The growth factor k in the normal grain growth law is important for any process where grain growth plays a role, and it is normally assumed to be a temperature-dependent material property. Here we show, using numerical simulations with the program Elle, that the factor k also incorporates the effect of the microstructure on grain growth. For example, a change in grain-size distribution from normal to log-normal in a thin section is found to correspond to an increase in k by a factor of 3.5.

Geology of the Cerro Quema Au-Cu deposit (Azuero Peninsula, Panama)

The Cerro Quema district, located on the Azuero Peninsula, Panama, is part of a large regional hydrothermal system controlled by regional faults striking broadly E-W, developed within the Rio Quema Formation. This formation is composed of volcanic, sedimentary and volcano-sedimentary rocks indicating a submarine depositional environment, corresponding to the fore-arc basin of a Cretaceous–Paleogene volcanic arc. The structures observed in the area and their tectono-stratigraphic relationship with the surrounding formations suggest a compressive and/or transpressive tectonic regime, at least during Late Cretaceous–Oligocene times. The igneous rocks of the Rio Quema Formation plot within the calc-alkaline field with trace and rare earth element (REE) patterns of volcanic arc affinity. This volcanic arc developed on the Caribbean large igneous province during subduction of the Farallon Plate. Mineralization consists of disseminations of pyrite and enargite as well as a stockwork of pyrite and barite with minor sphalerite, galena and chalcopyrite, hosted by a subaqueous dacitic lava dome of the Rio Quema Formation. Gold is present as submicroscopic grains and associated with pyrite as invisible gold. A hydrothermal alteration pattern with a core of advanced argillic alteration (vuggy silica with alunite, dickite, pyrite and enargite) and an outer zone of argillic alteration (kaolinite, smectite and illite) has been observed. Supergene oxidation overprinted the hydrothermal alteration resulting in a thick cap of residual silica and iron oxides. The ore minerals, the alteration pattern and the tectono-volcanic environment of Cerro Quema are consistent with a high sulfidation epithermal system developed in the Azuero peninsula during pre-Oligocene times

Insights to controls on dolomitization by means of reactive transport models applied to the Benicàssim case study (Maestrat Basin, eastern Spain)

Partially dolomitized carbonate rocks of the Middle East and North America host large hydrocarbon reserves. The origin of some of these dolomites has been attributed to a hydrothermal mechanism. The Benicàssim area (Maestrat Basin, eastern Spain) constitutes an excellent field analogue for fault-controlled stratabound hydrothermal dolomitization: dolostone geobodies are well exposed and extend over several kilometres away from seismic-scale faults. This work investigates the main controls on the formation of stratabound versus massive dolomitization in carbonate sequences by means of two-dimensional (2D) reactive transport models applied to the Benicàssim case study. Simulation results suggest that the dolomitization capacity of Mg-rich fluids reaches a maximum at temperatures around 100 °C and a minimum at 25 °C (studied temperature range: 25–150 °C). It takes of the order of hundreds of thousands to millions of years to completely dolomitize kilometre-long limestone sections, with solutions flowing laterally through strata at velocities of metres per year (m/a). Permeability differences of two orders of magnitude between layers are required to form stratabound dolomitization. The kilometre-long stratabound dolostone geobodies of Benicàssim must have formed under a regime of lateral flux greater than metres per year over about a million years. As long-term dolomitization tends to produce massive dolostone bodies not seen at Benicàssim, the dolomitizing process there must have been limited by the availability of fluid volume or the flow-driving mechanism. Reactive transport simulations have proven a useful tool to quantify aspects of the Benicàssim genetic model of hydrothermal dolomitization.

Strain and deformation history in a syntectonic pluton. The case of the Roses granodiorite (Cap de Creus, Eastern Pyrenees)

Abstract The Roses granodiorite is a Variscan stock with well developed syn- and post-magmatic deformation structures that crops out in the Pyrenean Axial Zone. Analysis of structures reveals a continuous deformation history during and after magma cooling. The deformation history is divided on the basis of mechanical behaviour into two stages: an early one with the development of magmatic structures and a late stage with the development of mylonitic fabrics along shear zones. Both stages are separated in time by the emplacement of aplite-pegmatite dykes. Time of dyke emplacement is thought to coincide with a sudden change in rheology of the granodiorite. The abundance of quartz dioritic enclaves permits the use of shape analysis to characterize the magmatic fabric as a homogeneous deformation. Later solid-state deformation led to the development of an inhomogeneous deformation pattern with different sizes of anastomosing shear zones wrapping around lozenge-shaped domains. The displacement/width ratio measured in shear zones ranges between one and two orders of magnitude. The Roses granodiorite is thought to be a synkinematically emplaced stock which records a continuous deformatioinal history with two distinct deformation stages, both recording bulk finite strains of similar order of magnitude but with a marked difference in finite strain distribution.

What happens to deformed rocks after deformation? A refined model for recovery based on numerical simulations

Abstract Deformation, in large parts of the middle crust, results in strained rocks consisting of grains with variable dislocation densities and microstructures which are characterized by gradual distortion and subgrain structures. Post-deformation residence of these rocks at elevated temperatures results in microstructural adjustments through static recovery and recrystallization. Here, we employ a numerical technique to simulate intragrain recovery at temperatures at or below the deformation temperature. The simulation is based on minimization of the stored energy, related to misorientation through local rotation of physical material points relative to their immediate environment. Three temperature- and/or deformation-geometry-dependent parameters were systematically varied: (1) deformation-induced dislocation types, (2) dislocation mobility and (3) size of dislocation interaction volume. Comparison with previously published in situ experiments shows consistency of numerical and experimental results. They show temperature- and dislocation-type-dependent small-scale fluctuations in subgrain-boundary misorientations and orientation variation within subgrains. These can be explained by the combined effect of increase in dislocation interaction volume and activation of climb. Our work shows microstructure can be significantly modified even if the post-deformational temperature is at or below the deformation temperature: a scenario relevant for most deformed rocks.

Porphyroblast rotation versus nonrotation: Conflict resolution! COMMENT

[Fay et al. (2008)][1] present a few numerical simulations (with unspecified code or boundary conditions) as evidence for nonrotation of porphyroblasts during non-coaxial flow. However, the large body of work on the behavior of porphyroblasts in deforming ductile rocks is significantly undervalued.

Dynamic recrystallization during deformation of polycrystalline ice: insights from numerical simulations

The flow of glaciers and polar ice sheets is controlled by the highly anisotropic rheology of ice crystals that have hexagonal symmetry (ice lh). To improve our knowledge of ice sheet dynamics, it is necessary to understand how dynamic recrystallization (DRX) controls ice microstructures and rheology at different boundary conditions that range from pure shear flattening at the top to simple shear near the base of the sheets. We present a series of two-dimensional numerical simulations that couple ice deformation with DRX of various intensities, paying special attention to the effect of boundary conditions. The simulations show how similar orientations of c-axis maxima with respect to the finite deformation direction develop regardless of the amount of DRX and applied boundary conditions. In pure shear this direction is parallel to the maximum compressional stress, while it rotates towards the shear direction in simple shear. This leads to strain hardening and increased activity of non-basal slip systems in pure shear and to strain softening in simple shear. Therefore, it is expected that ice is effectively weaker in the lower parts of the ice sheets than in the upper parts. Strain-rate localization occurs in all simulations, especially in simple shear cases. Recrystallization suppresses localization, which necessitates the activation of hard, non-basal slip systems. This article is part of the themed issue ‘Microdynamics of ice’.

The Application of In Situ 3D X-Ray Diffraction in Annealing Experiments: First Interpretation of Substructure Development in Deformed NaCl

n-situ 3D X-ray diffraction (3DXRD) annealing experiments were conducted at the ID-11 beamline at the European Synchrotron Radiation Facility in Grenoble. This allowed us to non-destructively document and subsequently analyse the development of substructures during heating, without the influence of surface effects. A sample of deformed single crystal halite was heated to between 260-400 °C. Before and after heating a volume of 500 by 500 by 300 μm was mapped using a planar beam, which was translated over the sample volume at intervals of 5-10 µm in the vertical dimension. In the following we present partially reconstructed orientation maps over one layer before and after heating for 240min at 260 °C. Additional small syn-heating maps over a constrained sample rotation of 12-30º. The purpose of this was to illuminate a few reflections from 1 or 2 subgrains and follow their evolution during heating. Preliminary results show that significant changes occurred within the sample volume, for which, surface effects can be excluded. Results show a number of processes, including: i) change in subgrain boundary misorientation angle and ii) subgrain subdivision into areas of similar lattice orientation with new subgrain boundary formation. These results demonstrate that 3DXRD coupled with in-situ heating is a successful non-destructive technique for examining real-time post-deformational annealing in strongly deformed crystalline materials with complicated microstructures.