Geology topics

Fred F. Pollitz

Publications and source records attributed to Fred F. Pollitz.

11 recordsLinked to original sources

USGS2015

Interpretation of S waves generated by near-surface chemical explosions at SAFOD

A series of near-surface chemical explosions conducted at the San Andreas Fault Observatory at Depth (SAFOD) were recorded by high-frequency downhole receiver arrays in separate experiments in November 2003 and May 2005. The 2003 experiment involved &sim;100  kg shots detonated along a 46-km-long line (Hole–Ryberg line) centered on SAFOD and recorded by 32 three-component geophones in the pilot hole between 0.8 and 2.0 km depth. The 2005 experiment involved &sim;36  kg shots detonated at Parkfield Area Seismic Observatory (PASO) stations (at &sim;1–8  km offset) recorded by 80 three-component geophones in the main hole between the surface and 2.4 km depth. These data sample the downgoing seismic wavefield and constrain the shallow velocity and attenuation structure, as well as the first-order characteristics of the source. Using forward modeling on a velocity structure designed for the near field, both observed P - and S -wave energy for the PASO shots are identified with the travel times expected for direct and/or reflected phases. Larger-offset recordings from shots along the Hole–Ryberg line reveal substantial SV and SH energy, especially southwest of SAFOD from the source as indicated by P -to- S amplitude ratios. The generated SV energy is interpreted to arise chiefly from P -to- S conversions at subhorizontal discontinuities. This provides a simple mechanism for often-observed low P -to- S amplitude ratios from nuclear explosions in the far field, as originating from strong near-field wave conversions.

Bulletin of the Seismological Society of America↗

USGS2015

The Mw6.0 24 August 2014 South Napa earthquake

The M w 6.0 South Napa earthquake, which occurred at 10:20 UTC 24 August 2014 was the largest earthquake to strike the greater San Francisco Bay area since the M w 6.9 1989 Loma Prieta earthquake. The rupture from this right‐lateral earthquake propagated mostly unilaterally to the north and up‐dip, directing the strongest shaking toward the city of Napa, where peak ground accelerations (PGAs) between 45% g and 61% g were recorded and modified Mercalli intensities (MMIs) of VII–VIII were reported. Tectonic surface rupture with dextral slip of up to 46 cm was observed on a 12.5 km long segment, some of which was along a previously mapped strand of the West Napa fault system, although the rupture extended to the north of the mapped Quaternary strand. Modeling of seismic and geodetic data suggests an average coseismic slip of 50 cm, with a maximum slip of about 1 m at depths of 10–11 km. We observed up to 35 cm of afterslip along the surface trace in the week following the mainshock, primarily along the southern half of the surface rupture that experienced relatively little coseismic offset. Relocation of the sparse aftershock sequence suggests en echelon southwest‐ and northeast‐dipping fault planes, reflective of the complex fault geometry in this region. The Napa basin and historic and late Holocene alluvial flood deposits in downtown Napa amplified the ground motions there. Few ground failures were mapped, reflecting the dry season (as well as a persistent drought that had lowered the groundwater table) and the short duration of strong shaking in the epicentral area.

California↗

USGS2014

Key recovery factors for the August 24, 2014, South Napa Earthquake

Through discussions between the Federal Emergency Management Agency (FEMA) and the U.S. Geological Survey (USGS) following the South Napa earthquake, it was determined that several key decision points would be faced by FEMA for which additional information should be sought and provided by USGS and its partners. This report addresses the four tasks that were agreed to. These tasks are (1) assessment of ongoing fault movement (called afterslip) especially in the Browns Valley residential neighborhood, (2) assessment of the shaking pattern in the downtown area of the City of Napa, (3) improvement of information on the fault hazards posed by the West Napa Fault System (record of past earthquakes and slip rate, for example), and (4) imagery acquisition and data processing to provide overall geospatial information support to FEMA.

California↗

USGS2012

A comparison among observations and earthquake simulator results for the allcal2 California fault model

In order to understand earthquake hazards we would ideally have a statistical description of earthquakes for tens of thousands of years. Unfortunately the ∼100‐year instrumental, several 100‐year historical, and few 1000‐year paleoseismological records are woefully inadequate to provide a statistically significant record. Physics‐based earthquake simulators can generate arbitrarily long histories of earthquakes; thus they can provide a statistically meaningful history of simulated earthquakes. The question is, how realistic are these simulated histories? This purpose of this paper is to begin to answer that question. We compare the results between different simulators and with information that is known from the limited instrumental, historic, and paleoseismological data. As expected, the results from all the simulators show that the observational record is too short to properly represent the system behavior; therefore, although tests of the simulators against the limited observations are necessary, they are not a sufficient test of the simulators’ realism. The simulators appear to pass this necessary test. In addition, the physics‐based simulators show similar behavior even though there are large differences in the methodology. This suggests that they represent realistic behavior. Different assumptions concerning the constitutive properties of the faults do result in enhanced capabilities of some simulators. However, it appears that the similar behavior of the different simulators may result from the fault‐system geometry, slip rates, and assumed strength drops, along with the shared physics of stress transfer. This paper describes the results of running four earthquake simulators that are described elsewhere in this issue of Seismological Research Letters. The simulators ALLCAL (Ward, 2012), VIRTCAL (Sachs et al., 2012), RSQSim (Richards‐Dinger and Dieterich, 2012), and ViscoSim (Pollitz, 2012) were run on our most recent all‐California fault model, allcal2. With the exception of ViscoSim, which ran for 10,000 years, all the simulators ran for 30,000 years. Presentations containing content similar to this paper can be found at http://scec.usc.edu/research/eqsims/.

