Laura A. Stern

Mantle phase changes and deep-earthquake faulting in subducting lithosphere

Inclined zones of earthquakes are the primary expression of lithosphere subduction. A distinct deep population of subduction-zone earthquakes occurs at depths of 350 to 690 kilometers. At those depths ordinary brittle fracture and frictional sliding, the faulting processes of shallow earthquakes, are not expected. A fresh understanding of these deep earthquakes comes from developments in several areas of experimental and theoretical geophysics, including the discovery and characterization of transformational faulting, a shear instability connected with localized phase transformations under nonhydrostatic stress. These developments support the hypothesis that deep earthquakes represent transformational faulting in a wedge of olivine-rich peridotite that is likely to persist metastably in coldest plate interiors to depths as great as 690 km. Predictions based on this deep structure of mantle phase changes are consistent with the global depth distribution of deep earthquakes, the maximum depths of earthquakes in individual subductions zones, and key source characteristics of deep events.

Peculiarities of Methane Clathrate Hydrate Formation and Solid-State Deformation, Including Possible Superheating of Water Ice

Slow, constant-volume heating of water ice plus methane gas mixtures forms methane clathrate hydrate by a progressive reaction that occurs at the nascent ice/liquid water interface. As this reaction proceeds, the rate of melting of metastable water ice may be suppressed to allow short-lived superheating of ice to at least 276 kelvin. Plastic flow properties measured on clathrate test specimens are significantly different from those of water ice; under nonhydrostatic stress, methane clathrate undergoes extensive strain hardening and a process of solid-state disproportionation or exsolution at conditions well within its conventional hydrostatic stability field.

Effects of dispersed particulates on the rheology of water ice at planetary conditions

We have investigated the effects of initial grain size and hard particulate impurities on the transient and steady state flow of water ice I at laboratory conditions selected to provide more quantitative constraints on the thermomechanical evolution of the giant icy moons of the outer solar system. Our samples were molded with particulate volume fractions, ϕ, of 0.001 to 0.56 and particle sizes of 1 to 150 μm. Deformation experiments were conducted at constant shortening rates of 4.4 × 10−7 to 4.9 × 10−4 s−1 at pressures of 50 and 100 MPa and temperatures 77 to 223 K. For the pure ice samples, initial grain sizes were 0.2–0.6 mm, 0.75–1.75 mm, and 1.25–2.5 mm. Stress-strain curves of pure ice I under these conditions display a strength maximum σu at plastic strains e ≤ 0.01 after initial yield, followed by strain softening and achievement of steady state levels of stress, σss, at e = 0.1 to 0.2. Finer starting grain size in pure ice generally raises the level of σu. Petrography indicates that the initial transient flow behavior is associated with the nucleation and growth of recrystallized ice grains and the approach to σss evidently corresponds to the development of a steady state grain texture. Effects of particulate concentrations ϕ < 0.1 are slight. At these concentrations, a small but significant reduction in σu with respect to that for pure water ice occurs. Mixed-phase ice with ϕ ≥ 0.1 is significantly stronger than pure ice; the strength of samples with ϕ = 0.56 approaches that of dry confined sand. The magnitude of the strengthening effect is far greater than expected from homogeneous strain-rate enhancement in the ice fraction or from pinning of dislocations (Orowan hardening). This result suggests viscous drag occurs in the ice as it flows around the hard particulates. Mixed-phase ice is also tougher than pure ice, extending the range of bulk plastic deformation versus faulting to lower temperatures and higher strain rates. The high-pressure phase ice II formed in ϕ = 0.56 mixed-phase ice during deformation at high stresses. Bulk planetary compositions of ice + rock (ϕ = 0.4–0.5) are roughly 2 orders of magnitude more viscous than pure ice, promoting the likelihood of thermal instability inside giant icy moons and possibly explaining the retention of crater topography on icy planetary surfaces.

Creep of water ices at planetary conditions: A compilation

