Geology topics

Stephen H. Kirby

Publications and source records attributed to Stephen H. Kirby.

51 records · Page 3Linked to original sources

USGS2005

The calcite → aragonite transformation in low-Mg marble: Equilibrium relations, transformations mechanisms, and rates

Experimental transformation of a rather pure natural calcite marble to aragonite marble did not proceed via the expected straightforward polymorphic replacement. Instead, the small amount of Mg in the starting material (0.36 wt %) was excluded from the growing aragonite and diffused preferentially into the remaining calcite grains, producing Mg-rich calcite rods that persisted as relicts. Nucleation of aragonite occurred exclusively on grain boundaries, with aragonite [001] oriented subparallel to calcite [0001]. The aragonite crystals preferentially consumed the calcite crystal on which they nucleated, and the reaction fronts developed preferentially along the {010} and {110} planes of aragonite. Each aragonite neoblast that grew was nearly free of Mg (typically <0.1 wt %). The excess Mg was taken up by the calcite grains in between, stabilizing them and causing a few volume percent rodlike relicts of Mg-enriched calcite (up to 10 wt % MgO) to be left behind by the advancing reaction front. The aragonite growth rates are approximately linear and range from &sim;3 × 10 −11 m s −1 at 600°C to &sim;9 × 10 −9 m s −1 at 850°C, with an apparent activation enthalpy of 166 ± 91 kJ mol −1 . This reaction mechanism and the resultant texture are akin to cellular precipitation reactions in metals. Similar transformation textures have been reported from high-Mg marbles in Japan and China that disproportionated to low-Mg calcite and dolomite.

Journal of Geophysical Research B: Solid Earth↗

USGS2003

Introduction to physical properties and elasticity models

Estimating the in situ methane hydrate volume from seismic surveys requires knowledge of the rock physics relations between wave speeds and elastic moduli in hydrate/sediment mixtures. The elastic moduli of hydrate/sediment mixtures depend on the elastic properties of the individual sedimentary particles and the manner in which they are arranged. In this chapter, we present some rock physics data currently available from literature. The unreferenced values in Table I were not measured directly, but were derived from other values in Tables I and II using standard relationships between elastic properties for homogeneous, isotropic material. These derivations allow us to extend the list of physical property estimates, but at the expense of introducing uncertainties due to combining property values measured under different physical conditions. This is most apparent in the case of structure II (sII) hydrate for which very few physical properties have been measured under identical conditions.

Book chapter↗

USGS2002

The Cascadia Subduction Zone and related subduction systems: Seismic structure, intraslab earthquakes and processes, and earthquake hazards

The following report is the principal product of an international workshop titled “Intraslab Earthquakes in the Cascadia Subduction System: Science and Hazards” and was sponsored by the U.S. Geological Survey, the Geological Survey of Canada and the University of Victoria. This meeting was held at the University of Victoria’s Dunsmuir Lodge, Vancouver Island, British Columbia, Canada on September 18–21, 2000 and brought 46 participants from the U.S., Canada, Latin America and Japan. This gathering was organized to bring together active research investigators in the science of subduction and intraslab earthquake hazards. Special emphasis was given to “warm-slab” subduction systems, i.e., those systems involving young oceanic lithosphere subducting at moderate to slow rates, such as the Cascadia system in the U.S. and Canada, and the Nankai system in Japan. All the speakers and poster presenters provided abstracts of their presentations that were a made available in an abstract volume at the workshop. Most of the authors subsequently provided full articles or extended abstracts for this volume on the topics that they discussed at the workshop. Where updated versions were not provided, the original workshop abstracts have been included. By organizing this workshop and assembling this volume, our aim is to provide a global perspective on the science of warm-slab subduction, to thereby advance our understanding of internal slab processes and to use this understanding to improve appraisals of the hazards associated with large intraslab earthquakes in the Cascadia system. These events have been the most frequent and damaging earthquakes in western Washington State over the last century. As if to underscore this fact, just six months after this workshop was held, the magnitude 6.8 Nisqually earthquake occurred on February 28th, 2001 at a depth of about 55 km in the Juan de Fuca slab beneath the southern Puget Sound region of western Washington. The Governor’s Office of the State of Washington estimated damage at more than US$2 billion, making it among the costliest earthquakes in U.S. history.

California, Oregon, Washington↗

USGS2001

The radiated seismic energy and apparent stress of interplate and intraplate earthquakes at subduction zone environments; implications for seismic hazard estimation

