Low-frequency and high-frequency expressions for the reflection and transmission coefficients of seismic waves for transition layers

Summary Low-frequency and high-frequency expressions are derived for the reflection and transmission coefficients of SH waves on transition layers, in which the modulus of torsion and density vary with depth only. Both the expressions are exact for any finite frequencies. However, the use of the former formulae is more convenient for low frequencies and, likewise, the latter formulae for high frequencies. Explicit expressions for several first terms of the expansions in terms of and 1/ are given.     

Traveling planetary waves in the stratosphere

This paper reviews the theory and observations of some traveling planetary waves in the stratosphere. Two categories of waves which appear prominently in the literature are discussed: westward propagating waves of periods in the range 3–7 days (the 5-day wave) and in the range 10–20 days (the 16-day wave). Although the observations seem to indicate that these waves are waves of the Rossby type (planetary waves), the evidence is less clear regarding (1) the question of whether these waves are forced internal waves or free (resonant) external waves, and (2) the identification of the observed waves with specific theoretical waves of the Rossby type. When recent observations are compared with theory, the evidence seems to favor the notion that the 5-day and 16-day waves of longitudinal wave number 1 may be identified, respectively, with the gravest and next gravest symmetric free Rossby modes. However, the observational evidence seems to be less clear regarding the nature of the 16-day wave than the 5-day wave.     

Computation of travel times and rays in a medium of two-dimensional velocity distribution

A ray-tracing computer program is described for a two-dimensional velocity distribution defined by the local wave velocitiesv _{i, j} in points at arbitrary depthv _{i, j} below the surface points with the horizontal coordinatesx _i . The velocity variation is assumed to be linear in the triangles formed by three neighbouring points. Travel times and rays are then calculated after the exact analytic formulae for any position of the source within the model. No assumptions other than of a piecewise linear velocity structure are made. A first-order discontinuity can be approximated by a thin layer with a strong velocity gradient and refracted waves or wide-angle reflections obtained in this way. As an example,P-wave rays were computed for section No. 05 of the Alpine Longitudinal Profile. The model includes a low-velocity channel which is cut off on the eastern side, first-order discontinuities and a sediment basin.Paper presented at the ESC-Workshop Meeting Seismic Waves in Laterally Inhomogeneous Media, Liblice, 1978.     

On wave motions in a tropical spectral experiment

An attempt is made to use the barotropic vorticity equation in spectral form in order to study barotropic instability when the basic current has east-west quasi-stationary asymmetries on the scale of long waves. This is done by expressing the spectral equations in three different ways. In the first experiment a 9-component system is integrated and the long waves are allowed to propagate freely. In the second experiment the long waves are constrained to propagate slowly and in the third experiment, they are removed altogether.The motivation behind the present investigation is due to observations made in the motion fields of the tropics. These are characterized by quasi-stationary long waves and very energetically active and propagating short waves.The presence of quasi-stationary long waves seems to enhance the energy exchanges between the short waves and the mean zonal current and also allow for larger energy values for the short waves. Long term integration (90 days) shows a 6-day mode in thev time spectra at latitude 10°N and a 15-day mode in theu time spectra at individual grid points in the equatorial latitudes. Any possible connection between this peak and the observed peak of 15 days in the completely different physics of Kelvin waves is left as a conjecture.     

Low frequency waves in plasmas with spatially varying electron temperature

Low frequency electrostatic waves are studied in magnetized plasmas with an electron temperature which varies with position in a direction perpendicular to the magnetic field. For wave frequencies below the ion cyclotron frequency, the waves need not follow any definite dispersion relation. Instead a band of phase velocities is allowed, with a range of variation depending on the maximum and minimum values of the electron temperature. Simple model equations are obtained for the general case which can be solved to give the spatial variation of a harmonically time varying potential. A simple analytical model for the phenomenon is presented and the results are supported by numerical simulations carried out in a 2.5-dimensional particle-in-cell numerical simulation. We find that when the electron temperature is striated along B₀ and low frequency waves (_ci) are excited in this environment, then the intensity of these low frequency waves will be striated in a manner following the electron temperature striations. High frequency ion acoustic waves (_ci) will on the other hand have a spatially more uniform intensity distribution.     

Identification of low-frequency kinetic wave modes in the Earth’s ion foreshock

