Search

Archeomagnetism focuses on exploring high-resolution variations of the geomagnetic field over hundreds to thousands of years. In this study, we carried out a comprehensive study of chronology, absolute and relative paleointensity on a late Neolithic site in central China. Ages of the samples are constrained to be ~3,500–3,000 BCE, a period when available paleointensity data are sparse. We present a total of 64 high-fidelity absolute paleointensities, demonstrating the field varied quickly from ~55 to ~90 ZAm2 between ~3,500–3,000 BCE. Our results record a new archeomagnetic jerk around 3,300 BCE, which is probably non-dipolar origin. The new results provide robust constraints on global geomagnetic models. We calculated a revised Chinese archeointensity reference curve for future application. The variations of absolute and relative paleointensity versus depth show good consistency, reinforcing the reliability of our results. This new attempt of combining absolute and relative paleointenstiy provides a useful tool for future archeomagnetic research.

This paper addresses one of the critical questions of scientific inquiry: How do we know when a given data set is representative of the phenomenon being examined? For paleomagnetists, the question is often whether a particular data set sufficiently averaged paleosecular variation (PSV). To this aim, we updated an existing PSV data set that now comprises 2,441 site mean directions from 94 individual studies (PSV10‐24). Minimal filtering for data quality resulted in 1,619 sites from 90 publications. Fitting PSV10‐24 with two newly defined parameters as well as two existing ones form the basis of a Giant Gaussian Process field model (THG24) consistent with the data. Drawing directions from THG24 yields directional distributions predicted for a given latitude allowing a comparison between empirical distributions and the cumulative distribution function generated by the model. This tests whether the observed data adequately averaged out PSV according to THG24. We applied these tests to five data sets from Large Igneous Provinces from the last billion years and find that they are consistent with the THG24 model as well. Sedimentary data sets that may have experienced inclination shallowing can be corrected using an (un)flattening factor that yields directions satisfying THG24 in a newly‐defined, four‐parameter space. This approach builds on the Elongation‐Inclination (E/I) method of Tauxe and Kent (2004), https://doi.org/10.1029/145gm08, so the approach introduced here is called SVEI. We show examples of the use of SVEI and explain how to use this newly developed Python code that is publicly available in the PmagPy GitHub repository. Plain Language Summary: Paleomagnetic vectors recorded by rocks and archeological materials yield information on the structure of the magnetic field through Earth's history. The geomagnetic field mostly resembles a geocentric dipole aligned with the spin axis, but the directions at any given time and place generally deviate from this simple model. Data sets produced over several decades help define global magnetic field behavior. Here, we updated a compilation of paleomagnetic data and used it to establish a new field model representative of the last 10 million years. Prior or emerging data sets can be tested against the model to see whether the two agree in terms of field structure and variability. We tested data sets as old as 1.1 billion years and found them compatible with the model. Moreover, this model can also assess, and potentially correct data from sedimentary rocks that may have suffered from inclination shallowing. Although other correction methods exist, our approach employs a more complete description of the geometry of directional data, thereby allowing a more quantitative comparison of empirical and predicted distributions. Corrected inclinations and their uncertainties define paleolatitudes more accurately, key for plate tectonic and paleoclimatic reconstructions. Key Points: Update of the paleosecular variation database from lava flows (PSV10‐24) for the last 10 MaDefinition of a new Giant Gaussian Process (GGP) model (THG24), compatible with PSV10‐24 matching VGP scatter including data sets as old as 1.1 GaDescription of the SVEI method (an update of E/I) to assess secular variation records and to correct for inclination shallowing in sediments [ABSTRACT FROM AUTHOR]

