Ioan Lascu

FORCulator: A micromagnetic tool for simulating first‐order reversal curve diagrams

We describe a method for simulating first-order reversal curve (FORC) diagrams of interacting single-domain particles. Magnetostatic interactions are calculated in real space, allowing simulations to be performed for particle ensembles with arbitrary geometry. For weakly interacting uniaxial particles, the equilibrium magnetization at each field step is obtained by direct solution of the Stoner-Wohlfarth model, assuming a quasi-static distribution of interaction fields. For all other cases, the equilibrium magnetization is calculated using an approximate iterated solution to the Landau-Lifshitz-Gilbert equation. Multithreading is employed to allow multiple curves to be computed simultaneously, enabling FORC diagrams to be simulated in reasonable time using a standard desktop computer. Statistical averaging and post processing lead to simulated FORC diagrams that are comparable to their experimental counterparts. The method is applied to several geometries of relevance to rock and environmental magnetism, including densely packed random clusters and partially collapsed chains. The method forms the basis of FORCulator, a freely available software tool with graphical user interface that will enable FORC simulations to become a routine part of rock magnetic studies.

Magnetic unmixing of first-order reversal curve diagrams using principal component analysis

We describe a quantitative magnetic unmixing method based on principal component analysis (PCA) of first-order reversal curve (FORC) diagrams. For PCA we resample FORC distributions on grids that capture diagnostic signatures of single-domain (SD), pseudo-single-domain (PSD), and multi-domain (MD) magnetite, as well as of minerals such as hematite. Individual FORC diagrams are recast as linear combinations of end-member (EM) FORC diagrams, located at user-defined positions in PCA space. The EM selection is guided by constraints derived from physical modeling and imposed by data scatter. We investigate temporal variations of two EMs in bulk North Atlantic sediment cores collected from the Rockall Trough and the Iberian Continental Margin. Sediments from each site contain a mixture of magnetosomes and granulometrically distinct detrital magnetite. We also quantify the spatial variation of three EM components (a coarse silt-sized MD component, a fine silt-sized PSD component, and a mixed clay-sized component containing both SD magnetite and hematite) in surficial sediments along the flow path of the North Atlantic Deep Water (NADW). These samples were separated into granulometric fractions, which helped constrain EM definition. PCA-based unmixing reveals systematic variations in EM relative abundance as a function of distance along NADW flow. Finally, we apply PCA to the combined dataset of Rockall Trough and NADW sediments, which can be recast as a four-EM mixture, providing enhanced discrimination between components. Our method forms the foundation of a general solution to the problem of unmixing multi-component magnetic mixtures, a fundamental task of rock magnetic studies. This article is protected by copyright. All rights reserved.

Age of the Laschamp excursion determined by U-Th dating of a speleothem geomagnetic record from North America

This project was funded by NSF-EAR grant 1316385, a University of Minnesota McKnight Land Grant Professorship to JMF, and ERC grant 320750. Confocal microscopy was performed at the University of Minnesota Imaging Centers. We are grateful to John Geissman, Brad Singer, and James Channell for their constructive reviews. This is Institute for Rock Magnetism contribution 1506.

Magnetic record of deglaciation using FORC‐PCA, sortable‐silt grain size, and magnetic excursion at 26 ka, from the Rockall Trough (NE Atlantic)

Research supported by US NSF grants 0850413 and 1014506, and the European R12esearch Council under the European Unions Seventh Framework Programme (FP/2007-2013)/ERC Grant Agreement No. 320750. The UK NERC and BGS funded the recovery of Core MD04-2822.

Deconvolution of u channel magnetometer data: Experimental study of accuracy, resolution, and stability of different inversion methods

We explore the effects of sampling density, signal/noise ratios, and position-dependent measurement errors on deconvolution calculations for u channel magnetometer data, using a combination of experimental and numerical approaches. Experiments involve a synthetic sample set made by setting hydraulic cement in a 30-cm u channel and slicing the hardened material into ∼2-cm lengths, and a natural lake sediment u channel sample. The cement segments can be magnetized and measured individually, and reassembled for continuous u channel measurement and deconvolution; the lake sediment channel was first measured continuously and then sliced into discrete samples for individual measurement. Each continuous data set was deconvolved using the ABIC minimization code of Oda and Shibuya (1996) and two new approaches that we have developed, using singular-value decomposition and regularized least squares. These involve somewhat different methods to stabilize the inverse calculations and different criteria for identifying the optimum solution, but we find in all of our experiments that the three methods converge to essentially identical solutions. Repeat scans in several experiments show that measurement errors are not distributed with position-independent variance; errors in setting/determining the u channel position (standard deviation ∼0.2 mm) translate in regions of strong gradients into measurement uncertainties much larger than those due to instrument noise and drift. When we incorporate these depth-dependent measurement uncertainties into the deconvolution calculations, the resulting models show decreased stability and accuracy compared to inversions assuming depth-independent measurement errors. The cement experiments involved varying directions and uniform intensities downcore, and very good accuracy was obtained using all of the methods when the signal/noise ratio was greater than a few hundred and the sampling interval no larger than half the length scale of magnetization changes. Addition of synthetic noise or reduction of sampling density decreased the resolution and accuracy of all the methods equally. The sediment-core experiment involved uniform (axial) magnetization direction and strongly varying intensities downcore. Intensity variations are well resolved and directions are accurate to within about 5 degrees, with errors attributable to omission and/or inaccurate calibration of cross terms in the instrument response function.