In this work we use ion and magnetic field data from the AMPTE-UKS mission to study the characteristics of low frequency (_{r p}) waves observed upstream of the Earths bow shock. We test the application of various plasma-field correlations and magnetic ratios derived from linear Vlasov theory to identify the modes in this region. We evaluate (for a parameter space consistent with the ion foreshock) the Alfvén ratio, the parallel compressibility, the crosshelicity, the noncoplanar ratio, the magnetic compression and the polarization for the two kinetic instabilities that can be generated in the foreshock by the interaction of hot diffuse ions with the solar wind: the left-hand resonant and the right-hand resonant ion beam instabilities. Comparison of these quantities with the observed plasma-field correlations and various magnetic properties of the waves observed during 10 intervals on 30 October 1984, where the waves are associated with diffuse ions, allows us to identify regions with Alfvénic waves and regions where the predominant mode is the right-hand resonant instability. In all the cases the waves are transverse, propagating at angles 33° and are elliptically polarized. Our results suggest that while the observed Alfvén waves are generated locally by hot diffuse ions, the right-handed waves may result from the superposition of waves generated by two different types of beam distribution (i.e. cold beam and diffuse ions). Even when there was good agreement between the values of observed transport ratios and the values given by the theory, some discrepancies were found. This shows that the observed waves are different from the theoretical modes and that mode identification based only on polarization quantities does not give a complete picture of the waves characteristics and can lead to mode identification of waves whose polarization may agree with theoretical predictions even when other properties can diverge from those of the theoretical modes.     

Approximate expression for the hilbert transform of a certain class of functions and their application to the ray theory of seismic waves

Summary Approximate expressions for the Hilbert transform of the functionf(t)=exp(- ₀ ² t ²/²) cos( ₀ t+v) are determined. This function, given a suitable choice of the three parameters ₀, and v, approximates a wide class of seismic signals very well. The approximate expressions for the Hilbert transform enable very simple formulae to be given for the elementary seismograms of the individual seismic body waves (in the zero approximation of the ray theory). This accelerates the computation of ray theoretical seismograms considerably.     

Scattering ofP andS waves by a spherically symmetric inclusion

Scattering of an arbitrary elastic wave incident upon a spherically symmetric inclusion is considered and solutions are developed in terms of the spherical vector system of Petrashen, which produces results in terms of displacements rather than displacement potentials and in a form suitable for accurate numerical computations. Analytical expressions for canonical scattering coefficients are obtained for both the cases of incidentP waves and incidentS waves. Calculations of energy flux in the scattered waves lead to elastic optical theorems for bothP andS waves, which relate the scattering cross sections to the amplitude of the scattered fields in the forward direction. The properties of the solutions for a homogeneous elastic sphere, a sphere filled by fluid, and a spherical cavity are illustrated with scattering cross sections that demonstrate important differences between these types of obstacles. A general result is that the frequency dependence of the scattering is defined by the wavelength of the scattered wave rather than the wavelength of the incident wave. This is consistent with the finding that the intensity of thePS scattering is generally much stronger than theSP scattering. When averaged over all scattering angles, the mean intensity of thePS converted waves is2V _p ² /V _s ⁴ times the mean intensity of theSP converted waves, and this ratio is independent of frequency. The exact solutions reduce to simple and easily used expressions in the case of the low frequency (Rayleigh) approximation and the low contrast (Rayleigh-Born) approximation. The case of energy absorbing inclusions can also be obtained by assigning complex values to the elastic parameters, which leads to the result that an increase in attenuation within the inclusion causes an increased scattering cross section with a marked preference for scatteredS waves. The complete generality of the results is demonstrated by showing waves scattered by the earth's core in the time domain, an example of high-frequency scattering that reveals a very complex relationship between geometrical arrivals and diffracted waves.     

Gravity Anomalies of Arbitrary 3D Polyhedral Bodies with Horizontal and Vertical Mass Contrasts

