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.     

Lander radioscience for obtaining the rotation and orientation of Mars

The paper presents the concept, the objectives, the approach used, and the expected performances and accuracies of a radioscience experiment based on a radio link between the Earth and the surface of Mars. This experiment involves radioscience equipment installed on a lander at the surface of Mars. The experiment with the generic name lander radioscience (LaRa) consists of an X-band transponder that has been designed to obtain, over at least one Martian year, two-way Doppler measurements from the radio link between the ExoMars lander and the Earth (ExoMars is an ESA mission to Mars due to launch in 2013). These Doppler measurements will be used to obtain Mars’ orientation in space and rotation (precession and nutations, and length-of-day variations). More specifically, the relative position of the lander on the surface of Mars with respect to the Earth ground stations allows reconstructing Mars’ time varying orientation and rotation in space.Precession will be determined with an accuracy better by a factor of 4 (better than the 0.1% level) with respect to the present-day accuracy after only a few months at the Martian surface. This precession determination will, in turn, improve the determination of the moment of inertia of the whole planet (mantle plus core) and the radius of the core: for a specific interior composition or even for a range of possible compositions, the core radius is expected to be determined with a precision decreasing to a few tens of kilometers.A fairly precise measurement of variations in the orientation of Mars’ spin axis will enable, in addition to the determination of the moment of inertia of the core, an even better determination of the size of the core via the core resonance in the nutation amplitudes. When the core is liquid, the free core nutation (FCN) resonance induces a change in the nutation amplitudes, with respect to their values for a solid planet, at the percent level in the large semi-annual prograde nutation amplitude and even more (a few percent, a few tens of percent or more, depending on the FCN period) for the retrograde ter-annual nutation amplitude. The resonance amplification depends on the size, moment of inertia, and flattening of the core. For a large core, the amplification can be very large, ensuring the detection of the FCN, and determination of the core moment of inertia.The measurement of variations in Mars’ rotation also determines variations of the angular momentum due to seasonal mass transfer between the atmosphere and ice caps. Observations even for a short period of 180 days at the surface of Mars will decrease the uncertainty by a factor of two with respect to the present knowledge of these quantities (at the 10% level).The ultimate objectives of the proposed experiment are to obtain information on Mars’ interior and on the sublimation/condensation of CO₂ in Mars’ atmosphere. Improved knowledge of the interior will help us to better understand the formation and evolution of Mars. Improved knowledge of the CO₂ sublimation/condensation cycle will enable better understanding of the circulation and dynamics of Mars’ atmosphere.     

Towards a better understanding of the Oulmes hydrogeological system (Mid-Atlas,Morocco)

Located in the Mid-Atlas (Morocco), the Oulmes plateau is famous for its mineral water springs “Sidi Ali” and “Lalla Haya” commercialised by the company “Les Eaux minérales d’Oulmès S.A”. Additionally, groundwater of the Oulmes plateau is intensively exploited for irrigation. The objective of this study, essentially performed from data collected during isotopic (summer 2004) and piezometric and hydrogeochemical field campaigns (spring 2007), is to improve the understanding of the Oulmes hydrogeological system. Analyses and interpretation of these data lead to the statement that this system is constituted by a main deep aquifer of large extension and by minor aquifers in a perched position. However, these aquifers interact enough to be in total equilibrium during the cold and wet period. As highlighted by isotopes, the origin of groundwater is mainly infiltration water except a small part of old groundwater with dissolved gas rising up from the granite through the schists.     

Identification of groundwater quality trends in a chalk aquifer threatened by intensive agriculture in Belgium

The European Union (EU) has adopted directives requiring that Member States take measures to reach a “good” chemical status of water resources by the year 2015 (Water Framework Directive: WFD). In order to achieve the environmental objectives for groundwater, the identification and reversal of significant upward trends in pollutant concentrations are required. A very detailed dataset available for the Hesbaye chalk aquifer in Belgium is used to evaluate tools and to propose efficient methodologies for identifying and quantifying nitrate trends in groundwater. Results indicate that the parametric linear regression and the non-parametric Mann-Kendall tests are robust; however, the latter test seems more adequate as it does not require verification of the normality of the dataset and it provides calculated nitrate trends very comparable to those obtained using linear regression. From a hydrogeological point of view, results highlight a general upward trend in the whole groundwater basin. The extrapolation of the trend analysis results indicates that measures have to be taken urgently in order to avoid further major degradation of groundwater quality within the next 10–70 years. However, a good groundwater quality status cannot be expected in the Hesbaye aquifer for the 2015 EU WFD deadline.     

