64 Matching Results

Search Results

Advanced search parameters have been applied.

Seismic Data for Evaluation of Ground Motion Hazards in Las Vegas in Support of Test Site Readiness Ground Motion

Description: In this report we describe the data sets used to evaluate ground motion hazards in Las Vegas from nuclear tests at the Nevada Test Site. This analysis is presented in Rodgers et al. (2005, 2006) and includes 13 nuclear explosions recorded at the John Blume and Associates network, the Little Skull Mountain earthquake and a temporary deployment of broadband station in Las Vegas. The data are available in SAC format on CD-ROM as an appendix to this report.
Date: January 16, 2008
Creator: Rodgers, A
Partner: UNT Libraries Government Documents Department


Description: This project seeks to compute ground motions for large (M>6.5) scenario earthquakes on the Hayward Fault using realistic pseudodynamic ruptures, the USGS three-dimensional (3D) velocity model and anelastic finite difference simulations on parallel computers. We will attempt to bound ground motions by performing simulations with suites of stochastic rupture models for a given scenario on a given fault segment. The outcome of this effort will provide the average, spread and range of ground motions that can be expected from likely large earthquake scenarios. The resulting ground motions will be based on first-principles calculations and include the effects of slip heterogeneity, fault geometry and directivity, however, they will be band-limited to relatively low-frequency (< 1 Hz).
Date: January 9, 2008
Creator: Rodgers, A & Xie, X
Partner: UNT Libraries Government Documents Department

Structure of the Lithosphere and Upper Mantle Across the Arabian Peninsula

Description: Analysis of modern broadband (BB) waveform data allows for the inference of seismic velocity structure of the crust and upper mantle using a variety of techniques. This presentation will report inferences of seismic structure of the Arabian Plate using BB data from various networks. Most data were recorded by the Saudi Arabian National Digital Seismic Network (SANDSN) which consists of 38 (26 BB, 11 SP) stations, mostly located on the Arabian Shield. Additional data were taken from the 1995-7 Saudi Arabian IRIS-PASSCAL Deployment (9 BB stations) and other stations across the Peninsula. Crustal structure, inferred from teleseismic P-wave receiver functions, reveals thicker crust in the Arabian Platform (40-45 km) and the interior of the Arabian Shield (35-40 km) and thinner crust along the Red Sea coast. Lithospheric thickness inferred from teleseismic S-wave receiver functions reveals very thin lithosphere (40-80 km) along the Red Sea coast which thickens rapidly toward the interior of the Arabian Shield (100-120 km). We also observe a step of 20-40 km in lithospheric thickness across the Shield-Platform boundary. Seismic velocity structure of the upper mantle inferred from teleseismic P- and S-wave travel time tomography reveals large differences between the Shield and Platform, with the Shield being underlain by slower velocities, {+-}3% for P-waves and {+-}6% for S-waves. Seismic anisotropy was inferred from shear-wave splitting, using teleseismic SKS waveforms. Results reveal a splitting time of approximately 1.4 seconds, with the fast axis slightly east of north. The shear-wave splitting results are consistent across the Peninsula, with a slight clockwise rotation parallel for stations near the Gulf of Aqaba. In summary, these results allow us to make several conclusions about the tectonic evolution and current state of the Arabian Plate. Lithospheric thickness implies that thinning near the Red Sea has accompanied the rupturing of the Arabian-Nubian continental lithosphere. ...
Date: January 5, 2007
Creator: Al-Amri, A & Rodgers, A
Partner: UNT Libraries Government Documents Department

Simulations of the 1906 San Francisco Earthquake

Description: Simulations of the Great 1906 San Francisco earthquake are being performed as part of the event's centenary. LLNL is participating in a large effort to study this event and possible consequences if the event were to happen today. This document is meant to describe our efforts to others participating in the project.
Date: December 16, 2005
Creator: Rodgers, A; Petersson, A & Tkalcic, H
Partner: UNT Libraries Government Documents Department

A Model-Based Signal Processing Approach to Nuclear Explosion Monitoring

Description: This report describes research performed under Laboratory Research and Development Project 05-ERD-019, entitled ''A New Capability for Regional High-Frequency Seismic Wave Simulation in Realistic Three-Dimensional Earth Models to Improve Nuclear Explosion Monitoring''. A more appropriate title for this project is ''A Model-Based Signal Processing Approach to Nuclear Explosion Monitoring''. This project supported research for a radically new approach to nuclear explosion monitoring as well as allowed the development new capabilities in computational seismology that can contribute to NNSA/NA-22 Programs.
Date: March 14, 2007
Creator: Rodgers, A; Harris, D & Pasyanos, M
Partner: UNT Libraries Government Documents Department