During the last 15 years, more attention has been paid to derive analytic formulae for the gravitational potential and field of polyhedral mass bodies with complicated polynomial density contrasts, because such formulae can be more suitable to approximate the true mass density variations of the earth (e.g., sedimentary basins and bedrock topography) than methods that use finer volume discretization and constant density contrasts. In this study, we derive analytic formulae for gravity anomalies of arbitrary polyhedral bodies with complicated polynomial density contrasts in 3D space. The anomalous mass density is allowed to vary in both horizontal and vertical directions in a polynomial form of \(\lambda =ax^m+by^n+cz^t\), where m, n, t are nonnegative integers and a, b, c are coefficients of mass density. First, the singular volume integrals of the gravity anomalies are transformed to regular or weakly singular surface integrals over each polygon of the polyhedral body. Then, in terms of the derived singularity-free analytic formulae of these surface integrals, singularity-free analytic formulae for gravity anomalies of arbitrary polyhedral bodies with horizontal and vertical polynomial density contrasts are obtained. For an arbitrary polyhedron, we successfully derived analytic formulae of the gravity potential and the gravity field in the case of \(m\le 1\), \(n\le 1\), \(t\le 1\), and an analytic formula of the gravity potential in the case of \(m=n=t=2\). For a rectangular prism, we derive an analytic formula of the gravity potential for \(m\le 3\), \(n\le 3\) and \(t\le 3\) and closed forms of the gravity field are presented for \(m\le 1\), \(n\le 1\) and \(t\le 4\). Besides generalizing previously published closed-form solutions for cases of constant and linear mass density contrasts to higher polynomial order, to our best knowledge, this is the first time that closed-form solutions are presented for the gravitational potential of a general polyhedral body with quadratic density contrast in all spatial directions and for the vertical gravitational field of a prismatic body with quartic density contrast along the vertical direction. To verify our new analytic formulae, a prismatic model with depth-dependent polynomial density contrast and a polyhedral body in the form of a triangular prism with constant contrast are tested. Excellent agreements between results of published analytic formulae and our results are achieved. Our new analytic formulae are useful tools to compute gravity anomalies of complicated mass density contrasts in the earth, when the observation sites are close to the surface or within mass bodies.     

Contributions to the spectral theory of seismic sources and determination of the parameters of earthquake foci

A spectral theory of seismic sources is developed which treats earthquake foci as band-pass filters. This concept is applied to focal mechanisms of the type double-couple with moment, and the obtained spectral function of the source is compared with those obtained by other authors. Finally some formulae for the determination of geometric and dynamic source parameters from the bandpass parameters are presented.     

Boundary integral solutions of geoelectric potential problems for the perturbing body of extreme conductivity

Summary Theoretical formulae for boundary integral solutions of the geoelectric potential fields for a perfectly conducting perturbing 3D body, as well as for a perfectly insulating body are presented. It is shown that the boundary integral equation for the double-layer density, can be solved for all physically possible conductivities of the perturbing body.Dedicated to the Memory of K. P     

Non-linear dynamo waves in an incompressible medium when the turbulence dissipative coefficients depend on temperature

Non-linear - dynamo waves existing in an incompressible medium with the turbulence dissipative coefficients depending on temperature are studied in this paper. We investigate of - solar non-linear dynamo waves when only the first harmonics of magnetic induction components are included. If we ignore the second harmonics in the non-linear equation, the turbulent magnetic diffusion coefficient increases together with the temperature, the coefficient of turbulent viscosity decreases, and for an interval of time the value of dynamo number is greater than 1. In these conditions a stationary solution of the non-linear equation for the dynamo waves amplitude exists; meaning that the magnetic field is sufficiently excited. The amplitude of the dynamo waves oscillates and becomes stationary. Using these results we can explain the existence of Maunders minimum.     

Analysis of various large-scale disturbances and the associated zonal mean motions in the atmosphere

Summary In this article, we present a scale analysis of planetary waves, extended long waves, and long waves. (We mean the extended long waves to be the disturbances whose east-west length is of order 10⁶ m and north-south extension 10⁷ m). We find for the extended long waves the two terms, the interaction between kinetic and available potential energy of the disturbances, and the interaction between the zonal mean available potential energy, and the eddy available potential energy, are of two orders of magnitude larger than the kinetic energy interaction between the disturbances and the associated zonal mean flow. This theoretical result concerning the relative importance of the various interaction terms may be of use in explaining the observational findings thus far available.It is also shown theoretically that the kinetic energy interaction between the planetary waves, the horizontal size of which is 10⁷ m, and the long waves, whose horizontal size is 10⁶ m, is of the same order as the interaction of kinetic energy between the zonal mean motion and the disturbances. This agrees fairly well with the observational estimates thus far obtained.     

A note on the general concept of wave breaking for Rossby and gravity waves

