Investigation of geomagnetic-field secular variation around the Sverdlovsk magnetic observatory

summary . The magnetic observatories at Vysokaya Dubrava and at Arti in Sverdlovsk region have operated in parallel since 1972 January. In spite of the rather small distance between them (160 km), it has been established that during that time, there have been differences in the geomagnetic field variations over a wide range of frequencies. The secular variation also differs appreciably. Over a period of four years the difference in the horizontal component between the observatories increased by 3 nT, the difference in the vertical component by 4 nT and the difference in the total intensity by 6 nT, while the difference in declination decreased by 0.6 min. Besides the normal geomagnetic secular variation, governed by processes in the Earth's core, secular-variation anomalies (SVA) linked with processes in the lithosphere may also take place. The intensity of the local SVAs reaches 10–2 nT/yr, whilst their size does not exceed 15 km. It is undesirable to locate an observatory in a region where an SVA exists since this will distort the normal secular-Variation pattern. It is therefore important to know to what extent observations at the observatory reflect the mean value of the secular variation for that area. For this purpose we would recommend performing a highly accurate survey, for example the total intensity in a network of radial lines in a region around the observatory of radius 50–100 km. Such investigations have been performed around the Arti observatory with an accuracy of the repeated measurements of ±0.3 nT in 1974 and ±0.2 nT in 1975. Ten local SVAs with intensities from 2–14 nT/yr were revealed around the observatory. The observatory of Arti is situated in a quiet region.     

Geomagnetic secular variation and secular acceleration in the Red Sea area

Summary. Estimates of the secular variation in the Red Sea over the period 1959 to 1972 have been obtained from an analysis of marine magnetic data. A total of 318 crossings of ships' tracks were used to determine the mean secular variation for the intervals 1959–72, 1959–65 and 1965–72. The mean secular variation 1959–72 shows a marked northward increase from approximately -10 nT/yr at 13°N to +27 nT/yr at 24°N. North of this, the data suggest a small decrease to + 25 nT/yr at 27° N. These values are consistent with the secular variation recorded at the nearby geomagnetic observatory at Helwan, Egypt, but less than those predicted by the 1965 IGRF for the same period.
Comparison of the mean secular variations for 1959–65 and 1965–72 yields a rough estimate of the secular acceleration of - 1.5 nT yr⁻². Analysis of the cross-over information, corrected for the latitude dependence of the secular variation, shows a regular decrease in the secular change over the period 1959–72 at all latitudes of about -1 nTyr⁻². This secular acceleration makes a substantial contribution to the overall secular change in the Red Sea and as such must be included in the correction of magnetic data covering more than a few years.     

Anomalous variations of the geomagnetic field in East Fergana — magnetic precursor of the May earthquake with M= 7.0 (1978 November 2)

Summary. A secular variation anomaly has been discovered at the north-east part of the Fergana vdey by repeated measurements every year or less. The change of total field Δ F at the 'magnetic epicentre' was 9 nT in 1977 and 16 nT in 1978 relative to the level of 1973. In 1977 an anomalous region was recognized, where according to the data from 25 observation points Δ F increased in the northern part up to 5.2 nT, and decreased by 4.7nTin the southern part according to a further 22 points. Permanent observations were begun at the epicentre in 1978 October. We normally observed variations of Δ F differences with magnitude ± 2–3 nT, which were not correlated with worldwide magnetic activity. Anomalous variations appeared on October 26 and rose to a maximum value of + 23 nT on October 30. The decrease of this anomalous field began on October 31. This made it possible to predict a potential earthquake. The Alay earthquake with M = 7.0 occurred on November 2 six hours after the prediction was issued; Δ F then returned to the initial level. Thus, using the geomagnetic field variations in the Fergana region, geophysicists were able to predict the moment of a strong earthquake.     

Monopoly

Summary. A model for the geomagnetic secular variation field is given, consisting of a series of magnetic monopoles at the surface of the Earth's core. These are distributed according to the density of the data to allow more detailed representation in areas where the density of observations is high, without introducing spurious detail where data are sparse. A monopole model is calculated from observatory secular change data for the epoch 1957.5–1962.5 and its usefulness assessed.     

Centennial to millennial geomagnetic secular variation

