Geology topics

R.A. Harris

Publications and source records attributed to R.A. Harris.

At least 19 recordsLinked to original sources

USGS2025

Paleo-scours within the layered sulfate-bearing unit at Gale crater, Mars: Evidence for intense wind erosion

The surface of modern Mars is largely shaped by wind, but the influence of past wind activity is less well constrained. Sedimentary rocks exposed in the lower foothills of Aeolis Mons, the central mound within Gale crater, record a transition from predominantly lacustrine deposition in the Murray formation to aeolian deposition in the Mirador formation. Here, we report a series of enigmatic decameter-wide, concave-up scour-and-fill structures within the Mirador formation and discuss their formation mechanisms. Using panoramic images of stratigraphy exposed in cliff faces acquired by the Curiosity rover, we map the extent, distribution and orientation of the scour-and-fill structures and document the sedimentary facies within and surrounding these structures. The scours are grouped into two classes: (A) scours with a simple, symmetric morphology and light-toned, draping infill; and (B) scours with lateral pinching and dark-toned infill. We find that the scour-enclosing environment is composed of planar, even-in-thickness laminations with a pin-stripe pattern which we interpret as wind-ripple strata formed within an aeolian sandsheet environment. Class B contains cm-scale cross-bedding and a wing-shaped feature making this scour-and-fill structure consistent with fluvial processes. We interpret scour fill of class A as an aeolian infill due to similarities with the surrounding sandsheet strata. The broad morphologies and distribution of class A are also consistent with the geometry of blowout structures formed by localized, enhanced wind deflation. These paleo-blowout structures occur clustered within the same stratigraphic interval, which may imply that they record an interval of intensified wind activity at Gale crater.

JGR Planets↗

USGS2012

Detection of a locked zone at depth on the Parkfield, California, segment of the San Andreas fault

The Parkfield, California, segment of the San Andreas fault is transitional in character between the creeping segment of the fault to the northwest and the locked Carrizo Plain segment to the southeast. The rate of shallow fault slip decreases from 25–30 mm/yr northwest of the epicenter of the 1966 Parkfield earthquake to zero at the southeastern end of the 1966 rupture zone. Data from a network of trilateration lines spanning the San Andreas fault near Parkfield and extending to the Pacific coast near San Luis Obispo shed light on the rate of fault slip at depth since the 1966 earthquake. In this study, average rates of line length change and shallow fault slip were inverted to determine the slip rate at depth on the Parkfield fault segment. The fault is taken to be a vertical surface with unknown distribution of strike-slip displacement in an elastic half-space. A striking result of the inversions is that all solutions providing acceptable fits to the data exhibit a locked zone essentially coincident with the rupture surface of the 1966 Parkfield earthquake. The data require that the locked zone extend nearly as far north as the 1966 epicenter; however, the vertical extent of the locked zone is not well resolved. Over much of the Parkfield segment the fault is slipping faster at the earth's surface than it is at seismogenic depths. In order to fit the trilateration measurements it is necessary to include a component of contraction normal to the trend of the San Andreas. The inversion results suggest a spatially uniform normal strain of −0.06 μstrain/yr. The orientation of the contraction is compatible with geologic and seismic evidence of active folding and reverse faulting in the region. The magnitude of the contraction is consistent with convergence rates inferred from global plate motion models.

Journal of Geophysical Research Solid Earth↗

USGS2011

Verifying a computational method for predicting extreme ground motion

In situations where seismological data is rare or nonexistent, computer simulations may be used to predict ground motions caused by future earthquakes. This is particularly practical in the case of extreme ground motions, where engineers of special buildings may need to design for an event that has not been historically observed but which may occur in the far-distant future. Once the simulations have been performed, however, they still need to be tested. The SCEC-USGS dynamic rupture code verification exercise provides a testing mechanism for simulations that involve spontaneous earthquake rupture. We have performed this examination for the specific computer code that was used to predict maximum possible ground motion near Yucca Mountain. Our SCEC-USGS group exercises have demonstrated that the specific computer code that was used for the Yucca Mountain simulations produces similar results to those produced by other computer codes when tackling the same science problem. We also found that the 3D ground motion simulations produced smaller ground motions than the 2D simulations.