Description: Comprehensive test ban monitoring in terms of location and discrimination has progressed significantly in recent years. However, the characterization of sources and the estimation of low yields remains a particular challenge. As the recent Korean shot demonstrated, we can probably expect to have a small set of teleseismic, far-regional and high-frequency regional data to analyze in estimating the yield of an event. Since stacking helps to bring signals out of the noise, it becomes useful to conduct comparable analyses on neighboring events, earthquakes in this case. If these auxiliary events have accurate moments and source descriptions, we have a means of directly comparing effective source strengths. Although we will rely on modeling codes, 1D, 2D, and 3D, we will also apply a broadband calibration procedure to use longer periods (P>5s) waveform data to calibrate short-period (P between .5 to 2 Hz) and high-frequency (P between 2 to 10 Hz) as path specify station corrections from well-known regional sources. We have expanded our basic Cut-and-Paste (CAP) methodology to include not only timing shifts but also amplitude (f) corrections at recording sites. The name of this method was derived from source inversions that allow timing shifts between 'waveform segments' (or cutting the seismogram up and re-assembling) to correct for crustal variation. For convenience, we will refer to these f-dependent refinements as CAP+ for (SP) and CAP++ for still higher frequency. These methods allow the retrieval of source parameters using only P-waveforms where radiation patterns are obvious as demonstrated in this report and are well suited for explosion P-wave data. The method is easily extended to all distances because it uses Green's function although there may be some changes required in t* to adjust for offsets between local vs. teleseismic distances. In short, we use a mixture of model-dependent and empirical corrections to ...
Date: July 16, 2007
Creator: Helmberger, D; Tromp, J & Rodgers, A
Partner: UNT Libraries Government Documents Department

Regionalization and calibration of seismic discriminants, path effects and signal-to-noise for station ABKT (Alibek, Turkmenistan)

Description: We report measurements and analysis of regional seismic phase amplitude ratios and signal-to-noise for earthquakes observed at the International Monitoring System primary station ABKT (Alibek, Turkmenistan). We measured noise and phase amplitudes of the regional phases Pn, Pg, Sn, and Lg in four frequency bands between 0.75-9.0 Hz. Measurements were made in both the time and frequency domains. The spatial variation of amplitude ratios (e.g., Pn/Lg, Pg/Lg, Pn/Sn, Pg/Sn) and signal-to-noise (phase/noise) reveal significant path effect differences between the Hindu Kush, Kazahk Platform, Iranian Plateau and Caspian Sea. In order to represent this behavior, we have investigated several techniques for characterizing the data. These techniques are: 1) correlation with along-path distance and waveguide properties; 2) sector analysis; and 3) spatial averaging. Along-path waveguide properties, such as mean elevation and rms topographic slope are found to be the strongest factors related to Pg/Lg amplitude ratios at the lowest frequencies (<3.0 Hz). Other path properties such as mean crustal thickness and basement depth are not strongly correlated with Pg/Lg ratios. For sector analysis we divided the data into four (4) azimuthal sectors and characterized the data within each sector by a distance trend. Sectors were chosen based on the behavior of Pn/Lg, Pg/Lg and Pn/Sn amplitude ratios as well as topographic and tectonic character. Results reveal significant reduction (up to a factor of two) in the scatter of the Pn/Lg and Pg/Lg amplitude ratios for the sectorized data compared to the entire data set from all azimuths. Spatial averaging involves smoothing and interpolation for the ratios projected at the event location. Methods such as cap averaging and kriging will be presented at the meeting. 7 refs., 6 figs.
Date: July 1, 1997
Creator: Rodgers, A.J. & Walter, W.R.
Partner: UNT Libraries Government Documents Department

Geophysical Model Applications for Monitoring

