Search

The Hebrew Bible and other ancient Near Eastern texts describe Egyptian, Aramean, Assyrian, and Babylonian military campaigns to the Southern Levant during the 10th to sixth centuries BCE. Indeed, many destruction layers dated to this period have been unearthed in archaeological excavations. Several of these layers are securely linked to specific campaigns and are widely accepted as chronological anchors. However, the dating of many other destruction layers is often debated, challenging the ability to accurately reconstruct the different military campaigns and raising questions regarding the historicity of the biblical narrative. Here, we present a synchronization of the historically dated chronological anchors and other destruction layers and artifacts using the direction and/or intensity of the ancient geomagnetic field recorded in mud bricks from 20 burnt destruction layers and in two ceramic assemblages. During the period in question, the geomagnetic field in this region was extremely anomalous with rapid changes and high-intensity values, including spikes of more than twice the intensity of today’s field. The data are useful in the effort to pinpoint these short-term variations on the timescale, and they resolve chronological debates regarding the campaigns against the kingdoms of Israel and Judah, the relationship between the two kingdoms, and their administrations.

The domal structure in the core of the Troodos ophiolite exposes lower crustal (gabbro suite) and mantle rocks (ultramafic province). This structure is part of a fossil ridge-transform intersection (RTI), where an extinct spreading axis meets the fossil oceanic transform, namely the Arakapas transform. A major feature in the RTI system is the Troodos Forest-Amiandos Fault (TAF), an off-axis and axis parallel fault that was active during the Cretaceous seafloor spreading. Here we investigate the deformation across the TAF by measuring paleomagnetic vectors from 34 sites in the gabbro suite around the domal ultramafic core. Special emphasize was along an E-W transect that crosses the TAF south of the sheeted dike complex and north of the ultramafic province. We also compiled dike dips along an E-W strip (6 km wide) north of the gabbro suite. All results were compared to previous paleomagnetic studies from the sheeted dikes and the gabbro suite. Accordingly, we have found that rotations in the gabbro are very similar to those in the sheeted dikes, suggesting coupling of the upper and the lower oceanic crust during axial deformation of seafloor spreading. All rotation axes were horizontal and parallel to the dike strikes, i.e., parallel to the extinct spreading axis. Rotations increase gradually towards the TAF from both sides, eastward in the footwall and westward in the hanging wall. The most plausible scenario is an upward and downward deflection in the footwall and the hanging wall, respectively, similarly described theoretically for the early stages of detachment development. The orientations of the rotation axes of all paleomagnetic vectors indicate spreading-related deformation. This suggests that the relative uplift of the deep-seated rocks was by the development of a young detachment during seafloor spreading rather than serpentinite diapirism. The detachment occurrence in the outside-corner is explained here by the shift from orthogonal to curved axis, inferred from sheeted dike orientations.

Magnetic properties of marine sediments are dictated not only by the detrital mineralogy, but also by diagenetic processes that can start instantaneously after deposition and may proceed for a long time as sediments are buried. Early diagenesis encompass a range of biochemical reactions associated with bacterial respirations. These may include oxidation at the water sediment interface, iron reduction, sulfate reduction, and anaerobic oxidation of methane (AOM) if methane is present. To investigate the link between diagenesis, sedimentary magnetic properties and quality of paleomagnetic recording, we collected eight 6m-long piston-cores from the Holocene Eastern Mediterranean continental shelf in four locations. Two locations are characterized by high concentration of methane and detectable sulfate-methane transition zone (SMTZ) at depth of 1-4 m. In the other two locations, organoclastic sulfate reduction is dominant throughout the entire cores. Sedimentation rates in this region range between 1 - 5 mm/year. In all cores, concentrations of sulfate, methane and ferrous iron were measured from the pore water. The geochemical data were compared to the mineral magnetic profile that include a range of parameters calculated from IRM, ARM, low- and high- field susceptibility, hysteresis, and FORCs. Paleomagnetic time-series of declination and inclination were obtained from demagnetization experiments carried out in 2 cm resolution. Age models were constructed from radiocarbon dating of carefully collected foraminifera. The results show a consistent link between the diagenetic zones and the magnetic mineralogy: Increase of magnetic properties in the shallow ferruginous zone, decay of magnetic parameters in the sulfate reduction zones, rapid decrease at the SMTZ and stabilization at the methanogenic zone. XHR-SEM analysis show multiple generations of greigite and pyrite framboids at all depths and unaltered detrital titanomagnetites. We find that except a short time interval below the SMTZ of one core, the paleomagnetic directions in these sediments do not represent the expected directions of the geomagnetic field. We conclude that continuous organoclastic sulfate reduction in marine sediments might have a profound effect on the quality of paleomagnetic recording, but AOM at the SMTZ may help stabilize the magnetic phase.

