Shaar R, Ben-Yosef E. paleomagnetic geochronology of Quaternary sequences in the Levant. In: Quaternary Environments, Climate Change, and Humans in the Levant. ; 2017. paleomagnetic_quaternary_levant_4_0_withfigures.pdf
Tauxe L, Shaar R, Jonestrask L, Swanson‐Hysell NL, Minnett R, Koppers AAP, Constable CG, Jarboe N, Gaastra K. PmagPy: Software package for paleomagnetic data analysis and a bridge to the Magnetics Information Consortium (MagIC) Database." Geochemistry, Geophysics, Geosystems (2016). eochemistry, Geophysics, Geosystems [Internet]. 2016. Publisher's Version
Shaar R, Tauxe L, Ron H, Ebert Y, Zuckerman S, Finkelstein I, Agnon A. Large geomagnetic field anomalies revealed in Bronze to Iron Age archeomagnetic data from Tel Megiddo and Tel Hazor, Israel. Earth and Planetray Science Letters [Internet]. 2016;442 :173-85. Publisher's Version
Shaar R, Tauxe L. Instability of thermoremanence and the problem of estimating the ancient geomagnetic field strength from non-single-domain recorders. Proceedings of the National Academy of Sciences of the United States of America [Internet]. 2015. Publisher's VersionAbstract

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.

