Late Weichselian local ice dome configuration and chronology in Northwestern Svalbard: early thinning, late retreat

The chronology and configuration of the Svalbard Barents Sea Ice Sheet (SBSIS) during the Late Weichselian (LW) are based on few and geographically scattered data. Thus, the timing and configuration of the SBSIS has been a subject of extensive debate. We present provenance data of erratic boulders and cosmogenic 10Be ages of bedrock and boulders from Northwest Spitsbergen (NWS), Svalbard to determine the thickness, configuration and chronology during the LW. We sampled bedrock and boulders of mountain summits and summit slopes, along with erratic boulders from coastal locations around NWS. We suggest that a local ice dome over central NWS during LW drained radially in all directions. Provenance data from erratic boulders from northern coastal lowland Reinsdyrflya suggest northeastward ice flow through Liefdefjorden. 10Be ages of high-elevation erratic boulders in central NWS (687–836 m above sea level) ranging from 18.3 ± 1.3 ka to 21.7 ± 1.4 ka, indicate that the centre of a local ice dome was at least 300 m thicker than at present. 10Be ages of all high-elevation erratics (>400 m above sea level, central and coastal locations) indicate the onset of ice dome thinning at 25–20 ka. 10Be ages from erratic boulders on Reinsdyrflya ranging from 11.1 ± 0.8 ka to 21.4 ± 1.7 ka, indicate an ice cover over the entire Reinsdyrflya during LW and a complete deglaciation prior to the Holocene, but apparently later than the thinning in the mountains. Lack of moraine deposits, but the preservation of beach terraces, suggest that the ice covering this peninsula possibly was cold-based and that Reinsdyrflya was part of an inter ice-stream area covered by slow-flowing ice, as opposed to the adjacent fjord, which possibly was filled by a fast-flowing ice stream. Despite the early thinning of the ice sheet (25–20 ka) we find a later timing of deglaciation of the fjords and the distal lowlands. Several bedrock samples (10Be) from vertical transects in the central mountains of NWS pre-date the LW, and suggest either ice free or pervasive cold-based ice conditions. Our reconstruction is aligned with the previously suggested hypothesis that a complex multi-dome ice-sheet-configuration occupied Svalbard and the Barents Sea during LW, with numerous drainage basins feeding fast ice streams, separated by slow flowing, possibly cold-based, inter ice-stream areas.

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.

A reconstruction of radiocarbon production and total solar irradiance from the Holocene ¹⁴C and CO₂ records: implications of data and model uncertainties

Radiocarbon production, solar activity, total solar irradiance (TSI) and solar-induced climate change are reconstructed for the Holocene (10 to 0 kyr BP), and TSI is predicted for the next centuries. The IntCal09/SHCal04 radiocarbon and ice core CO2 records, reconstructions of the geomagnetic dipole, and instrumental data of solar activity are applied in the Bern3D-LPJ, a fully featured Earth system model of intermediate complexity including a 3-D dynamic ocean, ocean sediments, and a dynamic vegetation model, and in formulations linking radiocarbon production, the solar modulation potential, and TSI. Uncertainties are assessed using Monte Carlo simulations and bounding scenarios. Transient climate simulations span the past 21 thousand years, thereby considering the time lags and uncertainties associated with the last glacial termination. Our carbon-cycle-based modern estimate of radiocarbon production of 1.7 atoms cm−2 s−1 is lower than previously reported for the cosmogenic nuclide production model by Masarik and Beer (2009) and is more in-line with Kovaltsov et al. (2012). In contrast to earlier studies, periods of high solar activity were quite common not only in recent millennia, but throughout the Holocene. Notable deviations compared to earlier reconstructions are also found on decadal to centennial timescales. We show that earlier Holocene reconstructions, not accounting for the interhemispheric gradients in radiocarbon, are biased low. Solar activity is during 28% of the time higher than the modern average (650 MeV), but the absolute values remain weakly constrained due to uncertainties in the normalisation of the solar modulation to instrumental data. A recently published solar activity–TSI relationship yields small changes in Holocene TSI of the order of 1 W m−2 with a Maunder Minimum irradiance reduction of 0.85 ± 0.16 W m−2. Related solar-induced variations in global mean surface air temperature are simulated to be within 0.1 K. Autoregressive modelling suggests a declining trend of solar activity in the 21st century towards average Holocene conditions.

