USGS Science⌕ Search

Geology topics

Leif Wennerberg

Publications and source records attributed to Leif Wennerberg.

11 recordsLinked to original sources
USGS2012

Influence of welded boundaries in anelastic media on energy flow, and characteristics of P, S-I, and S-II waves: Observational evidence for inhomogeneous body waves in low-loss solids

A general computer code, developed to calculate anelastic reflection-refraction coefficients, energy flow, and the physical characteristics for general P , S -I, and S -II waves, quantitatively describes physical characteristics for wave fields in anelastic media that do not exist in elastic media. Consideration of wave fields incident on boundaries between anelastic media shows that scattered wave fields experience reductions in phase and energy speeds, increases in maximum attenuation and Q −1 , and directions of maximum energy flow distinct from phase propagation. Each of these changes in physical characteristics are shown to vary with angle of incidence. Finite relaxation times for anelastic media result in energy flow due to interaction of superimposed radiation fields and contribute to energy flow across anelastic boundaries for all angles of incidence. Agreement of theoretical and numerical results with laboratory measurements argues for the validity of the theoretical and numerical formulations incorporating inhomogeneous wave fields. The agreement attests to the applicability of the model and helps confirm the existence of inhomogeneous body waves and their associated set of distinct physical characteristics in the earth. The existence of such body waves in layered, low-loss anelastic solids implies the need to reformulate some seismological models of the earth. The exact anelastic formulation for a liquid-solid interface with no low-loss approximations predicts the existence of a range of angles of incidence or an anelastic Rayleigh window, through which significant amounts of energy are transmitted across the boundary. The window accounts for the discrepancy apparent between measured reflection data presented in early textbooks and predictions based on classical elasticity theory. Characteristics of the anelastic Rayleigh window are expected to be evident in certain sets of wide-angle, ocean-bottom reflection data and to be useful in estimating Q −1 for some ocean bottom reflectors.

Journal of Geophysical Research B: Solid Earth↗
USGS1997

Bulk-friction modeling of afterslip and the modified Omori law

Afterslip data from the Superstition Hills fault in southern California, a creep event on the same fault, the modified Omori law, and cumulative moments from aftershocks of the 1957 Aleutian Islands earthquake all indicate that the original formulation by Dieterich (1981) [Constitutive properties of faults with simulated gouge. AGU, Geophys. Monogr. 24, 103–120] for friction evolution is more appropriate for systems far from instability than the commonly used approximation developed by Ruina (1983) [Slip instability and state variable friction laws. J. Geophys. Res. 88, 10359–10370] to study instability. The mathematical framework we use to test the friction models is a one-dimensional, massless spring-slider under the simplifying assumption, proposed by Scholz (1990) [The Mechanics of Earthquakes and Faulting. Cambridge University Press] and used by Marone et al. (1991) [On the mechanics of earthquake afterslip. J. Geophys. Res., 96: 8441–8452], that the state variable takes on its velocity-dependent steady-state value throughout motion in response to a step in stress. This assumption removes explicit state-variable dependence from the model, obviating the need to consider state-variable evolution equations. Anti-derivatives of the modified Omori law fit our data very well and are very good approximate solutions to our model equations. A plausible friction model with Omori-law solutions used by Wesson (1988) [Dynamics of fault creep. J. Geophys. Res. 93, 8929–8951] to model fault creep and generalized by Rice (1983) [Constitutive relations for fault slip and earthquake instabilities. Pure Appl. Geophys. 121, 443–475] to a rate-and-state variable friction model yields exactly Omori's law with exponents greater than 1, but yields unstable solutions for Omori exponents less than 1. We estimate from the Dieterich formulation the dimensionless parameter a∗ which is equal to the product of the nominal coefficient of friction and the more commonly reported friction parameter a . We find that a∗ is typically positive, qualitatively consistent with laboratory observations, although our observations are considerably larger than laboratory values. However, we also find good model fits for a∗ < 0 when data correspond to Omori exponents less than 1. A modification of the stability analysis by Rice and Ruina (1983) [Stability of steady frictional slipping. J. Appl. Mech. 50, 343–349] indicates that a∗ < 0 is not a consequence of our assumption regarding state-variable evolution. A consistent interpretation of a∗ < 0 in terms of laboratory models appears to be that the data are from later portions of processes better characterized by two-state-variable friction models. a∗ < 0 is explained by assuming that our data cannot resolve the co-seismic evolution of a short-length-scale state variable to a velocity-weakening state; our parameterization leads to an apparent negative instantaneous viscosity. We estimate the largest critical slip distance associated with afterslip to be ∼1–10 cm, consistent with other estimates for near-surface materials. We assume that our observed large values for a∗ reflect the fact that our model ignores the geometrical complexities of three-dimensional stresses in fractured crustal materials around a fault zone with frictional stresses that vary on a fault surface. Our one-dimensional model parameters reflect spatially averaged, bulk , stress and frictional properties of a fault zone, where we clearly cannot specify the details of the averaging process. Our analysis of Omori's law suggests that bulk-frictional properties of a fault zone are well described by our simple laboratory-based models, but they would need to change during the seismic cycle for a mainshock instability to recur, unless a mainshock-aftershock sequence were characterized by a process similar to the arrested instabilities possible in two-state-variable systems.

