High Precision Sr‐Nd‐Hf‐Pb Isotopic Compositions of USGS Reference Material BCR‐2

The Lamont‐Doherty Earth Observatory radiogenic isotope group has been systematically measuring Sr‐Nd‐Pb‐Hf isotopes of USGS reference material BCR‐2 (Columbia River Basalt 2), as a chemical processing and instrumental quality control monitor for isotopic measurements. BCR‐2 is now a widely used geochemical inter‐laboratory reference material (RM), with its predecessor BCR‐1 no longer available. Recognising that precise and accurate data on RMs is important for ensuring analytical quality and for comparing data between different laboratories, we present a compilation of multiple digestions and analyses made on BCR‐2 during the first author's dissertation research. The best estimates of Sr, Nd and Hf isotope ratios and measurement reproducibilities, after filtering at the 2s level for outliers, were ⁸⁷Sr/⁸⁶Sr = 0.705000 ± 11 (2s, 16 ppm, n = 21, sixteen digestions, one outlier), ¹⁴³Nd/¹⁴⁴Nd = 0.512637 ± 13 (2s, 25 ppm, n = 27, thirteen digestions, one outlier) and ¹⁷⁶Hf/¹⁷⁷Hf = 0.282866 ± 11 (2s, 39 ppm, n = 25, thirteen digestions, no outliers). Mean Nd and Hf values were within error of those reported by Weis et al. (2006, 2007) in their studies of RMs; mean Sr values were just outside the 2s uncertainty range of both laboratories. Moreover, a survey of published Sr‐Nd‐Hf data shows that our results fall within the range of reported values, but with a smaller variability. Our Pb isotope results on acid leached BCR‐2 aliquots (n = 26, twelve digestions, two outliers) were ²⁰⁶Pb/²⁰⁴Pb = 18.8029 ± 10 (2s, 55 ppm), ²⁰⁷Pb/²⁰⁴Pb = 15.6239 ± 8 (2s, 52 ppm), ²⁰⁸Pb/²⁰⁴Pb = 38.8287 ± 25 (2s, 63 ppm). We confirm that unleached BCR‐2 powder is contaminated with Pb, and that sufficient leaching prior to digestion is required to achieve accurate values for the uncontaminated Pb isotopic compositions.     

Measurement of 34S/32S Ratios of NBS 120c and BCR 32 Phosphorites Using Purge and Trap EA‐IRMS Technology

Measurements of sulfur stable isotope ratios (³⁴S/³²S) have suffered from technical difficulties in analysing low‐S materials reducing their use despite their undeniable scientific interest. The measurement of ³⁴S/³²S ratios is a powerful tool for deciphering problems such as determining the sources of environmental pollutants, to detect adulteration, tracking the evolution of the redox state of the oceans and quantifying the role of the bacterial activity in sulfide minerals genesis. We have used a high‐precision method of sulfur isotope determination using a new type of elemental analyser based on ‘purge and trap’ technology. This new technique demonstrates the high quality of ³⁴S/³²S measurements for samples with S concentrations lower than 1% m/m. International calibrated references of diverse sulfur‐bearing materials were used to calibrate two low (< 1%) S‐bearing phosphorites used as compositional reference material for future use as isotopic references: BCR 32 and NBS 120c. δ³⁴S_CDT values of, respectively, 18.2‰ (1s = 0.3; n = 23) and 18.3‰ (1s = 0.4; n = 20) are proposed for these. Calibration of both phosphorites with international reference materials led to calculation of a mean standard error close to 0.4‰. The demonstration of a capability to reliably measure S isotope ratios in low‐S phosphate minerals or rocks opens up new fields of palaeoenvironmental reconstructions.     

Matrix Corrections and Error Analysis in High‐Precision SIMS 18O/16O Measurements of Ca–Mg–Fe Garnet