The vast number of well-dated archaeological sites in the Levant and Mesopotamia, provide a unique opportunity for detailed archaeomagnetic research. In an effort to exploit the full archaeomagnetic potential of this region, we assemble an archaeointensity compilation based on two separate datasets which were constructed independently in different regions using different methods. The Mesopotamian dataset, comprised of 123 groups of samples (fragments) collected from Syria and north-western Iraq was analyzed for the most part using the Triaxe method. The Levantine dataset comprised of more than 150 groups of samples collected from Israel was analyzed using Thellier-IZZI-MagIC method. Together, these partially overlapped datasets compose a continuous record that spans over more than six millennia - from ~6300 BCE to 300 CE. Some intervals in this time span, such as 2000 BCE - 500 BCE and 300 BCE- 300 CE, include at least one context per century with a secure absolute age, which is based on radiocarbon or clear historically-dated events, providing sub-centennial temporal resolution. With these data we construct the second generation of the Levantine Archaeomagnetic Curve (LAC.v.2.0) that enables robust high-resolution analysis of the rates and amplitudes of geomagnetic secular variations. The LAC exhibits four geomagnetic spikes associated with the Levantine Iron Age Anomaly (LIAA), between 1050 and 600 BCE, with virtual axial dipole moment (VADM) reaching values of 155-162 ZAm2 - higher than predicted by currently reconstructed geomagnetic field models. Rates of VADM change associated with the four spikes are ~0.7-1.1 ZAm2/year – at least twice the maximum rate inferred from direct observations, which began ~190 years ago. The increase leading to the first spike from ~1750 BCE to ~1050 BCE depicts the Holocene largest change in field intensity. The highest value in the LAC is more than 3 times greater than the lowest value providing conservative bounds for the local variability of the geomagnetic field. Finally, we present a new compilation of archaeomagnetic directions obtained from 70 structures in Israel which, together with the intensity curve, provide a further step towards a full-vector description of the field.

This study demonstrates the feasibility of speleothem magnetism as a paleo-hydrology proxy in speleothems growing in semi-arid conditions. Soil-derived magnetic particles in speleothems retain valuable information on the physicochemical conditions of the overlying soil, and changes in bedrock hydrology. Yet, the link between magnetic and isotopic proxies of speleothems has been only partly established. We reveal strong coupling between the inflow of magnetic particles (quantified using the magnetic flux index, IRMflux) and δ13C in two Holocene speleothems from Soreq Cave (Israel). The stalagmite record spans from ca. 9.7 to ca. 5.4 ka, capturing the warm-humid conditions associated with the early Holocene and the transition to mid-Holocene wet-dry cycles. Extremely low IRMflux during the early Holocene, indicating minimal contribution from the overlying soil, is accompanied by anomalously high δ13C (approaching bedrock values) hypothesized to be caused by high rainfall and soil erosion. By contrast, IRMflux during the mid-Holocene covaries with the saw-tooth cyclicity of δ13C and δ18O, interpreted as rapid fluctuations in rainfall amount. The peaks in IRMflux precede the negative (wet) δ13C peaks by ~60–120 yr. The apparent lag is explained as a rapid physical translocation of overlying soil particles via groundwater (high IRMflux) as a response to increasing rainfall, compared with slower soil organic matter turnover rates (10–102 yr).

Shaar et al. (2022)

Apart from magnetostratigraphy, archaeomagnetism is rarely used in Middle and Late Pleistocene sites. Here we present detailed palaeomagnetic analyses of cemented hearths and burnt chert items from Amud Cave, Israel (68–55 ka) - two types of materials common in Levantine Middle Palaeolithic cave sites. Both materials are shown to be recorders of the geomagnetic field and were used to reconstruct either the ancient field direction (for cemented hearths) or intensity (palaeointensity) (for chert) at the time of the last burning or shortly after. We test the utility of palaeomagnetic data to further our understanding of temporal aspects of occupations in the cave by comparing the dispersion of the palaeomagnetic data to the known characteristics of geomagnetic secular variation in the Holocene. We show that divergent palaeointensities can help identify diachronic burning events, suggesting different activity patterns in two areas of the cave. Additionally, we used palaeomagnetic directional vectors to distinguish between a well-preserved hearth and one that had been mixed prior to cementation. Using rock magnetic investigations, we demonstrate that magnetic methods can be used as a relatively fast and inexpensive method to identify burning of cherts in antiquity above 500 °C. The palaeomagnetic results are in agreement with results of previous studies at Amud Cave, obtained by other independent methods. This study shows that palaeomagnetic methods can serve as a powerful tool in the study of Palaeolithic sites.

