Small quantity but large effect — How minor phases control strain localization in upper mantle shear zones

Low viscosity domains such as localized shear zones exert an important control on the geodynamics of the uppermost mantle. Grain size reduction and subsequent strain localization related to a switch from dislocation to diffusion creep is one mechanism to form low viscosity domains. To sustain strain localization, the grain size of mantle minerals needs to be kept small over geological timescales. One way to keep olivine grain sizes small is by pinning of mobile grain boundaries during grain growth by other minerals (second phases). Detailed microstructural studies based on natural samples from three shear zones formed at different geodynamic settings, allowed the derivation of the olivine grain-size dependence on the second-phase content. The polymineralic olivine grain-size evolution with increasing strain is similar in the three shear zones. If the second phases are to pin the mobile olivine grain boundary the phases need to be well mixed before grain growth. We suggest that melt-rock and metamorphic reactions are crucial for the initial phase mixing in mantle rocks. With ongoing deformation and increasing strain, grain boundary sliding combined with mass transfer processes and nucleation of grains promotes phase mixing resulting in fine-grained polymineralic mixtures that deform by diffusion creep. Strain localization due to the presence of volumetrically minor minerals in polymineralic mantle rocks is only important at high strain deformation (ultramylonites) at low temperatures (<~800°C). At smaller strain and stress conditions and/or higher temperatures other parameters like overall energy available to deform a given rock volume, the inheritance of mechanical anisotropies or the presence of water or melts needs to be considered to explain strain localization in the upper mantle.

Terrestrial Radar Interferometric Measurement of Hillslope Deformation and Atmospheric Disturbances in the Illgraben Debris-Flow Catchment, Switzerland

We describe a method for rapid identification and precise quantification of slope deformation using a portable radar interferometer. A rockslide with creep-like behavior was identified in the rugged and inaccessible headwaters of the Illgraben debris-flow catchment, located in the Central Swiss Alps. The estimated volume of the moving rock mass was approximately 0.5 x 10(6) m(3) with a maximum daily (3-D) displacement rate of 3 mm. Fast scene acquisition in the order of 6 s/scene led to uniquely precise mapping of spatial and temporal variability of atmospheric phase delay. Observations led to a simple qualitative model for prediction of atmospheric disturbances using a simple model for solar radiation, which can be used for advanced campaign planning for short observation periods (hours to days).

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.

Redistribution of particles across the nucleus of comet 67P/Churyumov-Gerasimenko

Context. We present an investigation of the surface properties of areas on the nucleus of comet 67P/Churyumov-Gerasimenko. Aims. We aim to show that transport of material from one part of the cometary nucleus to another is a significant mechanism that influences the appearance of the nucleus and the surface thermal properties. Methods. We used data from the OSIRIS imaging system onboard the Rosetta spacecraft to identify surface features on the nucleus that can be produced by various transport mechanisms. We used simple calculations based on previous works to establish the plausibility of dust transport from one part of the nucleus to another. Results. We show by observation and modeling that "airfall" as a consequence of non-escaping large particles emitted from the neck region of the nucleus is a plausible explanation for the smooth thin deposits in the northern hemisphere of the nucleus. The consequences are also discussed. We also present observations of aeolian ripples and ventifacts. We show by numerical modeling that a type of saltation is plausible even under the rarified gas densities seen at the surface of the nucleus. However, interparticle cohesive forces present difficulties for this model, and an alternative mechanism for the initiation of reptation and creep may result from the airfall mechanism. The requirements on gas density and other parameters of this alternative make it a more attractive explanation for the observations. The uncertainties and implications are discussed.

Crevassing and Calving of Glacial Ice

Calving of ice is a relatively new area of research in the still young field of glaciology. In the short time that calving has been studied, it has been mainly treated as an afterthought, with the predominant mode of thinking being that it will happen so to concern oneself with why is not important. Many studies dealt with observations of calving front positions over time vs. ice velocity in an attempt to quantify the calving rate as the difference between the two, while others have attempted to deduce some empirical relationship between calving rate and variables such as water depth or temperature. This study instead addresses the question of why, where, and when ice will first become crevassed, which is an obviously necessary condition for a later calving event to occur. Previous work examining the causes of calving used ideas put forth from a variety of fields, including civil engineering, materials science, and results from basic physics and mechanics. These theories are re-examined here and presented as part of a larger whole. Important results from the field of fracture mechanics are utilized frequently, and these results can be used as a predictor of ice behavior and intrinsic properties of ice, as well as properties like back stresses induced by local pinning points and resistive shears along glacial ice boundaries. A theory of fracture for a material experiencing creep is also presented with applications to ice shelves and crevasse penetration. Finally, a speculative theory regarding large scale iceberg formation is presented. It is meant mainly as an impetus to further discussion on the topic, with the hope that a model relating crevasse geometries to flow parameters can result in crevasse spacings that could produce the tabular icebergs which are so newsworthy. The primary focus of this thesis is to move away from the "after the fact" studies that are so common in calving research, and instead devote energy to determining what creates the conditions that drive the calving of ice in the first place.

