Long period seismic waves from large, near-surface nuclear explosions

1963 ◽  
Vol 53 (1) ◽  
pp. 109-149 ◽  
Author(s):  
Paul W. Pomeroy

Abstract Seismic surface waves were well-recorded from the larger explosions of the U.S. test series detonated in the Marshall Islands during the spring and summer of 1958 and of the U. S. S. R. test series detonated at the Novaya Zemlaya test site during October 1958. In addition to waves of the fundamental Rayleigh mode, some unusual or unexpected waves were identified. These include: 1.) Love waves with lengths as great as 90 km at Hong Kong from the U. S. S. R. explosions, 2.) Love waves at Agra, India, from the U. S. S. R. explosions, 3.) waves of the 1st shear mode at Agra and Uppsala, Sweden, from the U. S. S. R. explosions, and 4.) Love, 1st shear and PL waves at Guam from the U. S. explosions. Group velocity data were derived for many paths and, for the U. S. series, these data are very accurate because they are based on exact knowledge of locations and origin times. For the U. S. S. R. explosions, revised locations and origin times based on a limited number of P-wave observations were used to determine group velocities. Using records from Hong Kong and Honolulu for two U. S. tests, one at Eniwetok and one at Bikini, oceanic phase velocities in the period range of 15 to 40 seconds were measured for the path between these islands. Using Brune's method for initial phase determination with a known phase velocity, an initial phase, φ0, lying between +3π2 and +π2 was determined for the U.S. explosion, Oak. According to Aki (1960), these values of initial phase are associated, respectively, with forcing functions of a downward impulse and an explosive impulse. Seismic magnitudes of 4.7 and 4.8 were assigned to the U. S. Oak and Poplar events on the basis of the surface wave data, while magnitudes of the larger U. S. S. R. tests range from about 4.0 to 4.5. The ratio of seismic energy as computed from the surface wave magnitude to the total explosive energy available (yield) is apparently greater for the U. S. explosions, assuming equal yield for the two shots compared. This suggests a higher altitude of detonation for the Russian events. The Fourier amplitude spectra of the surface wave trains are generally single peaked where a significant portion of the path traversed is oceanic, but for Uppsala, where the path is short and continental, a second peak occurs at periods of about 10 to 13 seconds. The ratios of the predicted amplitudes, based on the data of one station and on a reasonable dissipation factor, to the actual recorded amplitudes vary with azimuth by as much as a factor of 5. These variations may be explained by asymmetry at the source but might also be explained by instrumental and geologic factors. Long period components of P and S and multiples thereof were sometimes recorded from these events.

1990 ◽  
Vol 80 (5) ◽  
pp. 1205-1231
Author(s):  
Jiajun Zhang ◽  
Thorne Lay

Abstract Determination of shallow earthquake source mechanisms by inversion of long-period (150 to 300 sec) Rayleigh waves requires epicentral locations with greater accuracy than that provided by routine source locations of the National Earthquake Information Center (NEIC) and International Seismological Centre (ISC). The effects of epicentral mislocation on such inversions are examined using synthetic calculations as well as actual data for three large Mexican earthquakes. For Rayleigh waves of 150-sec period, an epicentral mislocation of 30 km introduces observed source spectra phase errors of 0.6 radian for stations at opposing azimuths along the source mislocation vector. This is larger than the 0.5-radian azimuthal variation of the phase spectra at the same period for a thrust fault with 15° dip and 24-km depth. The typical landward mislocation of routinely determined epicenters of shallow subduction zone earthquakes causes source moment tensor inversions of long-period Rayleigh waves to predict larger fault dip than indicated by teleseismic P-wave first-motion data. For dip-slip earthquakes, inversions of long-period Rayleigh waves that use an erroneous source location in the down-dip or along-strike directions of a nodal plane, overestimate the strike, dip, and slip of that nodal plane. Inversions of strike-slip earthquakes that utilize an erroneous location along the strike of a nodal plane overestimate the slip of that nodal plane, causing the second nodal plane to dip incorrectly in the direction opposite to the mislocation vector. The effects of epicentral mislocation for earthquakes with 45° dip-slip fault mechanisms are more severe than for events with other fault mechanisms. Existing earth model propagation corrections do not appear to be sufficiently accurate to routinely determine the optimal surface-wave source location without constraints from body-wave information, unless extensive direct path (R1) data are available or empirical path calibrations are performed. However, independent surface-wave and body-wave solutions can be remarkably consistent when the effects of epicentral mislocation are accounted for. This will allow simultaneous unconstrained body-wave and surface-wave inversions to be performed despite the well known difficulties of extracting the complete moment tensor of shallow sources from fundamental modes.


