Parker, R. L. | Search NASA

NASA NTRS1997

The Influence of Correlated Crustal Signals in Modelling the Main Geomagnetic Field

Algorithms used in geomagnetic main-field modelling have for the most part treated the noise in the field measurements as if it were white. A major component of the noise consists of the field due to magnetization in the crust and it has been realized for some time that such signals are highly correlated at satellite altitude. Hence approximation by white noise, while of undoubted utility, is of unknown validity. In this paper we study two plausible statistical models for the crustal magnetization, in which the magnetization is a realization of a stationary, isotropic, random process. At a typical satellite altitude the associated fields exhibit significant correlation over ranges as great as 15 deg. or more, which introduces off-diagonal elements into the covariance matrix, elements that have usually been neglected in modelling procedures. Dealing with a full covariance matrix for a large data set would present a formidable computational challenge, but fortunately most of the entries in the covariance matrix are so small that they can be replaced by zeros. The resultant matrix comprises only about 3 per cent non-zero entries and thus we can take advantage of efficient sparse matrix techniques to solve the numerical system. We construct several main-field models based on vertical-component data from a selected 5 deg. by 5 deg. data set derived from the Magsat mission. Models with and without off-diagonal terms are compared.

Rygaard-Hjalsted, C.↗

NASA NTRS1988

Statistics of the geomagnetic secular variation for the past 5 m.y

A new statistical model is proposed for the geomagnetic secular variation over the past 5Ma. Unlike previous models, the model makes use of statistical characteristics of the present day geomagnetic field. The spatial power spectrum of the non-dipole field is consistent with a white source near the core-mantle boundary with Gaussian distribution. After a suitable scaling, the spherical harmonic coefficients may be regarded as statistical samples from a single giant Gaussian process; this is the model of the non-dipole field. The model can be combined with an arbitrary statistical description of the dipole and probability density functions and cumulative distribution functions can be computed for declination and inclination that would be observed at any site on Earth's surface. Global paleomagnetic data spanning the past 5Ma are used to constrain the statistics of the dipole part of the field. A simple model is found to be consistent with the available data. An advantage of specifying the model in terms of the spherical harmonic coefficients is that it is a complete statistical description of the geomagnetic field, enabling us to test specific properties for a general description. Both intensity and directional data distributions may be tested to see if they satisfy the expected model distributions.

Constable, C. G.↗

NASA NTRS1986

Statistics of the geomagnetic secular variation for the past 5Ma

A new statistical model is proposed for the geomagnetic secular variation over the past 5Ma. Unlike previous models, the model makes use of statistical characteristics of the present day geomagnetic field. The spatial power spectrum of the non-dipole field is consistent with a white source near the core-mantle boundary with Gaussian distribution. After a suitable scaling, the spherical harmonic coefficients may be regarded as statistical samples from a single giant Gaussian process; this is the model of the non-dipole field. The model can be combined with an arbitrary statistical description of the dipole and probability density functions and cumulative distribution functions can be computed for declination and inclination that would be observed at any site on Earth's surface. Global paleomagnetic data spanning the past 5Ma are used to constrain the statistics of the dipole part of the field. A simple model is found to be consistent with the available data. An advantage of specifying the model in terms of the spherical harmonic coefficients is that it is a complete statistical description of the geomagnetic field, enabling us to test specific properties for a general description. Both intensity and directional data distributions may be tested to see if they satisfy the expected model distributions.

Constable, C. G.↗

NASA NTRS1985

A preliminary harmonic spline model from Magsat data

A preliminary main field model for 1980 derived from a carefully selected subset of Magsat vector measurements, using the method of harmonic splines, is presented. This model (PHS /80/) for preliminary harmonic splines is the smoothest model (in the sense that the rms radial field at the core surface is minimum) consistent with the measurements (with an rms misfit of 10 nT to account for crustal and external fields as well as noise in the measurement procedure). Therefore PHS (80) is more suitable for studies of the core than models derived with the traditional least squares approach (e.g., GSFC /9/80/). A comparison is conducted of the characteristics of the harmonic spline spectrum, topology of the core field and especially the null-flux curves (loci where the radial field is zero) and the flux through patches bounded by such curves. PHS (80) is less complex than GSFC (9/80) and is therefore more representative of that part of the core field that the data constrain.