Elastic and magnetoelastic relaxation behaviour of multiferroic (ferromagnetic + ferroelectric + ferroelastic) Pb(Fe0.5Nb0.5)O3 perovskite.

Resonant Ultrasound Spectroscopy has been used to characterize elastic and anelastic anomalies in a polycrystalline sample of multiferroic Pb(Fe(0.5)Nb(0.5))O(3) (PFN). Elastic softening begins at ~550 K, which is close to the Burns temperature marking the development of dynamical polar nanoregions. A small increase in acoustic loss at ~425 K coincides with the value of T(*) reported for polar nanoregions starting to acquire a static or quasi-static component. Softening of the shear modulus by ~30-35% through ~395-320 K, together with a peak in acoustic loss, is due to classical strain/order parameter coupling through the cubic → tetragonal → monoclinic transition sequence of ferroelectric/ferroelastic transitions. A plateau of high acoustic loss below ~320 K is due to the mobility under stress of a ferroelastic microstructure but, instead of the typical effects of freezing of twin wall motion at some low temperature, there is a steady decrease in loss and increase in elastic stiffness below ~85 K. This is attributed to freezing of a succession of strain-coupled defects with a range of relaxation times and is consistent with a report in the literature that PFN develops a tweed microstructure over a wide temperature interval. No overt anomaly was observed near the expected Néel point, ~145 K, consistent with weak/absent spin/lattice coupling but heat capacity measurements showed that the antiferromagnetic transition is actually smeared out or suppressed. Instead, the sample is weakly ferromagnetic up to ~560 K, though it has not been possible to exclude definitively the possibility that this could be due to some magnetic impurity. Overall, evidence from the RUS data is of a permeating influence of static and dynamic strain relaxation effects which are attributed to local strain heterogeneity on a mesoscopic length scale. These, in turn, must have a role in determining the magnetic properties and multiferroic character of PFN.

Variable ecosystem response to climate change during the Holocene in northern Minnesota, USA

Both terrestrial and aquatic ecosystems should respond to abrupt climate changes such as those that affected the midcontinent of North America during the Holocene. A variety of paleorecords indicate that the onset of prolonged dry conditions in this region occurred as rapidly as 300 yr during the early Holocene, with a subsequent increase in moisture occurring rapidly in the late Holocene. Here, we report a 9500 yr multiproxy sediment record from Deming Lake, Minnesota, USA, that demonstrates only subtle dynamics during rapid climate changes that caused vegetation in the catchment to shift among pine forest, open grassland, and deciduous forest. The most substantial changes in ecosystem properties immediately followed deglaciation of the landscape, formation of the lake, and initial development of pine forests. In contrast, there were only muted responses to pronounced mid-Holocene climate changes that caused vegetation in the catchment to switch from pine forest to open grassland (ca. 8000 yr B.P.), and then deciduous forest (ca. 5400 yr B.P.). The flux of organic and inorganic terrigenous material, the processing of carbon, and catchment erosion changed rapidly during the most recent shift at 3300 yr B.P. to the modern pine forest. These changes coincided with the onset of meromictic conditions that influenced internal lake dynamics. However, the terrestrial influence on the lake sedimentary record gradually diminished over time, indicating a trajectory of increasing catchment stabilization that was relatively impervious to dramatic regional climate changes. The relative complacency of the Deming Lake record during the late Holocene indicates relative resistance to abrupt climate change at later stages of ecosystem development.

A comparison of magnetic susceptibility measurement techniques and ferrimagnetic component analysis from recent sediments in Lake Pepin (USA)