A recently proposed general definition of wave breaking is further discussed, in order to deal with some points on which misunderstanding appears to have arisen. As with surface and internal gravity waves, the classification of Rossby waves into breaking and not breaking is a generic classification based on dynamical considerations, and not a statement about any unique signature or automatically recognizable shape. Nor is it a statement about passive tracers uncorrelated with potential vorticity on isentropic surfaces. A strong motivation for the definition is that proofs of the nonacceleration theorem of wave, mean-flow interaction theory rely, explicitly or implicitly, on a hypothesis that the waves do not break in the sense envisaged.The general definition refers to the qualitative behaviour of a certain set of material contours, namely those, and only those, which would undulate reversibly, with small slopes, under the influence of the waves' restoring mechanism, in those circumstances for which linearized, nondissipative wave theory is a self-consistent approximation to nonlinear reality. The waves' restoring mechanism depends upon the basic-state vertical potential density gradient in the case of gravity waves, and upon the basic-state isentropic gradient of potential vorticity in the case of Rossby waves. In the usual linearized theory of planetary scale Rossby waves on a zonal shear flow, the relevant material contours lie along latitude circles when undisturbed.     

Wave propagation in a uniformly rotating fluid

Summary An unsteady flow generated by a harmonically oscillating pressure distribution of frequency acting on the paraboloidal free surface of an inviscid, incompressible fluid rotating with uniform angular velocity has been investigated. It is shown that case (i), >2 , corresponds to the usual surface waves, and case (ii), <2 , in contrast to the surface waves, corresponds to the inertial waves which are originated entirely due to rotation and have no counterpart in a non-rotating fluid motion. An explicit solution of the problem related to the above cases are obtained by the joint Laplace and Hankel transforms treatment in conjunction with asymptotic methods. The significant effects of the Coriolis force and the curvature of the free surface on the wave motions have been investigated. A comparison is made between the solutions of the problems with the horizontal and the paraboloidal free surface curvature. The analysis is concluded by exihibiting the characteristic features of the wave motions.     

M m: Extension to Love waves of the concept of a variable-period mantle magnitude

We extend to Love waves the concept of the mantle magnitudeM _mintroduced recently for Rayleigh waves. Spectral amplitudesX() of Love waves in the 50–300 s period range are measured on broad-band records from major events. A distance correctionC _D, regionalized to reflect the influence of different tectonic paths, and a source correctionC _S, compensating for the variation of excitation with period are effected; the exact geometry and depth of the event are however ignored. The resulting expression

Falling sphere observations of anisotropic gravity wave motions in the upper stratosphere over Australia

A study of inertial scale gravity wave motions in the region of the atmosphere between 30 and 60 km has been undertaken, using wind and temperature data derived from rocket-borne falling sphere density experiments performed over Woomera, Australia between 1962 nad 1976. The gross features of the wave field compare favorably with those found in similar northern hemispheric studies. Wave propagation is found to be both vertically and horizontally anisotropic. A rotary spectral analysis indicates predominately upgoing wave energy, suggesting that the majority of sources of these waves lie below 30 km. A detailed statistical investigation of the waves, made using the Stokes parameters technique, reveals that phase progression is also highly directional in the horizontal, with a significant zonal component in summer, but with a strong meridional component in winter. Propagation towards the southeast is inferred in summer, with the waves possibly emanating from tropospheric sources in equatorial regions to the north of Australia. The technique also shows that, on average, the waves appear to have mean ellipse eccentricities (=f/) around 0.4–0.45. Indirect estimates of a number of important wave parameters are made. In particular,v andw flux estimates are made over several height intervals. The vertical gradient of density weighted flux implies wave-induced mean flow accelerations of the order 0.1–1 ms^–1day^–1. This suggests that dissipating gravity waves are a significant source of the momentum residuals that are encountered in studies of satellite data from this region.     

Reflection of waves from the boundary of a random elastic semi-infinite medium

In this paper the smooth perturbation technique is employed to investigate the problem of reflection of waves incident on the plane boundary of a semi-infinite elastic medium with randomly varying inhomogeneities. Amplitude ratios have been obtained for various types of incident and reflected waves. It has been shown that an incidentSH orSV type of wave gives rise to reflectedSH, P andSV waves, the main components beingSH andP, SV in the respective cases. The reflected amplitudes have been calculated depending upon the randomness of the medium to the square of the small quantity , where measures the deviation of the medium from homogeneity. An incidentP-type wave produces mainly aP component and also a weakSH component to the order of ². The reflected amplitudes obtainable for elastic media are also altered by terms of the same order. The direction of the reflected wave is influenced by randomness in some cases.