Alteration geochemistry of the Nopal I uranium deposit (Sierra Peña Blanca,Mexico), a natural analogue for a radioactive waste repository in volcanic tuffs

Natural analogues provide an approach to characterize and test the long‐term modelling of a repository performance. This article presents geochemical information about the alteration conditions of the Nopal I uranium deposit, Mexico, an analogue for the proposed Yucca Mountain radioactive waste repository. Mineralization and hydrothermal alteration of volcanic tuffs are contemporaneous, according to petrographic observations. Trace element geochemistry (U, Th, REE) provides evidence for local mobilization of uranium under oxidizing conditions and further precipitation under reducing conditions. O‐ and H‐isotope geochemistry of kaolinite, smectite, opal and calcite suggests that argillic alteration proceeded at shallow depth with meteoric water at 25–75 °C, a low‐temperature context, unusual for volcanic‐hosted uranium deposits. This temperature range is compatible with some post‐closure evolution models of the proposed Yucca Mountain repository.     

A Rhaetian 40Ar/39Ar age for the Rochechouart impact structure (France) and implications for the latest Triassic sedimentary record

Abstract– ⁴⁰Ar/³⁹Ar dating of potassium feldspar (primary spherulitic‐blocky and secondary idiomorphic K‐feldspar) separated from impact‐metamorphosed gneiss found near Videix in the western central part of the Rochechouart impact structure (NW Massif Central, France) yielded a Rhaetian combined age of 201 ± 2 Ma (2σ), indistinguishable within uncertainty from the age of the Triassic/Jurassic boundary. Ballen quartz intergrown with the primary K‐feldspar indicates post‐shock temperatures exceeding approximately 1000 °C that affected the precursor gneiss. Geochemically, both feldspar types represent essentially pure potassium end‐members. Apart from the approximately 15 km diameter impact deposit area, the youngest crystallization age known for basement rocks in this part of the Massif Central is approximately 300 Ma. No endogenic magmatic‐thermal events are known to have occurred later in this region. The K‐feldspar recrystallized from local feldspar melts and superimposed post‐shock hydrothermal crystallization, probably within some thousands of years after the impact. It is, therefore, suggested that the ⁴⁰Ar/³⁹Ar age for the Videix gneiss (as a potassic “impact metasomatite”) dates the Rochechouart impact, in consistence with evidence for K‐metasomatism in the Rochechouart impactites. The new age value is distinctly younger than the previously obtained Karnian–Norian age for Rochechouart and, thus, contradicts the Late Triassic multiple impact theory postulated some years ago. In agreement with the paleogeographic conditions in the western Tethys domain around the Triassic/Jurassic boundary, the near‐coastal to shallow marine Rochechouart impact is compatible with the formation of seismites and tsunami deposits in the latest Triassic of the British Isles and possible related deposits in other parts of Europe.     

Stochastic Hydrogeology: What Professionals Really Need?

Quantitative hydrogeology celebrated its 150th anniversary in 2006. Geostatistics is younger but has had a very large impact in hydrogeology. Today, geostatistics is used routinely to interpolate deterministically most of the parameters that are required to analyze a problem or make a quantitative analysis. In a small number of cases, geostatistics is combined with deterministic approaches to forecast uncertainty. At a more academic level, geostatistics is used extensively to study physical processes in heterogeneous aquifers. Yet, there is an important gap between the academic use and the routine applications of geostatistics. The reasons for this gap are diverse. These include aspects related to the hydrogeology consulting market, technical reasons such as the lack of widely available software, but also a number of misconceptions. A change in this situation requires acting at different levels. First, regulators must be convinced of the benefit of using geostatistics. Second, the economic potential of the approach must be emphasized to customers. Third, the relevance of the theories needs to be increased. Last, but not least, software, data sets, and computing infrastructure such as grid computing need to be widely available.     

Mobility of potentially harmful metals in latosols impacted by the municipal solid waste deposit of Londrina,Brazil

