Imaging Strong Lateral Heterogeneities with USArray using Body-to-Surface Wave Scattering

NASA Astrophysics Data System (ADS)

Yu, C.; Zhan, Z.; Hauksson, E.; Cochran, E. S.

2017-12-01

Seismic scattering is commonly observed and results from wave propagation in heterogeneous medium. Yet, deterministic characterization of scatterers remains challenging. In this study, we analyze broadband waveforms recorded by the USArray across the entire conterminous US. With array analysis, we observe strong scattered surface waves following the arrival of teleseismic body waves over several hundreds of kilometers. We use back-projection to locate the body-to-surface scattering sources, and detect strong scatterers both around and within the conterminous US. For the former, strong scattering is associated with pronounced bathymetric relief, such as the Patton Escarpment in the Southern California Continental Borderland. For the latter, scatterers are consistent with sharp lateral heterogeneities, such as near the Yellowstone hotspot and Southern California fault zones. We further model the body-to-surface wave scattering using finite-difference simulations. As an example, in the Southern California Continental Borderland a simplified 2-D bathymetric and crustal model are able to predict the arrival times and amplitudes of major scatterers. The modeling also suggests a relatively low shear wave velocity in the Continental Borderland. These observation of strong body-to-surface wave scattering and waveform modeling not only helps us image sharp heterogeneities but also are useful for assessing seismic hazard, including the calibration and refinement of seismic velocity models used to locate earthquakes and simulate strong ground motions.

A surface wave reflector in Southwestern Japan

NASA Astrophysics Data System (ADS)

Mak, S.; Koketsu, K.; Miyake, H.; Obara, K.; Sekine, S.

2009-12-01

Surface waves at short periods (<35s) are affected severely by heterogeneities in the crust and the uppermost mantle. When the scale of heterogeneity is sufficiently large, its effect can be studied in a deterministic way using conventional concepts of reflection and refraction. A well-known example is surface wave refraction at continental margin. We present a case study to investigate the composition of surface wave coda in a deterministic approach. A long duration of surface wave coda with a predominant period of 20s is observed during various strong earthquakes around Japan. The coda shows an unambiguous propagation direction, implying a deterministic nature. Beamforming and particle motion analysis suggest that the surface wave later arrivals could be explained by Love wave reflections by a point reflector located at offshore southeast to Kyushu. The reflection demonstrates a seemingly incidence-independent favorable azimuth in emitting strength. In additional to beamforming, we use a new regional crustal velocity model to perform a grid-search ray-tracing with the assumption of point reflector to further constrain to location of coda generation. Because strong velocity anomalies exist near the zone of interest, we decide to use a network shortest-path ray-tracing method, instead of analytical methods like shooting and bending, to avoid the problems like convergence, shadow zone, and smooth model assumption. Two geological features are found to be related to the formation of the coda. The primary one is the intersection between the Kyushu-Palau Ridge and the Nankai Trough at offshore southeast to Kyushu (hereafter referred as "KPR-NT"), which may act as a point reflector. There is a strong Love wave phase velocity anomaly at KPR-NT but not other parts of the ridge, implying that topography is irrelevant. Rayleigh wave phase velocity does not experience a strong anomaly there, which is consistent to the absence of Rayleigh wave reflections implied by the observed particle motions. The secondary one is a low phase velocity (<2km/s for T=20s) at the accretionary wedge of the Nankai Trough due to the thick sediment. Such a long and narrow low velocity zone, with its southwest tip at KPR-NT, is a potential wave-guide to channel waves towards KPR-NT. The longer duration of deterministic later arrivals than the direct arrival is partially explained by multi-pathing due to the wave-guide. The surface wave coda is observable for earthquakes whose propagation path does not include the accretionary wedge, implying that the wedge is an enhancer but not indispensable of the formation of the observed coda.

Seismic effects from major basin formation on the Moon and Mercury

NASA Technical Reports Server (NTRS)

Schultz, P. H.; Gault, D. E.

1974-01-01

Grooved and hilly terrains are reported which occur at the antipode of major basins on the Moon (Imbrium, Orientale) and Mercury (Caloris). Order-of-magnitude calculations, for an Imbrium-size impact on the Moon, indicate P-wave-induced surface displacements of 10 m at the basin antipode that would arrive prior to secondary ejecta. Comparable surface waves are reported which would arrive subsequent to secondary ejecta impacts and would increase in magnitude as they converge at the antipode. Other seismically induced surface features include: subdued, furrowed crater walls produced by landslides and concomitant secondary impacts; emplacement and leveling of light plains units owing to seismically induced "fluidization" of slide material; knobby, pitted terrain around old basins from enhancement of seismic waves in ancient ejecta blankets; and the production and enhancement of deep-seated fractures that led to the concentration of farside lunar maria in the Apollo-Ingenii region.

Seismic interferometry of railroad induced ground motions: body and surface wave imaging

NASA Astrophysics Data System (ADS)

Quiros, Diego A.; Brown, Larry D.; Kim, Doyeon

2016-04-01

Seismic interferometry applied to 120 hr of railroad traffic recorded by an array of vertical component seismographs along a railway within the Rio Grande rift has recovered surface and body waves characteristic of the geology beneath the railway. Linear and hyperbolic arrivals are retrieved that agree with surface (Rayleigh), direct and reflected P waves observed by nearby conventional seismic surveys. Train-generated Rayleigh waves span a range of frequencies significantly higher than those recovered from typical ambient noise interferometry studies. Direct P-wave arrivals have apparent velocities appropriate for the shallow geology of the survey area. Significant reflected P-wave energy is also present at relatively large offsets. A common midpoint stack produces a reflection image consistent with nearby conventional reflection data. We suggest that for sources at the free surface (e.g. trains) increasing the aperture of the array to record wide angle reflections, in addition to longer recording intervals, might allow the recovery of deeper geological structure from railroad traffic. Frequency-wavenumber analyses of these recordings indicate that the train source is symmetrical (i.e. approaching and receding) and that deeper refracted energy is present although not evident in the time-offset domain. These results confirm that train-generated vibrations represent a practical source of high-resolution subsurface information, with particular relevance to geotechnical and environmental applications.

Joint inversion of seismic and gravity data for imaging seismic velocity structure of the crust and upper mantle beneath Utah, United States

NASA Astrophysics Data System (ADS)

Syracuse, E. M.; Zhang, H.; Maceira, M.

2017-10-01

We present a method for using any combination of body wave arrival time measurements, surface wave dispersion observations, and gravity data to simultaneously invert for three-dimensional P- and S-wave velocity models. The simultaneous use of disparate data types takes advantage of the differing sensitivities of each data type, resulting in a comprehensive and higher resolution three-dimensional geophysical model. In a case study for Utah, we combine body wave first arrivals mainly from the USArray Transportable Array, Rayleigh wave group and phase velocity dispersion data, and Bouguer gravity anomalies to invert for crustal and upper mantle structure of the region. Results show clear delineations, visible in both P- and S-wave velocities, between the three main tectonic provinces in the region. Without the inclusion of the surface wave and gravity constraints, these delineations are less clear, particularly for S-wave velocities. Indeed, checkerboard tests confirm that the inclusion of the additional datasets dramatically improves S-wave velocity recovery, with more subtle improvements to P-wave velocity recovery, demonstrating the strength of the method in successfully recovering seismic velocity structure from multiple types of constraints.

Impacts of short-time scale water column variability on broadband high-frequency acoustic wave propagation

NASA Astrophysics Data System (ADS)

Eickmeier, Justin

Acoustical oceanography is one way to study the ocean, its internal layers, boundaries and all processes occurring within using underwater acoustics. Acoustical sensing techniques allows for the measurement of ocean processes from within that logistically or financially preclude traditional in-situ measurements. Acoustic signals propagate as pressure wavefronts from a source to a receiver through an ocean medium with variable physical parameters. The water column physical parameters that change acoustic wave propagation in the ocean include temperature, salinity, current, surface roughness, seafloor bathymetry, and vertical stratification over variable time scales. The impacts of short-time scale water column variability on acoustic wave propagation include coherent and incoherent surface reflections, wavefront arrival time delay, focusing or defocusing of the intensity of acoustic beams and refraction of acoustic rays. This study focuses on high-frequency broadband acoustic waves, and examines the influence of short-time scale water column variability on broadband high-frequency acoustics, wavefronts, from 7 to 28 kHz, in shallow water. Short-time scale variability is on the order of seconds to hours and the short-spatial scale variability is on the order of few centimeters. Experimental results were collected during an acoustic experiment along 100 m isobaths and data analysis was conducted using available acoustic wave propagation models. Three main topics are studied to show that acoustic waves are viable as a remote sensing tool to measure oceanographic parameters in shallow water. First, coherent surface reflections forming striation patterns, from multipath receptions, through rough surface interaction of broadband acoustic signals with the dynamic sea surface are analyzed. Matched filtered results of received acoustic waves are compared with a ray tracing numerical model using a sea surface boundary generated from measured water wave spectra at the time of signal propagation. It is determined that on a time scale of seconds, corresponding to typical periods of surface water waves, the arrival time of reflected acoustic signals from surface waves appear as striation patterns in measured data and can be accurately modelled by ray tracing. Second, changes in acoustic beam arrival angle and acoustic ray path influenced by isotherm depth oscillations are analyzed using an 8-element delay-sum beamformer. The results are compared with outputs from a two-dimensional (2-D) parabolic equation (PE) model using measured sound speed profiles (SSPs) in the water column. Using the method of beamforming on the received signal, the arrival time and angle of an acoustic beam was obtained for measured acoustic signals. It is determined that the acoustic ray path, acoustic beam intensity and angular spread are a function of vertical isotherm oscillations on a time scale of minutes and can be modeled accurately by a 2-D PE model. Third, a forward problem is introduced which uses acoustic wavefronts received on a vertical line array, 1.48 km from the source, in the lower part of the water column to infer range dependence or independence in the SSP. The matched filtering results of received acoustic wavefronts at all hydrophone depths are compared with a ray tracing routine augmented to calculate only direct path and bottom reflected signals. It is determined that the SSP range dependence can be inferred on a time scale of hours using an array of hydrophones spanning the water column. Sound speed profiles in the acoustic field were found to be range independent for 11 of the 23 hours in the measurements. A SSP cumulative reconstruction process, conducted from the seafloor to the sea surface, layer-by-layer, identifies critical segments in the SSP that define the ray path, arrival time and boundary interactions. Data-model comparison between matched filtered arrival time spread and arrival time output from the ray tracing was robust when the SSP measured at the receiver was input to the model. When the SSP measured nearest the source (at the same instant in time) was input to the ray tracing model, the data-model comparison was poor. It was determined that the cumulative sound speed change in the SSP near the source was 1.041 m/s greater than that of the SSP at the receiver and resulted in the poor data-model comparison. In this study, the influences on broadband acoustic wave propagation in the frequency range of 7 to 28 kHz of spatial and temporal changes in the oceanography of shallow water regions are addressed. Acoustic waves can be used as remote sensing tools to measure oceanographic parameters in shallow water and data-model comparison results show a direct relationship between the oceanographic variations and acoustic wave propagations.

Rapid Estimation of Building Damage by Conventional Weapons

DTIC Science & Technology

2014-09-01

impulse applied to reflect surface 17 iso = impulse of incident wave t0 = wave duration Lw = wavelength of incident wave ta = arrival time...floor area = 3 x 100 x 90 = 27000 ft2 Tributary area support per column = 27000 /30 = 900 ft2 Therefore, the tributary area supported by a column is

Acoustic resonances of fluid-immersed elastic cylinders and spheroids: Theory and experiment

NASA Astrophysics Data System (ADS)

Niemiec, Jan; Überall, Herbert; Bao, X. L.

2002-05-01

Frequency resonances in the scattering of acoustic waves from a target object are caused by the phase matching of surface waves repeatedly encircling the object. This is exemplified here by considering elastic finite cylinders and spheroids, and the phase-matching condition provides a means of calculating the complex resonance frequencies of such objects. Tank experiments carried out at Catholic University, or at the University of Le Havre, France by G. Maze and J. Ripoche, have been interpreted using this approach. The experiments employed sound pulses to measure arrival times, which allowed identification of the surface paths taken by the surface waves, thus giving rise to resonances in the scattering amplitude. A calculation of the resonance frequencies using the T-matrix approach showed satisfactory agreement with the experimental resonance frequencies that were either measured directly (as at Le Havre), or that were obtained by the interpretation of measured arrival times (at Catholic University) using calculated surface wave paths, and the extraction of resonance frequencies therefrom, on the basis of the phase-matching condition. Results for hemispherically endcapped, evacuated steel cylinders obtained in a lake experiment carried out by the NSWC were interpreted in the same fashion.

