Publications

Export 22 results:
Sort by: Author Title Type [ Year  (Desc)]
2019
Korte, M, Brown MC, Gunnarson SR, Nilsson A, Panovska S, Wardinski I, Constable CG.  2019.  Refining Holocene geochronologies using palaeomagnetic records. Quaternary Geochronology. 50:47-74.   10.1016/j.quageo.2018.11.004   AbstractWebsite

The aperiodic nature of geomagnetic field variations, both in intensity and direction, can aid in dating archaeological artefacts, volcanic rocks, and sediment records that carry a palaeomagnetic signal. The success of palaeomagnetic dating relies upon our knowledge of past field variations at specific locations. Regional archaeo- and palaeomagnetic reference curves and predictions from global geomagnetic field models provide our best description of field variations through the Holocene. State-of-the-art palaeomagnetic laboratory practices and accurate independent age controls are prerequisites for deriving reliable reference curves and models from archaeological, volcanic, and sedimentary palaeomagnetic data. In this review paper we give an overview of these prerequisites and the available reference curves and models, discuss techniques for palaeomagnetic dating, and outline its limitations. In particular, palaeomagnetic dating on its own cannot give unique results, but rather serves to refine or confirm ages obtained by other methods. Owing to the non-uniform character of magnetic field variations in different regions, care is required when choosing a palaeomagnetic dating curve, so that the distance between the dating curve and the record to be dated is not too large. Accurate reporting and incorporation of new, independently dated archaeo- and palaeomagnetic results into databases will help to improve reference curves and global models for all regions on Earth.

2018
Korte, M, Constable CG.  2018.  Archeomagnetic intensity spikes: Global or regional geomagnetic field features? Frontiers in Earth Science. 6   10.3389/feart.2018.00017   AbstractWebsite

Variations of the geomagnetic field prior to direct observations are inferred from archeo- and paleomagnetic experiments. Seemingly unusual variations not seen in the present-day and historical field are of particular interest to constrain the full range of core dynamics. Recently, archeomagnetic intensity spikes, characterized by very high field values that appear to be associated with rapid secular variation rates, have been reported from several parts of the world. They were first noted in data from the Levant at around 900 BCE. A recent re-assessment of previous and new Levantine data, involving a rigorous quality assessment, interprets the observations as an extreme local geomagnetic high with at least two intensity spikes between the 11th and 8th centuries BCE. Subsequent reports of similar features from Asia, the Canary Islands and Texas raise the question of whether such features might be common occurrences, or whether they might even be part of a global magnetic field feature. Here we use spherical harmonic modeling to test two hypotheses: firstly, whether the Levantine and other potential spikes might be associated with higher dipole field intensity than shown by existing global field models around 1,000 BCE, and secondly, whether the observations from different parts of the world are compatible with a westward drifting intense flux patch. Our results suggest that the spikes originate from intense flux patches growing and decaying mostly in situ, combined with stronger and more variable dipole moment than shown by previous global field models. Axial dipole variations no more than 60% higher than observed in the present field, probably within the range of normal geodynamo behavior, seem sufficient to explain the observations.

2017
Panovska, S, Constable CG.  2017.  An activity index for geomagnetic paleosecular variation, excursions, and reversals. Geochemistry Geophysics Geosystems. 18:1366-1375.   10.1002/2016gc006668   AbstractWebsite

Magnetic indices provide quantitative measures of space weather phenomena that are widely used by researchers in geomagnetism. We introduce an index focused on the internally generated field that can be used to evaluate long term variations or climatology of modern and paleomagnetic secular variation, including geomagnetic excursions, polarity reversals, and changes in reversal rate. The paleosecular variation index, P-i, represents instantaneous or average deviation from a geocentric axial dipole field using normalized ratios of virtual geomagnetic pole colatitude and virtual dipole moment. The activity level of the index, sigma P-i, provides a measure of field stability through the temporal standard deviation of P-i. P-i can be calculated on a global grid from geomagnetic field models to reveal large scale geographic variations in field structure. It can be determined for individual time series, or averaged at local, regional, and global scales to detect long term changes in geomagnetic activity, identify excursions, and transitional field behavior. For recent field models, P-i ranges from less than 0.05 to 0.30. Conventional definitions for geomagnetic excursions are characterized by P-i exceeding 0.5. Strong field intensities are associated with low P-i unless they are accompanied by large deviations from axial dipole field directions. sigma P-i provides a measure of geomagnetic stability that is modulated by the level of PSV or frequency of excursional activity and reversal rate. We demonstrate uses of P-i for paleomagnetic observations and field models and show how it could be used to assess whether numerical simulations of the geodynamo exhibit Earth-like properties.

2016
Constable, C, Korte M, Panovska S.  2016.  Persistent high paleosecular variation activity in southern hemisphere for at least 10,000 years. Earth and Planetary Science Letters. 453:78-86.   10.1016/j.epsl.2016.08.015   AbstractWebsite

