Dynamics of supercooled water in highly compacted clays studied by neutron scattering

The freezing behavior of water confined in compacted charged and uncharged clays (montmorillonite in Na-and Ca-forms, illite in Na-and Ca-forms, kaolinite and pyrophyllite) was investigated by neutron scattering. Firstly, the amount of frozen (immobile) water was measured as a function of temperature at the IN16 backscattering spectrometer, Institute Laue-Langevin (ILL). Water in uncharged, partly hydrophobic (kaolinite) and fully hydrophobic (pyrophyllite) clays exhibited a similar freezing and melting behavior to that of bulk water. In contrast, water in charged clays which are hydrophilic could be significantly supercooled. To observe the water dynamics in these clays, further experiments were performed using quasielastic neutron scattering. At temperatures of 250, 260 and 270 K the diffusive motion of water could still be observed, but with a strong reduction in the water mobility as compared with the values obtained above 273 K. The diffusion coefficients followed a non-Arrhenius temperature dependence well described by the Vogel-Fulcher-Tammann and the fractional power relations. The fits revealed that Na-and Ca-montmorillonite and Ca-illite have similar Vogel-Fulcher-Tammann temperatures (T-VFT, often referred to as the glass transition temperature) of similar to 120 K and similar temperatures at which the water undergoes the 'strong-fragile' transition, T-s similar to 210 K. On the other hand, Na-illite had significantly larger values of T-VFT similar to 180 K and T-s similar to 240 K. Surprisingly, Ca-illite has a similar freezing behavior of water to that of montmorillonites, even though it has a rather different structure. We attribute this to the stronger hydration of Ca ions as compared with the Na ions occurring in the illite clays.

Strategies for reducing fumigant loss to the atmosphere

A model is developed to describe transport and loss of methyl bromide (MeBr) in soil following application as a soil fumigant. The model is used to investigate the effect of soil and management factors on MeBr volatilization. Factors studied include depth of injection, soil water content, presence or absence of tarp, depth to downward barrier, and irrigation after injection. Of these factors, the most important was irrigation after injection followed by covering with the tarp, which increased the diffusive resistance of the soil and prevented early loss of MeBr. The model offers an explanation for the apparently contradictory observations of earlier field studies of MeBr volatilization from soils under different conditions. The model was also used to calculate the concentration-time index for various management alternatives, showing that the irrigation application did not make the surface soil more difficult to fumigate, except at very early times. Therefore, irrigation shows promise for reducing fumigant loss while at the same time permitting control of target organisms during fumigation.

Transport of 234U in the Opalinus Clay on centimetre to decimetre scales