Seismological Research Letters↗

USGS2009

The SCEC/USGS dynamic earthquake rupture code verification exercise

Numerical simulations of earthquake rupture dynamics are now common, yet it has been difficult to test the validity of these simulations because there have been few field observations and no analytic solutions with which to compare the results. This paper describes the Southern California Earthquake Center/U.S. Geological Survey (SCEC/USGS) Dynamic Earthquake Rupture Code Verification Exercise, where codes that simulate spontaneous rupture dynamics in three dimensions are evaluated and the results produced by these codes are compared using Web-based tools. This is the first time that a broad and rigorous examination of numerous spontaneous rupture codes has been performed—a significant advance in this science. The automated process developed to attain this achievement provides for a future where testing of codes is easily accomplished. Scientists who use computer simulations to understand earthquakes utilize a range of techniques. Most of these assume that earthquakes are caused by slip at depth on faults in the Earth, but hereafter the strategies vary. Among the methods used in earthquake mechanics studies are kinematic approaches and dynamic approaches. The kinematic approach uses a computer code that prescribes the spatial and temporal evolution of slip on the causative fault (or faults). These types of simulations are very helpful, especially since they can be used in seismic data inversions to relate the ground motions recorded in the field to slip on the fault(s) at depth. However, these kinematic solutions generally provide no insight into the physics driving the fault slip or information about why the involved fault(s) slipped that much (or that little). In other words, these kinematic solutions may lack information about the physical dynamics of earthquake rupture that will be most helpful in forecasting future events. To help address this issue, some researchers use computer codes to numerically simulate earthquakes and construct dynamic, spontaneous rupture (hereafter called “spontaneous rupture”) solutions. For these types of numerical simulations, rather than prescribing the slip function at each location on the fault(s), just the friction constitutive properties and initial stress conditions are prescribed. The subsequent stresses and fault slip spontaneously evolve over time as part of the elasto-dynamic solution. Therefore, spontaneous rupture computer simulations of earthquakes allow us to include everything that we know, or think that we know, about earthquake dynamics and to test these ideas against earthquake observations.

Seismological Research Letters↗

USGS2006

Introduction to the special issue on the 2004 Parkfield earthquake and the Parkfield earthquake prediction experiment

The 28 September 2004 M 6.0 Parkfield earthquake, a long-anticipated event on the San Andreas fault, is the world's best recorded earthquake to date, with state-of-the-art data obtained from geologic, geodetic, seismic, magnetic, and electrical field networks. This has allowed the preearthquake and postearthquake states of the San Andreas fault in this region to be analyzed in detail. Analyses of these data provide views into the San Andreas fault that show a complex geologic history, fault geometry, rheology, and response of the nearby region to the earthquake-induced ground movement. Although aspects of San Andreas fault zone behavior in the Parkfield region can be modeled simply over geological time frames, the Parkfield Earthquake Prediction Experiment and the 2004 Parkfield earthquake indicate that predicting the fine details of future earthquakes is still a challenge. Instead of a deterministic approach, forecasting future damaging behavior, such as that caused by strong ground motions, will likely continue to require probabilistic methods. However, the Parkfield Earthquake Prediction Experiment and the 2004 Parkfield earthquake have provided ample data to understand most of what did occur in 2004, culminating in significant scientific advances.

Bulletin of the Seismological Society of America↗

USGS2005

Material contrast does not predict earthquake rupture propagation direction

Earthquakes often occur on faults that juxtapose different rocks. The result is rupture behavior that differs from that of an earthquake occurring on a fault in a homogeneous material. Previous 2D numerical simulations have studied simple cases of earthquake rupture propagation where there is a material contrast across a fault and have come to two different conclusions: 1) earthquake rupture propagation direction can be predicted from the material contrast, and 2) earthquake rupture propagation direction cannot be predicted from the material contrast. In this paper we provide observational evidence from 70 years of earthquakes at Parkfield, CA, and new 3D numerical simulations. Both the observations and the numerical simulations demonstrate that earthquake rupture propagation direction is unlikely to be predictable on the basis of a material contrast. Copyright 2005 by the American Geophysical Union.

