t* for S waves with a continental ray path

The purpose of this study was to determine t* for S waves with ray paths under the continental United States. The data set consists of long- and short-period body waves from the Borrego Mountain earthquake as observed in the northeastern U.S. The P wave forms are dominated by the sP phase and the SH wave forms by the sS. It is assumed that there are no losses in pure compression so that the relative attenuation rate of P and S waves is known. The initial source radiation is determined from the sP phase and the value of t_β* from the spectral content of the S wave. The results indicate that t_β* is 5.2 ± 0.7 sec along this ray path. Long- and short-period body waves from some deep South American events are used to test for lateral asymmetry of the Q distribution under the U.S. No lateral amplitude variation exists in this data, but this result is difficult to correlate with many previous results. The t_ β* value for a 600-km deep earthquake appears to be about 3 sec. A comparison of these values with values computed from current models of the Earth's Q distribution indicates that the models are slightly too high in Q overall and that more of the total body-wave attenuation occurs above 600 km than is indicated by the models

Modeling crustal structure through the use of converted phases in teleseismic body-wave forms

By comparing records of the radial component of motion of teleseismic P waves to records of the vertical component, it is possible to identify S phases within the P wave form. These phases are generated by the mechanism of P to S conversion at discontinuities in velocity under the receiving station. Similar phases of the S to P converted type appear as precursors to the direct SV arrival. Models for the crustal structure can be tested by generating synthetic seismograms for both components of motion of both the P and SV waves and comparing with the data. The technique has been used to model the crustal structure at WWSSN stations CAR and COR. It has also been used to check a recently proposed model for the crustal structure in eastern Canada which contains a large low-shear-velocity zone at the base of the crust. This study indicates that the crustal structure in eastern Canada is highly non-uniform with perhaps few features common to the whole region. Finally, the technique is used to identify several stations in the WWSSN which appear to be located on highly anomalous structure

A practical autonomous path planner for turn-of-the-century planetary microrovers

With the success of Mars Pathfinder's Sojourner rover, a new era of planetary exploration has opened, with demand for highly capable mobile robots. These robots must be able to traverse long distances over rough, unknown terrain autonomously, under severe resource constraints. Based on the authors' firsthand experience with the Mars Pathfinder mission, this paper reviews issues which are critical for successful autonomous navigation of planetary rovers. No currently proposed methodology addresses all of these issues. We next report on the 'Wedgebug' algorithm, which is applicable to planetary rover navigation in SE(2). The Wedgebug algorithm is complete, correct, requires minimal memory for storage of its worked model, and uses only on-board sensors, which are guided by the algorithm to efficiently senses only the data needed for motion planning. The implementation of a version of Wedgebug on the Rocky7 Mars Rover prototype at the Jet Propulsion Laboratory is described, and experimental results from operation in simulated martian terrain are presented

Inversion of the body waves from the Borrego Mountain earthquake to the source mechanism

The generalized linear inverse technique has been adapted to the problem of determining an earthquake source model from body-wave data. The technique has been successfully applied to the Borrego Mountain earthquake of April 9, 1968. Synthetic seismograms computed from the resulting model match in close detail the first 25 sec of long-period seismograms from a wide range of azimuths. The main shock source-time function has been determined by a new simultaneous short period-long period deconvolution technique as well as by the inversion technique. The duration and shape of this time function indicate that most of the body-wave energy was radiated from a surface with effective radius of only 8 km. This is much smaller than the total surface rupture length or the length of the aftershock zone. Along with the moment determination of M_o = 11.2 × 10^(25) dyne-cm, this radius implies a high stress drop of about 96 bars. Evidence in the amplitude data indicates that the polarization angle of shear waves is very sensitive to lateral structure

Time functions appropriate for nuclear explosions

The source-time function of megaton class nuclear explosions has been determined by modeling teleseismic short- and long-period body waves with synthetic seismograms. A simple analytic expression for the time function was used to closely match observations from both Novaya Zemlya and the U.S. test site at Amchitka. It was found that the time functions of all the events have a substantial overshoot. It was also found that, although the durations of the time functions did appear to depend on yield, the effect was very difficult to observe even in short-period records. All synthetics were computed by assuming a simple point source in a layered elastic half-space. It was not necessary to appeal to any nonlinear processes in the source region to explain the observations. Numerical calculations are presented to show that tectonic release triggered by earthquakes does not have a substantial effect on the P waves unless the long-period level of the tectonic event is as large or larger than the long-period level of the explosion. The pS wave, on the other hand, is shown to be very sensitive to even a moderate amount of tectonic release

Time functions appropriate for deep earthquakes

The seismic signatures of isolated body phases from many deep-focus earthquakes were analyzed in the time domain. Most shocks were found to be multiple events when examined in detail. The time history derived from P waves for single events predict synthetic S-wave shapes that match the observations, indicating compatibility with shear dislocation theory. Several other features of source functions in the time domain have been brought to light