Diagenetic processes in anoxic sedimentary environments influence sediment magnetic properties mainly through dissolution of detrital magnetite and precipitation of authigenic greigite. Recently exposed late Holocene Dead Sea sediments provide an opportunity to study the processes governing greigite formation and preservation, and their relation to different hydrological settings. Magnetic data and pore-fluid compositions were obtained from three Holocene sections along a N-S transect on the western Dead Sea shore: Og, Ein-Feshkha (EF), and Ein-Gedi. The northern sections are closer to the major freshwater source to the Dead Sea-the Jordan River. Detrital titanomagnetite is present at all sections, but greigite is the dominant magnetic phase at Og and EF. Bulk rock magnetic data vary between and within the sections by over 3 orders of magnitude, where higher values indicate higher greigite concentrations. At the three sites, pore fluids have similar or lower salinity than the modern and Holocene Dead Sea brine, with variable and dissolved iron (Fe2+) and sulfate (SO42−). Magnetic property changes are reflected by iron and/or sulfate microbial reduction that controlled sedimentary greigite formation. We propose that the N-S greigite decrease suggests that anoxic microbial activity was controlled by labile organic matter and/or reactive iron brought by, or formed as a result of, freshwater influx from the Jordan River. Hence, greigite concentration changes depended on past freshwater input to the hypersaline lake and proximity to the freshwater source. The apparent relationship between hydrological conditions and magnetic properties provides a new method to trace past hydrological changes in the Dead Sea.

Paleomagnetic analysis of archaeological materials is crucial for understanding the behavior of the geomagnetic field in the past. As it is often difficult to accurately date the acquisition of magnetic information recorded in archaeological materials, large age uncertainties and discrepancies are common in archaeomagnetic datasets, limiting the ability to use these data for geomagnetic modeling and archaeomagnetic dating. Here we present an accurately dated reconstruction of the intensity and direction of the field in Jerusalem in August, 586 BCE, the date of the city's destruction by fire by the Babylonian army, which marks the end of the Iron Age in the Levant. We analyzed 54 floor segments, of unprecedented construction quality, unearthed within a large monumental structure that had served as an elite or public building and collapsed during the conflagration. From the reconstructed paleomagnetic directions, we conclude that the tilted floor segments had originally been part of the floor of the second story of the building and cooled after they had collapsed. This firmly connects the time of the magnetic acquisition to the date of the destruction. The relatively high field intensity, corresponding to virtual axial dipole moment (VADM) of 148.9 ± 3.9 ZAm2, accompanied by a geocentric axial dipole (GAD) inclination and a positive declination of 8.3°, suggests instability of the field during the 6th century BCE and redefines the duration of the Levantine Iron Age Anomaly. The narrow dating of the geomagnetic reconstruction enabled us to constrain the age of other Iron Age finds and resolve a long archaeological and historical discussion regarding the role and dating of royal Judean stamped jar handles. This demonstrates how archaeomagnetic data derived from historically-dated destructions can serve as an anchor for archaeomagnetic dating and its particular potency for periods in which radiocarbon is not adequate for high resolution dating.

Archaeomagnetic dating is a firmly established dating technique applicable to a wide variety of heat-treated anthropological materials and is advantageous for sites that lack materials suitable for radiocarbon dating. To correct recent misinterpretations of the method, we provide examples of how archaeomagnetic dating curves are calibrated and show how, in some instances, the technique can provide superior results. We emphasize that no single dating technique is capable of resolving the challenging chronology controversies in the Levant, and instead argue that multiple dating methods must be integrated in order to achieve the highest possible temporal resolution.

Author(s): Behar, Nicole; Shaar, Ron; Tauxe, Lisa; Asefaw, Hanna; Ebert, Yael; Heimann, Ariel; Koppers, Anthony AP; Ron, Hagai

Cave sediments pose dating challenges due to complex depositional and post-depositional processes that operate during their transport and accumulation. Here, we confront these challenges and investigate the stratified sedimentary sequence from Wonderwerk Cave, which is a key site for the Earlier Stone Age (ESA) in Southern Africa. The precise ages of the Wonderwerk sediments are crucial for our understanding of the timing of critical events in hominin biological and cultural evolution in the region, and its correlation with the global paleontological and archaeological records. We report new constraints for the Wonderwerk ESA chronology based on magnetostratigraphy, with 178 samples passing our rigorous selection criteria, and fourteen cosmogenic burial ages. We identify a previously unrecognized reversal within the Acheulean sequence attributed to the base of the Jaramillo (1.07 Ma) or Cobb Mtn. subchrons (1.22 Ma). This reversal sets an early age constraint for the onset of the Acheulean, and supports the assignment of the basal stratum to the Olduvai subchron (1.77–1.93 Ma). This temporal framework offers strong evidence for the early establishment of the Oldowan and associated hominins in Southern Africa. Notably, we found that cosmogenic burial ages of sediments older than 1 Ma are underestimated due to changes in the inherited 26Al/10Be ratio of the quartz particles entering the cave. Back calculation of the inherited 26Al/10Be ratios using magnetostratigraphic constraints reveals a decrease in the 26Al/10Be ratio of the Kalahari sands with time. These results imply rapid aeolian transport in the Kalahari during the early Pleistocene which slowed during the Middle Pleistocene and enabled prolonged and deeper burial of sand while transported across the Kalahari Basin.

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.

