Paleointensity Experiments Research Articles

On the early postclassic collapse of Plazuelas (Northern Border of Mesoamerica): Archaeomagnetic evidence based on full geomagnetic vectors

We report a detailed rock-magnetic, archaeomagnetic and absolute intensity investigation on burned in situ remains of the Plazuelas archaeological site (Northern Border of Mesoamerica). Fifty-six specimens were subjected to step-wise alternating field demagnetizations to retrieve characteristic remanence directions. All samples yielded evidence of stable, single-component magnetization carried by magnetite or Ti-poor titanomagnetite. Sixty-four specimens were selected for paleointensity experiments following the Thellier double heating technique, including remanence anisotropy correction. The archaeomagnetic dating, based on full geomagnetic vectors, was performed using the global geomagnetic models and the regional Paleosecular Variation Curves available for Mesoamerica. The most probable age intervals obtained for the burned floor of Site 1 (Casas Tapadas North) are within the interval 972AD – 744 AD. In the case of the burned floor of Site 2 (Casas Tapadas West), the ages obtained fall within 971 AD - 752 AD. In the case of the burned floor of Site 3 (South Sector), the ages obtained fall within the interval 752 AD – 969 AD. The calculation of the most probable age intervals, carried out using the Bootstrap resampling method, provided the dates of 940 ± 6 AD, 896 ± 14 AD and 883 ± 14 AD, respectively. These ages agree well with various archaeological hypotheses on the sudden collapse of Plazuelas as a result of intentional destruction, mutilation and generalized firing due to the local war and/or social unrest.

First Full‐Vector Archeomagnetic Data From Central Asia (3 BCE to 15 CE Centuries): Evidence for a Large Non‐Dipole Field Contribution Around the First Century BCE

AbstractUnraveling the short‐term behavior of the Earth's past geomagnetic field at regional scales is crucial for understanding its global behavior and, thus, the dynamics of the deep Earth. In this context, obtaining accurate full‐vector geomagnetic field records from regions where archeomagnetic data are absent becomes essential. Here, we present the first full‐vector archeomagnetic data from Central Asia, derived from the analysis of nine archeological kilns sampled in South Uzbekistan, dating back to the period between 200 BCE and 1429 CE. To obtain these new data, we conducted thermal and alternating field demagnetization procedures, along with Thellier‐Thellier paleointensity experiments, including partial thermoremanent magnetization checks, thermoremanent magnetization anisotropy and cooling rate corrections. The comparison between the new data, previous selected data from Central Asia, and available global models reveals important differences between approximately 400 BCE and 400 CE, especially concerning the geomagnetic field intensity element. In order to investigate this in detail, we have developed a regional update of the SHAWQ global models family by incorporating, for the first time, high‐quality data from Central Asia. The results suggest that this deviation is linked to non‐dipolar sources of the geomagnetic field in Central Asia reaching a maximum contribution around the first century BCE. According to the updated global paleoreconstruction, this non‐dipole feature, manifested at the Earth's surface as low intensities, is associated with the presence of a reversed flux patch at the core‐mantle boundary beneath this region.

Paleomagnetism, rock magnetism and age determination of effusive and explosive Holocene volcanism in the Momotombo-Managua-Masaya region, Nicaragua

The Central America Volcanic Arc shows an intense and continuous volcanic activity spanning the last 10 ky because of the subduction processes parallel to the Middle America Trench. Several volcanic events have occurred during the Plio-Holocene, along the coastal Nicaraguan region, characterized either by explosive or effusive activity. Volcanic deposits were emplaced at high temperatures and have similar basaltic-andesitic composition suggesting a long-lived magma chamber but emplaced and cooled by different mechanisms. We present paleodirectional and paleointensity results of four lava flows, and the paleointensity recorded on different pyroclastic deposits in the Managua-Masaya region, western Nicaragua. For each, we characterize the rock magnetic properties as a function of temperature and magnetic field. Absolute paleointensity experiments estimated a mean range from 38.6 to 52.9 μT; and a mean virtual axial dipole moment from 9.4 × 1022 Am2 to 12.8 × 1022 Am2. Two paleomagnetic ages were determined for two cooling units: (1) the Masaya Caldera Walls with an age range of 8580–8535 BP, and (2) a dating for a lava flow from Momotombo Volcanic Complex estimated of 1500–1600 CE. These results are consistent with the expected values according to the actual paleosecular variation global models, and with the regional stratigraphy and historical ages. The high reliability of paleointensity determination from these pyroclastic deposits and lava flows suggest that both are suitable for paleomagnetic dating. Finally, an accurate chronology of the volcanic events, will improve the understanding of the evolution of the volcanism and their implications for volcanic risk on the region. This knowledge will help to decrease the considerable hazard for the population on the surrounding areas.