The radiated seismic energies (ES) of 980 shallow subduction-zone earthquakes with magnitudes ? 5.8 are used to examine global patterns of energy release and apparent stress. In contrast to traditional methods which have relied upon empirical formulas, these energies are computed through direct spectral analysis of broadband seismic waveforms. Energy gives a physically different measure of earthquake size than moment. Moment, being derived from the low-frequency asymptote of the displacement spectra, is related to the final static displacement. Thus, moment is crucial to the long-term tectonic implication of an earthquake. In contrast, energy, being derived from the velocity power spectra, is more a measure of seismic potential for damage to anthropogenic structures. There is considerable scatter in the plot of ES-M0 for worldwide earthquakes. For any given M0, the ES can vary by as much as an order of magnitude about the mean regression line. The global variation between ES and M0, while large, is not random. When subsets of ES-M0 are plotted as a function of seismic region, tectonic setting and faulting type, the scatter in data is often substantially reduced. There are two profound implications for the estimation of seismic and tsunamic hazard. First, it is now feasible to characterize the apparent stress for particular regions. Second, a given M0 does not have a unique ES. This means that M0 alone is not sufficient to describe all aspects of an earthquake. In particular, we have found examples of interplate thrust-faulting earthquakes and intraslab normal-faulting earthquakes occurring in the same epicentral region with vastly different macroseismic effects. Despite the gross macroseismic disparities, the MW?s in these examples were identical. However, the Me?s (energy magnitudes) successfully distinguished the earthquakes that were more damaging.

Open-File Report↗

USGS2000

Taking the temperature of slabs

No abstract available.

Nature↗

USGS1999

Steady-state flow of solid CO2: Preliminary results

To help answer the question of how much solid CO 2 exists in the Martian south polar cap, we performed a series of laboratory triaxial deformation experiments at constant displacement rate in compression on jacketed cylinders of pure, polycrystalline CO 2 . Test conditions were temperatures 150 < T < 190 K, hydrostatic confining pressures 5≤ P ≤40 MPa, and strain rates 4.5×10 −8 ≤ ε ≤4.3×10 −4 s −1 . Most of the measurements follow a constitutive law of the form ε = Aσ n exp(−Q/RT), where σ is the applied differential stress, R is the gas constant, and the other constants have values as follows: A = 10 3 86 MPa −n s −1 , n = 5.6, and Q = 33 kJ/mol. Solid CO 2 is markedly weaker than water ice. Our results suggest that the south polar cap on Mars is unlikely to be predominately solid CO 2 , because the elevation and estimated age of the cap is difficult to reconcile with the very weak rheology of the material.

Geophysical Research Letters↗

USGS1996

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.

Science↗

USGS1995

Frequency-moment distribution of deep earthquakes; Implications for the seismogenic zone at the bottom of slabs

We present a systematic investigation of the variation with depth of the frequency of earthquake occurrence vs. seismic moment based on 16 years of Harvard Centroid Moment Tensor (CMT) solutions. We analyze depth variations of earthquake size distribution in terms of variations in the absolute value of the slope of the regression of the logarithm of the population vs. seismic moment, a quantity known as the β parameter. The shallowest earthquakes (0–50 km depth) exhibit a well-defined and robust size distribution regime characterized by a discontinuous increase in β with increasing moment. Others have shown that this increase probably represents the effects of a physical limit in the dimensions of the area of seismogenic slip of shallow earthquake sources. The population of deep earthquakes in the depth interval 500–600 km shows two markedly different distributions. The deep earthquakes in the Tonga region feature an initially high β value (0.92) at small moments and a lower β value (0.41) at high moments. In contrast, the size distribution of non-Tonga deep events shows the reverse of those changes (β = 0.41 at low moment and β = 1.17 at higher moment). To help explain these observations, we propose a model of deep seismogenesis that assumes three-dimensional earthquake source regions that vary principally in their transverse dimensions. The two-β segment behavior in the Tonga region and other subduction zones is thought to represent, in part, constraints owing to the threshold of completeness of the CMT catalog and to its short time interval of sampling. We interpret the differences between Tonga and other deep Wadati-Benioff zones as being a consequence of Tonga's markedly higher subduction rate and, hence, its colder thermal structure and presumably thicker region of seismogenesis. We interpret the critical moments at which β values change in terms of variations in the transverse thickness of deep seismogenic zones and estimate that it is about 11 km for the Tonga region and about 3 km for other zones at depths of 500–600 km. These results are generally consistent with deep earthquakes being restricted to wedge-shaped regions of peridotite persisting metastably to as deep as 700 km in old, rapidly descending and hence cold slabs. Failure is thought to occur in metastable peridotite by transformational faulting. Great deep earthquakes present special challenges to any theory of deep earthquakes based on slab thermal structure. For example, a continuing question is how such large events can fit in a thermally controlled seismogenic zone that is diminishing in its transverse dimensions with increasing depth. The very concept of a scale-invariant earthquake size distribution may be inappropriate for these rare events and the unusual settings in which they are found.

