Sampling and Characterizing Rare Earth Elements in Groundwater in Deep-Lying Fractures in Granitoids Under In Situ High-Pressure and Low-Redox Conditions

Several countries are preparing to dispose of radioactive nuclear waste deep underground in crystalline rock. This type of bedrock is commonly extensively fractured and consequently carries groundwater that serves as a medium for transporting metals and radionuclides. A group of metals of particular interest in this context is the rare earth elements (REEs), because they are analogues of actinides contained within radioactive waste and are tracers of hydrological pathways and geochemical processes. Concentrations of REEs are commonly low in these groundwaters, leading to values below detection limits of standard monitoring methods, particularly for the heavy REEs. We present a new technical set-up for monitoring REEs (and other trace metals) in groundwater in fractured crystalline rock. The technique consists of passing the fracture groundwater, commonly under high pressure and containing reduced chemical species, through a device that maintains the physicochemical character of the groundwater. Within the device, diffusive gradient in thin-film (DGT) discs are installed in triplicate. With this set-up, we studied REEs in groundwater in fractures at depths of approximately ?144, ?280, and ?450 m in granitoids in the Äspö Hard Rock Laboratory in southern Sweden. The entire REE suite was detected (concentrations down to 0.1 ng L^?1) and was differently fractionated among the groundwaters. The shallowest groundwater, composed of dilute modern Baltic Sea water, was enriched in the heavy REEs, whereas the deeper groundwaters, dominated by old saline water, were depleted in the heavy REEs. Deployment periods varying from 1 to 4 weeks delivered similar REE concentrations, indicating stability and reproducibility of the experimental set-up. The study finds that 1 week of deployment may be enough. However, if the overall setting and construction allow for longer deployment times, 2–3 weeks will be optimal in terms of reaching reliable REE concentrations well above the detection limit while maintaining the performance of the DGT samplers.     

Recharge and sustainability of a coastal aquifer in northern Albania

The River Mati in Albania has formed a coastal plain with Holocene and Pleistocene sediments. The outer portion of the plain is clay, with three underlying aquifers that are connected to an alluvial fan at the entry of the river into the plain. The aquifers supply water for 240,000 people. Close to the sea the aquifers are brackish. The brackish water is often artesian and found to be thousands of years old. Furthermore, the salinity, supported by δ¹⁸O results, does not seem to be due to mixing with old seawater but due to diffusion from intercalated clay layers. Heavy metals from mines in the upstream section of River Mati are not an immediate threat, as the pH buffering of the river water is good. Moreover, the heavy metals are predominantly found in suspended and colloidal phases. Two sulphur isotope signatures, one mirroring seawater sulphate in the brackish groundwater (δ³⁴S >21 ‰) and one showing the influence of sulphide in the river and the fresh groundwater (δ³⁴S <10 ‰), indicate that the groundwater in the largest well field is recharged from the river. The most serious threat is gravel extraction in the alluvial fan, decreasing the hydraulic head necessary for recharge and causing clogging of sediments.     

New Methodology to Characterize Shear Behavior of Joints by Combination of Direct Shear Box Testing and Numerical Simulations

A new procedure is presented, which combines big shear box tests on rocks and corresponding numerical simulations with explicit consideration of joint roughness to get deeper insight into the shear behavior of rock joints. The procedure consists of three parts: (1) constant normal load- or CNS-shear box tests with registration of shear- and normal-components of stress and displacements and deduction of basis rock mechanical parameters; (2) high resolution 3D-scanning of joint surface to deduce joint topography; and (3) set-up, run and evaluation of 3-dimensional numerical model with explicit duplication of joint roughness as back-analysis of shear box tests. The numerical back-analysis provides deeper insight into the joint behavior at the micro-scale. Several parameters can be deduced, like micro-slope angle distribution, aperture size distribution, local normal stress distribution and detailed analysis of dilation in relation to shear direction. The potential of the new procedure is illustrated exemplary by shear box tests on slate.     

