Secular variation in North America during historical times

Summary. The secular variation of the declination, inclination and total force of the geomagnetic field has been plotted for 74 locations in North America. A comparison of the occurrences of maxima and minima in the SV curves from different stations shows little evidence of drift in North America. Although a declination maximum exhibits westward drift up to 1915, all other extrema in declination, inclination and total intensity occur almost simultaneously over a wide area. The major feature of SV in North America appears to be a 4000 nT decrease in the total field since 1850, which may be due to a decrease in the dipole moment coupled with the decay of a large non-dipole anomaly situated under the continent. Short-period changes in the rate of decrease are possibly jerks of the magnetic field. Maximum entropy spectral analysis of all three components of the field indicates periods of 102 and 53 yr.     

Geomagnetic secular variation since 1901

Summary. All available annual means, from the world-wide network of magnetic observatories, of north intensity ( X ), east intensity ( Y ) and vertical intensity (Z) from 1901 to 1977 are subjected to spherical harmonic analysis to obtain 38 models of the Earth's geomagnetic field at two-year intervals. These models are differenced to give 37 models of secular variation at two-year intervals from 1903.5 to 1975.5. The results show the decreasing trend of the dipole moment and are analysed for possible information on the westward drift of the magnetic field.     

Polarity Reversals and Oxidation State in Rocks

The correlation between reversed polarity in rocks and strong oxidation of magnetic minerals is briefly re-examined, with especial reference to a well-documented Columbia Plateau basalt suite, recently described in this journal. The correlation appears to vary in quality from place to place and some reasons for this are suggested. Attention is particularly drawn to recent studies which suggest that for palaeomagnetic purposes there may be two groups of lavas, divisible on the basis of whether their magnetic minerals became oxidized under high or low temperature conditions.
A sizable body of opinion now holds (not necessarily correctly) that self-reversals are a negligible factor in the correlation between oxidation and reversed polarity, and that some compositional effect is most probably involved. Virtually the only known petrographic criterion which can (statistically) distinguish reversely from normally magnetized rocks is their (statistically) higher oxidation states. An obvious implication is that the water content of magmas, and perhaps of high-grade metamorphic rocks, is somehow increased during episodes of reversed polarity, so that a slightly larger proportion of rocks is more strongly oxidized.
Since motions of the Earth's core, upper mantle processes, and field reversals are apparently all related, compositional changes originating in the upper mantle could well be also related to field reversals. Alternating linear zones of more and less altered baslat recently reported from one place near the mid-Atlantic ridge might have arisen in such a way.
The importance of relating palaeomagnetic measurements to detailed mineralogical studies of magnetic minerals is again stressed.     

Anisotropic turbulence in weakly stratified rotating magnetoconvection

Numerical simulations of the 3-D MHD-equations that describe rotating magnetoconvection in a Cartesian box have been performed using the code NIRVANA. The characteristics of averaged quantities like the turbulence intensity and the turbulent heat flux that are caused by the combined action of the small-scale fluctuations are computed. The correlation length of the turbulence significantly depends on the strength and orientation of the magnetic field and the anisotropic behavior of the turbulence intensity induced by Coriolis and Lorentz force is considerably more pronounced for faster rotation. The development of isotropic behavior on the small scales—as it is observed in pure rotating convection—vanishes even for a weak magnetic field which results in a turbulent flow that is dominated by the vertical component. In the presence of a horizontal magnetic field the vertical turbulent heat flux slightly increases with increasing field strength, so that cooling of the rotating system is facilitated. Horizontal transport of heat is always directed westwards and towards the poles. The latter might be a source of a large-scale meridional flow whereas the first would be important in global simulations in case of non-axisymmetric boundary conditions for the heat flux.     

Equation for the field lines of an axisymmetric magnetic multipole