Joint inversion of seismic and gravity data for imaging seismic velocity structure of the crust and upper mantle beneath Utah, United States

DOE PAGES

Syracuse, Ellen Marie; Zhang, Haijiang; Maceira, Monica

2017-07-11

Here, we present a method for using any combination of body wave arrival time measurements, surface wave dispersion observations, and gravity data to simultaneously invert for three-dimensional P- and S-wave velocity models. The simultaneous use of disparate data types takes advantage of the differing sensitivities of each data type, resulting in a comprehensive and higher resolution three-dimensional geophysical model. In a case study for Utah, we combine body waves first arrivals mainly from the USArray Transportable Array, Rayleigh wave group and phase velocity dispersion data, and Bouguer gravity anomalies to invert for crustal and upper mantle structure of the region.more » Results show clear delineations, visible in both P- and S-wave velocities, between the three main tectonic provinces in the region. In conclusion, without the inclusion of the surface wave and gravity constraints, these delineations are less clear, particularly for S-wave velocities. Indeed, checkerboard tests confirm that the inclusion of the additional datasets dramatically improves S-wave velocity recovery, with more subtle improvements to P-wave velocity recovery, demonstrating the strength of the method in successfully recovering seismic velocity structure from multiple types of constraints.« less

Joint inversion of seismic and gravity data for imaging seismic velocity structure of the crust and upper mantle beneath Utah, United States

DOE Office of Scientific and Technical Information (OSTI.GOV)

Syracuse, Ellen Marie; Zhang, Haijiang; Maceira, Monica

Here, we present a method for using any combination of body wave arrival time measurements, surface wave dispersion observations, and gravity data to simultaneously invert for three-dimensional P- and S-wave velocity models. The simultaneous use of disparate data types takes advantage of the differing sensitivities of each data type, resulting in a comprehensive and higher resolution three-dimensional geophysical model. In a case study for Utah, we combine body waves first arrivals mainly from the USArray Transportable Array, Rayleigh wave group and phase velocity dispersion data, and Bouguer gravity anomalies to invert for crustal and upper mantle structure of the region.more » Results show clear delineations, visible in both P- and S-wave velocities, between the three main tectonic provinces in the region. In conclusion, without the inclusion of the surface wave and gravity constraints, these delineations are less clear, particularly for S-wave velocities. Indeed, checkerboard tests confirm that the inclusion of the additional datasets dramatically improves S-wave velocity recovery, with more subtle improvements to P-wave velocity recovery, demonstrating the strength of the method in successfully recovering seismic velocity structure from multiple types of constraints.« less

Earthquake site response in Santa Cruz, California

USGS Publications Warehouse

Carver, D.; Hartzell, S.H.

1996-01-01

Aftershocks of the 1989 Loma Prieta, California, earthquake are used to estimate site response in a 12-km2 area centered on downtown Santa Cruz. A total of 258 S-wave records from 36 aftershocks recorded at 33 sites are used in a linear inversion for site-response spectra. The inversion scheme takes advantage of the redundancy of the large data set for which several aftershocks are recorded at each site. The scheme decomposes the observed spectra into source, path, and site terms. The path term is specified before the inversion. The undetermined degree of freedom in the decomposition into source and site spectra is removed by specifying the site-response factor to be approximately 1.0 at two sites on crystalline bedrock. The S-wave site responses correlate well with the surficial geology and observed damage pattern of the mainshock. The site-response spectra of the floodplain sites, which include the heavily damaged downtown area, exhibit significant peaks. The largest peaks are between 1 and 4 Hz. Five floodplain sites have amplification factors of 10 or greater. Most of the floodplain site-response spectra also have a smaller secondary peak between 6 and 8 Hz. Residential areas built on marine terraces above the flood-plain experienced much less severe damage. Site-response spectra for these areas also have their largest peaks between 1 and 4 Hz, but the amplification is generally below 6. Several of these sites also have a secondary peak between 6 and 8 Hz. The response peaks seen at nearly all sites between 1 and 4 Hz are probably caused by the natural resonance of the sedimentary rock column. The higher amplifications at floodplain sites may be caused by surface waves generated at the basin margins. The secondary peak between 6 and 8 Hz at many sites may be a harmonic of the 1- to 4-Hz peaks. We used waveforms from a seven-station approximately linear array located on the floodplain to calculate the apparent velocity and azimuth of propagation of coherent arrivals within moving windows of the S-wave codas. The initial windows give results that are consistent with direct S-wave arrivals. The apparent velocities are high (greater than 4.0 km/sec), and azimuths are from the source. Waves arriving later than 2 sec after the direct S waves have apparent velocities of less than 1 km/ sec, indicating that they are surface waves, and arrive from divergent azimuths. This analysis indicates that after the direct S-wave arrival, long-duration shaking comes from surface waves that are generated at the basin margin and reverberate in the floodplain sediments.

Triggered Seismicity in Utah from the November 3, 2002, Denali Fault Earthquake

NASA Astrophysics Data System (ADS)

Pankow, K. L.; Nava, S. J.; Pechmann, J. C.; Arabasz, W. J.

2002-12-01

Coincident with the arrival of the surface waves from the November 3, 2002, Mw 7.9 Denali Fault, Alaska earthquake (DFE), the University of Utah Seismograph Stations (UUSS) regional seismic network detected a marked increase in seismicity along the Intermountain Seismic Belt (ISB) in central and north-central Utah. The number of earthquakes per day in Utah located automatically by the UUSS's Earthworm system in the week following the DFE was approximately double the long-term average during the preceding nine months. From these preliminary data, the increased seismicity appears to be characterized by small magnitude events (M = 3.2) and concentrated in five distinct spatial clusters within the ISB between 38.75°and 42.0° N. The first of these earthquakes was an M 2.2 event located ~20 km east of Salt Lake City, Utah, which occurred during the arrival of the Love waves from the DFE. The increase in Utah earthquake activity at the time of the arrival of the surface waves from the DFE suggests that these surface waves triggered earthquakes in Utah at distances of more than 3,000 km from the source. We estimated the peak dynamic shear stress caused by these surface waves from measurements of their peak vector velocities at 43 recording sites: 37 strong-motion stations of the Advanced National Seismic System and six broadband stations. (The records from six other broadband instruments in the region of interest were clipped.) The estimated peak stresses ranged from 1.2 bars to 3.5 bars with a mean of 2.3 bars, and generally occurred during the arrival of Love waves of ~15 sec period. These peak dynamic shear stress estimates are comparable to those obtained from recordings of the 1992 Mw 7.3 Landers, California, earthquake in regions where the Landers earthquake triggered increased seismicity. We plan to present more complete analyses of UUSS seismic network data, further testing our hypothesis that the DFE remotely triggered seismicity in Utah. This hypothesis is important to investigate because well-documented evidence for triggering of seismicity by distant earthquakes comes primarily from areas characterized by recent volcanic or geothermal activity. The regions of apparent triggered seismicity from the DFE in Utah fall into neither of these two categories.

A repeating source of infrasound from the Wells, Nevada earthquake sequence

DOE Office of Scientific and Technical Information (OSTI.GOV)

Arrowsmith, Stephen J.; Whitaker, Rod; Randall, George

2009-01-01

The Wells, Nevada earthquake of February 21, 2008, generated a complex seismoacoustic wakefield. In addition to epicentral infrasound, the earthquake triggered a secondary source of infrasound, which was also initiated by subsequent aftershocks. By applying simple constraints on the propagation of seismic and infrasound waves, we show that the secondary source is an isolated peak that appears to efficiently generate infrasound through the interaction with seismic surface waves. By measuring peak-to-peak amplitudes of epicentral and secondary arrivals and correcting them for the effects of distance and winds, we find that epicentral arrivals lit with empirical relationships of Mutschlecner and Whitakermore » (2005) and Le Pichon et al. (2006), which form the basis for a proposed infrasound discriminant (Anderson et al., Pers. Comm.). In contrast, the secondary arrivals are much higher in amplitude, highlighting the importance of being able to separate epicentral and secondary arrivals for infrasonic event discrimination.« less

Multichannel analysis of surface-waves and integration of downhole acoustic televiewer imaging, ultrasonic Vs and Vp, and vertical seismic profiling in an NEHRP-standard classification, South of Concordia, Kansas, USA

NASA Astrophysics Data System (ADS)

Raef, Abdelmoneam; Gad, Sabreen; Tucker-Kulesza, Stacey

2015-10-01

Seismic site characteristics, as pertaining to earthquake hazard reduction, are a function of the subsurface elastic moduli and the geologic structures. This study explores how multiscale (surface, downhole, and laboratory) datasets can be utilized to improve "constrained" average Vs30 (shear-wave velocity to a 30-meter depth). We integrate borehole, surface and laboratory measurements for a seismic site classification based on the standards of the National Earthquake Hazard Reduction Program (NEHRP). The seismic shear-wave velocity (Vs30) was derived from a geophysical inversion workflow that utilized multichannel analysis of surface-waves (MASW) and downhole acoustic televiewer imaging (DATI). P-wave and S-wave velocities, based on laboratory measurements of arrival times of ultrasonic-frequency signals, supported the workflow by enabling us to calculate Poisson's ratio, which was incorporated in building an initial model for the geophysical inversion of MASW. Extraction of core samples from two boreholes provided lithology and thickness calibration of the amplitudes of the acoustic televiewer imaging for each layer. The MASW inversion, for calculating Vs sections, was constrained with both ultrasonic laboratory measurements (from first arrivals of Vs and Vp waveforms at simulated in situ overburden stress conditions) and the downhole acoustic televiewer (DATV) amplitude logs. The Vs30 calculations enabled categorizing the studied site as NEHRP-class "C" - very dense soil and soft rock. Unlike shallow fractured carbonates in the studied area, S-wave and P-wave velocities at ultrasonic frequency for the deeper intact shale core-samples from two boreholes were in better agreement with the corresponding velocities from both a zero-offset vertical seismic profiling (VSP) and inversion of Rayleigh-wave velocity dispersion curves.

Dense Seismic Recordings of Two Surface-Detonated Chemical Explosions

NASA Astrophysics Data System (ADS)

Koper, K. D.; Hale, J. M.; Burlacu, R.; Goddard, K. J.; Trow, A.; Linville, L. M.; Stein, J. R.; Drobeck, D.; Leidig, M.

2015-12-01

In the summer of 2015 two controlled chemical explosions were carried out near Dugway, Utah. The 2 June 2015 explosion consisted of 30,000 lbs of ammonium nitrate fuel oil (ANFO) and the 22 July 2015 explosion consisted of 60,000 lbs of ANFO. The explosion centroids were 1-2 m above the Earth's surface and both created significant craters in the soft desert alluvium. To better understand the seismic source associated with surface explosions, we deployed an array of wireless, three-component, short-period (5 Hz corner frequency) seismometers for several days around each shot. For the first explosion, 46 receivers were deployed in a "lollipop" geometry that had a sparse ring at a radius of 1 km, and a dense stem with 100 m spacing for distances of 0.5-4.5 km. For the second explosion, 48 receivers were deployed similarly, but with a dense ring spaced in azimuthal increments of 10 degrees at a distance of 1 km, and a sparse stem (~500 m spacing) that extended to a distance of nearly 6 km. A rich variety of phases were recorded including direct P waves, refracted and reflected P waves, nearly monochromatic air-coupled Rayleigh waves, normally dispersed fundamental mode Rayleigh waves (Rg), primary airblast arrivals, some secondary airblast arrivals, and possibly tertiary airblast arrivals. There is also evidence of converted S waves on the radial components and possibly direct S energy on the radial and transverse components, although the transverse energy does not always possess a simple, coherent move-out with distance, implying that it might have a scattering origin. To aid in the phase identification, especially of the apparent SH and Love energy, we are currently performing tau-p, f-k, and particle motion analysis.

Joint Inversion of Body-Wave Arrival Times and Surface-Wave Dispersion Data in the Wavelet Domain Constrained by Sparsity Regularization

NASA Astrophysics Data System (ADS)

Zhang, H.; Fang, H.; Yao, H.; Maceira, M.; van der Hilst, R. D.

2014-12-01