Direct observations of the geomagnetic field show that secular variation is strong in the Atlantic hemisphere, and comparatively reduced in the Pacific region. The dipole has been decaying since at least 1840 AD, driven by growth and migration of reverse flux patches in the southern hemisphere. We investigate whether anything like this modern pattern of geomagnetic secular variation persists and can be detected in global paleomagnetic field models. Synthesis of results from two new time-varying spherical harmonic models shows that geographically distinct geomagnetic secular variation extends to at least 10000 BP. The models use the same database but differ in methodology, leading to some regional differences in results. Consistent large-scale surface features include strong average fields in the northern hemisphere and weaker fields with greater overall variability in the south. Longitudinal structure is present, with weaker average fields in the western Pacific than in the east, and prominent negative inclination anomalies extending beneath Indonesia, across Africa and to Brazil, but weaker anomalies in the central Pacific. Marginally positive inclination anomalies occur west of the Americas. Paleosecular variation activity peaks at high southern latitudes, and there is a pattern of reduced activity at equatorial and mid-latitudes beneath the Pacific. Although the dipole has exhibited both growth and decay over the interval 0-10 000 BP, our results show that geomagnetic paleosecular variation is preferentially focused in similar geographic regions to secular variation seen in the modern field. (C) 2016 The Authors. Published by Elsevier B.V.

2015
Ziegler, LB, Constable CG.  2015.  Testing the geocentric axial dipole hypothesis using regional paleomagnetic intensity records from 0 to 300 ka. Earth and Planetary Science Letters. 423:48-56.   10.1016/j.epsl.2015.04.022   AbstractWebsite

Absolute and relative geomagnetic paleointensity records reveal variations in geomagnetic dipole strength, either via averaging time series of virtual axial dipole moments, or through formal inversion strategies like the penalized maximum likelihood (PML) method used for the PADM2M (Paleomagnetic Axial Dipole Moment for 0-2 Ma) model. However, departures from the most basic geocentric axial dipole (GAD) structure are obvious on centennial to millennial time scales, and paleomagnetic records from igneous rocks suggest small deviations persist on million year time scales. Spatial variations in heat flow at the core-mantle boundary (inferred from large low shear velocity provinces, LLSVPs) are widely suspected to influence both the average geomagnetic field and its regional secular variation. Long term departures from a GAD configuration should be visible from regional differences in paleointensity reconstructions. We use a PML method to construct time-varying models of regional axial dipole moment (RADMs) from a combined set of absolute and relative palebintensity data, and compare results from the last 300 kyr. RADMs are created from sediment records selected from specific latitude and longitude bands. We also test whether grouping records lying above each of the 2 major LLSVPs (centered on Africa and the Pacific) produce RADMs that are distinct from those above regions lacking anomalous seismic structure. Systematic differences appear in the various regional results. In the most recent part of the record regional differences are broadly similar to the Holocene, CALS10k.1b, time-varying geomagnetic field model spanning 0-10 ka. However, lack of Southern hemisphere records prevents direct confirmation of the hemispheric asymmetry present in CALS10k.1b in both average virtual axial dipole moment and its variability. As expected, the 300 kyr RADMs exhibit greater overall temporal field variability than is seen over 0-10 ka. Average RADM is higher in the Pacific and in Equatorial regions than in the Atlantic and in mid-high latitude northern hemisphere regions. Higher average RADMs are associated with lower overall field variability and less pronounced excursional signatures. Notably, the lower variability in the Pacific sector seen here (defined by either longitude band or LLSVP location) suggests that the modern low paleosecular variation there extends over at least the past few hundred thousand years. RADMs identified with LLSVPs show systematic deviations from the non-LLSVP group of records, with distinct characteristics for the African and Pacific provinces. The African LLSVP generates more pronounced RADM minima associated with geomagnetic excursions, and in general paleointensity decreases associated with excursions occur first in the Atlantic longitude sector and over the African LLSVP. (C) 2015 Elsevier B.V. All rights reserved.

Panovska, S, Korte M, Finlay CC, Constable CG.  2015.  Limitations in paleomagnetic data and modelling techniques and their impact on Holocene geomagnetic field models. Geophysical Journal International. 202:402-418.   10.1093/gji/ggv137   AbstractWebsite