The Opalinus Clay formation in North Switzerland is a potential host rock for a deep underground radioactive waste repository. The distribution of U-238, U-234 and Th-230 was studied in rock samples of the Opalinus Clay from an exploratory borehole at Benken (Canton of Zurich) using MC-ICP-MS. The aim of U-234 was to assess the in situ, long-term migration behaviour in this rock. Very low hydraulic conductivities of the Opalinus Clay, reducing potential of the pore water and its chemical equilibrium with the host rock are expected to render both U-238 and Th-230 immobile. If U is heterogeneously distributed in the Opalinus Clay, gradients in the supply of U-234 from the rock matrix to the pore water by the decay of U-238 will be established. Diffusive redistribution separates U-234 from its immobile parent U-238 resulting in bulk rock U-234/U-238 activity disequilibria. These may provide a means of estimating the mobility of U-234 in the rock if the diffusion rate of U-234 is significant compared to its decay rate. Sampling was carried out on two scales. Drilling of cm-spaced samples from the drill-core was done to study mobility over short distances and elucidate possible small-scale lithological control. Homogenized 25-cm-long portions of a 2-m-long drill-core section were prepared to provide information on transport over a longer distance. Variations in U and/or Th content on the cm-scale between clays and carbonate-sandy layers are revealed by beta-scanning, which shows that the (dominant) clay is richer in both elements. Samples were digested using aqua regia followed by total HF dissolution, yielding two fractions. in all studied samples U was found to be concentrated in the HF digestion fraction. It has a high U/Th ratio and a study by SEM-EDS points to sub-mu m up to several mu m in size zircon grains as the main U-rich phase. This fraction consistently has U-234/U-238 activity ratios below unity. The minute zircon grains constitute the major reservoir of U in the rock and act as constant rate suppliers of U-234 into the rock matrix and the pore water. The aqua regia leach fraction was found to be enriched in Th, and complementary to the HF fraction, having U-234/U-238 activity ratios above unity. It is believed that these U activity ratios reflect the surplus of having U-234 delivered from the zircon grains. Some cm-spaced samples show bulk rock U-234/U-238 activity ratios that are markedly out of equilibrium. In most of them a striking negative correlation between the total U content and the bulk rock U-234/U-238 activity ratios is observed. This is interpreted to indicate net U-234 transfer from regions of higher supply of U-234 towards those of lower supply which is, in most cases, equivalent to transfer from clayey towards carbonate/sandy portions of the rock. In contrast, the 25 cm averaged samples all have uniform bulk rock U-234/U-238 activity ratios in equilibrium, indicating U immobility in the last 1-1.5 Ma on this spatial scale. It is concluded that the small-scale lithological variations which govern U spatial distribution in the Opalinus Clay are the major factor determining U-234 in situ supply rates, regulating its diffusive fluxes and controlling the observed bulk rock U-234/U-238 activity ratios. A simple box-model is presented to simulate the measured bulk rock U-234/U-238 activity ratios and to give an additional insight into the studied system. (C) 2008 Elsevier Ltd. All rights reserved.

Light noble gases in 12 meteorites from the Omani desert, Australia, Mauritania, Canada, and Sweden

We measured the concentrations and isotopic compositions of He, Ne, and Ar in 14 fragments from 12 different meteorites: three carbonaceous chondrites, six L chondrites (three most likely paired), one H chondrite, one R chondrite, and one ungrouped chondrite. The data obtained for the CV3 chondrites Ramlat as Sahmah (RaS) 221 and RaS 251 support the hypothesis of exposure age peaks for CV chondrites at approximately 9 Ma and 27 Ma. The exposure age for Shişr 033 (CR chondrite) of 7.3 Ma is also indicative of a possible CR chondrite exposure age peak. The three L chondrites Jiddat al Harasis (JaH) 091, JaH 230, and JaH 296, which are most likely paired, fall together with Hallingeberg into the L chondrite exposure age peak of approximately 15 Ma. The two L chondrites Shelburne and Lake Torrens fall into the peaks at approximately 40 Ma and 5 Ma, respectively. The ages for Bassikounou (H chondrite) and RaS 201 (R chondrite) are approximately 3.5 Ma and 5.8 Ma, respectively. Six of the studied meteorites show clear evidence for 3He diffusive losses, the deficits range from approximately 17% for one Lake Torrens aliquot to approximately 45% for RaS 211. The three carbonaceous chondrites RaS 221, RaS 251, and Shişr 033 all have excess 4He, either of planetary or solar origin. However, very high 4He/20Ne ratios occur at relatively low 20Ne/22Ne ratios, which is unexpected and needs further study. The measured 40Ar ages fit well into established systematics. They are between 2.5 and 4.5 Ga for the carbonaceous chondrites, older than 3.6 Ga for the L and H chondrites, and about 2.4 Ga for the R chondrite as well as for the ungrouped chondrite. Interestingly, none of our studied L chondrites has been degassed in the 470 Ma break-up event. Using the amount of trapped 36Ar as a proxy for noble gas contamination due to terrestrial weathering we are able to demonstrate that the samples studied here are not or only very slightly affected by terrestrial weathering (at least in terms of their noble gas budget).