Recently, Zhang et al. (2014, Pure and Appiled Geophysics) have developed a joint inversion code incorporating body-wave arrival times and surface-wave dispersion data. The joint inversion code was based on the regional-scale version of the double-difference tomography algorithm tomoDD. The surface-wave inversion part uses the propagator matrix solver in the algorithm DISPER80 (Saito, 1988) for forward calculation of dispersion curves from layered velocity models and the related sensitivities. The application of the joint inversion code to the SAFOD site in central California shows that the fault structure is better imaged in the new model, which is able to fit both the body-wave and surface-wave observations adequately. Here we present a new joint inversion method that solves the model in the wavelet domain constrained by sparsity regularization. Compared to the previous method, it has the following advantages: (1) The method is both data- and model-adaptive. For the velocity model, it can be represented by different wavelet coefficients at different scales, which are generally sparse. By constraining the model wavelet coefficients to be sparse, the inversion in the wavelet domain can inherently adapt to the data distribution so that the model has higher spatial resolution in the good data coverage zone. Fang and Zhang (2014, Geophysical Journal International) have showed the superior performance of the wavelet-based double-difference seismic tomography method compared to the conventional method. (2) For the surface wave inversion, the joint inversion code takes advantage of the recent development of direct inversion of surface wave dispersion data for 3-D variations of shear wave velocity without the intermediate step of phase or group velocity maps (Fang et al., 2014, Geophysical Journal International). A fast marching method is used to compute, at each period, surface wave traveltimes and ray paths between sources and receivers. We will test the new joint inversion code at the SAFOD site to compare its performance over the previous code. We will also select another fault zone such as the San Jacinto Fault Zone to better image its structure.

Velocity Structure Determination Through Seismic Waveform Modeling and Time Deviations

NASA Astrophysics Data System (ADS)

Savage, B.; Zhu, L.; Tan, Y.; Helmberger, D. V.

2001-12-01

Through the use of seismic waveforms recorded by TriNet, a dataset of earthquake focal mechanisms and deviations (time shifts) relative to a standard model facilitates the investigation of the crust and uppermost mantle of southern California. The CAP method of focal mechanism determination, in use by TriNet on a routine basis, provides time shifts for surface waves and Pnl arrivals independently relative to the reference model. These shifts serve as initial data for calibration of local and regional seismic paths. Time shifts from the CAP method are derived by splitting the Pnl section of the waveform, the first arriving Pn to just before the arrival of the S wave, from the much slower surface waves then cross-correlating the data with synthetic waveforms computed from a standard model. Surface waves interact with the entire crust, but the upper crust causes the greatest effect. Whereas, Pnl arrivals sample the deeper crust, upper mantle, and source region. This natural division separates the upper from lower crust for regional calibration and structural modeling and allows 3-D velocity maps to be created using the resulting time shifts. Further examination of Pnl and other arrivals which interact with the Moho illuminate the complex nature of this boundary. Initial attempts at using the first 10 seconds of the Pnl section to determine upper most mantle structure have proven insightful. Two large earthquakes north of southern California in Nevada and Mammoth Lakes, CA allow the creation of record sections from 200 to 600 km. As the paths swing from east to west across southern California, simple 1-D models turn into complex structure, dramatically changing the waveform character. Using finite difference models to explain the structure, we determine that a low velocity zone is present at the base of the crust and extends to 100 km in depth. Velocity variations of 5 percent of the mantle in combination with steeply sloping edges produces complex waveform variations. Characteristics of this complex propagation appear from the southern Sierra Nevada Mountains, in the west, to Death Valley in the east. The structure does not cross the Garlock fault to the south, but we are unsure of the structures northern extent.

A practical implementation of wave front construction for 3-D isotropic media

NASA Astrophysics Data System (ADS)

Chambers, K.; Kendall, J.-M.

2008-06-01

Wave front construction (WFC) methods are a useful tool for tracking wave fronts and are a natural extension to standard ray shooting methods. Here we describe and implement a simple WFC method that is used to interpolate wavefield properties throughout a 3-D heterogeneous medium. Our approach differs from previous 3-D WFC procedures primarily in the use of a ray interpolation scheme, based on approximating the wave front as a `locally spherical' surface and a `first arrival mode', which reduces computation times, where only first arrivals are required. Both of these features have previously been included in 2-D WFC algorithms; however, until now they have not been extended to 3-D systems. The wave front interpolation scheme allows for rays to be traced from a nearly arbitrary distribution of take-off angles, and the calculation of derivatives with respect to take-off angles is not required for wave front interpolation. However, in regions of steep velocity gradient, the locally spherical approximation is not valid, and it is necessary to backpropagate rays to a sufficiently homogenous region before interpolation of the new ray. Our WFC technique is illustrated using a realistic velocity model, based on a North Sea oil reservoir. We examine wavefield quantities such as traveltimes, ray angles, source take-off angles and geometrical spreading factors, all of which are interpolated on to a regular grid. We compare geometrical spreading factors calculated using two methods: using the ray Jacobian and by taking the ratio of a triangular area of wave front to the corresponding solid angle at the source. The results show that care must be taken when using ray Jacobians to calculate geometrical spreading factors, as the poles of the source coordinate system produce unreliable values, which can be spread over a large area, as only a few initial rays are traced in WFC. We also show that the use of the first arrival mode can reduce computation time by ~65 per cent, with the accuracy of the interpolated traveltimes, ray angles and source take-off angles largely unchanged. However, the first arrival mode does lead to inaccuracies in interpolated angles near caustic surfaces, as well as small variations in geometrical spreading factors for ray tubes that have passed through caustic surfaces.

Integration of an Acoustic Modem onto a Wave Glider Unmanned Surface Vehicle

DTIC Science & Technology

2012-06-01

of the wave and ωτ represents the phase of the wave. After some amount of math and taking the limit as ω →∞ , we arrive at a form of the eikonal ...the phase front. (5.5) 22 0A Aτ τ∇ ⋅∇ + ∇ = The transport equation and the eikonal equation can be solved by using multiple methods to give

Thermal-hydraulic behaviors of vapor-liquid interface due to arrival of a pressure wave

DOE Office of Scientific and Technical Information (OSTI.GOV)

Inoue, Akira; Fujii, Yoshifumi; Matsuzaki, Mitsuo

In the vapor explosion, a pressure wave (shock wave) plays a fundamental role for triggering, propagation and enhancement of the explosion. Energy of the explosion is related to the magnitude of heat transfer rate from hot liquid to cold volatile one. This is related to an increasing rate of interface area and to an amount of transient heat flux between the liquids. In this study, the characteristics of transient heat transfer and behaviors of vapor film both on the platinum tube and on the hot melt tin drop, under same boundary conditions have been investigated. It is considered that theremore » exists a fundamental mechanism of the explosion in the initial expansion process of the hot liquid drop immediately after arrival of pressure wave. The growth rate of the vapor film is much faster on the hot liquid than that on the solid surface. Two kinds of roughness were observed, one due to the Taylor instability, by rapid growth of the explosion bubble, and another, nucleation sites were observed at the vapor-liquid interface. Based on detailed observation of early stage interface behaviors after arrival of a pressure wave, the thermal fragmentation mechanism is proposed.« less

Assessing Multiple Methods for Determining Active Source Travel Times in a Dense Array

NASA Astrophysics Data System (ADS)

Parker, L.; Zeng, X.; Thurber, C. H.; Team, P.

2016-12-01

238 three-component nodal seismometers were deployed at the Brady Hot Springs geothermal field in Nevada to characterize changes in the subsurface as a result of changes in pumping conditions. The array consisted of a 500 meter by 1600 meter irregular grid with 50 meter spacing centered in an approximately rectangular 1200 meter by 1600 meter grid with 200 meter spacing. A large vibroseis truck (T-Rex) was deployed as an active seismic source at 216 locations. Over the course of 15 days, the truck occupied each location up to four times. At each location a swept-frequency source between 5 and 80 Hz over 20 seconds was produced using three vibration modes: longitudinal S-wave, transverse S-wave, and P-wave. Seismic wave arrivals were identified using three methods: cross-correlation, deconvolution, and Wigner-Ville distribution (WVD) plus the Hough Transform (HT). Surface wave arrivals were clear for all three modes of vibration using all three methods. Preliminary tomographic models will be presented, using the arrivals of the identified phases. This analysis is part of the PoroTomo project: Poroelastic Tomography by Adjoint Inverse Modeling of Data from Seismology, Geodesy, and Hydrology; http://geoscience.wisc.edu/feigl/porotomo.

Sub-micron elastic property characterization of materials using a near-field scanning optical microscope

NASA Astrophysics Data System (ADS)

Blodgett, David W.; Spicer, James B.

2001-12-01

The ability to characterize the sub-surface mechanical properties of a bulk or thin film material at the sub-micron level has applications in the microelectronics and thin film industries. In the microelectronics industry, with the decrease of line widths and the increase of component densities, sub-surface voids have become increasingly detrimental. Any voids along an integrated circuit (IC) line can lead to improper electrical connections between components and can cause failure of the device. In the thin film industry, the detection of impurities is also important. Any impurities can detract from the film's desired optical, electrical, or mechanical properties. Just as important as the detection of voids and impurities, is the measurement of the elastic properties of a material on the nanometer scale. These elastic measurements provide insight into the microstructural properties of the material. We have been investigating a technique that couples the high-resolution surface imaging capabilities of the apertureless near-field scanning optical microscope (ANSOM) with the sub-surface characterization strengths of high-frequency ultrasound. As an ultrasonic wave propagates, the amplitude decreases due to geometrical spreading, attenuation from absorption, and scattering from discontinuities. Measurement of wave speeds and attenuation provides the information needed to quantify the bulk or surface properties of a material. The arrival of an ultrasonic wave at or along the surface of a material is accompanied with a small surface displacement. Conventional methods for the ultrasound detection rely on either a contact transducer or optical technique (interferometric, beam deflection, etc.). However, each of these methods is limited by the spatial resolution dictated by the detection footprint. As the footprint size increases, variations across the ultrasonic wavefront are effectively averaged, masking the presence of any nanometer-scale sub-surface or surface mechanical property variations. The use of an ANSOM for sensing ultrasonic wave arrivals reduces the detection footprint allowing any nanometer scale variations in the microstructure of a material to be detected. In an ANSOM, the ultrasonic displacement is manifested as perturbations on the near-field signal due to the small variations in the tip-sample caused by the wave arrival. Due to the linear dependence of the near-field signal on tip-sample separation, these perturbations can be interpreted using methods identical to those for conventional ultrasonic techniques. In this paper, we report results using both contact transducer (5 MHz) and laser-generated ultrasound.

Joint Inversion of Body-Wave Arrival Times and Surface-Wave Dispersion Data for Three-Dimensional Seismic Velocity Structure Around SAFOD

NASA Astrophysics Data System (ADS)

Zhang, H.; Thurber, C. H.; Maceira, M.; Roux, P.

2013-12-01

The crust around the San Andreas Fault Observatory at depth (SAFOD) has been the subject of many geophysical studies aimed at characterizing in detail the fault zone structure and elucidating the lithologies and physical properties of the surrounding rocks. Seismic methods in particular have revealed the complex two-dimensional (2D) and three-dimensional (3D) structure of the crustal volume around SAFOD and the strong velocity reduction in the fault damage zone. In this study we conduct a joint inversion using body-wave arrival times and surface-wave dispersion data to image the P-and S-wave velocity structure of the upper crust surrounding SAFOD. The two data types have complementary strengths - the body-wave data have good resolution at depth, albeit only where there are crossing rays between sources and receivers, whereas the surface waves have very good near-surface resolution and are not dependent on the earthquake source distribution because they are derived from ambient noise. The body-wave data are from local earthquakes and explosions, comprising the dataset analyzed by Zhang et al. (2009). The surface-wave data are for Love waves from ambient noise correlations, and are from Roux et al. (2011). The joint inversion code is based on the regional-scale version of the double-difference (DD) tomography algorithm tomoDD. The surface-wave inversion code that is integrated into the joint inversion algorithm is from Maceira and Ammon (2009). The propagator matrix solver in the algorithm DISPER80 (Saito, 1988) is used for the forward calculation of dispersion curves from layered velocity models. We examined how the structural models vary as we vary the relative weighting of the fit to the two data sets and in comparison to the previous separate inversion results. The joint inversion with the 'optimal' weighting shows more clearly the U-shaped local structure from the Buzzard Canyon Fault on the west side of SAF to the Gold Hill Fault on the east side.

Explosion source strong ground motions in the Mississippi embayment

USGS Publications Warehouse

Langston, C.A.; Bodin, P.; Powell, C.; Withers, M.; Horton, S.; Mooney, W.

2006-01-01