Summary. An exact equation is derived for the magnetic field lines of the general axisymmetric magnetic multipole of arbitrary degree ( n ). This new result has important applications in studies of the possible nature of solarterrestrial physics during geomagnetic polarity reversals. In the limiting case of a magnetic dipole ( n=1 ), the equation for the magnetic field lines of the general axisymmetric magnetic multipole simplifies correctly to the well-known dipolar form, which is used extensively in geomagnetism, magnetospheric physics and cosmic-ray physics as a first-order approximation to the actual configuration of the geomagnetic field.
It is also shown theoretically that suites of similar magnetic field lines of the general axisymmetric multipole attain their maximum radial distances from the origin on a set of circular conical surfaces, with coincident vertices at the centre of the Earth; this set includes the equatorial plane if the degree ( n ) of the multipole is odd. The magnetic field is horizontal everywhere on all these surfaces.
Palaeomagnetic studies have suggested that during some polarity reversals the magnetic field in the inner magnetosphere can be represented approximately either by a single, non-dipolar, low-degree (2 < n < 4), axisymmetric magnetic multipole or by a linear combination of such multipoles. In this situation, the equation for the field lines of an axisymmetric magnetic multipole of low degree (2 < n < 4) would be as fundamental to a proper understanding of magnetospheric, ionospheric and cosmic-ray physics during polarity reversals as is the equation for dipolar field lines in the case of the contemporary geomagnetic field.     

Late Mesozoic and Cenozoic palaeomagnetism of Australia – III. Bias-corrected pole paths for Australia, Antarctica and India

Summary. Palaeomagnetic results from Part I of this study and their analysis in Part II are combined to eliminate bias from the Cenozoic apparent polar wander path for Australia – a bias due to non-dipole components in past geomagnetic fields or, for poles calculated from hot-spot data, due to the motion of hot spots relative to the Earth's rotational axis. This path is extended in approximately bias-free form to the late Mesozoic, and indicates a significant change in the drift direction of the continent between 26 and about 60 Ma.
The bias-corrected Australian path is used, first, with seafloor spreading data for the Southern Ocean to derive a corresponding late Mesozoic–Cenozoic pole path for Antarctica. The latter shows that the Antarctic drift direction reversed in the early Tertiary. It is suggested that the early Tertiary directional changes of both Australia and Antarctica are part of a global reorganization of plates during the Eocene, postulated by Rona & Richardson, Cande & Mutter and Patriat & Achache.
Next, the Australian path is compared with hot-spot data from the African and Australian plates, indicating a movement of the hot spots relative the Earth's rotational axis during the Cenozoic. The direction of this movement is found to be consistent with previous results from other parts of the world.
Finally, the Australian path is used together with non-dipole components in the geomagnetic field to explain a prominent westward displacement of the mid- and late Cenozoic poles of India relative to those of Australia.
Because of uncertainties in the original poles and in the analysis, the present results are likely to contain appreciable errors. Nevertheless, their consistency with independent findings supports the dipole-quadrupole model of Part II for mid- and late Cenozoic geomagnetic fields.     

11 million years of Oligocene geomagnetic field behaviour

An 11 million year long record of the Oligocene geomagnetic field has been obtained from pelagic sediments of DSDP Hole 522 An average sample spacing of 4 cm yielded approximately one specimen per 4 to 8 kyr. The rock magnetics are remarkabh consistent across the entire interval. Previous work demonstrated a magnetic mineralogy dominated by magnetically stable magnetite. The natural remanent magnetism (NRM) carries an Oligocene polarity timescale that is in excellent agreement with the Oligocene reversal record as determined from marine magnetic anomalies (MMAs), including many of the so-called 'crypto-chrons'. Normalized NRM intensities from the undisturbed portions of the record yield a time series of variations with features consistent with a number of other palaeointensity time series derived from both sedimentary and lava sequences. These features include consistent, major decreases in palaeointensity (DIPs) at reversal boundaries, and occasional DIPs between reversal boundaries that could correspond to lineated 'tiny wiggles' in the MMA records. The data set suggests that the overall field strength was 40 per cent higher in the first half of the Oligocene when the average reversal frequency was 1.6 Myr^-1 than in the second half when the reversal frequency was 4 Myr^-1. There is also a weak dependence of average field strength on length of polarity interval. Finally, in the three cores suited to spectral analysis (of coherent polarity and relative intensity independent of lithological contamination), there is a persistent ca. 30–50ka periodicity in the variations of the relative intensity, suggesting that the geomagnetic field 'pulses' at about this frequency, not only during the Brunhes (as demonstrated by Tauxe & Shackleton 1994), but in the Oligocene as well.     