The Magnetics Information Consortium (MagIC) database provides an archive with a flexible data model for paleomagnetic and rock magnetic data. The PmagPy software package is a cross-platform and open-source set of tools written in Python for the analysis of paleomagnetic data that serves as one interface to MagIC, accommodating various levels of user expertise. PmagPy facilitates thorough documentation of sampling, measurements, data sets, visualization, and interpretation of paleomagnetic and rock magnetic experimental data. Although not the only route into the MagIC database, PmagPy makes preparation of newly published data sets for contribution to MagIC as a byproduct of normal data analysis and allows manipulation as well as reanalysis of data sets downloaded from MagIC with a single software package. The graphical user interface (GUI), Pmag GUI enables use of much of PmagPy's functionality, but the full capabilities of PmagPy extend well beyond that. Over 400 programs and functions can be called from the command line interface mode, or from within the interactive Jupyter notebooks. Use of PmagPy within a notebook allows for documentation of the workflow from the laboratory to the production of each published figure or data table, making research results fully reproducible. The PmagPy design and its development using GitHub accommodates extensions to its capabilities through development of new tools by the user community. Here we describe the PmagPy software package and illustrate the power of data discovery and reuse through a reanalysis of published paleointensity data which illustrates how the effectiveness of selection criteria can be tested.

Studies of early Paleogene climate suffer from the scarcity of well-dated sedimentary records from the southern Pacific Ocean, the largest ocean basin during this time. We present a new magnetostratigraphic record from marine sediments that outcrop along the mid-Waipara River, South Island, New Zealand. Fully oriented samples for paleomagnetic analyses were collected along 45 m of stratigraphic section, which encompasses magnetic polarity Chrons from C23n to C21n (~51.5-47 Ma). These results are integrated with foraminiferal, calcareous nannofossil, and dinoflagellate cyst (dinocyst) biostratigraphy from samples collected in three different expeditions along a total of ~80 m of section. Biostratigraphic data indicates relatively continuous sedimentation from the lower Waipawan to the upper Heretaungan New Zealand stages (i.e., lower Ypresian to lower Lutetian, 55.5 to 46 Ma). We provide the first magnetostratigraphically-calibrated age of 48.88 Ma for the base of the Heretaungan New Zealand stage (latest early Eocene). To improve the correlation of the climate record in this section with other Southern Ocean records, we reviewed the magnetostratigraphy of Ocean Drilling Program (ODP) Site 1172 (East Tasman Plateau) and Integrated Ocean Drilling Program (IODP) Site U1356 (Wilkes Land Margin, Antarctica). A paleomagnetic study of discrete samples could not confirm any reliable magnetic polarity reversals in the early Eocene at Site 1172. We use the robust magneto-biochronology of a succession of dinocyst bioevents that are common to mid-Waipara, Site 1172, and Site U1356 to assist correlation between the three records. A new integrated chronology offers new insights into the nature and completeness of the southern high-latitude climate histories derived from these sites.

Hematite-bearing red beds are renowned for their chemical remanent magnetization (CRM). If the CRM was acquired substantially later than the sediment was formed, this severely compromises paleomagnetic records. To improve our interpretation of the natural remanent magnetization, the intricacies of the CRM acquisition process must be understood. Here, we contribute to this issue by synthesizing hematite under controlled 'Earth-like' field conditions (≲100 μT). CRM was imparted in 90 oriented samples with varying inclinations. The final synthesis product appeared to be dominated by hematite with traces of ferrimagnetic iron oxides. When the magnetic field intensity is ≳40 μT, the CRM records the field direction faithfully. However, for field intensities ≲40 μT, the CRM direction may deviate considerably from that of the applied field during synthesis. The CRM intensity normalized by the isothermal remanent magnetization (CRM/IRM@2.5 T) increases linearly with the intensity of growth field, implying that CRM could potentially be useful for relative paleointensity studies if hematite particles of chemical origins have consistent properties. CRM in hematite has a distributed unblocking temperature spectrum from ~200 to ~650 °C, while hematite with a depositional remanent magnetization (DRM) has a more confined spectrum from ~600to680°C because it is usually coarser-grained and more stoichiometric. Therefore, the thermal decay curves of CRM with their concave shape are notably different from their DRM counterparts which are convex. These differences together are suggested to be a potential discriminator of CRM from DRM carried by hematite in natural red beds, and of significance for the interpretation of paleomagnetic studies on red beds.

The Arabian Sea in the northern Indian Ocean pre-serves regional sedimentary records of rifting, tectonic subsidence, and paleoceanographic history, and also provides archives of long-term erosion of the Himalaya since the start of collision between In?dia and Eurasia. Investigations reveal that drilling in this region can provide erosion records through analyses of the sediment cores, along with providing age control for the regional seismic stratigraphy. It is only by quantifying the volume of sedi?ment deposited in the fan that researchers can mass balance the volume of bedrock eroded from the mountains, constrained by thermochronology, with the volume of eroded rock deposited in the offshore and in the foreland basin.

