Publications

Abstract not provided.

Muons are subatomic particles capable of penetrating the earth's crust several kilometers. Muons have been used to image the Pyramid of Khafre of Giza, various volcanoes, and smaller targets like cargo. For objects like a volcano, the detector is placed at the volcano's base and muon fluxes for paths through the volcano are recorded for many days to weeks.

Abstract not provided.

Abstract not provided.

Seismic attenuation is defined as the loss of the seismic wave amplitude as the wave propagates excluding losses strictly due to geometric spreading. Information gleaned from seismic waves can be utilized to solve for the attenuation properties of the earth. One method of solving for earth attenuation properties is called t*. This report will start by introducing the basic theory behind t* and delve into inverse theory as it pertains to how the algorithm called tstarTomog inverts for attenuation properties using t* observations. This report also describes how to use the tstarTomog package to go from observed data to a 3-D model of attenuation structure in the earth.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

A Rayleigh wave propagates laterally without dispersion in the vicinity of the plane stress-free surface of a homogeneous and isotropic elastic halfspace. The phase speed is independent of frequency and depends only on the Poisson ratio of the medium. However, after temporal and spatial discretization, a Rayleigh wave simulated by a 3D staggered-grid finite-difference (FD) seismic wave propagation algorithm suffers from frequency- and direction-dependent numerical dispersion. The magnitude of this dispersion depends critically on FD algorithm implementation details. Nevertheless, proper gridding can control numerical dispersion to within an acceptable level, leading to accurate Rayleigh wave simulations. Many investigators have derived dispersion relations appropriate for body wave propagation by various FD algorithms. However, the situation for surface waves is less well-studied. We have devised a numerical search procedure to estimate Rayleigh phase speed and group speed curves for 3D O(2,2) and O(2,4) staggered-grid FD algorithms. In contrast with the continuous time-space situation (where phase speed is obtained by extracting the appropriate root of the Rayleigh cubic), we cannot develop a closed-form mathematical formula governing the phase speed. Rather, we numerically seek the particular phase speed that leads to a solution of the discrete wave propagation equations, while holding medium properties, frequency, horizontal propagation direction, and gridding intervals fixed. Group speed is then obtained by numerically differentiating the phase speed with respect to frequency. The problem is formulated for an explicit stress-free surface positioned at two different levels within the staggered spatial grid. Additionally, an interesting variant involving zero-valued medium properties above the surface is addressed. We refer to the latter as an implicit free surface. Our preliminary conclusion is that an explicit free surface, implemented with O(4) spatial FD operators and positioned at the level of the compressional stress components, leads to superior numerical dispersion performance. Phase speeds measured from fixed-frequency synthetic seismograms agree very well with the numerical predictions.

An efficient numerical algorithm for treating earth models composed of fluid and solid portions is obtained via straightforward modifications to a 3D time-domain finite-difference algorithm for simulating isotropic elastic wave propagation.

Abstract not provided.

P- and S-body wave travel times collected from stations in and near the state of Nevada were inverted for P-wave velocity and the Vp/Vs ratio. These waves consist of Pn, Pg, Sn and Sg, but only the first arriving P and S waves were used in the inversion. Travel times were picked by University of Nevada Reno colleagues and were culled for inclusion in the tomographic inversion. The resulting tomographic model covers the entire state of Nevada to a depth of {approx}90 km; however, only the upper 40 km indicate relatively good resolution. Several features of interest are imaged including the Sierra Nevada, basin structures, and low velocities at depth below Yucca Mountain. These velocity structure images provide valuable information to aide in the interpretation of geothermal resource areas throughout the state on Nevada.

Abstract not provided.

Dispersion and attenuation relations are derived for both the continuous and discrete velocity-memory-stress systems governing 3D anelastic wave propagation in a standard linear solid. Phase speed and attenuation factor curves extracted from these relations enable optimal selection of spatial and temporal gridding intervals to achieve finite-difference algorithm efficiency, while simultaneously minimizing numerical inaccuracy.

Abstract not provided.

Stable and accurate numerical modeling of seismic wave propagation in the vicinity of high-contrast interfaces is achieved with straightforward modifications to the conventional, rectangular-staggered-grid, finite-difference (FD) method. Improvements in material parameter averaging and spatial differencing of wavefield variables yield high-quality synthetic seismic data.