Calibration of cosmogenic noble gas production in ordinary chondrites based on 36 Cl- 36 Ar ages. Part 1: Refined produced rates for cosmogenic 21 Ne and 38 Ar

We measured the concentrations and isotopic compositions of He, Ne, and Ar in bulk samples and metal separates of 14 ordinary chondrite falls with long exposure ages and high metamorphic grades. In addition, we measured concentrations of the cosmogenic radionuclides 10Be, 26Al, and 36Cl in metal separates and in the nonmagnetic fractions of the selected meteorites. Using cosmogenic 36Cl and 36Ar measured in the metal separates, we determined 36Cl-36Ar cosmic-ray exposure (CRE) ages, which are shielding-independent and therefore particularly reliable. Using the cosmogenic noble gases and radionuclides, we are able to decipher the CRE history for the studied objects. Based on the correlation 3He/21Ne versus 22Ne/21Ne, we demonstrate that, among the meteorites studied, only one suffered significant diffusive losses (about 35%). The data confirm that the linear correlation 3He/21Ne versus 22Ne/21Ne breaks down at high shielding. Using 36Cl-36Ar exposure ages and measured noble gas concentrations, we determine 21Ne and 38Ar production rates as a function of 22Ne/21Ne. The new data agree with recent model calculations for the relationship between 21Ne and 38Ar production rates and the 22Ne/21Ne ratio, which does not always provide unique shielding information. Based on the model calculations, we determine a new correlation line for 21Ne and 38Ar production rates as a function of the shielding indicator 22Ne/21Ne for H, L, and LL chondrites with preatmospheric radii less than about 65 cm. We also calculated the 10Be/21Ne and 26Al/21Ne production rate ratios for the investigated samples, which show good agreement with recent model calculations.

Influence of the background wind on the local soil moisture-precipitation feedback

The importance of soil moisture anomalies on airmass convection over semiarid regions has been recognized in several studies. The underlying mechanisms remain partly unclear. An open question is why wetter soils can result in either an increase or a decrease of precipitation (positive or negative soil moisture–precipitation feedback, respectively). Here an idealized cloud-resolving modeling framework is used to explore the local soil moisture–precipitation feedback. The approach is able to replicate both positive and negative feedback loops, depending on the environmental parameters. The mechanism relies on horizontal soil moisture variations, which may develop and intensify spontaneously. The positive expression of the feedback is associated with the initiation of convection over dry soil patches, but the convective cells then propagate over wet patches where they strengthen and preferentially precipitate. The negative feedback may occur when the wind profile is too weak to support the propagation of convective features from dry to wet areas. Precipitation is then generally weaker and falls preferentially over dry patches. The results highlight the role of the midtropospheric flow in determining the sign of the feedback. A key element of the positive feedback is the exploitation of both low convective inhibition (CIN) over dry patches (for the initiation of convection) and high CAPE over wet patches (for the generation of precipitation).

Thermal Convection in Ice Sheets: We Look But Do Not See

Thermal convection in the Antarctic and Greenland ice sheets has been dismissed on the grounds that radio-echo stratigraphy is undisturbed for long distances. However, the undisturbed stratigraphy lies, for the most part, above the density inversion in polar ice sheets and therefore does not disprove convection. An echo-free zone is widespread below the density inversion, yet nobody has cited this as a strong indication that convection is indeed present at d�pth. A generalized Rayleigh criterion for thermal convection in e1astic-viscoplastic polycrystalline solids heated from below is developed and applied to ice-sheet convection. An infinite Rayleigh number at the onset of primary creep decreases with time and becomes constant when secondary creep dominates, suggesting that any thermal buoyancy stress can initiate convection but convection cannot be sustained below a buoyancy stress of about 3 kPa. An analysis of the temperature profile down the Byrd Station core hole suggests that about 1000 m of ice below the density inversion will sustain convection. Creep along the Byrd Station strain network, radar sounding in East Antarctica, and seismic sounding in West Antarctica are examined for evidence of convective creep superimposed on advective creep. It is concluded that the evidence for convection is there, if we look for it with the intention offinding it.