Characterization of geomagnetic field behaviour on timescales of centuries to millennia is necessary to understand the mechanisms that sustain the geodynamo and drive its evolution. As Holocene paleomagnetic and archeomagnetic data have become more abundant, strategies for regularized inversion of modern field data have been adapted to produce numerous time-varying global field models. We evaluate the effectiveness of several approaches to inversion and data handling, by assessing both global and regional properties of the resulting models. Global Holocene field models cannot resolve Southern hemisphere regional field variations without the use of sediments. A standard data set is used to construct multiple models using two different strategies for relative paleointensity calibration and declination orientation and a selection of starting models in the inversion procedure. When data uncertainties are considered, the results are similar overall regardless of whether we use iterative calibration and reorientation, or co-estimation of the calibration and orientation parameters as part of the inversion procedure. In each case the quality of the starting model used for initial relative paleointensity calibration and declination orientation is crucial and must be based on the best absolute information available. Without adequate initial calibration the morphology of dipole moment variations can be recovered but its absolute value will be correlated with the initial intensity calibrations, an effect that might be mitigated by ensuring an appropriate fit to enough high quality absolute intensity data with low uncertainties. The declination reorientation mainly impacts regional field structure and in the presence of non-zonal fields will result in a non-zero local average. The importance of declination orientation is highlighted by inconsistencies in the West Pacific and Australian sediment records in CALS10k.1b model. Great care must also be taken to assess uncertainties associated with both paleomagnetic and age data and to evaluate the effects of poor data distribution. New consistently allocated uncertainty estimates for sediment paleomagnetic records highlight the importance of adequate uncertainties in the inversion process, as they determine the relative weighting among the data and overall normalized misfit levels which in turn influence the complexity of the inferred field models. Residual distributions suggest that the most appropriate misfit measure is the L-1 norm (minimum absolute deviation) rather than L-2 (least squares), but this seems to have relatively minor impact on the overall results. For future Holocene field modelling we see a need for comprehensive methods to assess uncertainty in individual archeomagnetic data so that these data or models derived from them can be used for reliable initial relative paleointensity calibration and declination orientation in sediments. More work will be needed to assess whether co-estimation or an iterative approach to inversion is more efficient overall. This would be facilitated by realistic and globally consistent data and age uncertainties from the paleomagnetic community.

Brown, MC, Donadini F, Korte M, Nilsson A, Korhonen K, Lodge A, Lengyel SN, Constable CG.  2015.  GEOMAGIA50.v3: 1. general structure and modifications to the archeological and volcanic database. Earth Planets and Space. 67:1-31.   10.1186/s40623-015-0232-0   AbstractWebsite

Background: GEOMAGIA50.v3 is a comprehensive online database providing access to published paleomagnetic, rock magnetic, and chronological data from a variety of materials that record Earth's magnetic field over the past 50 ka. Findings: Since its original release in 2006, the structure and function of the database have been updated and a significant number of data have been added. Notable modifications are the following: (1) the inclusion of additional intensity, directional and metadata from archeological and volcanic materials and an improved documentation of radiocarbon dates; (2) a new data model to accommodate paleomagnetic, rock magnetic, and chronological data from lake and marine sediments; (3) a refinement of the geographic constraints in the archeomagnetic/volcanic query allowing selection of particular locations; (4) more flexible methodological and statistical constraints in the archeomagnetic/volcanic query; (5) the calculation of predictions of the Holocene geomagnetic field from a series of time varying global field models; (6) searchable reference lists; and (7) an updated web interface. This paper describes general modifications to the database and specific aspects of the archeomagnetic and volcanic database. The reader is referred to a companion publication for a description of the sediment database. Conclusions: The archeomagnetic and volcanic part of GEOMAGIA50.v3 currently contains 14,645 data (declination, inclination, and paleointensity) from 461 studies published between 1959 and 2014. We review the paleomagnetic methods used to obtain these data and discuss applications of the data within the database. The database continues to expand as legacy data are added and new studies published. The web-based interface can be found at http://geomagia.gfz-potsdam.de

2014
Davies, CJ, Constable CG.  2014.  Insights from geodynamo simulations into long-term geomagnetic field behaviour. Earth and Planetary Science Letters. 404:238-249.   10.1016/j.epsl.2014.07.042   AbstractWebsite

Detailed knowledge of the long-term spatial configuration and temporal variability of the geomagnetic field is lacking because of insufficient data for times prior to 10 ka. We use realisations from suitable numerical simulations to investigate three important questions about stability of the geodynamo process: is the present field representative of the past field; does a time-averaged field actually exist; and, supposing it exists, how long is needed to define such a field. Numerical geodynamo simulations are initially selected to meet existing criteria for morphological similarity to the observed magnetic field. A further criterion is introduced to evaluate similarity of long-term temporal variations. Allowing for reasonable uncertainties in the observations, observed and synthetic axial dipole moment frequency spectra for time series of order a million years in length should be fit by the same power law model. This leads us to identify diffusion time as the appropriate time scaling for such comparisons. In almost all simulations, intervals considered to have good morphological agreement between synthetic and observed field are shorter than those of poor agreement. The time needed to obtain a converged estimate of the time-averaged field was found to be comparable to the length of the simulation, even in non-reversing models, suggesting that periods of stable polarity spanning many magnetic diffusion times are needed to obtain robust estimates of the mean dipole field. Long term field variations are almost entirely attributable to the axial dipole; nonzonal components converge to long-term average values on relatively short timescales (15-20 kyr). In all simulations, the time-averaged spatial power spectrum is characterised by a zigzag pattern as a function of spherical harmonic degree, with relatively higher power in odd degrees than in even degrees. We suggest that long-term spatial characteristics of the observed field may emerge on averaging times that are within reach for the next generation of global time-varying paleomagnetic field models. (C) 2014 Elsevier B.V. All rights reserved.

2011
Ziegler, LB, Constable CG, Johnson CL, Tauxe L.  2011.  PADM2M: a penalized maximum likelihood model of the 0-2 Ma palaeomagnetic axial dipole moment. Geophysical Journal International. 184:1069-1089.   10.1111/j.1365-246X.2010.04905.x   AbstractWebsite

