Spherical harmonic representation of the gravitational potential of discrete spherical mass elements

b
We present expressions in a spherical harmonic framework for the gravitational potential of discrete point, surface, and volume mass elements located at any depth within a sphere. Through analysis of the spherical harmonic spectrum, insight is gained into the properties of the potentials arising from a variety of mass distributions. A point mass at the surface of a sphere displays the richest harmonic spectrum in all degrees; spectra become increasingly reddened as the source mass is distributed through larger elements of area or volume, or is located at greater depths below the surface of the reference sphere. The spectra of dipolar distributions, useful in representing compensated masses, are depressed, especially in the low harmonic degrees, relative to the spectra of monopole elements.     

Assessment of regional geomagnetic field modelling methods using a standard data set: spherical cap harmonic analysis

Various methods that take account of the potential nature of the field have been proposed for modelling geomagnetic data on a regional scale. Several of these have been applied to a standard data set based on annual mean values from observatories in Europe. Here, we examine some of the properties of spherical cap harmonic analysis when applied to this data set, and compare the quality of fit with that of the other models. It is found that, for this data set, rectangular polynomial analysis provides a compact fit to main field data, but that in most other cases, for both main field and anomaly data, spherical cap harmonic analysis provides the better fit. Although relatively insensitive to chosen cap size, spherical cap harmonic analysis deteriorates more rapidly than the other methods when the number of coefficients is reduced.     

Anomalous Love-wave phase velocities in the Pacific: sequential pure-path and spherical harmonic inversion

Summary. We have determined the lateral distribution of Love-wave phase velocities in the Pacific for the periods 40, 67, 91 and 125 s. Application of the pure-path and spherical harmonic representation methods indicates that the velocities are primarily a function of the age of the seafloor. A comparison of the results from these two techniques indicates inherent modelling constraints in both methods. The pure-path method is limited by its a priori nature while the spherical harmonic approach is unsuitable in describing sharp lateral velocity gradients. To circumvent these limitations, we propose the sequential application of the pure-path and spherical harmonic methods. The sequential inversion separates the velocity distribution into two separate components; velocity as a function of the age of the oceanic plate and variations superimposed on this relationship. Application of this method demonstrates the presence of velocity anomalies which cannot be modelled by an age–velocity relationship. These anomalies are tentatively correlated with regions of anomalous seafloor depths and/or the presence of active hot-spots. In the central south Pacific, an area with numerous active hot-spots coincides roughly with a region of anomalously slow Love wave velocities. A method for determining the errors associated with the slowness distributions calculated by the spherical harmonic method is presented and provides a means for determining the resolvability of these features.     

极区电离层TEC经验模型的建立及适用性分析

利用南极地区40多个GPS跟踪站2010年全年的实测数据,实现了极区电离层TEC建模,对多项式模型、广义三角级数函数模型、低阶球谐函数模型、改进的球谐函数模型以及球冠谐函数模型等五种电离层经验模型进行了比较,并评估了其在极区的适用性情况。结果表明,各个模型在极区都可以取得比较好的拟合精度,残差均值在0.1TECU以内,均方根误差在2 TECU以内。     

Rayleigh-wave dispersion function for a transversely isotropic layered spherical earth using the Thomson-Haskell matrix method

We consider the two coupled differential equations of the two radial functions appearing in the displacement components of spheroidal oscillations for a transversely isotropic (TI) medium in spherical coordinates. Elements of the layer matrix have been explicitly written—perhaps for the first time—to extend the use of the Thomson-Haskell matrix method to the derivation of the dispersion function of Rayleigh waves in a transversely isotropic spherical layered earth. Furthermore, an earth-flattening transformation (EFT) is found and effectively used for spheroidal oscillations. The exponential function solutions obtained for each layer give the dispersion function for TI spherical media the same form as that on a flat earth. This has been achieved by assuming that the five elastic parameters involved vary as r ^p and that the density varies as r ^p-2, where p is an arbitrary constant and r is the radial distance. A numerical illustration with p = - 2 shows that, in spite of the inhomogeneity assumed within layers, the results for spherical harmonic degree n , versus time period T , obtained here for the Primary Reference Earth Model (PREM), agree well with those obtained earlier by other authors using numerical integration or variational methods. The results for isotropic media derived here are also in agreement with previous results. The effect of transverse isotropy on phase velocity for the first two modes of Rayleigh waves in the period range 20 to 240 s is calculated and discussed for continental and oceanic models.     