Using Preisach Theory to Evaluate Chemical Remanent Magnetization and Its Behavior During Thellier‐Thellier‐Coe Paleointensity Experiments

AbstractThe behavior of grain‐growth chemical remanent magnetizations (gCRM) are investigated for different coercivity and magnetostatic‐interaction‐field distributions and acquisition conditions using a thermally activated Preisach model for assemblages of interacting single‐domain grains. A new growth‐rate dependent equation was derived, from which it was found that gCRM intensity is over 10% more sensitive to growth rate than previously modeled. We compare the gCRM results to the behavior of thermoremanences (TRM). gCRMs are two times more sensitive to changes in coercivity distribution, whereas TRMs are four times more sensitive to changes in magnetostatic interactions. The Thellier‐Thellier‐Coe paleointensity protocol was simulated in Preisach space, and gCRMs were found to produce concave‐up Arai plots with pTRM checks which plot to the left of the Arai plot and positive partial‐TRM tail checks that increase with magnetostatic interactions. This often leads to the failure of selection criteria, but high‐temperature segments can pass the criteria for weakly interacting gCRMs; these estimates can underestimate the field by up to 66%.

Rock Magnetic Characterization of Returned Samples From Asteroid (162173) Ryugu: Implications for Paleomagnetic Interpretation and Paleointensity Estimation

AbstractIn this study, systematic rock magnetic measurements and saturation isothermal remanent magnetization (SIRM) paleointensity calibration experiments were conducted for the returned samples from C‐type asteroid (162173) Ryugu and two carbonaceous chondrites (Orgueil and Tagish Lake) to evaluate the remanence carriers of the Ryugu sample and its ability as a paleomagnetic recorder. Our magnetic measurements show that Ryugu samples exhibit signatures for framboidal magnetite, coarse‐grained magnetite, and pyrrhotite, and that framboidal magnetite is the dominant remanence carrier of Ryugu samples in the middle‐coercivity range. The SIRM paleointensity constant was obtained for two Ryugu samples, and the median value was 3,318 ± 1,038 μT, which is close to the literature's value based on the average among magnetite, titanomagnetite, pyrrhotite, and FeNi alloys and is widely used for SIRM paleointensity experiments. The paleointensity values estimated using the obtained SIRM paleointensity constant indicate a strong magnetic field of the protoplanetary disk, suggesting that Sun's protoplanetary disk existed at the disk location of Ryugu's parent planetesimal when framboidal magnetite precipitated from the aqueous fluid.

From Toltecs to Franciscans: The history of El Cerrito pyramid (Middle Mesoamerica) revealed by the fire

In order to obtain the first absolute chronology data, a detailed rock-magnetic, archaeomagnetic and radiocarbon investigation was performed on some in situ features of the El Cerrito pyramid. Archaeomagnetic sampling was carried out on two rectangular fire pits and one burned soil located at the piedmont of the pyramid. Marked irreversibility between heating and cooling cycles and consequently thermal instability is observed on continuous thermomagnetic experiments which discouraged to perform paleointensity experiments. Alternating field demagnetization revealed essentially single component magnetization with negligible secondary overprint probably of viscous origin. Archaeomagnetic age intervals were obtained using available global models and local curves with rather concordant results. Calibrated radiocarbon age agrees well with archaeomagnetic date in one fire pit analyzed while slight differences are observed for the burned soil and the second fire pit. This may be interpreted as reuse and continuous use of these ceremonial areas involving fire as an essential element. The evidence of possible colonial occupation indicates that El Cerrito persisted as a ceremonial center over the centuries, as also attested by Franciscan historical documents.

Absolute Paleointensity Experiments on Aged Thermoremanent Magnetization: Assessment of Reliability of the Tsunakawa‐Shaw and Other Methods With Implications for “Fragile” Curvature

