Search USGS⌕ Search

Geology topics

R.W. Simpson

Publications and source records attributed to R.W. Simpson.

At least 19 recordsLinked to original sources
USGS2014

Northward migration of the Cascadia forearc in the northwestern U.S. and implications for subduction deformation

Geologic and paleomagnetic data from the Cascadia forearc indicate long-term northward migration and clockwise rotation of an Oregon coastal block with respect to North America. Paleomagnetic rotation of coastal Oregon is linked by a Klamath Mountains pole to geodetically and geologically determined motion of the Sierra Nevada block to derive a new Oregon Coast—North America (OC-NA) pole of rotation and velocity field. This long-term velocity field, which is independent of Pacific Northwest GPS data, is interpreted to be the result of Basin-Range extension and Pacific-North America dextral shear. The resulting Oregon Coast pole compares favorably to those derived solely from GPS data, although uncertainties are large. Subtracting the long-term motion from forearc GPS velocities reveals ENE motion with respect to an OC reference frame that is parallel to the direction of Juan de Fuca-OC convergence and decreases inland. We interpret this to be largely the result of subduction-related deformation. The adjusted mean GPS velocities are generally subparallel to those predicted from elastic dislocation models for Cascadia, but more definitive interpretations await refinement of the present large uncertainty in the Sierra Nevada block motion.

British Columbia, Washington↗
USGS2012

Introduction to special section on the California-Arizona crustal transect: CACTIS, Part 3

The CACTIS (California‐Arizona Crustal Transect Interim Synthesis) workshop in May 1988 brought researchers together in Flagstaff, Arizona, to discuss the geologic evolution and crustal structure of the southern Cordillera between the San Andreas fault in southeastern California and the Colorado Plateau in Arizona [ Sass et al ., 1988]. The first set of papers resulting from the workshop appeared in the Journal of Geophysical Research (JGR) as the special CACTIS 1 section in January 1990 [ Howard et al ., 1990], and a second as CACTIS 2 in November 1990 [ Haxel et al ., 1990]. Twelve papers in this issue form part 3 of the CACTIS series. A final collection of papers will appear in the future as a joint special section with the California Consortium for Crustal Studies (CALCRUST).

Journal of Geophysical Research B: Solid Earth↗
USGS2012

Introduction to special section on the California-Arizona crustal transect: Part II

In May 1988, a multidisciplinary conference was held in Flagstaff, Arizona, to synthesize recent research on the crustal structural and evolution of a transect across the southern Cordillera. This transect extends from the southern Colorado Plateau southwestward across the Arizona Transition Zone, through the Basin and Range province in the Mojave and Sonoran deserts, to the San Andreas fault system, and thence westward across the Imperial Valley and Peninsular Ranges and into the southern California Borderland. Most discussion concerned the eastern part of the transect, from the Colorado Plateau to the San Andreas fault [ Sass et al ., 1988; Howard et al ., 1990]. For this meeting, the name CACTIS (California‐Arizona Crustal Transect‐Interim Synthesis) was coined by J. Sass.

Journal of Geophysical Research↗
USGS2009

Geodetically inferred coseismic and postseismic slip due to the M 5.4 31 October 2007 Alum Rock earthquake

On 31 October 2007 the M 5.4 Alum Rock earthquake occurred near the junction between the Hayward and Calaveras faults in the San Francisco Bay Area, producing coseismic and postseismic displacements recorded by 10 continuously operating Global Positioning System (GPS) instruments. The cumulative postseismic displacements over the four months following the earthquake are linearly related to the cumulative number of aftershocks and are comparable in magnitude to the coseis mic displacements. The postseismic signal suggests that, in addition to afterslip at seismogenic depths, localized right-lateral/reverse slip occurred on dipping shallow fault surfaces southwest of the Calaveras. The spatial distribution of slip inferred by inverting the GPS data is compatible with a model in which moderate Calaveras fault earthquakes rupture locked patches surrounded by areas of creep, afterslip, and microseismicity (Oppenheimer et al., 1990). If this model and existing Calaveras fault slip rate estimates are correct, a slip deficit remains on the 2007 Alum Rock rupture patch that may be made up by aseismic slip or slip in larger earthquakes. Recent studies (e.g., Manaker et al., 2005) suggest that at depth the Hayward and central Calaveras faults connect via a simple continuous surface illuminated by the Mission Seismic Trend (MST), implying that a damaging earthquake rupture could involve both faults (Graymer et al., 2008). If this geometry is correct, the combined coseismic and postseismic slip we infer for the 2007 Alum Rock event predicts static Coulomb stress increases of ???0:6 bar on the MST surface and on the northern Calaveras fault ???5 km northwest of the Alum Rock hypocenter.