P>We present a new time-varying model for palaeomagnetic axial dipole moment (PADM) for the past 2 Myr and compare it with earlier virtual axial dipole moment (VADM) reconstructions which have been based on stacking and averaging scaled relative palaeointensity records. The PADM is derived from both absolute and relative palaeointensity data and constructed using a new penalized maximum likelihood (PML) approach to recover a cubic B-spline representation of axial-dipole field variations on million year timescales. The PML method is explicitly intended to reduce bias in estimating the true axial dipole moment that arises in average VADM reconstructions. We apply the PML method to a set of 96 032 published data (1800 palaeointensities from igneous rocks, 3300 archaeointensities and 86 relative palaeointensity time-series of variable lengths and resolutions). Two models are discussed: PADM2Mp is a trial model based on a subset of the nine longest available sedimentary records; PADM2M uses a comprehensive data set (76 records, 81 446 data; 10 records were eliminated) and is our preferred model. PADM2M has a lower mean than existing VADM reconstructions but shows similarities in long-period variability. Some differences in timing, amplitude and resolution of certain features can be attributed to variations in age assignments. Others result from our more comprehensive data set and a reduction in bias attributable to PML modelling. PADM2M has an average axial dipole moment over 0-2 Ma of 5.3 x 1022 Am2 with a standard deviation of 1.5 x 1022 Am2. The Brunhes chron average (6.2 x 1022 Am2) is higher than for earlier epochs of Matuyama (4.8 x 1022 Am2), as seen in some previous studies. The power spectrum for our model agrees with previous estimates of the global palaeomagnetic power spectrum for frequencies up to about 102 Myr-1. We see no distinctive evidence in the power spectrum for orbital forcing of geodynamo behaviour.

Korte, M, Constable C, Donadini F, Holme R.  2011.  Reconstructing the Holocene geomagnetic field. Earth and Planetary Science Letters. 312:497-505.   10.1016/j.epsl.2011.10.031   AbstractWebsite

Knowledge of the Holocene evolution of Earth's magnetic field is important for understanding geodynamo processes in the core, is necessary for studying long-term solar-terrestrial relationships, and can provide useful age constraints for archeologicaland stratigraphic applications. Continuous time-varying global field models based on archeo- and paleomagnetic data are useful tools in this regard. We use a comprehensive data compilation and recently refined modelling strategies to produce CALS10k.1b, the first time-varying spherical harmonic geomagnetic field model spanning 10 ky. The model is an average obtained from bootstrap sampling to take account of uncertainties in magnetic components and ages in the data (and hence has version number 1b instead of 1). This model shows less spatial and temporal resolution than earlier versions for 0-3 ka, and particularly aims to provide a robust representation of the large-scale field at the core-mantle boundary (CMB). We discuss the geomagnetic dipole evolution and changes in Holocene magnetic field morphology at the CMB as shown by the new reconstruction. The results are compatible with earlier models (CALS3k.3 and CALS3k.4) for 0-3 ka, but reveal some clear deficiencies in the 0-7 ka CALS7K.2 model prior to 3 ka. CALS10k.1b is able to resolve mobile and structurally-evolving high latitude radial field flux lobes at the CMB in both hemispheres, as well as persistent non-zonal structure, in the 10 ky average. Contributions to the average field from time-varying structures in the equatorial Indonesian-Australian region are particularly striking. (C) 2011 Elsevier B.V. All rights reserved.

Amit, H, Korte M, Aubert J, Constable C, Hulot G.  2011.  The time-dependence of intense archeomagnetic flux patches. Journal of Geophysical Research-Solid Earth. 116   10.1029/2011jb008538   AbstractWebsite

The long-term temporal behavior of intense geomagnetic flux patches at the core-mantle boundary and the relation with lower mantle lateral heterogeneity are under debate. We apply an algorithm to detect centers of intense flux patches and track their time-evolution in a recent archeomagnetic field model in order to study the kinematics of such intense magnetic flux patches on millennial timescale. We find that most intense flux patches appear near the edge of the tangent cylinder. Quasi-stationary periods with small oscillations of patches occur more than drifts. Detailed comparison of the archeomagnetic patches' behavior with that seen in numerical dynamos with tomographic heat flux boundary conditions suggests that core-mantle thermal coupling could be the cause of a statistical preference for some longitudes on the long term, which does not exclude significant time spent away from the preferred longitudes. This could explain the roughly coincident locations of high-latitude patches in the historical geomagnetic field with that of the time-average paleomagnetic field together with the much weaker patches intensity in the latter. Alternating eastward and westward drifts are also observed. The drifts are more westward than eastward, especially in the southern hemisphere, indicating that the time-average zonal core flow may also be driven by core-mantle thermal coupling. An average patch lifetime of similar to 300 years is found, which we hypothesize may indicate the vortex lifetime in the outer core.

2010
Donadini, F, Korte M, Constable C.  2010.  Millennial Variations of the Geomagnetic Field: from Data Recovery to Field Reconstruction. Space Science Reviews. 155:219-246.   10.1007/s11214-010-9662-y   AbstractWebsite