2015 ◽  
Vol 58 (5) ◽  
Author(s):  
Sankar N. Bhattacharya

<p>Sensitivity kernels or partial derivatives of phase velocity (<em>c</em>) and group velocity (<em>U</em>) with respect to medium parameters are useful to interpret a given set of observed surface wave velocity data. In addition to phase velocities, group velocities are also being observed to find the radial anisotropy of the crust and mantle. However, sensitivities of group velocity for a radially anisotropic Earth have rarely been studied. Here we show sensitivities of group velocity along with those of phase velocity to the medium parameters <em>V<sub>SV</sub>, V<sub>SH </sub>, V<sub>PV</sub>, V<sub>PH , </sub></em><em>h</em><em> </em>and density in a radially anisotropic spherical Earth. The peak sensitivities for <em>U</em> are generally twice of those for <em>c</em>; thus <em>U</em> is more efficient than <em>c</em> to explore anisotropic nature of the medium. Love waves mainly depends on <em>V<sub>SH</sub></em> while Rayleigh waves is nearly independent of <em>V<sub>SH</sub></em> . The sensitivities show that there are trade-offs among these parameters during inversion and there is a need to reduce the number of parameters to be evaluated independently. It is suggested to use a nonlinear inversion jointly for Rayleigh and Love waves; in such a nonlinear inversion best solutions are obtained among the model parameters within prescribed limits for each parameter. We first choose <em>V<sub>SH</sub></em>, <em>V<sub>SV </sub></em>and <em>V<sub>PH</sub></em> within their corresponding limits; <em>V<sub>PV</sub></em> and <em>h</em> can be evaluated from empirical relations among the parameters. The density has small effect on surface wave velocities and it can be considered from other studies or from empirical relation of density to average P-wave velocity.</p>


2021 ◽  
Vol 40 (8) ◽  
pp. 567-575
Author(s):  
Myrto Papadopoulou ◽  
Farbod Khosro Anjom ◽  
Mohammad Karim Karimpour ◽  
Valentina Laura Socco

Surface-wave (SW) tomography is a technique that has been widely used in the field of seismology. It can provide higher resolution relative to the classical multichannel SW processing and inversion schemes that are usually adopted for near-surface applications. Nevertheless, the method is rarely used in this context, mainly due to the long processing times needed to pick the dispersion curves as well as the inability of the two-station processing to discriminate between higher SW modes. To make it efficient and to retrieve pseudo-2D/3D S-wave velocity (VS) and P-wave velocity (VP) models in a fast and convenient way, we develop a fully data-driven two-station dispersion curve estimation, which achieves dense spatial coverage without the involvement of an operator. To handle higher SW modes, we apply a dedicated time-windowing algorithm to isolate and pick the different modes. A multimodal tomographic inversion is applied to estimate a VS model. The VS model is then converted to a VP model with the Poisson's ratio estimated through the wavelength-depth method. We apply the method to a 2D seismic exploration data set acquired at a mining site, where strong lateral heterogeneity is expected, and to a 3D pilot data set, recorded with state-of-the-art acquisition technology. We compare the results with the ones retrieved from classical multichannel analysis.


Geophysics ◽  
2017 ◽  
Vol 82 (6) ◽  
pp. WB47-WB61 ◽  
Author(s):  
Sylvain Pasquet ◽  
Ludovic Bodet