Reduction in Northward Incursions of the South Asian Monsoon Since ~1400 AD Inferred from a Mt. Everest Ice Core

A highly resolved Mt. Everest ice core reveals a decrease in marine and increase in continental air masses related to relatively high summer surface pressure over Mongolia, and reduction in northward incursions of the summer South Asian monsoon since similar to 1400 AD. Previously published proxy records from lower sites south of the Himalayas indicate strengthening of the monsoon since this time. These regional differences are consistent with a south north seesaw in convective activity in the Asian monsoon region, and reflect a southward shift in the mean summer position of the monsoon trough since similar to 1400 AD. The change in monsoonal circulation at 1400 AD is synchronous with a reduction in solar irradiance and the onset of the LIA. This demonstrates a hemispheric scale circulation reorganization at this time, and the potential for future large shifts in monsoonal circulation.

Temperature-Time History of Subducted Continental-Crust, Mount Olympos Region, Greece

In the Mt. Olympos region of northeastern Greece, continental margin strata and basement rocks were subducted and metamorphosed under blueschist facies conditions, and thrust over carbonate platform strata during Alpine orogenesis. Subsequent exposure of the subducted basement rocks by normal faulting has allowed an integrated study of the timing of metamorphism, its relationship to deformation, and the thermal history of the subducted terrane. Alpine low-grade metamorphic assemblages occur at four structural levels. Three thrust sheets composed of Paleozoic granitic basement and Mesozoic metasedimentary cover were thrust over Mesozoic carbonate rocks and Eocene flysch; thrusting and metamorphism occurred first in the highest thrust sheets and progressed downward as units were imbricated from NE to SW. 40Ar/39Ar spectra from hornblende, white mica, and biotite samples indicate that the upper two units preserve evidence of four distinct thermal events: (1) 293–302 Ma crystallization of granites, with cooling from >550°C to <325°C by 284 Ma; (2) 98–100 Ma greenschist to blueschist-greenschist transition facies metamorphism (T∼350–500°C) and imbrication of continental thrust sheets; (3) 53–61 Ma blueschist facies metamorphism and deformation of the basement and continental margin units at T<350–400°C; (4) 36–40 Ma thrusting of blueschists over the carbonate platform, and metamorphism at T∼200–350°C. Only the Eocene and younger events affected the lower two structural packages. A fifth event, indicated by diffusive loss profiles in microcline spectra, reflects the beginning of uplift and cooling to T<100–150°C at 16–23 Ma, associated with normal faulting which continued until Quaternary time. Incomplete resetting of mica ages in all units constrains the temperature of metamorphism during continental subduction to T≤350°C, the closure temperature for Ar in muscovite. The diffusive loss profiles in micas and K-feldspars enable us to “see through” the younger events to older events in the high-T parts of the release spectra. Micas grown during earlier metamorphic events lost relatively small amounts of Ar during subsequent high pressure-low temperature metamorphism. Release spectra from phengites grown during Eocene metamorphism and deformation record the ages of the Ar-loss events. Alpine deformation in northern Greece occurred over a long time span (∼90 Ma), and involved subduction and episodic imbrication of continental basement before, during, and after the collision of the Apulian and Eurasian plates. Syn-subduction uplift and cooling probably combined with intermittently higher cooling rates during extensional events to preserve the blueschist facies mineral assemblages as they were exhumed from depths of >20 km. Extension in the Olympos region was synchronous with extension in the Mesohellenic trough and the Aegean back-arc, and concurrent with westward-progressing shortening in the external Hellenides.

Water–montmorillonite systems: Neutron scattering and tracer throughdiffusion studies