We report technical and data treatment methods for making accurate, high‐precision measurements of ¹⁸O/¹⁶O in Ca–Mg–Fe garnet utilising the Cameca IMS 1280 multi‐collector ion microprobe. Matrix effects were similar to those shown by previous work, whereby Ca abundance is correlated with instrumental mass fractionation (IMF). After correction for this effect, there appeared to be no significant secondary effect associated with Mg/Fe²⁺ for routine operational conditions. In contrast, investigation of the IMF associated with Mn‐ or Cr‐rich garnet showed that these substitutions are significant and require a more complex calibration scheme. The Ca‐related calibration applied to low‐Cr, low‐Mn garnet was reproducible across different sample mounts and under a range of instrument settings and therefore should be applicable to similar instruments of this type. The repeatability of the measurements was often better than ± 0.2‰ (2s), a precision that is similar to the repeatability of bulk techniques. At this precision, the uncertainties due to spot‐to‐spot repeatability were at the same magnitude as those associated with matrix corrections (± 0.1–0.3‰) and the uncertainties in reference materials (± 0.1–0.2‰). Therefore, it is necessary to accurately estimate and propagate uncertainties associated with these parameters – in some cases, uncertainties in reference materials or matrix corrections dominate the uncertainty budget.     

Eclogite and pyroxenite xenoliths from off-craton kimberlites near the Kaapvaal Craton, South Africa

Mantle xenoliths brought to the surface by kimberlite magmas along the south-western margin of the Kaapvaal craton in South Africa can be subdivided into eclogites sensu stricto, kyanite eclogites and orthopyroxene eclogites, all containing omphacite, and garnet clinopyroxenites and garnet websterites characterised by diopside. Texturally, chemically (major elements) and thermally, we observe an evolution from garnet websterites (T_EG = 742–781 °C) towards garnet clinopyroxenites (T_EG = 715–830 °C) and to eclogites (T_EG = 707–1056 °C, mean value of 913 °C). Pressures calculated for orthopyroxene-bearing samples suggest upper mantle conditions of equilibration (P = 16–33 kb for the garnet websterites, 18 kb for a garnet clinopyroxenite and 23 kb for an opx-bearing eclogite). The overall geochemical similarity between the two groups of xenoliths (omphacite-bearing and diopside-bearing) as well as the similar trace element patterns of clinopyroxenes and garnet suggest a common origin for these rocks. Recently acquired oxygen isotope data on garnet (δ¹⁸O_gnt = 5.25–6.78 ‰ for eclogites, δ¹⁸O_gnt = 5.24–7.03 ‰ for garnet clinopyroxenites) yield values ranging from typical mantle values to other interpreted as resulting from low-temperature alteration or precursors sea-floor basalts and associated rocks. These rocks could then represent former magmatic oceanic rocks that crystallised from a same parental magma as plagioclase free diopside-bearing and plagioclase-bearing crustal rocks. During subduction, these oceanic rock protoliths equilibrated at mantle depth, with the plagioclase-bearing rocks converting to omphacite and garnet-bearing lithologies (eclogites sensu largo), whereas the plagioclase-free diopside-bearing rocks converted to diopside and garnet-bearing lithologies (garnet websterites and garnet clinopyroxenites).     

Holocene relative sea-level changes and vertical movements along the Italian and Istrian coastlines

Published and new data exist for relative sea-level change for 105 locations (127 samples) during the late Holocene, along the Italian (and Istrian) coasts. These data, compared with predictions (derived from two different models associated with the last glacial cycle) allowed the calculation of the tectonic vertical movements. They are based on precise measures of geomorphological and archaeological markers between 0.4 and 12.6 ka cal. BP, sampled at elevations between +7 and −51 m. In order to decipher the broad pattern of Holocene tectonic vertical movements along the Italian coastline, these data were compared with predicted sea-level curves using the most recent models published for the Mediterranean sea. Tectonic rates varied from −4.85 mm/a, in a core at Sybaris, to 5 mm/a, in the volcanic areas of Pozzuoli and Pantelleria. New MIS 5.5 (125 ka) data, mostly from the Venetian plain, are reported. In particular the depth of the base of MIS 5.5 paralic deposits found in four cores near Venezia provides a mean subsidence of 0.62 mm/a. New, precise mass spectrometer U-Th analyses on Cladocora layers from the bottom of a long core (named ENEA), indicate older ages (195.7 ± 1.6 and 161.2 ± 1.2 ka, respectively), relative to the published MIS 5.5 ages, which were based on alpha-counting U-Th data.Instrumental data obtained from tide gauges and repeated levelling measurements from the NE Adriatic and Sicily are also considered. These methods have one great advantage with respect to continuous GPS measurements and the satellite altimetric observations, in that a much greater time span is available. Although the altimetric measurements are available for 16 years, and the GPS for less than a decade, repeated levelling lines cover up to 50 years and tide gauge observations in some cases to 100 years or more. The greater time span allows for more stable differential rate estimates. The repeated levelling shows that the plain east of Mestre is subsiding (to −4 mm/a). The Messina tidal gauge demonstrates a total coseismic and post-seismic subsidence of 77 cm associated with the event of 1908, the post-seismic phase lasting for at least 13 years. The Reggio Calabria tidal station points to an uplift of this station relative to Palermo in the order of 1–2 mm/a.     