Description: Geophysical models constitute an important component of calibration for nuclear explosion monitoring. We will focus on four major topics and their applications: (1) surface wave models, (2) receiver function profiles, (3) regional tomography models, and (4) stochastic geophysical models. First, we continue to improve upon our surface wave model by adding more paths. This has allowed us to expand the region to all of Eurasia and into Africa, increase the resolution of our model, and extend results to even shorter periods (7 sec). High-resolution models exist for the Middle East and the YSKP region. The surface wave results can be inverted either alone, or in conjunction with other data, to derive models of the crust and upper mantle structure. One application of the group velocities is to construct phase-matched filters in combination with regional surface-wave magnitude formulas to improve the mb:Ms discriminant and extend it to smaller magnitude events. Next, we are using receiver functions, in joint inversions with the surface waves, to produce profiles directly under seismic stations throughout the region. In the past year, we have been focusing on deployments throughout the Middle East, including the Arabian Peninsula and Turkey. By assembling the results from many stations, we can see how regional seismic phases are affected by complicated upper mantle structure, including lithospheric thickness and anisotropy. The next geophysical model item, regional tomography models, can be used to predict regional travel times such as Pn and Sn. The times derived by the models can be used as a background model for empirical measurements or, where these don't exist, simply used as is. Finally, we have been exploring methodologies such as Markov Chain Monte Carlo (MCMC) to generate data-driven stochastic models. We have applied this technique to the YSKP region using surface wave dispersion data, body wave travel time ...
Date: July 11, 2005
Creator: Pasyanos, M; Walter, W; Tkalcic, H; Franz, G; Gok, R & Rodgers, A
Partner: UNT Libraries Government Documents Department

Pre-Shot Simulations of Far-Field Ground Motions for the Source Physics Experiment (SPE) Explosions at the Climax Stock, Nevada National Security Site

Description: The Source Physics Experiment (SPE) will involve a series of explosions in various geologic and emplacement conditions to validate numerical simulation methods to predict behavior of seismic wave excitation and propagation for nuclear test monitoring. The first SPE's currently underway involve explosions in the Climax Stock (granitic geology) at the Nevada National Security Site (NNSS). Detailed geologic data and published material properties for the major lithologic units of the NNSS and surrounding region were used to build three-dimensional models for seismic wave propagation simulations. The geologic structure near the SPE shot point is quite varied including granitic, carbonate, tuff and alluvium lithologies. We performed preliminary ground motion simulations for a near-source domain covering 8 km x 8 km at the surface centered on the shot point to investigate various source and propagation effects using WPP, LLNL's anelastic seismic wave finite difference code. Simulations indicate that variations in wave propagation properties of the sub-surface will generate strongly path-dependent response once the energy has left the relatively small granitic geology of the near-surface Climax Stock near the SPE shot point. Rough topography to the north and west of SPE shot point causes additional complexity in the signals including energy on the transverse components. Waves propagate much faster through the granitic and carbonate formations and slower through the tuff and alluvium. Synthetic seismograms for a pure explosion source in a 3D geologic structure show large amplitudes on transverse component. For paths to the south sampling the granite, tuff and alluvium lithologies transverse component amplitudes are as high as 50% of that on the vertical and radial components.
Date: November 7, 2010
Creator: Rodgers, A J; Wagoner, J; Petersson, N A & Sjogreen, B
Partner: UNT Libraries Government Documents Department

Tomographic Imaging of Upper Mantle P- and S-wave Velocity Heterogeneity Beneath the Arabian Peninsula

Description: We report the estimates of three-dimensional P- and S-wave velocity structure beneath the Arabian Peninsula estimated from travel time delay tomography. We have completed travel time measurements and inversion of a partial data set provided by King Abdulaziz City for Science and Technology (KACST). This study builds on previous work by Benoit et al. (2003) following the methods of VanDecar and Crosson (1990) and VanDecar (1991). Data were collected from the Saudi Arabian National Digital Seismic Network (SANDSN) operated by KACST. The network consists of 38 stations (27 broadband and 11 short-period). We augmented the KACST data with delay times measured from permanent Incorporated Research Institutions for Seismology (IRIS) stations in the region (RAYN, EIL and MRNI) and the 1996 Saudi Arabian PASSCAL Experiment. This study shows the inverted P- and S-wave models computed with the combined data with all three different seismic networks (KASCST, IRIS, and the 1996 Saudi Arabian PASSCAL experiment) with best coverage beneath the Arabian Shield. Tomographic images reveal low velocity features in the upper mantle along a north-south line from the southern Asir region to the northeastern portion of the Arabian Shield.
Date: August 30, 2005
Creator: Park, Y; Nyblade, A; Rodgers, A & Al-Amri, A
Partner: UNT Libraries Government Documents Department

The March 11, 2002 Masafi, United Arab Emirates Earthquake: Insights into the Seismotectonics of the Northern Oman Mountains

