In situ monitoring of karst springs in Wallonia (southern Belgium)

The implementation of the Water Framework Directive, as well as the forthcoming entry into force of other European water policies focusing on water resources, require a thorough reorganization of groundwater resources monitoring in most European countries. In Wallonia, monitoring programs were initiated in 2005 to control the quantitative and chemical statuses of groundwater bodies, as well as the effectiveness of management plans aimed at achieving environmental objectives. In karst aquifers, springs are preferential monitoring targets even if, as a result of the high heterogeneity of this type of aquifer, the interpretation of time series and spatial data remains a challenge. Since 2006, a “springs” monitoring network has been progressively set up in the southern part of Belgium. Currently, in situ measurements of discharge, temperature, conductivity, turbidity and fluorescence of water are conducted at nine karstic outlets in the Devonian and the Carboniferous limestone aquifers of Wallonia. The main objective of this ongoing research is to evaluate the suitability and the robustness of such standard measurements to discriminate anthropogenic impacts from natural variations, at both short- and long-term scales. Discharge and temperature data revealed being useful in the assessment of regional hydro-climatic trends. Groundwater proved to be a better natural filter for the assessment of long-term climatic changes than mathematical filtering of noisy surface signals. Investigations of such data at a European scale could therefore give additional insight into the impacts of climate change on groundwater resources. Conductivity, turbidity and fluorescence data characterize the particulate, mineral and organic health of water and a better understanding of their natural dynamics could help in the early detection of anthropogenic deviations. However, the current reproducibility of these measurements is too low to ascertain observed trends and fluctuations, and further research is still needed.     

Dissolution–precipitation processes governing the carbonation and silicification of the serpentinite sole of the New Caledonia ophiolite

The weathering of mantle peridotite tectonically exposed to the atmosphere leads commonly to natural carbonation processes. Extensive cryptocrystalline magnesite veins and stock-work are widespread in the serpentinite sole of the New Caledonia ophiolite. Silica is systematically associated with magnesite. It is commonly admitted that Mg and Si are released during the laterization of overlying peridotites. Thus, the occurrence of these veins is generally attributed to a per descensum mechanism that involves the infiltration of meteoric waters enriched in dissolved atmospheric CO₂. In this study, we investigate serpentinite carbonation processes, and related silicification, based on a detailed petrographic and crystal chemical study of serpentinites. The relationships between serpentine and alteration products are described using an original method for the analysis of micro-X-ray fluorescence images performed at the centimeter scale. Our investigations highlight a carbonation mechanism, together with precipitation of amorphous silica and sepiolite, based on a dissolution–precipitation process. In contrast with the per descensum Mg/Si-enrichment model that is mainly concentrated in rock fractures, dissolution–precipitation process is much more pervasive. Thus, although the texture of rocks remains relatively preserved, this process extends more widely into the rock and may represent a major part of total carbonation of the ophiolite.     

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.     

Construction of three-phase data to model multiphase flow in porous media: Comparing an optimization approach to the finite element approach

Multiphase flow modelling is a major issue in the assessment of groundwater pollution. Three-phase flows are commonly governed by mathematical models that associate a pressure equation with two saturation equations. These equations involve a number of secondary variables that reflect the fluid behaviour in a porous medium. To improve the computational efficiency of multiphase flow simulators, several simplified reformulations of three-phase flow equations have been proposed. However, they require the construction of new secondary variables adapted to the reformulated flow equations. In this article, two different approaches are compared to quantify these variables. A numerical example is given for a typical fine sand.     

Adaptation of prestack migration to multi-offset ground-penetrating radar (GPR) data

In adapting the prestack migration technique used in seismic imaging to the inversion of ground‐penetrating radar (GPR) from time‐ to depth‐sections, we show that the theoretical integral formulation of the inversion can be applied to electromagnetic problems, albeit with three assumptions. The first two assumptions concern the electromagnetic characteristics of the medium, primarily that the medium must be perfectly resistive and non‐dispersive, and the third concerns the antennae radiation pattern, which is taken to be 2D. The application of this adaptation of the inversion method is confirmed by migrating actual GPR measurements acquired on the test site of the Laboratoire Central des Ponts et Chaussées. The results show good agreement with the geometry of the structures in the medium and confirm that the possible departure from the assumption of a purely resistive medium has no visible effect on the information concerning the geometry of scattering and reflecting structures. The field experiments also show that prestack migration processing is sufficiently robust with regard to the assumption of a non‐dispersive medium. The assumption of a 2D antennae radiation pattern, however, produces artefacts that could be significant for laterally heterogeneous media. Nevertheless, where the medium is not highly laterally heterogeneous, the migration gives a clear image of the scattering potential due to the geometry of structural contrasts in the medium; the scatterers are well focused from diffraction hyperbolae and well localized. Spatial geometry has limited dimensional accuracy and positions are located with a maximum error equal to the minimum wavelength of the signal bandpass. Objects smaller than one wavelength can nevertheless be detected and well focused if their dielectric contrasts are sufficiently high, as in the case of iron or water in gneiss gravels. Furthermore, the suitability of multi‐offset protocols to estimate the electromagnetic propagating velocity and to decrease the non‐coherent noise level of measurements is confirmed. Our velocity estimation is based on the semblance calculation of multi‐offset migrated images, and we confirmed the relevance of this quantification method using numerical data. The signal‐to‐noise ratio is improved by summing multi‐offset results after the addition of random noise on measurements. Thus the adaptation of prestack migration to multi‐offset radar measurements significantly improves the resolution of the scattering potential of the medium. Limitations associated with the methods used here suggest that 3D algorithms should be applied to strongly laterally heterogeneous media and further studies concerning the waveform inversion are necessary to obtain information about the electric nature of the medium.     

