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.     

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.     

An investigation of the magnetic variation within two recent lava flows

Samples collected in vertical sections from two recent (1169 AD and 1971 AD) lava flows show significant variations in both the palaeomagnetic vector and the rock magnetic characteristics. The age and the short cooling time of the flows (both are less than 3 m thick) are such that the observed range of palaeomagnetic variation is more than can reasonably be explained by secular variation. The variations in palaeomagnetic and rock magnetic characteristics appear in part to be a response to the relative proportions of two coexisting ferrimagnetic phases, one that displays a primary coarsegrained and Ti-rich composition and a second that has finer-grained, Ti-poor characteristics. The use of AF demagnetization in this study means that the latter phase will have a disproportionate influence on the palaeomagnetic results. It has been suggested elsewhere that in such a magnetic system the Ti-poor phase may represent a CRM formed at some temperature below its final Curie temperature; if so, the results of this study imply that not only does it affect the palaeointensity estimate (which requires a pure TRM) but it appears also to influence the palaeomagnetic direction to an extent that far exceeds the normal secular change that might be expected to occur during delayed acquisition. Such a result has important implications for palaeomagnetic studies that use lavas, both in the selection of a suitable sampling strategy and in the geomagnetic interpretation of palaeomagnetic variation within single flow units.     

Palaeomagnetism and radiometric dating of the British Tertiary Igneous Province: Muck and Eigg

Summary. The lavas, dykes and other intrusive bodies of Muck and Eigg have been sampled for palaeomagnetic purposes. All the bodies are reversely magnetized except for four dykes of Muck which are normally magnetized. The lavas show variations of direction of magnetism believed to be a record of secular variation.
The radiometric dates support the geological evidence that the lavas of Eigg and Muck predate parts of the Rhum central igneous complex, the lavas of Rhum and Canna and some of those of Skye. The combined evidence shows that most of the activity took place during an R-N-R-N sequence between 63 and 58 Ma. The younger R-N sequence is the same as the first part of the R-N-R sequence found in Mull but the Sgurr pitchstone of Eigg and parts of the Eastern Redhills of Skye belong to later reversed intervals and correspond to Chron 26r and Chron 26. The igneous activity occurred over a 10 Myr period from about 63 Ma and can probably be associated with Chrons 27r to 23r of the Harland et al. time-scale.     

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.     

Geomagnetic secular variation in historical lavas from the Canary Islands

Summary. The age of historical lava flows on the Canary Islands can be traced back at least to 1585 ad . During the last 400 years eruptions occurred at fairly regular intervals. The palaeomagnetic directions of these lavas record the ambient field direction with high precision so that an historical secular variation curve (SVC) can be constructed for the Canaries which closely resembles the curve registered at Paris. The SVC can be used as a dating tool for historical eruptions of unknown age and can probably be extended into the last 4000–5000 yr of recent volcanic activity on the Canaries.     

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.     

Magnetostratigraphy of Palaeocene basalts from the Vaigat Formation of West Greenland

A palaeomagnetic study comprising the directional results from 289 individual lava flows, sampled along eight sections in the Palaeocene basalts of West Greenland, is reported. The eight individual sections are correlated using lithostratigraphical marker horizons to form a single composite profile. Generally, the lithological correlation is in good agreement with the record of geomagnetic secular variation.
  The total composite palaeomagnetic profile represents a stratigraphic thickness of 1.6  km through the Vaigat Formation, which is the lowermost of the two volcanic formations formed during the main stage of plateau volcanism. Only two polarity zones are found in the composite profile, suggesting a very short duration for the West Greenland main plateau-building volcanism. ⁴⁰Ar/³⁹Ar dates support a high extrusion rate and also indicate that the lower normal polarity zone is Chron C27n and that the upper reverse polarity zone is Chron C26r.
  The C27n–C26r transition is fully recorded along one of the sections (Nuusap Qaqqarsua), with intermediate directions covering a 200  m thick succession of lavas. A combined palaeomagnetic, field and geochemical study along this profile showed good agreement; that is, geochemically and geologically derived single magmatic events show groupings of the palaeomagnetic directions. Supposing a duration for the geomagnetic transition of 5000 years, the eruption frequency during this period was as high as one flow every 80 years.     