Geomagnetic field measurements from the past few centuries show heightened secular variation activity in the southern hemisphere associated with the south Atlantic anomaly (SAA). It is uncertain whether geomagnetic anomalies at a similar scale have existed in the past owing to limited coverage and uncertainties in the paleomagnetic database. Here we provide new evidence from archaeological sources in the Levant suggesting a large positive northern hemisphere anomaly, similar in magnitude to the SAA during the 9th–8th centuries BCE, called “Levantine Iron Age anomaly”. We also report an additional geomagnetic spike in the 8th century. The new dataset comprises 73 high precision paleointensity estimates from ca. 3000 BCE to 732 BCE, and five directional measurements between the 14th and the 9th centuries BCE. Well-dated pottery and cooking ovens were collected from twenty archaeological strata in two large contemporaneous stratigraphical mounds (tells) in Israel: Tel Megiddo and Tel Hazor. The new data are combined with previously published data and interpreted automatically using the PmagPy Thellier GUI program. The Tel Megiddo and Tel Hazor data sets demonstrate excellent internal consistency and remarkable agreement with published data from Mesopotamia (Syria). The data illustrate the evolution of an extreme geomagnetic high that culminated in at least two spikes between the 11th and the 8th centuries BCE (Iron Age in the Levant). The paleomagnetic directional data of the 9th century BCE show positive inclination anomalies, and deviations of up to 22° from the averaged geocentric axial dipole (GAD) direction. From comparison of the Levantine archaeomagnetic data with IGRF model for 2015 we infer the “Levantine Iron Age anomaly” between the 10th and the 8th centuries BCE is a local positive anomaly. The eastward extent of the anomaly is currently unknown.

Accumulations of ice and dust mixtures may acquire magnetization during deposition in a manner analogous to sedimentary rocks. Here, we consider the process of particles descending through an atmosphere and depositing in a preferential orientation that serves to record the ambient magnetic field during emplacement. We use a simple model for the settling and reorientation of ice particles with magnetic inclusions that includes magnetic torque, aerodynamic forces and gravity, to investigate the parameter space governing the process. For fields in the range of 10's – 100's μ T we find that ice particles of sizes up to ∼100 μ m which contain smaller magnetic grains as nuclei will produce a deposit indeed magnetized in the direction aligned with the applied field, but with a moment that is independent of the field strength. For particles in the 100's μ m range, the magnetic moment increases with the field strength. To demonstrate the effect experimentally, we performed a suite of laboratory deposition simulations followed by measurements of the magnetic moment of the samples. We show that in the idealized laboratory conditions dusty ice magnetizes in the direction of the applied field, with the alignment increasing with its intensity. For the chosen conditions, the magnetization increases rapidly with field intensity in the range 10 – 200 μ T , and approaches a maximal value above that. For a mixture with dust/ice ratio of 5 × 10 − 3 we obtained maximal magnetization values in the range 1.6 × 10 − 5 – 3 × 10 − 3 A m2/kg, depending on the distribution of particle sizes. We show that magnetic particle concentration in the ice determines the level of magnetic remanence, and conclude that the remanent magnetization of natural ice deposit in various settings may be measurable (if unobscured by post-depositional, wind, or other effects) and thus could provide a new paleomagnetic record on Earth and other planetary objects.

Marine and lacustrine sediments represent an important source of global paleomagnetic data. Although it is usually assumed that detrital iron oxides record most of the magnetic signal in sediments, iron sulfides—which form during bacterial sulfate reduction—can also represent a significant source of sedimentary magnetism. Knowing how sulfate reduction impacts sedimentary magnetism is critical to the interpretation of paleomagnetic records. Here, we show that three distinct types of magnetic particles can be produced by bacterial sulfate reduction, each of which impacts the bulk sediment magnetism in a distinct way. We combined magnetic force microscopy and electron probe microanalysis to image magnetic mineral extracts from Dead Sea sediments from a glacial period and an interglacial period. In sediments from the dry interglacial period, during which bacterial sulfate reduction was suppressed, we found greigite framboids (Fe3S4) with strong intergrain magnetic interactions. Contrastingly, in sediments from the wet glacial period, which experienced extensive sulfate reduction, pyrite (FeS2) is the dominant sulfide phase. High-resolution magnetic imaging of glacial pyrite reveals that greigite is present as single-domain particles within the pyrite. We also found that as titanomagnetite grains undergo bacterially mediated alteration to form pyrite, the original magnetic grains become divided into smaller regions, which potentially facilitates acquisition of secondary magnetization by the reorganization of these magnetic domains. Our results provide a previously undocumented mechanism by which bacterially mediated alteration can overwrite primary detrital magnetic records.

