MSP-tool: A VBA-based software tool for the analysis of multispecimen paleointensity data

The multispecimen protocol (MSP) is a method to estimate the Earth's magnetic field's past strength from volcanic rocks or archeological materials. By reducing the amount of heating steps and aligning the specimens parallel to the applied field, thermochemical alteration and multi-domain effects are minimized. We present a new software tool, written for Microsoft Excel 2010 in Visual Basic for Applications (VBA), that evaluates paleointensity data acquired using this protocol. In addition to the three ratios (standard, fraction-corrected, and domain-state-corrected) calculated following Dekkers and Böhnel (2006) and Fabian and Leonhardt (2010) and a number of other parameters proposed by Fabian and Leonhardt (2010), it also provides several reliability criteria. These include an alteration criterion, whether or not the linear regression intersects the y axis within the theoretically prescribed range, and two directional checks. Overprints and misalignment are detected by isolating the remaining natural remanent magnetization (NRM) and the partial thermoremanent magnetization (pTRM) gained and comparing their declinations and inclinations. The NRM remaining and pTRM gained are then used to calculate alignment-corrected multispecimen plots. Data are analyzed using bootstrap statistics. The program was tested on lava samples that were given a full TRM and that acquired their pTRMs at angles of 0, 15, 30, and 90° with respect to their NRMs. MSP-Tool adequately detected and largely corrected these artificial alignment errors

The performance of various palaeointensity techniques as a function of rock magnetic behaviour - A case study for La Palma

Three different palaeointensity methods were applied to six historical and three carbon-dated flows from the island of La Palma (Spain); in total fifteen sites were processed. The two 20th-century flows were sampled at multiple locations as their obtained directions and intensities can be compared directly to those from the International Geomagnetic Reference Field (IGRF). After determination of the declinations and inclinations of the natural remanent magnetisation (NRM) by thermal and alternating-field demagnetisation, the samples were subjected to standard rock magnetic analyses to determine their Curie and alteration temperatures. Based on these characteristics, the sites were allocated to one of four rock magnetic groups labelled L*, L, C, and H, a division primarily based on the temperature-dependent behaviour of the low-field susceptibility that has been used in studies of other volcanic edifices. Scanning electron microscope (SEM) observations revealed little oxidation and exsolution (oxidation classes I to III). Palaeointensities were determined using the classic Thellier-Thellier method (Aitken and IZZI protocols), the microwave method and the domain-state-corrected multispecimen method. Thellier-Thellier and microwave results were analysed using the ThellierTool A and B sets of selection criteria as modified by Paterson et al. (2014). Their combined success rate was around 40%. Of the eight IGRF sites, two yielded average intensities within 10% of the IGRF value. For the microwave method, three sites reproduced the IGRF intensity within 10%. In the domain-state-corrected multispecimen protocol, just one site (site 9, 1971) passed the 'ARM-test' (applied in retrospect) and showed less than 3% progressive alteration. Its multispecimen result reproduced the palaeofield within error. The other IGRF sites over- or underestimated the palaeofield by up to 50%. The seven older sites produced plausible palaeointensities, generally within a few μT of model data, and if multiple methods were successful, the results were within error of each other. For all three PI methods, it seems that sites with low Curie temperatures (<150°C; group L*), are more likely to pass all selection criteria while substantially over- or underestimating the palaeofield. It is hypothesised that time-dependent processes after cooling of the lava would be a prime reason for this discrepancy: PI experiments with a laboratory thermoremanent magnetisation (TRM), imparted at a temperature above the site's dominant Curie temperature but below its alteration temperature, yielded the correct intensity of the laboratory-imparted TRM. When two or three methods agree to within a few μT, the obtained palaeointensity is close to the palaeofield. Multi-method consistency provides an additional palaeointensity reliability check