Two strong-motion arrays were deployed for the October 2002 Embayment Seismic Excitation Experiment to study the spatial variation of strong ground motions in the deep, unconsolidated sediments of the Mississippi embayment because there are no comparable strong-motion data from natural earthquakes in the area. Each linear array consisted of eight three-component K2 accelerographs spaced 15 m apart situated 1.2 and 2.5 kin from 2268-kg and 1134-kg borehole explosion sources, respectively. The array data show distinct body-wave and surface-wave arrivals that propagate within the thick, unconsolidated sedimentary column, the high-velocity basement rocks, and small-scale structure near the surface. Time-domain coherence of body-wave and surface-wave arrivals is computed for acceleration, velocity, and displacement time windows. Coherence is high for relatively low-frequency verticalcomponent Rayleigh waves and high-frequency P waves propagating across the array. Prominent high-frequency PS conversions seen on radial components, a proxy for the direct S wave from earthquake sources, lose coherence quickly over the 105-m length of the array. Transverse component signals are least coherent for any ground motion and appear to be highly scattered. Horizontal phase velocity is computed by using the ratio of particle velocity to estimates of the strain based on a plane-wave-propagation model. The resulting time-dependent phase-velocity map is a useful way to infer the propagation mechanisms of individual seismic phases and time windows of three-component waveforms. Displacement gradient analysis is a complementary technique for processing general spatial-array data to obtain horizontal slowness information.

Estimates of velocity structure and source depth using multiple P waves from aftershocks of the 1987 Elmore Ranch and Superstition Hills, California, earthquakes

USGS Publications Warehouse

Mori, J.

1991-01-01

Event record sections, which are constructed by plotting seismograms from many closely spaced earthquakes recorded on a few stations, show multiple free-surface reflections (PP, PPP, PPPP) of the P wave in the Imperial Valley. The relative timing of these arrivals is used to estimate the strength of the P-wave velocity gradient within the upper 5 km of the sediment layer. Consistent with previous studies, a velocity model with a value of 1.8 km/sec at the surface increasing linearly to 5.8 km/sec at a depth of 5.5 km fits the data well. The relative amplitudes of the P and PP arrivals are used to estimate the source depth for the aftershock distributions of the Elmore Ranch and Superstition Hills main shocks. Although the depth determination has large uncertainties, both the Elmore Ranch and Superstition Hills aftershock sequencs appear to have similar depth distribution in the range of 4 to 10 km. -Author

Air-ground interface: Surface waves, surface impedance and acoustic-to-seismic coupling coefficient

NASA Technical Reports Server (NTRS)

Daigle, Gilles; Embleton, Tony

1990-01-01

In atmospheric acoustics, the subject of surface waves has been an area of discussion for many years. The existence of an acoustic surface wave is now well established theoretically. The mathematical solution for spherical wave propagation above an impedance boundary includes the possibility of a contribution that possesses all the standard properties for a surface wave. Surface waves exist when the surface is sufficiently porous, relative to its acoustical resistance, that it can influence the airborne particle velocity near the surface and reduce the phase velocity of sound waves in air at the surface. This traps some of the sound energy in the air to remain near the surface as it propagates. Above porous grounds, the existence of surface waves has eluded direct experimental confirmation (pulse experiments have failed to show a separate arrival expected from the reduced phase speed) and indirect evidence for its existence has appeared contradictory. The experimental evidence for the existence of an acoustical surface wave above porous boundaries is reviewed. Recent measurements including pulse experiments are also described. A few years ago the acoustic impedance of a grass-covered surface was measured in the frequency range 30 to 300 Hz. Here, further measurements on the same site are discussed. These measurements include core samples, a shallow refractive survey to determine the seismic velocities, and measurements of the acoustic-to-seismic coupling coefficient.

Multichannel analysis of surface waves

USGS Publications Warehouse

Park, C.B.; Miller, R.D.; Xia, J.

1999-01-01

The frequency-dependent properties of Rayleigh-type surface waves can be utilized for imaging and characterizing the shallow subsurface. Most surface-wave analysis relies on the accurate calculation of phase velocities for the horizontally traveling fundamental-mode Rayleigh wave acquired by stepping out a pair of receivers at intervals based on calculated ground roll wavelengths. Interference by coherent source-generated noise inhibits the reliability of shear-wave velocities determined through inversion of the whole wave field. Among these nonplanar, nonfundamental-mode Rayleigh waves (noise) are body waves, scattered and nonsource-generated surface waves, and higher-mode surface waves. The degree to which each of these types of noise contaminates the dispersion curve and, ultimately, the inverted shear-wave velocity profile is dependent on frequency as well as distance from the source. Multichannel recording permits effective identification and isolation of noise according to distinctive trace-to-trace coherency in arrival time and amplitude. An added advantage is the speed and redundancy of the measurement process. Decomposition of a multichannel record into a time variable-frequency format, similar to an uncorrelated Vibroseis record, permits analysis and display of each frequency component in a unique and continuous format. Coherent noise contamination can then be examined and its effects appraised in both frequency and offset space. Separation of frequency components permits real-time maximization of the S/N ratio during acquisition and subsequent processing steps. Linear separation of each ground roll frequency component allows calculation of phase velocities by simply measuring the linear slope of each frequency component. Breaks in coherent surface-wave arrivals, observable on the decomposed record, can be compensated for during acquisition and processing. Multichannel recording permits single-measurement surveying of a broad depth range, high levels of redundancy with a single field configuration, and the ability to adjust the offset, effectively reducing random or nonlinear noise introduced during recording. A multichannel shot gather decomposed into a swept-frequency record allows the fast generation of an accurate dispersion curve. The accuracy of dispersion curves determined using this method is proven through field comparisons of the inverted shear-wave velocity (??(s)) profile with a downhole ??(s) profile.Multichannel recording is an efficient method of acquiring ground roll. By displaying the obtained information in a swept-frequency format, different frequency components of Rayleigh waves can be identified by distinctive and simple coherency. In turn, a seismic surface-wave method is derived that provides a useful noninvasive tool, where information about elastic properties of near-surface materials can be effectively obtained.

Improving arrival time identification in transient elastography

NASA Astrophysics Data System (ADS)

Klein, Jens; McLaughlin, Joyce; Renzi, Daniel

2012-04-01

In this paper, we improve the first step in the arrival time algorithm used for shear wave speed recovery in transient elastography. In transient elastography, a shear wave is initiated at the boundary and the interior displacement of the propagating shear wave is imaged with an ultrasound ultra-fast imaging system. The first step in the arrival time algorithm finds the arrival times of the shear wave by cross correlating displacement time traces (the time history of the displacement at a single point) with a reference time trace located near the shear wave source. The second step finds the shear wave speed from the arrival times. In performing the first step, we observe that the wave pulse decorrelates as it travels through the medium, which leads to inaccurate estimates of the arrival times and ultimately to blurring and artifacts in the shear wave speed image. In particular, wave ‘spreading’ accounts for much of this decorrelation. Here we remove most of the decorrelation by allowing the reference wave pulse to spread during the cross correlation. This dramatically improves the images obtained from arrival time identification. We illustrate the improvement of this method on phantom and in vivo data obtained from the laboratory of Mathias Fink at ESPCI, Paris.

The incorporation of fault zone head wave and direct wave secondary arrival times and arrival polarizations into seismic tomography: Application to the Parkfield, California area

NASA Astrophysics Data System (ADS)

Bennington, N. L.; Thurber, C. H.; Peng, Z.; Zhao, P.

2012-12-01

We present a 3D P-wave velocity (Vp) model of the Parkfield region that utilizes existing P-wave arrival time data, including fault zone head waves (FZHW), plus new data from direct wave secondary arrivals (DWSA). The first-arrival and DWSA travel times are obtained as the global and local minimum travel time paths, respectively. The inclusion of DWSA results in as much as a 10% increase in the across-fault velocity contrast for the Vp model at Parkfield relative to Thurber et al. (2006). Viewed along strike, three pronounced velocity contrast regions are observed: a pair of strong positive velocity contrasts (SW fast), one NW of the 1966 Parkfield hypocenter and the other SE of the 2004 Parkfield hypocenter, and a strong negative velocity contrast (NE fast) between the two hypocenters. The negative velocity contrast partially to entirely encompasses peak coseismic slip estimated in several slip models for the 2004 earthquake, suggesting that the negative velocity contrast played a part in defining the rupture patch of the 2004 Parkfield earthquake. We expand on this work by modifying our seismic tomography algorithm to incorporate arrival polarizations (azimuths). Synthetic tests will be presented to demonstrate the improvements in velocity structure when arrival polarizations are incorporated. These tests will compare the synthetic model recovered when FZHW/DWSA arrivals as well as existing P-wave arrival time data are inverted to that recovered with the same dataset with the inclusion of arrival polarizations. We plan to extend this work to carry out a full scale seismic tomography/relocation inversion at Parkfield, CA utilizing arrival polarizations from all first-P arrivals, and FZHW/DWSA arrivals as well as existing P-wave arrival time data. This effort requires the determination of polarization data for all P-waves and FZHW's at Parkfield. To this end, we use changes in the arrival azimuth from fault normal to source-receiver direction to identify FZHW and DWSA arrivals. We also use an eigenvalue decomposition to determine the direction of the incoming wave field, and to measure the arrival azimuths. This work is supported by the USGS Earthquake Hazards Program under grant numbers G11AP20027 and G11AP20028.

Analyzing the Possibility of Dynamic Earthquake Triggering in Socorro, New Mexico

NASA Astrophysics Data System (ADS)

Morton, E.; Bilek, S. L.

2011-12-01

The release of energy during an earthquake changes the stress state and seismicity both locally and remotely. Far-field stress changes can lead to triggered earthquakes coinciding with the arrival of the surface waves. This dynamic triggering is found to occur in a variety of tectonic settings, but in particular magmatic regions. Here we test whether the Socorro Magma Body region in central New Mexico hosts triggered seismicity. Preliminary inspection of continuous network data in central New Mexico suggested a local triggered event with the passage of surface waves from an MW 6.9 event in 2009. For a more comprehensive view, we examine data from 379 earthquakes MW ≥ 6.0 between January 15, 2008 to March 13, 2010 recorded on the EarthScope USArray Transportable Network stations located within New Mexico and providing more dense coverage for better detectability. Waveforms from twenty EarthScope stations were windowed around the time of the large event, high-pass filtered at 5 Hz to remove low frequency signals and analyzed to detect high frequency triggered local earthquakes. For each possible trigger detected, waveforms from nine short-period stations in the Socorro Seismic Network were added to aid in locating the events. In the time period analyzed, twelve triggered events were detected. Only one of these events, on August 30, 2009, corresponded to the arrival of surface waves, occurring about a minute after their arrival. The majority of the triggered events occur well after the arrival of the surface waves, indicating that they are either independent of the main shock or the result of delayed dynamic triggering. Delayed dynamic triggering can occur hours or days after the passage of surface waves, and are marked by an increase in seismicity relative to background. Only one of the events, on September 18, 2009, occurred within the Socorro Magma Body area. The rest of these events occur spread throughout New Mexico. The widely spread distribution of possibly triggered events and the low ratio of triggers to main shocks indicates that the rifted magmatic region above the Socorro Magma Body is not particularly susceptible to dynamic triggering from remote main shocks. The lack of direct correspondence to a seismic phase can mean that the detected events may be independent (not triggered events), or the result of delayed dynamic triggering. A comparison to randomly chosen waveforms within the time period as background will reveal if the possible events are a result of delayed dynamic triggering or part of the background.

Measuring ionospheric movements using totally reflected radio waves

NASA Astrophysics Data System (ADS)

Sadler, Elaine M.; Whitehead, J. D.; From, W. R.

1988-02-01

It is shown that for radio waves of a particular frequency reflected totally from the ionosphere the effect of refraction as well as reflection can be simulated by an effective reflecting surface. This mirror-like surface will give the correct angle of arrival and Doppler shift for all radars operating at this frequency. It is theoretically possible for the effective reflecting surface to be folded back on itself, but this is unlikely except for F-region echoes refracted by sporadic E-clouds. If the surface is not folded and exists everywhere, it is always possible to describe its motion and change in terms of wave undulations. Experimental data for F-region echoes show that these wave undulations are very dispersive. However, the matching between the best fitting model and the experimental data is worse than expected for reasons we do not understand.

Measuring ionospheric movements using totally reflected radio waves

NASA Astrophysics Data System (ADS)

From, W. R.; Sadler, Elaine M.; Whitehead, J. D.

1988-02-01

It is shown that for radio waves of a particular frequency reflected totally from the ionosphere the effect of refraction as well as reflection can be simulated by an effective reflecting surface. This mirrorlike surface will give the correct angle of arrival and Doppler shift for all radars operating at this frequency. It is theoretically possible for the effective reflecting surface to be folded back on itself, but this is unlikely except for F-region echoes refracted by sporadic E-clouds. If the surface is not folded and exists everywhere, it is always possible to describe its motion and change in terms of wave undulations. Experimental data for F-region echoes show that these wave undulations are very dispersive. However, the matching between the best fitting model and the experimental data is worse than expected for reasons we do not understand.

Head wave correlations in ambient noise.

PubMed

Gebbie, John; Siderius, Martin

2016-07-01