Historical events are sometimes expressed in destruction layers. We present here a study in which aspects of construction, destruction, and chronostratigraphy of fired mud bricks were explored using archaeomagnetism, infrared spectroscopy, and micromorphology. We measured 88 oriented samples mostly collected from one stratum, dated ca. 1000 B.C.E., representing a destroyed late Canaanite (late Iron Age I) city in Tel Megiddo, Israel. Firing temperatures, evaluated from infrared spectroscopy, micromorphology, and high‐temperature magnetic susceptibility cycles, range between 300°C and 800°C. Samples studied in one archaeomagnetic site yield a single stable magnetization vector in demagnetization experiments. Archaeomagnetic site means of three standing walls are grouped near the expected direction of the ancient geomagnetic field. We propose that walls in the destruction layer were constructed from sun‐dried mud bricks that later burned during the destruction. Collapsed bricks and tilted walls show variable directions, diagnostic for the relative timing of collapse and cooling of bricks, during and following the destruction event. In addition, we attempt to assign stratigraphic affiliation based on archaeomagnetic considerations to standing walls, which are spatially disconnected from the studied destruction layer. Altogether, this study demonstrates the usefulness of archaeomagnetism to understanding site formation processes related to fire and destruction.

Earth’s magnetic field, one of the most enigmatic physical phenomena of the planet, is constantly changing on various time scales, from decades to millennia and longer. The reconstruction of geomagnetic field behavior in periods predating direct observations with modern instrumentation is based on geological and archaeological materials and has the twin challenges of ( i ) the accuracy of ancient paleomagnetic estimates and ( ii ) the dating of the archaeological material. Here we address the latter by using a set of storage jar handles (fired clay) stamped by royal seals as part of the ancient administrative system in Judah (Jerusalem and its vicinity). The typology of the stamp impressions, which corresponds to changes in the political entities ruling this area, provides excellent age constraints for the firing event of these artifacts. Together with rigorous paleomagnetic experimental procedures, this study yielded an unparalleled record of the geomagnetic field intensity during the eighth to second centuries BCE. The new record constitutes a substantial advance in our knowledge of past geomagnetic field variations in the southern Levant. Although it demonstrates a relatively stable and gradually declining field during the sixth to second centuries BCE, the new record provides further support for a short interval of extreme high values during the late eighth century BCE. The rate of change during this “geomagnetic spike” [defined as virtual axial dipole moment > 160 ZAm 2 (10 21 Am 2 )] is further constrained by the new data, which indicate an extremely rapid weakening of the field (losing ∼27% of its strength over ca. 30 y).

[1] 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.

