Palaeointenisty measurements are vital for understanding the deep Earth and the evolution of the geomagnetic field, however, they suffer from non-ideal behaviour that compromises the fidelity of the results. In an attempt to screen for the non-ideal behaviour that plagues palaeointensity experiments, selection criteria are used. One such non-ideal behaviour is partial thermoremanent magnetisation (pTRM) tails, which in the IZZI protocol result in a zig-zagging pattern in the Arai plot. While pTRM tail checks can be carried out, these are not routinely used in IZZI experiments. We set out a new framework with which to create, develop and re-assess selection criteria before applying it to our proposed new selection statistic, Ziggie, which is specific to detecting zig-zagging behaviour. Ziggie meets all the requirements set out; that is, it is able to objectively and robustly identify the target behaviour, consistently provide a quantifiable assessment of the target behaviour, and is predominately sensitive to the target behaviour while being free of undue influences. Other pre-existing statistics for zig-zagging in Arai plots, Z*, IZZI_MD, β and SCAT, are also tested, but fail to meet all requirements. We find that Ziggie should be adopted in selection criteria sets for the IZZI protocol with a criterion of Ziggie ≤ 0.1.

First-order reversal curve (FORC) diagrams are a standard rock magnetic tool for analyzing bulk magnetic hysteresis behaviors, which are used to estimate the magnetic mineralogies and magnetic domain states of grains within natural materials. However, the interpretation of FORC distributions is challenging due to complex domain-state responses, which introduce well-documented uncertainties and subjectivity. Here, we propose a neural network algorithm (FORCINN) to invert the size and aspect ratio distribution from measured FORC data. We trained and tested the FORCINN model using a dataset of synthetic numerical FORCs for single magnetite grains with various grain-sizes (45-400 nm) and aspect ratios (oblate and prolate grains). In addition to successfully testing FORCINN against synthetic datasets, we also tested FORCINN against FORC data measured on natural samples with accurately determined grain-size and aspect ratio distributions. FORCINN was found to provide good estimates of the grain-size distributions for basalt samples and marine sediments.

The minerals carrying the magnetic remanence in geological samples are commonly a solid solution series of iron-titanium spinels known as titanomagnetites. Despite the range of compositions within this series, micromagnetic studies that characterize the magnetic domain structures present in these minerals have typically focused on magnetite. No studies systematically comparing the domain-states present in titanomagnetites have been undertaken since the discovery of the single vortex (SV) structure and the advent of modern micromagnetism. The magnetic properties of the titanomagnetite series are known to vary with composition, which may influence the domain states present in these minerals, and therefore the magnetic stability of the samples bearing them. We present results from micromagnetic simulations of titanomagnetite ellipsoids of varying shape and composition to find the size ranges of the single domain (SD) and SV structures. These size ranges overlap, allowing for regions where the SD and SV structures are both available. These regions are of interest as they may lead to magnetic instability and “pTRM tails’ in paleointensity experiments. We find that although this SD+SV zone occupies a narrow range of sizes for equidimensional magnetite, it is widest for intermediate (TM30-40) titanomagnetite compositions, and increases for both oblate and prolate particles, with some compositions and sizes having an SD+SV zone up to 100s of nm wide. Our results help to explain the prevalence of pTRM tail-like behavior in paleointensity experiments. They also highlight regions of particles with unusual domain states to target for further investigation into the definitive mechanism behind paleointensity failure.