Variations of the geomagnetic field over past millennia can be determined from archeomagnetic data and paleomagnetic sediment records. The resolution and validity of any field reconstruction depends on the reliability of such indirect measurements of past field values. Considerable effort is invested to ensure that the magnetic minerals carrying the ancient magnetization are good, if not ideal, recorders of the magnetic field. This is achieved by performing a wide array of rock magnetic and microscopy investigations, many of which are outlined here. In addition to data quality, the spatial and temporal distributions of archeomagnetic and sediment records play a significant role in the accuracy of past field reconstruction. Global field reconstructions enable studies of dynamic processes in Earth's core. They rely on data compilations which ideally include information about the quality of a measurement and provide a useful archive for selecting data with the best characteristics. There is, however, a trade off between the total number of reliable data and the geographic or temporal coverage. In this review we describe the various types of paleomagnetic recorders, and the kind of measurements that are performed to gather reliable geomagnetic field information. We show which modeling strategies are most suitable, and the main features of the field that can be derived from the resulting models. Finally, we discuss prospects for progress in this kind of research.

2009
Donadini, F, Korte M, Constable CG.  2009.  Geomagnetic field for 0-3 ka: 1. New data sets for global modeling. Geochemistry Geophysics Geosystems. 10   10.1029/2008gc002295   AbstractWebsite

Paleomagnetic and archeomagnetic records are used in both regional and global studies of Earth's magnetic field. We present a description and assessment of five newly compiled data sets, also used in the companion paper by Korte et al. (2009) to produce a series of time-varying spherical harmonic models of the geomagnetic field for the last 3000 years. Data are drawn from our compilation of lake sediment records and from the online database, GEOMAGIA50v2. The five selections are available from the EarthRef Digital Archive at http://earthref.org/cgi-bin/erda.cgi?n=944. Data are grouped according to the source of material, and we conducted separate assessments of reliability for archeomagnetic artifacts and lava flows (the ARCH3k_dat data set) and for sediments (SED3k_dat). The overall number of data is 55% greater than in previous compilations. Constrained data sets were selected using different criteria for each group. Winnowing of archeological data was based on uncertainties supplied by the original data providers. The lake sediment data assessment relied on preassigned age uncertainties and one or more of the following: comparisons with archeomagnetic data from the same region, regional consistency among several lakes, and consistency with global archeomagnetic models. We discuss relative merits of a larger unconstrained data set or a smaller (possibly) more reliable one. The constrained data sets eliminate a priori up to 35% of the available data in each case and rely on potentially subjective assessments of data quality. Given the limited data available our analyses indicate that iterative rejection of a small number (1-1.5%) of outlying data during global field modeling is a preferable approach. Specific regional comparisons among the models and data support the conclusion that Korte et al.'s outlier-free CALS3k.3 model based on all available measurements from sediments and archeological artifacts currently provides the best global representation of the 0-3 ka field; the ARCH3k.1 model provides a better fit to the denser European archeomagnetic data and may be better in that region.

Korte, M, Donadini F, Constable CG.  2009.  Geomagnetic field for 0-3 ka: 2. A new series of time-varying global models. Geochemistry Geophysics Geosystems. 10   10.1029/2008gc002297   AbstractWebsite

Steadily increasing numbers of archeomagnetic and paleomagnetic data for the Holocene have allowed development of temporally continuous global spherical harmonic models of the geomagnetic field extending present and historical global descriptions of magnetic field evolution. The current work uses various subsets of improved data compilations, details of which are given in a companion paper by Donadini et al. (2009), and minor modifications of standard modeling strategies (using temporally and spatially regularized inversion of the data and cubic spline parameterizations for temporal variations) to produce five models with enhanced spatial and temporal resolution for 0-3 ka. Spurious end effects present in earlier models are eliminated by enforcing large-scale agreement with the gufm1 historical model for 1650-1990 A.D. and by extending the model range to accommodate data older than 3 ka. Age errors are not considered as a contribution to data uncertainties but are included along with data uncertainties in an investigation of statistical uncertainty estimates for the models using parametric bootstrap resampling techniques. We find common features but also significant differences among the various models, indicating intrinsic uncertainties in global models based on the currently available Holocene data. Model CALS3k.3 based on all available archeomagnetic and sediment data, without a priori quality selection, currently constitutes the best global representation of the past field. The new models have slightly higher dipole moments than our previous models. Virtual axial dipole moments (VADMs) calculated directly from the data are in good agreement with all corresponding model predictions of VADMs. These are always higher than the spherical harmonic dipole moment, indicating the limitations of using VADMs as a measure of geomagnetic dipole moments.

2008
Genevey, A, Gallet Y, Constable CG, Korte M, Hulot G.  2008.  ArcheoInt: An upgraded compilation of geomagnetic field intensity data for the past ten millennia and its application to the recovery of the past dipole moment. Geochemistry Geophysics Geosystems. 9   10.1029/2007gc001881   AbstractWebsite