AbstractAbsolute paleointensity (API) of the geomagnetic field can be estimated from volcanic rocks by comparing the natural remanent magnetization (NRM) to a laboratory‐induced thermoremanent magnetization (Lab‐TRM). Plots of NRM unblocking versus Lab‐TRM blocking from API experiments often exhibit nonideal curvature, which can result in biased estimates. Previous work showed that curvature can increase with age; however, selection criteria designed to eliminate such behavior yielded accurate estimates for two‐year‐aged specimens (70.3 ± 3.8 μT;N = 96 specimens out of 120 experiments). API can also be estimated in coercivity space. Here, we use the Tsunakawa‐Shaw (TS) method applied to 20 specimens aged in the laboratory field of 70.0 μT for 4 years, after acquisition of zero‐age (fresh) Lab‐TRM in the same field. Selection criteria for the TS experiment also yielded accurate results (68.5 ± 4.5 μT;N = 17 specimens). In thermal API experiments, curvature is related to internal structure with more single domain‐like behavior having the least curvature. Here we show that the fraction of anhysteretic remanent magnetization demagnetized by low‐temperature treatment was larger for samples with larger thermal curvatures suggesting a magnetocrystalline anisotropy source. We also tested experimental remedies that have been proposed to improve the accuracy of paleointensity estimates. In particular, we test the efficacy of the multi‐specimen approach and a strategy pretreating specimens with low field alternating field demagnetization prior to the paleointensity experiment. Neither yielded accurate results.

Investigating the utility of a high-temperature Thellier-style paleointensity experimental protocol

High-quality data are vital to the research field of paleointensity, which has long suffered from poor-quality and/or sparse data. Previous paleointensity work has established that repeatedly heating specimens increases the opportunity for thermochemical alteration to decrease the reliability of paleointensity data. In addition, recent work has shown that repeatedly heating specimens in paleointensity experiments can also exaggerate the effects of non-ideal, non-single domain grains. Arai plots resulting from paleointensity experiments containing such grains are often curvilinear (two-slope) across most of the specimen’s unblocking temperature spectrum, except in the temperature range nearest to the grains’ Curie temperature. This study tests the following strategy to mitigate these effects: that of performing a Thellier paleointensity experiment using fewer temperature steps that are also concentrated at higher temperatures. For this purpose, we use natural specimens with well-constrained rock magnetic data from the Hawaiian Scientific Observation Hole 1 (SOH1) drill core in paleointensity experiments with starting temperatures ranging from 200 °C to 500+ °C. Those experiments that focused in on the portion of the unblocking temperature spectrum nearest the Curie temperature of the specimen (HiTeCT) had an exceptionally low success rate, whereas those with initial temperatures at comparatively moderate temperatures (200–400 °C) had high success rates (~ 70%). Thermochemical alteration was minimized with a start temperature of 400 °C, but the curvature of the Arai plots had no clear dependance on start temperature. We conclude herein that increasing the start temperature can help avoid the effects of low-temperature alterations. Additionally, we found that the approach of focusing in on the highest temperature range is not a feasible one to apply in paleointensity experiments, in general, and consider this likely to be a result of a lack of intermediate-temperature checks for alteration and insufficient independence of temperature steps.Graphical abstract

New Paleointensity Data From Aniakchak Volcano, Alaska, USA

AbstractThis study presents the first set of Holocene, high‐quality absolute paleointensity data from Alaska, USA. Existing paleointensity data for the Holocene are generally located at mid‐northern latitudes in North America, Europe, the Middle East, and eastern Asia. Relatively few data are from the Alaska region. IZZI‐modified paleointensity experiments were conducted on glassy volcanic materials from Aniakchak volcano, a mid‐ to high‐latitude composite volcano on the Alaska‐Aleutian arc. The CCRIT selection criteria were applied to the paleointensity results. A total of 30 specimens from six samples with estimated ages ranging from 1931 CE to 2,300 years before present passed all selection criteria. The sample‐mean paleointensities ranged from 49.5 to 68.0 microTesla (μT). The sample‐mean paleointensity results are comparable to modeled intensities, however all except for one sample‐mean paleointensity are lower than those predicted by geomagnetic field models. The paleointensity estimate for the historical 1931 CE eruption was about 15 μT greater than the expected field strength. This overestimate may result from unrecognized alteration or non‐ideal remanence carriers in this sample. Further evaluation of samples from the 1931 CE eruption using a Bayesian estimation method resulted in a paleointensity estimate that encompasses the historical field strength within uncertainty. These new paleointensity results are a valuable contribution to the mid‐ to high‐northern latitude paleomagnetic data set. Incorporation of these data into future geomagnetic field models will improve the predictions of geomagnetic field behavior in the Alaska region.