Spectral analysis using orthonormal functions with a case study on the sea surface topography

b
Spherical harmonics are orthonormalized using the Gram-Schmidt process in a function space. The problem of linear dependence of spherical harmonics over the oceans is studied using the Gram matrices and consequently three sets of orthonormal (ON) functions have been constructed. For the process an efficient formula for computing inner products of spherical harmonics has been developed. Important spectral properties of the ON functions are addressed. The ON functions may be used for representing the sea surface topography (SST) in the analysis of satellite altimeter data. The geoid error can be transformed to a representation by the ON functions and hence the comparison of powers of the geoid error and the SST signal only over the oceans is possible, leading to a better way of determining the cut-off frequency of the SST in the simultaneous solution using satellite altimeter data. As a case study, the modified Levitus SST is expanded into the ON functions. The results show that 99.90 per cent of that signal's energy is contained within degree 24 of the orthonormal functions. Such expansions also render better spectral behaviour of oceanic signals as compared to that from spherical harmonic expansions. The study shows that these generalized Fourier functions are suitable for spectral analyses of oceanic signals and they can be applied to future altimetric mission where the geoid and the SST are to be recovered.     

The geomagnetic field over the past 5 million years

The modern geomagnetic field is usually expressed as a spherical harmonic expansion. Although the palaeomagnetic record is very incomplete in both space and time, sufficient data are available from a span of ages to generate time-averaged spherical harmonic field models with many degrees of freedom. Here three data sets are considered: directional measurements from lavas, inclination measurements from ocean sediments, and intensity measurements from lavas. Individual data are analysed, as well as site-averages, using the same methods that have been developed for the modern field, to give models for the past 5 Myr. The normal-polarity field model has an axial-dipole intensity similar to that of the modern-day field, whilst the equatorial-dipole component is very much smaller. The field is not axisymmetric, but shows flux concentrations at the core's surface under Canada and Siberia similar to those observed in the field over historical timescales. Tests on synthetic data show that it is unlikely that these similarities result from the overprinting of the palaeomagnetic field due to inadequate cleaning of the samples. The reverse-polarity field model does not show such obvious features, but this may be due to the sparsity of the data.
The patterns observed in the normal-polarity field, with persistent features in the northern hemisphere and a smooth southern hemisphere, could be explained if the present pattern of secular variation is typical of the past several million years. This would reveal itself as large variations over time in the direction of the magnetic vector in regions of high secular variation, with relatively little change over quieter regions. However, we have been unable to find any evidence for a geographical pattern of secular variation in the data.     

An improved geomagnetic data selection algorithm for global geomagnetic field modelling

A spherical harmonic degrees 60, global internal field model is described (called BGS/G/L/0706). This model includes a degree 15 core and piecewise-linear secular variation model and is derived from quiet-time Ørsted and Champ satellite data sampled between 2001.0 and 2005.0. For the satellite data selection, a wide range of geomagnetic index and other data selection filters have been used to best isolate suitably quiet magnetospheric and ionospheric conditions. Only a relatively simple, degree one spherical harmonic, external field model is then required. It is found that a new 'Vector Magnetic Disturbance' index ( VMD ), the existing longitude sector A indices, the auroral zone index IE , and the polar cap index PC are better than Kp and Dst at rejecting rapidly varying external field signals at low, middle, auroral and polar latitudes. The model quality is further enhanced by filling spatial and temporal gaps in the quiet data selection with a second selection containing slightly more disturbed data. It is shown that VMD provides a better parametrization than Dst of the large-scale, rapidly changing, external field. The lithospheric field model between degrees 16 and 50 is robust and displays good coherence with other recently published models for this epoch. BGS/G/L/0706 also shows crustal anomalies consistent with other studies, although agreement is poorer in the southern polar cap. Intermodel coherency reduces above about degree 40, most likely due to incompletely filtered signals from polar ionospheric currents and auroral field aligned currents. The absence of the PC index for the southern hemisphere for 2003 onwards is a particular concern.     