This paper presents a compilation of intensity data covering the past 10 millennia (ArcheoInt). This compilation, which upgrades the one of Korte et al. (2005), contains 3648 data and incorporates additional intensity and directional data sets. A large majority of these data (similar to 87%) were acquired on archeological artifacts, and the remaining similar to 13% correspond to data obtained from volcanic products. The present compilation also includes important metadata for evaluating the intensity data quality and providing a foundation to guide improved selection criteria. We show that similar to 50% of the data set fulfill reasonable reliability standards which take into account the anisotropic nature of most studied objects (potsherds), the stability of the magnetization, and the data dispersion. The temporal and geographical distributions of this sub-data set are similar to those of the main data set, with similar to 72% of the data dated from the past three millennia and similar to 76% obtained from western Eurasia. Approximately half of the selected intensity data are associated with at least an inclination value. To constrain the axial and full dipole evolution over the past three millennia requires that we avoid any overrepresentation of the western Eurasian data. We introduce a first-order regional weighting scheme based on the definition of eight widely distributed regions of 30 degrees width within which the selected data are numerous enough. The regional curves of virtual axial dipole moments (VADM) and of mixed VADM-virtual dipole moments (VDM) averaged over sliding windows of 200 years and 500 years testify for strong contributions from either equatorial dipole or nondipole components. The computation of global VADM and mixed VADM/VDM variation curves, assuming an equal weight for each region, yields a dipole evolution marked by a distinct minimum around 0 B.C./A.D. followed by a maximum around the third-fourth century A. D. A second minimum is present around the eighth century A. D. This variation pattern is compatible with the one deduced from earlier, more sophisticated analysis based on the inversion of both intensity and directional data. In particular, there is a good agreement among all VADMs and dipole moment estimates over the historical period, which further strengthens the validity of our weighting scheme.

2005
Constable, C, Johnson C.  2005.  A paleomagnetic power spectrum. Physics of the Earth and Planetary Interiors. 153:61-73.   10.1016/j.pepi.2005.03.015   AbstractWebsite

We construct a power spectrum of geomagnetic dipole moment variations or their proxies that spans the period range from some tens of million down to about 100 years. Empirical estimates of the spectrum are derived from the magnetostratigraphic time scale, from marine sediment relative paleointensity records, and from a time varying paleomagnetic field model for the past 7 kyr. The spectrum has the most power at long periods, reflecting the influence of geomagnetic reversals and in general decreases with increasing frequency (decreasing period). The empirical spectrum is compared with predictions from simple models. Discrepancies between the observed and predicted spectra are discussed in the context of: (i) changes in reversal rate, (ii) overall average reversal rate, (iii) cryptochrons, (iv) the time taken for a reversal to occur, and (v) long term paleosecular variations and average estimates of the field strength and variance from other sources. (c) 2005 Elsevier B.V. All rights reserved.

Korte, M, Genevey A, Constable CG, Frank U, Schnepp E.  2005.  Continuous geomagnetic field models for the past 7 millennia: 1. A new global data compilation. Geochemistry Geophysics Geosystems. 6   10.1029/2004gc000800   AbstractWebsite

A global data set of archeomagnetic and paleomagnetic data covering the past 7000 years has been compiled. It consists of 16,085 results of inclination, 13,080 of declination, and 3188 of intensity for the time span 5000 BC to 1950 AD. Declination and inclination data come partly from existing databases and partly from original literature. A new global compilation of intensity data for the millennial scale is included. Data and dating uncertainties are discussed as we attempted to obtain an internally coherent data set. The global distribution of the data is very inhomogeneous in both time and space. All the data are compared to predictions from the previous 3000 year global model, CALS3K.1. This collection of data will be useful for global secular variation studies and geomagnetic field modeling, although southern hemisphere data are still underrepresented. In particular, we will use it in a further study to update and extend the existing global model, CALS3K.1. The huge increase in data compared to the previous compilation will result in significant changes from current models. As we might have missed some suitable data, we encourage the reader to notify us about any data that have not been included yet and might fit in, as improving our global millennial scale models remains our aim for the future.

2003
Love, JJ, Constable CG.  2003.  Gaussian statistics for palaeomagnetic vectors. Geophysical Journal International. 152:515-565.   10.1046/j.1365-246X.2003.01858.x   AbstractWebsite