Underground pumped storage hydroelectricity using abandoned works (deep mines or open pits) and the impact on groundwater flow

Underground pumped storage hydroelectricity (UPSH) plants using open-pit or deep mines can be used in flat regions to store the excess of electricity produced during low-demand energy periods. It is essential to consider the interaction between UPSH plants and the surrounding geological media. There has been little work on the assessment of associated groundwater flow impacts. The impacts on groundwater flow are determined numerically using a simplified numerical model which is assumed to be representative of open-pit and deep mines. The main impact consists of oscillation of the piezometric head, and its magnitude depends on the characteristics of the aquifer/geological medium, the mine and the pumping and injection intervals. If an average piezometric head is considered, it drops at early times after the start of the UPSH plant activity and then recovers progressively. The most favorable hydrogeological conditions to minimize impacts are evaluated by comparing several scenarios. The impact magnitude will be lower in geological media with low hydraulic diffusivity; however, the parameter that plays the more important role is the volume of water stored in the mine. Its variation modifies considerably the groundwater flow impacts. Finally, the problem is studied analytically and some solutions are proposed to approximate the impacts, allowing a quick screening of favorable locations for future UPSH plants.     

Innovative Contaminant Mass Flux Monitoring in an Aquifer Subject to Tidal Effects

Exposure from groundwater contamination to aquatic receptors residing in receiving surface water is dependent upon the rate of contaminated groundwater discharge. Characterization of groundwater fluxes is challenging, especially in coastal environments where tidal fluctuations result in transient groundwater flows towards these receptors. This can also be further complicated by the high spatial heterogeneity of subsurface deposits enhanced by anthropogenic influences such as the mixing of natural sediments and backfill materials, the presence of subsurface built structures such as sheet pile walls or even occurrence of other sources of contaminant discharge. In this study, the finite volume point dilution method (FVPDM) was successfully used to characterize highly transient groundwater flows and contaminant mass fluxes within a coastal groundwater flow system influenced by marked tides. FVPDM tests were undertaken continuously for more than 48 h at six groundwater monitoring wells, in order to evaluate groundwater flow dynamics during several tide cycles. Contaminant concentrations were measured simultaneously which allowed calculating contaminant mass fluxes. The study highlighted the importance of the aquifer heterogeneity, with groundwater fluxes ranging from 10⁻⁷ to 10⁻³ m/s. Groundwater flux monitoring enabled a significant refinement of the conceptual site model, including the fact that inversion of groundwater fluxes was not observed at high tide. Results indicated that contaminant mass fluxes were particularly higher at a specific monitoring well, by more than three orders of magnitude, than at other wells of the investigated aquifer. This study provided crucial information for optimizing further field investigations and risk mitigation measures.     

An atlas of predicted exotic gravitational lenses

We present an investigation of the relationships between the radio properties of a giant radio galaxy MRC B0319−454 and the surrounding galaxy distribution with the aim of examining the influence of intergalactic gas and gravity associated with the large-scale structure on the evolution in the radio morphology. Our new radio continuum observations of the radio source, with high surface brightness sensitivity, images the asymmetries in the megaparsec-scale radio structure in total intensity and polarization. We compare these with the three-dimensional galaxy distribution derived from galaxy redshift surveys. Galaxy density gradients are observed along and perpendicular to the radio axis: the large-scale structure is consistent with a model wherein the galaxies trace the ambient intergalactic gas and the evolution of the radio structures are ram-pressure limited by this associated gas. Additionally, we have modelled the off-axis evolution of the south-west radio lobe as deflection of a buoyant jet backflow by a transverse gravitational field: the model is plausible if entrainment is small. The case study presented here is a demonstration that giant radio galaxies may be useful probes of the warm-hot intergalactic medium believed to be associated with moderately over dense galaxy distributions.     

Tracking the sources and fate of nitrate pollution by combining hydrochemical and isotopic data with a statistical approach