Discrimination between primary magmatic biotites, reequilibrated biotites and neoformed biotites

The ternary diagram TiO₂–FeO*–MgO (FeO* = FeO + MnO) is proposed as a quantitative objective tool for distinguishing between primary magmatic biotites and those that are more or less reequilibrated, or possibly neoformed, by or within a hydrothermal fluid. The limit of the domains of the primary magmatic biotites, the reequilibrated biotites and the neoformed biotites were determined on the basis of optical, paragenetic and chemical criteria. To cite this article: H. Nachit et al., C. R. Geoscience 337 (2005).     

Non‐Matrix‐Matched Calibration for the Multi‐Element Analysis of Geological and Environmental Samples Using 200 nm Femtosecond LA‐ICP‐MS: A Comparison with Nanosecond Lasers

LA‐ICP‐MS is one of the most promising techniques for in situ analysis of geological and environmental samples. However, there are some limitations with respect to measurement accuracy, in particular for volatile and siderophile/chalcophile elements, when using non‐matrix‐matched calibration. We therefore investigated matrix‐related effects with a new 200 nm femtosecond (fs) laser ablation system (NWRFemto200) using reference materials with different matrices and spot sizes from 10 to 55 μm. We also performed similar experiments with two nanosecond (ns) lasers, a 193 nm excimer (ESI NWR 193) and a 213 nm Nd:YAG (NWR UP‐213) laser. The ion intensity of the 200 nm fs laser ablation was much lower than that of the 213 nm Nd:YAG laser, because the ablation rate was a factor of about 30 lower. Our experiments did not show significant matrix dependency with the 200 nm fs laser. Therefore, a non‐matrix‐matched calibration for the multi‐element analysis of quite different matrices could be performed. This is demonstrated with analytical results from twenty‐two international synthetic silicate glass, geological glass, mineral, phosphate and carbonate reference materials. Calibration was performed with the certified NIST SRM 610 glass, exclusively. Within overall analytical uncertainties, the 200 nm fs LA‐ICP‐MS data agreed with available reference values.     

Zinc Isotopic Compositions of NIST SRM 683 and Whole‐Rock Reference Materials

In this study the homogeneity of the zinc isotopic composition in the NIST SRM 683 reference material was examined by measuring the Zn isotopic signature in microdrilled sample powders from two metal nuggets. Zinc was purified using AG MP‐1M resin and then measured by MC‐ICP‐MS. Instrumental mass bias was corrected using the “sample‐standard bracketing” method and empirical external normalisation with Cu doping. After evaluating the potential effects of varying acid mass fractions and different matrices, high‐precision Zn isotope data were obtained with an intermediate measurement precision better than ± 0.05‰ (δ⁶⁶Zn, 2s) over a period of 5 months. The δ⁶⁶Zn_JMC‐Lyon mean values of eighty‐four and fourteen drilled powders from two nuggets were 0.11 ± 0.02‰ and 0.12 ± 0.02‰, respectively, indicating that NIST SRM 683 is a good isotopic reference material with homogeneous Zn isotopes. The Zn isotopic compositions of seventeen rock reference materials were also determined, and their δ⁶⁶Zn values were in agreement with most previously published data within 2s. The δ⁶⁶Zn values of most of the rock reference materials analysed were in the range 0.22–0.36‰, except for GSP‐2 (1.07 ± 0.06‰, n = 12), NOD‐A‐1 (0.96 ± 0.03‰, n = 6) and NOD‐P‐1 (0.78 ± 0.03‰, n = 6). These comprehensive data should serve as reference values for quality assurance and interlaboratory calibration exercises.     

High‐Precision Measurement of Molybdenum Isotopic Compositions of Selected Geochemical Reference Materials