Bias Corrected Estimation of Paleointensity (BiCEP): An Improved Methodology for Obtaining Paleointensity Estimates

AbstractThe assumptions of paleointensity experiments are violated in many natural and archeological materials, leading to Arai plots which do not appear linear and yield inaccurate paleointensity estimates, leading to bias in the result. Recently, paleomagnetists have adopted sets of “selection criteria” that exclude specimens with nonlinear Arai plots from the analysis, but there is little consensus in the paleomagnetic community on which set to use. In this study, we present a statistical method we call Bias Corrected Estimation of Paleointensity (BiCEP), which assumes that the paleointensity recorded by each specimen is biased away from a true answer by an amount that is dependent a single metric of nonlinearity (the curvature parameter) on the Arai plot. We can use this empirical relationship to estimate the recorded paleointensity for a specimen where, that is, a perfectly straight line. We apply the BiCEP method to a collection of 30 sites for which the true value of the original field is well constrained. Our method returns accurate estimates of paleointensity, with similar levels of accuracy and precision to restrictive sets of paleointensity criteria, but accepting as many sites as permissive criteria. The BiCEP method has a significant advantage over using these selection criteria because it achieves these accurate results without excluding large numbers of specimens from the analysis. It yields accurate, albeit imprecise estimates from sites whose specimens all fail traditional criteria. BiCEP combines the accuracy of the strictest selection criteria with the low failure rates of the less reliable “loose” criteria.

New criteria for selecting reliable Thellier-type paleointensity results from the 1960 Kilauea lava flows, Hawaii

Thellier-type paleointensity experiments associated with partial thermal remanent magnetization checks have been widely used to determine paleointensity values from volcanic and archaeological media. However, previous studies have revealed that a substantial portion of paleointensity results with positive checks for historical lava samples largely fails to predict known Earth magnetic field intensity values. To determine the fidelity of paleointensity values, conventional Thellier-type paleointensity experiments were performed on Kilauea lava flows that erupted in 1960. The positive partial thermal remanent magnetization checks for our results range from 30.28 ± 1.38 µT to 52.94 ± 1.89 µT. This strongly indicates that conventional paleointensity checks cannot guarantee the fidelity of paleointensity results, especially when the unblocking temperatures for the newly formed magnetic particles are higher than the treatment temperature. Therefore, in this study, to check for thermal alteration during heating, the temperature dependence of the hysteresis parameter measured at room temperature for the thermally treated samples was also measured. Our new results show that nearly all biased paleointensity values correspond to a ratio of the coercivity of remanence to the magnetic coercivity of > 3 and a chemical alteration index > ~ 10%, which indicates the strong effect of the domain state and thermal alteration on the fidelity of the paleointensity results. Our study provides feasible criteria to further improve the fidelity of paleointensity estimations.

An Integrated Paleomagnetic, Multimethod‐Paleointensity, and Radiometric Study on Cretaceous and Paleogene Lavas From the Lesser Caucasus: Geomagnetic and Tectonic Implications

AbstractSixteen rhyolitic and dacitic Cretaceous and Paleocene‐Eocene lavas from the Lesser Caucasus have been subjected to paleomagnetic and multimethod paleointensity experiments to analyze the variations of the Earth's magnetic field. Paleointensity experiments were performed with two methods. Thellier‐type experiments with the IZZI method on 65 specimens (nine flows) yielded 15 successful determinations and experiments with the multispecimen method on 14 samples (seven flows) yielded two successful determinations. The joint analysis of the results obtained with both methods produced a mean FuK = (19.9 ± 3.7) µT for upper Cretaceous and FPg = (20.7 ± 3.3) µT for Paleogene sites. Low virtual axial dipole moments for the Cretaceous (3.4 × 1022 Am2) and Paleogene (3.5 × 1022 Am2) samples support the idea of a lower average dipole moment during periods of stable polarity of the Earth magnetic field. Mean flow paleomagnetic directions did not match expected upper Cretaceous to Paleogene directions calculated from the European Apparent Polar Wander Path. While inclination results roughly agreed with expected values, a group of sites showed nearly North‐South paleodeclinations (D = 1.1° ± 14.2°), and another group displayed eastward deviated paleodeclinations (D = 72.9° ± 26.6°). These results suggest the occurrence of nearly vertical‐axis rotations, probably as a result of continental collision since Oligocene. In addition to paleomagnetic and palaeointensity analyses, new K‐Ar absolute age determinations have been performed on three of the studied sites, yielding Late Cretaceous ages (78.7 ± 1.7, 79.7 ± 1.6, and 83.4 ± 1.8 Ma (2σ)).

