Solute transport in crystalline rocks at Äspö — II: Blind predictions, inverse modelling and lessons learnt from test STT1

Based on the results from detailed structural and petrological characterisation and on up-scaled laboratory values for sorption and diffusion, blind predictions were made for the STT1 dipole tracer test performed in the Swedish A¨ spo¨ Hard Rock Laboratory. The tracers used were nonsorbing, such as uranine and tritiated water, weakly sorbing 22Na+, 85Sr2 +, 47Ca2 +and more strongly sorbing 86Rb+, 133Ba2 +, 137Cs+. Our model consists of two parts: (1) a flow part based on a 2D-streamtube formalism accounting for the natural background flow field and with an underlying homogeneous and isotropic transmissivity field and (2) a transport part in terms of the dual porosity medium approach which is linked to the flow part by the flow porosity. The calibration of the model was done using the data from one single uranine breakthrough (PDT3). The study clearly showed that matrix diffusion into a highly porous material, fault gouge, had to be included in our model evidenced by the characteristic shape of the breakthrough curve and in line with geological observations. After the disclosure of the measurements, it turned out that, in spite of the simplicity of our model, the prediction for the nonsorbing and weakly sorbing tracers was fairly good. The blind prediction for the more strongly sorbing tracers was in general less accurate. The reason for the good predictions is deemed to be the result of the choice of a model structure strongly based on geological observation. The breakthrough curves were inversely modelled to determine in situ values for the transport parameters and to draw consequences on the model structure applied. For good fits, only one additional fracture family in contact with cataclasite had to be taken into account, but no new transport mechanisms had to be invoked. The in situ values for the effective diffusion coefficient for fault gouge are a factor of 2–15 larger than the laboratory data. For cataclasite, both data sets have values comparable to laboratory data. The extracted Kd values for the weakly sorbing tracers are larger than Swedish laboratory data by a factor of 25–60, but agree within a factor of 3–5 for the more strongly sorbing nuclides. The reason for the inconsistency concerning Kds is the use of fresh granite in the laboratory studies, whereas tracers in the field experiments interact only with fracture fault gouge and to a lesser extent with cataclasite both being mineralogically very different (e.g. clay-bearing) from the intact wall rock.

Coupled W–Os–Pt isotope systematics in IVB iron meteorites: In situ neutron dosimetry for W isotope chronology

Tungsten isotope compositions of magmatic iron meteorites yield ages of differentiation that are within ±2 Ma of the formation of CAIs, with the exception of IVB irons that plot to systematically less radiogenic compositions yielding erroneously old ages. Secondary neutron capture due to galactic cosmic ray (GCR) irradiation is known to lower the ε182W of iron meteorites, adequate correction of which requires a measure of neutron dosage which has not been available, thus far. The W, Os and Pt isotope systematics of 12 of the 13 known IVB iron meteorites were determined by MC-ICP-MS (W, Os, Pt) and TIMS (Os). On the same dissolutions that yield precise ε182W, stable Os and Pt isotopes were determined as in situ neutron dosimeters for empirical correction of the ubiquitous cosmic-ray induced burn-out of 182W in iron meteorites. The W isotope data reveal a main cluster with ε182W of ∼−3.6, but a much larger range than observed in previous studies including irons (Weaver Mountains and Warburton Range) that show essentially no cosmogenic effect on their ε182W. The IVB data exhibits resolvable negative anomalies in ε189Os (−0.6ε) and complementary ε190Os anomalies (+0.4ε) in Tlacotepec due to neutron capture on 189Os which has approximately the same neutron capture cross section as 182W, and captures neutrons to produce 190Os. The least irradiated IVB iron, Warburton Range, has ε189Os and ε190Os identical to terrestrial values. Similarly, Pt isotopes, which are presented as ε192Pt, ε194Pt and ε196Pt range from +4.4ε to +53ε, +1.54ε to −0.32ε and +0.73ε to −0.20ε, respectively, also identify Tlacotepec and Dumont as the most GCR-damaged samples. In W–Os and W–Pt isotope space, the correlated isotope data back-project toward a 0-epsilon value of ε192Pt, ε189Os and ε190Os from which a pre-GCR irradiation ε182W of −3.42±0.09 (2σ) is derived. This pre-GCR irradiation ε182W is within uncertainty of the currently accepted CAI initial ε182W. The Pt and Os isotope correlations in the IVB irons are in good agreement with a nuclear model for spherical irons undergoing GCR spallation, although this model over-predicts the change of ε182W by ∼2×, indicating a need for better W neutron capture cross section determinations. A nucleosynthetic effect in ε184W in these irons of −0.14±0.08 is confirmed, consistent with the presence of Mo and Ru isotope anomalies in IVB irons. The lack of a non-GCR Os isotope anomaly in these irons requires more complex explanations for the production of W, Ru and Mo anomalies than nebular heterogeneity in the distribution of s-process to r-process nuclides.