©2015. American Geophysical Union. All Rights Reserved. The Earth's magnetic field is assumed to be a geocentric axial dipole (GAD) when averaged over sufficient time (105-106years). Recent investigations of global paleosecular variation and time-averaged field behavior on million year timescales generally support a predominantly dipole field in the Northern Hemisphere, but unique field structures at high southern latitudes suggest the presence of a substantial g20 quadrupolar component. Average paleointensity results from Antarctica are approximately half the value predicted by a GAD field; this behavior has not been sufficiently investigated because there is a paucity of absolute paleointensity data from the high latitudes of the Arctic and Antarctic, so no adequate comparisons have been made between the two regions. We collected glassy volcanic material from 129 subaerial and subglacial volcanic units in Iceland in order to provide a suitable intensity data set at high northern latitudes. Forty-four sites met our very strict specimen and site level selection criteria. Four Holocene sites have a median intensity value of 55.8 ± 15.6μT (virtual axial dipole moment = 78.1 ± 22.0 ZAm2), consistent with the present-day field. Thirty-seven sites are between 11 ka and 3.35 Ma with a median intensity of 33.1 ± 8.3 μT (47.0 ± 11.6 ZAm2). This median intensity is indistinguishable from some long-term global field strength estimates. Reevaluation of existing high-latitude data suggests a general agreement with our Iceland results, but there are still too few Antarctic sites to adequately compare Arctic and Antarctic field behaviors

©2015. American Geophysical Union. All Rights Reserved. We have carried out an archaeomagnetic study on a late Neolithic locality (Liujiazhai) in Sichuan, southwestern China. We pull together various dating techniques, including radiocarbon analysis, optically stimulated luminescence dating, stratigraphic information as well as archaeological and archaeomagnetic estimations, to constrain the age of the studied samples. Rock magnetic results indicate thermally stable fine-grained magnetite or titanomagnetite as the dominant magnetic carriers. More than half of the specimens (141/246) in the paleointensity experiment pass the selection criteria and are considered to record robust intensity values. The virtual axial dipole moments range from approximately (2.8 to 7.8) × 10 < sup > 22 < /sup > Am < sup > 2 < /sup > with an average of 5.9 × 10 < sup > 22 < /sup > Am < sup > 2 < /sup > , indicating that the geomagnetic intensity around 3000 before the Common Era (B.C.E.) is overall lower than the present field intensity (9.8 × 10 < sup > 22 < /sup > Am < sup > 2 < /sup > ) of this area. The new results from Liujiazhai are generally consistent with the published data of similar age but deviate from the only available model of CALS10k.1b at certain time periods, making them important for future improvements of the model. Those data are significant for constraining the variation of geomagnetic field intensity between ∼3100 and 2600 B.C.E. and improving the regional model of eastern Asia.

Published paleointensity estimates derived from lavas extruded in known fields are highly variable and rarely recover the expected field strength within an accuracy of better than 10%. Inconsistent results on modern volcanic rocks lend even greater uncertainty to intensity experiments performed on lava flows emplaced during periods of unknown geomagnetic field strength. The majority of published paleointensity data are collected from the slowly cooled, massive centers of lava flows, where the magnetic grains are more likely to be multi-domain and produce non-ideal experimental results. Glassy volcanic material (found on subaerial lava flow tops and in sub-aqueous and subglacial environments), however is rapidly cooled, and therefore most likely of all volcanic materials to behave as single-domain particles demanded by Néel theory. We present a new paleointensity study of historic and modern Hawaiian lavas and test the viability of subaerially emplaced glassy basaltic material as an accurate recorder of magnetic field intensity. Six of eight lava flows sampled on the Big Island of Hawaii (1843, 1859, 1935, 1950, 1960, 1990 C.E.) produce well behaved Arai plots and recover an average intensity to within 2.7. μT of the expected field strength or better than 8% accuracy. We apply very strict selection criteria, including a minimum of three specimens per site, to prevent extraneous field estimates from affecting the final results. Individual volcanic glass results from the 1960 C.E. lava flow have a much lower variance than published data from the same volcanic unit. Glassy materials should therefore be collected wherever possible as they allow recovery of geomagnetic field strength with unprecedented accuracy.