Here we describe high‐precision molybdenum isotopic composition measurements of geological reference materials, performed using multi‐collector inductively coupled plasma‐mass spectrometry (MC‐ICP‐MS). Purification of Mo for isotopic measurements was achieved by ion exchange chromatography using Bio‐Rad AG^® 1‐X8 anion exchange resin. Instrumental mass bias was corrected using ¹⁰⁰Mo‐⁹⁷Mo double spiking techniques. The precision under intermediate measurement conditions (eighteen measurement sessions over 20 months) in terms of δ^98/95Mo was 0.10‰ (2s). The measurement output was approximately four times more efficient than previous techniques, with no compromise in precision. The Mo isotopic compositions of seven geochemical reference materials, seawater (IAPSO), manganese nodules (NOD‐P‐1 and NOD‐A‐1), copper‐molybdenum ore (HV‐2), basalt (BCR‐2) and shale (SGR‐1b and SCo‐1), were measured. δ^98/95Mo values were obtained for IAPSO (2.25 ± 0.09‰), NOD‐P‐1 (?0.66 ± 0.05‰), NOD‐A‐1 (?0.48 ± 0.05‰), HV‐2 (?0.23 ± 0.10‰), BCR‐2 (0.21 ± 0.07‰), SCo‐1 (?0.24 ± 0.06‰) and SGR‐1b (0.63 ± 0.02‰) by calculating δ^98/95Mo relative to NIST SRM 3134 (0.25‰, 2s). The molybdenum isotopic compositions of IAPSO, NOD‐A‐1 and NOD‐P‐1 obtained in this study are within error of the compositions reported previously. Molybdenum isotopic compositions for BCR‐2, SCo‐1 and SGR‐1b are reported for the first time.     

Magnesium Isotopic Compositions of International Geological Reference Materials

Magnesium isotopic compositions are reported for twenty‐four international geological reference materials including igneous, metamorphic and sedimentary rocks, as well as phlogopite and serpentine minerals. The long‐term reproducibility of Mg isotopic determination, based on 4‐year analyses of olivine and seawater samples, was ≤ 0.07‰ (2s) for δ²⁶Mg and ≤ 0.05‰ (2s) for δ²⁵Mg. Accuracy was tested by analysis of synthetic reference materials down to the quoted long‐term reproducibility. This comprehensive dataset, plus seawater data produced in the same laboratory, serves as a reference for quality assurance and inter‐laboratory comparison of high‐precision Mg isotopic data.     

Quasi-static slip on the plate boundary associated with the 2003 M8.0 Tokachi-oki and 2004 M7.1 off-Kushiro earthquakes, Japan

We analyzed small repeating earthquakes recorded over a 13-year period and GPS data recorded over an 8-month period to estimate interplate quasi-static slip associated with the 2003 Tokachi-oki earthquake (M8.0) and the 2004 off-Kushiro earthquake (M7.1). The repeating-earthquake analysis revealed that the slip rate near the source region of the Tokachi-oki earthquake was relatively low (< 5 cm/year) prior to the earthquake; however, in the last 3 years leading up to the event, a minor acceleration in slip occurred upon the deeper extension of the coseismic slip area of the earthquake. Repeating-earthquake and GPS data indicate that large amounts of afterslip occurred around the rupture area following the earthquake; the afterslip mainly propagated to the east of the coseismic slip area. We also infer that the occurrence of the 2004 off-Kushiro earthquake, located about 100 km northeast of the epicenter of the Tokachi-oki earthquake, was advanced by the afterslip associated with the Tokachi-oki earthquake.     

Dating Micrometre‐Thin Rims Using a LA‐ICP‐MS Depth Profiling Technique on Zircon from an Archaean Metasediment: Comparison with the SIMS Depth Profiling Method

Laser ablation‐inductively coupled plasma‐mass spectrometry (LA‐ICP‐MS) was examined as a tool for measuring isotopic variation as a function of ablation depth in unpolished zircon from an Archaean metasediment specimen. This technique was able to identify micrometre‐thin (> 3 μm) isotopically distinct mineral domains characterised by ca. 100 Myr younger ²⁰⁷Pb/²⁰⁶Pb ages associated with 2s age uncertainties as low ~ 0.2%, as well as elevated U content relative to grain interiors (up to an order of magnitude). Our calculated drilling rate suggests that each laser pulse excavated depths of ~ 0.06 μm. Ages resolved through the LA‐ICP‐MS methods overlap the 2s uncertainties of ²⁰⁷Pb/²⁰⁶Pb ages measured using SIMS depth profiling on the same zircon population. The rims were further evinced by the detection of relative enrichment (> 3 orders of magnitude) in REE in the outermost micrometres of the same zircon, measured using a different and independent LA‐ICP‐MS depth profiling technique. We propose a LA‐ICP‐MS U–Pb technique capable of quickly identifying and quantifying rims, which are indication of late, yet geologically significant, fluid events that are otherwise undefined.     