With the aim of treating the statistics of palaeomagnetic directions and intensities jointly and consistently, we represent the mean and the variance of palaeomagnetic vectors, at a particular site and of a particular polarity, by a probability density function in a Cartesian three-space of orthogonal magnetic-field components consisting of a single (unimodal) non-zero mean, spherically-symmetrical (isotropic) Gaussian function. For palaeomagnetic data of mixed polarities, we consider a bimodal distribution consisting of a pair of such symmetrical Gaussian functions, with equal, but opposite, means and equal variances. For both the Gaussian and bi-Gaussian distributions, and in the spherical three-space of intensity, inclination, and declination, we obtain analytical expressions for the marginal density functions, the cumulative distributions, and the expected values and variances for each spherical coordinate (including the angle with respect to the axis of symmetry of the distributions). The mathematical expressions for the intensity and off-axis angle are closed-form and especially manageable, with the intensity distribution being Rayleigh-Rician. In the limit of small relative vectorial dispersion, the Gaussian (bi-Gaussian) directional distribution approaches a Fisher (Bingham) distribution and the intensity distribution approaches a normal distribution. In the opposite limit of large relative vectorial dispersion, the directional distributions approach a spherically-uniform distribution and the intensity distribution approaches a Maxwell distribution. We quantify biases in estimating the properties of the vector field resulting from the use of simple arithmetic averages, such as estimates of the intensity or the inclination of the mean vector, or the variances of these quantities. With the statistical framework developed here and using the maximum-likelihood method, which gives unbiased estimates in the limit of large data numbers, we demonstrate how to formulate the inverse problem, and how to estimate the mean and variance of the magnetic vector field, even when the data consist of mixed combinations of directions and intensities. We examine palaeomagnetic secular-variation data from Hawaii and Reunion, and although these two sites are on almost opposite latitudes, we find significant differences in the mean vector and differences in the local vectorial variances, with the Hawaiian data being particularly anisotropic. These observations are inconsistent with a description of the mean field as being a simple geocentric axial dipole and with secular variation being statistically symmetrical with respect to reflection through the equatorial plane. Finally, our analysis of palaeomagnetic acquisition data from the 1960 Kilauea flow in Hawaii and the Holocene Xitle flow in Mexico, is consistent with the widely held suspicion that directional data are more accurate than intensity data.

1999
Constable, CG, Johnson CL.  1999.  Anisotropic paleosecular variation models: implications for geomagnetic field observables. Physics of the Earth and Planetary Interiors. 115:35-51.   10.1016/s0031-9201(99)00065-5   AbstractWebsite

We present a family of statistical models for paleosecular variation (PSV) of the geomagnetic field that are compatible with paleodirectional and paleointensity variations in lava flows sampling the last 5 Ma, and explore what paleomagnetic observables might be used to discriminate among the various family members. We distinguish statistical models with axial anisotropy, which provide a suitable description for an earth with homogeneous boundary conditions at the core-mantle interface from those with more general anisotropy corresponding to geographically heterogeneous boundary conditions. The models revise and extend earlier ones, which are themselves descendants of CP88, devised by Constable and Parker [Constable, C.G., Parker, R.L., 1988. Statistics of the geomagnetic secular variation for the past 5 m.y. J. Geophys, Res. 93, 11569-11581]. In CP88, secular variation is described by statistical variability of each Gauss coefficient in a spherical harmonic description of the geomagnetic field, with each coefficient treated as a normally distributed random variable: the Gauss coefficients of the non-dipole part of the field exhibit isotropic variability, and the variances are derived from the present field spatial power spectrum. The dipole terms have a special status in CP88, with a non-zero mean for the axial-dipole, and lower variance than predicted from the spatial power spectrum. All non-dipole terms have zero mean except the axial-quadrupole. CP88 is untenable for two reasons: it fails to predict the observed geographic dependence of directional variability in the magnetic field, and it grossly underpredicts the variance in paleointensity data. The new models incorporate large variance in the axial-dipole, and in the non-axial-quadrupole Gauss coefficients, g1/2: and h1/2:. The resulting variance in paleomagnetic observables depends only on latitude (zonal models), unless the variance in h1/2: is different from that in g1/2 (non-zonal models). Non-zonal (longitudinal) variations in PSV, such as the flux lobes seen in the historical magnetic field, are simulated using the non-zonal models. Both the zonal and non-zonal models fit summary statistics of the present dataset. We investigate the influence of persistent non-zonal influences in PSV on various paleomagnetic observables. It is shown that virtual geomagnetic pole (VGP) dispersion is rather insensitive to longitudinal variations in structure of PSV, and that inclination dispersion has the potential to be more informative given the right site distribution. There is also the possibility of using paleointensity and geographic variations in the frequency of occurrence of excursional directions to identify appropriate PSV models. (C) 1999 Elsevier Science B.V. All rights reserved.

1997
Johnson, CL, Constable CG.  1997.  The time-averaged geomagnetic field: global and regional biases for 0-5 Ma. Geophysical Journal International. 131:643-+.   10.1111/j.1365-246X.1997.tb06604.x   AbstractWebsite

Palaeodirectional data from lava flows and marine sediments provide information about the long-term structure and variability in the geomagnetic held. We present a detailed analysis of the internal consistency and reliability of global compilations of sediment and lava-flow data. Time-averaged field models are constructed for normal and reverse polarity periods for the past 5 Ma, using the combined data sets. Non-zonal models are required to satisfy the lava-flow data, but not those from sediments alone. This is in part because the sediment data are much noisier than those from lavas, but is also a consequence of the site distributions and the way that inclination data sample the geomagnetic field generated in the Earth's core. Different average held configurations for normal and reverse polarity periods are consistent with the palaeomagnetic directions; however, the differences are insignificant relative to the uncertainty in the average field models. Thus previous inferences of non-antipodal normal and reverse polarity field geometries will need to be re-examined using recently collected high-quality palaeomagnetic data. Our new models indicate that current global sediment and lava-flow data sets combined do not permit the unambiguous detection of northern hemisphere flux lobes in the 0-5 Ma time-averaged field, highlighting the need for the collection of additional high-latitude palaeomagnetic data. Anomalous time-averaged held structure is seen in the Pacific hemisphere centred just south of Hawaii. The location of the anomaly coincides with heterogeneities in the lower mantle inferred from seismological data. The seismic observations can be partly explained by lateral temperature variations; however, they also suggest the presence of lateral compositional variations and/or the presence of partial melt. The role of such heterogeneities in influencing the geomagnetic held observed at the Earth's surface remains an unresolved issue, requiring higher-resolution time-averaged geomagnetic field models, along with the integration of future results from seismology, mineral physics and numerical simulations.