Understanding Nonideal Paleointensity Recording in Igneous Rocks: Insights From Aging Experiments on Lava Samples and the Causes and Consequences of “Fragile” Curvature in Arai Plots

AbstractThe theory for recording of thermally blocked remanences predicts a quasilinear relationship between low fields like the Earth's in which rocks cool and acquire a magnetization. This serves as the foundation for estimating ancient magnetic field strengths. Addressing long‐standing questions concerning Earth's magnetic field requires a global paleointensity data set, but recovering the ancient field strength is complicated because the theory only pertains to uniformly magnetized particles. A key requirement of a paleointensity experiment is that a magnetization blocked at a given temperature should be unblocked by zero‐field reheating to the same temperature. However, failure of this requirement occurs frequently and the causes and consequences of failure are understood incompletely. Recent experiments demonstrate that the remanence in many samples typical of those used in paleointensity experiments is unstable, exhibiting an “aging” effect in which the (un)blocking temperature spectra can change over only a few years resulting in nonideal experimental behavior. While a fresh remanence may conform to the requirement of equality of blocking and unblocking temperatures, aged remanences may not. Blocking temperature spectra can be unstable (fragile), which precludes reproduction of the conditions under which the original magnetization was acquired. This limits our ability to acquire accurate and precise ancient magnetic field strength estimates because differences between known and estimated fields can be significant for individual specimens, with a low field bias. Fragility of unblocking temperature spectra may be related to grain sizes with lower energy barriers and may be detected by features observed in first‐order reversal curves.

First archeomagnetic data from Kenya and Chad: Analysis of iron furnaces from Mount Kenya and Guéra Massif

The variation of the Earth's magnetic field over the last millennia is poorly known in Africa, especially in Sub-Saharan Africa that represents less than 1% of the global archeomagnetic dataset. Fourteen iron furnaces from Kenya and Chad have been studied here for archeomagnetic purposes. These structures were dated by 14C from the 14th century in Kenya and from the 18-19th centuries in Chad. Samples, oriented with the plaster cap technique, have been subjected to thermal and alternating field demagnetization, rock magnetic analysis and paleointensity experiment, with the classical Thellier-Thellier technique corrected for anisotropy and cooling rate effects. New directional data (4 from Kenya and 10 from Chad) and one new paleointensity estimate from Kenya (43.5 ± 1.8 μT) were obtained. The Kenyan declinations are in good agreement with the geomagnetic global models, but the inclinations are lower and the intensity higher than predicted. For Chad, the new results prove the efficiency of the archeomagnetic dating, with a precision that is better than radiocarbon for the 18–19th centuries.

Paleomagnetism of Paleocene‐Maastrichtian (60–70 Ma) Lava Flows From Tian Shan (Central Asia): Directional Analysis and Paleointensities

AbstractThe Tuoyun volcanics from the Tian Shan range (Central Asia) give an opportunity to investigate the variability of the Earth's magnetic field during the Cretaceous and Early Paleogene. In the paper we focus on Maastrichtian‐Paleocene basalts and report new paleomagnetic results from 70 lava flows (respectively 45 directional groups) sampled near Tuoyun village (75.33°E; 40.18°N) from three distinct sections. Combined with previous results, our new data set yields a virtual geomagnetic pole at φS = 180.2°E, ϑS = 49.5°N with an angular standard deviation (for N = 93 directional groups). The mean inclination Is = 36.2° is indicative of an emplacement of the Tuoyun volcanics at a paleolatitude around 20°N, consistent with the high dispersion k = 14.6 of the directions. Such an inclination value is 10–15° lower than the predictions from the apparent polar wander paths for stable Europe and East Asia. This observation can be partly explained by crustal shortening within the East Asian plate during the India‐Asia convergence, but also needs to invoke a local field anomaly in Central Asia as previously proposed. Our absolute paleointensity experiments conducted on two lava flows and two baked sedimentary layers yield a mean virtual dipole moment of 58.9 ± 6.9 ZA·m2, consistent with the values at 60–70 Ma in the global paleointensity database. Compared to relative paleointensity results during the Cretaceous Normal Superchron and during the Plio‐Pleistocene, the relative variability in intensity—supposed to be a proxy for the geodynamo's activity—estimated here at 31 ± 5% (N = 15) from pseudo‐Thellier experiments possibly corroborates a correlation between geomagnetic reversal frequency and paleosecular variation rate during the past 120 Myr.