Ambient ocean noise is processed with a vertical line array to reveal coherent time-separated arrivals suggesting the presence of head wave multipath propagation. Head waves, which are critically propagating water waves created by seabed waves traveling parallel to the water-sediment interface, can propagate faster than water-only waves. Such eigenrays are much weaker than water-only eigenrays, and are often completely overshadowed by them. Surface-generated noise is different whereby it amplifies the coherence between head waves and critically propagating water-only waves, which is measured by cross-correlating critically steered beams. This phenomenon is demonstrated both experimentally and with a full wave simulation.

Amplification of seismic waves by the Seattle basin, Washington State

USGS Publications Warehouse

Pratt, T.L.; Brocher, T.M.; Weaver, C.S.; Creager, K.C.; Snelson, C.M.; Crosson, R.S.; Miller, K.C.; Trehu, A.M.

2003-01-01

Recordings of the 1999 Mw 7.6 Chi-Chi (Taiwan) earthquake, two local earthquakes, and five blasts show seismic-wave amplification over a large sedimentary basin in the U.S. Pacific Northwest. For weak ground motions from the Chi-Chi earthquake, the Seattle basin amplified 0.2- to 0.8-Hz waves by factors of 8 to 16 relative to bedrock sites west of the basin. The amplification and peak frequency change during the Chi-Chi coda: the initial S-wave arrivals (0-30 sec) had maximum amplifications of 12 at 0.5-0.8 Hz, whereas later arrivals (35-65 sec) reached amplifications of 16 at 0.3-0.5 Hz. Analysis of local events in the 1.0- to 10.0-Hz frequency range show fourfold amplifications for 1.0-Hz weak ground motion over the Seattle basin. Amplifications decrease as frequencies increase above 1.0 Hz, with frequencies above 7 Hz showing lower amplitudes over the basin than at bedrock sites. Modeling shows that resonance in low-impedance deposits forming the upper 550 m of the basin beneath our profile could cause most of the observed amplification, and the larger amplification at later arrival times suggests surface waves also play a substantial role. These results emphasize the importance of shallow deposits in determining ground motions over large basins.

Explosively Driven Shock Induced Damage in OFHC Copper

DOE Office of Scientific and Technical Information (OSTI.GOV)

Koller, D. D.; Hixson, R. S.; Gray, G. T. III

OFHC Cu samples were subjected to shock loading using plane wave HE lenses to produce a uniaxial Taylor wave profile (shock followed by immediate release). Upon arrival of the shock wave at the free surface of the sample, the wave is reflected and propagates back into the sample as a release wave. It is the interaction of initial and reflected release waves that place the material in a localized state of tension which can ultimately result in damage and possibly complete failure of the material. The peak tensile stress and its location in the material are determined by the wavemore » shape. Damage evolution processes and localized behavior are discussed based on results from time-resolved free surface velocity (VISAR) interferometry and post shock metallurgical analysis of the soft recovered samples.« less

A previously unreported type of seismic source in the firn layer of the East Antarctic Ice Sheet

NASA Astrophysics Data System (ADS)

Lough, Amanda C.; Barcheck, C. Grace; Wiens, Douglas A.; Nyblade, Andrew; Anandakrishnan, Sridhar

2015-11-01

We identify a unique type of seismic source in the uppermost part of the East Antarctic Ice Sheet recorded by temporary broadband seismic arrays in East Antarctica. These sources, termed "firnquakes," are characterized by dispersed surface wave trains with frequencies of 1-10 Hz detectable at distances up to 1000 km. Events show strong dispersed Rayleigh wave trains and an absence of observable body wave arrivals; most events also show weaker Love waves. Initial events were discovered by standard detection schemes; additional events were then detected with a correlation scanner using the initial arrivals as templates. We locate sources by determining the L2 misfit for a grid of potential source locations using Rayleigh wave arrival times and polarization directions. We then perform a multiple-filter analysis to calculate the Rayleigh wave group velocity dispersion and invert the group velocity for shear velocity structure. The resulting velocity structure is used as an input model to calculate synthetic seismograms. Inverting the dispersion curves yields ice velocity structures consistent with a low-velocity firn layer ~100 m thick and show that velocity structure is laterally variable. The absence of observable body wave phases and the relative amplitudes of Rayleigh waves and noise constrain the source depth to be less than 20 m. The presence of Love waves for most events suggests the source is not isotropic. We propose the events are linked to the formation of small crevasses in the firn, and several events correlate with shallow crevasse fields mapped in satellite imagery.

Automatic pickup of arrival time of channel wave based on multi-channel constraints

NASA Astrophysics Data System (ADS)

Wang, Bao-Li

2018-03-01

Accurately detecting the arrival time of a channel wave in a coal seam is very important for in-seam seismic data processing. The arrival time greatly affects the accuracy of the channel wave inversion and the computed tomography (CT) result. However, because the signal-to-noise ratio of in-seam seismic data is reduced by the long wavelength and strong frequency dispersion, accurately timing the arrival of channel waves is extremely difficult. For this purpose, we propose a method that automatically picks up the arrival time of channel waves based on multi-channel constraints. We first estimate the Jaccard similarity coefficient of two ray paths, then apply it as a weight coefficient for stacking the multichannel dispersion spectra. The reasonableness and effectiveness of the proposed method is verified in an actual data application. Most importantly, the method increases the degree of automation and the pickup precision of the channel-wave arrival time.

Calculating qP-wave traveltimes in 2-D TTI media by high-order fast sweeping methods with a numerical quartic equation solver

NASA Astrophysics Data System (ADS)

Han, Song; Zhang, Wei; Zhang, Jie

2017-09-01

A fast sweeping method (FSM) determines the first arrival traveltimes of seismic waves by sweeping the velocity model in different directions meanwhile applying a local solver. It is an efficient way to numerically solve Hamilton-Jacobi equations for traveltime calculations. In this study, we develop an improved FSM to calculate the first arrival traveltimes of quasi-P (qP) waves in 2-D tilted transversely isotropic (TTI) media. A local solver utilizes the coupled slowness surface of qP and quasi-SV (qSV) waves to form a quartic equation, and solve it numerically to obtain possible traveltimes of qP-wave. The proposed quartic solver utilizes Fermat's principle to limit the range of the possible solution, then uses the bisection procedure to efficiently determine the real roots. With causality enforced during sweepings, our FSM converges fast in a few iterations, and the exact number depending on the complexity of the velocity model. To improve the accuracy, we employ high-order finite difference schemes and derive the second-order formulae. There is no weak anisotropy assumption, and no approximation is made to the complex slowness surface of qP-wave. In comparison to the traveltimes calculated by a horizontal slowness shooting method, the validity and accuracy of our FSM is demonstrated.

On near-source earthquake triggering

USGS Publications Warehouse

Parsons, T.; Velasco, A.A.

2009-01-01

When one earthquake triggers others nearby, what connects them? Two processes are observed: static stress change from fault offset and dynamic stress changes from passing seismic waves. In the near-source region (r ??? 50 km for M ??? 5 sources) both processes may be operating, and since both mechanisms are expected to raise earthquake rates, it is difficult to isolate them. We thus compare explosions with earthquakes because only earthquakes cause significant static stress changes. We find that large explosions at the Nevada Test Site do not trigger earthquakes at rates comparable to similar magnitude earthquakes. Surface waves are associated with regional and long-range dynamic triggering, but we note that surface waves with low enough frequency to penetrate to depths where most aftershocks of the 1992 M = 5.7 Little Skull Mountain main shock occurred (???12 km) would not have developed significant amplitude within a 50-km radius. We therefore focus on the best candidate phases to cause local dynamic triggering, direct waves that pass through observed near-source aftershock clusters. We examine these phases, which arrived at the nearest (200-270 km) broadband station before the surface wave train and could thus be isolated for study. Direct comparison of spectral amplitudes of presurface wave arrivals shows that M ??? 5 explosions and earthquakes deliver the same peak dynamic stresses into the near-source crust. We conclude that a static stress change model can readily explain observed aftershock patterns, whereas it is difficult to attribute near-source triggering to a dynamic process because of the dearth of aftershocks near large explosions.

Three-dimensional waveform sensitivity kernels

NASA Astrophysics Data System (ADS)

Marquering, Henk; Nolet, Guust; Dahlen, F. A.

1998-03-01

The sensitivity of intermediate-period (~10-100s) seismic waveforms to the lateral heterogeneity of the Earth is computed using an efficient technique based upon surface-wave mode coupling. This formulation yields a general, fully fledged 3-D relationship between data and model without imposing smoothness constraints on the lateral heterogeneity. The calculations are based upon the Born approximation, which yields a linear relation between data and model. The linear relation ensures fast forward calculations and makes the formulation suitable for inversion schemes; however, higher-order effects such as wave-front healing are neglected. By including up to 20 surface-wave modes, we obtain Fréchet, or sensitivity, kernels for waveforms in the time frame that starts at the S arrival and which includes direct and surface-reflected body waves. These 3-D sensitivity kernels provide new insights into seismic-wave propagation, and suggest that there may be stringent limitations on the validity of ray-theoretical interpretations. Even recently developed 2-D formulations, which ignore structure out of the source-receiver plane, differ substantially from our 3-D treatment. We infer that smoothness constraints on heterogeneity, required to justify the use of ray techniques, are unlikely to hold in realistic earth models. This puts the use of ray-theoretical techniques into question for the interpretation of intermediate-period seismic data. The computed 3-D sensitivity kernels display a number of phenomena that are counter-intuitive from a ray-geometrical point of view: (1) body waves exhibit significant sensitivity to structure up to 500km away from the source-receiver minor arc; (2) significant near-surface sensitivity above the two turning points of the SS wave is observed; (3) the later part of the SS wave packet is most sensitive to structure away from the source-receiver path; (4) the sensitivity of the higher-frequency part of the fundamental surface-wave mode is wider than for its faster, lower-frequency part; (5) delayed body waves may considerably influence fundamental Rayleigh and Love waveforms. The strong sensitivity of waveforms to crustal structure due to fundamental-mode-to-body-wave scattering precludes the use of phase-velocity filters to model body-wave arrivals. Results from the 3-D formulation suggest that the use of 2-D and 1-D techniques for the interpretation of intermediate-period waveforms should seriously be reconsidered.

Modelling guided waves in the Alaskan-Aleutian subduction zone

NASA Astrophysics Data System (ADS)

Coulson, Sophie; Garth, Thomas; Reitbrock, Andreas

2016-04-01

Subduction zone guided wave arrivals from intermediate depth earthquakes (70-300 km depth) have a huge potential to tell us about the velocity structure of the subducting oceanic crust as it dehydrates at these depths. We see guided waves as the oceanic crust has a slower seismic velocity than the surrounding material, and so high frequency energy is retained and delayed in the crustal material. Lower frequency energy is not retained in this crustal waveguide and so travels at faster velocities of the surrounding material. This gives a unique observation at the surface with low frequency energy arriving before the higher frequencies. We constrain this guided wave dispersion by comparing the waveforms recorded in real subduction zones with simulated waveforms, produced using finite difference full waveform modelling techniques. This method has been used to show that hydrated minerals in the oceanic crust persist to much greater depths than accepted thermal petrological subduction zone models would suggest in Northern Japan (Garth & Rietbrock, 2014a), and South America (Garth & Rietbrock, in prep). These observations also suggest that the subducting oceanic mantle may be highly hydrated at intermediate depth by dipping normal faults (Garth & Rietbrock 2014b). We use this guided wave analysis technique to constrain the velocity structure of the down going ~45 Ma Pacific plate beneath Alaska. Dispersion analysis is primarily carried out on guided wave arrivals recorded on the Alaskan regional seismic network. Earthquake locations from global earthquake catalogues (ISC and PDE) and regional earthquake locations from the AEIC (Alaskan Earthquake Information Centre) catalogue are used to constrain the slab geometry and to identify potentially dispersive events. Dispersed arrivals are seen at stations close to the trench, with high frequency (>2 Hz) arrivals delayed by 2 - 4 seconds. This dispersion is analysed to constrain the velocity and width of the proposed waveguide. The velocity structure of this relatively young subducting plate is compared to the velocity structure resolved in the older oceanic lithosphere subducted beneath Northern Japan. We also use guided wave observations to investigate the thickness and low velocity structure of the subducting Yakutat terrain. Additionally we discuss the dependence of the inferred slab geometry on the earthquake catalogues that are used.

6C polarization analysis - seismic direction finding in coherent noise, automated event identification, and wavefield separation

NASA Astrophysics Data System (ADS)

Schmelzbach, C.; Sollberger, D.; Greenhalgh, S.; Van Renterghem, C.; Robertsson, J. O. A.

2017-12-01