Cosmogenic production rates and recoil loss effects in micrometeorites and interplanetary dust particles

We present a purely physical model to determine cosmogenic production rates for noble gases and radionuclides in micrometeorites (MMs) and interplanetary dust particles (IDPs) by solar cosmic-rays (SCR) and galactic cosmic-rays (GCR) fully considering recoil loss effects. Our model is based on various nuclear model codes to calculate recoil cross sections, recoil ranges, and finally the percentages of the cosmogenic nuclides that are lost as a function of grain size, chemical composition of the grain, and the spectral distribution of the projectiles. The main advantage of our new model compared with earlier approaches is that we consider the entire SCR particle spectrum up to 240 MeV and not only single energy points. Recoil losses for GCR-produced nuclides are assumed to be equal to recoil losses for SCR-produced nuclides. Combining the model predictions with Poynting-Robertson orbital lifetimes, we calculate cosmic-ray exposure ages for recently studied MMs, cosmic spherules, and IDPs. The ages for MMs and the cosmic-spherule are in the range <2.2–233 Ma, which corresponds, according to the Poynting-Robertson drag, to orbital distances in the range 4.0–34 AU. For two IDPs, we determine exposure ages of longer than 900 Ma, which corresponds to orbital distances larger than 150 AU. The orbital distance in the range 4–6 AU for one MM and the cosmic spherule indicate an origin either in the asteroid belt or release from comets coming either from the Kuiper Belt or the Oort Cloud. Three of the studied MMs have orbital distances in the range 23–34 AU, clearly indicating a cometary origin, either from short-period comets from the Kuiper Belt or from the Oort Cloud. The two IDPs have orbital distances of more than 150 AU, indicating an origin from Oort Cloud comets.

Chronostratigraphy of the 600,000 year old continental record of Lake Van

Lake Van sediment cores from the Ahlat Ridge and Northern Basin drill sites of the ICDP project PALEOVAN contain a wealth of information about past environmental processes. The sedimentary sequence was dated using climatostratigraphic alignment, varve chronology, tephrostratigraphy, argon-argon single-crystal dating, radiocarbon dating, magnetostratigraphy, and cosmogenic nuclides. Based on the lithostratigraphic framework, the different age constraints are compiled and a robust and precise chronology of the 600,000 year-old Lake Van record is constructed. Proxy records of total organic carbon content and sediment color, together with the calcium/potassium-ratios and arboreal pollen percentages of the 174-meter-long Ahlat Ridge record, mimic the Greenland isotope stratotype (NGRIP). Therefore, the proxy records are systematically aligned to the onsets of interstadials reflected in the NGRIP or synthesized Greenland ice-core stratigraphy. The chronology is constructed using 27 age control points derived from visual synchronization with the GICC05 timescale, an absolutely-dated speleothem record (e.g., Hulu, Sanbao, Linzhu cave) and the Epica Dome C timescale. In addition, the uppermost part of the sequence is complemented with four ages from Holocene varve chronology and two calibrated radiocarbon ages. Furthermore, nine argon-argon ages and a comparison of the relative paleointensity record of the magnetic field with reference curve PISO-1500 confirm the accuracy of the age model. Also the identification of the Laschamp event via measurements of 10Be in the sediment confirms the presented age model. The chronology of the Ahlat Ridge record is transferred to the 79-meter-long event-corrected composite record from the Northern Basin and supplemented by additional radiocarbon dating on organic marco-remains. The basal age of the Northern Basin record is estimated at ~90 ka. The variations of the time series of total organic carbon content, the Ca/K ratio, and the arboreal pollen percentages illustrate that the presented chronology and paleoclimate data are suited for reconstructions and modeling of the Quaternary and Pleistocene climate evolution in the Near East at millennial timescales. Furthermore, the chronology of the last 250 kyr can be used to test other dating techniques.