Late Mesozoic and Cenozoic palaeomagnetism of Australia–II. Implications for geomagnetism and true polar wander

Summary. New palaeomagnetic results from Australia indicate that throughout the Cenozoic era the continent lay further south than suggested by hot-spot data. Moreover, while hot spots give a uniform rate of drift during most of the Cenozoic, the drift rate obtained from apparent polar wander varies considerably.
The discrepancies between the palaeomagnetic and hot-spot results are analysed by comparing the Australian data with those of Europe and the central Pacific. The analysis suggests that the discrepancies are due to: (1) departures of the Earth's magnetic field from the geocentric axial dipole model, and (2), either true polar wander or a non-axial inclined dipole component. It is found that since the mid-Tertiary the dominant non-dipole component has been a quadrupole, and that during this period both the quadrupolar field and the true polar displacement/non-axial dipole component decreased progressively. During the Quaternary, and also at the earliest Tertiary, the non-dipole components appear to have been moderate or small.
The comparison of data sets demonstrates that considerable errors may be incurred when Cenozoic, and presumably earlier, poles from one geographic region are used to derive those of another, widely separated, region. The results also imply that absolute plate velocities estimated from palaeomagnetic data can contain substantial errors, and that hot-spot data may need significant adjustments for true polar wander to yield correct palaeolatitudes.
Finally, the new early Tertiary pole for Australia is used in conjunction with updated early Tertiary poles from other lithospheric plates to reapply the McKenzie test for true polar wander. The results indicate a small true polar displacement since the beginning of the Tertiary. The amount and direction of the displacement, however, differ from those generally obtained from hot-spot data.     

Rotation of the geomagnetic field about an optimum pole

Since 1693, when Halley proposed that secular change was the result of the westward drift of the main field, his simple model has undergone many refinements. These include different drift rates for dipole and non-dipole parts; separation into drifting and standing parts; latitudinal dependence of drift rate; northward drift of the dipole; and non-longitudinal rotations of the individual harmonics of the geomagnetic field. Here we re-examine the model of Malin and Saunders, in which the main field is rotated about an optimum pole which does not necessarily coincide with the geographical pole. The optimum pole and rotation angle are those that bring the main field for epoch T ₁ closest to that for T ₂ , as indicated by the coefficients of correlation between the spherical harmonic coefficients for the two epochs, after rotation. Malin and Saunders examined the pole positions and rates of rotation using data from 1910 to 1965, and noticed a number of trends. We show that these trends are confirmed by recent IGRF models, spanning the interval 1900–2000 and to degree and order 10. We also show that the effect of the level of truncation is small.     

Time series analysis of the 10 000 yr geomagnetic secular variation record from SE Australia

Summary. A 10 000 yr continuous secular variation record from intensively dated lake sediments in SE Australia has been subjected to periodogram and maximum entropy method analysis. Tests on synthetic data reveal some of the limitations of the latter method, particularly when applied to complex number series. Anticlockwise precession of the magnetic vector at a period of 5000 ± 1000 yr is tentatively ascribed to dipole precession, and clockwise precession at a period of about 2800 yr is probably due to westward drift of features of the non-dipole field.
The effect of calibrating the radiocarbon time-scale is important and results in periodicity shifts of up to 25 per cent. Even for well-dated lacustrine sequences power spectra are poorly constrained: it is thus possible that the geomagnetic secular variation on a time-scale of thousands of years is more uniform than often supposed. Mismatches between declination and inclination spectra can arise as a natural consequence of certain types of source mechanism and should not be simply attributed to figments of the analysis employed.     