California↗
USGS2008

Three-Dimensional Geologic Map of the Hayward Fault Zone, San Francisco Bay Region, California

A three-dimensional (3D) geologic map of the Hayward Fault zone was created by integrating the results from geologic mapping, potential field geophysics, and seismology investigations. The map volume is 100 km long, 20 km wide, and extends to a depth of 12 km below sea level. The map volume is oriented northwest and is approximately bisected by the Hayward Fault. The complex geologic structure of the region makes it difficult to trace many geologic units into the subsurface. Therefore, the map units are generalized from 1:24,000-scale geologic maps. Descriptions of geologic units and structures are offered, along with a discussion of the methods used to map them and incorporate them into the 3D geologic map. The map spatial database and associated viewing software are provided. Elements of the map, such as individual fault surfaces, are also provided in a non-proprietary format so that the user can access the map via open-source software. The sheet accompanying this manuscript shows views taken from the 3D geologic map for the user to access. The 3D geologic map is designed as a multi-purpose resource for further geologic investigations and process modeling.

Scientific Investigations Map↗
USGS2007

Why the sacramento delta area differs from other parts of the great valley: numerical modeling of thermal structure and thermal subsidence of forearc basins

Data on present-day heat flow, subsidence history, and paleotemperature for the Sacramento Delta region, California, have been employed to constrain a numerical model of tectonic subsidence and thermal evolution of forearc basins. The model assumes an oceanic basement with an initial thermal profile dependent on its age subjected to refrigeration caused by a subducting slab. Subsidence in the Sacramento Delta region appears to be close to that expected for a forearc basin underlain by normal oceanic lithosphere of age 150 Ma, demonstrating that effects from both the initial thermal profile and the subduction process are necessary and sufficient. Subsidence at the eastern and northern borders of the Sacramento Valley is considerably less, approximating subsidence expected from the dynamics of the subduction zone alone. These results, together with other geophysical data, show that Sacramento Delta lithosphere, being thinner and having undergone deeper subsidence, must differ from lithosphere of the transitional type under other parts of the Sacramento Valley. Thermal modeling allows evaluation of the rheological properties of the lithosphere. Strength diagrams based on our thermal model show that, even under relatively slow deformation (10 −17 s −1 ), the upper part of the delta crystalline crust (down to 20–22 km) can fail in brittle fashion, which is in agreement with deeper earthquake occurrence. Hypocentral depths of earthquakes under the Sacramento Delta region extend to nearly 20 km, whereas, in the Coast Ranges to the west, depths are typically less than 12–15 km. The greater width of the seismogenic zone in this area raises the possibility that, for fault segments of comparable length, earthquakes of somewhat greater magnitude might occur than in the Coast Ranges to the west.

California↗
USGS2005

Seismic reflection evidence for a northeast-dipping Hayward fault near Fremont, California: Implications for seismic hazard

A 1.6-km-long seismic reflection profile across the creeping trace of the southern Hayward fault near Fremont, California, images the fault to a depth of 650 m. Reflector truncations define a fault dip of about 70 degrees east in the 100 to 650 m depth range that projects upward to the creeping surface trace, and is inconsistent with a nearly vertical fault in this vicinity as previously believed. This fault projects to the Mission seismicity trend located at 4-10 km depth about 2 km east of the surface trace and suggests that the southern end of the fault is as seismically active as the part north of San Leandro. The seismic hazard implication is that the Hayward fault may have a more direct connection at depth with the Calaveras fault, affecting estimates of potential event magnitudes that could occur on the combined fault surfaces, thus affecting hazard assessments for the south San Francisco Bay region.

Geophysical Research Letters↗
USGS2002

Spatial and temporal deformation along the northern San Jacinto fault, southern California: Implications for slip rates