Folate Receptor Targeted Alpha-Therapy Using Terbium-149

Terbium-149 is among the most interesting therapeutic nuclides for medical applications. It decays by emission of short-range α-particles (Eα = 3.967 MeV) with a half-life of 4.12 h. The goal of this study was to investigate the anticancer efficacy of a 149Tb-labeled DOTA-folate conjugate (cm09) using folate receptor (FR)-positive cancer cells in vitro and in tumor-bearing mice. 149Tb was produced at the ISOLDE facility at CERN. Radiolabeling of cm09 with purified 149Tb resulted in a specific activity of ~1.2 MBq/nmol. In vitro assays performed with 149Tb-cm09 revealed a reduced KB cell viability in a FR-specific and activity concentration-dependent manner. Tumor-bearing mice were injected with saline only (group A) or with 149Tb-cm09 (group B: 2.2 MBq; group C: 3.0 MBq). A significant tumor growth delay was found in treated animals resulting in an increased average survival time of mice which received 149Tb-cm09 (B: 30.5 d; C: 43 d) compared to untreated controls (A: 21 d). Analysis of blood parameters revealed no signs of acute toxicity to the kidneys or liver in treated mice over the time of investigation. These results demonstrated the potential of folate-based α-radionuclide therapy in tumor-bearing mice.

Snow shielding factors for cosmogenic nuclide dating inferred from long-term neutron detector monitoring

The depth-dependent attenuation of the secondary cosmic-ray particle flux due to snow cover and its effects on production rates of cosmogenic nuclides constitutes a potential source of uncertainty for studies conducted in regions characterized by frequent seasonal snow burial. Recent experimental and numerical modelling studies have yielded new constraints on the effect of hydrogen-rich media on the production rates of cosmogenic nuclides by low- and high-energy neutrons (<10(-3) MeV and >10(2) MeV, respectively). Here we present long-term neutron-detector monitoring data from a natural setting that we use to quantify the effect of snow cover on the attenuation of fast neutrons (0.1-10 MeV), which are responsible for the production of Ne-21 from Mg and Cl-36 from K. We use data measured between July 2001 and May 2008 at seven stations located throughout the Ecrins-Pelvoux massif (French Western Alps) and its surroundings, at elevations ranging from 200 to 2500 m a.s.l. From the cosmic-ray fluxes recorded during summer, when snow is absent, we infer an apparent attenuation length of 148 g cm(-2) in the atmosphere at a latitude of similar to 45 degrees N and for altitudes ranging from similar to 200 to 2500 m a.s.l. Using snow water-equivalent (SWE) values obtained through snow-coring campaigns that overlap in time the neutron monitoring for five stations, we show that fast neutrons are much more strongly attenuated in snow than predicted by a conventional mass-shielding formulation and the attenuation length estimated in the atmosphere. We suggest that such strong attenuation results from boundary effects at the atmosphere/snow interface induced by the high efficiency of water as a neutron moderator. Finally, we propose an empirical model that allows calculating snow-shielding correction factors as a function of SWE for studies using Ne-21 and Cl-36 analyses in Mg- and K-rich minerals, respectively. This empirical model is of interest for studies with a focus on cosmic-ray exposure dating, particularly if the target rocks are made up of mafic to ultramafic units where seasonal snow-cover is a common phenomenon.

Decomposition studies of group 6 hexacarbonyl complexes. Part 1: Production and decomposition of Mo(CO)6 and W(CO)6