Geomagnetic effects of the Earth's ellipticity

We analyse the external field generated by a uniform distribution of magnetic susceptibility contained in an oblate spheroidal shell when it is magnetized by an internal magnetic field of arbitrary complexity. The situation is more relevant to the Earth than that of a spherical shell considered by Runcorn (1975a ) (in the context of lunar magnetism), because of the larger flattening of the Earth than that of the Moon. We find that, to first order in the susceptibility, each internal harmonic in a spheroidal harmonic expansion of the magnetic potential generates just one non-vanishing external field coefficient, unlike in the spherical case when all harmonics vanish identically. The field generated is proportional to the susceptibility, thickness of the shell and square of the Earth's eccentricity, and hence it appears that this field amplification mechanism will be very ineffective for the Earth.     

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.     

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.     

The geomagnetic power spectrum

Combining CHAMP satellite magnetic measurements with aeromagnetic and marine magnetic data, the global geomagnetic field has now been modelled to spherical harmonic degree 720. An important tool in field modelling is the geomagnetic power spectrum. It allows the comparison of field models estimated from different data sets and can be used to identify noise levels and systematic errors. A correctly defined geomagnetic power spectrum is flat (white) for an uncorrelated field, such as the Earth's crustal magnetic field at long wavelengths. It can be inferred from global spherical harmonic models as well as from regional grids. Marine and aeromagnetic grids usually represent the anomaly of the total intensity of the magnetic field. Appropriate corrections have to be applied in estimating the geomagnetic power spectrum from such data. The comparison of global and regional spectra using a consistently defined azimuthally averaged geomagnetic power spectrum facilitates quality control in field modelling and should provide new insights in magnetic anomaly interpretation.     

THE EARTH'S MAIN MAGNETIC FIELD—EPOCH 1955.0

Summary. A brief account of the material employed in the preparation of the Admiralty world charts for epoch 1955^.0 is given together with the results of a spherical harmonic analysis made to assist in estimating the field values in poorly observed areas. The deduced Gaussian coefficients are tabulated.     

Taking into account truncation problems and geomagnetic model accuracy in assessing computed flows at the core-mantle boundary

SUMMARY
Since the time Roberts & Scott (1965) first expressed the key 'frozen flux' hypothesis relating the secular variation of the geomagnetic field (SV) to the flow at the core surface, a large number of studies have been devoted to building maps of the flow and inferring its fundamental properties from magnetic observations at the Earth's surface. There are some well-known difficulties in carrying out these studies, such as the one linked to the non-uniqueness of the flow solution [if no additional constraint is imposed on the flow (Backus 1968)] which has been thoroughly investigated. In contrast little investigation has been made up to now to estimate the exact importance of other difficulties, although the different authors are usually well aware of their existence. In this paper we intend to make as systematic as possible a study of the limitations linked to the use of truncated spherical harmonic expansions in the computation of the flow. Our approach does not rely on other assumptions than the frozen flux, the insulating mantle and the large-scale flow assumptions along with some simple statistical assumptions concerning the flow and the Main Field. Our conclusions therefore apply to any (toroidal, steady or tangentially geostrophic) of the flow models that have already been produced; they can be summarized in the following way: first, because of the unavoidable truncation of the spherical harmonic expansion of the Main Field to degree 13, no information will ever be derived for the components of the flow with degree larger than 12; second, one may truncate the spherical harmonic expansion of the flow to degree 12 with only a small impact on the first degrees of the flow. Third, with the data available at the present day, the components of the flow with degree less than 5 are fairly well known whereas those with degree greater than 8 are absolutely unconstrained.     