The San Timoteo badlands is an area of uplift and erosional dissection that has formed as a result of late Quaternary uplift along a restraining bend in the San Jacinto fault, of the San Andreas fault system in southern California. This bend currently is located in a region where late Quaternary deposits and associated surfaces have formed in lower San Timoteo Canyon. We have used morphometric analysis of these surfaces, in conjunction with computer modeling of deformational patterns along the San Jacinto fault, to reconstruct spatial and temporal variations in uplift along the bend. Morphometric techniques used include envelope/subenvelope mapping, a gradient-length index along channels, and denudation values. Age control is determined using a combination of thermoluminescence (TL) and near infrared optical simulation luminescence dating (IROSL) and correlation of soil-development indices. These approaches are combined with an elastic half-space model used to determine the deformation associated with the fault bend. The region of modeled uplift has a similar distribution as that determined by morphometric techniques. Luminescence dates and soil-correlation age estimates generally agree. Based on soil development, surfaces within the study area were stabilized at approximately 300-700 ka for Q3, 43-67 ka for Q2, and 27.5-67 ka for Q1. Luminescence ages (both TL and IROSL) for the formation of the younger two surfaces are 58 to 94 ka for Q2 and 37 to 62 ka for Q1 (ages reported to 1?? uncertainty). Periods of uplift were determined for the surfaces in the study area, resulting in approximate uplift rates of 0.34 to 0.84 m/ka for the past 100 ka and 0.13 to 1.00 m/ka for the past 66 ka. Comparison of these rates of uplift to those generated by the model support a higher rate of lateral slip along the San Jacinto fault than commonly assumed (greater than 20 mm/yr, as compared to 8-12 mm/yr commonly cited). This higher slip rate supports the proposal that a greater amount of slip has transferred from the San Andreas fault to the San Jacinto fault than generally held. The San Jacinto fault may have accommodated a significant portion of the plate boundary slip during the past 100 ka.

California↗
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↗
USGS2002

Comments on potential geologic and seismic hazards affecting Mare Island, Solano County, California

This report was prepared in response to a written request from the City of Vallejo, California, to the U.S. Geological Survey (USGS). By letter of October 4, 2002, the City requested that the USGS "provide advice to the City’s LNG Health and Safety Committee on its review of a potential liquid natural gas project" on the southern portion of Mare Island. The City specifically requested that the USGS advise the committee on potential hazards including fault rupture, earthquake ground motion, soil failure during earthquakes, tsunami and seiche, and landslides. The City requested that the USGS: (1) comment on these hazards, (2) describe its degree of confidence in its opinions, and (3) describe the scope of additional studies that will be needed if the City enters into an agreement with project sponsors. Advice was also requested on the selection of the safe shutdown and operating basis earthquakes as specified in the NFPA 59A standard (NFPA, 2001). This review of published reports and other publicly available information indicates that all of the hazards on which the USGS was asked to comment should be considered for the proposed project on the southern portion of Mare Island. Available information differs greatly for each of these potential hazards, and adequate understanding for design will require detailed site-specific investigations.

California↗
USGS2001

Variations in creep rate along the Hayward Fault, California, interpreted as changes in depth of creep

Variations ill surface creep rate along the Hayward fault are modeled as changes in locking depth using 3D boundary elements. Model creep is driven by screw dislocations at 12 km depth under the Hayward and other regional faults. Inferred depth to locking varies along strike from 4-12 km. (12 km implies no locking.) Our models require locked patches under the central Hayward fault, consistent with a M6.8 earthquake in 1868, but the geometry and extent of locking under the north and south ends depend critically on assumptions regarding continuity and creep behavior of the fault at its ends. For the northern onshore part of the fault, our models contain 1.4-1.7 times more stored moment than the model of Bu??rgmann et al. [2000]; 45-57% of this stored moment resides in creeping areas. It is important for seismic hazard estimation to know how much of this moment is released coseismically or as aseismic afterslip.

Geophysical Research Letters↗
USGS2000

Pore fluid pressure, apparent friction, and Coulomb failure

Many recent studies of stress-triggered seismicity rely on a fault failure model with a single free parameter, the apparent coefficient of friction, presumed to be a material constant with possible values 0 ≤ μ′ ≤ 1. These studies may present a misleading view of fault strength and the role of pore fluid pressure in earthquake failure. The parameter μ′ is intended to incorporate the effects of both friction and pore pressure, but is a material constant only if changes in pore fluid pressure induced by changes in stress are proportional to the normal stress change across the potential failure plane. Although specific models of fault zones permit such a relation, neither is it known that fault zones within the Earth behave this way, nor is this behavior expected in all cases. In contrast, for an isotropic homogeneous poroelastic model the pore pressure changes are proportional to changes in mean stress, μ′ is not a material constant, and −∞ ≤ μ′ ≤ +∞. Analysis of the change in Coulomb failure stress for tectonically loaded reverse and strike-slip faults shows considerable differences between these two pore pressure models, suggesting that such models might be distinguished from one another using observations of triggered seismicity (e.g., aftershocks). We conclude that using the constant apparent friction model exclusively in studies of Coulomb failure stress is unwise and could lead to significant errors in estimated stress change and seismic hazard.

Journal of Geophysical Research B: Solid Earth↗
USGS1999

Stress sensitivity of fault seismicity: A comparison between limited-offset oblique and major strike-slip faults