Designs for deep geological respositories of nuclear waste include bentonite as a hydraulic and chemisorption buffer material to protect the biosphere from leakage of radionuclides. Bentonite is chosen because it is a cheap, naturally occurring material with the required properties. It consists essentially of montmorillonite, a swelling clay mineral. Upon contact with groundwater such clays can seal the repository by incorporating water in the interlayers of their crystalline structure. The intercalated water exhibits significantly different properties to bulk water in the surrounding interparticle pores, such as lower diffusion coefficients (González Sánchez et. al. 2008). This doctoral thesis presents water distribution and diffusion behavior on various time and space scales in montmorillonite. Experimental results are presented for Na- and Cs-montmorillonite samples with a range of bulk dry densities (0.8 to 1.7 g/cm3). The experimental methods employed were neutron scattering (backscattering, diffraction, time-of-flight), adsorption measurements (water, nitrogen) and tracer-through diffusion. For the tracer experiments the samples were fully saturated via the liquid phase under volume-constrained conditions. In contrast, for the neutron scattering experiments, the samples were hydrated via the vapor phase and subsequently compacted, leaving a significant fraction of interparticle pores unfilled with water. Owing to these differences in saturation, the water contents of the samples for neutron scattering were characterized by gravimetry whereas those for the tracer experiments were obtained from the bulk dry density. The amount of surface water in interlayer pores could be successfully discriminated from the amount of bulk-like water in interparticle pores in Na- and Csmontmorillonite using neutron spectroscopy. For the first time in the literature, the distribution of water between these two pore environments was deciphered as a function of gravimetric water content. The amount was compared to a geometrical estimation of the amount of interlayer and interparticle water determined by neutron diffraction and adsorption measurements. The relative abundances of the 1 to 4 molecular water layers in the interlayer were determined from the area ratios of the (001)-diffraction peaks. Depending on the characterization method, different fractions of surface water and interlayer water were obtained. Only surface and interlayer water exists in amontmorillonite with water contents up to 0.18 g/g according to spectroscopic measurements and up to 0.32 g/g according to geometrical estimations, respectively. At higher water contents, bulk-like and interparticle water also exists. The amounts increase monotonically, but not linearly, from zero to 0.33 g/g for bulk-like water and to 0.43 g/g for interparticle water. It was found that water most likely redistributes between the surface and interlayer sites during the spectroscopic measurements and therefore the reported fraction is relevant only below about -10 ºC (Anderson, 1967). The redistribution effect can explain the discrepancy in fractions between the methods. In a novel approach the fractions of water in different pore environments were treated as a fixed parameter to derive local diffusion coefficients for water from quasielastic neutron scattering data, in particular for samples with high water contents. Local diffusion coefficients were obtained for the 1 to 4 molecular water layers in the interlayer of 0.5·10–9, 0.9·10–9, 1.5·10–9 and 1.4·10–9 m²/s, respectively, taking account of the different water fractions (molecular water layer, bulk-like water). The diffusive transport of 22Na and HTO through Na-montmorillonite was measured on the laboratory experimental scale (i.e. cm, days) by tracer through-diffusion experiments. We confirmed that diffusion of HTO is independent of the ionic strength of the external solution in contact with the clay sample but dependent on the bulk dry density. In contrast, the diffusion of 22Na was found to depend on both the ionic strength of the pore solution and on the bulk dry density. The ratio of the pore and surface diffusion could be experimentally determined for 22Na from the dependence of the diffusion coefficient on the ionic strength. Activation energies were derived from the temperaturedependent diffusion coefficients via the Arrhenius relation. In samples with high bulk dry density the activation energies are slightly higher than those of bulk water whereas in low density samples they are lower. The activation energies as a function of ionic strengths of the pore solutions are similar for 22Na and HTO. The facts that (i) the slope of the logarithmic effective diffusion coefficients as a function of the logarithmic ionic strength is less than unity for low bulk dry densities and (ii) two water populations can be observed for high gravimetric water contents (low bulk dry densities) support the interlayer and interparticle porosity model proposed by Glaus et al. (2007), Bourg et al. (2006, 2007) and Gimmi and Kosakowski (2011).