Comparison of Thermal and Microwave Paleointensity Estimates in Specimens Displaying Non‐Ideal Behavior in Thellier‐Style Paleointensity Experiments

AbstractDetermining the strength of the ancient geomagnetic field is vital to our understanding of the core and geodynamo, but obtaining reliable measurements of the paleointensity is fraught with difficulties. Over a quarter of magnetic field strength estimates within the global paleointensity database from 0 to 5 Ma come from Hawaii. Two previous studies on the SOH1 drill core gave inconsistent, apparently method‐dependent paleointensity estimates, with an average difference of 30%. The paleointensity methods employed in the two studies differed both in demagnetization mechanism (thermal or microwave radiation) and Thellier‐style protocol (perpendicular and original Thellier protocols)—both variables that could cause the strong differences in the estimates obtained. Paleointensity experiments have therefore been conducted on 79 specimens using the previously untested combinations of thermal‐perpendicular and microwave‐original Thellier methods to analyze the effects of demagnetization mechanism and protocol in isolation. We find that, individually, neither demagnetization mechanism nor protocol entirely explains the differences in paleointensity estimates. Specifically, we found that non‐ideal multidomain‐like effects are enhanced using the original Thellier protocol (independent of demagnetization mechanism), often resulting in paleointensity overestimation. However, we also find evidence, supporting recent findings from the 1960 Kilauea lava flow, that microwave‐perpendicular experiments performed without partial thermal remanent magnetization checks can produce underestimates of the paleointensity due to unaccounted‐for sample alteration at higher microwave powers. Together, these findings support that the true paleointensities fall between the estimates previously published and emphasize the need for future studies (thermal or microwave) to use protocols with both partial thermal remanent magnetization checks and a means of detecting non‐ideal grain effects.

Archeointensity of the Four Corners Region of the American Southwest

AbstractNew paleointensity experiments were conducted using the IZZI protocol (a variation on the Thellier‐Thellier method) on 289 specimens from 54 baked pottery fragments collected from four archeological sites in the American Southwest. Anisotropy experiments were conducted to correct for anisotropy of remanence, a common problem in archeological material. Additionally, the effect of cooling rate was evaluated and the calculated paleointensity values were adjusted accordingly. Using the Thellier graphical user interface program, the specimen results were analyzed, averaged by sample (i.e., pottery fragment), filtered for the highest quality, and converted to Virtual Axial Dipole Moments (VADMs). Stylistic evidence, historical documentation, and dendrochronology analyses provide age constraints with up to decade resolution for the VADM results. We compared these new results with the highest quality previously published paleointensity values from the Four Corners region of the American Southwest—defined here as the four states of New Mexico, Arizona, Utah, and Colorado. None of the previously published data were corrected for cooling rate and only some were corrected for anisotropy, resulting in a systematic bias between data sets. To accommodate this difference, an estimated cooling rate correction was applied to all the previously published data. No correction can be made for anisotropy as this is specimen specific. Our estimated cooling rate correction and the empirical correction applied to the new data both require an estimation of the historical cooling rate of the pottery. Experimental pottery firings using one ancient technique and outfitted with thermocouples were conducted to determine the historical cooling rate.

Absolute Paleointensity Study of Miocene Tiva Canyon Tuff, Yucca Mountain, Nevada: Role of Fine‐Particle Grain‐Size Variations

