Model of a strong volcanic blast and a method of estimating the mass ejected

A hydrodynamic model of a strong volcanic blast associated with large (caldera-forming) ash–gas-flow eruptions is presented. We show that the mass erupted from a vent influences considerably the damping of the shock wave. If the ejected mass flux is constant, the shock amplitude as a function of distance will decrease slower (in accordance with the inverse 4/3 power law) than that of an atomic explosion (the inverse cube power law). From the modelling, we extract a relation connecting the effective energy of a blast with the mass ejected during the blast phase. The mass released during the largest explosion of the 1883 eruption of Krakatau, and the lateral blast energy at Mount St. Helens on 1980 May 18 are estimated on the basis of the equation obtained.     

Elastic scattered waves from a continuous and heterogeneous layer

Elastic scattering from a continuous and laterally unbounded heterogeneous layer has been formulated using the Born approximation. A general solution of the scattered wave equation for the above-stated medium has been given in terms of a Fourier integral over plane waves. Far-field asymptotic expressions for weak elastic scattering by a finite, continuous and inhomogeneous layer have been presented which agree with earlier results. For perturbations of the two elastic parameters and the density having the same form of spatial variation, the spectrum of plane waves scattered from a heterogeneous layer is expressed as a product of an 'elastic scattering factor'and a 'distribution factor'. As in earlier results for small-scale heterogeneity, the scattering pattern depends on various combinations of perturbations of elastic parameters and density. In order to show the general characteristics of the elastic wave scattering, some scattering patterns have been given.     

Extended formulation for Rayleigh-wave computations at high frequencies in a spherical layered earth

Summary. The earlier formulation of dispersion function, amplitude response and surface ellipticity for Rayleigh waves in a spherical layered earth has been extended by the reduced-delta matrix to avoid computational instability at high frequencies. We derive a modified formulation, where the layer thickness is kept separated during matrix multiplications. The present formulations are valid for an unlimited frequency range. A numerical experiment with single precision (six decimal digits accuracy) on an IBM 360/44 computer is performed at frequencies from 0.4 to 12.5 Hz on a Shield structure down to a depth of 1090 km below which a complete sphere is assumed. The experiment shows that both reduced-delta matrix and modified formulations generate roots of the dispersion equation without any reduction of layers. However, to avoid unnecessary computation, a layer-reduction procedure has been derived from a modified formulation when the contribution of deep layers becomes insignificant on surface waves.     

近45年全球海域波浪能资源研究及等级区划

基于来自ECMWF的资料将1957年09月---2002年08月风浪、涌浪分离的ERA-40 wave reanalysis,对全球海域的波浪能资源进行重新审视．充分依据涌浪具有能量大、稳定性好等优点,利于波浪能的采集与转换,从提高波浪能资源有效利用率的角度出发,综合考虑能流密度的大小、能级频率、能流密度的稳定性和长期变化趋势等,对全球海域的波浪能资源进行系统性研究,并构建一套波浪能资源评估系统,对全球海域的波浪能资源进行功能区划,为海浪发电、海水淡化等波浪能资源开发工作提供科学依据．     

南极绕极波的行波与驻波共存系统分析

南极绕极波(ACW)是近10年来海洋和大气科学领域研究的重要内容之一,其特征表现为一种海气耦合的年际信号在南大洋的明显东向传播。文中依据波动的基本原理,利用非线性拟合方法得到了绕极波行波的平均波速、波数和角频率和印度洋驻波的表达形式,得出行波和驻波具有相同的角频率。针对绕极波纬圈传播的不连续现象,提出了南极绕极波是行波和驻波共存系统的新见解。行波与驻波共存系统由三部分组成:在南大洋的太平洋西侧发生SLP振动,与发生在太平洋中低纬度的过程有关;在太平洋和大西洋存在传播的行波,行波在太平洋最强,到大西洋明显减弱;在印度洋SLP系统表现为局地变化,体现明显的驻波特征。三部分有共同的周期,表明三部分是统一的动力学系统,证实了行波与驻波共存系统可以很好地表达南极绕极波的时空特性和传播特征。     

2.5-dimensional crosshole acoustic response: a variable density approach

This paper presents the development of a 2.5-D simulation technique for acoustic wave propagation in media with variable density and velocity. A comparative study of the 2-D and 2.5-D responses of a model reveals the spatially and temporally damped nature of the 2.5-D acoustic wave equations. The simulated results for constant and variable density models show that the density variation affects only the reflectivity of the layer. The computational cost for variable density models is 2.17 and 2.26 times that for constant density models for the 2.5-D and 2-D cases, respectively. Furthermore, the 2.5-D computational cost in the time domain is only about 10–15 per cent more than that for two dimensions, so this modest increase in computational cost can avoid the exorbitant 3-D computational cost.
Snapshots for a crosshole geometry were computed at various times in order to study the effect of heterogeneity on the amplitude and shape of the wave front. Extensive analysis of an oil-bearing reservoir with and without the inclusion of a gas zone was performed using a point source as well as multiple sources. In addition, the effects of the thickness of a low-velocity layer (oil-bearing) and of the location of the source have been studied. It is concluded from the numerical response that the waveguide action of the low-velocity layer depends on its thickness in terms of the dominant wavelength. Trapping of waves was not observed when the source was outside the low-velocity layer. Furthermore, the presence of heterogeneity in the low-velocity layer contributes considerably to the leakage of energy in the adjacent layers due to scattering/diffraction. It was found that, in the 2.5-D numerical simulation, the stability condition and the requirement of the number of grid points per wavelength to avoid grid dispersion are the same as for the 2-D case.     

Teleseismic waveform modelling with a one-way wave equation

It is now widely accepted that elastic properties of the continental lithosphere and the underlying sublithospheric mantle are both anisotropic and laterally heterogeneous at a range of scales. To fully exploit modern three-component broad-band array data sets requires the use of comprehensive modelling tools. In this work, we investigate the use of a wide-angle, one-way wave equation to model variations in teleseismic 3-D waveforms due to 2-D elastic heterogeneity and anisotropy. The one-way operators are derived based on a high-frequency approximation of the square-root operator and include the effects of wave propagation as well as multiple scattering. Computational cost is reduced through a number of physically motivated approximations. We present synthetic results from simple 1-D (layer over a half-space) and 2-D (subduction zone) models that are compared with reference solutions. The algorithm is then used to model data from an array of broad-band seismograph stations deployed in northwestern Canada as part of the IRIS-PASSCAL/LITHOPROBE CANOE experiment. In this region radial-component receiver functions show a clear continental Moho and the presence of crustal material dipping into the mantle at the suture of two Palaeo-Proterozoic terranes. The geometry of the suture is better defined on the transverse component where subduction is associated with a ∼10 km thick layer exhibiting strong elastic anisotropy. The modelling reproduces the main features of the receiver functions, including the effects of anisotropy, heterogeneity and finite-frequency scattering.