The contamination of soils by metals issuing from municipal solid waste (MSW) disposal in tropical environments has hardly been studied with regard to the particular problems associated with them, i.e., generally a high permeability of soils despite the abundance of clay, and the role of reactive Fe compounds. From a previous geotechnical and chemical survey, three latosol profiles differently affected by MSW leachates in the region of Londrina (Paraná, Brazil) were selected. The aims were to evaluate the extent of their contamination, to better understand the fate of potentially harmful metals in tropical soils and rank the determining factors. Samples between 0.5 and 7 m depth were analyzed for their physical, mineralogical and chemical properties, and their micro-morphology was described by optical and transmission electron microscopy. Two steps of a sequential extraction procedure helped to assess the mobility of elements and to better discriminate between metals originating from pedogenesis and issued from MSW. These combined approaches showed that exposed soil profiles have been impacted at various depths, down to 7 m, through increased metal content, especially enhanced mobility of Zn, Co, Mn, Cu and Fe, and through increased salinity and organic matter. The mobility of potentially harmful metals should decrease with pH, which significantly increased in some impacted horizons, but other factors can reverse this trend.     

Définition et datation de la Formation de Venaco (Corse) : dépôt d’origine gravitaire d’âge Priabonien

Priabonian age is highlighted for the first time in Corsica in the Venaco Formation using the presence of specific microfauna (in particular some representatives of Turborotalia cerroazulensis lineage). This silicoclastic formation is mainly represented by coarse facies. It is composed of three members from south to north and from oldest towards youngest: member of Uboli, Cardo and Orsu. The sedimentologic analysis reveals a gravity depositional environment, involving different type of currents. Sedimentologic and chronologic characteristics make the Venaco Formation and the Annot Formation (p.p.) equivalent. Dating the Venaco Fm. brings confirmation that the green schist metamorphism of the Variscan batholith and the related deformation are from the pre-Priabonian period.     

Géomorphologie de Tromelin,océan Indien

Tromelin is a small coral reef island (1 km²) located in the Indian Ocean, approximately 440 km east of Madagascar and 580 km north of La Réunion. Despite the presence of a permanent Météo France weather station on Tromelin since the 1950 s, a detailed geomorhological study of the island has never been undertaken. In this paper, we describe results from the 2008 field season, which enabled us to map and describe seven geomorphological zones on the atoll. (1) Tromelin's bathymetry is characterised by a pronounced submarine slope, which attains depths of 1000 m at only 2.5 km from the island's coastline. This geomorphology is typical of volcanic hotspots. (2) A coral platform borders the emerged crown of the island. This platform is exposed at low tide and has been fashioned in an ancient Pleistocene substratum. The substratum has been eroded by abrasive marine action to its present shape since the stabilisation of relative sea level around 6000 years ago. Erosion of the reef yields most of the biogenic sediment supply to the island. (3) The upper foreshore is characterised by outcrops of beachrock. These formations show a stratfication in slabs and dip gently towards the sea in the same manner as the contemporary clastic sediments. (4) Tromelin's beaches are divided into three sub-units, which reflect the different energy dynamics around the island, in addition to the transfer of sediment from the windward to the leeward side: (i) the southern storm beaches are characterised by coral blocks reworked during episodic high-energy events; they form boulder ramparts. Four perched rampart ridges have been recorded at ∼1 m, ∼3 m, ∼4.5 m and ∼7 m above mean sea level; (ii) transitional beaches are observed between the south and the north, constituting storm blocks and coarse sands and gravels; and (iii) sandy beaches are noted on the northern leeward side, comprising rounded gravels and coarse sands. The northern tip of the island is characterised by a sand spit (∼125 × ∼225 m), whose geometry varies based on the seasonal and meteo-marine contexts. (5) On Tromelin, a number of dune formations are observed. These are best developed in the north of the island due to the south-north translation of clastic sediments by aeolian wind action, and the trapping of sands by the well-developed vegetation. The height of the dunes varies between 10 cm (the southern micro-dunes) to 250 cm in the Northeast of Tromelin. (6) Exceptional waves have deposited storm tracts up to 250 m from the foreshore zone, comprising coarse gravels and coral blocks. The storm tract deposits are most prevalent on the windward side of the island, due not only to the high-energy dynamics, but also to the absence of vegetation cover in this area. (7) The southern depression covers an area of ∼40,000 m². At its lowest point it is only 1.2 m above sea level.     