A time-varying spherical harmonic model of the palaeomagnetic field for 0–7 ka is used to investigate large-scale global geomagnetic secular variation on centennial to millennial scales. We study dipole moment evolution over the past 7 kyr, and estimate its rate of change using the Gauss coefficients of degree 1 (dipole coefficients) from the CALS7K.2 field model and by two alternative methods that confirm the robustness of the predicted variations. All methods show substantial dipole moment variation on timescales ranging from centennial to millennial. The dipole moment from CALS7K.2 has the best resolution and is able to resolve the general decrease in dipole moment seen in historical observations since about 1830. The currently observed rate of dipole decay is underestimated by CALS7K.2, but is still not extraordinarily strong in comparison to the rates of change shown by the model over the whole 7 kyr interval. Truly continuous phases of dipole decrease or increase are decadal to centennial in length rather than longer-term features. The general large-scale secular variation shows substantial changes in power in higher spherical harmonic degrees on similar timescales to the dipole. Comparisons are made between statistical variations calculated directly from CALS7K.2 and longer-term palaeosecular variation models: CALS7K.2 has lower overall variance in the dipole and quadrupole terms, but exhibits an imbalance between dispersion in   g ¹₂  and   h ¹₂  , suggestive of long-term non-zonal structure in the secular variations.     

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.     

Secular variation of the poloidal magnetic field at the core–mantle boundary

A region of enhanced conductivity at the base of the mantle is modelled by an infinitesimally thin sheet of uniform effective conductance adjacent to the core–mantle boundary. Currents induced in this sheet by the temporally varying magnetic field produced by the geodynamo give rise to a discontinuity in the horizontal components of the poloidal magnetic field on crossing the sheet, while the radial component is continuous across the sheet. Treating the rest of the mantle as an insulator, the horizontal components of the poloidal magnetic field and their secular variation at the top of the core are determined from geomagnetic field, secular variation and secular acceleration models. It is seen that for an assumed effective conductance of the sheet of 10⁸  S, which may be not unrealistic, the changes produced in the horizontal components of the poloidal field at the top of the core are usually ≤10 per cent, but corrections to the secular variation in these components at the top of the core are typically 40 per cent, which is greater than the differences that exist between different secular variation models for the same epoch. Given the assumption that all the conductivity of the mantle is concentrated into a thin shell, the present method is not restricted to a weakly conducting mantle. Results obtained are compared with perturbation solutions.     

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.     

Daily Variation and Secular Variation of the Geomagnetic Field from Shipboard Observations in the Gulf of Aden

Summary The problems of reducing geomagnetic observations from ships at sea in areas influenced by the effect of the equatorial electrojet are discussed. In particular, observations within the Gulf of Aden have been corrected for daily variation and secular variation for the purposes of constructing a contoured magnetic anomaly chart.
An empirical formula is given with which the range of daily variation at different latitudes within the Gulf was estimated for the purpose of correcting the data for daily variation. The observed secular variation, which was used to correct the data, is—11 γ/yr. which differs from the secular variation of +19 γ/yr. in the Gulf of Aden given by the recently adopted International Geomagnetic Reference Field (Zmuda 1969).     

Maximum entropy regularization of time-dependent geomagnetic field models

We incorporate a maximum entropy image reconstruction technique into the process of modelling the time-dependent geomagnetic field at the core–mantle boundary (CMB). In order to deal with unconstrained small lengthscales in the process of inverting the data, some core field models are regularized using a priori quadratic norms in both space and time. This artificial damping leads to the underestimation of power at large wavenumbers, and to a loss of contrast in the reconstructed picture of the field at the CMB. The entropy norm, recently introduced to regularize magnetic field maps, provides models with better contrast, and involves a minimum of a priori information about the field structure. However, this technique was developed to build only snapshots of the magnetic field. Previously described in the spatial domain, we show here how to implement this technique in the spherical harmonic domain, and we extend it to the time-dependent problem where both spatial and temporal regularizations are required. We apply our method to model the field over the interval 1840–1990 from a compilation of historical observations. Applying the maximum entropy method in space—for a fit to the data similar to that obtained with a quadratic regularization—effectively reorganizes the magnetic field lines in order to have a map with better contrast. This is associated with a less rapidly decaying spectrum at large wavenumbers. Applying the maximum entropy method in time permits us to model sharper temporal changes, associated with larger spatial gradients in the secular variation, without producing spurious fluctuations on short timescales. This method avoids the smearing back in time of field features that are not constrained by the data. Perspectives concerning future applications of the method are also discussed.     

Geomagnetic disturbance fields: an analysis of observatory monthly means