Ultimate and failure strengths of normal muds and serpentinite muds of ODP Leg 125 samples

Serpentinite seamounts in the Mariana forearc have been explained as diapirs rising from the Benioff zone. This hypothesis predicts that the serpentinites should have low strengths as well as low densities relative to the surrounding rocks. Drilling during Leg 125 showed that the materials forming Conical Seamount in the Mariana forearc and Torishima Forearc Seamount in the Izu-Bonin forearc are water-charged serpentinite muds of density <2 g/cm**3. Wykeham-Farrance torsion-vane tests showed that they are plastic solids with a rheology that bears many similarities to the idealized Cam clay soil model and is well described by critical-state soil mechanics. The serpentinite muds have ultimate strengths of 1.3 to 273.7 kPa and yield strengths of approximately 1.0 to 50 kPa. These muds thus are orders of magnitude weaker than salt and are, in fact, comparable in density and strength to common deep-sea clay muds. Such weak and low-density materials easily become diapiric. Serpentinite muds from the summit of Conical Seamount are weaker and more ductile than those on its flanks or on Torishima Forearc Seamount. Moreover, the summit muds do not contain the pronounced pure- and simple-shear fabrics that characterize those on the seamount flanks. The seamounts are morphologically similar to shield volcanoes, and anastomosing serpentinite debris flows descending from their summits are similar in map view to pahoehoe flows. These morphologic features, together with the physical properties of the muds and their similarities to other oceanic muds and the geochemistry of the entrained waters, suggest that many forearc serpentinite seamounts are gigantic (10-20 km wide, 1.5-2.0 km high) mud volcanoes that formed by the eruption of highly liquid serpentinite muds. Torishima Forearc Seamount, which is blanketed by more ìnormalî pelagic/volcaniclastic sediment, has probably been inactive since the Miocene. Conical Seamount, which seems to consist entirely of serpentinite mud and is venting fresh water of unusual chemistry from its summit, is presently active. Muds from the flanks of Conical Seamount are stronger and more brittle than those from the summit site, and muds from Torishima Forearc Seamount are stronger yet; this suggests that the serpentinite debris flows are compacted and dewatered as they mature. The shear fabrics probably result from downslope creep and flow, but may also be inherited.

Benthic foraminiferal composition at the inner wall of the Middle America Trench, Costa Rica

The sediments of Deep Sea Drilling Project Site 565 and University of Texas Marine Science Institute Cores IG-24-7-38 to -42 taken on the landward slope of the Middle America Trench exhibit characteristics of material subject to reworking during downslope mass flow. These characteristics include a generally homogeneous texture, lack of sedimentary structures, pervasive presence of a penetrative scaly fabric, and presence of transported benthic foraminifers. Although these features occur throughout the sediments examined, trends in bulk density, porosity, and water content, and abrupt shifts in these index physical properties and in sediment magnetic properties at Site 565 indicate that downslope sediment creep is presently most active in the upper 45 to 50 m of sediment. It cannot be determined whether progressive dewatering of sediment has brought the material at this depth to a plastic limit at which sediment can no longer flow (thus resulting in its accretion to the underlying sediments) or whether this depth represents a surface along which slumping has occurred. We suspect both are true in part, that is, that mass movements and downslope reworking accumulate sediments in a mobile layer of material that is self-limiting in thickness.

Measurements of velocity-dependent and slip-dependent frictional strength in laboratory friction experiments on natural samples from IODP Expedition 316 and ODP Hole 190-1174B

Slowslip forms part of the spectrum of fault behaviour between stable creep and destructive earthquakes. Slow slip occurs near the boundaries of large earthquake rupture zones and may sometimes trigger fast earthquakes. It is thought to occur in faults comprised of rocks that strengthen under fast slip rates, preventing rupture as a normal earthquake, or on faults that have elevated pore-fluid pressures. However, the processes that control slow rupture and the relationship between slow and normal earthquakes are enigmatic. Here we use laboratory experiments to simulate faulting in natural rock samples taken from shallow parts of the Nankai subduction zone, Japan, where very low-frequency earthquakes - a form of slow slip - have been observed.We find that the fault rocks exhibit decreasing strength over millimetre-scale slip distances rather than weakening due to increasing velocity. However, the sizes of the slip nucleation patches in our laboratory simulations are similar to those expected for the very lowfrequency earthquakes observed in Nankai. We therefore suggest that this type of fault-weakening behaviour may generate slow earthquakes. Owing to the similarity between the expected behaviour of slow earthquakes based on our data, and that of normal earthquakes during nucleation, we suggest that some types of slow slip may represent prematurely arrested earthquakes.

Mean annual surface displacement of the Laurichard rock glacier, French Alps