Shure, L.↗

NASA NTRS1985

Separation of core and crustal magnetic field sources

Fluid motions in the electrically conducting core and magnetized crustal rocks are the two major sources of the magnetic field observed on or slightly above the Earth's surface. The exact separation of these two contributions is not possible without imposing a priori assumptions about the internal source distribution. Nonetheless models like these were developed for hundreds of years Gauss' method, least squares analysis with a truncated spherical harmonic expansion was the method of choice for more than 100 years although he did not address separation of core and crustal sources, but rather internal versus external ones. Using some arbitrary criterion for appropriate truncation level, we now extrapolate downward core field models through the (approximately) insulating mantle. Unfortunately our view can change dramatically depending on the degree of truncation for describing core sources.

Shure, L.↗

NASA NTRS1985

Does the geoid drift west?

In 1970 Hide and Malin noted a correlation of about 0.8 between the geoid and the geomagnetic potential at the Earth's surface when the latter is rotated eastward in longitude by about 160 degrees and the spherical harmonic expansions of both functions are truncated at degree 4. From a century of magnetic observatory data, Hide and Malin inferred an average magnetic westward drift rate of about 0.27 degrees/year. They attributed the magnetic-gravitational correlation to a core event at about 1350 A.D. which impressed the mantle's gravity pattern at long wavelengths onto the core motion and the resulting magnetic field. The impressed pattern was then carried westward 160 degrees by the nsuing magnetic westward drift. An alternative possibility is some sort of steady physical coupling between the magnetic and gravitational fields (perhaps migration of Hide's bumps on the core-mantle interface). This model predicts that the geoid will drift west at the magnetic rate. On a rigid earth, the resulting changes in sea level would be easily observed, but they could be masked by adjustment of the mantle if it has a shell with viscosity considerably less than 10 to the 21 poise. However, steady westward drift of the geoid also predicts secular changes in g, the local acceleration of gravity, at land stations. These changes are now ruled out by recent independent high-accuracy absolute measurements of g made by several workers at various locations in the Northern Hemisphere.

Backus, G. E.↗

NASA NTRS1974

Materials science and manufacturing in space: The NBS program of ground based research

National Bereau of Standards work for NASA in support of NASA's Space Processing Program is described. The objectives of the NBS program are to perform ground-based studies of those aspects of space that could possibly provide a unique environment for making materials more perfect or more pure. The approach taken deals primarily with experimental and theoretical studies of the possible effects of the absence of gravitational forces on those materials preparation processes where the presence of gravity may be important in reducing perfection or purity. The materials preparation processes studied comprise 6 tasks in the areas of crystal growth, purification and chemical processing, and the preparation of composites. They are: (1) crystal perfection in Czochralski growth; (2) evaporative preparation of ultrahigh purity materials; (3) vacuum effects in the preparation of composite materials; (4) melt shape in weightless crystal growth; (5) vapor transport synthesis and crystal growth of oxides; and (6) surface traction and other surface phenomena.

Passaglia, E.↗

NASA NTRS1972

The Backup-Gilbert method and its application to the electrical conductivity problem

The theory of Backus and Gilbert gives a technique for solving the general linear inverse problem. Observational error and lack of data are shown to reduce the reliability of the solution in different ways: the former introduces statistical uncertainties in the model, while the latter smooths out the detail. Precision can be improved by sacrificing resolving power, and vice versa, so that some compromise may be made between the two in choosing the best model. Nonlinear inverse problems can be brought into the domain of the theory by linearizing about a typical solution. The inverse problem of electrical conductivity in the mantle is used to illustrate the Backus-Gilbert technique; an example of the tradeoff diagram is given.

Parker, R. L.↗