The simultaneous estimation of 2D pressure (P-) and S-wave velocities ([Formula: see text] and [Formula: see text], respectively) is a promising approach for imaging subsurface mechanical properties. It can be performed with a single acquisition setup by combining P-wave refraction and surface-wave (SW) analysis. Although SW methods are commonly applied for the 1D estimation of [Formula: see text], 2D profiling requires the implementation of specific processing and inversion tools not yet widely available in the community. We have developed an open-source MATLAB-based package that performs SW inversion and profiling (SWIP) so as to retrieve 1D to 2D variations of [Formula: see text] from any kind of linear active-source near-surface seismic data. Each step of the workflow involves up-to-date processing and inversion techniques and provides ready-to-use outputs with quality control tools. First, windowing and stacking techniques are implemented to enhance the signal-to-noise ratio and extract local dispersion images along the line. Then, dispersion curves are picked for each window with an uncertainty range in the phase velocity including higher uncertainties at low frequency. These curves are next inverted using a Monte Carlo approach with various parameterizations (e.g., user defined, refraction based). The best models are finally selected according to their fit to the data to build an average final model with a suggested investigation depth. As an example, we used SWIP to process data collected at a Yellowstone hydrothermal system. Our results show the benefits of estimating [Formula: see text] and [Formula: see text] from a single seismic setup to highlight subsurface gas pathways.


Geophysics ◽  
2018 ◽  
Vol 83 (6) ◽  
pp. B339-B353 ◽  
Author(s):  
Bojan Brodic ◽  
Alireza Malehmir ◽  
André Pugin ◽  
Georgiana Maries

We deployed a newly developed 3C microelectromechanical system-based seismic land streamer over porous glacial sediments to delineate water table and bedrock in Southwestern Finland. The seismic source used was a 500 kg vertical impact drop hammer. We analyzed the SH-wave component and interpreted it together with previously analyzed P-wave component data. In addition to this, we examined the land streamer’s potential for multichannel analysis of surface waves and delineated the site’s stratigraphy with surface-wave-derived S-wave velocities and [Formula: see text] ratios along the entire profile. These S-wave velocities and [Formula: see text] ratios complement the interpretation conducted previously on P-wave stacked section. Peculiarly, although the seismic source used is of a vertical-type nature, the data inspection indicated clear bedrock reflection on the horizontal components, particularly the transverse component. This observation led us to scrutinize the horizontal component data through side-by-side inspection of the shot records of all the three components and particle motion analysis to confirm the S-wave nature of the reflection. Using the apparent moveout velocity of the reflection, as well as the known depth to bedrock based on drilling, we used finite-difference synthetic modeling to further verify its nature. Compared with the P-wave seismic section, bedrock is relatively well delineated on the transverse component S-wave section. Some structures connected to the kettle holes and other stratigraphic units imaged on the P-wave results were also notable on the S-wave section, and particularly on the surface-wave derived S-wave velocity model and [Formula: see text] ratios. Our results indicate that P-, SV-, and SH-wave energy is generated simultaneously at the source location itself. This study demonstrates the potential of 3C seismic for characterization and delineation of the near-surface seismics.


1961 ◽  
Vol 51 (1) ◽  
pp. 1-12
Author(s):  
Jack Oliver

Abstract S and multiple S phases at moderate to large epicentral distances are frequently followed by normally-dispersed, long-period, wave trains for which surface particle motion is elliptical and progressive and in the plane of propagation of the SV wave. The character of such phases can be explained as the result of coupling between the incident shear waves and dispersive PL waves in the near-surface wave guide. A detailed study of shocks in Mexico and in Montana recorded at Resolute, and less detailed studies of other data support this hypothesis.


2019 ◽  
Vol 113 (4) ◽  
pp. 818-821

In July of 2019, the U.S. Senate gave advice and consent to protocols updating tax treaties with Spain, Switzerland, Japan, and Luxembourg, after a nearly decade-long period during which no tax treaties were approved by the Senate. This drought was primarily due to the privacy concerns of a single senator, Rand Paul of Kentucky, who deployed the Senate's procedural rules to increase the difficulty of the advice and consent process. Tax treaties with Hungary, Chile, and Poland, as well as a protocol to a multilateral tax convention, remained pending in the Senate Foreign Relations Committee as of mid-August of 2019.


2019 ◽  
Vol 24 (1) ◽  
pp. 101-120
Author(s):  
Kajetan Chrapkiewicz ◽  
Monika Wilde-Piórko ◽  
Marcin Polkowski ◽  
Marek Grad