A spectrophotometric measurement of soil cation exchange capacity based on cobaltihexamine chloride absorbance

Among numerous methods for cation exchange capacity (CEC) determination for soils and sediments, the cobaltihexamine chloride method is frequently used due to its ability to measure CEC at soil pH. After exchange with Co(NH₃)₆³⁺ ions, CEC is estimated via the measurement of the Co remaining in solution. The modified method proposed allows a more rapid determination of CEC based on the measurement of the absorbance at 472 nm of the cobaltihexamine chloride solution before and after exchange. This method has been applied to various soil's horizons from four sites, selected to cover a wide range of CEC and pH values. The model obtained allows one to calculate CEC from absorbance at 472 nm with 95% confidence intervals. As CEC is of relevant meaning in agronomical and environmental purposes, and more recently in ecotoxicological studies, this modified method can be proposed as a rapid test for CEC evaluation.     

The Alleret maar (Massif Central, France): A new lacustrine sequence of the early Middle Pleistocene in western Europe

The Alleret maar (Massif Central, France) provides a long lacustrine sequence (40.6 m) attributed to the early Middle Pleistocene. Sediment, pollen and diatoms analysis of its upper part (AL2 core, 14.6 m) indicates two temperate phases marked by high lake levels, forest development and vegetation expansion. They are separated by a cold period during which lake level drops, coarse sediment input increases and steppic and xerophilous plants develop. Pollen data suggests that this sequence belongs to the upper part of the Cromerian complex. These results are in agreement with the 557 ± 3 ka (±12 ka, including all errors) ⁴⁰Ar/³⁹Ar age obtained from an interbedded tephra layer emitted by the Mont-Dore/Sancy strato-volcano and establish that this sequence probably covers the MIS 15 substages.     

Alluvial architecture of mid-channel fluvial–tidal barforms: The mesotidal Lower Columbia River,Oregon/Washington,USA

Barforms of mesotidal to macrotidal fluvial–tidal transitions, regardless of fluvial-discharge, are currently thought to display a sedimentary architecture dominated by tidal signatures. Due to the scarcity of observations from modern mesotidal fluvial–tidal transitions, especially those of multi-channelled large-rivers (mean annual discharge ≥7000 m³ s⁻¹ and peak discharges ≥15 000 m³ s⁻¹) with mid-channel bars, this concept remains unproven. The present study analyses data produced by a combination of high-resolution ground penetrating radar and coupled shallow vibracores (<5 m depth), collected from modern fluvial–tidal mid-channel bars of the mesotidal multi-channelled Lower Columbia River, Washington/Oregon, USA, which can experience peak discharges ≥18 000 m³ s⁻¹. These data were used alongside time-sequenced aerial imagery to characterize the spatio-temporal sedimentological evolution of these barforms in singular flows or combined flows consisting of river, tidal and/or wind-wave oscillatory, current components operating in unique fluvial–tidal transition regimes. Results indicate that ca 75% of the Lower Columbia River fluvial–tidal transition produces braid-bars with basal to bar-top sedimentological architectures that are indistinguishable from fluvial-only braid-bars recorded in the literature. Barform stratal characteristics within the fluvial–tidal transitions of mesotidal large-rivers are therefore more likely to be dominated by downstream-oriented currents. Furthermore, a new style of low-angle (<5°) inclined heterolithic stratification found in bar-top accretion-sets within upper-mixed tidal–fluvial regime braid-bars is observed. This common stratification is created by combined-flows characterized by intrabasinal wind-wave oscillatory-currents and bidirectional tidal-currents. This inclined heterolithic stratification marks the initial downstream fluvial–tidal crossover point from Lower Columbia River up-dip fully-fluvial braid-bar architectures, to those possessing bar-top facies produced by the hydraulic-sedimentation response of combined intrabasinal wind-wave and tidal influence. When preserved, this form of mid-channel bar inclined heterolithic stratification provides a unique sedimentological signature of multi-channelled fluvial–tidal transitions that possess an open-water lower basin with intrabasinal wind-waves.     