© 2015. American Geophysical Union. All Rights Reserved. Geomagnetic models based on direct observations since the 1830s show that the averaged relative change in field intensity on Earth's surface over the past 170 years is less than 4.8% per decade. It is unknown if these rates represent the typical behavior of secular variations due to insufficient temporal resolution of archaeomagnetic records from earlier periods. To address this question, we investigate two ancient slag mounds in Cyprus-Skouriotissa Vouppes (SU1, fourth to fifth centuries CE, 21 m in height), and Mitsero Kokkinoyia (MK1, seventh to fifth centuries BCE, 8 m in height). The mounds are multilayered sequences of slag and charcoals that accumulated near ancient copper production sites. We modeled the age-height relation of the mounds using radiocarbon dates, and estimated paleointensities using Thellier-type IZZI experiments with additional anisotropy, cooling rate, and nonlinear TRM assessments. To screen out ambiguous paleointensity interpretations, we applied strict selection criteria at the specimen/sample levels. To ensure objectivity, consistency, and robust error estimation, we employed an automatic interpretation technique and put the data available in the MagIC database. The analyses yielded two independent subcentury-scale paleointensity time series. The MK1 data indicate relatively stable field at the time the mound accumulated. In contrast, the SU1 data demonstrate changes that are comparable in magnitude to the fastest changes inferred from geomagnetic models. We suggest that fast changes observed in the published archaeomagnetic data from the Levant are driven by two longitudinally paired regions, the Middle East and South Africa, that show unusual activity in geomagnetic models.

© 2014 Macmillan Publishers Limited. All rights reserved. The Pliocene and Early Pleistocene, between 5.3 and 0.8 million years ago, span a transition from a global climate state that was 2-3 °C warmer than present with limited ice sheets in the Northern Hemisphere to one that was characterized by continental-scale glaciations at both poles. Growth and decay of these ice sheets was paced by variations in the Earth's orbit around the Sun. However, the nature of the influence of orbital forcing on the ice sheets is unclear, particularly in light of the absence of a strong 20,000-year precession signal in geologic records of global ice volume and sea level. Here we present a record of the rate of accumulation of iceberg-rafted debris oshore from the East Antarctic ice sheet, adjacent to the Wilkes Subglacial Basin, between 4.3 and 2.2 million years ago. We infer that maximum iceberg debris accumulation is associated with the enhanced calving of icebergs during ice-sheet margin retreat. In the warmer part of the record, between 4.3 and 3.5 million years ago, spectral analyses show a dominant periodicity of about 40,000 years. Subsequently, the powers of the 100,000-year and 20,000-year signals strengthen. We suggest that, as the Southern Ocean cooled between 3.5 and 2.5 million years ago, the development of a perennial sea-ice field limited the oceanic forcing of the ice sheet. After this threshold was crossed, substantial retreat of the East Antarctic ice sheet occurred only during austral summer insolation maxima, as controlled by the precession cycle.

In this study, we have carried out paleointensity experiments on 918 specimens spanning the last ~7 kyr, including pottery fragments, baked clay and slag, collected from Shandong, Liaoning, Zhejiang and Hebei Provinces in China. Approximately half of the specimens yielded results that passed strict data selection criteria and give high-fidelity paleointensities. The virtual axial dipole moments (VADMs) of our sites range from ~2×1022to ~13×1022Am2. At ~2250 BCE our results suggest a paleointensity low of ~2×1022Am2, which increases to a high of ~13×1022Am2by ~1300 BCE. This rapid (less than 1000 yrs) six-fold change in the paleointensity may have important implications for the dynamics of core flow at this time. Our data from the last ~3 kyr are generally in good agreement with the ARCH3k.1 model, but deviate significantly at certain time periods from the CALS3k.4 and CALS10k.1b model, which is likely due to differences in the data used to constrain these models. At ages older than ~3 ka, where only the CALS10k.1b model is available for comparison, our data deviate significantly from the model. Combining our new results with the published data from China and Japan, we provide greatly improved constraints for the regional model of Eastern Asia. When comparing the variations of geomagnetic field in three global representative areas of Eastern Asia, the Middle East and Southern Europe, a common general trend of sinusoidal variations since ~8 ka is shown, likely dominated by the dipole component. However, significant disparities are revealed as well, which we attribute to non-dipolar components caused by movement of magnetic flux patches at the core-mantle boundary. © 2014 Elsevier B.V.