Geophysical Research Letters↗

USGS2005

The wrinkle-like slip pulse is not important in earthquake dynamics

A particular solution for slip on an interface between different elastic materials, the wrinkle-like slip pulse, propagates in only one direction with reduced normal compressive stress. More general solutions, and natural earthquakes, need not share those properties. In a 3D dynamic model with a drop in friction and heterogeneous initial stress, the wrinkle-like slip pulse is only a small part of the solution. Rupture propagation is determined primarily by the potential stress drop, not by the wrinkle-like slip pulse. A 2D calculation with much finer resolution shows that energy loss to friction might not be significantly reduced in the wrinkle-like slip pulse. Copyright 2005 by the American Geophysical Union.

Geophysical Research Letters↗

USGS2004

Numerical simulations of large earthquakes: Dynamic rupture propagation on heterogeneous faults

Our current conceptions of earthquake rupture dynamics, especially for large earthquakes, require knowledge of the geometry of the faults involved in the rupture, the material properties of the rocks surrounding the faults, the initial state of stress on the faults, and a constitutive formulation that determines when the faults can slip. In numerical simulations each of these factors appears to play a significant role in rupture propagation, at the kilometer length scale. Observational evidence of the earth indicates that at least the first three of the elements, geometry, material, and stress, can vary over many scale dimensions. Future research on earthquake rupture dynamics needs to consider at which length scales these features are significant in affecting rupture propagation. ?? Birkha??user Verlag, Basel, 2004.

Pure and Applied Geophysics↗

USGS2004

Inverse kinematic and forward dynamic models of the 2002 Denali fault earthquake, Alaska

We perform inverse kinematic and forward dynamic models of the M 7.9 2002 Denali fault, Alaska, earthquake to shed light on the rupture process and dynamics of this event, which took place on a geometrically complex fault system in central Alaska. We use a combination of local seismic and Global Positioning System (GPS) data for our kinematic inversion and find that the slip distribution of this event is characterized by three major asperities on the Denali fault. The rupture nucleated on the Susitna Glacier thrust fault, and after a pause, propagated onto the strike-slip Denali fault. Approximately 216 km to the east, the rupture abandoned the Denali fault in favor of the more southwesterly directed Totschunda fault. Three-dimensional dynamic models of this event indicate that the abandonment of the Denali fault for the Totschunda fault can be explained by the Totschunda fault's more favorable orientation with respect to the local stress field. However, a uniform tectonic stress field cannot explain the complex slip pattern in this event. We also find that our dynamic models predict discontinuous rupture from the Denali to Totschunda fault segments. Such discontinuous rupture helps to qualitatively improve our kinematic inverse models. Two principal implications of our study are (1) a combination of inverse and forward modeling can bring insight into earthquake processes that are not possible with either technique alone, and (2) the stress field on geometrically complex fault systems is most likely not due to a uniform tectonic stress field that is resolved onto fault segments of different orientations; rather, other forms of stress heterogeneity must be invoked to explain the observed slip patterns.

Alaska↗

USGS2003

The 2002 Denali fault earthquake, Alaska: A large magnitude, slip-partitioned event

The MW (moment magnitude) 7.9 Denali fault earthquake on 3 November 2002 was associated with 340 kilometers of surface rupture and was the largest strike-slip earthquake in North America in almost 150 years. It illuminates earthquake mechanics and hazards of large strike-slip faults. It began with thrusting on the previously unrecognized Susitna Glacier fault, continued with right-slip on the Denali fault, then took a right step and continued with right-slip on the Totschunda fault. There is good correlation between geologically observed and geophysically inferred moment release. The earthquake produced unusually strong distal effects in the rupture propagation direction, including triggered seismicity.