Thermal Variability in a Tidal River

In this paper, we discuss observations of temperature variability in the tidal portion of the San Joaquin River in California. The San Joaquin River makes up the southern portion of the Sacramento San Joaquin Delta, the eastern end of San Francisco Bay. Observations made in August 2004 and August 2005 show significant diurnal variations in temperature in response to surface heat exchange. However, to account for observed changes in heat content a sizeable downstream heat flux (approximately 100 W m⁻²) must be added to the surface heat flux. To account for this flux via Fickian dispersion, a flow-dependent dispersion coefficient varying from 500 to 4,000 m² s⁻¹ is needed. These values are much larger than would be predicted for a river of this size, suggesting that the complex topology of the Delta greatly enhances longitudinal dispersion. Building on these observations, we present a simple theory that explores how the subtidal temperature field varies in response to changes in flow rate, dispersion, and heat exchange.     

Lithospheric flexure, uplift and expected horizontal strain rate in the Pannonian Carpathian region

The Apuseni Mountains are located between the Pannonian Basin and the Transylvanian Basin along a direction of SE convergence with the Carpathian belt. A flexural model based on the cylindrical bending of a semi-infinite, isostatically supported, thin elastic plate is here examined with the Apuseni playing the role of flexural bulge, and under the assumption that the plate is deforming under the action of a vertical shear force and a bending moment applied at the end of the plate, beneath the Carpathians. The model yields estimates of the plate thickness ranging between 13 and 14.5 km, depending on the assumed density contrast between crust/sediments and mantle providing buoyancy. The vertical shear force which is necessary to bend the plate is in the range between 60 and 300 × 10¹¹ N m^− 1, depending on the assumed density contrast. This force is shown to be modelled by a gravitational ‘slab pull’ force, using model parameters derived from seismic tomography. If the height of the flexural bulge, after correction for erosion, is allowed to increase, the model yields an estimate of the horizontal strain rate at the top of the bulge. For example, 5 mm/yr vertical change of the flexural bulge of a 14 km thick plate results in a horizontal deformation rate of approximately 7 nanostrain/yr at the top of the bulge, a value which is at the threshold of sensitivity of continuous GPS measurements. Different vertical rates will change the horizontal strain rate almost proportionally.     

Iron Distribution in Size‐Resolved Aerosols Generated by UV‐Femtosecond Laser Ablation: Influence of Cell Geometry and Implications for In Situ Isotopic Determination by LA‐MC‐ICP‐MS

The influence of ablation cell geometry (Frames single‐ and HelEx two‐volume cells) and laser wavelength (198 and 266 nm) on aerosols produced by femtosecond laser ablation (fs‐LA) were evaluated. Morphologies, iron mass distribution (IMD) and ⁵⁶Fe/⁵⁴Fe ratios of particles generated from magnetite, pyrite, haematite and siderite were studied. The following two morphologies were identified: spherules (10–200 nm) and agglomerates (5–10 nm). Similarity in IMD and ablation rate at 198 and 266 nm indicates similar ablation mechanisms. ⁵⁶Fe/⁵⁴Fe ratios increased with aerodynamic particle size as a result of kinetic fractionation during laser plasma plume expansion, cooling and aerosol condensation. The HelEx cell produces smaller particles with a larger range of ⁵⁶Fe/⁵⁴Fe ratios (1.85‰) than particles from the Frames cell (1.16‰), but the bulk aerosol matches the bulk substrate for both cells, demonstrating stoichiometric fs‐LA sampling. IMD differences are the result of faster wash out of the HelEx cell allowing less time for agglomeration of small, low‐δ ⁵⁶Fe particles with larger, high‐δ ⁵⁶Fe particles in the cell. Even with a shorter ablation time, half the total Fe ion intensity, and half the ablation volume, the HelEx cell produced Fe isotope determinations for magnetite that were as precise as the Frames cell, even when the latter included an aerosol‐homogenising mixing chamber. The HelEx cell delivered a more constant stream of small particles to the ICP, producing a more stable Fe ion signal (0.7% vs. 1.5% RSE for ⁵⁶Fe in a forty‐cycle single analysis), constant instrumental mass bias and thus a more precise measurement.