Seismological Research Letters↗

USGS2012

The 11 April 2012 east Indian Ocean earthquake triggered large aftershocks worldwide

Large earthquakes trigger very small earthquakes globally during passage of the seismic waves and during the following several hours to days 1, 2, 3, 4, 5, 6, 7, 8, 9, 10 , but so far remote aftershocks of moment magnitude M ≥5.5 have not been identified 11 , with the lone exception of an M =6.9 quake remotely triggered by the surface waves from an M =6.6 quake 4,800 kilometres away 12 . The 2012 east Indian Ocean earthquake that had a moment magnitude of 8.6 is the largest strike-slip event ever recorded. Here we show that the rate of occurrence of remote M ≥5.5 earthquakes (>1,500 kilometres from the epicentre) increased nearly fivefold for six days after the 2012 event, and extended in magnitude to M ≥7. These global aftershocks were located along the four lobes of Love-wave radiation; all struck where the dynamic shear strain is calculated to exceed 10 -7 for at least 100 seconds during dynamic-wave passage. The other M ≥8.5 mainshocks during the past decade are thrusts; after these events, the global rate of occurrence of remote M ≥5.5 events increased by about one-third the rate following the 2012 shock and lasted for only two days, a weaker but possibly real increase. We suggest that the unprecedented delayed triggering power of the 2012 earthquake may have arisen because of its strike-slip source geometry or because the event struck at a time of an unusually low global earthquake rate, perhaps increasing the number of nucleation sites that were very close to failure.

Indian Ocean;Sumatra↗

USGS2012

Illumination of rheological mantle heterogeneity by the M7.2 2010 El Mayor-Cucapah earthquake

Major intracontinental strike-slip faults tend to mark boundaries between lithospheric blocks of contrasting mechanical properties along much of their length. Both crustal and mantle heterogeneities can form such boundaries, but the role of crustal versus mantle strength contrasts for localizing strain sufficiently to generate major faults remains unclear. Using the crustal velocity field observed through the Global Positioning System (GPS) in the epicentral area of the M7.2 2010 El Mayor-Cucapah earthquake, Baja California, we find that transient deformation observed after the event is anomalously small in areas of relatively high seismic velocity in the shallow upper mantle (∼50 km depth). This pattern is best explained with a laterally heterogeneous viscoelastic structure that mimics the seismic structure. The mantle of the Southern Colorado River Desert (SCRD) and Peninsular Ranges (PR), which bound the fault system to its east and west, respectively, have anomalously high viscosity and seismic velocity. We hypothesize that compared with the rest of the San Andreas fault (SAF) system to its north, the strike-slip fault system in northern Baja California is narrow because of the presence of the PR and SCRD high-viscosity regions which bound it.

Arizona, Baja California, California, Sonora↗

USGS2011

High-frequency Born synthetic seismograms based on coupled normal modes

High-frequency and full waveform synthetic seismograms on a 3-D laterally heterogeneous earth model are simulated using the theory of coupled normal modes. The set of coupled integral equations that describe the 3-D response are simplified into a set of uncoupled integral equations by using the Born approximation to calculate scattered wavefields and the pure-path approximation to modulate the phase of incident and scattered wavefields. This depends upon a decomposition of the aspherical structure into smooth and rough components. The uncoupled integral equations are discretized and solved in the frequency domain, and time domain results are obtained by inverse Fourier transform. Examples show the utility of the normal mode approach to synthesize the seismic wavefields resulting from interaction with a combination of rough and smooth structural heterogeneities. This approach is applied to an ∼4 Hz shallow crustal wave propagation around the site of the San Andreas Fault Observatory at Depth (SAFOD).

Geophysical Journal International↗

USGS2011

Epistemic uncertainty in California-wide synthetic seismicity simulations