The Laurichard active rock glacier is the permafrost-related landform with the longest record of monitoring in France, including an annual geodetic survey, repeated geoelectrical campaigns from 1979 onwards and continuous recording of ground temperature since 2003. These data were used to examine changes in creep rates and internal structure from 1986 to 2006. The control that climatic variables exert on rock glacier kinematics was investigated over three time scales. Between the 1980s and the early 2000s, the main observed changes were a general increase in surface velocity and a decrease in internal resistivity. At a multi-year scale, the high correlation between surface movement and snow thickness in the preceding December appears to confirm the importance of snow cover conditions in early winter through their influence on the ground thermal regime. A comparison of surface velocities, regional climatic datasets and ground sub-surface temperatures over six years suggests a strong relation between rock glacier deformation and ground temperature, as well as a role for liquid water due to melt of thick snow cover. Finally, unusual surface lowering that accompanied peak velocities in 2004 may be due to a general thaw of the top of the permafrost, probably caused both by two successive snowy winters and by high energy inputs during the warm summer of 2003.

Results from Purgatory Conglomerate investigations

Quartz Crystallographic Preferred Orientation (CPO) patterns are most commonly a result of deformation by dislocation creep. We investigated whether Dissolution-Precipitation Creep (DPC) a process that occur at lower differential stresses and temperatures, may result in CPO in quartz. Within the Purgatory Conglomerate, DPC led to quartz dissolution along cobble surfaces perpendicular to the shortening direction, and quartz precipitation in overgrowths at the ends of the cobbles (strain shadows), parallel to the maximum extension direction. The Purgatory Conglomerate is part of the SE Narragansett basin where strain intensity increases from west to east and is associated with top-to-the-west transport and folding during the Alleghanian orogeny. Quartz c-axis orientations as revealed by Electron Backscatter Diffraction (EBSD) methods, were random in all analyzed domains within the cobbles and strain shadows irrespective of the intensity of strain or metamorphic grade of the sample. Quartz dissolution probably occurred exclusively along the cobbles' margins, leaving the remaining grains unaffected by DPC. The fact that quartz precipitated in random orientations may indicate that the strain shadows were regions of little or no differential stress.

A numerical discrete element study on the influence of an embedded viscoelastic-plastic layer at different viscosity values

The major aim of this study was to examine the influence of an embedded viscoelastic-plastic layer at different viscosity values on accretionary wedges at subduction zones. To quantify the effects of the layer viscosity, we analysed the wedge geometry, accretion mode, thrust systems and mass transport pattern. Therefore, we developed a numerical 2D 'sandbox' model utilising the Discrete Element Method. Starting with a simple pure Mohr Coulomb sequence, we added an embedded viscoelastic-plastic layer within the brittle, undeformed 'sediment' package. This layer followed Burger's rheology, which simulates the creep behaviour of natural rocks, such as evaporites. This layer got thrusted and folded during the subduction process. The testing of different bulk viscosity values, from 1 × 10**13 to 1 × 10**14 (Pa s), revealed a certain range where an active detachment evolved within the viscoelastic-plastic layer that decoupled the over- and the underlying brittle strata. This mid-level detachment caused the evolution of a frontally accreted wedge above it and a long underthrusted and subsequently basally accreted sequence beneath it. Both sequences were characterised by specific mass transport patterns depending on the used viscosity value. With decreasing bulk viscosities, thrust systems above this weak mid-level detachment became increasingly symmetrical and the particle uplift was reduced, as would be expected for a salt controlled forearc in nature. Simultaneously, antiformal stacking was favoured over hinterland dipping in the lower brittle layer and overturning of the uplifted material increased. Hence, we validated that the viscosity of an embedded detachment strongly influences the whole wedge mechanics, both the respective lower slope and the upper slope duplex, shown by e.g. the mass transport pattern.

Lithospheric thickness of Pangea

Surface wave tomography, using the fundamental Rayleigh wave velocities and those of higher modes between 1 and 4 and periods between 50 and 160 s, is used to image structures with a horizontal resolution of ~250 km and a vertical resolution of ~50 km to depths of ~300 km in the mantle. A new model, PM_v2_2012, obtained from 3×10**6 seismograms, agrees well with earlier lower resolution models. It is combined with temperature estimates from oceanic plate models and with pressure and temperature estimates from the mineral compositions of garnet peridotite nodules to generate a number of estimates of SV(P,T) based on geophysical and petrological observations alone. These are then used to estimate the unrelaxed shear modulus and its derivatives with respect to pressure and temperature, which agree reasonably with values from laboratory experiments. At high temperatures relaxation occurs, causing the shear wave velocity to depend on frequency. This behaviour is parameterised using a viscosity to obtain a Maxwell relaxation time. The relaxation behaviour is described using a dimensionless frequency, which depends on an activation energy E and volume Va. The values of E and Va obtained from the geophysical models agree with those from laboratory experiments on high temperature creep. The resulting expressions are then used to determine the lithospheric thickness from the shear wave velocity variations. The resolution is improved by about a factor of two with respect to earlier models, and clearly resolves the thick lithosphere beneath active intracontinental belts that are now being shortened. The same expressions allow the three dimensional variations of the shear wave attenuation and viscosity to be estimated.