AbstractFine‐grained, Ti‐poor titanomagnetite in the ~12.7 Ma Tiva Canyon (TC) Tuff systematically increases in grain size from superparamagnetic (SP) at the flow base to single domain (SD) at a few meters height. This allows us to examine the role of grain‐size variation on paleointensity, within the transition from SP to stable SD. We present magnetic properties from two previously unreported sections of the TC Tuff, as well as Thellier‐type paleointensity estimates from the lowermost ~7.0 m of the flow. Magnetic hysteresis, frequency‐dependent susceptibility, and thermomagnetic data show that sample grain‐size distribution is dominated by SP in the lower ~3.6 m, transitioning upwards to mostly stable SD. Paleointensity results are closely tied to stratigraphic height and to magnetic properties linked to domain state. SD samples have consistent absolute paleointensity values of 28.5 ± 1.94 μT (VADM of 51.3 ZAm2) and behaved ideally during paleointensity experiments. The samples including a significant SP fraction have consistently higher paleointensities and less ideal behavior but would likely pass many traditional quality‐control tests. We interpret the SD remanence to be a primary thermal remanent magnetization but discuss the possibility of a partial thermal‐chemical remanent magnetization if microcrystal growth continued at T < Tc and/or the section is affected by post‐emplacement vapor‐phase alteration. The link between paleointensity and domain state is stronger than correlations with water content or other evidence of alteration and suggests that the presence of a significant SP population may adversely impact paleointensity results, even in the presence of a stable SD fraction.

Experimental Evaluation of Remanence Carriers Using the Microcoercivity‐Unblocking Temperature Diagram

AbstractKnowledge of remanence carriers in recording media is crucial for obtaining reliable paleomagnetic data. In this paper, we report a newly developed technique to determine remanence carriers in natural rock samples that uses anhysteretic remanent magnetization (ARM) spectra in the microcoercivity (hc) and unblocking temperature (TUB) planes. The distribution of ARM carriers was determined via systematic stepwise demagnetization experiments of ARM combining thermal demagnetization and alternating field demagnetization treatments and calculating the change in ARM intensity on the hc‐TUB plane (hc‐TUB diagram). Samples from three basaltic lavas and two welded tuffs were analyzed. Combining the hc‐TUB diagram with the change in ARM coercivity spectra, both decreases and increases in magnetic minerals due to the thermal alteration during laboratory heating could be identified. From experimental samples, ARM carrier components were identified as being original rather than produced as a result of thermal alteration. In the case of the Izu‐Oshima 1986 A lava flow (1986La) sample, it was confirmed from the hc‐TUB diagram that high‐hc components can be separated using stepwise alternating field demagnetization treatment but not stepwise thermal demagnetization treatment. Paleointensity experiments were also conducted on Izu‐Oshima 1986 A lava flow samples. Two of the three paleointensity estimates were greater than the expected field intensity, suggesting a certain rock‐magnetic property that affects the paleointensity estimation. The hc‐TUB diagrams of the Izu‐Oshima 1986 A lava flow samples suggest that the differences in the paleointensity estimates were likely caused by different degrees of the high‐temperature oxidation state.

Investigating the Accuracy, Precision, and Cooling Rate Dependence of Laboratory‐Acquired Thermal Remanences During Paleointensity Experiments

AbstractWe examine the behavior of natural basaltic and trachytic samples during paleointensity experiments on both the original and laboratory‐acquired thermal remanences and characterize the samples using proxies for domain state including curvature (k) and the bulk domain stability parameters of Paterson (2011,https://doi.org/10.1029/2011JB008369) and Paterson et al. (2017,https://doi.org/10.1073/pnas.1714047114), respectively. A curvature value of 0.164 (suggested by Paterson, 2011,https://doi.org/10.1029/2011JB008369) as a critical threshold that separates single‐domain‐like remanences from multidomain‐like remanances on the original paleointensity data was used to separate samples into “straight” (single‐domain‐like) and “curved” (multidomain‐like) groups. Specimens from the two sample sets were given a “fresh” thermal remanent magnetization in a 70 μT field and subjected to an infield‐zerofield, zerofield‐infield (IZZI)‐type (Yu et al., 2004,https://doi.org/10.1029/2003GC000630) paleointensity experiment. The straight sample set recovered the laboratory field with high precision while the curved set had much more scattered results (70.5 ± 1.5 and 71.9 ± 5.2 μT, respectively). The average intensity of both sets for straight and curved was quite close to the laboratory field of 70 μT, however, suggesting that if experiments contain a sufficient number of specimens, there does not seem to be a large bias in the field estimate. We found that the dependence of the laboratory thermal remanent magnetization on cooling rate was significant in most samples and did not depend on domain states inferred from proxies based on hysteresis measurements and should be estimated for all samples whose cooling rates differ from that used in the laboratory.