Description: A moderate (M{approx}5) earthquake struck the northeastern United Arab Emirates (UAE) and northern Oman on March 11, 2002. The event was felt over a wide area of the northern Emirates and was accompanied by smaller (felt) events before and after the March 11 main shock. The event was large enough to be detected and located by global networks at teleseismic distances. We estimated focal mechanism and depth from broadband complete regional waveform modeling. We report a normal mechanism with a slight right-lateral strike-slip component consistent with the large-scale tectonics. The normal component suggests relaxation of obducted crust of the Semail Ophilite (specifically, the Khor Fakkan Block) while the right-lateral strike-slip component of the mechanism is consistent with shear across the Oman Line. Felt earthquakes are rare in the region, however no regional seismic network exists in the UAE to determine local seismicity. This event offers a unique opportunity to study the active tectonics of the region as well as inform future studies of seismic hazard in the UAE and northern Oman.
Date: April 26, 2005
Creator: Rodgers, A; Fowler, A; Al-Amri, A & Al-Enezi, A
Partner: UNT Libraries Government Documents Department

Seismicity and Improved Velocity Structure in Kuwait

Description: The Kuwait National Seismic Network (KNSN) began operation in 1997 and consists of nine three-component stations (eight short-period and one broadband) and is operated by the Kuwait Institute for Scientific Research. Although the region is largely believed to be aseismic, considerable local seismicity is recorded by KNSN. Seismic events in Kuwait are clustered in two main groups, one in the south and another in the north. The KNSN station distribution is able to capture the southern cluster within the footprint of the network but the northern cluster is poorly covered. Events tend to occur at depths ranging from the free surface to about 20 km. Events in the northern cluster tend to be deeper than those in south, however this might be an artifact of the station coverage. We analyzed KNSN recordings of nearly 200 local events to improve understanding of seismic events and crustal structure in Kuwait, performing several analyses with increasing complexity. First, we obtained an optimized one-dimensional (1D) velocity model for the entire region using the reported KNSN arrival times and routine locations. The resulting model is consistent with a recently obtained model from the joint inversion of receiver functions and surface wave group velocities. Crustal structure is capped by the thick ({approx} 7 km) sedimentary rocks of the Arabian Platform underlain by normal velocities for stable continental crust. Our new model has a crustal thickness of 44 km, constrained by an independent study of receiver functions and surface wave group velocities by Pasyanos et al (2006). Locations and depths of events after relocation with the new model are broadly consistent with those reported by KISR, although a few events move more than a few kilometers. We then used a double-difference tomography technique (tomoDD) to jointly locate the events and estimate three-dimensional (3D) velocity structure. TomoDD is ...
Date: January 26, 2006
Creator: Gok, R M; Rodgers, A J & Al-Enezi, A
Partner: UNT Libraries Government Documents Department

Crustal Structure of Iraq from Receiver Functions and Surface Wave Dispersion

Description: We report the crustal structure of Iraq, located in the northeastern Arabian plate, estimated by joint inversion of P-wave receiver functions and surface wave group velocity dispersion. Receiver functions were computed from teleseismic recordings at two temporary broadband seismic stations in Mosul (MSL) and Baghdad (BHD), separated by approximately 360 km. Group velocity dispersion curves at the sites were derived from continental-scale tomography of Pasyanos (2006). The inversion results show that the crustal thicknesses are 39 km at MSL and 43 km at BHD. Both sites reveal low velocity surface layers consistent with sedimentary thickness of about 3 km at station MSL and 7 km at BHD, agreeing well with the existing models. Ignoring the sediments, the crustal velocities and thicknesses are remarkably similar between the two stations, suggesting that the crustal structure of the proto-Arabian Platform in northern Iraq was uniform before subsidence and deposition of the sediments in the Cenozoic. Deeper low velocity sediments at BHD are expected to result in higher ground motions for earthquakes.
Date: August 31, 2006
Creator: Gok, R; Mahdi, H; Al-Shukri, H & Rodgers, A J
Partner: UNT Libraries Government Documents Department

Model-Based Signal Processing: Correlation Detection With Synthetic Seismograms