Shaar R, Tauxe L, Ben-Yosef E, Kassianidou V, Lorentzen B, Feinberg JM, Levy TE. Decadal-scale variations in geomagnetic field intensity from ancient Cypriot slag mounds. Geochemistry Geophysics Geosystems [Internet]. 2015;16 :195-214. Publisher's VersionAbstract
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.
Ben-Yosef E, Shaar R, Tauxe L, Levy TE, Kassianidou V, Lorentzen B. Archaeomagnetic investigations of copper slag deposits in Cyprus. In: Kassianidou V, Papasavvas G Studies in Mediterranean Archaeology, Introduction to the Archaeometallurgy of Cyprus. Uppsala: Astrom Editions ; 2014.
Paterson GA, Tauxe L, Biggin AJ, Shaar R, Jonestrask LC. On improving the selection of Thellier-type paleointensity data. Geochemistry Geophysics Geosystems [Internet]. 2014;15 :1180-1192. Publisher's VersionAbstract
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.
Shaar R, Feinberg JM. Rock magnetic properties of dendrites: insights from MFM imaging and implications for paleomagnetic studies. Geochemistry Geophysics Geosystems [Internet]. 2013;14 :407-421. Publisher's VersionAbstract
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 (< 1 mu m) and few, if any, octahedra. The tips of these dendrites consist of closely intergrown, rounded, acicular branches that behave as pseudo single-domain grains due to interactions between the branches. The largest, thickest dendrites show MD behavior and MFM images show their magnetic domain walls to be arranged in crystallographically controlled, geometrically repetitive patterns within elongated branches, which give rise to their anisotropy characteristics.
Shaar R, Tauxe L. Thellier GUI: An integrated tool for analyzing paleointensity data from Thellier-type experiments. Geochemistry Geophysics Geosystems [Internet]. 2013;14 :677-692. Publisher's VersionAbstract
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.
Ben-Yosef E, Shaar R, Tauxe L, Ron H. A New Chronological Framework for Iron Age Copper Production at Timna (Israel). Bulletin of the American Schools of Oriental Research. 2012;366 :1-41.
Zilberman E, Ron H, Shaar R. Evaluating the potential seismic hazards of the Ahihud Ridge fault system by paleomagnetic and morphological analyses of calcretes. 2011.
Shaar R, Ben-Yosef E, Ron H, Tauxe L, Agnon A, Kessel R. Geomagnetic field intensity: How high can it get? How fast can it change? Constraints from Iron Age copper slag. Earth and Planetary Science Letters [Internet]. 2011;301 :297-306. Publisher's VersionAbstract
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 C-14-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 C-14 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). The geomagnetic spikes rise and fall over a period of less than 30 years and are associated with VADM fluctuations of at least 70 ZAm2. Thus, the Levantine archaeomagnetic record places new constraints on maximum geomagnetic intensity as well as for its rate of change. Yet, it is not clear whether the geomagnetic spikes are local non-dipolar features or a geomagnetic dipolar phenomenon. (C) 2010 Elsevier B.V. All rights reserved.
Shaar R, Ron N, Tauxe L, Kessel R, Agnon A. Paleomagnetic field intensity derived from non-SD: Testing the Thellier IZZI technique on MD slag and a new bootstrap procedure. Earth and Planetary Science Letters [Internet]. 2011;310 :213-224. Publisher's VersionAbstract
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 (B-TRM) and the field in which the NRM was acquired (B-NRM). 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 B-TRM is parallel and significantly stronger than B-NRM. 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 B-TRM and B-NRM. 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. (C) 2011 Elsevier BM. All rights reserved.
Shaar R, Ben-Yosef E. Earth's magnetic field: new discoveries from archaeomagnetic study of copper slag from Arava valley. In: Starinsky A Melakh Haaretz (Salt of the Earth). Vol. 5. Jerusalem: The Hebrew University Magnes Press Ltd. ; 2010.
Ebert Y, Kessel R, Shaar R, Agnon A, Ron H. Petrology and rock magnetism of the gabbro of Troodos ophiolite. Physics of the Earth and Planetary Interiors [Internet]. 2010;183 :413-420. Publisher's VersionAbstract
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, <350 degrees C, and the high component, <580 degrees C, are the contribution of the pyrrhotite and magnetite, respectively. 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 degrees 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. (C) 2010 Elsevier B.V. All rights reserved.
Shaar R, Ron H, Tauxe L, Kessel R, Agnon A, Ben-Yosef E, Feinberg JM. Testing the accuracy of absolute intensity estimates of the ancient geomagnetic field using copper slag material. Earth and Planetary Science Letters [Internet]. 2010;290 :201-213. Publisher's VersionAbstract
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 1721 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 submicronelongated 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 nonlinear 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 similar to 5% after applying the required corrections. (C) 2009 Elsevier B.V. All rights reserved.
Matmon A, Katz O, Shaar R, Ron H, Porat N, Agnon A. Timing of relay ramp growth and normal fault linkage, Upper Galilee, northern Israel. Tectonics [Internet]. 2010;29. Publisher's VersionAbstract
We investigate a kilometer-scale steep relay ramp structure in the Galilee, northern Israel. The data indicate an asymmetric fold above the buried tip of the southern boundary fault as well as at the lower part of the ramp. Structural analysis suggests that the buried tip of the bounding fault is only a few hundred meters below the surface. A sequence of colluvial wedges, exposed at the base of the southern flank of the relay ramp over the blind tip of the boundary fault, presents successively decreasing angles from 67 at the bottom to 30 at the top. The lower three wedges rest at angles greater than the angle of repose, suggesting tilting during and after deposition. We suggest that the increasing angle of the wedges is the result of the development of the monocline over the buried tip of the boundary fault, which in turn reflects the development of the entire relay ramp structure. Paleomagnetic measurements test this hypothesis. Samples from three of the four lower wedges and the carbonate fill in the fractured bedrock yield northerly declination and positive inclination. Inclination anomalies are 19.4 +/- 3.3 and 11.3 +/- 1.8 for the fractured bedrock and the wedges, respectively. These anomalies reflect part of the tilt. Paleomagnetic measurements and optically stimulated luminescence (OSL) constrain the time of initial deformation of the monocline above the southern boundary fault between 780 ka and 461 +/- 75 ka (OSL age of wedge 1) and the termination of deformation to 176 +/- 22 ka (OSL age of wedge 5). These bounding ages suggest tilting rates that range between 1 degrees/16 kyr and 1 degrees/8 kyr. Termination of tilting along the southern flank of the relay ramp and the fact that the tip of the buried normal fault is only a few hundreds of meters below the surface suggest that the relay ramp may be close to being breached by a connecting fault. Citation: Matmon, A., O. Katz, R. Shaar, H. Ron, N. Porat, and A. Agnon (2010), Timing of relay ramp growth and normal fault linkage, Upper Galilee, northern Israel, Tectonics, 29, TC2016, doi: 10.1029/2009TC002510.
Ben-Yosef E, Tauxe L, Levy TE, Shaar R, Ron H, Najjar M. Geomagnetic intensity spike recorded in high resolution slag deposit in Southern Jordan. Earth and Planetary Science Letters [Internet]. 2009;287 :529-539. Publisher's VersionAbstract
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 mu T (compared to the current field strength of 44 mu 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. (C) 2009 Elsevier B.V. All rights reserved.