SUMMARY
This work quantifies the extent to which disturbance phenomena contribute to the observed geomagnetic field on the time-scale of months to years. A deterministic approach was adopted in which geomagnetic monthly means were analysed in the time domain. By presupposing the time dependence of the external and induced field variations, which we assume to vary as the aa geomagnetic index, we arrived at an estimate of the relative amplitude of the disturbance to each component of the field at a number of observatories. In all some 58 000 geomagnetic monthly means from 59 observatories operating this century were analysed. Special attention was paid to the selection and validation of the data as otherwise the numerical computations would have been unduly compromised.
The results consist of 177 estimates of the relative amplitude of disturbance, one for each of three orthogonal components north ( X ), east ( Y ) and vertically down ( Z ) at each observatory. The disturbance to the X component was found to be consistently negative over the whole of the Earth's surface with an intensification in auroral regions. The disturbances to the Z component were found to be smaller than that for X except in high latitudes. Mean disturbances to the Y component were smaller still. Results were in general consistent with the dipole field of the magnetospheric ring current, aligned with, though of opposite polarity to, the Earth's main field. Typical amplitudes of the mean disturbance field from month to month were of the order 10 nT.
The results can be used to estimate the variation of the disturbance field. Subtracting this from both monthly and annual means yields an improved estimate of the field originating in the core and its secular variation. Some examples are presented.     

The Canary Islands swell: a coherence analysis of bathymetry and gravity

The Canary Archipelago is an intraplate volcanic chain, located near the West African continental margin, emplaced on old oceanic lithosphere of Jurassic age, with an extended volcanic activity since Middle Miocene. The adjacent seafloor does not show the broad oceanic swell usually observed in hotspot-generated oceanic islands. However, the observation of a noticeable depth anomaly in the basement west of the Canaries might indicate that the swell is masked by a thick sedimentary cover and the influence of the Canarian volcanism. We use a spectral approach, based on coherence analysis, to determine the swell and its compensation mechanism. The coherence between gravity and topography indicates that the swell is caused by a subsurface load correlated with the surface volcanic load. The residual gravity/geoid anomaly indicates that the subsurface load extends 600 km SSW and 800 km N and NNE of the islands. We used computed depth anomalies from available deep seismic profiles to constrain the extent and amplitude of the basement uplift caused by a relatively low-density anomaly within the lithospheric mantle, and coherence analysis to constrain the elastic thickness of the lithosphere ( T_e ) and the compensation depth of the swell. Depth anomalies and coherence are well simulated with T_e =28–36 km, compensation depth of 40–65 km, and a negative density contrast within the lithosphere of ∼33 kg m⁻³. The density contrast corresponds to a temperature increment of ∼325°C, which we interpret to be partially maintained by a low-viscosity convective layer in the lowermost lithosphere, and which probably involves the shallower parts of the asthenosphere. This interpretation does not require a significant rejuvenation of the mechanical properties of the lithosphere.     

Geomagnetic evidence for fluid upwelling at the core-mantle boundary

Summary. Previous studies, both geomagnetic and seismic, have been unable to show conclusively whether or not there is fluid upwelling at the core-mantle boundary. Here a new method is developed, in which an attempt is made to invert geomagnetic secular variation data measured at the Earth's surface for a frozen-flux purely toroidal core-mantle boundary (CMB) velocity field, under the assumption that the mantle is electrically insulating and flux is frozen in at the CMB. These data have previously been inverted for the core-mantle boundary radial secular variation, from which the appropriate fit between model and data is known. Two different main field models were used to assess the effect of uncertainty in its radial component at the CMB. The conclusions were the same in both cases: frozen-flux purely toroidal motions provide a poor fit. A statistical test allows very firm rejection of the hypothesis that the residuals are not significantly larger, whereas there is no statistical difference between the residuals of inversions for radial secular variation and frozen-flux velocity fields at the CMB if upwelling and down-welling is included. The inherent non-uniqueness in the velocity field obtained is not of concern, since only their statistical properties are utilized and no physical significance is attached to the flows obtained.     

Discrete scale invariance connects geodynamo timescales

The geodynamo exhibits a bewildering gamut of time-dependent fluctuations, on timescales from years to at least hundreds of millions of years. No framework yet exists that comprises all and relates each to all others in a quantitative sense. The technique of bootstrapped discrete scale invariance quantifies characteristic timescales of a process, based upon log-periodic fits of modulated power-law scaling of size-ranked event durations. Four independent geomagnetic data sets are analysed therewith, each spanning different timescales: the sequence of 332 known dipole reversal intervals (0–161 Ma); dipole intensity fluctuations (0–2 Ma); archeomagnetic secular variation (5000 B.C.–1950 A.D.); and historical secular variation (1590–1990 A.D.).
Six major characteristic timescales are empirically attested: circa 1.43 Ma, 56 Ka, and 763, 106, 21 and 3 yr. Moreover, all detected wavelengths and phases of the detected scaling signatures are highly similar, suggesting that a single process underlies all. This hypothesis is reinforced by extrapolating the log-periodic scaling signal of any particular data set to higher timescales than observed, through which predictions are obtained for characteristic scales attested elsewhere. Not only do many confirm one another, they also predict the typical duration of superchrons and geomagnetic jerks. A universal scaling bridge describes the complete range of geodynamo fluctuation timescales with a single power law.     