[1] Dendrites are crystals that grow in branches that diverge along crystallographically defined directions. Despite the importance of dendrites in paleomagnetic research, little is known about how dendrites act as magnetic recorders, because they exhibit complicated magnetic domain structures. In this study, we experimentally examine how textures and sizes of dendrites affect their magnetic domain structure and magnetic properties. We study two basaltic glass samples and three synthetic slag samples, which collectively define a wide range of dendritic morphologies. We use electron microscopy to characterize the morphology of the dendrites and magnetic force microscopy (MFM) to observe their magnetic domain structure. We characterize the dendrites’ bulk properties by firs-order reversal curve distributions, Thellier-style paleointensity experiments, anisotropy of remanence, and anisotropy of susceptibility. The samples with the thinnest dendrites have high coercivity, stable single-domain (SD) – pseudo single-domain magnetization, and yield ideal Arai plots. By contrast, the sample with the thickest dendrites has the lowest coercivity and shows the most extreme multidomain (MD) behavior. All samples except one, exhibit significant remanence and susceptibility anisotropy. MFM observations show that dendrites built from branches of interconnected octahedra, typical for basaltic glass, have a stable, high coercivity, SD-like magnetization despite the fact that their overall dimensions exceed the SD-MD threshold. Their stability is likely due to interactions between the octahedra and their narrow rod-like interconnections. Dendrites that crystallize in faster cooling environments, such as in archaeological slag, display finer branch thicknesses (

Recent archaeomagnetic data from ancient Israel revealed the existence of a so-called “Levantine Iron Age geomagnetic anomaly” (LIAA) which spanned the first 350 years of the first millennium before the Common Era (B.C.E.) and was characterized by a high averaged geomagnetic field (virtual axial dipole moments, VADM > 140 Z Am2, nearly twice of today's field), short decadal-scale geomagnetic spikes (VADM of 160–185 Z Am2), fast field variations, and substantial deviation from dipole field direction. The geographic constraints of the LIAA have remained elusive due to limited high-quality paleointensity data in surrounding locations. Here we report archaeointensity data from Georgia showing high field values (VADM > 150 Z Am2) in the tenth or ninth century B.C.E., low field values (VADM < 60 Z Am2) in the twelfth century B.C.E., and fast field variation in the fifth and fourth centuries B.C.E. High field values in the time frame of LIAA have been observed so far only in three localities near the Levant: Eastern Anatolia, Turkmenistan, and now Georgia, all located east of longitude 30°E. West of this, in the Balkans, field values in the same time are moderate to low. These constraints put geographic limits on the extent of the LIAA and support the hypothesis of an unusually intense regional geomagnetic anomaly during the beginning of the first half of the first millennium B.C.E., comparable in area and magnitude (but of opposite sign) to the presently active South Atlantic anomaly. ©2017. American Geophysical Union. All Rights Reserved.

Magnetization in rock samples is crucial for paleomagnetometry research, as it harbors valuable geological information on long term processes, such as tectonic movements and the formation of oceans and continents. Nevertheless, current techniques are limited in their ability to measure high spatial resolution and high-sensitivity quantitative vectorial magnetic signatures from individual minerals and micrometer scale samples. As a result, our understanding of bulk rock magnetization is limited, specifically for the case of multi-domain minerals. In this work, we use a newly developed nitrogen-vacancy magnetic microscope, capable of quantitative vectorial magnetic imaging with optical resolution. We demonstrate direct imaging of the vectorial magnetic field of a single, multi-domain dendritic magnetite, as well as the measurement and calculation of the weak magnetic moments of an individual grain on the micron scale. These results pave the way for future applications in paleomagnetometry and for the fundamental understanding of magnetization in multi-domain samples.

This paper presents the results of the 2009 excavations at Site 30 in the Timna Valley, Israel. The results, coupled with a suite of 11 new radiocarbon dates, fix the chronology of the site between...

Experimental techniques to determine paleomagnetic field intensity are based on a theoretical framework that is valid only for single-domain (SD) ferromagnetic particles. Yet, most of the available materials exhibit distinctly non-SD properties. Designing the optimal paleointensity methodology for non-SD is, therefore, a fundamental challenge in paleomagnetism. The objective of this study is to experimentally test the IZZI Thellier absolute paleointensity method on small MD recorders. The test has two purposes: 1) to describe the characteristic non-SD patterns occurring in Arai plots, and 2) to identify the optimal approach in interpreting non-SD behavior. We carried out paleointensity experiments on 40 specimens from 4 synthetic re-melted slag samples with identical magnetic properties (mineralogy, texture, and non-SD state) produced under different field intensities. We ran three batches of IZZI experiments using different conditions that allow for a detailed characterization of the non-SD behavior. We find that the curvature of the Arai plot is systematically dependent on the angle and the proportion between the field used in the paleointensity experiment (BTRM) and the field in which the NRM was acquired (BNRM). Straight-line Arai plot occur when the two fields are parallel and equal, and seems to always give the ‘true’ slope. Convex curves occur when BTRM is parallel and significantly stronger than BNRM. Concave curves occur in all the other cases and yield two end-case slopes that are always different than the ‘true’ slope. In addition, zigzagged patterns increase with the angle the proportion between BTRM and BNRM. We test the accuracy of the ‘best fitting’ line approach and conclude that ‘best fitting’ line in curved plots cannot provide robust paleointensity estimates. Yet, the two ‘end-case’ slopes in concave curves provide adequate constraints for the true value. We introduce a new procedure to calculate a 95% confidence interval of the paleointensity from curved plots using bootstrap statistics. We substantiate the new procedure by conducting two independent tests. The first uses synthetic re-melted slag produced under known field intensities — 3 SD samples and 4 non-SD samples. The second compares paleointensity determinations from archeological slag samples of the same age — 34 SD samples and 10 non-SD samples. The two tests demonstrate that the bootstrap technique may be the optimal approach for non-ideal dataset.

article The intensity of the geomagnetic field varies over different time scales. Yet, constraints on the maximum intensity of the field as well as for its maximum rate of change are inadequate due to poor temporal resolution and large uncertainties in the geomagnetic record. The purpose of this study is to place firm limits on these fundamental properties by constructing a high-resolution archaeointensity record of the Levant from the 11th century to the early 9th century BCE, a period over which the geomagnetic field reached its maximum intensity in Eurasia over the past 50,000 years. We investigate a 14 C-dated sequence of ten layers of slag material, which accumulated within an ancient industrial waste mound of an Iron Age copper-smelting site in southern Israel. Depositional stratigraphy constrains relative ages of samples analyzed for paleointensity, and 14 C dates from different horizons of the mound constrain the age of the whole sequence. The analysis yielded 35 paleointenisty data points with accuracy better than 94% and precision better than 6%, covering a period of less than 350 years, most probably 200 years. We construct a new high-resolution quasi-continuous archaeointensity curve of the Levant that displays two dramatic spikes in geomagnetic intensity, each corresponding to virtual axial dipole moment (VADM) in excess of 200 ZAm 2