Tectonophysics↗
USGS1997

The Loma Prieta, California, earthquake of October 17, 1989: Aftershocks and postseismic effects

While the damaging effects of the earthquake represent a significant social setback and economic loss, the geophysical effects have produced a wealth of data that have provided important insights into the structure and mechanics of the San Andreas Fault system. Generally, the period after a large earthquake is vitally important to monitor. During this part of the seismic cycle, the primary fault and the surrounding faults, rock bodies, and crustal fluids rapidly readjust in response to the earthquake's sudden movement. Geophysical measurements made at this time can provide unique information about fundamental properties of the fault zone, including its state of stress and the geometry and frictional/rheological properties of the faults within it. Because postseismic readjustments are rapid compared with corresponding changes occurring in the preseismic period, the amount and rate of information that is available during the postseismic period is relatively high. From a geophysical viewpoint, the occurrence of the Loma Prieta earthquake in a section of the San Andreas fault zone that is surrounded by multiple and extensive geophysical monitoring networks has produced nothing less than a scientific bonanza. The reports assembled in this chapter collectively examine available geophysical observations made before and after the earthquake and model the earthquake's principal postseismic effects. The chapter covers four broad categories of postseismic effect: (1) aftershocks; (2) postseismic fault movements; (3) postseismic surface deformation; and (4) changes in electrical conductivity and crustal fluids.

California↗
USGS1996

Site amplification at five locations in San Francisco, California: A comparison of S waves, codas, and microtremors