Evidence for delayed poleward expansion of North Atlantic surface waters during the last interglacial (MIS 5e)

The Last Interglacial (Marine Isotopic Stage or MIS 5e) surface ocean heat flux from the Rockall Basin (NE Atlantic) towards the Arctic Ocean was reconstructed by analysing dinoflagellate cyst (dinocyst) assemblages in four sediment cores. Together with records of stable isotopes and ice-rafted detritus, the assemblage data reflect the northward retreat of ice(berg)-laden waters and the gradual development towards interglacial conditions at the transition from the Saalian deglaciation (Termination II) into MIS 5e. At the Rockall Basin, this onset of the Last Interglacial is soon followed by the appearance of the thermophilic dinocyst species Spiniferites mirabilis, with relative abundances higher than those observed at present in the area. North of the Iceland-Scotland Ridge, however, S. mirabilis only appears in significant numbers during late MIS 5e, between ～118 and 116.5 ka. Hence, fully marine Last Interglacial conditions with most intense Atlantic surface water influence occurred during late MIS 5e in the Nordic seas, and consequently also farther north in the Arctic Ocean, and at times when northern hemisphere summer insolation was already significantly decreased. The stratigraphic position of this Late Interglacial optimum is supported by planktic foraminifers and contrasts with the timing of the early Holocene climatic optimum in this area. We interpret the delayed northward expansion of Atlantic waters towards the polar latitudes as a result of the Saalian ice sheet deglaciation and its specific impact on the subsequent water mass evolution in this region.     

A mass‐balance approach to estimate in‐stream processes in a large river

A mass‐balance approach was used to estimate in‐stream processes related to inorganic nitrogen species (NH₄⁺, NO₂^? and NO₃^?) in a large river characterized by highly variable hydrological conditions, the Garonne River (south‐west France). Studies were conducted in two consecutive reaches of 30 km located downstream of the Toulouse agglomeration (population 760 000, seventh order), impacted by modification of discharge regime and high nitrogen concentrations. The mass‐balance was calculated by two methods: the first is based on a variable residence time (VRT) simulated by a one‐dimensional (1‐D) hydraulic model; the second is a based on a calculation using constant residence time (CRT) evaluated according to hydrographic peaks. In the context of the study, removal of dissolved inorganic nitrogen (DIN) for a reach of 30 km is underestimated by 11% with the CRT method. In sub‐reaches, the discrepancy between the two methods led to a 50% overestimation of DIN removal in the upper reach (13 km) and a 43% underestimation in the lower reach (17 km) using the CRT method. The study highlights the importance of residence time determination when using modelling approaches in the assessment of whole stream processes in short‐duration mass‐balance for a large river under variable hydrological conditions. Copyright © 2007 John Wiley & Sons, Ltd.     

Orbit Determination of Binary Asteroids

In addition to the detection of an asteroid moon or a binary asteroid, the knowledge of the satellite’s true orbit is of high importance to derive fundamental physical parameters of the binary system such as its mass and to shed light on its possible formation history and dynamical evolution (prograde/retrograde orbit, large/small eccentricity or inclination, etc.). A new methodology for preliminary orbit determination of binary asteroids – and visual binaries in general – is proposed. It is based on Thiele–Innes method combined with a ‘trial and error’ Monte-Carlo technique. This method provides the full set of solutions (bundle of orbits, with the 7 orbital elements) even for a reduced number of observations. The mass is a direct by-product of this orbit determination, from which one can next infer the bulk-density and porosity. In addition to the bundle of orbits, the method provides the marginal probability densities of the foreseen parameters. Such error analysis – since it avoids linear approximation – can be of importance for the prediction of the satellite’s position in the plane-of-sky during future stellar occultations or subsequent observations, but also for the analysis of the orbit’s secular evolution. After briefly describing the method, we present the algorithm and its application to some practical cases, with particular emphasis on asteroids binaries and applications on orbital evolution.