The generation of seismicity catalogs on synthetic fault networks holds the promise of providing key inputs into probabilistic seismic-hazard analysis, for example, the coefficient of variation, mean recurrence time as a function of magnitude, the probability of fault-to-fault ruptures, and conditional probabilities for foreshock–mainshock triggering. I employ a seismicity simulator that includes the following ingredients: static stress transfer, viscoelastic relaxation of the lower crust and mantle, and vertical stratification of elastic and viscoelastic material properties. A cascade mechanism combined with a simple Coulomb failure criterion is used to determine the initiation, propagation, and termination of synthetic ruptures. It is employed on a 3D fault network provided by Steve Ward (unpublished data, 2009) for the Southern California Earthquake Center (SCEC) Earthquake Simulators Group. This all-California fault network, initially consisting of 8000 patches, each of ∼12 square kilometers in size, has been rediscretized into Graphic patches, each of ∼1 square kilometer in size, in order to simulate the evolution of California seismicity and crustal stress at magnitude M∼5–8. Resulting synthetic seismicity catalogs spanning 30,000 yr and about one-half million events are evaluated with magnitude-frequency and magnitude-area statistics. For a priori choices of fault-slip rates and mean stress drops, I explore the sensitivity of various constructs on input parameters, particularly mantle viscosity. Slip maps obtained for the southern San Andreas fault show that the ability of segment boundaries to inhibit slip across the boundaries (e.g., to prevent multisegment ruptures) is systematically affected by mantle viscosity.

California↗

USGS2011

Coseismic slip distribution of the February 27, 2010 Mw 8.9 Maule, Chile earthquake

[1] Static offsets produced by the February 27, 2010 Mw = 8.8 Maule, Chile earthquake as measured by GPS and InSAR constrain coseismic slip along a section of the Andean megathrust of dimensions 650 km (in length) × 180 km (in width). GPS data have been collected from both campaign and continuous sites sampling both the near-field and far field. ALOS/PALSAR data from several ascending and descending tracks constrain the near-field crustal deformation. Inversions of the geodetic data for distributed slip on the megathrust reveal a pronounced slip maximum of order 15 m at ∼15–25 km depth on the megathrust offshore Lloca, indicating that seismic slip was greatest north of the epicenter of the bilaterally propagating rupture. A secondary slip maximum appears at depth ∼25 km on the megathrust just west of Concepción. Coseismic slip is negligible below 35 km depth. Estimates of the seismic moment based on different datasets and modeling approaches vary from 1.8 to 2.6 × 1022 N m. Our study is the first to model the static displacement field using a layered spherical Earth model, allowing us to incorporate both near-field and far-field static displacements in a consistent manner. The obtained seismic moment of 1.97 × 1022 N m, corresponding to a moment magnitude of 8.8, is similar to that obtained by previous seismic and geodetic inversions.

Geophysical Research Letters↗

USGS2010

On the resolution of shallow mantle viscosity structure using post-earthquake relaxation data: Application to the 1999 Hector Mine, California, earthquake

Most models of lower crust/mantle viscosity inferred from postearthquake relaxation assume one or two uniform-viscosity layers. A few existing models possess apparently significant radially variable viscosity structure in the shallow mantle (e.g., the upper 200 km), but the resolution of such variations is not clear. We use a geophysical inverse procedure to address the resolving power of inferred shallow mantle viscosity structure using postearthquake relaxation data. We apply this methodology to 9 years of GPS-constrained crustal motions after the 16 October 1999 M = 7.1 Hector Mine earthquake. After application of a differencing method to isolate the postearthquake signal from the “background” crustal velocity field, we find that surface velocities diminish from ∼20 mm/yr in the first few months to ≲2 mm/yr after 2 years. Viscoelastic relaxation of the mantle, with a time-dependent effective viscosity prescribed by a Burgers body, provides a good explanation for the postseismic crustal deformation, capturing both the spatial and temporal pattern. In the context of the Burgers body model (which involves a transient viscosity and steady state viscosity), a resolution analysis based on the singular value decomposition reveals that at most, two constraints on depth-dependent steady state mantle viscosity are provided by the present data set. Uppermost mantle viscosity (depth ≲ 60 km) is moderately resolved, but deeper viscosity structure is poorly resolved. The simplest model that explains the data better than that of uniform steady state mantle viscosity involves a linear gradient in logarithmic viscosity with depth, with a small increase from the Moho to 220 km depth. However, the viscosity increase is not statistically significant. This suggests that the depth-dependent steady state viscosity is not resolvably different from uniformity in the uppermost mantle.

California↗

USGS2009

A Viscoelastic earthquake simulator with application to the San Francisco Bay region

Earthquake simulation on synthetic fault networks carries great potential for characterizing the statistical patterns of earthquake occurrence. I present an earthquake simulator based on elastic dislocation theory. It accounts for the effects of interseismic tectonic loading, static stress steps at the time of earthquakes, and postearthquake stress readjustment through viscoelastic relaxation of the lower crust and mantle. Earthquake rupture initiation and termination are determined with a Coulomb failure stress criterion and the static cascade model. The simulator is applied to interacting multifault systems: one, a synthetic two-fault network, and the other, a fault network representative of the San Francisco Bay region. The faults are discretized both along strike and along dip and can accommodate both strike slip and dip slip. Stress and seismicity functions are evaluated over 30,000 yr trial time periods, resulting in a detailed statistical characterization of the fault systems. Seismicity functions such as the coefficient of variation and a- and b-values exhibit systematic patterns with respect to simple model parameters. This suggests that reliable estimation of the controlling parameters of an earthquake simulator is a prerequisite to the interpretation of its output in terms of seismic hazard.

California↗