In order to determine the nature of remanence carriers of the layered gabbro of the Troodos ophiolite, Cyprus, we report optical and electron microscopic observations, together with rock magnetic and paleomagnetic experiments. Above all, the study aims to understand and clarify the time of magnetic acquisition relative to the brittle deformation of the oceanic crust manifested by the ridge-transform intersection (Solea graben and the Arakapas transform). Petrographic examination of pyroxene grains revealed isolated magnetite inclusions ranging in size from single-domain (SD) to multi-domain (MD) and in addition, MD pyrrhotite inclusions residing in veins and cracks. Thermal demagnetization and thermomagnetic procedures indicate two components, low and high temperature. We argue that the low temperature component, The SD magnetite inclusions exsolved in pyroxene are the dominant carriers of magnetic remanence in Troodos Gabbro. The initial formation of these inclusions occurred via exsolution reaction at temperatures between 520 and 850 °C, above the Curie temperature of pure magnetite during the solidification of the magma. Therefore, acquisition of remanent magnetization of the Troodos gabbro took place during the earliest stages of crustal accretion, before any brittle deformation associated with the spreading ridge and the transform fault occurred.

The Middle-Eastern copper slag is a promising new material for studying intensity variations in the geomagnetic field with high resolution and precision. The purpose of this study is to test the accuracy of archaeointensity estimates determined using copper slag by addressing two questions: 1) “Does slag material display the magnetic properties required for valid Thellier experiments?” and 2) “What is the accuracy of the archaeointensity estimates derived from Thellier-style experiments on optimal samples?” We address the first question through a comprehensive microscopic and magnetic study of representative archaeological slag samples in order to identify the properties responsible for optimal behavior in Thellier experiments. To address the second question, we reproduced slag samples in the laboratory under controlled magnetic fields and analyzed them using the same IZZI paleointensity technique used for the ancient slag. Microscopic analyses of the archaeological slag show that ferromagnetic phases occur as three-dimensional dendritic structures whose branches consist of submicron-elongated particles. Magnetic analyses show that these dendrites behave as an assemblage of shape-controlled, single-domain-like particles and that their magnetization is thermoremanent. We conclude that slag material can be magnetically suitable for valid Thellier experiments. The laboratory-produced slag material demonstrated similar magnetic and mineralogical properties as the archaeological slag. IZZI experiments showed that non-linear TRM acquisition, even at field strengths similar to Earth's, and TRM anisotropy are important factors to monitor during paleointensity studies of slag material. Anisotropy and non-linearity are probably related to the dendritic shape of the oxide grains. Intensity estimates derived from three laboratory-produced slag samples demonstrated accuracy to within ∼ 5% after applying the required corrections.

article i nfo In paleomagnetism, periods of high field intensity have been largely ignored in favor of the more spectacular directional changes associated with low field intensity periods of excursions and reversals. Hence, questions such as how strong the field can get and how fast changes occur are still open. In this paper we report on data obtained from an archaeometallurgical excavation in the Middle East, designed specifically for archaeomagnetic sampling. We measured 342 specimens from 72 samples collected from a 6.1 m mound of well stratified copper production debris at the early Iron Age (12th-9th centuries BCE) site of Khirbat en- Nahas in Southern Jordan. Seventeen samples spanning 200 yr yielded excellent archaeointensity results that demonstrate rapid changes in field intensity in a period of overall high field values. The results display a remarkable spike in field strength, with sample mean values of over 120 μT (compared to the current field strength of 44 μT). A suite of 13 radiocarbon dates intimately associated with our samples, tight control of sample location and relative stratigraphy provide tight constraints on the rate and magnitude of changes in archaeomagnetic field intensities.

Data on the past intensity of Earth's magnetic field (paleointensity) are essential for understanding Earth's deep interior, climatic modeling, and geochronology applications, among other items. Here we demonstrate the possibility that much of available paleointensity data could be biased by instability of thermoremanent magnetization (TRM) associated with non-single-domain (SD) particles. Paleointensity data are derived from experiments in which an ancient TRM, acquired in an unknown field, is replaced by a laboratory-controlled TRM. This procedure is built on the assumption that the process of ancient TRM acquisition is entirely reproducible in the laboratory. Here we show experimental results violating this assumption in a manner not expected from standard theory. We show that the demagnetization-remagnetization relationship of non-SD specimens that were kept in a controlled field for only 2 y show a small but systematic bias relative to sister specimens that were given a fresh TRM. This effect, likely caused by irreversible changes in micromagnetic structures, leads to a bias in paleointensity estimates.

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.