The magnitude of the palaeomagnetic field in Iceland between 2 and 6 Myr ago

Summary. Sixty-eight palaeomagnetic field magnitude values have been determined from a sequence of Icelandic lavas, ranging from 2 to 6 Myr in age. The results indicate large and rapid changes in the palaeomagnetic dipole field and provide a mean value of the palaeomagnetic field magnitude in Iceland for this period.     

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.     

A physical model for palaeosecular variation

Summary. A new model to describe the latitude dependence of the angular dispersion of the palaeomagnetic field (palaeosecular variation) is developed following previous models, but with crucial differences. It is shown that if the probability distribution of virtual geomagnetic poles (VGPs) is circularly symmétric about the rotation axis then the geométry of the distribution of field directions is latitude dependent. This has a significant effect on the latitude dependence of dispersion and is accounted for in the model. The dipole and non-dipole parts of the field are not artificially separated but are intimately linked through an observationally determined relation that the time averaged intensity of the non-dipole field is dependent upon the intensity of the dipole field. It is shown that a consequence of this relation is that no knowledge of the probability distribution of the geomagnetic dipole moment is required. This is a fundamental improvement over previous models.
The model provides excellent fits to the palaeodata and, unlike previous models, is not inconsistent with the latitude variation of the non-dipole field dispersion determined from the present field. For the past 5 Ma the point estimate of the VGP dispersion due to dipole wobble is 7.2° and of the VGP dispersion at the equator due to variation in the non-dipole field is 10.6°. This estimate of the dispersion due to variation in the non-dipole field is in excellent agreement with the value predicted from an analysis of the variation in field intensities over the same period. Fits of the model to data from earlier periods indicate that dispersion due to variation in the non-dipole field is essentially independent of the geomagnetic reversal rate while dipole wobble is positively correlated with reversal rate.     

A definitive model of the geomagnetic field and its secular variation for 1965 — I. Derivation of model and comparison with the IGRF

Summary. Using a very large body of post-1955 data, a spherical harmonic model of the geomagnetic field and its secular variation is derived for 1965.0. This model is compared with the original International Geomagnetic Reference Field (IGRF) and with individual models used, or proposed for use, in producing the IGRF. Positions of the dip-poles, the geomagnetic poles and the eccentric dipole are derived from the model, together with their rates of change, and comparisons are made with other estimates of these positions.     

Palaeosecular variation over the past 5 Myr based on a new generalized database

A generalized database of global palaeomagnetic data from 3719 lava flows and thin dykes of age 0–5 Ma has been constructed for use with a relational database management system. The database includes all data whose virtual geomagnetic poles (VGP) lie within 45 of the spin axis and can be used for studies of palaeosecular variation and for geomagnetic field modelling. Because many of these data were collected and processed more than 15–20 years ago, each result has been characterized according to the demagnetization procedures carried out. Analysis of these data in terms of the latitude variation of the angular dispersion of VGPs (palaeosecular variation from lavas) strongly suggests that careful data selection is required and that many of the older studies may need to be redone using more modern methods. Differences between the angular dispersions for separate normal- and reverse-polarity data sets confirm that many older studies have not been adequately cleaned magnetically. Therefore, the use of the database for geomagnetic field modelling should be carried out with some caution. Using a VGP cut-off angle that varies with latitude, the best data set consists of 2636 results that show a smooth increase of VGP angular dispersion with latitude. Model G for palaeosecular variation, which is based on modelling of the antisymmetric (dipole) and symmetric (quadrupole) dynamo families, provides a good fit to these results.     

Palaeomagnetism of A Core Sediment From the Inland Sea, Japan (Seto Naikai)