Physics of the Earth and Planetary Interiors↗

USGS1991

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.

Science↗

USGS1987

Inelastic properties of several high pressure crystalline phases of H2O: Ices II, III, and V

We have performed deformation experiments on cylinders of polycrystalline H 2 O at temperatures from 178 to 257 K at pressures to 500 MPa in the stability fields of ices II, III, and V. Ice II is the strongest of the phases, having a strength under laboratory conditions roughly comparable to that of ice I h . Ice V is somewhat weaker than ice II. Ice III is extremely weak and over geologic times must behave essentially as a liquid bounded below by ice V and above by ice II or I h . Phase relationships are complicated by a number of phase metastabilities, the most important of which is the existence of ice III in the ice II field for extended periods of time. Even under deformation at temperatures as low as 211 K (over 30 K below the ice III field), the transformations from III to II can not be made to happen in the laboratory.

Journal de Physique Colloques↗

USGS1987

Ionic conductivity of quartz: DC time dependence and transition in charge carriers

The time dependence of DC electrical conductivity in the c-axis direction of quartz can be accounted for by a transition in charge carriers from interstitial alkali impurities to interstitial H. The diffusive transport rates of Li, Na, and K are rapid parallel to c and have been shown to be responsible for the highly anisotropic electrical conductivity measured at short times. With increasing time, however, conductivities parallel to c decrease progressively to values that are roughly equal to those measured perpendicular to c. Comparison of these ultimate, nearly isotropic conductivities with those derived from recent measurements of H diffusion parallel and perpendicular to c suggests that H interstitials are the principal charge carriers at long times. The transient decrease in conductivities parallel to c is interpreted to result from depletion of initial alkali impurities, whereas the steady-state conductivities measured at long times may be sustained by the steady supply of H by the dissociation of atmospheric water vapor. The mobility of H along the c axis is anomalously low and at variance with the trend of increasing mobility with decreasing ionic radius exhibited by Cs, Rb, K, Na, and Li. Although the elastic lattice distortions required for H transport are insignificant in comparison with those required by the larger alkali impurities, the strong association of H interstitials with Al substitutions for Si may be responsible for the relatively low H mobilities.

American Mineralogist↗

USGS1985

Rock mechanics observations pertinent to the rheology of the continental lithosphere and the localization of strain along shear zones

Emphasized in this paper are the deformation processes and rheologies of rocks at high temperatures and high effective pressures, conditions that are presumably appropriate to the lower crust and upper mantle in continental collision zones. Much recent progress has been made in understanding the flexure of the oceanic lithosphere using rock-mechanics-based yield criteria for the inelastic deformations at the top and base. At mid-plate depths, stresses are likely to be supported elastically because bending strains and elastic stresses are low. The collisional tectonic regime, however, is far more complex because very large permanent strains are sustained at mid-plate depths and this requires us to include the broad transition between brittle and ductile flow. Moreover, important changes in the ductile flow mechanisms occur at the intermediate temperatures found at mid-plate depths. Two specific contributions of laboratory rock rheology research are considered in this paper. First, the high-temperature steady-state flow mechanisms and rheology of mafic and ultramafic rocks are reviewed with special emphasis on olivine and crystalline rocks. Rock strength decreases very markedly with increases in temperature and it is the onset of flow by high temperature ductile mechanisms that defines the base of the lithosphere. The thickness of the continental lithosphere can therefore be defined by the depth to a particular isotherm Tc above which (at geologic strain rates) the high-temperature ductile strength falls below some arbitrary strength isobar (e.g., 100 MPa). For olivine Tc is about 700??-800??C but for other crustal silicates, Tc may be as low as 400??-600??C, suggesting that substantial decoupling may take place within thick continental crust and that strength may increase with depth at the Moho, as suggested by a number of workers on independent grounds. Put another way, the Moho is a rheological discontinuity. A second class of laboratory observations pertains to the general phenomenon of ductile faulting in which ductile strains are localized into shear zones. Ductile faults have been produced in experiments of five different rock types and is generally expressed as strain softening in constant-strain-rate tests or as an accelerating-creep-rate stage at constant differential stress. A number of physical mechanisms have been identified that may be responsible for ductile faulting, including the onset of dynamic recrystallization, phase changes, hydrothermal alteration and hydrolytic weakening. Microscopic evidence for these processes as well as larger-scale geological and geophysical observations suggest that ductile faulting in the middle to lower crust and upper mantle may greatly influence the distribution and magnitudes of differential stresses and the style of deformation in the overlying upper continental lithosphere.