Polarization analysis of standard three-component (3C) seismic data is an established tool to determine the propagation directions of seismic waves recorded by a single station. A major limitation of seismic direction finding methods using 3C recordings, however, is that a correct propagation-direction determination is only possible if the wave mode is known. Furthermore, 3C polarization analysis techniques break down in the presence of coherent noise (i.e., when more than one event is present in the analysis time window). Recent advances in sensor technology (e.g., fibre-optical, magnetohydrodynamic angular rate sensors, and ring laser gyroscopes) have made it possible to accurately measure all three components of rotational ground motion exhibited by seismic waves, in addition to the conventionally recorded three components of translational motion. Here, we present an extension of the theory of single station 3C polarization analysis to six-component (6C) recordings of collocated translational and rotational ground motions. We demonstrate that the information contained in rotation measurements can help to overcome some of the main limitations of standard 3C seismic direction finding, such as handling multiple arrivals simultaneously. We show that the 6C polarisation of elastic waves measured at the Earth's free surface does not only depend on the seismic wave type and propagation direction, but also on the local P- and S-wave velocities just beneath the recording station. Using an adaptation of the multiple signal classification algorithm (MUSIC), we demonstrate how seismic events can univocally be identified and characterized in terms of their wave type. Furthermore, we show how the local velocities can be inferred from single-station 6C data, in addition to the direction angles (inclination and azimuth) of seismic arrivals. A major benefit of our proposed 6C method is that it also allows the accurate recovery of the wave type, propagation directions, and phase velocities of multiple, interfering arrivals in one time window. We demonstrate how this property can be exploited to separate the wavefield into its elastic wave-modes and to isolate or suppress waves arriving from specific directions (directional filtering), both in a fully automated fashion.

P-wave fault-plane solutions and the generation of surface waves by earthquakes in the western United States

NASA Astrophysics Data System (ADS)

Patton, Howard J.

1985-08-01

Surface waves recorded at regional distances are used to study the source mechanisms of seven earthquakes in the western United States with magnitudes between 4.3 and 5.5. The source mechanisms of events in or on the margins of the Basin and Range show T-axis with an azimuth of N85°W +/- 16° and a plunge of 12° +/- 16°. Of the seven events, four have P-wave solutions that are inconsistent with surface-wave observations. Azimuths of the T-axis obtained from the surface-wave mechanisms and from the P-wave solutions differ by up to 45°. These events have dip-slip or oblique-slip mechanisms, and the source depths for three of the events are 5 km or less. Their source mechanisms and small magnitudes make identification of the P-wave first motion difficult due to poor signal-to-noise ratio of the initial P-wave and close arrivals of pP or sP with significant amplitude. We suggest that mis-identification of the P-wave first motion and distortion of the body-wave ray paths due to non-planar structure were sources of error in determining the nodal planes for these events.

Constraints for the subsurface structure at the Abydos site on 67P/Churyumov-Gerasimenko resulting from CASSE listening to the MUPUS insertion phase

NASA Astrophysics Data System (ADS)

Knapmeyer, M.; Fischer, H. H.; Joerg, K.; Seidensticker, K. J.

2016-12-01

During the more than 3 hours of the MUPUS PEN insertion phase at Abydos, the Comet Acoustic Surface Sounding Experiment (CASSE), part of SESAME, recorded hammer strokes of MUPUS with all three accelerometers, which are housed in Philae's feet. Stroke times recorded in the MUPUS housekeeping data, the identification of two adjacent strokes in a single recording, and laboratory experiments concerning the properties of sound transmission through Philae's legs and feet foster the identification of the recorded signals as hammer strokes, and as elastic waves transmitted through the comet. A hammer stroke exerted on a surface dominantly excites surface waves of the Rayleigh type. By comparison of arrival times at the individual feet, we estimate the propagation velocity of these Rayleigh waves to be at least 80 m/s. With the bulk density of 533 ± 6 kg/m3 as derived from tracking Rosetta (Pätzold et al., Nature, vol. 530, 2016), this velocity translates into a shear modulus of the comet material of at least 3.2 MPa. Shear modulus scales with velocity squared, so when taking into account the formal uncertainties arising from the arrival time inversion, the shear modulus may easily be as large as 10 MPa. This is still low compared to solid rock or monocrystalline ice, but is compatible with highly porous materials. The recorded signals are only weakly dispersive: Energy at frequencies below approx. 200 Hz arrives slightly later than at higher frequencies. This indicates the presence of a surface layer, to which the above propagation velocity and shear modulus apply, on top of a material with an even lower propagation velocity and shear modulus. The boundary between the two materials is likely more than 20 cm below the surface. We conclude that the results of CASSE listening to MUPUS support the hypothesis of surficial sintering at least for the Abydos site on 67P/Churyumov-Gerasimenko.

Optical fiber repeatered transmission systems utilizing SAW filters

NASA Astrophysics Data System (ADS)

Rosenberg, R. L.; Ross, D. G.; Trischitta, P. R.; Fishman, D. A.; Armitage, C. B.

1983-05-01

Baseband digital transmission-line systems capable of signaling rates of several hundred to several thousand Mbit/s are presently being developed around the world. The pulse regeneration process is gated by a timing wave which is synchronous with the symbol rate of the arriving pulse stream. Synchronization is achieved by extracting a timing wave from the arriving pulse stream, itself. To date, surface acoustic-wave (SAW) filters have been widely adopted for timing recovery in the in-line regenerators of high-bit-rate systems. The present investigation has the objective to acquaint the SAW community in general, and SAW filter suppliers in particular, with the requirements for timing recovery filters in repeatered digital transmission systems. Attention is given to the system structure, the timing loop function, the system requirements affecting the timing-recovery filter, the decision process, timing jitter accumulation, the filter 'ringing' requirement, and aspects of reliability.

Forward volume and surface magnetostatic modes in an yttrium iron garnet film for out-of-plane magnetic fields: Theory and experiment

NASA Astrophysics Data System (ADS)

Lim, Jinho; Bang, Wonbae; Trossman, Jonathan; Amanov, Dovran; Ketterson, John B.

2018-05-01

We present experimental and theoretical results on the propagation of magnetostatic spin waves in a film of yttrium iron garnet (YIG) for out-of-plane magnetic fields for which propagation in opposite directions is nonreciprocal in the presence of a metal layer. The plane studied is defined by the film normal n and n × k where k is the wave vector of the mode. Spin waves in this setting are classified as forward volume waves or surface waves and display non-reciprocity in the presence of an adjacent metal layer except for when H//n. The measurements are carried out in a transmission geometry, and a microwave mixer is used to measure the change of phase, and with it the evolution of wavevector, of the arriving spin wave with external magnetic field.

S-wave triggering of tremor beneath the Parkfield, California, section of the San Andreas fault by the 2011 Tohoku, Japan earthquake: observations and theory

USGS Publications Warehouse

Hill, David P.; Peng, Zhigang; Shelly, David R.; Aiken, Chastity

2013-01-01

The dynamic stresses that are associated with the energetic seismic waves generated by the Mw 9.0 Tohoku earthquake off the northeast coast of Japan triggered bursts of tectonic tremor beneath the Parkfield section of the San Andreas fault (SAF) at an epicentral distance of ∼8200  km. The onset of tremor begins midway through the ∼100‐s‐period S‐wave arrival, with a minor burst coinciding with the SHSH arrival, as recorded on the nearby broadband seismic station PKD. A more pronounced burst coincides with the Love arrival, followed by a series of impulsive tremor bursts apparently modulated by the 20‐ to 30‐s‐period Rayleigh wave. The triggered tremor was located at depths between 20 and 30 km beneath the surface trace of the fault, with the burst coincident with the S wave centered beneath the fault 30 km northwest of Parkfield. Most of the subsequent activity, including the tremor coincident with the SHSH arrival, was concentrated beneath a stretch of the fault extending from 10 to 40 km southeast of Parkfield. The seismic waves from the Tohoku epicenter form a horizontal incidence angle of ∼14°, with respect to the local strike of the SAF. Computed peak dynamic Coulomb stresses on the fault at tremor depths are in the 0.7–10 kPa range. The apparent modulation of tremor bursts by the small, strike‐parallel Rayleigh‐wave stresses (∼0.7  kPa) is likely enabled by pore pressure variations driven by the Rayleigh‐wave dilatational stress. These results are consistent with the strike‐parallel dynamic stresses (δτs) associated with the S, SHSH, and surface‐wave phases triggering small increments of dextral slip on the fault with a low friction (μ∼0.2). The vertical dynamic stresses δτd do not trigger tremor with vertical or oblique slip under this simple Coulomb failure model.

Nonlinear interaction of strong S-waves with the rupture front in the shallow subsurface

NASA Astrophysics Data System (ADS)

Sleep, N. H.

2017-12-01

Shallow deformation in moderate to large earthquakes is sometimes distributed rather than being concentrated on a single fault plane. Strong high-frequency S-waves interact with the rupture front to produce this effect. For strike-slip faults, the rupture propagation velocity is a fraction of the S-wave velocity. The rupture propagation vector refracts essentially vertically in the low (S-wave) velocity shallow subsurface. So does the propagation direction of S-waves. The shallow rupture front is essentially mode 3 near the surface. Strong S-waves arrive before the rupture front. They continue to arrive for several seconds in a large event. There are simple scaling relationships. The dynamic Coulomb stress ratio of horizontal stress on horizontal planes from S-waves is the normalized acceleration in g's. For fractured rock and gravel, frictional failure occurs when the normalized acceleration exceeds the effective coefficient of friction. Acceleration tends to saturate at that level as the anelastic strain rate increases rapidly with stress. For muddy materials, failure begins at a low normalized acceleration but increases slowly with dynamic stress. Dynamic accelerations sometimes exceed 1 g. In both cases, the rupture tip finds the shallow subsurface already in nonlinear failure down to a few to tens of meters depth. The material does not distinguish between S-wave and rupture tip stresses. Both stresses add to the stress invariant and hence to the anelastic strain rate tensor. Surface anelastic strain from fault slip is thus distributed laterally over a distance scaling to the depth of nonlinearity from S-waves. The environs of the fault anelastically accommodate the fault slip at depth. This process differs from blind faults where the shallow coseismic strain is mostly elastic and interseismic anelastic processes accommodate the long-term shallow deformation.

Observations of basin ground motions from a dense seismic array in San Jose, California

USGS Publications Warehouse

Frankel, A.; Carver, D.; Cranswick, E.; Bice, T.; Sell, R.; Hanson, S.

2001-01-01

We installed a dense array of 41 digital seismographs in San Jose, California, to evaluate in detail the effects of a deep sedimentary basin and shallow sedimentary deposits on earthquake ground motions. This urban array is located near the eastern edge of the Santa Clara Valley and spans the Evergreen sedimentary basin identified by gravity data. Average station spacing is 1 km, with three stations initially spaced 110 m apart. Despite the high-noise urban environment, the stations of the array successfully triggered on and recorded small local earthquakes (M 2.5-2.8 at 10-25 km distance) and larger regional events such as the M 5.0 Bolinas earthquake (90 km distance), M 4.6-5.6 earthquakes near Mammoth Lakes (270 km distance), M 4.9-5.6 events in western Nevada (420 km distance) and the M 7.1 Hector Mine earthquake (590 km distance). Maps of spectral ratios across the array show that the highest amplitudes in all frequency bands studied (0.125-8 Hz) are generally observed at stations farther from the eastern edge of the Santa Clara Valley. Larger spectral amplitudes are often observed above the western edge of the Evergreen Basin. Snapshots of the recorded wavefield crossing the array for regional events to the east reveal that large, low-frequency (0.125-0.5 Hz) arrivals after the S-wave travel from south to north across the array. A moving-window, cross-correlation analysis finds that these later arrivals are surface waves traveling from the south. The timing and propagation direction of these arrivals indicates that they were likely produced by scattering of incident S waves at the border of the Santa Clara Valley to the south of the array. It is remarkable that the largest low-frequency phases at many of the valley sites for regional events to the east are basin surface waves coming from a direction about 70 degrees different from that of the epicenters. Basin surface waves emanating from the eastern edge of the valley are also identified by the cross-correlation analysis.

HF Surface Wave Radar Tests at the Eastern China Sea

NASA Astrophysics Data System (ADS)

Wu, Xiong Bin; Cheng, Feng; Wu, Shi Cai; Yang, Zi Jie; Wen, Biyang; Shi, Zhen Hua; Tian, Jiansheng; Ke, Hengyu; Gao, Huotao

2005-01-01