We present a new three-dimensional inventory of the southern San Francisco Bay area faults and use it to calculate stress applied principally by the 1989 M = 7.1 Loma Prieta earthquake and to compare fault seismicity rates before and after 1989. The major high-angle right-lateral faults exhibit a different response to the stress change than do minor oblique (right-lateral/thrust) faults. Seismicity on oblique-slip faults in the southern Santa Clara Valley thrust belt increased where the faults were undamped. The strong dependence of seismicity change on normal stress change implies a high coefficient of static friction. In contrast, we observe that faults with significant offset (>50–100 km) behave differently; microseismicity on the Hayward fault diminished where right-lateral shear stress was reduced and where it was undamped by the Loma Prieta earthquake. We observe a similar response on the San Andreas fault zone in southern California after the Landers earthquake sequence. Additionally, the offshore San Gregorio fault shows a seismicity rate increase where right-lateral/oblique shear stress was increased by the Loma Prieta earthquake despite also being clamped by it. These responses are consistent with either a low coefficient of static friction or high pore fluid pressures within the fault zones. We can explain the different behavior of the two styles of faults if those with large cumulative offset become impermeable through gouge buildup; coseismically pressurized pore fluids could be trapped and negate imposed normal stress changes, whereas in more limited offset faults, fluids could rapidly escape. The difference in behavior between minor and major faults may explain why frictional failure criteria that apply intermediate coefficients of static friction can be effective in describing the broad distributions of aftershocks that follow large earthquakes, since many of these events occur both inside and outside major fault zones.

Journal of Geophysical Research B: Solid Earth↗
USGS1998

Strain accumulation rates in the San Francisco Bay area, 1972–1989

Maps of the strain accumulation rate in the San Francisco Bay area have been constructed from trilateration observations extending from about 1972 until the Loma Prieta earthquake in late 1989. The observations were corrected to remove offsets imposed by shallow fault creep and by four M ∼6 earthquakes that occurred in the Bay area during that time interval. The Bay area was divided into 32 contiguous polygons, and the uniform (in both space and time) strain rates that best explain the changes in the corrected (earthquake and shallow fault creep offsets removed) distances within each polygon were calculated. In a coordinate system with the 1 axis directed N58°E and the 2 axis N32°W (perpendicular and parallel to the local tangent to the small circle drawn about the Pacific-Sierra Nevada pole of rotation) the averages of these 32 strain rates (each weighted by the area of the polygon) are ε11=9.2±7.4, ε12=−160.7±4.6, and ε22=8.2±6.2 nanostrain/yr, where extension is reckoned positive and quoted uncertainties are standard deviations. As expected from the Pacific-Sierra Nevada relative plate motion, the overall strain rate is predominantly right-lateral shear across a vertical plane striking N32°W. The net increase in the 12,225 km2 area of the trilateration network is only 212±110 m2/yr, which arises from almost equal extensions in the N32°W and N58°E directions. Within the network the strain rates vary from polygon to polygon. Significant areal dilatation rates are observed in almost ⅓ of the individual polygons and the N32°W extension rates tend to be negative to the west of the Hayward-Rodgers Creek fault trend and positive east of it. The N58E extension rate is positive in 22 out of 32 polygons, a proportion that is significantly larger than would be expected by chance if the N58°E extension rate were zero or negative. The pre-1989 strain accumulation across the eventual site of the Loma Prieta rupture involves fault normal contraction as well as right-lateral shear, consistent with the rupture mechanism.

California↗
USGS1997

Stressing of the New Madrid seismic zone by a lower crust detachment fault

A new mechanical model for the cause of the New Madrid seismic zone in the central United States is analyzed. The model contains a subhorizontal detachment fault which is assumed to be near the domed top surface of locally thickened anomalous lower crust (“rift pillow”). Regional horizontal compression induces slip on the fault, and the slip creates a stress concentration in the upper crust above the rift pillow dome. In the coseismic stage of the model earthquake cycle, where the three largest magnitude 7–8 earthquakes in 1811–1812 are represented by a single model mainshock on a vertical northeast trending fault, the model mainshock has a moment equivalent to a magnitude 8 event. During the interseismic stage, corresponding to the present time, slip on the detachment fault exerts a right-lateral shear stress on the locked vertical fault whose failure produces the model mainshock. The sense of shear is generally consistent with the overall sense of slip of 1811–1812 and later earthquakes. Predicted rates of horizontal strain at the ground surface are about 10 −7 year −1 and are comparable to some observed rates. The model implies that rift pillow geometry is a significant influence on the maximum possible earthquake magnitude.

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↗