Chemical studies of superheavy elements require fast and efficient techniques, due to short half-lives and low production rates of the investigated nuclides. Here, we advocate for using a tubular flow reactor for assessing the thermal stability of the Sg carbonyl complex – Sg(CO)6. The experimental setup was tested with Mo and W carbonyl complexes, as their properties are established and supported by theoretical predictions. The suggested approach proved to be effective in discriminating between the thermal stabilities of Mo(CO)6 and W(CO)6. Therefore, an experimental verification of the predicted Sg–CO bond dissociation energy seems to be feasible by applying this technique. By investigating the effect of 104,105Mo beta-decay on the formation of 104,105Tc carbonyl complex, we estimated the lower reaction time limit for the metal carbonyl synthesis in the gas phase to be more than 100 ms. We examined further the influence of the wall material of the recoil chamber, the carrier gas composition, the gas flow rate, and the pressure on the production yield of 104Mo(CO)6, so that the future stability tests with Sg(CO)6 can be optimized accordingly.

Radiochemical determination of 129I and 36Cl in MEGAPIE, a proton irradiated lead-bismuth eutectic spallation target

The concentrations of the long-lived nuclear reaction products 129I and 36Cl have been measured in samples from the MEGAPIE liquid metal spallation target. Samples from the bulk target material (lead-bismuth eutectic, LBE), from the interface of the metal free surface with the cover gas, from LBE/steel interfaces and from noble metal absorber foils installed in the cover gas system were analysed using Accelerator Mass Spectrometry at the Laboratory of Ion beam Physics at ETH Zürich. The major part of 129I and 36Cl was found accumulated on the interfaces, particularly at the interface of LBE and the steel walls of the target container, while bulk LBE samples contain only a minor fraction of these nuclides. Both nuclides were also detected on the absorber foils to a certain extent (≪ 1% of the total amount). The latter number is negligible concerning the radio-hazard of the irradiated target material; however it indicates a certain affinity of the absorber foils for halogens, thus proving the principle of using noble metal foils for catching these volatile radionuclides. The total amounts of 129I and 36Cl in the target were estimated from the analytical data by averaging within the different groups of samples and summing up these averages over the total target. This estimation could account for about half of the amount of 129I and 36Cl predicted to be produced using nuclear physics modelling codes for both nuclides. The significance of the results and the associated uncertainties are discussed.

Cosmic-ray exposure ages of pallasites

We analyzed cosmogenic nuclides in metal and/or silicate (primarily olivine) separated from the main-group pallasites Admire, Ahumada, Albin, Brahin, Brenham, Esquel, Finmarken, Glorieta Mountain, Huckitta, Imilac, Krasnojarsk, Marjalahti, Molong, Seymchan, South Bend, Springwater, and Thiel Mountains and from Eagle Station. The metal separates contained an olivine fraction which although small, <1 wt% in most cases, nonetheless contributes significantly to the budgets of some nuclides (e.g., up to 35% for Ne-21 and Al-26). A correction for olivine is therefore essential and was made using model calculations and/or empirical relations for the production rates of cosmogenic nuclides in iron meteoroids and/or measured elemental concentrations. Cosmic-ray exposure (CRE) ages for the metal phases of the main-group pallasites range from 7 to 180 Ma, but many of the ages cluster around a central peak near 100 Ma. These CRE ages suggest that the parent body of the main-group pallasites underwent a major break-up that produced most of the meteorites analyzed. The CRE age distribution for the pallasites overlaps only a small fraction of the distribution for the IIIAB iron meteorites. Most pallasites and IIIAB irons originated in different collisions, probably on different parent bodies; a few IIIABs and pallasites may have come out of the same collision but a firm conclusion requires further study. CRE ages calculated from noble gas and radionuclide data of the metal fraction are higher on average than the Ne-21 exposure ages obtained for the olivine samples. As the metal and olivine fractions were taken in most cases from different specimens, the depth-dependency of the production rate ratio Be-10/Ne-21 in metal, not accounted for in our calculations, may explain the difference.