We compare microtremor data to weak-motion S -wave and coda recordings at sites in San Francisco in order to clarify the range of applicability of microtremor data to ground-motion prediction. We also compare S -wave results to coda results. For each type of data, we compute spectral ratios of motions from two soil/rock station pairs and from an uphole/downhole pair in the Marina district. We compute horizontal/vertical ratios (Nakamura's method) at a soil site, a rock site, and the surface and borehole instruments. In the station-pair analyses, microtremor data show amplifications at the same fundamental frequency as S waves, but the frequencies of other peaks do not agree. The amplification at frequencies higher than 2 Hz is greater in the microtremor data. Station-pair ratios of coda data generally show spectral peaks occurring at the same frequencies, but with levels varying from one to four times the amplification from S -wave ratios. Nakamura's method of analyzing microtremors agrees better with S -wave station-pair results than the microtremor station-pair method over a limited frequency band that varies from station to station.

California↗
USGS1992

The Loma Prieta, California, earthquake of October, 17, 1989: Marina District

During the earthquake, a total land area of about 4,300 km2 was shaken with seismic intensities that can cause significant damage to structures. The area of the Marina District of San Francisco is only 4.0 km2--less than 0.1 percent of the area most strongly affected by the earthquake--but its significance with respect to engineering, seismology, and planning far outstrips its proportion of shaken terrain and makes it a centerpiece for lessons learned from the earthquake. The Marina District provides perhaps the most comprehensive case history of seismic effects at a specific site developed for any earthquake. The reports assembled in this chapter, which provide an account of these seismic effects, constitute a unique collection of studies on site, as well as infrastructure and societal, response that cover virtually all aspects of the earthquake, ranging from incoming ground waves to the outgoing airwaves used for emergency communication. The Marina District encompasses the area bounded by San Francisco Bay on the north, the Presidio on the west, and Lombard Street and Van Ness Avenue on the south and east, respectively. Nearly all of the earthquake damage in the Marina District, however, occurred within a considerably smaller area of about 0.75 km2, bounded by San Francisco Bay and Baker, Chestnut, and Buchanan Streets. At least five major aspects of earthquake response in the Marina District are covered by the reports in this chapter: (1) dynamic site response, (2) soil liquefaction, (3) lifeline performance, (4) building performance, and (5) emergency services.

California↗
USGS1990

Stochastic summation of empirical Green's functions

Two simple strategies are presented that use random delay times for repeatedly summing the record of a relatively small earthquake to simulate the effects of a larger earthquake. The simulations do not assume any fault plane geometry or rupture dynamics, but realy only on the ω−2 spectral model of an earthquake source and elementary notions of source complexity. The strategies simulate ground motions for all frequencies within the bandwidth of the record of the event used as a summand. The first strategy, which introduces the basic ideas, is a single-stage procedure that consists of simply adding many small events with random time delays. The probability distribution for delays has the property that its amplitude spectrum is determined by the ratio of ω−2 spectra, and its phase spectrum is identically zero. A simple expression is given for the computation of this zero-phase scaling distribution. The moment rate function resulting from the single-stage simulation is quite simple and hence is probably not realistic for high-frequency (>1 Hz) ground motion of events larger than ML ∼ 4.5 to 5. The second strategy is a two-stage summation that simulates source complexity with a few random subevent delays determined using the zero-phase scaling distribution, and then clusters energy around these delays to get an ω−2 spectrum for the sum. Thus, the two-stage strategy allows simulations of complex events of any size for which the ω−2 spectral model applies. Interestingly, a single-stage simulation with too few ω−2 records to get a good fit to an ω−2 large-event target spectrum yields a record whose spectral asymptotes are consistent with the ω−2 model, but that includes a region in its spectrum between the corner frequencies of the larger and smaller events reasonably approximated by a power law trend. This spectral feature has also been discussed as reflecting the process of partial stress release (Brune, 1970), an asperity failure (Boatwright, 1984), or the breakdown of ω−2 scaling due to rupture significantly longer than the width of the seismogenic zone (Joyner, 1984).

Bulletin of the Seismological Society of America↗
USGS1989

Rupture process of the Ms 6.6 Superstition Hills, California, earthquake determined from strong-motion recordings: application of tomographic source inversion

We analyze strong-motion recordings of the M s 6.6 Superstition Hills earthquake to determine the timing, location, spatial extent, and rupture velocity of the subevents that produced the bulk of the high-frequency (0.5 to 4 Hz) seismic energy radiated by this shock. The earthquake can be characterized by three principal subevents, the largest ones occurring about 3 and 10 sec after initiation of rupture. Timing relationships between pulses on the seismograms indicate that the three subevents are located within 8 km of each other along the northern portion of the Superstition Hills fault. The two largest subevents display different directivity effects. We apply a tomographic source inversion to the integrated accelerograms to determine the slip acceleration on the fault as a function of time and distance, based on a one-dimensional fault model. The azimuthal distribution of amplitudes for the second subevent can be largely explained by a rupture that propagated about 2 km to the southeast along the Superstition Hills fault at a velocity about equal to the P -wave velocity. An alternative model with rupture propagating to the northeast along a conjugate fault plane can also account for the observed directivity of this subevent, but it is not supported by the aftershock distribution. The third subevent ruptured to the southeast along an 8-km long portion of the Superstition Hills fault at about the shear-wave velocity. This rupture propagation caused the relatively large accelerations and velocities observed in strong-motion records for stations southeast of the hypocenter. The long time intervals between the subevents and their relative proximity to each other indicate a very slow component to the rupture development. The southern half of the Superstition Hills fault did not generate significant high-frequency strong ground motion, although it showed substantial co-seismic surface displacement. The subevents are situated along the same northern portion of the fault where most of the aftershocks are located. The locations of the subevents appear to be controlled by bends in the fault mapped at the surface and by changes in basement structure at depth.

California↗
USGS1989

Microearthquake spectra from the Anza, California, seismic network: Site response and source scaling

We analyzed spectra of local microearthquakes recorded by the Anza, California, seismic network to isolate the effects of site response and to investigate the scaling of source parameters for small earthquakes. Spectra of microearthquakes (M < 2; Mo < 1019 dyne-cm) at Anza have shapes characteristic of the receiver sites and are generally independent of the source region. Thus, the site response is a major conditioner of the observed spectral shape. To remove the effects of site response from the spectra of a M ∼ 3 event and isolate its source spectrum, we divided by the spectra of an adjacent aftershock used as an empirical Green's function event. The spectral ratios indicate that the apparent corner frequencies of small earthquakes (Mo < 1019dyne-cm) observed at even the high-fmax stations on hard rock are much lower than the source corner frequencies. The spectral ratios are consistent with stress drop remaining constant with decreasing seismic moment, for events with moments as small as 1018 dyne-cm. The spectral ratios display remarkable agreement between sites which showed vast differences in their original spectra, indicating that the spectral division effectively removed the site response. The source spectrum of the M ∼ 3 event has a high-frequency spectral fall-off of about ω−2. An apparent dependence of high-frequency fall-off with seismic moment in the original spectra can also be explained by the effects of site response. The difference between the P- and S-wave corner frequencies and high-frequency roll-offs in the observed spectra for these events is the result of the site response and is not a source property. The shapes of the spectra of microearthquakes at Anza can largely be explained by attenuation at shallow depth with a frequency-independent Q. For some sites, near-surface resonances are also apparent in the spectra of microearthquakes. It is indicated by t* values determined for each site that Qp ∼ Qs for the shallow low-Q layer. Further evidence of low near-surface Q is observed in the anomalous spectra of an unusually shallow earthquake (source depth ∼ 1 km) in the network. The spectra from this shallow event are depleted in high-frequency energy at most stations, relative to those of deeper events. This observation can be explained by a low-Q surficial zone. For stations of the network situated on alluvium, this low-Q layer has a maximum thickness of about 3 km and maximum P- and S-wave Q values of 30 to 50.

California↗
USGS1989

On the similarity of theories of anelastic and scattering attenuation

We point out basic parallels between theories of anelastic and scattering attenuation. We consider approximations to scattering effects presented by O'Doherty and Anstey (1971), Sato (1982), and Wu (1982). We use the linear theory of anelasticity. We note that the frequency dependence of Q can be related to a distribution of scales of physical properties of the medium. The frequency dependence of anelastic Q is related to the distribution of relaxation times in exactly the same manner as the frequency dependence of scattering Q is related to the distribution of scatterer sizes. Thus, the well-known difficulty of separating scattering from intrinsic attenuation is seen from this point of view as a consequence of the fact that certain observables can be interpreted by identical equations resulting from either of two credible physical theories describing fundamentally different processes.

Bulletin of the Seismological Society of America↗
USGS1985

General P, type-I S, and type-II S waves in anelastic solids; inhomogeneous wave fields in low-loss solids

The physical characteristics for general plane-wave radiation fields in an arbitrary linear viscoelastic solid are derived. Expressions for the characteristics of inhomogeneous wave fields, derived in terms of those for homogeneous fields, are utilized to specify the characteristics and a set of reference curves for general P and S wave fields in arbitrary viscoelastic solids as a function of wave inhomogeneity and intrinsic material absorption. The expressions show that an increase in inhomogeneity of the wave fields causes the velocity to decrease, the fractional-energy loss ( Q −1 ) to increase, the deviation of maximum energy flow with respect to phase propagation to increase, and the elliptical particle motions for P and type-I S waves to approach circularity. Q −1 for inhomogeneous type-I S waves is shown to be greater than that for type-II S waves, with the deviation first increasing then decreasing with inhomogeneity. The mean energy densities (kinetic, potential, and total), the mean rate of energy dissipation, the mean energy flux, and Q −1 for inhomogeneous waves are shown to be greater than corresponding characteristics for homogeneous waves, with the deviations increasing as the inhomogeneity is increased for waves of fixed maximum displacement amplitude. For inhomogeneous wave fields in low-loss solids, only the tilt of the particle motion ellipse for P and type-I S waves is independent to first order of the degree of inhomogeneity. Quantitative estimates for the characteristics of inhomogeneous plane body waves in layered low-loss solids are derived and guidelines established for estimating the effect of inhomogeneity on seismic body waves and a Rayleigh-type surface wave in low-loss media.

Bulletin of the Seismological Society of America↗