Magnetic viscosity, diffusion after-effect, and disaccommodation in natural and synthetic samples

Summary. The magnetic viscosity of some oceanic basalts from IPOD site 417D is investigated. The samples are fine-grained pillow basalts and coarse-grained massive flows. Experiments on the viscous behaviour of these samples include (1) acquisition and decay of weak-field viscous remanent magnetization (VRM) and (2) the influence on VRM produced by the zero field storage time following demagnetization by an alternating magnetic field. In all cases the aquisition of VRM is found to be more rapid than its subsequent decay. The time dependence of the intensity of VRM is observed to be significantly non-logarithmic during acquisition, but less so during decay. The length of time a sample remains in zero field following demagnetization produces a decrease in the aquisition and decay coefficients of VRM. These results cannot be reconciled entirely with predictions based on existing models of magnetic viscosity due to thermal fluctuations. Available evidence also suggest that the effect on VRM produced by (1) zero field storage following demagnetization or produced by (2) different initial states of magnetization, is attributed to susceptibility disaccommodation. Unfortunately, the available data are insufficient fully to ascertain the importance of diffusion after-effect and disaccommodation to long-term viscous behaviour. It is suggested, however, that diffusional processes may contribute to VRM in young pillow basalts (< 0.1 Ma) near mid-ocean ridges where the oxidation of titanomagnetite produces an increase in the number of octahedral lattice site vacancies and the elevated temperature of the crust enhances the diffusive reordering of these vacancies and ferrous ions. In order to gain some additional insights into VRM and attempt to explain some of these observations, theoretical and experimental results pertaining to diffusion after-effect and disaccommodation in magnetic materials are reviewed.     

Palaeointensity studies of a late Permian lava succession from Guizhou Province, South China: implications for post-Kiaman dipole field behaviour

Palaeomagnetic results are presented from a volcanic sequence in Zhijin County, Guizhou Province, Southern China. The lavas and associated volcanic breccias comprising the sequence represent a southern extension of the Emeishan volcanic province. Biostratigraphic dating of interbedded limestone units and stratigraphic constraints indicate that the section formed during the late Permian (∼263–255  Ma), and is thus somewhat older than the Emei stratotype section in Sichuan Province, and close in age to reported estimates of the termination of the Permo-Carboniferous (Kiaman) reverse superchron. Rock magnetic analyses and reflected light microscopy indicate that the magnetic mineralogy of the lava units is dominated by fresh, primary magnetites containing a significant fraction of single-domain grains. Thermal demagnetization behaviour of the breccia units is poor, but most lava samples have one or two components of remanence above 250 °C. The normal polarity characteristic remanence held by the lavas implies a post-Kiaman age for this succession and suggests that the termination of the Kiaman occurred prior to 263  Ma, supporting recently published estimates. According to standard criteria, Thellier palaeointensity results from the lavas are of good quality and reveal that the dipole field strength was comparatively low shortly after the termination of the superchron. 80 per cent of samples record relative VDM values in the range 42–52 per cent of the present-day value, supporting recent studies of mid-Kiaman field intensity. This suggests that a low-energy dipole existed at least between 300 and 255  Ma and does not appear to have been confined to the stable reverse polarity interval.     

Ordovician magnetostratigraphy: Llanvirn-Caradoc limestones of the Baltic platform

