A three-dimensional spherical mesh generator

A new spherical mesh generator is described. It represents an efficient, deterministic packing of tetrahedra into a solid sphere, a spherical shell, or both. The mesh can be used for finite-element solutions to a wide variety of global numerical modelling problems in the geosciences. The nodes within the mesh are distributed uniformly, and long, thin tetrahedra are avoided. The method proposed here offers several advantages over 3-D Delaunay algorithms for finite-element mesh generation. For the related problem of trivariate scattered data interpolation, which is not considered here, the 3-D Delaunay algorithms are the method of choice.     

The growing spherical seismic source

Summary. A solution is found for the seismic radiation from an arbitrarily growing spherical source in an inhomogeneously prestressed elastic medium. The general problem of the growing seismic source in a prestressed medium is formulated as a boundary value problem. For the special case of the growing spherical source, an expansion in vector spherical harmonics reduces the problem to a set of one-dimensional Volterra integral equations. The equations can be easily formed through the use of Bessel function recursion relations. The integral equations for a growing spherical cavity are solved numerically. Waveforms are then computed for homogeneous and inhomogeneous stress fields for several growth histories. The resulting waveforms are similar to the waveforms of the corresponding instantaneous problem, but stretched out in time and reduced in amplitude. The effects of diffraction and the overshoot of equilibrium are reduced with a reduction in growth rate. The effects caused by inhomogeneity of the stress field are quite strong for the growing as well as for the instantaneous seismic source.     

Bending of spherical lithosphere–axisymmetric case

Effects of sphericity are commonly ignored in the lithospheric bending problem. In order to examine its effects, I solve a simple axisymmetric spherical-shell model. The full solution and the asymptotic solution are derived from the basic equations, and their relationship to the flat-plate solution is examined. For displacement, effects of sphericity are small, and use of the flat-plate solution produces results that are numerically indistinguishable from those of the spherical solution. The most significant effect of sphericity appears in the stress, in particular the normal stress along the strike direction of the trench. This stress is approximately given by Eu_r/R , where E is Young's modulus, u_r is the vertical deformation of the shell and R is its radius of curvature. If the shell (lithosphere) is bent downwards and reaches 30 km, this stress can become about 5 kbar in the Earth. While plastic behaviour may set in under such high pressure conditions and analysis beyond elasticity theory may be required, sphericity may be a cause of large compressive stress in the trench strike direction. This stress may play an important role in forming the overall shape of the Earth's subduction zones.     

World population in a grid of spherical quadrilaterals

"We report on a project that converted subnational population data to a raster of cells on the earth. We note that studies using satellites as collection devices yield results indexed by latitude and longitude. Thus it makes sense to assemble the terrestrial arrangement of people in a compatible manner. This alternative is explored here, using latitude/longitude quadrilaterals as bins for population information.... The results to date of putting world boundary coordinates together with estimates of the number of people are described. The estimated 1994 population of 219 countries, subdivided into 19,032 polygons, has been assigned to over six million five minute by five minute quadrilaterals covering the world."     

Parametrizing surface wave tomographic models with harmonic spherical splines

We present a mathematical framework and a new methodology for the parametrization of surface wave phase-speed models, based on traveltime data. Our method is neither purely local, like block-based approaches, nor is it purely global, like those based on spherical harmonic basis functions. Rather, it combines the well-known theory and practical utility of the spherical harmonics with the spatial localization properties of spline basis functions. We derive the theoretical foundations for the application of harmonic spherical splines to surface wave tomography and summarize the results of numerous numerical tests illustrating the performance of a practical inversion scheme based upon them. Our presentation is based on the notion of reproducing-kernel Hilbert spaces, which lends itself to the parametrization of fully 3-D tomographic earth models that include body waves as well.     

Normal-mode representations of multiple-ray reflections in a spherical earth

Summary. In a spherically symmetric, isotropic earth model the duality between seismic rays and modes can be established completely by application of the principle of stationary phase to the summed normal-mode representation of the time signal. The requirement of stationary phase must be applied not only on the sum over the angular order but also over the radial-order summation.
This approach is illustrated by using asymptotic approximations to the equations for toroidal oscillations. In this way the travel-time formulae for rays from a surface source in a sufficiently smooth earth model are easily derived, and the distribution of modes between the different rays can be found. The result may be used as a selection criterion to reduce the number of modes that must be summed to construct synthetic seismograms at a certain distance and within a certain time window.     

Extension of the Thomson-Haskell method to non-homogeneous spherical layers

Summary. Amplitude spectra of Rayleigh and Love waves in a layered non-gravitating spherical earth have been obtained using as a source, displacement and stress discontinuities. In each layer elastic parameters and density follow specified functions of radial distance and the solutions of the equations of motion are obtained in terms of exponential functions. The Thomson—Haskell method is extended to this case. The problem reduces to simple calculations as in a plane-layered medium. Numerical results of phase and group velocities up to periods of 300 s in various earth models when compared with earlier results (obtained by numerical integration) show that the present method can be used with sufficient accuracy. The differences in phase velocity, group velocity and amplitude (also surface ellipticity in the case of Rayleigh waves) between spherical- and flat-earth models have been investigated in the range 20–300–s period and expressed in polynomials in the period.     

On the efficient calculation of ordinary and generalized spherical harmonics

Algorithms for the stable computation of generalized and ordinary spherical harmonics are presented. The algorithms are fast and have the useful property that they can compute harmonics for isolated harmonic degrees. fortran and C programs implementing these algorithms are available from the authors.     

Polarization and amplitude attributes of reflected plane and spherical waves

The characteristics of a reflected spherical wave at a free surface are investigated by numerical methods; in particular, the polarization angles and amplitude coefficients of a reflected spherical wave are studied. The classical case of the reflection of a plane P wave from a free surface is revisited in order to establish our terminology, and the classical results are recast in a way which is more suited for the study undertaken. The polarization angle of a plane P wave, for a given angle of incidence, is shown to be 90° minus twice the angle of reflection of the reflected S wave. For a Poisson's ratio less than 1/3, there is a non-normal incident angle for which both amplification coefficients are 2 precisely; for this incident angle the direction of the particle motion at the free surface is also the direction of the incident wave. For a wave emanating from a spherical source, the polarization angle, for all angles of incidence, is always less than, or equal to, the polarization angle of a plane P wave. The vector amplification coefficient of a spherical wave, for all angles of incidence, is always greater than the vector amplification coefficient of a plane P wave. As expected, the results for a spherical wave approach the results for a plane P wave in the far field. Furthermore, there was a good agreement between the theoretical modelling and the numerical modelling using the dynamic finite element method (DFEM).