Lake surface temperatures in a changing climate: a global sensitivity analysis

We estimate the effects of climatic changes, as predicted by six climate models, on lake surface temperatures on a global scale, using the lake surface equilibrium temperature as a proxy. We evaluate interactions between different forcing variables, the sensitivity of lake surface temperatures to these variables, as well as differences between climate zones. Lake surface equilibrium temperatures are predicted to increase by 70 to 85 % of the increase in air temperatures. On average, air temperature is the main driver for changes in lake surface temperatures, and its effect is reduced by ~10 % by changes in other meteorological variables. However, the contribution of these other variables to the variance is ~40 % of that of air temperature, and their effects can be important at specific locations. The warming increases the importance of longwave radiation and evaporation for the lake surface heat balance compared to shortwave radiation and convective heat fluxes. We discuss the consequences of our findings for the design and evaluation of different types of studies on climate change effects on lakes.

Analytical Models of Exoplanetary Atmospheres. II. Radiative Transfer via the Two-Stream Approximation

We present a comprehensive analytical study of radiative transfer using the method of moments and include the effects of non-isotropic scattering in the coherent limit. Within this unified formalism, we derive the governing equations and solutions describing two-stream radiative transfer (which approximates the passage of radiation as a pair of outgoing and incoming fluxes), flux-limited diffusion (which describes radiative transfer in the deep interior) and solutions for the temperature-pressure profiles. Generally, the problem is mathematically under-determined unless a set of closures (Eddington coefficients) is specified. We demonstrate that the hemispheric (or hemi-isotropic) closure naturally derives from the radiative transfer equation if energy conservation is obeyed, while the Eddington closure produces spurious enhancements of both reflected light and thermal emission. We concoct recipes for implementing two-stream radiative transfer in stand-alone numerical calculations and general circulation models. We use our two-stream solutions to construct toy models of the runaway greenhouse effect. We present a new solution for temperature-pressure profiles with a non-constant optical opacity and elucidate the effects of non-isotropic scattering in the optical and infrared. We derive generalized expressions for the spherical and Bond albedos and the photon deposition depth. We demonstrate that the value of the optical depth corresponding to the photosphere is not always 2/3 (Milne's solution) and depends on a combination of stellar irradiation, internal heat and the properties of scattering both in optical and infrared. Finally, we derive generalized expressions for the total, net, outgoing and incoming fluxes in the convective regime.

The isotopic composition of pore-water from Tournemire argillite (France): An inter-comparison study.

The determination of stable isotope contents of pore-water from consolidated argillaceous rocks remains a critical issue. In order to understand the processes involved in techniques developed for acquiring stable isotope compositions of pore-water, a comparative study between different methods was based on core samples of the Tournemire argillite. It concerns two water extraction techniques based on vacuum distillation and two pore-water equilibration techniques (radial diffusion in liquid phase and diffusive exchange in vapor phase). The water-content values obtained from vacuum distillation at 50 °C are always the lowest, on average 8% lower than the values obtained by heating at 105 °C and 17% lower than the values obtained by heating at 150 °C. The amounts of pore-water estimated from vacuum distillation at 105 °C and 150 °C and from radial diffusion method are in good agreement with those determined by heating. On the contrary, the vapor exchange method provides the highest values of water contents. Concerning stable isotope data, a good agreement was found between those obtained by equilibration techniques and those of fracture water, especially for 2H. Vacuum distillation at high temperature (particularly at 150 °C) also provided results consistent with data of fracture fluids. On the other hand, distillation at 50 °C provides a systematic depletion in heavy isotopes (about –20‰ for 2H and –2.7‰ for 18O) that can be modelled by an incomplete Rayleigh-type distillation process.

A multi-model comparison of Atlantic multidecadal variability

A multi-model analysis of Atlantic multidecadal variability is performed with the following aims: to investigate the similarities to observations; to assess the strength and relative importance of the different elements of the mechanism proposed by Delworth et al. (J Clim 6:1993–2011, 1993) (hereafter D93) among coupled general circulation models (CGCMs); and to relate model differences to mean systematic error. The analysis is performed with long control simulations from ten CGCMs, with lengths ranging between 500 and 3600 years. In most models the variations of sea surface temperature (SST) averaged over North Atlantic show considerable power on multidecadal time scales, but with different periodicity. The SST variations are largest in the mid-latitude region, consistent with the short instrumental record. Despite large differences in model configurations, we find quite some consistency among the models in terms of processes. In eight of the ten models the mid-latitude SST variations are significantly correlated with fluctuations in the Atlantic meridional overturning circulation (AMOC), suggesting a link to northward heat transport changes. Consistent with this link, the three models with the weakest AMOC have the largest cold SST bias in the North Atlantic. There is no linear relationship on decadal timescales between AMOC and North Atlantic Oscillation in the models. Analysis of the key elements of the D93 mechanisms revealed the following: Most models present strong evidence that high-latitude winter mixing precede AMOC changes. However, the regions of wintertime convection differ among models. In most models salinity-induced density anomalies in the convective region tend to lead AMOC, while temperature-induced density anomalies lead AMOC only in one model. However, analysis shows that salinity may play an overly important role in most models, because of cold temperature biases in their relevant convective regions. In most models subpolar gyre variations tend to lead AMOC changes, and this relation is strong in more than half of the models.

Revisiting the Nd-142 deficits in the 1.48 Ga Khariar alkaline rocks, India

The 146Sm–142Nd system plays a central role in tracing the silicate differentiation of the Earth prior to 4.1 Ga. After this time, given its initial abundance, the 146Sm can be considered to be extinct. Upadhyay et al. (2009) reported unexpected negative 142Nd anomalies in 1.48 Ga rocks of the Khariar nepheline syenite complex (India) and inferred that an early enriched, low-Sm/Nd reservoir must have contributed to the mantle source rocks of the Khariar complex. As 146Sm had been effectively extinct for about 2.6 billion years before the crystallisation of the Khariar samples, this Nd signature should have remained isolated from the convective mantle for at least that long. It was thus suggested that the source rock of Khariar samples had been sequestered in the lithospheric root of the Indian craton. Using a different chemical separation method, and a different Thermal Ionization Mass Spectrometry (TIMS) analysis protocol, the present study attempted to replicate these negative 142Nd anomalies, but none were found. To determine which data set is correct, we investigated three possible sources of bias between them: imperfect cancellation of Faraday collector efficiencies during multidynamic TIMS analysis, rapid sample fractionation between the sequential measurement of 146Nd/144Nd and 142Nd/144Nd, and non-exponential law behaviour resulting from so-called “domain mixing.” Incomplete cancellation of collector efficiencies was found unlikely to cause resolvable biases at the estimated level of variation among collector efficiencies. Even in the case of highly variable efficiency and resolvable biases, there is no reason to suspect that they would reproducibly affect only four rocks out of 10 analysed by Upadhyay et al. (2009). Although domain mixing may explain apparent “reverse” fractionation trends observed in some TIMS analyses, it cannot be the cause of the apparent negative anomalies in the study of Upadhyay et al. (2009). It was determined that rapid mass fractionation during the course of a multidynamic TIMS analysis can bias all measured Nd ratios. After applying an approximate correction for this effect, only one rock from Upadhyay et al. (2009) retained an apparent negative 142Nd anomaly. This, in conjunction with our new, anomaly-free data set measured at fractionation rates too low to cause bias, leads to the conclusion that the anomalies reported by Upadhyay et al. (2009) are a subtle and reproducible analytical artefact. The absence of negative 142Nd anomalies in these rocks relaxes the need for a mechanism (other than crust formation) that can isolate a Nd reservoir from the convective mantle for billions of years.