Micromagnetic modelling allows the systematic study of the effects of particle size and shape on the first-order reversal curve (FORC) magnetic hysteresis response for magnetite particles in the single-domain (SD) and pseudo-single domain (PSD) particle size range. The interpretation of FORCs, though widely used, has been highly subjective.  Here, we use micromagnetics to model randomly oriented distributions of particles to allow more physically meaningful interpretations.  We show that one commonly found type of PSD particle - namely single vortex (SV) particles - has far more complex signals than SD particles, with multiple peaks and troughs in the FORC distribution, where the peaks have higher switching fields for larger SV particles. Particles in the SD to SV transition zone have the lowest switching fields. Symmetrical and prolate particles display similar behavior, with distinctive peaks forming near the vertical axis of the FORC diagram. In contrast, highly oblate particles produce `butterfly' structures, suggesting that these are potentially diagnostic of particle morphology. We also consider FORC diagrams for distributions of particle sizes and shapes and produce an online application that users can use to build their own FORC distributions.  There is good agreement between the model predictions for distributions of particle sizes and shapes, and the published experimental literature.

Magnetic hysteresis measurements are routinely made in the Earth and planetary sciences to identify geologically meaningful  magnetic recorders, and to study variations in present and past environments. Interpreting magnetic hysteresis data in terms of domain state (particle size)  and paleomagnetic stability are major motivations behind undertaking these measurements, but the interpretations remain  fraught with challenges and ambiguities. To shed new light on these ambiguities, we have undertaken a systematic micromagnetic study to quantify the magnetic hysteresis behavior of room-temperature magnetite as a function of particle size (50-195 nm; equivalent spherical volume diameter) and shape (oblate, prolate and equant); our models span uniformly magnetized single domain (SD) to non-uniformly magnetized single vortex (SV) states. Within our models the reduced magnetization  marks a clear boundary between SD (≥0.5) and SV (<0.5) magnetite. We further identify particle sizes and shapes with unexpectedly low coercivity and coercivity of remanence. These low coercivity regions correspond to magnetite particles that typically have multiple possible magnetic domain states, which has been previously linked to a zone of unstable magnetic recorders. Of all hysteresis parameters investigated, transient hysteresis is most sensitive to particles that exhibit such domain state multiplicity, leading us to suggest that transient behavior be more routinely measured during rock magnetic investigations.

The ability of rocks to hold a reliable record of the ancient geomagnetic field depends on the structure and stability of magnetic domain-states contained within the rock’s magnetic particles. In paleomagnetic studies, the Day plot is an easily constructed graph of magnetic hysteresis parameters that is frequently used (and mis-used) to estimate the likely magnetic recording stability of samples. Often samples plot in the region of the Day plot attributed to so-called pseudo-single-domain (PSD) particles with little under standing of the implications for domain-states or recording fidelity. Here we use micromagnetic models to explore the hysteresis parameters of magnetite particles with idealized prolate and oblate truncated-octahedral geometries containing single domain (SD), single-vortex (SV) and occasionally multi-vortex (MV) states. We show that these do main states exhibit a well-defined trend in the Day plot that extends from the SD region well into the multi-domain (MD) region, all of which are likely to be stable remanence carriers. We suggest that although the interpretation of the Day plot and its vari33 ants might be subject to ambiguities, if the magnetic mineralogy is known, it can still provide some useful insights about paleomagnetic specimens’ dominant domain state, average particle sizes and, consequently, their paleomagnetic stability.

Thermal demagnetization furnaces are routine facilities that underpin countless paleomagnetic studies by allowing the progressive removal of naturally acquired magnetic remanence or the imparting of well controlled laboratory magnetization. The ideal thermal demagnetizer should maintain “zero” magnetic field during thermal treatments. However, magnetic field noise, including residual magnetic fields of material used to construct the furnace and induced fields caused by the heating current in the furnace are always present. As technology advances allowing the measurement of ever weaker magnetic remanences, it is essential that high-performance demagnetization furnaces are developed to reduce these sources of magnetic field noise. By combining efficient demagnetization of shielding and a new structure of heating wire, we have developed a new demagnetization furnace with low magnetic field noise. Repeated progressive thermal demagnetization experiments using specimens that were previously completely thermally demagnetized above their Curie temperature were carried out to explore the effects of fields within various types of furnace during demagnetization. These experiments confirm that magnetic field noise in various furnaces can have an observable and detrimental impact on demagnetization behavior and that this is reduced with our new design. The new heating element design and procedure for reducing magnetic field noises represent a significant improvement in the design of thermal demagnetizers and allows for extremely weak specimens to be successfully measured.

Understanding the evolution of Earth’s magnetic field can provide insights into core processes and can constrain plate tectonics and atmospheric shielding. The absolute paleointensity database PINT provides a curated repository of site mean, (i.e., cooling unit), estimates of the strength of the magnetic field. We present a minor update to the PINT database to version 8.1.0 by adding 248 records from 31 studies. The PINT database is used to define a continuous model of the dipole field, using an approach combining non-parametric and Monte Carlo resampling termed MCADAM. Three dipole field strength models spanning 50 ka to 3.7-4.2 Ga (MCADAM.1a-c) are presented, reflecting three tiers of increasingly more stringent data selection. The MCADAM models allow for the estimation of the magnetic standoff distance, constraining the shielding of Earth’s atmosphere against solar wind erosion provided by the geodynamo.

Magnetic hysteresis loops are an important tool in theoretical and applied rock magnetism with applications to paleointensities, paleoenvironmental analysis, and tectonic studies, among many others. Hence, information derived from these data is amongst the most ubiquitous rock magnetic data used by the Earth science community. Despite their prevalence, there are no general guidelines to aid scientists in obtaining the best possible data and no widely available software to allow the efficient analysis of hysteresis loop data using the most advanced and appropriate methods. Here we provide an outline of detrimental factors and simple approaches to measuring better hysteresis loops as well as introducing a new MATLAB software package called Hysteresis Loop analysis box (HystLab) for processing and analyzing loop data. This graphical user interface software is capable of reading the wide range of data formats that are generated by the multiple types of equipment typically used to measure hysteresis loops. HystLab provides an easy-to-use interface allowing users to visualize their data and perform advanced processing, including loop centering, drift correction, linear and approach to saturation high-field slope corrections, as well as loop fitting to improve the results from noisy specimens. A large number of hysteresis loop properties and statistics are calculated by HystLab and can be exported to text files for further analysis. All plots generated by HystLab are customizable and user preferences can be saved for future use. In addition, all plots can be exported to encapsulated postscript (EPS) files that are publication ready with little or no adjustment, greatly enhancing workflow productivity when processing and analyzing large data sets. HystLab is freely available for download at https://github.com/greigpaterson/HystLab and in combination with our simple measurement guide should help the paleo- and rock magnetic communities get the most from their hysteresis data.