Palaeomagnetic measurements on a giant core sample 20 cm in diameter and 7.38 m long collected from Mizushima-Nada, the Inland Sea. Japan (Seto Naikai) provide evidence of post-depositional magnetization. the geomagnetic secular variation from about 4000 to 8000 yr BP is characterized by a long period of westerly declination before 6600 yr BP. the maximum deflection is beyond 50°W during this period. the palaeomagnetic record further demonstrates that there is a hiatus over 3000 yr in sedimentation due 10 the sea-level change around 6.85 m from sea bottom, and that the lower limit of the sea-level around 8000 yr BP is 18.5 m beneath the present sea-level.     

Geomagnetic secular variations in Greece through the last 6000 years obtained from lake sediment studies

Summary. Geomagnetic declination and inclination secular variation curves have been constructed back to 6000 yr bp from palaeomagnetic measurements of sediments cored in Lakes Volvi, Trikhonis and Begoritis, Greece. They have been dated by palynological studies and by comparison with observatory and archaeomagnetic records. Radiocarbon ages are all affected by the'old carbon'factor. Cores from all three lakes contain a thin ash layer which we associate with the Pompeii eruption of Vesuvius in ad 79. Averaged curves of secular variations of declination and inclination have been computed for the whole study. The broader scale pattern of secular variations is basically similar to lake sediment records from the UK, Poland and Switzerland, though finer detail in the structure of some of the principal features is resolved.     

Geomagnetic field analysis - I. Stochastic inversion

Summary. Most of the Earth's magnetic field and its secular change originate in the core. Provided the mantle can be treated as an electrical insulator, stochastic inversion enables surface observations to be analysed for the core field. A priori information about the variation of the field at the core boundary leads to very stringent conditions at the Earth's surface. The field models are identical with those derived from the method of harmonic splines (Shure, Parker & Backus) provided the a priori information is specified appropriately.
The method is applied to secular variation data from 106 magnetic observatories. Model predictions for fields at the Earth's surface have error estimates associated with them that appear realistic. For plausible choices of a priori information the error of the field at the core is unbounded, but integrals over patches of the core surface can have finite errors. The hypothesis that magnetic fields are frozen to the core fluid implies that certain integrals of the secular variation vanish. This idea is tested by computing the integrals and their standard and maximum errors. Most of the integrals are within one standard deviation of zero, but those over the large patches to the north and south of the magnetic equator are many times their standard error, because of the dominating influence of the decaying dipole. All integrals are well within their maximum error, indicating that it will be possible to construct core fields, consistent with frozen flux, that satisfy the observations.     

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.     

Palaeosecular Variation of the Geomagnetic Field in the Aleutian Islands, Alaska

Sequences of late Tertiary to Quaternary lava flows at six Aleutian Island sites were sampled for palaeosecular variation (PSV) measurement using palaeomagnetic methods. Due to limitations in probable time span represented by each sequence of flows and statistical constraints, data from only two sites are considered in themselves to be reliable indicators of PSV. When all the data are combined to give greater statistical significance, an angular standard deviation of palaeofield directions of = 10.8° is obtained. This is compatible with other Pacific Ocean PSV data and points to the persistence of the present secular variation low in this region at least into Quaternary time.     

Annual variation of the geomagnetic field

Summary. Horizontal and vertical intensity data, obtained between 1957.0 and 1961.0 at 69 observatories, are analysed to determine the worldwide distribution of the annual variation of the geomagnetic field. Only data observed near local midnight are used, to avoid the small, but significant contamination from Sq. Over most of the world the variation is found to be small, with a clear dependence on latitude, but near the poles it is larger and more erratic. The non-polar variation is subjected to spherical harmonic analysis and separated into parts of internal and external origin. The polar variations are shown to be consistent with a north—south oscillation of the mean position of the auroral electrojets during the year. It is suggested that, with the exception of the polar effect, the annual variation is not due to ionospheric currents (as was hitherto believed), but results from an annual variation in the latitude of the ring current.     

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.