Science↗

USGS2002

The 1999 Izmit, Turkey, earthquake: A 3D dynamic stress transfer model of intraearthquake triggering

Before the August 1999 Izmit (Kocaeli), Turkey, earthquake, theoretical studies of earthquake ruptures and geological observations had provided estimates of how far an earthquake might jump to get to a neighboring fault. Both numerical simulations and geological observations suggested that 5 km might be the upper limit if there were no transfer faults. The Izmit earthquake appears to have followed these expectations. It did not jump across any step-over wider than 5 km and was instead stopped by a narrower step-over at its eastern end and possibly by a stress shadow caused by a historic large earthquake at its western end. Our 3D spontaneous rupture simulations of the 1999 Izmit earthquake provide two new insights: (1) the west- to east-striking fault segments of this part of the North Anatolian fault are oriented so as to be low-stress faults and (2) the easternmost segment involved in the August 1999 rupture may be dipping. An interesting feature of the Izmit earthquake is that a 5-km-long gap in surface rupture and an adjacent 25° restraining bend in the fault zone did not stop the earthquake. The latter observation is a warning that significant fault bends in strike-slip faults may not arrest future earthquakes.

Izmit↗

USGS2002

The 1999 Mw 7.1 Hector Mine, California, earthquake: A test of the stress shadow hypothesis?

We test the stress shadow hypothesis for large earthquake interactions by examining the relationship between two large earthquakes that occurred in the Mojave Desert of southern California, the 1992 M w 7.3 Landers and 1999 M w 7.1 Hector Mine earthquakes. We want to determine if the 1999 Hector Mine earthquake occurred at a location where the Coulomb stress was increased (earthquake advance, stress trigger) or decreased (earthquake delay, stress shadow) by the previous large earthquake. Using four models of the Landers rupture and a range of possible hypocentral planes for the Hector Mine earthquake, we discover that most scenarios yield a Landers-induced relaxation (stress shadow) on the Hector Mine hypocentral plane. Although this result would seem to weigh against the stress shadow hypothesis, the results become considerably more uncertain when the effects of a nearby Landers aftershock, the 1992 M L 5.4 Pisgah earthquake, are taken into account. We calculate the combined static Coulomb stress changes due to the Landers and Pisgah earthquakes to range from -0.3 to +0.3 MPa (- 3 to +3 bars) at the possible Hector Mine hypocenters, depending on choice of rupture model and hypocenter. These varied results imply that the Hector Mine earthquake does not provide a good test of the stress shadow hypothesis for large earthquake interactions. We use a simple approach, that of static dislocations in an elastic half-space, yet we still obtain a wide range of both negative and positive Coulomb stress changes. Our findings serve as a caution that more complex models purporting to explain the triggering or shadowing relationship between the 1992 Landers and 1999 Hector Mine earthquakes need to also consider the parametric and geometric uncertainties raised here.

California↗

USGS2000

Earthquake stress triggers, stress shadows, and seismic hazard

Many aspects of earthquake mechanics remain an enigma at the beginning of the twenty-first century. One potential bright spot is the realization that simple calculations of stress changes may explain some earthquake interactions, just as previous and ongoing studies of stress changes have begun to explain human- induced seismicity. This paper, which is an update of Harris1, reviews many published works and presents a compilation of quantitative earthquake-interaction studies from a stress change perspective. This synthesis supplies some clues about certain aspects of earthquake mechanics. It also demonstrates that much work remains to be done before we have a complete story of how earthquakes work.

Current Science↗

USGS1999

Dynamic 3D simulations of earthquakes on en echelon faults

One of the mysteries of earthquake mechanics is why earthquakes stop. This process determines the difference between small and devastating ruptures. One possibility is that fault geometry controls earthquake size. We test this hypothesis using a numerical algorithm that simulates spontaneous rupture propagation in a three-dimensional medium and apply our knowledge to two California fault zones. We find that the size difference between the 1934 and 1966 Parkfield, California, earthquakes may be the product of a stepover at the southern end of the 1934 earthquake and show how the 1992 Landers, California, earthquake followed physically reasonable expectations when it jumped across en echelon faults to become a large event. If there are no linking structures, such as transfer faults, then strike-slip earthquakes are unlikely to propagate through stepover s >5 km wide.