Preliminary observations from the use of US-Soviet Joint Seismic Program data to model upper mantle triplications beneath Asia

New short-period waveform data from the US-Soviet Joint Seismic Program (JSP) make possible investigations of Asian upper mantle structure. the goal of this paper is to explore the potential use of the newly available JSP data to gain a qualitative view of upper mantle structure beneath Asia, and to facilitate more detailed future detailed future upper mantle studies. In a reconnaissance approach, waveform upper mantle studies. In a reconnaissance approach, waveform predictions from upper mantle P-wave velocity models of previous studies are compared to the JSP data to investigate regional differences in the central Asian upper mantle. Data coverage brackets the upper mantle triplications with excellent multi-source-to-stations sections. the abundance of data for controlled source-receiver geometries and the impulsive nature of the arrivals enable us to stack seismograms to improve signal-to-noise ratio. Arrivals from the 400 and 670 km discontinuities are apparent in the data and are compared to predictions of the mantle models. the principal result is that, for the regions studied, paths through cratonic regions of Asia are compatible with shield-type models, while paths through highly deformed regions of Asia are compatible with models derived for tectonically active regions, suggesting large lateral variations beneath the Eurasian continent. Use of the JSP data in a comparative approach is fast and simple, and proves effective in obtaining a first-order understanding of the Asian upper mantle. This result also presents the potential for qualitative studies elsewhere with digital portable stations

The source mechanism of the August 7, 1966 El Golfo earthquake

The El Golfo earthquake of August 7, 1966 (m_b = 6.3, M_S = 6.3) occurred near the mouth of the Colorado River at the northern end of the Gulf of California. Synthetic seismograms for this event were computed for both the body waves and the surface waves to determine the source parameters of the earthquake. The body-wave model indicated the source was a right lateral, strike-slip source with a depth of 10 km and a far-field time function 4 sec in duration. The body-wave moment was computed to be 5.0 × 10^(25) dyne-cm. The surface-wave radiation pattern was found to be consistent with that of the body waves with a surface-wave moment of 6.5 × 10^(25) dyne-cm. The agreement of the two different moments indicates that the earthquake had a simple source about 4 sec long. A comparison of this earthquake source with the Borrego Mountain and Truckee events demonstrates that all three of these earthquakes behaved as high stress-drop events. El Golfo was shown to be different from the low stress-drop, plate-boundary events which were located on the Gibbs fracture zone in 1967 and 1974

Time Domain Regional Discriminants

The time and frequency domains are equivalent displays of seismic trace, information, though some qualities of the signal are more easily observed in one domain than the other. The relative frequency excitation of Lg, for instance, is most easily viewed in the frequency domain, but such waveform qualities as the sequence in which pulses arrive in the wave train or the sharpness of pulse onset are most easily studied in the time domain (Murphy and Bennett, 1982, Blandford, 1981). Because of the tremendous complexity of high frequency regional data, most attempts at using it for discrimination purposes have involved analysis of the frequency content of the various arrivals either through transforming selected windows or through multiple bandpass filtering. We report here on our initial attempts to explore the alternative and to discriminate events using those waveform characteristics most easily observed in the time domain. A second advantage of time domain analysis approaches is that they permit a deeper insight into the physical processes creating a seismic signal's character. For this reason, they can be more e3silv used to evaluate the transportabilty of a discriminant to varying geophysical and tectonic regimes. This is an especially important feature in the development of regional discriminants. The most prominent and successful spectral regional discriminants have been empirically developed. This means that they must be redeveloped and reverified in each new area. As we shall show in the following, through rigorous time domain analysis such features as regional depth phases can be identified and used to discriminate. Discriminants based on such simple physical features as source depth should be transportable anywhere. In work recently completed under the treaty verification program, we have proved that such time domain discriminants do exist. In analyzing a test discrimination data set from the western U. S., we have discovered that the onset of P_n is always very similar for explosions and that few earthquakes have this unique waveform character. This information can be constructed into a simple discrimination scheme by testing the correlation of observed P_n waveform onsets with average waveforms observed from explosions. High correlations indicate explosions and low correlations earthquakes. We have also discovered that the regional phase P_g is actually composed of a sequence of sub-arrivals which correspond to successively higher orders of reverberation in the crust. In realistic crust models, the depth phases play an important role in the waveshapes of these sub-arrivals. By selecting an appropriate frequency band to analyze, we have been able to accurately model this type of data from explosions in the western United States. Over the very relevant regional distance ranges of 200 to 600 km, it appears that a discrimination procedure very similar to the one which is known to work for P_n will also be effective for P_g. We are investigating whether similar discriminants can be constructed based on the phases S_n and S_g in areas where those phases are prominent arrivals