Early contamination of European flounder (Platichthys flesus) by PCDD/Fs and dioxin-like PCBs in European waters

Contamination levels and profiles of 7 polychlorinated-p-dioxins, 10 polychlorinated furans (PCDD/Fs) and 12 dioxin-like polychlorinated biphenyls (dl-PCBs) were investigated in juvenile European flounder (Platichthys flesus) captured in different nursery areas in the northeastern Atlantic coast across its geographical distribution range. The toxic equivalent concentrations (WHO-TEQ_fish) were also determined in order to evaluate which P. flesus population was more exposed to dioxin-like toxicity. Juveniles caught in the Sørfjord (Norway) showed the lowest WHO-TEQ_fish concentration (0.052 pg WHO-TEQ_fish g⁻¹ wet weight) whereas the highest value was observed in fish from the Wadden Sea (The Netherlands; 0.291 pg WHO-TEQ_fish g⁻¹ ww), mainly due to the greater contribution of 2,3,7,8-tetrachlorodibenzo-p-dioxin, the most toxic congener. Nonetheless, when comparing the results with existent tissue residue-based toxicity benchmarks, no adverse effects resulting from PCDD/Fs and dl-PCBs are expected to occur in flounder from the studied systems.     

Effects of magnetic interactions in anisotropy of magnetic susceptibility: Models, experiments and implications for igneous rock fabrics quantification

Besides granites of the ilmenite series, in which the anisotropy of magnetic susceptibility (AMS) is mainly controlled by paramagnetic minerals, the AMS of igneous rocks is commonly interpreted as the result of the shape-preferred orientation of unequant ferromagnetic grains. In a few instances, the anisotropy due to the distribution of ferromagnetic grains, irrespective of their shape, has also been proposed as an important AMS source. Former analytical models that consider infinite geometry of identical and uniformly magnetized and coaxial particles confirm that shape fabric may be overcome by dipolar contributions if neighboring grains are close enough to each other to magnetically interact. On these bases we present and experimentally validate a two-grain macroscopic numerical model in which each grain carries its own magnetic anisotropy, volume, orientation and location in space. Compared with analytical predictions and available experiments, our results allow to list and quantify the factors that affect the effects of magnetic interactions. In particular, we discuss the effects of (i) the infinite geometry used in the analytical models, (ii) the intrinsic shape anisotropy of the grains, (iii) the relative orientation in space of the grains, and (iv) the spatial distribution of grains with a particular focus on the inter-grain distance distribution. Using documented case studies, these findings are summarized and discussed in the framework of the generalized total AMS tensor recently introduced by Cañon-Tapia (Cañon-Tapia, E., 2001. Factors affecting the relative importance of shape and distribution anisotropy in rocks: theory and experiments. Tectonophysics, 340, 117–131.). The most important result of our work is that analytical models far overestimate the role of magnetic interaction in rock fabric quantification. Considering natural rocks as an assemblage of interacting and non-interacting grains, and that the effects of interaction are reduced by (i) the finite geometry of the interacting clusters, (ii) the relative orientation between interacting grains, (iii) their heterogeneity in orientation, shape and bulk susceptibility, and (iv) their inter-distance distribution, we reconcile analytical models and experiments with real case studies that minimize the role of magnetic interaction onto the measured AMS. Limitations of our results are discussed and guidelines are provided for the use of AMS in geological interpretation of igneous rock fabrics where magnetic interactions are likely to occur.     

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.     

A roadmap for a dedicated Earth Science Grid platform

Due to its intensive data processing and highly distributed organization, the multidisciplinary Earth Science applications community is uniquely positioned for the uptake and exploitation of Grid technologies. Currently Enabling Grids for E-sciencE, and other large Grid infrastructures are already deployed and capable of operational services. So far however, the adoption and exploitation of Grid technology throughout the Earth Science community has been slower than expected. The Dissemination and Exploitation of GRids in Earth sciencE project, proposed by the European Commission to assist and accelerate this process in a number of different ways, had between its main goals the creation of a roadmap towards Earth Science Grid platform. This paper presents the resulting roadmap.