Geophysical Research Letters↗

USGS1998

Introduction to special section: Stress triggers, stress shadows, and implications for seismic hazard

Many aspects of earthquake mechanics remain an enigma as we enter the closing years of the twentieth century. One potential bright spot is the realization that simple calculations of stress changes may explain some earthquake interactions, just as previous and on going studies of stress changes have begun to explain human-induced seismicity. This paper, which introduces the special section “Stress Triggers, Stress Shadows, and Implications for Seismic Hazard,” reviews many published works and presents a compilation of quantitative earthquake interaction studies from a stress change perspective. This synthesis supplies some clues about certain aspects of earthquake mechanics. It also demonstrates that much work remains before we can understand the complete story of how earthquakes work.

Journal of Geophysical Research B: Solid Earth↗

USGS1996

In the shadow of 1857-the effect of the great Ft. Tejon earthquake on subsequent earthquakes in southern California

The great 1857 Fort Tejon earthquake is the largest earthquake to have hit southern California during the historic period. We investigated if seismicity patterns following 1857 could be due to static stress changes generated by the 1857 earthquake. When post-1857 earthquakes with unknown focal mechanisms were assigned strike-slip mechanisms with strike and rake determined by the nearest active fault, 13 of the 13 southern California M≥5.5 earthquakes between 1857 and 1907 were encouraged by the 1857 rupture. When post-1857 earthquakes in the Transverse Ranges with unknown focal mechanisms were assigned reverse mechanisms and all other events were assumed strike-slip, 11 of the 13 earthquakes were encouraged by the 1857 earthquake. These results show significant correlations between static stress changes and seismicity patterns. The correlation disappears around 1907, suggesting that tectonic loading began to overwhelm the effect of the 1857 earthquake early in the 20th century.

Geophysical Research Letters↗

USGS1995

Influence of static stress changes on earthquake locations in southern California

Earthquakes induce changes in static stress on neighbouring faults that may delay, hasten or even trigger subsequent earthquakes 1–10 . The length of time over which such effects persist has a bearing on the potential contribution of stress analyses to earthquake hazard assessment, but is presently unknown. Here we use an elastic half-space model 11 to estimate the static stress changes generated by damaging (magnitude M ≥5) earthquakes in southern California over the past 26 years, and to investigate the influence of these changes on subsequent earthquake activity. We find that, in the 1.5-year period following a M ≥5 earthquake, any subsequent nearby M ≥5 earthquake almost always ruptures a fault that is loaded towards failure by the first earthquake. After this period, damaging earthquakes are equally likely to rupture loaded and relaxed faults. Our results suggest that there is a short period of time following a damaging earthquake in southern California in which simple Coulomb failure stress models could be used to identify regions of increased seismic hazard.

California↗

USGS1992

Changes in static stress on southern California faults after the 1992 Landers earthquake

The magnitude 7.5 Landers earthquake of 28 June 1992 was the largest earthquake to strike California in 40 years. The slip that occurs in such an earthquake would be expected to induce large changes in the static stress on neighbouring faults; these changes in stress should in turn affect the likelihood of future earthquakes. Stress changes that load faults towards failure have been cited as the cause of small 1–5 , moderate 6 and large 7 earthquakes; conversely, those that relax neighbouring faults have been related to a decrease in seismicity 5 . Here we use an elastic half-space model 8 to estimate the stress changes produced by the Landers earthquake on selected southern California faults, including the San Andreas. We find that the estimated stress changes are consistent with the triggering of four out of the five aftershocks with magnitude greater than 4.5, and that the largest changes (1–10 bar), occurring on part of the San Bernardino segment of the San Andreas fault, may have decreased the time to the next magnitude 8 earthquake by about 14 years.

California↗