Ice Elevation Near the West Antarctic Ice Sheet Divide During the Last Glaciation

Interior ice elevations of the West Antarctic Ice Sheet (WAIS) during the last glaciation, which can serve as benchmarks for ice-sheet models, are largely unconstrained. Here we report past ice elevation data from the Ohio Range, located near the WAIS divide and the onset region of the Mercer Ice Stream. Cosmogenic exposure ages of glacial erratics that record a WAIS highstand similar to 125 m above the present surface date to similar to 11.5 ka. The deglacial chronology prohibits an interior WAIS contribution to meltwater pulse 1A. Our observational data of ice elevation changes compare well with predictions of a thermomechanical ice-sheet model that incorporates very low basal shear stress downstream of the present day grounding line. We conclude that ice streams in the Ross Sea Embayment had thin, low-slope profiles during the last glaciation and interior WAIS ice elevations during this period were several hundred meters lower than previous reconstructions.

Terrestrial and submarine evidence for the extent and timing of the Last Glacial Maximum and the onset of deglaciation on the maritime-Antarctic and sub-Antarctic islands

This paper is the maritime and sub–Antarctic contribution to the Scientific Committee for Antarctic Research (SCAR) Past Antarctic Ice Sheet Dynamics (PAIS) community Antarctic Ice Sheet reconstruction. The overarching aim for all sectors of Antarctica was to reconstruct the Last Glacial Maximum (LGM) ice sheet extent and thickness, and map the subsequent deglaciation in a series of 5000 year time slices. However, our review of the literature found surprisingly few high quality chronological constraints on changing glacier extents on these timescales in the maritime and sub–Antarctic sector. Therefore, in this paper we focus on an assessment of the terrestrial and offshore evidence for the LGM ice extent, establishing minimum ages for the onset of deglaciation, and separating evidence of deglaciation from LGM limits from those associated with later Holocene glacier fluctuations. Evidence included geomorphological descriptions of glacial landscapes, radiocarbon dated basal peat and lake sediment deposits, cosmogenic isotope ages of glacial features and molecular biological data. We propose a classification of the glacial history of the maritime and sub–Antarctic islands based on this assembled evidence. These include: (Type I) islands which accumulated little or no LGM ice; (Type II) islands with a limited LGM ice extent but evidence of extensive earlier continental shelf glaciations; (Type III) seamounts and volcanoes unlikely to have accumulated significant LGM ice cover; (Type IV) islands on shallow shelves with both terrestrial and submarine evidence of LGM (and/or earlier) ice expansion; (Type V) Islands north of the Antarctic Polar Front with terrestrial evidence of LGM ice expansion; and (Type VI) islands with no data. Finally, we review the climatological and geomorphological settings that separate the glaciological history of the islands within this classification scheme.

Continental degassing of 4He by surficial discharge of deep groundwater

Radiogenic He is produced by the decay of uranium and thorium in the Earth’s mantle and crust. From here, it is degassed to the atmosphere and eventually escapes to space. Assuming that all of the 4He produced is degassed, about 70% of the total He degassed from Earth comes from the continental crust. However, the outgoing flux of crustal He has not been directly measured at the Earth’s surface and the migration pathways are poorly understood. Here we present measurements of helium isotopes and the long-lived cosmogenic radio-isotope Kr in the deep, continental-scale Guarani aquifer in Brazil and show that crustal He reaches the atmosphere primarily by the surficial discharge of deep groundwater. We estimate that He in Guarani groundwater discharge accounts for about 20% of the assumed global flux from continental crust, and that other large aquifers may account for about 33%. Old groundwater ages suggest that He in the Guarani aquifer accumulates over half- to one-millionyear timescales. We conclude that He degassing from the continents is regulated by groundwater discharge, rather than episodic tectonic events, and suggest that the assumed steady state between crustal production and degassing of He, and its resulting atmospheric residence time, should be re-examined.

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.

Transient sediment supply in a high-altitude Alpine environment evidenced through a 10Be budget of the Etages catchment (French Western Alps)