The selection of paleointensity data is a challenging, but essential step for establishing data reliability. There is, however, no consensus as to how best to quantify paleointensity data and which data selection processes are most effective. To address these issues, we begin to lay the foundations for a more unified and theoretically justified approach to the selection of paleointensity data. We present a new compilation of standard definitions for paleointensity statistics to help remove ambiguities in their calculation. We also compile the largest-to-date data set of raw paleointensity data from historical locations and laboratory control experiments with which to test the effectiveness of commonly used sets of selection criteria. Although most currently used criteria are capable of increasing the proportion of accurate results accepted, criteria that are better at excluding inaccurate results tend to perform poorly at including accurate results and vice versa. In the extreme case, one widely used set of criteria, which is used by default in the ThellierTool software (v4.22), excludes so many accurate results that it is often statistically indistinguishable from randomly selecting data. We demonstrate that, when modified according to recent single domain paleointensity predictions, criteria sets that are no better than a random selector can produce statistically significant increases in the acceptance of accurate results and represent effective selection criteria. The use of such theoretically derived modifications places the selection of paleointensity data on a more justifiable theoretical foundation and we encourage the use of the modified criteria over their original forms. © 2014. American Geophysical Union. All Rights Reserved.

Paleointensity data from the Atlantic Ocean are rare. We present new paleointensity data from São Miguel (Azores Islands, Portugal) based on 20 paleomagnetic sites from 13 lava flows emplaced over the last 3000years. Ten lava flows are radiocarbon dated, whereas three flows were paleomagnetically dated and one site was dated using stratigraphic relations. All the samples, previously investigated to recover paleodirections, were subjected to IZZI experiments. Importantly, the new data are internally consistent, agree with Moroccan and European datasets, and offer new constraints for global geomagnetic field models. Some of the ages of the paleomagnetically dated lava flows have been revised based on the intensity data presented here. The inferred Virtual Axial Dipole Moments (VADMs) range from 68.2 to 163.5ZAm2. A peak in field strength with an estimated age of around 600 BC is well supported by two sites from the same flow (Furna), and is comparable to the high intensity values found in Portugal for the same age and the earlier field peak at about 1000 BC in the Levant. A gradient in VADM values with latitude from northwestern Africa and across Europe between 100 and 1000 AD is confirmed as well as its absence from between 0 to 100 AD. © 2014 Elsevier B.V.

Tholeiite of the oldest oceanic crust was drilled during ODP Legs 129 and 185 at Hole 801C in the western Pacific. Fresh appearing submarine basaltic glass (SBG) was recovered from the tholeiites (∼167 Ma) which has been shown to be nearly ideal for determining absolute paleointensity. Paleointensities of the younger, off-axis, alkalic basalts (∼160 Ma), overlying the tholeiites, had been studied earlier. Here we report results from the older tholeiitic (on-axis) sequence. We subjected a total of 73 specimens from 17 cooling units to absolute paleointensity experiments. Of these, 30 specimens and six cooling unit averages met our strictest reliability criteria, yielding an average of 11.9 ± 3.9 μT. The bulk of evidence suggests a paleolatitude of the site of 14°S (with an uncertainty of 10°). This translates the intensity to a value for the virtual axial dipole moment of 28 ZAm2, slightly lower than values determined from the plagioclase crystals in the three cooling units of the younger alkalic basalts overlying the tholeiites (Tarduno & Cottrell, 2005). Our value is low when compared to the long-term median value of the field of 42 ZAm2. Our results and those of the published literature therefore support the contention of a low magnetic field strength in the Jurassic (average of 28 ± 14 ZAm2; N = 138 individual estimates), as initially suggested by Prévot et al. (1990). Our interpretation of the body of available data argue for low field strengths for the entire Jurassic extending into the early Cretaceous. Key Points We review the literature of Jurassic paleoinensity We present new paleointensity data for the Jurassic We support the Mesozoic Dipole Low hypothesis ©2013. American Geophysical Union. All Rights Reserved.