Lunar and solar tidal components in the occurrence of earthquakes in Italy

Summary. A method for the determination of lunar and solar daily components of earthquakes is proposed. Its application to the most intense earthquakes with epicentres in Italy (from 1900 to 1983) provides for the first time large and significant lunar and solar tidal terms. The absence of significant results from the analysis of 9157 world-wide large shocks processed according to the proposed method and to routine spherical harmonic techniques indicates that the Italian earthquakes are locally tidally triggered.     

The Sq current system during the International Geophysical Year

Summary. The daily variation for the mean of the eight most quiet days of the IGY is represented as a function of latitude, longitude and UT by means of a spherical harmonic model. The most appropriate choice of Sq baseline is investigated, on the assumption that the external electric currents deduced from the model should be negligible at local midnight. The properties of the final model are discussed, particularly with regard to the external current foci, and are compared with those of previous models.     

Spherical Slepian functions and the polar gap in geodesy

The estimation of potential fields such as the gravitational or magnetic potential at the surface of a spherical planet from noisy observations taken at an altitude over an incomplete portion of the globe is a classic example of an ill-posed inverse problem. We show that this potential-field estimation problem has deep-seated connections to Slepian's spatiospectral localization problem which seeks bandlimited spherical functions whose energy is optimally concentrated in some closed portion of the unit sphere. This allows us to formulate an alternative solution to the traditional damped least-squares spherical harmonic approach in geodesy, whereby the source field is now expanded in a truncated Slepian function basis set. We discuss the relative performance of both methods with regard to standard statistical measures such as bias, variance and mean squared error, and pay special attention to the algorithmic efficiency of computing the Slepian functions on the region complementary to the axisymmetric polar gap characteristic of satellite surveys. The ease, speed, and accuracy of our method make the use of spherical Slepian functions in earth and planetary geodesy practical.     

Gravitational and magnetic fields of some simple solids of revolution

Summary. Exact spherical harmonic expansions are given for calculating the gravitational and magnetic fields associated with certain uniform solids of revolution. The figures are those made by rotating a conic section about one of its principal axes. The coefficients in the expansions can be computed accurately and efficiently and this approach leads to a very satisfactory method for calculating the fields of geological bodies with approximate circular symmetry about a vertical axis. A complete theory of convergence is given for the expansions. Somewhat unexpectedly, the sphere of convergence is determined by the location of a number of equivalent point or line sources that lie within the body or on its edges.     

Geomagnetic secular acceleration

Summary. A spherical harmonic model of the second time-derivative of the geomagnetic field is determined, for the first time, directly from measures of the secular acceleration based on observatory annual mean data. The data span the interval 1964.5–1975.5, and 165 observatories are included. The model comprises the 32 coefficients of degree and order up to 6 that are significant at the 5 per cent level. Its primary purpose is to aid in the reduction of data to epoch for the 1980 series of navigational charts. The model is compared with earlier estimates of secular acceleration, derived by less direct methods.     

Numerical aspects of a spline-based multiresolution recovery of the harmonic mass density out of gravity functionals

We show the numerical applicability of a multiresolution method based on harmonic splines on the 3-D ball which allows the regularized recovery of the harmonic part of the Earth's mass density distribution out of different types of gravity data, for example, different radial derivatives of the potential, at various positions which need not be located on a common sphere. This approximated harmonic density can be combined with its orthogonal anharmonic complement, for example, determined out of the splitting function of free oscillations, to an approximation of the whole mass density function. The applicability of the presented tool is demonstrated by several test calculations based on simulated gravity values derived from EGM96. The method yields a multiresolution in the sense that the localization of the constructed spline basis functions can be increased which yields in combination with more data a higher resolution of the resulting spline. Moreover, we show that a locally improved data situation allows a highly resolved recovery in this particular area in combination with a coarse approximation elsewhere which is an essential advantage of this method, for example, compared to polynomial approximation.