This study contributes to identifying and spatializing the different types of nitrate sources by combining hydrogeochemical and isotopic data with principal component analysis (PCA) and t-distributed stochastic neighbor embedding (t-SNE) multicriteria statistical methods. The methodology is applied to the strategic Mons Basin chalk aquifer (Belgium). The results are based on a whole dataset containing 72 water samples with analyses of the hydrogeochemical parameters (temperature, pH, electrical conductivity (EC), redox potential, dissolved O₂), alkalinity, total organic carbon (TOC), silica (SiO₂), major and minor ions (NO₃^–, NH⁴⁺, Ca²⁺, dissolved Fe and Mn, K⁺, Mg²⁺, Na⁺, Sr²⁺, Cl^–, F^–, SO₄^–_, B) and multiple stable isotope ratios (δ¹¹B, δ¹⁵N–NO₃^–, δ¹⁸O–NO₃^–). Compared to classical PCA, the recently developed t-SNE method, which considers nonlinear relationships between variables and preserves local-scale similarities in a low-dimensional space, showed much better performance in discriminating different groups of samples and related zones in the aquifer. t-SNE results combined with isotope ratios highlighted four zones in the aquifer (grouped as A–D) and the presence of denitrification fronts. Group A presents a manure signature (δ¹⁵N–NO₃^– – mean (μ) +12.78‰, standard deviation (σ) 6.48‰; δ¹¹B – μ 29.96‰, σ 6.91‰). Group B exhibits both manure and inorganic fertilizer signatures (δ¹⁵N–NO₃^– – μ 6.27‰, σ 2.55‰; δ¹¹B – μ 15.86‰, σ 9.69‰). Group C shows a contamination by sewage (δ¹⁵N–NO₃^– – μ 12.67‰, σ 5.60‰; δ¹¹B – μ 9.97‰, σ 7.08‰). Group D presents a mixed signature (δ¹⁵N–NO₃^– – μ 9.25‰, σ 2.94‰; δ¹¹B – μ 20.00‰, σ 6.70‰).

     

The path towards high-contrast imaging with the VLTI: the Hi-5 project

The development of high-contrast capabilities has long been recognized as one of the top priorities for the VLTI. As of today, the VLTI routinely achieves contrasts of a few 10^??3 in the near-infrared with PIONIER (H band) and GRAVITY (K band). Nulling interferometers in the northern hemisphere and non-redundant aperture masking experiments have, however, demonstrated that contrasts of at least a few 10^??4 are within reach using specific beam combination and data acquisition techniques. In this paper, we explore the possibility to reach similar or higher contrasts on the VLTI. After reviewing the state-of-the-art in high-contrast infrared interferometry, we discuss key features that made the success of other high-contrast interferometric instruments (e.g., integrated optics, nulling, closure phase, and statistical data reduction) and address possible avenues to improve the contrast of the VLTI by at least one order of magnitude. In particular, we discuss the possibility to use integrated optics, proven in the near-infrared, in the thermal near-infrared (L and M bands, 3-5 \(\upmu \)m), a sweet spot to image and characterize young extra-solar planetary systems. Finally, we address the science cases of a high-contrast VLTI imaging instrument and focus particularly on exoplanet science (young exoplanets, planet formation, and exozodiacal disks), stellar physics (fundamental parameters and multiplicity), and extragalactic astrophysics (active galactic nuclei and fundamental constants). Synergies and scientific preparation for other potential future instruments such as the Planet Formation Imager are also briefly discussed. This project is called Hi-5 for High-contrast Interferometry up to 5 μm.     

Simulation of spatial and temporal trends in nitrate concentrations at the regional scale in the Upper Dyle basin,Belgium

Models are the only tools capable of predicting the evolution of groundwater systems at a regional scale, by taking into account a large amount of information. This study presents the association of a water balance model (WetSpass) with a groundwater flow and solute transport model (SUFT3D, saturated and unsaturated flow and transport in 3D) in order to simulate the present and future groundwater quality in terms of nitrate in the Upper Dyle basin (439 km²) Belgium. The HFEMC (hybrid finite element mixing cell) method implemented in the SUFT3D code is used to model groundwater flow and nitrate transport. Spatially distributed recharge, modelled with WetSpass, is considered for prescribing the recharge to the groundwater flow model. The feasibility of linking the WetSpass model with the finite-elements SUFT3D code is demonstrated. Time evolution and distribution of nitrate concentration are then simulated using the calibrated model. Nitrate inputs are spatially distributed according to land use. The spatial simulations and temporal trends are compared with previously published data on this aquifer and show good results.