1995
Johnson, CL, Constable CG.  1995.  The Time-Averaged Geomagnetic-Field As Recorded By Lava Flows Over The Past 5 Million-Years. Geophysical Journal International. 122:489-519.   10.1111/j.1365-246X.1995.tb07010.x   AbstractWebsite

A recently compiled lava flow data base spanning the last 5 million years is used to investigate properties of the time-averaged geomagnetic field. More than 90 per cent of the power in the palaeofield can be accounted for by a geocentric axial dipole; however, there are significant second-order structures in the held. Declination and inclination anomalies for the new data base indicate that the main second-order signal is the 'far-sided' effect, and there is also evidence for non-zonal structure. VGP (virtual geomagnetic pole) latitude distributions indicate that, over the last 5 million years, normal and reverse polarity morphologies are different, and that any changes in the normal polarity field morphology are undetectable, given the present data distribution. Regularized non-linear inversions of the palaeomagnetic directions support all these observations. We test the hypothesis that zonal models for the time-averaged field are adequate to describe the data and find that they are not. Non-zonal models are needed to fit the data to within the required tolerance level. Normal and reverse polarity held models obtained are significantly different. Field models obtained for the Brunhes epoch data alone are much smoother than those obtained from combining an the normal polarity data; simulations indicate that these differences can be explained by the less extensive data distribution for the Brunhes epoch. The field model for all of the normal polarity data (LN1) contains features observed in the historical field maps, although the details differ. LN1 suggests that, although the two northern hemisphere flux lobes observed in the historical field are stationary to a first-order approximation, they do show changes in position and amplitude. A. third, less pronounced flux lobe is observed in LN1 over central Europe. The lack of structure ih the southern hemisphere is due in part to the paucity of data. Jackknife estimates of the field models for different subsets of the data suggest that a few sites contribute significant structure to the final field models. More conservative estimates of the time-averaged field morphology are obtained by removing these sites.

1993
Constable, CG, Parker RL, Stark PB.  1993.  Geomagnetic-Field Models Incorporating Frozen-Flux Constraints. Geophysical Journal International. 113:419-433.   10.1111/j.1365-246X.1993.tb00897.x   AbstractWebsite

Techniques for modelling the geomagnetic field at the surface of Earth's core often penalize contributions at high spherical harmonic degrees to reduce the effect of mapping crustal fields into the resulting field model at the core-mantle boundary (CMB). Ambiguity in separating the observed field into crustal and core contributions makes it difficult to assign error bounds to core field models, and this makes it hard to test hypotheses that involve pointwise values of the core field. The frozen-flux hypothesis, namely that convective terms dominate diffusive terms in the magnetic-induction equation, requires that the magnetic flux through every patch on the core surrounded by a zero contour of the radial magnetic field remains constant, although the shapes, areas and locations (but not the topology) of these patches may change with time. Field models exactly satisfying the conditions necessary for the hypothesis have not yet been constructed for the early part of this century. We show that such models must exist, so testing the frozen-flux hypothesis becomes the question of whether the models satisfying it are geophysically unsatisfactory on other grounds, for example because they are implausibly rough or complicated. We introduce an algorithm to construct plausible fleld models satisfying the hypothesis, and present such models for epochs 1945.5 and 1980. Our algorithm is based on a new parametrization of the field in terms of its radial component B(r) at the CMB. The model consists of values of B(r) at a finite set of points on the CMB, together with a rule for interpolating the values to other points. The interpolation rule takes the specified points to be the vertices of a spherical triangle tessellation of the CMB, with B(r) varying linearly in the gnomonic projections of the spherical triangles onto planar triangles in the planes tangent to the centroids of the spherical triangles. This parametrization of B(r) provides a direct means of constraining the integral invariants required by the frozen-flux hypothesis. Using this parametrization, we have constructed field models satisfying the frozen-flux hypothesis for epochs 1945.5 and 1980, while fitting observatory and survey data for 1945.5 and Magsat data for 1980. We use the better constrained 1980 CMB field model as a reference for 1945.5: we minimize the departure of the 1945.5 CMB field model from a regularized 1980 CMB field model, while constraining the 1945.5 model to have the same null-flux curves and flux through those curves as the 1980 model. The locations, areas and shapes of the curves are allowed to change. The resulting 1945.5 CMB field model is nearly as smooth as that for 1980, fits the data adequately, and satisfies the conditions necessary for the frozen-flux hypothesis.