Paleomagnetic results from globally distributed lava flows have been collected and analyzed under the time-averaged field initiative (TAFI), a multi-institutional collaboration started in 1996 and designed to improve the geographic and temporal coverage of the 0-5 Ma paleomagnetic database for studying both the time-averaged field and its very long-term secular variations. Paleomagnetic samples were collected from 35 volcanic units, either lava flows or ignimbrites, in Costa Rica in December 1998 and February 2000 from the Cordilleras Central and Guanacaste, the underlying Canas, Liberia and Bagaces formations and from Volcano Arenal. Age estimates range from approximately 40 ka to slightly over 6 Ma. Although initial results from these sites were used in a global synthesis of TAFI data by Johnson et al. (2008), a full description of methodology was not presented. This paper documents the definitive collection of results comprising 28 paleomagnetic directions (24 normal, 4 reversed), with enhanced precision and new geological interpretations, adding two paleointensity estimates and 19 correlated 40Ar/39Ar radiometric ages. The average field direction is consistent with that of a geocentric axial dipole and dispersion of virtual geomagnetic poles (17.3 ± 4.6°) is in general agreement with predictions from several statistical paleosecular variation models. Paleointensity estimates from two sites give an average field strength of 26.3 μT and a virtual axial dipole moment of 65 ZAm2. The definitive results provide a useful augmentation of the global database for the longer term goal of developing new statistical descriptions of paleomagnetic field behavior. © 2013. American Geophysical Union. All Rights Reserved.

Thellier-type experiments are a method used to estimate the intensity of the ancient geomagnetic field from samples carrying thermoremanent magnetization. The analysis of Thellier-type experimental data is conventionally done by manually interpreting data from each specimen individually. The main limitations of this approach are: (1) manual interpretation is highly subjective and can be biased by misleading concepts, (2) the procedure is time consuming, and (3) unless the measurement data are published, the final results cannot be reproduced by readers. These issues compound when trying to combine together paleointensity data from a collection of studies. Here, we address these problems by introducing the Thellier GUI: a comprehensive tool for interpreting Thellier-type experimental data. The tool presents a graphical user interface, which allows manual interpretation of the data, but also includes two new interpretation tools: (1) Thellier Auto Interpreter: an automatic interpretation procedure based on a given set of experimental requirements, and 2) Consistency Test: a self-test for the consistency of the results assuming groups of samples that should have the same paleointensity values. We apply the new tools to data from two case studies. These demonstrate that interpretation of non-ideal Arai plots is nonunique and different selection criteria can lead to significantly different conclusions. Hence, we recommend adopting the automatic interpretation approach, as it allows a more objective interpretation, which can be easily repeated or revised by others. When the analysis is combined with a Consistency Test, the credibility of the interpretations is enhanced. We also make the case that published paleointensity studies should include the measurement data (as supplementary files or as a contributions to the MagIC database) so that results based on a particular data set can be reproduced and assessed by others. ©2013. American Geophysical Union. All Rights Reserved.

Warm intervals within the Pliocene epoch (5.33-2.58 million years ago) were characterized by global temperatures comparable to those predicted for the end of this century and atmospheric CO 2 concentrations similar to today. Estimates for global sea level highstands during these times imply possible retreat of the East Antarctic ice sheet, but ice-proximal evidence from the Antarctic margin is scarce. Here we present new data from Pliocene marine sediments recovered offshore of Adélie Land, East Antarctica, that reveal dynamic behaviour of the East Antarctic ice sheet in the vicinity of the low-lying Wilkes Subglacial Basin during times of past climatic warmth. Sedimentary sequences deposited between 5.3 and 3.3 million years ago indicate increases in Southern Ocean surface water productivity, associated with elevated circum-Antarctic temperatures. The geochemical provenance of detrital material deposited during these warm intervals suggests active erosion of continental bedrock from within the Wilkes Subglacial Basin, an area today buried beneath the East Antarctic ice sheet. We interpret this erosion to be associated with retreat of the ice sheet margin several hundreds of kilometres inland and conclude that the East Antarctic ice sheet was sensitive to climatic warmth during the Pliocene. © 2013 Macmillan Publishers Limited. All rights reserved.

The Earth's inner core grows by the freezing of liquid iron at its surface. The point in history at which this process initiated marks a step-change in the thermal evolution [...]