Improved magnetic anomalies of the Antarctic lithosphere from satellite and near-surface data

The Antarctic magnetic anomaly map compiled marine and airborne surveys collected south of 60°S through 1999 and used Magsat data to help fill in the regional gaps between the surveys. Ørsted and CHAMP satellite magnetic observations with greatly improved measurement accuracies and temporal and spatial coverage of the Antarctic, have now supplanted the Magsat data. We combined the new satellite observations with the near-surface survey data for an improved magnetic anomaly map of the Antarctic lithosphere. Specifically, we separated the crustal from the core and external field components in the satellite data using crustal thickness variations estimated from the terrain and the satellite-derived free-air gravity observations. Regional gaps in the near-surface surveys were then filled with predictions from crustal magnetization models that jointly satisfied the near-surface and satellite crustal anomalies. Comparisons in some of the regional gaps that also considered newly acquired aeromagnetic data demonstrated the enhanced anomaly estimation capabilities of the predictions over those from conventional minimum curvature and spherical harmonic geomagnetic field models. We also noted that the growing number of regional and world magnetic survey compilations involve coverage gaps where these procedures can contribute effective near-surface crustal anomaly estimates.     

t*— an unsuitable parameter to characterize anelastic attenuation in the Eastern Carpathians

Investigation of teleseismic P -wave recordings at a temporary network in the Eastern Carpathians, equipped with predominantly short-period sensors, is compared with synthetic modelling of anelastic attenuation of teleseismic waves in the upper mantle. Using the t * approach, we examine variations of amplitude decrease over frequency for teleseismic recordings in the frequency band 0.5–1.5 Hz. The results reveal a consistent pattern of increased t * values in the centre of the network, in the Vrancea region at the bend of the Carpathian Arc, although the magnitude of the observed variation in t * is much higher than expected. Synthetic t * parameter computations for the same event-receiver configurations reproduce the observed pattern in terms of relative variations. However, the amplitude of the synthetic t * values explains only 10–20 per cent of the observed variation in t *. t * is not a direct measure for anelastic attenuation but rather for a combination of anelastic and other attenuating effects such as scattering and amplitude fluctuation related to velocity inhomogeneities. If regional amplitude variations are solely attributed to anelastic attenuation, all other effects are mapped into Q . We discuss the role of anelastic attenuation and other effects in the case of the Eastern Carpathians and conclude that t * is an unsuitable parameter to characterize anelastic attenuation in the Eastern Carpathians.     

Apparent polar wander of the mean-lithosphere reference frame

Apparent polar wander in the mean-lithosphere (= no-net-rotation = no-net-torque uniform drag) reference frame is compared with apparent polar wander in the hotspot reference frame over the past 100 Myr. Palaeo-magnetic poles and plate rotations previously used to determine an apparent polar wander path for the hotspot reference frame are here used to determine an apparent polar wander path in the mean-lithosphere reference frame. We find that the two paths are similar, especially for Late Cretaceous time, when a 10°–20° shift of the pole occurred. To first-order the hotspots and lithosphere (as a whole) moved in unison relative to the palaeomagnetic axis during Late Cretaceous time. A non-dipole field explanation for the apparent shift can probably be excluded. However, either motion of the time-averaged geomagnetic axis relative to the spin axis or polar wandering could have caused this shift, the latter being the more likely explanation.     

Lake sediment record of the geomagnetic secular variation in Britain during Holocene times

Summary. A palaeomagnetic record of the geomagnetic secular variation during the last 10000 years has been obtained from 10 cores of sediment from Loch Lomond, Scotland, Lake Windermere, North England, and Llyn Geirionydd, North Wales. A time-scale is provided by 30 radiocarbon age determinations and pollen analyses on several of the cores. The main swings and much fine detail of both declination and inclination records repeat well between cores and between lakes, and the overall record is much more detailed than previous European records.
The new record shows that neither declination nor inclination swings have been periodic over the past 10000 years, but that the main swings have become progressively shorter in duration during that time. Each swing is characterized by fine detail which enables use of the record as a secondary method of dating other European sediments.
The motion of the geomagnetic vector has been predominantly clockwise throughout the time period spanned, but confirms a period of anticlockwise motion from 1100 to 600 bp first discovered by British archaeomagnetic investigations. The record agrees with British and Czechoslovakian archaeomagnetic records, but not with Japanese archaeomagnetic or North American lake sediment records. This suggests that the secular changes are controlled by local growing and decaying, drifting sources, rather than by wobbling of the main geomagnetic dipole.     

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.