We investigate a kilometer‐scale steep relay rampstructure in the Galilee, northern Israel. The data indi-cate an asymmetric fold above the buried tip of thesouthern boundary fault as well as at the lower partof the ramp. Structural analysis suggests that the buriedtip of the bounding fault is only a few hundred metersbelow the surface. A sequence of colluvial wedges,exposed at the base of the southern flank of the relayramp over the blind tip of the boundary fault, presentssuccessively decreasing angles from 67° at the bottomto 30° at the top. The lower three wedges rest at anglesgreater than the angle of repose, suggesting tilting dur-ing and after deposition. We suggest that the increasingangle of the wedges is the result of the development ofthe monocline over the buried tip ofthe boundaryfault,which in turn reflects the development of the entirerelay ramp structure. Paleomagnetic measurements testthis hypothesis. Samples from three of the four lowerwedges and the carbonate fill in the fractured bedrockyield northerly declination and positive inclination.Inclination anomalies are 19.4° ± 3.3° and 11.3° ±1.8° for the fractured bedrock and the wedges, respec-tively. These anomalies reflect part of the tilt. Paleo-magnetic measurements and optically stimulatedluminescence (OSL) constrain the time of initial defor-mation of the monocline above the southern boundaryfault between 780 ka and 461 ± 75 ka (OSL age ofwedge 1) and the termination of deformation to 176 ±22 ka (OSL age of wedge 5). These bounding agessuggest tilting rates that range between 1°/16 kyrand 1°/8 kyr. Termination of tilting along the south-ern flank of the relay ramp and the fact that the tipof the buried normal fault is only a few hundreds ofmeters below the surface suggest that the relay rampmay be close to being breached by a connecting fault.

Archeomagnetic records are an important source of information on the past behavior of the geomagnetic field. Frequently, however, coeval archeomagnetic intensity (archeointensity) datasets from nearby locations display significant discrepancies, hampering precise reconstruction of high-resolution secular variation curve. This is the case for the time interval between the later phase of the Early Bronze and the early phase of the Late Bronze Ages (23rd–15th centuries BCE) in the Levant and Mesopotamia. We address the problem by cross-correlating archeointensity datasets from four major multilayered archeological sites in the southern Levant (Hazor and Megiddo), northern Levant (Ebla), and western Upper Mesopotamia (Mari). We report new archeointensity data, obtained using the Thellier-IZZI-MagIC and the Triaxe methods, from six strata at Hazor and four radiocarbon-dated strata at Megiddo. From 39 pottery fragments, 199 specimens passed our selection criteria, from which we calculated the mean archeointensity for each stratum. To strengthen the comparison of these data with previously published data from Mari and Ebla, obtained using the Triaxe method, we conducted a blind test of the methods that resulted in indistinguishable results or a difference of less than 1 μT. The synchronized compilation, constrained by radiocarbon data from Megiddo, displays a V-shaped pattern with a prominent minimum of at least 200 years centered around the 18th century BCE. The study highlights the importance of stacking archeomagnetic data obtained by different archeointensity methods only after cross-testing the methods and ensuring that archeological samples were dated in a consistent manner.

Archeomagnetic data from the Levant revealed periods within the Holocene with fast and extreme changes in the geomagnetic field. Yet, the availability of the archeomagnetic data is sporadic and the correlation with the available sedimentary records from the region is rather poor. To further explore decadal variations in the directions of the field, we investigate three outcrops of the late Holocene Dead Sea that are exposed along the western retreating shores of the modern lake. The sediments were deposited under spatially varying limnological-environmental conditions, influencing their magnetic properties. The southern section, located near Ein-Gedi Spa (EG section) is dominated by detrital titanomagnetite whereas the northern sections, Nahal Og (Og section) and Ein-Feshkha (EF section), are dominated by authigenic greigite. The chronology of the sections was established by radiocarbon dating of short lived organic debris. The magnetic data were obtained in a 2 cm resolution. The EF section, spanning the time interval from ca. 2,500 cal yr BP to ca. 1,000 cal yr BP, is dominated by greigite and thus providing the most robust geomagnetic record with precise paleomagnetic directions. Greigite forms very early in the sediment and the effects of smoothing and the inclination shallowing are negligible. The new data reveal a maximal deviation of 20° from the geocentric axial dipole field between 2,400 to 2,200 cal yr BP accompanied with a fast swing in inclination from 60° to 35° over about a century. This suggests high geomagnetic field activity associated with the Levantine Iron Age geomagnetic anomaly.