Accuracy and efficiency of time integration methods for 1D diffusive wave equation

The diffusive wave approximation of the Saint-Venant equations is commonly used in hydrological models to describe surface flow processes. Numerous numerical approaches can be used to solve this highly nonlinear equation. Nonlinear time integration schemes—also called methods of lines (MOL)—were proven very efficient to solve other nonlinear problems in geosciences but were never considered to deal with surface flow modeling with the diffusive wave equation. In this paper, we study the relative performance of different time and space integration schemes by comparing the results obtained with classical approaches and with nonlinear time integration approaches. The results show that (i) the integration method with a higher order in space shows high accuracy regarding an integrated indicator such as the global mass balance error but is less accurate regarding local indicators, and (ii) nonlinear time integration techniques perform better than classical ones. Overall, it seems that integration techniques combining nonlinear time integration and a low spatial order need to be considered when developing hydrological modeling tools owing to their simplicity of implementation and very good performance.     

Peierls dislocation modelling in perovskite (CaTiO<Subscript>3</Subscript>): comparison with tausonite (SrTiO<Subscript>3</Subscript>) and MgSiO<Subscript>3</Subscript> perovskite

We present here a numerical modelling study of dislocations in perovskite CaTiO₃. The dislocation core structures and properties are calculated through the Peierls–Nabarro model using the generalized stacking fault (GSF) results as a starting model. The GSF are determined from first-principles calculations using the VASP code. The dislocation properties such as collinear, planar core spreading and Peierls stresses are determined for the following slip systems: [100](010), [100](001), [010](100), [010](001), [001](100), [001](010), and All dislocations exhibit lattice friction, but glide appears to be easier for [100](010) and [010](100). [001](010) and [001](100) exhibit collinear dissociation. Comparing Peierls stresses among tausonite (SrTiO₃), perovskite (CaTiO₃) and MgSiO₃ perovskite demonstrates the strong influence of orthorhombic distortions on lattice friction. However, and despite some quantitative differences, CaTiO₃ appears to be a satisfactory analogue material for MgSiO₃ perovskite as far as dislocation glide is concerned.     

Using Multiple-Point Geostatistics for Tracer Test Modeling in a Clay-Drape Environment with Spatially Variable Conductivity and Sorption Coefficient

This study investigates the effect of fine-scale clay drapes on tracer transport. A tracer test was performed in a sandbar deposit consisting of cross-bedded sandy units intercalated with many fine-scale clay drapes. The heterogeneous spatial distribution of the clay drapes causes a spatially variable hydraulic conductivity and sorption coefficient. A fluorescent tracer (sodium naphthionate) was injected in two injection wells and ground water was sampled and analyzed from five pumping wells. To determine (1) whether the fine-scale clay drapes have a significant effect on the measured concentrations and (2) whether application of multiple-point geostatistics can improve interpretation of tracer tests in media with complex geological heterogeneity, this tracer test is analyzed with a local three-dimensional ground-water flow and transport model in which fine-scale sedimentary heterogeneity is modeled using multiple-point geostatistics. To reduce memory needs and calculation time for the multiple-point geostatistical simulation step, this study uses the technique of direct multiple-point geostatistical simulation of edge properties. Instead of simulating pixel values, model cell edge properties indicating the presence of irregularly shaped surfaces are simulated using multiple-point geostatistical simulations. Results of a sensitivity analysis show under which conditions clay drapes have a significant effect on the concentration distribution. Calibration of the model against measured concentrations from the tracer tests reduces the uncertainty on the clay-drape parameters. The calibrated model shows which features of the breakthrough curves can be attributed to the geological heterogeneity of the aquifer and which features are caused by other processes.