Beryllium-10 and aluminium-26 concentration and production rates in glacial deposits from the Quartermain Mountains

Many glacial deposits in the Quartermain Mountains, Antarctica present two apparent contradictions regarding the degradation of unconsolidated deposits. The glacial deposits are up to millions of years old, yet they have maintained their meter-scale morphology despite the fact that bedrock and regolith erosion rates in the Quartermain Mountains have been measured at 0.1-4.0 m/Ma. Additionally, ground ice persists in some Miocene-aged soils in the Quartermain Mountains even though modeled and measured sublimation rates of ice in Antarctic soils suggest that without any recharge mechanisms ground ice should sublimate in the upper few meters of soil on the order of 10**3 to 10**5 years. This paper presents results from using the concentration of cosmogenic nuclides beryllium-10 (10Be) and aluminum-26 (26Al) in bulk sediment samples from depth profiles of three glacial deposits in the Quartermain Mountains. The measured nuclide concentrations are lower than expected for the known ages of the deposits, erosion alone does not always explain these concentrations, and deflation of the tills by the sublimation of ice coupled with erosion of the overlying till can explain some of the nuclide concentration profiles. The degradation rates that best match the data range 0.7-12 m/Ma for sublimation of ice with initial debris concentrations ranging 12-45% and erosion of the overlying till at rates of 0.4-1.2 m/Ma. Overturning of the tills by cryoturbation, vertical mixing, or soil creep is not indicated by the cosmogenic nuclide profiles, and degradation appears to be limited to within a few centimeters of the surface. Erosion of these tills without vertical mixing may partially explain how some glacial deposits in the Quartermain Mountains maintain their morphology and contain ground ice close to the surface for millions of years.

Long term low latitude and high elevation cosmogenic 3He production rate inferred from a 107 ka-old lava flow in northern Chile; 22°S-3400 m a.s.l

Available geological calibration sites used to estimate the rate at which cosmogenic 3He is produced at the Earth’s surface are mostly clustered in medium to high latitudes. Moreover, most of them have exposure histories shorter than tens of thousands of years. This lack of sites prevents a qualitative assessment of available production models used to convert cosmogenic 3He concentrations into exposure ages and/or denudation rates. It thus limits our ability to take into account the atmospheric, geomagnetic and solar modulation conditions that might have affected the production of cosmogenic nuclides in the past for longer exposure histories and in low latitude regions. We present the cosmogenic 3He production rate inferred from a new geological calibration site located in northern Chile. Five samples were collected on the surface of the largest and best-preserved lava flow of the San Pedro volcano (21.934°S-68.510°W- 3390 m a.s.l), which displays pristine crease-structure features. 40Ar/39Ar dating yield a reliable plateau age of 107±12 ka for the eruption of this lava flow. Eight pyroxene aliquots separated from the surface samples yield a weighted average cosmogenic 3He concentration of 99.3±1.2 Mat.g-1 from which a local cosmogenic 3He production rate of 928±101 at.g-1.yr-1 is calculated. The local production rate is then scaled to a sea level high latitude (SLHL) reference position using different combinations of geographic spatialization schemes, atmosphere models and geomagnetic field reconstructions, yielding SLHL production rates between 103±11 and 130±14 at.g-1.yr-1 consistent with the most recent estimates available from the literature. Finally, we use the same scaling frameworks to re-evaluate the mean global-scale cosmogenic 3He production rate in olivine and pyroxene minerals at 120±16 at.g-1.yr-1 from the compilation of previously published calibration datasets.

Uncertainty methods inactivation and inventory calculations

For a number of important nuclides, complete activation data libraries with covariance data will be produced, so that uncertainty propagation in fuel cycle codes (in this case ACAB,FISPIN, ...) can be developed and tested. Eventually, fuel inventory codes should be able to handle the complete set of uncertainty data, i.e. those of nuclear reactions (cross sections, etc.), radioactive decay and fission yield data. For this, capabilities will be developed both to produce covariance data and to propagate the uncertainties through the inventory calculations.