Sustained Performance of 10+ Teraflop/s in Simulation on Seismic Waves Using 507 Nodes of the Earth Simulator

Author(s):  
Seiji Tsuboi
2001 ◽  
Vol 09 (04) ◽  
pp. 1407-1416 ◽  
Author(s):  
GIULIANA ROSSI ◽  
ALDO VESNAVER

Converted waves can play a basic role in the traveltime inversion of seismic waves. The sought velocity fields of P and S waves are almost decoupled, when considering pure P and S arrivals: their only connection are the possible common reflecting interfaces in the Earth. Converted waves provide new equations in the linear system to be inverted, which directly relates the two velocity fields. Since the new equations do not introduce additional unknowns, they increase the system rank or its redundancy, so making its solutions better constrained and robust.


Author(s):  
John A. Adam

This chapter focuses on the underlying mathematics of seismic rays. Seismic waves caused by earthquakes and explosions are used in seismic tomography to create computer-generated, three-dimensional images of Earth's interior. If the Earth had a uniform composition and density, seismic rays would travel in straight lines. However, it is broadly layered, causing seismic rays to be refracted and reflected across boundaries. In order to calculate the speed along the wave's ray path, the time it takes for a seismic wave to arrive at a seismic station from an earthquake needs to be determined. Arrival times of different seismic waves allow scientists to define slower or faster regions deep in the Earth. The chapter first presents the relevant equations for seismic rays before discussing how rays are propagated in a spherical Earth. The Wiechert-Herglotz inverse problem is considered, along with the properties of X in a horizontally stratified Earth.


Author(s):  
Olga Popova

The asteroid impact near the Russian city of Chelyabinsk on February 15, 2013, was the largest airburst on Earth since the 1908 Tunguska event, causing a natural disaster in an area with a population exceeding 1 million. On clear morning at 9:20 a.m. local time, an asteroid about 19 m in size entered the Earth atmosphere near southern Ural Mountains (Russia) and, with its bright illumination, attracted the attention of hundreds of thousands of people. Dust trail in the atmosphere after the bolide was tens of kilometers long and was visible for several hours. Thousands of different size meteorites were found in the areas south-southwest of Chelyabinsk. A powerful airburst, which was formed due to meteoroid energy deposition, shattered thousands of windows and doors in Chelyabinsk and wide surroundings, with flying glass injuring many residents. The entrance and destruction of the 500-kt Chelyabinsk asteroid produced a number of observable effects, including light and thermal radiation; acoustic, infrasound, blast, and seismic waves; and release of interplanetary substance. This unexpected and unusual event is the most well-documented bolide airburst, and it attracted worldwide attention. The airburst was observed globally by multiple instruments. Analyses of the observational data allowed determination of the size of the body that caused the superbolide, its velocity, its trajectory, its behavior in the atmosphere, the strength of the blast wave, and other characteristics. The entry of the 19-m-diameter Chelyabinsk asteroid provides a unique opportunity to calibrate the different approaches used to model meteoroid entry and to calculate the damaging effects. The recovered meteorite material was characterized as brecciated LL5 ordinary chondrite, in which three different lithologies can be distinguished (light-colored, dark-colored, and impact-melt). The structure and properties of meteorites demonstrate that before encountering Earth, the Chelyabinsk asteroid had experienced a very complex history involving at least a few impacts with other bodies and thermal metamorphism. The Chelyabinsk airburst of February 15, 2013, was exceptional because of the large kinetic energy of the impacting body and the damaging airburst that was generated. Before the event, decameter-sized objects were considered to be safe. With the Chelyabinsk event, it is possible, for the first time, to link the damage from an impact event to a well-determined impact energy in order to assess the future hazards of asteroids to lives and property.


Geophysics ◽  
1985 ◽  
Vol 50 (11) ◽  
pp. 1684-1688 ◽  
Author(s):  
S. E. Pullan ◽  
J. A. Hunter

One aspect of the modification of seismic waves on passage through the Earth is the partitioning of energy at subsurface interfaces as described by the Zoeppritz equations. These equations have been applied to simple two‐ and three‐layer models to determine the variations in the amplitude and phase of a reflection signal at nonnormal angles of incidence. Synthetic seismograms have been produced to illustrate the effect of these variations on a seismic wavelet. It is found that the phase variations can lead to substantial changes in the character of the reflected wavelet as the source‐geophone distance increases. These changes are dependent on the angle of incidence and on the elastic properties of the subsurface layers. In particular, for models approximating an overburden over bedrock situation, the reflected pulse is predicted to “change phase” when the velocity contrast between the two layers is relatively small. This effect has been observed on field records. Geophysicists should be aware of this phenomenon, both in terms of the problems it may cause in observing a bedrock reflection and in terms of the potential it has for indicating subsurface elastic properties.


2016 ◽  
Vol 16 (04) ◽  
pp. 1650009
Author(s):  
Josselin Garnier ◽  
Knut Sølna

Waves propagating through heterogeneous media experience scattering that can convert a coherent pulse into small incoherent fluctuations. This may appear as attenuation for the transmitted front pulse. The classic O’Doherty–Anstey theory describes such a transformation for scalar waves in finely layered media. Recent observations for seismic waves in the earth suggest that this theory can explain a significant component of seismic attenuation. An important question to answer is then how the O’Doherty–Anstey theory generalizes to seismic waves when several wave modes, possibly with the same velocity, interact. An important aspect of the O’Doherty–Anstey theory is the statistical stability property, which means that the transmitted front pulse is actually deterministic and depends only on the statistics of the medium but not on the particular medium realization when the medium is modeled as a random process. It is shown in this paper that this property generalizes in the case of elastic waves in a nontrivial way: the energy of the transmitted front pulse, but not the pulse shape itself, is statistically stable. This result is based on a separation of scales technique and a diffusion-approximation theorem that characterize the transmitted front pulse as the solution of a stochastic partial differential equation driven by two Brownian motions.


2014 ◽  
Vol 6 (1) ◽  
pp. 1371-1400 ◽  
Author(s):  
J. M. Carcione ◽  
F. Poletto ◽  
B. Farina ◽  
A. Craglietto

Abstract. The Earth crust presents two dissimilar rheological behaviours depending on the in-situ stress-temperature conditions. The upper, cooler, part is brittle while deeper zones are ductile. Seismic waves may reveal the presence of the transition but a proper characterization is required. We first obtain a stress–strain relation including the effects of shear seismic attenuation and ductility due to shear deformations and plastic flow. The anelastic behaviour is based on the Burgers mechanical model to describe the effects of seismic attenuation and steady-state creep flow. The shear Lamé constant of the brittle and ductile media depends on the in-situ stress and temperature through the shear viscosity, which is obtained by the Arrhenius equation and the octahedral stress criterion. The P- and S-wave velocities decrease as depth and temperature increase due to the geothermal gradient, an effect which is more pronounced for shear waves. We then obtain the P-S and SH equations of motion recast in the velocity-stress formulation, including memory variables to avoid the computation of time convolutions. The equations correspond to isotropic anelastic and inhomogeneous media and are solved by a direct grid method based on the Runge–Kutta time stepping technique and the Fourier pseudospectral method. The algorithm is tested with success against known analytical solutions for different shear viscosities. A realistic example illustrates the computation of surface and reverse-VSP synthetic seismograms in the presence of an abrupt brittle-ductile transition.


Sign in / Sign up

Export Citation Format

Share Document