Many constitutive laws for the flow of ice have been published since the advent of the Voyager explorations of the outer solar system. Conflicting data have occasionally come from different laboratories, and refinement of experimental techniques has led to the publication of laws that supersede earlier ones. In addition, there are unpublished data from ongoing research that also amend the constitutive laws. Here we compile the most current laboratory-derived flow laws for water ice phases I, II, III, V, and VI, and ice I mixtures with hard particulates. The rheology of interest is mainly that of steady state, and the conditions reviewed are the pressures and temperatures applicable to the surfaces and interiors of icy moons of the outer solar system. Advances in grain-size-dependent creep in ices I and II as well as in phase transformations and metastability under differential stress are also included in this compilation. At laboratory strain rates the several ice polymorphs are rheologically distinct in terms of their stress, temperature, and pressure dependencies but, with the exception of ice III, have fairly similar strengths. Hard particulates strengthen ice I significantly only at high particulate volume fractions. Ice III has the potential for significantly affecting mantle dynamics because it is much weaker than the other polymorphs and its region of stability, which may extend metastably well into what is nominally the ice II field, is located near likely geotherms of large icy moons.

Rheology of ice I at low stress and elevated confining pressure

Triaxial compression testing of pure, polycrystalline water ice I at conditions relevant to planetary interiors and near-surface environments (differential stresses 0.45 to 10 MPa, temperatures 200 to 250 K, confining pressure 50 MPa) reveals that a complex variety of rheologies and grain structures may exist for ice and that rheology of ice appears to depend strongly on the grain structures. The creep of polycrystalline ice I with average grain size of 0.25 mm and larger is consistent with previously published dislocation creep laws, which are now extended to strain rates as low as 2×10−8 s−1. When ice I is reduced to very fine and uniform grain size by rapid pressure release from the ice II stability field, the rheology changes dramatically. At 200 and 220 K the rheology matches the grain-size-sensitive rheology measured by Goldsby and Kohlstedt [1997, this issue] at 1 atm. This finding dispels concerns that the Goldsby and Kohlstedt results were influenced by mechanisms such as microfracturing and cavitation, processes not expected to operate at elevated pressures in planetary interiors. At 233 K and above, grain growth causes the fine-grained ice to become more creep resistant. Scanning electron microscopy investigation of some of these deformed samples shows that grains have markedly coarsened and the strain hardening can be modeled by normal grain growth and the Goldsby and Kohlstedt rheology. Several samples also displayed very heterogeneous grain sizes and high aspect ratio grain shapes. Grain-size-sensitive creep and dislocation creep coincidentally contribute roughly equal amounts of strain rate at conditions of stress, temperature, and grain size that are typical of terrestrial and planetary settings, so modeling ice dynamics in these settings must include both mechanisms.

Scanning Electron Microscopy investigations of laboratory-grown gas clathrate hydrates formed from melting ice, and comparison to natural hydrates

Abstract Scanning electron microscopy (SEM) was used to investigate grain texture and pore structure development within various compositions of pure sI and sII gas hydrates synthesized in the laboratory, as well as in natural samples retrieved from marine (Gulf of Mexico) and permafrost (NW Canada) settings. Several samples of methane hydrate were also quenched after various extents of partial reaction for assessment of mid-synthesis textural progression. All laboratory-synthesized hydrates were grown under relatively high-temperature and high-pressure conditions from rounded ice grains with geometrically simple pore shapes, yet all resulting samples displayed extensive recrystallization with complex pore geometry. Growth fronts of mesoporous methane hydrate advancing into dense ice reactant were prevalent in those samples quenched after limited reaction below and at the ice point. As temperatures transgress the ice point, grain surfaces continue to develop a discrete “rind” of hydrate, typically 5 to 30 μm thick. The cores then commonly melt, with rind microfracturing allowing migration of the melt to adjacent grain boundaries where it also forms hydrate. As the reaction continues under progressively warmer conditions, the hydrate product anneals to form dense and relatively pore-free regions of hydrate grains, in which grain size is typically several tens of micrometers. The prevalence of hollow, spheroidal shells of hydrate, coupled with extensive redistribution of reactant and product phases throughout reaction, implies that a diffusion-controlled shrinking-core model is an inappropriate description of sustained hydrate growth from melting ice. Completion of reaction at peak synthesis conditions then produces exceptional faceting and euhedral crystal growth along exposed pore walls. Further recrystallization or regrowth can then accompany even short-term exposure of synthetic hydrates to natural ocean-floor conditions, such that the final textures may closely mimic those observed in natural samples of marine origin. Of particular note, both the mesoporous and highly faceted textures seen at different stages during synthetic hydrate growth were notably absent from all examined hydrates recovered from a natural marine-environment setting.

Flow of ices in the ammonia‐water system