AbstractNon-linear inverse problems arising in seismology are usually addressed either by linearization or by Monte Carlo methods. Neither approach is flawless. The former needs an accurate starting model; the latter is computationally intensive. Both require careful tuning of inversion parameters. An additional challenge is posed by joint inversion of data of different sensitivities and noise levels such as receiver functions and surface wave dispersion curves. We propose a generic workflow that combines advantages of both methods by endowing the linearized approach with an ensemble of homogeneous starting models. It successfully addresses several fundamental issues inherent in a wide range of inverse problems, such as trapping by local minima, exploitation of a priori knowledge, choice of a model depth, proper weighting of data sets characterized by different uncertainties, and credibility of final models. Some of them are tackled with the aid of novel 1D checkerboard tests—an intuitive and feasible addition to the resolution matrix. We applied our workflow to study the south-western margin of the East European Craton. Rayleigh wave phase velocity dispersion and P-wave receiver function data were gathered in the passive seismic experiment “13 BB Star” (2013–2016) in the area of the crust recognized by previous borehole and refraction surveys. Final models of S-wave velocity down to 300 km depth beneath the array are characterized by proximity in the parameter space and very good data fit. The maximum value in the mantle is higher by 0.1–0.2 km/s than reported for other cratons.


Geophysics ◽  
2013 ◽  
Vol 78 (3) ◽  
pp. D169-D179 ◽  
Author(s):  
Zijian Zhang ◽  
De-hua Han ◽  
Daniel R. McConnell

Hydrate-bearing sands and shallow nodular hydrate are potential energy resources and geohazards, and they both need to be better understood and identified. Therefore, it is useful to develop methodologies for modeling and simulating elastic constants of these hydrate-bearing sediments. A gas-hydrate rock-physics model based on the effective medium theory was successfully applied to dry rock, water-saturated rock, and hydrate-bearing rock. The model was used to investigate the seismic interpretation capability of hydrate-bearing sediments in the Gulf of Mexico by computing elastic constants, also known as seismic attributes, in terms of seismic interpretation, including the normal incident reflectivity (NI), Poisson’s ratio (PR), P-wave velocity ([Formula: see text]), S-wave velocity ([Formula: see text]), and density. The study of the model was concerned with the formation of gas hydrate, and, therefore, hydrate-bearing sediments were divided into hydrate-bearing sands, hydrate-bearing sands with free gas in the pore space, and shallow nodular hydrate. Although relations of hydrate saturation versus [Formula: see text] and [Formula: see text] are different between structures I and II gas hydrates, highly concentrated hydrate-bearing sands may be interpreted on poststack seismic amplitude sections because of the high NI present. The computations of elastic constant implied that hydrate-bearing sands with free gas could be detected with the crossplot of NI and PR from prestack amplitude analysis, and density may be a good hydrate indicator for shallow nodular hydrate, if it can be accurately estimated by seismic methods.


1990 ◽  
Vol 80 (6A) ◽  
pp. 1677-1695 ◽  
Author(s):  
Ik Bum Kang ◽  
George A. McMechan

Abstract Full wave field modeling of wide-aperture data is performed with a pseudospectral implementation of the elastic wave equation. This approach naturally produces three-component stress and two-component particle displacement, velocity, and acceleration seismograms for compressional, shear, and Rayleigh waves. It also has distinct advantages in terms of computational requirements over finite-differencing when data from large-scale structures are to be modeled at high frequencies. The algorithm is applied to iterative two-dimensional modeling of seismograms from a survey performed in 1985 by The University of Texas at El Paso and The University of Texas at Dallas across the Anadarko basin and the Wichita Mountains in southwestern Oklahoma. The results provide an independent look at details of near-surface structure and reflector configurations. Near-surface (&lt;3 km deep) structure and scattering effects account for a large percentage (&gt;70 per cent) of the energy in the observed seismograms. The interpretation of the data is consistent with the results of previous studies of these data, but provides considerably more detail. Overall, the P-wave velocities in the Wichita Uplift are more typical of the middle crust than the upper crust (5.3 to 7.1 km/sec). At the surface, the uplift is either exposed as weathered outcrop (5.0 to 5.3 km/sec) or is overlain with sediments of up to 0.4 km in thickness, ranging in velocity from 2.7 to 3.4 km/sec, generally increasing with depth. The core of the uplift is relatively seismically transparent. A very clear, coherent reflection is observed from the Mountain View fault, which dips at ≈40° to the southwest, to at least 12 km depth. Velocities in the Anadarko Basin are typical of sedimentary basins; there is a general increase from ≈2.7 km/sec at the surface to ≈5.9 km/sec at ≈16 km depth, with discontinuous reflections at depths of ≈8, 10, 12, and 16 km.


Sign in / Sign up

Export Citation Format

Share Document