Abstract Magnetic susceptibility (MS) is commonly measured on lake sediments and used as a proxy for clastic input and soil erosion, for correlating among cores in the same lake basin and aligning successive overlapping drives of the same core. There are several common techniques for measuring MS, each with its own advantages. Here we compare three such techniques measured on a sediment core from Lake Pepin: (a) loop-sensor MS logging on wet sediment of the intact core; (b) point-sensor MS logging on wet sediment of a split (lengthwise) core; and (c) discrete MS measurements of dried subsamples using a susceptibility bridge. To obtain further information about the origin of downcore MS variability, additional magnetic measurements were performed on discrete samples to model ferrimagnetic sedimentary components. Overall trends and individual features in MS curves agree reasonably well between techniques; however, the amplitude of local minima and maxima varies according to the technique used. All three MS techniques captured distinct events c. 1900 and 1940, attributed to increases in allochthonous ferrimagnetic components. The ferrimagnetic particle flux has declined over the past half-century despite increasing sediment accumulation in Lake Pepin, suggesting a possible shift in sediment sources from fields to stream banks.

Elastic and anelastic relaxation behaviour of perovskite multiferroics II: PbZr0.53Ti0.47O3 (PZT)–PbFe0.5Ta0.5O3 (PFT)

Elastic and anelastic properties of ceramic samples of multiferroic perovskites with nominal compositions across the binary join PbZr0.53Ti0.47O3–PbFe0.5Ta0.5O3 (PZT–PFT) have been assembled to create a binary phase diagram and to address the role of strain relaxation associated with their phase transitions. Structural relationships are similar to those observed previously for PbZr0.53Ti0.47O3–PbFe0.5Nb0.5O3 (PZT–PFN), but the magnitude of the tetragonal shear strain associated with the ferroelectric order parameter appears to be much smaller. This leads to relaxor character for the development of ferroelectric properties in the end member PbFe0.5Ta0.5O3. As for PZT–PFN, there appear to be two discrete instabilities rather than simply a reorientation of the electric dipole in the transition sequence cubic–tetragonal–monoclinic, and the second transition has characteristics typical of an improper ferroelastic. At intermediate compositions, the ferroelastic microstructure has strain heterogeneities on a mesoscopic length scale and, probably, also on a microscopic scale. This results in a wide anelastic freezing interval for strain-related defects rather than the freezing of discrete twin walls that would occur in a conventional ferroelastic material. In PFT, however, the acoustic loss behaviour more nearly resembles that due to freezing of conventional ferroelastic twin walls. Precursor softening of the shear modulus in both PFT and PFN does not fit with a Vogel–Fulcher description, but in PFT there is a temperature interval where the softening conforms to a power law suggestive of the role of fluctuations of the order parameter with dispersion along one branch of the Brillouin zone. Magnetic ordering appears to be coupled only weakly with a volume strain and not with shear strain but, as with multiferroic PZT–PFN perovskites, takes place within crystals which have significant strain heterogeneities on different length scales.

Elastic and anelastic relaxation behaviour of perovskite multiferroics II: PbZr0.53Ti0.47O3 (PZT)-PbFe0.5Ta0.5O3 (PFT)

Elastic and anelastic properties of ceramic samples of multiferroic perovskites with nominal compositions across the binary join PbZr0.53Ti0.47O3–PbFe0.5Ta0.5O3 (PZT–PFT) have been assembled to create a binary phase diagram and to address the role of strain relaxation associated with their phase transitions. Structural relationships are similar to those observed previously for PbZr0.53Ti0.47O3–PbFe0.5Nb0.5O3 (PZT–PFN), but the magnitude of the tetragonal shear strain associated with the ferroelectric order parameter appears to be much smaller. This leads to relaxor character for the development of ferroelectric properties in the end member PbFe0.5Ta0.5O3. As for PZT–PFN, there appear to be two discrete instabilities rather than simply a reorientation of the electric dipole in the transition sequence cubic–tetragonal–monoclinic, and the second transition has characteristics typical of an improper ferroelastic. At intermediate compositions, the ferroelastic microstructure has strain heterogeneities on a mesoscopic length scale and, probably, also on a microscopic scale. This results in a wide anelastic freezing interval for strain-related defects rather than the freezing of discrete twin walls that would occur in a conventional ferroelastic material. In PFT, however, the acoustic loss behaviour more nearly resembles that due to freezing of conventional ferroelastic twin walls. Precursor softening of the shear modulus in both PFT and PFN does not fit with a Vogel–Fulcher description, but in PFT there is a temperature interval where the softening conforms to a power law suggestive of the role of fluctuations of the order parameter with dispersion along one branch of the Brillouin zone. Magnetic ordering appears to be coupled only weakly with a volume strain and not with shear strain but, as with multiferroic PZT–PFN perovskites, takes place within crystals which have significant strain heterogeneities on different length scales.