Tectonophysics↗

USGS1977

Mechanical twinning in diopside Ca(Mg,Fe)Si2O6: Structural mechanism and associated crystal defects

iopside twins mechanically on two planes, (100) and (001), and the associated macroscopic twinning strains are identical (Raleigh and Talbot, 1967). An analysis based on crystal structural arguments predicts that both twin mechanisms involve shearing of the (100) octahedral layers (containing Ca 2+ , Mg 2+ and Fe 2+ ions) by a magnitude of c /2. Small adjustments or shuffles occur in the adjacent layers containing the [SiO 4 ] 4− tetrahedral chains. While the (100) twins are conventional with shear parallel to the composition plane, this analysis predicts that (001) twins form by a mechanism closely related to kinking. A polycrystalline diopside specimen was compressed 8% at a temperature of 400° C, a pressure of 16 kilobars, and a compressive strain rate of about 10 −4 /s. Transmission electron microscopy on this specimen has revealed four basic lamellar features: 1) (100) mechanical twin lamellae; 2) (100) glide bands containing unit dislocations; 3) (001) twin lamellae; 4) (101) lamellar features, not as yet identified. The (001) twins often contain remnant (100) lamellae of untwinned host. Twinning dislocations occur in these (100) lamellae and in the (001) twin boundaries with very high densities. Diffraction contrast experiments indicate that the twinning dislocations associated with both twin laws glide on (100) with Burgers vector b = X [001] where X is probably equal to 1/2 on the basis of the structural analysis. Parallels are drawn between mechanical twinning in clinopyroxenes and clinoamphiboles. The exclusive natural occurrence of basal twins in shock-loaded clinopyroxenes and of analogous ( 1 ¯ "> 1 ¯ 1¯ 01) twins in clinoamphiboles is given a simple explanation in terms of the relative difficulty of the “kinking” mechanism as compared to direct glide parallel to the composition plane.

Physics and Chemistry of Minerals↗

USGS1975

The orthoenstatite to clinoenstatite transformation by shearing and reversion by annealing: Mechanism and potential applications

Clinoenstatite (CE) was produced by deforming single-crystal specimens of ortho-enstatite (OE) in several different sorta of experiments. Examination with light and trans-mission electron microscopes shows that the transformation is coherent and involves a macroscopic shear on (100) [001] through an angle of 12.8±1.3 °, in good agreement with the theoretically expected value of 13.3 °, and that the transformation is accomplished by glide on (100) of partial dislocations with b = 0.83[001]. Structural analysis provides further insight into the transformation mechanism. Reversion occurs in specimens annealed under a variety of conditions, and thin lamellae of CE in unconstrained, low-strain specimens recover their original shape during transformation back to OE. Our experiments and thermodynamic estimates both suggest that the equilibrium transition temperature is raised roughly 300 ° C per kilobar of shear stress on (100) [001]. This provides the basis of a method by which it may be possible to determine the magnitude as well as the orientation of the principal stresses that produce CE in nature during deformation of enstatite-bearing rocks.

Contributions to Mineralogy and Petrology↗

USGS1973

Mechanisms of high-temperature, solid-state flow in minerals and ceramics and their bearing on the creep behavior of the mantle

The problem of applying laboratory silicate-flow data to the mantle, where conditions can be vastly different, is approached through a critical review of high-temperature flow mechanisms in ceramics and their relation to empirical flow laws. The intimate association of solid-state diffusion and high-temperature creep in pure metals is found to apply to ceramics as well. It is shown that in ceramics of moderate grain size, compared on the basis of self-diffusivity and elastic modulus, normalized creep rates compare remarkably well. This comparison is paralleled by the near universal occurrence of similar creep-induced structures, and it is thought that the derived empirical flow laws can be associated with dislocation creep. Creep data in fine-grained ceramics, on the other hand, are found to compare poorly with theories involving the stress-directed diffusion of point defects and have not been successfully correlated by self-diffusion rates. We conclude that these fine-grained materials creep primarily by a quasi-viscous grain-boundary sliding mechanism which is unlikely to predominate in the earth's deep interior. Creep predictions for the mantle reveal that under most conditions the empirical dislocation creep behavior predominates over the mechanisms involving the stress-directed diffusion of point defects. The probable role of polymorphic transformations in the transition zone is also discussed.

Tectonophysics↗