The HF surface wave radar system OSMAR2000 adopts Frequency Modulated Interrupted Continuous Waveform (FMICW) and its 120m-antenna array is transmitting/receiving co-used. MUSIC and MVM are applied to obtain sea echo's direction of arrival (DOA) when extracting currents information. Verification tests of OSMAR2000 ocean surface dynamics detection against in-situ measurements had been accomplished on Oct. 23~29, 2000. Ship detection test was carried out on Dec.24, 2001. It shows that OSMAR2000 is capable of detecting 1000 tons ships with a wide beam out to 70 km. This paper introduces the radar system and the applied DOA estimation methods in the first, and then presents ship detection results and some sea state measurement results of surface currents and waves. The results indicate the validity of the developed radar system and the effectiveness of the applied signal processing methods.

Effects of upper ocean sound-speed structure on deep acoustic shadow-zone arrivals at 500- and 1000-km range.

PubMed

Van Uffelen, Lora J; Worcester, Peter F; Dzieciuch, Matthew A; Rudnick, Daniel L; Colosi, John A

2010-04-01

Deep acoustic shadow-zone arrivals observed in the late 1990s in the North Pacific Ocean reveal significant acoustic energy penetrating the geometric shadow. Comparisons of acoustic data obtained from vertical line arrays deployed in conjunction with 250-Hz acoustic sources at ranges of 500 and 1000 km from June to November 2004 in the North Pacific, with simulations incorporating scattering consistent with the Garrett-Munk internal-wave spectrum, are able to describe both the energy contained in and vertical extent of deep shadow-zone arrivals. Incoherent monthly averages of acoustic timefronts indicate that lower cusps associated with acoustic rays with shallow upper turning points (UTPs), where sound-speed structure is most variable and seasonally dependent, deepen from June to October as the summer thermocline develops. Surface-reflected rays, or those with near-surface UTPs, exhibit less scattering due to internal waves than in later months when the UTP deepens. Data collected in November exhibit dramatically more vertical extension than previous months. The depth to which timefronts extend is a complex combination of deterministic changes in the depths of the lower cusps as the range-average profiles evolve with seasonal change and of the amount of scattering, which depends on the mean vertical gradients at the depths of the UTPs.

Teleseismic and regional data analysis for estimating depth, mechanism and rupture process of the 3 April 2017 MW 6.5 Botswana earthquake and its aftershock (5 April 2017, MW 4.5)

NASA Astrophysics Data System (ADS)

Letort, J.; Guilhem Trilla, A.; Ford, S. R.; Sèbe, O.; Causse, M.; Cotton, F.; Campillo, M.; Letort, G.

2017-12-01

We constrain the source, depth, and rupture process of the Botswana earthquake of April 3, 2017, as well as its largest aftershock (5 April 2017, Mw 4.5). This earthquake is the largest recorded event (Mw 6.5) in the East African rift system since 1970, making one important case study to better understand source processes in stable continental regions. For the two events an automatic cepstrum analysis (Letort et al., 2015) is first applied on respectively 215 and 219 teleseismic records, in order to detect depth phase arrivals (pP, sP) in the P-coda. Coherent detections of depth phases for different azimuths allow us to estimate the hypocentral depths respectively at 28 and 23 km, suggesting that the events are located in the lower crust. A same cepstrum analysis is conducted on five other earthquakes with mb>4 in this area (from 2002 to 2017), and confirms a deep crustal seismicity cluster (around 20-30 km). The source mechanisms are then characterized using a joint inversion method by fitting both regional long-period surface-waves and teleseismic high-frequency body-waves. Combining regional and teleseismic data (as well as systematic comparisons between theoretical and observed regional surface-waves dispersion curves prior to the inversion) allows us to decrease epistemic uncertainties due to lack of regional data and poor knowledge about the local velocity structure. Focal mechanisms are both constrained as normal faulting with a northwest trending, and hypocentral depths are confirmed at 28 and 24 km. Finally, in order to study the mainshock rupture process, we originally apply a kymograph analysis method (an image processing method, commonly used in the field of cell biology for identifying motions of molecular motors, e.g. Mangeol et al., 2016). Here, the kymograph allows us to better identify high-frequency teleseismic P-arrivals inside the P-coda by tracking both reflected depth phase and direct P-wave arrivals radiated from secondary sources during the faulting process. Secondary P-arrivals are thus identified with a significant azimuthal variation of their arrival times (until 4s), allowing the localization of the source that generated these secondary waves. This analysis shows that the mainshock is probably a mix of at least two events, the second being 20-30 km further northwest along the fault.

Seismic anisotropy in the upper mantle beneath the MAGIC array, mid-Atlantic Appalachians: Constraints from SKS splitting and quasi-Love wave propagation

NASA Astrophysics Data System (ADS)

Aragon, J. C.; Long, M. D.; Benoit, M. H.; Servali, A.

2016-12-01

North America's eastern passive continental margin has been modified by several cycles of supercontinent assembly. Its complex surface geology and distinct topography provide evidence of these events, while also raising questions about the extent of deformation in the continental crust, lithosphere, and mantle during past episodes of rifting and mountain building. The Mid-Atlantic Geophysical Integrative Collaboration (MAGIC) is an EarthScope and GeoPRISMS-funded project that involves a collaborative effort among seismologists, geodynamicists, and geomorphologists. One component of the project is a broadband seismic array consisting of 28 instruments in a linear path from coastal Virginia to western Ohio, which operated between October 2013 and October 2016. A key science question addressed by the MAGIC project is the geometry of past lithospheric deformation and present-day mantle flow beneath the Appalachians, which can be probed using observations of seismic anisotropy Here we present observations of SKS splitting and quasi-Love wave arrivals from stations of the MAGIC array, which together constrain seismic anisotropy in the upper mantle. SKS splitting along the array reveals distinct regions of upper mantle anisotropy, with stations in and to the west of the range exhibiting fast directions parallel to the strike of the mountains. In contrast, weak splitting and null SKS arrivals dominate eastern stations in the coastal plain. Documented Love-to-Rayleigh wave scattering for surface waves originating the magnitude 8.3 Illapel, Chile earthquakes in September 2015 provides complementary constraints on anisotropy. These quasi-Love wave arrivals suggest a pronounced change in upper mantle anisotropy at the eastern edge of present-day Appalachian topography. Together, these observations increase our understanding of the extent of lithospheric deformation beneath North America associated with Appalachian orogenesis, as well as the pattern of present-day mantle flow beneath the passive margin.

Passive characterization of hydrofracture properties using signals from the hydraulic pumps

DOE Office of Scientific and Technical Information (OSTI.GOV)

Rector, J.W. III; Dong, Qichen

1995-12-31

In this study we utilize conical shear wave arrivals recorded in geophone observation wells to characterize a hydrofracture performed in the South Belridge Diatomite oil field. The conical wave arrivals are initially created by the hydraulic pumps on the surface, which send tube waves down the treatment borehole. Since the tube wave velocity in the Diatomite is greater than the shear formation velocity (the shear velocity in the diatomite is about 2,200 ft/s) cortical shear waves are radiated into the formation by the tube waves traveling down the treatment borehole. We use the decrease in amplitude of the tube wavemore » as it passes through the fracture zone to image changes in hydraulic conductivity of the fracture. By combining this information with estimates of the fracture height we obtain estimates of fracture width changes over time using the model of Tang and Cheng (1993). We find an excellent qualitative agreement between tube wave attenuation and pump pressure over time. Fracture widths estimated from the Tang and Cheng model appear to be consistent with the volume of injected fluid and the known length of the hydrofracture. Provided a monitor well can be instrumented, this technique holds potential for obtaining a relatively inexpensive real-time characterization of hydrofracs.« less

Auto Correlation Analysis of Coda Waves from Local Earthquakes for Detecting Temporal Changes in Shallow Subsurface Structures: the 2011 Tohoku-Oki, Japan Earthquake

NASA Astrophysics Data System (ADS)

Nakahara, Hisashi

2015-02-01

For monitoring temporal changes in subsurface structures I propose to use auto correlation functions of coda waves from local earthquakes recorded at surface receivers, which probably contain more body waves than surface waves. Use of coda waves requires earthquakes resulting in decreased time resolution for monitoring. Nonetheless, it may be possible to monitor subsurface structures in sufficient time resolutions in regions with high seismicity. In studying the 2011 Tohoku-Oki, Japan earthquake (Mw 9.0), for which velocity changes have been previously reported, I try to validate the method. KiK-net stations in northern Honshu are used in this analysis. For each moderate earthquake normalized auto correlation functions of surface records are stacked with respect to time windows in the S-wave coda. Aligning the stacked, normalized auto correlation functions with time, I search for changes in phases arrival times. The phases at lag times of <1 s are studied because changes at shallow depths are focused. Temporal variations in the arrival times are measured at the stations based on the stretching method. Clear phase delays are found to be associated with the mainshock and to gradually recover with time. The amounts of the phase delays are 10 % on average with the maximum of about 50 % at some stations. The deconvolution analysis using surface and subsurface records at the same stations is conducted for validation. The results show the phase delays from the deconvolution analysis are slightly smaller than those from the auto correlation analysis, which implies that the phases on the auto correlations are caused by larger velocity changes at shallower depths. The auto correlation analysis seems to have an accuracy of about several percent, which is much larger than methods using earthquake doublets and borehole array data. So this analysis might be applicable in detecting larger changes. In spite of these disadvantages, this analysis is still attractive because it can be applied to many records on the surface in regions where no boreholes are available.

Shear wave arrival time estimates correlate with local speckle pattern.

PubMed

Mcaleavey, Stephen A; Osapoetra, Laurentius O; Langdon, Jonathan

2015-12-01

We present simulation and phantom studies demonstrating a strong correlation between errors in shear wave arrival time estimates and the lateral position of the local speckle pattern in targets with fully developed speckle. We hypothesize that the observed arrival time variations are largely due to the underlying speckle pattern, and call the effect speckle bias. Arrival time estimation is a key step in quantitative shear wave elastography, performed by tracking tissue motion via cross-correlation of RF ultrasound echoes or similar methods. Variations in scatterer strength and interference of echoes from scatterers within the tracking beam result in an echo that does not necessarily describe the average motion within the beam, but one favoring areas of constructive interference and strong scattering. A swept-receive image, formed by fixing the transmit beam and sweeping the receive aperture over the region of interest, is used to estimate the local speckle pattern. Metrics for the lateral position of the speckle are found to correlate strongly (r > 0.7) with the estimated shear wave arrival times both in simulations and in phantoms. Lateral weighting of the swept-receive pattern improved the correlation between arrival time estimates and speckle position. The simulations indicate that high RF echo correlation does not equate to an accurate shear wave arrival time estimate-a high correlation coefficient indicates that motion is being tracked with high precision, but the location tracked is uncertain within the tracking beam width. The presence of a strong on-axis speckle is seen to imply high RF correlation and low bias. The converse does not appear to be true-highly correlated RF echoes can still produce biased arrival time estimates. The shear wave arrival time bias is relatively stable with variations in shear wave amplitude and sign (-20 μm to 20 μm simulated) compared with the variation with different speckle realizations obtained along a given tracking vector. We show that the arrival time bias is weakly dependent on shear wave amplitude compared with the variation with axial position/ local speckle pattern. Apertures of f/3 to f/8 on transmit and f/2 and f/4 on receive were simulated. Arrival time error and correlation with speckle pattern are most strongly determined by the receive aperture.

Shear Wave Arrival Time Estimates Correlate with Local Speckle Pattern

PubMed Central

McAleavey, Stephen A.; Osapoetra, Laurentius O.; Langdon, Jonathan

2016-01-01

We present simulation and phantom studies demonstrating a strong correlation between errors in shear wave arrival time estimates and the lateral position of the local speckle pattern in targets with fully developed speckle. We hypothesize that the observed arrival time variations are largely due to the underlying speckle pattern, and call the effect speckle bias. Arrival time estimation is a key step in quantitative shear wave elastography, performed by tracking tissue motion via cross correlation of RF ultrasound echoes or similar methods. Variations in scatterer strength and interference of echoes from scatterers within the tracking beam result in an echo that does not necessarily describe the average motion within the beam, but one favoring areas of constructive interference and strong scattering. A swept-receive image, formed by fixing the transmit beam and sweeping the receive aperture over the region of interest, is used to estimate the local speckle pattern. Metrics for the lateral position of the speckle are found to correlate strongly (r>0.7) with the estimated shear wave arrival times both in simulations and in phantoms. Lateral weighting of the swept-receive pattern improved the correlation between arrival time estimates and speckle position. The simulations indicate that high RF echo correlation does not equate to an accurate shear wave arrival time estimate – a high correlation coefficient indicates that motion is being tracked with high precision, but the location tracked is uncertain within the tracking beam width. The presence of a strong on-axis speckle is seen to imply high RF correlation and low bias. The converse does not appear to be true – highly correlated RF echoes can still produce biased arrival time estimates. The shear wave arrival time bias is relatively stable with variations in shear wave amplitude and sign (−20 μm to 20 μm simulated) compared to the variation with different speckle realizations obtained along a given tracking vector. We show that the arrival time bias is weakly dependent on shear wave amplitude compared to the variation with axial position/local speckle pattern. Apertures of f/3 to f/8 on transmit and f/2 and f/4 on receive were simulated. Arrival time error and correlation with speckle pattern are most strongly determined by the receive aperture. PMID:26670847