We have fabricated in the laboratory and subsequently deformed crystalline hydrates and partial melts of the water-rich end of the NH3-H2O system, with the aim of improving our understanding of physical processes occurring in icy moons of the outer solar system. Deformation experiments were carried out at constant strain rate. The range of experimental variables was strain rate 3.5 × 10−7 < ϵe < 3.5 × 10−4 s−1, temperature 132 < T < 220 K, pressure 50 ≤ P ≤ 100 MPa, and mole fraction NH3 0 ≤ xNH3 ≤ 0.295. Phase relationships in the NH3-H2O system indicate that water ice and ammonia dihydrate, NH3 · 2H2O, are the stable phases under our experimental conditions. X ray diffraction of our samples usually revealed these as the dominant phases, but we have also observed an amorphous phase (in unpressurized samples only) and occasionally significant ammonia monohydrate, NH3 · H2O. The onset of partial melting at the peritectic temperature at about 176 K appeared as a sharp transition in strength observed in samples of xNH3 = 0.15 and 0.295. In samples of xNH3 = 0.05 and 0.01, the effect of melt was less pronounced. For any given water ice + dihydrate alloy in the subsolidus region, we observed one rheological law over the entire temperature range from 176 K to about 140 K. Below 140 K, a shear instability similar to that occurring in pure water ice under the same conditions limited our ability to measure ductile flow. The rheological laws for the several alloys vary systematically from that of pure ice to that of dihydrate. Pure dihydrate is about 4 orders of magnitude less viscous than water ice just below the peritectic temperature, but because of a very pronounced temperature dependence in dihydrate (100 kJ/mol versus 43 kJ/mol for water ice) the viscosity of dihydrate equals or exceeds that of water ice at T < 140 K. The large variation in viscosity of dihydrate with relatively small changes in temperature may be helpful in explaining the rich variety of tectonic and volcanic features seen on the surfaces of icy moons in the outer solar system.

Laboratory measurements of compressional and shear wave speeds through methane hydrate

Abstract: Simultaneous measurements of compressional and shear wave speeds through polycrystalline methane hydrate have been made. Methane hydrate, grown directly in a wave speed measurement chamber, was uniaxially compacted to a final porosity below 2%. At 277 K, the compacted material had a compressional wave speed of 3,650 ± 50 m/s. The shear wave speed, measured simultaneously, was 1,890 ± 30 m/s. From these wave speed measurements, we derive Vp/Vs, Poissons ratio, bulk, shear, and Youngs moduli.

Grain‐size‐induced weakening of H2O ices I and II and associated anisotropic recrystallization

Grain-size-dependent flow mechanisms tend to be favored over dislocation creep at low differential stresses and can potentially influence the rheology of low-stress, low-strain rate environments such as those of planetary interiors. We experimentally investigated the effect of reduced grain size on the solid-state flow of water ice I, a principal component of the asthenospheres of many icy moons of the outer solar system, using techniques new to studies of this deformation regime. We fabricated fully dense ice samples of approximate grain size 2±1 μm by transforming “standard” ice I samples of 250±50 μm grain size to the higher-pressure phase ice II, deforming them in the ice II field, and then rapidly releasing the pressure deep into the ice I stability field. At T≤200 K, slow growth and rapid nucleation of ice I combine to produce a fine grain size. Constant-strain rate deformation tests conducted on these samples show that deformation rates are less stress sensitive than for standard ice and that the fine-grained material is markedly weaker than standard ice, particularly during the transient approach to steady state deformation. Scanning electron microscope examination of the deformed fine-grained ice samples revealed an unusual microstructure dominated by platelike grains that grew normal to the compression direction, with c axes preferentially oriented parallel to compression. In samples tested at T≥220 K the elongation of the grains is so pronounced that the samples appear finely banded, with aspect ratios of grains approaching 50:1. The anisotropic growth of these crystallographically oriented neoblasts likely contributes to progressive work hardening observed during the transient stage of deformation. We have also documented remarkably similar microstructural development and weak mechanical behavior in fine-grained ice samples partially transformed and deformed in the ice II field.

Laboratory synthesis of pure methane hydrate suitable for measurement of physical properties and decomposition behavior

Gas hydrates are an intriguing class of nonstoichiometric compounds that have significant commercial and scientific applications both as an energy resource and as a manufactured material. The last half-century has witnessed a marked escalation in the scope of experimental research on gas hydrates, particularly directed towards the determination of their phase equilibria, formation kinetics, crystallographic and structural properties, transport and thermal properties, effects of inhibitors, and a number of related geochemical topics.