Description: Recent applications of correlation methods to seismological problems illustrate the power of coherent signal processing applied to seismic waveforms. Examples of these applications include detection of low amplitude signals buried in ambient noise and cross-correlation of sets of waveforms to form event clusters and accurately measure delay times for event relocation and/or earth structure. These methods rely on the exploitation of the similarity of individual waveforms and have been successfully applied to large sets of empirical observations. However, in cases with little or no empirical event data, such as aseismic regions or exotic event types, correlation methods with observed seismograms will not be possible due to the lack of previously observed similar waveforms. This study uses model-based signals computed for three-dimensional (3D) Earth models to form the basis for correlation detection. Synthetic seismograms are computed for fully 3D models estimated from the Markov Chain Monte-Carlo (MCMC) method. MCMC uses stochastic sampling to fit multiple seismological data sets. Rather than estimate a single ''optimal'' model, MCMC results in a suite of models that sample the model space and incorporates uncertainty through variability of the models. The variability reflects our ignorance of Earth structure, due to limited resolution, data and modeling errors, and produces variability in the seismic waveform response. Model-based signals are combined using a subspace method where the synthetic signals are decomposed into an orthogonal basis by singular-value decomposition (SVD) and the observed waveforms are represented with a linear combination of a sub-set of eigenvectors (signals) associated with the most significant eigenvalues. We have demonstrated the method by modeling long-period (80-10 seconds) regional seismograms for a moderate (M{approx}5) earthquake near the China-North Korea border. Synthetic seismograms are computed with the Spectral Element Method for a suite of long-wavelength (2 degree) seismic velocity models based on the MCMC method. We are ...
Date: August 30, 2006
Creator: Rodgers, A; Harris, D; Pasyanos, M; Blair, S & Matt, R
Partner: UNT Libraries Government Documents Department

Combined Plate Motion and Density Driven Flow in the Asthenosphere beneath Saudi Arabia: Evidence from Shearwave Splitting and Seismic Anisotropy

Description: Mantle anisotropy along the Red Sea and across the Arabian Peninsula was analyzed using shear-wave splitting recorded by stations from three different seismic networks: the largest, most widely distributed array of stations examined across the Arabian Peninsula to date. Stations near the Gulf of Aqaba display fast orientations aligned parallel to the Dead Sea Transform Fault, most likely related to the strike-slip motion between Africa and Arabia However, most of our observations across Arabia are statistically the same (at a 95% confidence level), with north-south oriented fast directions and delay times averaging about 1.4 s. Since end-member models of fossilized anisotropy and present-day asthenospheric flow do not adequately explain these observations, we interpret them as a combination of plate and density driven flow in the asthenosphere. Combining northeast oriented flow associated with absolute plate motion with northwest oriented flow associated with the channelized Afar upwelling along the Red Sea produces a north-south resultant that matches the observations and supports models of active rifting.
Date: September 8, 2006
Creator: Hansen, S; Schwartz, S; Al-Amri, A & Rodgers, A
Partner: UNT Libraries Government Documents Department

Modeling the conversion of hydroacoustic to seismic energy at island and continental margins: preliminary analysis of Ascension Island data

Description: Seismic stations at islands and continental margins will be an essential component of the International Monitoring System (IMS) for event location and identification in support of Comprehensive Nuclear-Test-Ban Treaty (CTBT) monitoring. Particularly important will be the detection and analysis of hydroacoustic-to-seismic converted waves (T-phases) at island or continental margins. Acoustic waves generated by sources in or near the ocean propagate for long distances very efficiently due to the ocean sound speed channel (SOFAR) and low attenuation. When ocean propagating acoustic waves strike an island or continental margin they are converted to seismic (elastic) waves. We are using a finite difference code to model the conversion of hydroacoustic T-waves at an island or continental margin. Although ray-based methods are far more efficient for modeling long-range (&gt; 1000 km) high-frequency hydroacoustic propagation, the finite difference method has the advantage of being able to model both acoustic and elastic wave propagation for a broad range of frequencies. The method allows us to perform simulations of T-phases to relatively high frequencies ({&gt;=}10 Hz). Of particular interest is to identify factors that affect the efficiency of T-phase conversion, such as the topographic slope and roughness at the conversion point and elastic velocity structure within the island or continent. Previous studies have shown that efficient T-phase conversion occurs when the topographic slope at the conversion point is steep (Cansi and Bethoux, 1985; Talandier and Okal, 1998). Another factor impacting T-phase conversion may be the near-shore structure of the sound channel. It is well known that the depth to the sound channel axis decreases in shallow waters. This can weaken the channeled hydroacoustic wave. Elastic velocity structure within the island or continent will impact how the converted seismic wave is refracted to recording stations at the surface and thus impact the T-phase amplitudes. For this paper we ...
Date: July 26, 1999
Creator: Harben, P. & Rodgers, A.
Partner: UNT Libraries Government Documents Department