Full wave field recording of the vertical strain at SAFOD from local, regional and teleseismic earthquakes

NASA Astrophysics Data System (ADS)

Ellsworth, W. L.; Karrenbach, M. H.; Zumberge, M. A.

2017-12-01

The main borehole at the San Andreas Fault Observatory at Depth (SAFOD) contains optical fibers cemented in place in between casing strings from the surface to just below the top of the basement. The fibers are under tension of approximately 1 N and are housed in a 0.9 mm diameter stainless steel tube. Earth strain is transmitted to the fiber by frictional contact with the tube wall. One fiber has been in use as a vertical strainmeter since 2005, measuring the total strain between 9 and 740 m by laser interferometry. In June 2017 we attached an OptaSense Distributed Acoustic Sensing (DAS) system, model ODH3.1, to a second fiber that terminates at 864 m depth. The DAS laser interrogator measures the strain over a gauge length with a set spacing between gauge intervals. For this experiment we set the gauge length to 10 m with 1 m spacing between gauges. Including the surface run of the fiber, this gives us 936 channels measuring the vertical strain at a sample interval of 0.4 msec (2500 samples/s). Continuous recording of the string produces approximately 1 TB/day. During one month of data collection, we recorded local, regional and teleseismic earthquakes. With this recording geometry, the DAS system captures the full vertical wavefield between the basement interface and free surface, revealing direct, converted and refracted waves. Both P- and S- strain waves are clearly visible in the data, even for 10 km deep earthquakes located almost directly below the well (see figure). The incident and surface reflected wavefields can be separated by frequency-wavenumber filtering due to the large-aperture and fine spatial and temporal sampling. Up- and downgoing strain waves illuminate the subsurface within the sensor array's depth range. Accurate arrival time determinations of the initial arrival phase are possible due to consistent wave forms recorded at 1 m spatial intervals that can be used for fine-scale shallow velocity model estimation.

Regional Wave Propagation in Southeastern United States

NASA Astrophysics Data System (ADS)

Jemberie, A. L.; Langston, C. A.

2003-12-01

Broad band seismograms from the April 29, 2003, M4.6 Fort Payne, Alabama earthquake are analyzed to infer mechanisms of crustal wave propagation, crust and upper mantle velocity structure in southeastern United States, and source parameters of the event. In particular, we are interested in producing deterministic models of the distance attenuation of earthquake ground motions through computation of synthetic seismograms. The method first requires constraining the source parameters of an earthquake and then modeling the amplitude and times of broadband arrivals within the waveforms to infer appropriate layered earth models. A first look at seismograms recorded by stations outside the Mississippi Embayment (ME) show clear body phases such P, sP, Pnl, Sn and Lg. The ME signals are qualitatively different from others because they have longer durations and large surface waves. A straightforward interpretation of P wave arrival times shows a typical upper mantle velocity of 8.18 km/s. However, there is evidence of significantly higher P phase velocities at epicentral distances between 400 and 600km, that may be caused by a high velocity upper mantle anomaly; triplication of P-waves is seen in these seismograms. The arrival time differences between regional P and the depth phase sP at different stations are used to constrain the depth of the earthquake. The source depth lies between 9.5 km and 13km which is somewhat more shallow than the network location that was constrained to 15km depth. The Fort Payne earthquake is the largest earthquake to have occurred within the Eastern Tennessee Seismic Zone.

Cancellation of spurious arrivals in Green's function extraction and the generalized optical theorem

USGS Publications Warehouse

Snieder, R.; Van Wijk, K.; Haney, M.; Calvert, R.

2008-01-01

The extraction of the Green's function by cross correlation of waves recorded at two receivers nowadays finds much application. We show that for an arbitrary small scatterer, the cross terms of scattered waves give an unphysical wave with an arrival time that is independent of the source position. This constitutes an apparent inconsistency because theory predicts that such spurious arrivals do not arise, after integration over a complete source aperture. This puzzling inconsistency can be resolved for an arbitrary scatterer by integrating the contribution of all sources in the stationary phase approximation to show that the stationary phase contributions to the source integral cancel the spurious arrival by virtue of the generalized optical theorem. This work constitutes an alternative derivation of this theorem. When the source aperture is incomplete, the spurious arrival is not canceled and could be misinterpreted to be part of the Green's function. We give an example of how spurious arrivals provide information about the medium complementary to that given by the direct and scattered waves; the spurious waves can thus potentially be used to better constrain the medium. ?? 2008 The American Physical Society.

Towards Simulating a Realistic Planetary Seismic Wavefield: The Contribution of the Megaregolith and Low-Velocity Waveguides

NASA Technical Reports Server (NTRS)

Schmerr, Nicholas C.; Weber, Renee C.; Lin, Pei-Ying Patty; Thorne, Michael Scott; Garnero, Ed J.

2011-01-01

Lunar seismograms are distinctly different from their terrestrial counterparts. The Apollo lunar seismometers recorded moonquakes without distinct P- or S-wave arrivals; instead waves arrive as a diffuse coda that decays over several hours making the identification of body waves difficult. The unusual character of the lunar seismic wavefield is generally tied to properties of the megaregolith: it consists of highly fractured and broken crustal rock, the result of extensive bombardment of the Moon. The megaregolith extends several kilometers into the lunar crust, possibly into the mantle in some regions, and is covered by a thin coating of fine-scale dust. These materials possess very low seismic velocities that strongly scatter the seismic wavefield at high frequencies. Directly modeling the effects of the megaregolith to simulate an accurate lunar seismic wavefield is a challenging computational problem, owing to the inherent 3-D nature of the problem and the high frequencies (greater than 1 Hz) required. Here we focus on modeling the long duration code, studying the effects of the low velocities found in the megaregolith. We produce synthetic seismograms using 1-D slowness integration methodologies, GEMINI and reflectivity, and a 3-D Cartesian finite difference code, Wave Propagation Program, to study the effect of thin layers of low velocity on the surface of a planet. These codes allow us generate seismograms with dominant frequencies of approximately 1 Hz. For background lunar seismic structure we explore several models, including the recent model of Weber et al., Science, 2011. We also investigate variations in megaregolithic thickness, velocity, attenuation, and seismogram frequency content. Our results are compared to the Apollo seismic dataset, using both a cross correlation technique and integrated envelope approach to investigate coda decay. We find our new high frequency results strongly support the hypothesis that the long duration of the lunar seismic codes is generated by the presence of the low velocity megaregolith, and that the diffuse arrivals are a combination of scattered energy and multiple reverberations within this layer. The 3-D modeling indicates the extreme surface topography of the Moon adds only a small contribution to scattering effects, though local geology may play a larger role. We also study the effects of the megaregolith on core reflected and converted phases and other body waves. Our analysis indicates detection of core interacting arrivals with a polarization filter technique is robust and lends the possibility of detecting other body waves from the Moon.

Three-dimensional numerical modeling of bottom-diffracted surface-reflected arrivals in the North Pacific

NASA Astrophysics Data System (ADS)

Stephen, R. A.; Udovydchenkov, I. A.; Bolmer, T.; Komatitsch, D.; Tromp, J.; Casarotti, E.; Xie, Z.; Worcester, P. F.

2015-12-01

Bottom-diffracted surface-reflected (BDSR) arrivals were first identified in the 2004 Long-range Ocean Acoustic Propagation Experiment (Stephen et al, 2013, JASA, v.134, p.3307-3317). The BDSR mechanism provides a means for acoustic signals and noise from distant sources to appear with significant strength on the deep seafloor. At depths greater than the conjugate depth ambient noise and PE- predicted arrivals are sufficiently quiet that BDSR paths, scattered from small seamounts, can be the largest amplitude arrivals observed. The Ocean Bottom Seismometer Augmentation in the North Pacific (OBSANP) Experiment in June-July 2013 was designed to further define the characteristics of the BDSRs and to understand the conditions under which BDSRs are excited and propagate. An example of BDSR arrivals is shown on the record section in the figure. Other arrivals are the direct water wave and first and second multiples. The reciprocal of the BDSR mechanism also plays a role in T-phase excitation. To further understand the BDSR mechanism, the SPECFEM3D code was extended to handle high-frequency, deep water bottom scattering problems with actual bathymetry and a typical sound speed profile in the water column. The model size is 38km x 27km x 6.5km. The source is centered at 10Hz with a 5Hz bandwidth. Work supported by NSF and ONR.

On vertical seismic profile processing

DOE Office of Scientific and Technical Information (OSTI.GOV)

Tariel, P.; Michon, D.

1984-10-01

From the wealth of information which can be deduced from VSP, the information most directly comparable to well logs is considered: P-wave and S-wave interval velocity, acoustic impedance, and the velocity ratio ..gamma.. = V /SUB s/ /V /SUB p/ . This information not only allows better interpretation of surface seismic sections but also improves processing. For these results to be usable a number of precautions must be taken during acquisition and processing; the sampling in depth should be chosen in such a way that aliasing phenomena do not unnecessarily limit the spectra during the separation of upwards and downwardsmore » travelling waves. True amplitudes should be respected and checked by recording of signatures, and the interference of upwards and downwards travelling waves should be taken into account for the picking of first arrivals. The different steps in processing and the combination of results in the interpretation of surface seismic results are described with actual records.« less

First seismic shear wave velocity profile of the lunar crust as extracted from the Apollo 17 active seismic data by wavefield gradient analysis

NASA Astrophysics Data System (ADS)

Sollberger, David; Schmelzbach, Cedric; Robertsson, Johan O. A.; Greenhalgh, Stewart A.; Nakamura, Yosio; Khan, Amir

2016-04-01

We present a new seismic velocity model of the shallow lunar crust, including, for the first time, shear wave velocity information. So far, the shear wave velocity structure of the lunar near-surface was effectively unconstrained due to the complexity of lunar seismograms. Intense scattering and low attenuation in the lunar crust lead to characteristic long-duration reverberations on the seismograms. The reverberations obscure later arriving shear waves and mode conversions, rendering them impossible to identify and analyze. Additionally, only vertical component data were recorded during the Apollo active seismic experiments, which further compromises the identification of shear waves. We applied a novel processing and analysis technique to the data of the Apollo 17 lunar seismic profiling experiment (LSPE), which involved recording seismic energy generated by several explosive packages on a small areal array of four vertical component geophones. Our approach is based on the analysis of the spatial gradients of the seismic wavefield and yields key parameters such as apparent phase velocity and rotational ground motion as a function of time (depth), which cannot be obtained through conventional seismic data analysis. These new observables significantly enhance the data for interpretation of the recorded seismic wavefield and allow, for example, for the identification of S wave arrivals based on their lower apparent phase velocities and distinct higher amount of generated rotational motion relative to compressional (P-) waves. Using our methodology, we successfully identified pure-mode and mode-converted refracted shear wave arrivals in the complex LSPE data and derived a P- and S-wave velocity model of the shallow lunar crust at the Apollo 17 landing site. The extracted elastic-parameter model supports the current understanding of the lunar near-surface structure, suggesting a thin layer of low-velocity lunar regolith overlying a heavily fractured crust of basaltic material showing high (>0.4 down to 60 m) Poisson's ratios. Our new model can be used in future studies to better constrain the deep interior of the Moon. Given the rich information derived from the minimalistic recording configuration, our results demonstrate that wavefield gradient analysis should be critically considered for future space missions that aim to explore the interior structure of extraterrestrial objects by seismic methods. Additionally, we anticipate that the proposed shear wave identification methodology can also be applied to the routinely recorded vertical component data from land seismic exploration on Earth.

Role of Compressibility on Tsunami Propagation

NASA Astrophysics Data System (ADS)

Abdolali, Ali; Kirby, James T.

2017-12-01

In the present paper, we aim to reduce the discrepancies between tsunami arrival times evaluated from tsunami models and real measurements considering the role of ocean compressibility. We perform qualitative studies to reveal the phase speed reduction rate via a modified version of the Mild Slope Equation for Weakly Compressible fluid (MSEWC) proposed by Sammarco et al. (2013). The model is validated against a 3-D computational model. Physical properties of surface gravity waves are studied and compared with those for waves evaluated from an incompressible flow solver over realistic geometry for 2011 Tohoku-oki event, revealing reduction in phase speed.