b
Palaeomagnetic sampling has been performed covering 43 stratigraphic levels within the Baltoscandian Ordovician carbonates. After removing a ubiquitous Permo-Carboniferous (287 ± 14 Ma) remagnetization between 200 and 500 C, a Llanvirn-Caradoc reversal stratigraphy is delineated by components with maximum unblocking temperatures up to 550-580 C. Three reversed (SE, down) and three normal (NW. up) antipodal polarity intervals have been recognized. A primary/early diagenetic remanence age is therefore inferred for the stratigraphically linked polarity chrons. Primary magnetizations are resident in detrital/biogenic or early diagenetically formed single- and pseudo-single domain magnetite phases and subordinate early diagenetic pigmentary haematite.
The recognition of a primary remanence within these well-dated Ordovician carbonates has the following important tectonic and magnetostratigraphic consequences.
(1) Accurate time-calibration of the Baltic APW path implies that rapid counterclockwise rotation took place in late Tremadoc and Llandeilo times. The Arenig-Llanvirn epochs are characterized by a 'still stand'. Baltica occupied intermediate to high Southerly latitudes during the early Ordovician (Tremadoc-Llanvirn). Systematic northward drift is recognized in post-Llanvirn times.
(2) A time-calibrated Ordovician reversal stratigraphy is proposed. Presently available data suggest the geomagnetic field was predominantly reversely polarized during Tremadoc and Arenig times. Two normal polarity zones of short duration are identified within mid-Llanvirn and mid-Llandeilo strata. Discontinuities within the succession may mask other short-period events. Late Llandeilo to mid-Caradoc times were then characterized by a normal polarity field.     

Secular variation of geomagnetic field intensity during polarity transitions

Summary. The palaeointensity of the geomagnetic field during some ancient polarity transitions as well as and during the Jaramillo event was determined by the method of alternating-field remagnetization. The periods of the secular variation during these transitions were also determined.     

Upper Cretaceous palaeomagnetic stratigraphy at Moria (Umbrian Apennines, Italy): verification of the Gubbio section

Summary. The Gubbio section, in the Umbrian Apennines, yielded a sequence of Upper Cretaceous polarity zones precisely corresponding to the sea-floor magnetic anomalies of that age; the Gubbio polarity zones have been dated with planktonic foraminifera. As a test of the validity of the Gubbio reversal sequence, we have studied the section at Moria, 16 km to the north. The polarity sequence at Moria closely agrees with the sequence at Gubbio; only the two thinnest polarity subzones were not found. Comparison of polarity zone thicknesses in the two sections gives an excellent correlation. Declinations in the upper part of the Moria section are somewhat erratic, suggesting the possibility of detachment surfaces and relative structural rotations at these levels. Remanent magnetic intensity varies in a cyclical pattern that is unmistakably the same in the Gubbio and Moria sections. The results from the Moria section provide strong confirmation of the validity of the magnetic polarity sequence established at Gubbio.     

Use of the frozen flux approximation in the interpretation of archaeomagnetic and palaeomagnetic data

Summary. The frozen flux approximation of Roberts & Scott is a constraint on the core field that can be used to aid interpretation of the very sparse datasets that palaeomagnetism and archaeomagnetism provide. It gives bounds on the size of the components of the magnetic field at a point, of the Gauss coefficients, and, if valid over such long time periods, limits the shape of the field during transitions between normal and reversed polarities. The maximum intensity at a point, consistent with the present flux, is 281 μT or 4 times the maximum field observed today. The present dipole is about 50 per cent of its upper bound. Polarity reversal is impossible if the transition field is purely axisymmetric. None of the measurements we consider violate the frozen flux approximation.     

Geomagnetic long-term secular variations in Italian Lower Cretaceous shallow-water carbonates