Although beryllium-10 (10Be) concentrations in stream sediments provide useful synoptic views of catchment-wide erosion rates, little is known on the relative contributions of different sediment supply mechanisms to the acquisition of their initial signature in the headwaters. Here we address this issue by conducting a 10Be-budget of detrital materials that characterize the morphogenetic domains representative of high-altitude environments of the European Alps. We focus on the Etages catchment, located in the Ecrins-Pelvoux massif (southeast France), and illustrate how in situ 10Be concentrations can be used for tracing the origin of the sand fraction from the bedload in the trunk stream. The landscape of the Etages catchment is characterized by a geomorphic transient state, high topographic gradients, and a large variety of modern geomorphic domains ranging from glacial environments to scarcely vegetated alluvial plains. Beryllium-10 concentrations measured in the Etages catchment vary from similar to 1 x 104 to 4.5 x 105 atoms per gram quartz, while displaying consistent 10Be signatures within each representative morphogenetic unit. We show that the basic requirements for inferring catchment-wide denudation from 10Be concentration measurements are not satisfied in this small, dynamic catchment. However, the distinct 10Be signature observed for the geomorphic domains can be used as a tracer. We suggest that a terrestrial cosmogenic nuclide (TCN) budget approach provides a valuable tool for the tracing of material origin in basins where the let nature do the averaging' principles may be violated.

Timing of retreat of the Reuss Glacier (Switzerland) at the end of the Last Glacial Maximum

We used cosmogenic 10Be and 36Cl to establish the timing of the onset of deglaciation after the Last Glacial Maximum of the Reuss Glacier, one of the piedmont lobes of the Alpine ice cap that reached the northern Alpine foreland in Switzerland. In this study, we sampled erratic boulders both at the frontal position in the foreland (Lenzburg and Wohlen, canton Aargau) and at the lateral Alpine border position (Seeboden moraine, Rigi, canton Schwyz). The minimum age for the beginning of retreat is 22.2 ± 1.0 ka at the frontal (terminal) position and 20.4 ± 1.0 ka at the lateral position. These ages are directly comparable with exposure ages from the other piedmont lobes in the northern Alpine foreland. Our data from the mountain called Rigi, do not support the hypothesis that boulders located external to the Seeboden moraine were deposited prior to the last glacial cycle. We present a first exposure age from an erratic boulder in a retreat position in the Alpine foreland. The Reuss Glacier was approximately 12 km behind the maximal extent no later than at 18.6 ± 0.9 ka.

The first major incision of the Swiss Deckenschotter landscape

The Swiss Deckenschotter (“cover gravels”) is the oldest Quaternary units in the northern Swiss Alpine Foreland. They are a succession of glaciofluvial gravel layers intercalated with glacial and/or overbank deposits. This lithostratigraphic sequence is called Deckenschotter because it “covers” Molasse or Mesozoic bedrock and forms mesa-type hill-tops. Deckenschotter occurs both within and beyond the extent of the Last Glacial Maximum glaciers. The Swiss Deckenschotter consist of two sub-units: Höhere (Higher) and Tiefere (Lower) Deckenschotter. Although the Höhere Deckenschotter sub-unit (HDS) is topographically higher than the Tiefere Deckenschotter, it is older. The only available age for the Swiss Deckenschotter is 2.5–1.8 Ma based on mammal remains found in HDS at the Irchel site. In this study, we present an exposure age for the topographically lowest HDS, calculated from a cosmogenic 10Be depth-profile. Our results show that the first phase of the Deckenschotter glaciations in the Swiss Alps terminated at least 1,020+80−120 ka ago, which is indicated by a significant fluvial incision. This line of evidence seems to be close to synchronous with the beginning of the Mid-Pleistocene Revolution, when the frequency of the glacial-interglacial cyclicity changed from 41 to 100 ka and the amplitude from low to high, between marine isotope stages 23 and 22.

Minor inheritance inhibits the calibration of the 10Be production rate from the AD 1717 Val Ferret rock avalanche, European Alps

To calibrate the in situ 10Be production rate, we collected surface samples from nine large granitic boulders within the deposits of a rock avalanche that occurred in AD 1717 in the upper Ferret Valley, Mont Blanc Massif, Italy. The 10Be concentrations were extremely low and successfully measured within 10% analytical uncertainty or less. The concentrations vary from 4829 ± 448 to 5917 ± 476 at g−1. Using the historical age exposure time, we calculated the local and sea level-high latitude (i.e. ≥60°) cosmogenic 10Be spallogenic production rates. Depending on the scaling schemes, these vary between 4.60 ± 0.38 and 5.26 ± 0.43 at g−1 a−1. Although they correlate well with global values, our production rates are clearly higher than those from more recent calibration sites. We conclude that our 10Be production rate is a mean and an upper bound for production rates in the Massif region over the past 300 years. This rate is probably influenced by inheritance and will yield inaccurate (e.g. too young) exposure ages when applied to surface-exposure studies in the area. Other independently dated rock-avalanche deposits in the region that are approximately 103 years old could be considered as possible calibration sites.

