Geomagnetic field analysis—IV. Testing the frozen-flux hypothesis

Summary. We present a model of the magnetic field at the core–mantle boundary, for epoch 1959.5, based on a large set of observatory and survey measurements. Formal error estimates for the radial field at the core are 50 μT, compared with 30 and 40 μT for our previous MAGSAT (1980) and POGO (1970) models.
Current work on the determination of the velocity of the core fluid relies on the assumption that the core behaves as a perfect conductor, so that the field lines remain frozen to the fluid at the core surface. This frozen-flux condition requires that the integrated flux over patches of the core surface bounded by contours of zero radial field remain constant in time. A new method is presented for constructing core fields that satisfy these frozen-flux constraints. The constraints are non-linear when applied to main field data, unlike the case of secular variation which was considered in an earlier paper. The method is applied to datasets from epochs 1969.5 and 1959.5 to produce fields with the same flux integrals as the 1980 model.
The frozen-flux hypothesis is tested by comparing the changes in the flux integrals between 1980/1969.5, 1969.5/1959.5 and 1980/1959.5 with their errors. We find that the hypothesis can be rejected with 95 per cent confidence. The main evidence for flux diffusion is in the South Atlantic region, where a new null flux curve appears between 1960 and 1970, and continues to grow at a rapid rate from 1970 to 1980. However, the statistical result depends critically on our error estimates for the field at the core surface, which are difficult to assess with any certainty; indeed, doubling the error estimates negates the statistical argument. The conclusion is therefore, at this stage, tentative, and requires further evidence, either from older data, if good enough, or from future satellite measurements.     

Geomagnetic field analysis—III. Magnetic fields on the core—mantle boundary

Summary. The method of stochastic inversion, previously applied to secular variation data, is applied to main field data. Adaptations to the method are required: non-linear, as well as linear, data are used; allowance is made for crustal components in the observatory data; and the prior information is specified differently. The requirement that the models should satisfy a finite lower bound on the Ohmic heating in the core provides strong prior information and gives finite error estimates at the core—mantle boundary.
The new method is applied to data from the epochs 1969.5 and 1980.0. The resulting field models are very much more complex than other models, such as the IGRF models extrapolated to the core, and show considerable small-scale detail which, on the basis of the error analysis, can be believed.
The flux integral over the northern hemisphere is computed at each epoch; the difference between the two epochs is approximately one standard deviation, suggesting that the question as to whether the decay of the dipole is consistent with the frozen-flux hypothesis has been resolved in favour of the hypothesis.     

Geomagnetic lunar analysis by least-squares

Summary. A simple and accurate method is presented for the determination of lunar and solar harmonic terms present in a series of geomagnetic data. It is based on least-squares and incorporates a direct means for determining the confidence limits. The method is tested by applying it to a series of artificially generated data, and comparing the results both with the true answers and with results obtained using existing methods.     

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.     

Geomagnetic excursion in the late Cretaceous

Summary. Two sedimentary cores from the western Pacific display a palaeo-magnetic record of the late Cretaceous long normal interval and the boundary reversed interval corresponding to seafloor spreading anomalies 33–34. Near the young end of this reversed interval, a systematic excursion of inclinations is observed in both cores. Samples are very stable to both alternating field and thermal demagnetization. Blocking temperatures and Curie points suggest that the remanence is carried primarily by magnetite, but with an additional contribution from hematite. Approximate sedimentation rates derived from biostratigraphy suggest that the excursion had a duration of between 46 000 and 54 000yr and occurred about 236000–303000 yr before the succeeding polarity reversal. The excursion, thus, may represent an aborted geomagnetic field reversal.