The remanent magnetic properties of an 88  m bore core are unrelated to either the dolomite content or the sedimentological textures and are considered to be carried primarily by biogenetic magnetite that was cemented in during very early diagenesis. Individual readings represent time intervals of c . 720 ± 32  yr and, after 40  mT partial demagnetization, they provide an almost continuous record of averaged geomagnetic secular variations over a period of some 3.17  Myr. The magnitude of directional secular variation is twice that of the present day, despite being smoothed, and the secular variations appear to grade into polarity transitions, suggesting no difference in their mechanisms. The rates of change in direction between subjacent levels in the core have a log-normal distribution which extends smoothly beyond 90° and has a median value of 13°/700  yr, the same as for unsmoothed European secular variation during the last 2000  yr. The intensity of remanence, after 40  mT partial demagnetization, appears to provide a reasonable approximation to geomagnetic field intensity. This tends to be weaker when the direction is moving faster, reflecting averaging, but is unrelated to the distance of the vector from the mean direction; that is, it depends on the rate of change and not on the virtual pole latitude. The virtual poles, after correction for tectonic rotations about horizontal and vertical axes, have latitudes that form a strongly platykurtic Fisherian distribution, while their longitudes have a circular distribution on which are superimposed two Gaussian peaks, 180° apart. This bore core thus provides detailed information of smoothed geomagnetic secular variation in the Lower Cretaceous (127 ± 3  Ma) which shows clear regularities in behaviour, some related to changes in the Earth's orbital parameters.     

Palaeointensity In the Pelagic Realm: Marine Sediment Data Compared With Archaeomagnetic and Lake Sediment Records

Summary. Four box cores collected from the Ontong—Java plateau during the Eurydice expedition have been used to make relative geomagnetic palaeo-intensity measurements. Rock magnetic measurements on the sediments show that they are characterized by a uniform magnetic mineralogy, and that they are suitable for relative intensity estimates. These are obtained by normalizing the NRM by an ARM imparted in a low DC bias field. the palaeoceanographic event known as the preservation spike is used to establish a crude time-scale for the record so that it may be compared with other data from the same region, and also with global palaeointensity estimates. the marine sediment data are quite similar to Australian intensity data from lake sediments and archaeomagnetic sources, but as might be expected exhibit some obvious differences from the global record.     

On the interpretation of frequency response functions for oceanic gravity and bathymetry

Summary. In 1972–78, the late N. D. Watkins and others carried out ajoint field programme of geological mapping in the Mio-Pliocene flood basalts of Iceland, including sampling for K-Ar dating and palaeomagnetic research. The major part of the palaeomagnetic sampling is represented by 2462 lavas in five long composite sections through the lava pile. This paper deals with various statistical properties of this data set.
It is concluded that geomagnetic reversals occur more frequently than is assumed in the current ocean-floor polarity time-scale. There is no evidence for significant asymmetries between normal and reverse polarity states of the field, neither as regards chron lengths, secular variation, or virtual dipole moment magnitude. Intensities of remanence in these lavas are shown to be well approximated in terms of a hyperbolic distribution. The latitude distribution of virtual magnetic poles can be fitted with a Bingham function having k' ∼ 4.5, and low-latitude poles do not occur preferentially in any particular longitude interval.     

A study on the multi-component substrom current

A new technique of eigen mode analysis, Method of Natural Orthogonal Components (MNOC) is used to analyze the ionospheric equivalent current systems obtained on the basis of magnetic data at six meridian magnetometer chains in the northern hemisphere during March 17 19, 1978. The results show that the whole current pattern for any given instant consists of a few eigen modes with different intensities. The first eigen mode exhibits a two cell current construction, characterizing the large scale magnetospheric convection and directly driven process, while the second eigen mode shows a concentrated westward electrojet at midnight sector, characterizing the substorm current wedge and the loading unloading process. The first mode consistently exists whenever during quiet periods or at substorms, and its intensity increases from the beginning of the growth phase of substorms, then quickly intensifies in the expansion phase, followed by a gradual decrease in the recovery phase. On the other hand, the intensity of the second mode remains to be near zero during both quiet time and the growth phase of substorms. Its rapid enhancement occurs in the expansion phase. These characteristics in the current patterns and the intensity variations coincide with the defined physical processes of the directly driven and loading unloading components.