36Cl production rate from K-spallation in the European Alps (Chironico landslide, Switzerland)

The abundant production of in situ cosmogenic 36Cl from potassium renders 36Cl measurements in K-rich rocks or minerals, such as K-feldspars, potentially useful for precisely dating rock surfaces, either in single-nuclide or in multi-nuclide studies, for example combined with 10Be measurements in quartz. However, significant discrepancies in experimentally calibrated 36Cl production rates from spallation of potassium (36PK-sp), referenced to sea-level/high-latitude (SLHL), limit the accuracy of 36Cl dating from K-rich lithologies. We present a new 36Cl calibration using K-feldspars, in which K-spallation is the most dominant 36Cl production pathway (>92% of total 36Cl), thus minimizing uncertainties from the complex multi-pathway 36Cl production systematics. The samples are derived from boulders of an ∼13.4 ka-old landslide in the Swiss Alps (∼820 m, 46.43°N, 8.85°E). We obtain a local 36PK-sp of 306 ± 16 atoms 36Cl (g K)−1 a−1 and an SLHL 36PK-sp of 145.5 ± 7.7 atoms 36Cl (g K)−1 a−1, when scaled with a standard scaling protocol (‘Lm’). Applying this SLHL 36PK-sp to determine 36Cl exposure ages of K-feldspars from 10Be-dated moraine boulders yields excellent agreement, confirming the validity of the new SLHL 36PK-sp for surface exposure studies, involving 36Cl in K-feldspars, in the Alps.

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.

Crete fault scarp point cloud data

This paper assesses the along strike variation of active bedrock fault scarps using long range terrestrial laser scanning (t-LiDAR) data in order to determine the distribution behaviour of scarp height and the subsequently calculate long term throw-rates. Five faults on Cretewhich display spectacular limestone fault scarps have been studied using high resolution digital elevation model (HRDEM) data. We scanned several hundred square metres of the fault system including the footwall, fault scarp and hanging wall of the investigated fault segment. The vertical displacement and the dip of the scarp were extracted every metre along the strike of the detected fault segment based on the processed HRDEM. The scarp variability was analysed by using statistical and morphological methods. The analysis was done in a geographical information system (GIS) environment. Results show a normal distribution for the scanned fault scarp's vertical displacement. Based on these facts, the mean value of height was chosen to define the authentic vertical displacement. Consequently the scarp can be divided into above, below and within the range of mean (within one standard deviation) and quantify the modifications of vertical displacement. Therefore, the fault segment can be subdivided into areas which are influenced by external modification like erosion and sedimentation processes. Moreover, to describe and measure the variability of vertical displacement along strike the fault, the semi-variance was calculated with the variogram method. This method is used to determine how much influence the external processes have had on the vertical displacement. By combining of morphological and statistical results, the fault can be subdivided into areas with high external influences and areas with authentic fault scarps, which have little or no external influences. This subdivision is necessary for long term throw-rate calculations, because without this differentiation the calculated rates would be misleading and the activity of a fault would be incorrectly assessed with significant implications for seismic hazard assessment since fault slip rate data govern the earthquake recurrence. Furthermore, by using this workflow areas with minimal external influences can be determined, not only for throw-rate calculations, but also for determining samples sites for absolute dating techniques such as cosmogenic nuclide dating. The main outcomes of this study include: i) there is no direct correlation between the fault's mean vertical displacement and dip (R² less than 0.31); ii) without subdividing the scanned scarp into areas with differing amounts of external influences, the along strike variability of vertical displacement is ±35%; iii) when the scanned scarp is subdivided the variation of the vertical displacement of the authentic scarp (exposed by earthquakes only) is in a range of ±6% (the varies depending on the fault from 7 to 12%); iv) the calculation of the long term throw-rate (since 13 ka) for four scarps in Crete using the authentic vertical displacement is 0.35 ± 0.04 mm/yr at Kastelli 1, 0.31 ± 0.01 mm/yr at Kastelli 2, 0.85 ± 0.06 mm/yr at the Asomatos fault (Sellia) and 0.55 ± 0.05 mm/yr at the Lastros fault.