Application of multiple-point geostatistics on modelling groundwater flow and transport in a cross-bedded aquifer (Belgium)

Sedimentological processes often result in complex three-dimensional subsurface heterogeneity of hydrogeological parameter values. Variogram-based stochastic approaches are often not able to describe heterogeneity in such complex geological environments. This work shows how multiple-point geostatistics can be applied in a realistic hydrogeological application to determine the impact of complex geological heterogeneity on groundwater flow and transport. The approach is applied to a real aquifer in Belgium that exhibits a complex sedimentary heterogeneity and anisotropy. A training image is constructed based on geological and hydrogeological field data. Multiple-point statistics are borrowed from this training image to simulate hydrofacies occurrence, while intrafacies permeability variability is simulated using conventional variogram-based geostatistical methods. The simulated hydraulic conductivity realizations are used as input to a groundwater flow and transport model to investigate the effect of small-scale sedimentary heterogeneity on contaminant plume migration. Results show that small-scale sedimentary heterogeneity has a significant effect on contaminant transport in the studied aquifer. The uncertainty on the spatial facies distribution and intrafacies hydraulic conductivity distribution results in a significant uncertainty on the calculated concentration distribution. Comparison with standard variogram-based techniques shows that multiple-point geostatistics allow better reproduction of irregularly shaped low-permeability clay drapes that influence solute transport.     

A multifractal formalism for vector-valued random fields based on wavelet analysis: application to turbulent velocity and vorticity 3D numerical data

Extreme atmospheric events are intimately related to the statistics of atmospheric turbulent velocities. These, in turn, exhibit multifractal scaling, which is determining the nature of the asymptotic behavior of velocities, and whose parameter evaluation is therefore of great interest currently. We combine singular value decomposition techniques and wavelet transform analysis to generalize the multifractal formalism to vector-valued random fields. The so-called Tensorial Wavelet Transform Modulus Maxima (TWTMM) method is calibrated on synthetic self-similar 2D vector-valued multifractal measures and monofractal 3D vector-valued fractional Brownian fields. We report the results of some application of the TWTMM method to turbulent velocity and vorticity fields generated by direct numerical simulations of the incompressible Navier–Stokes equations. This study reveals the existence of an intimate relationship between the singularity spectra of these two vector fields which are found significantly more intermittent than previously estimated from longitudinal and transverse velocity increment statistics.

---

Forecasting the energy produced by a windmill on a yearly basis

The objective of this article is to study as extensively as possible the uncertainties affecting the annual energy produced by a windmill. In the literature, the general approach is to estimate the mean annual energy from a transformation of a Weibull distribution law. Then the issue is reduced to estimating the coefficients of this distribution. This is obtained by classical statistical methods. Therefore, the uncertainties are mostly limited to those resulting from the statistical procedures. But in fact, the real uncertainty of the random variable which represents the annual energy cannot been reduced to the uncertainty on its mean and to the uncertainties induced from the estimation procedure. We propose here a model, which takes advantage of the fact that the annual energy production is the sum of many random variables representing the 10?min energy production during the year. Under some assumptions, we make use of the central limit theorem and show that an intrinsic uncertainties of wind power, usually not considered, carries an important risk. We also explain an observation coming from practice that the forecasted annual production is always overestimated, which creates a risk of reducing the profitability of the operation.     

Nutrient fluxes between water column and sediments: Potential influence of the pearl oyster culture

This study quantifies benthic nutrient fluxes and sedimentation rates in the Ahe Atoll lagoon (French Polynesia), in two stations located under pearl oyster frames, and two control stations away from the pearl culture facility. Dissolved inorganic nitrogen fluxes ranged between 2 and 35μmolNm(-2)h(-1) and Soluble Reactive Phosphorus varied between -3 and 8.2μmolPm(-2)h(-1). Particulate sedimentation rates beneath the oysters were approximately five times higher than in the control zone and the percentage of small particles (?63μm) were about the twice. In contrast, sediment composition was similar under and outside the direct influence of oyster frames. In this ecosystem, where primary production is dependent on the available nitrogen, our study revealed that, while highly variable, benthic fluxes could sometimes contribute up to 28% of the nitrogen demand for primary production.     

Use of a Sodar to Improve the Forecast of Fogs and Low Clouds on Airports

A sodar was deployed at Roissy–Charles de Gaulle airport near Paris, France, in 2008 with the aim of improving the forecast of low visibility conditions there. During the winter of 2008–2009, an experiment was conducted that showed that the sodar can effectively detect and locate the top of fog layers which is signaled by a strong peak of acoustic reflectivity. The peak is generated by turbulence activity in the inversion layer that contrasts sharply with the low reflectivity recorded in the fog layer below. A specific version of the 1D-forecast model deployed at Roissy for low visibility conditions (COBEL-ISBA) was developed in which fogs’ thicknesses are initialized by the sodar measurements rather than the information derived from the down-welling IR fluxes observed on the site. It was tested on data archived during the winters of 2008–2009 and 2009–2010 and compared to the version of the model presently operational. The results show a significant improvement—dissipation times of fogs are better predicted.     

Temporal variability of suspended sediment sources in an alpine catchment combining river/rainfall monitoring and sediment fingerprinting

Influence of the rainfall regime on erosion and transfer of suspended sediment in a 905‐km² mountainous catchment of the southern French Alps was investigated by combining sediment monitoring, rainfall data, and sediment fingerprinting (based on geochemistry and radionuclide concentrations). Suspended sediment yields were monitored between October 2007 and December 2009 in four subcatchments (22–713 km²). Automatic sediment sampling was triggered during floods to trace the sediment origin in the catchment. Sediment exports at the river catchment outlet (330 ± 100 t km^‐2 yr^‐1) were mainly driven (80%) by widespread rainfall events (long duration, low intensities). In contrast, heavy, local and short duration storms, generated high peak discharges and suspended sediment concentrations in small upstream torrents. However, these upstream floods had generally not the capacity to transfer the sediment down to the catchment outlet and the bulk of this fine sediment deposited along downstream sections of the river. This study also confirmed the important contribution of black marls (up to 70%) to sediment transported in rivers, although this substrate only occupies c. 10% of the total catchment surface. Sediment exports generated by local convective storms varied significantly at both intra‐ and inter‐flood scales, because of spatial heterogeneity of rainfall. However, black marls/marly limestones contribution remained systematically high. In contrast, widespread flood events that generate the bulk of annual sediment supply at the outlet were characterized by a more stable lithologic composition and by a larger contribution of limestones/marls, Quaternary deposits and conglomerates, which corroborates the results of a previous sediment fingerprinting study conducted on riverbed sediment. Copyright © 2012 John Wiley & Sons, Ltd.     

Monitoring Soil Water Content by Vertical Temperature Variations

The availability of high sensitivity temperature sensors (0.001 K sensitivity platinum resistors), which can be positioned at intervals of a few centimeters along a vertical profile in the unsaturated zone, allows short‐term in situ determinations—one day or even less—of the thermal diffusivity. The development of high data storage capabilities also makes this possible over long periods and the relative variations in thermal diffusivity allow the monitoring of the variations in water content. The processing of temperature measurements recorded at different depths is achieved by solving the heat equation, using the finite elements method, with both conductive and convective heat transfers. A first set of measurements has allowed this approach to be validated. Water content variations derived from thermal diffusivity values are in excellent agreement with TDR measurements carried out on the experimental site at Boissy‐le‐Châtel (Seine et Marne, France).     

Physically Based Groundwater Vulnerability Assessment Using Sensitivity Analysis Methods

A general physically based method is presented to assess the vulnerability of groundwater to external pressures by numerical simulation of groundwater flow. The concept of groundwater vulnerability assessment considered here is based on the calculation of sensitivity coefficients for a user‐defined groundwater state for which we propose several physically based indicators. Two sensitivity analysis methods are presented: the sensitivity equation method and the adjoint operator method. We show how careful selection of a method can significantly minimize the computational effort. An illustration of the general methodology is presented for the Herten aquifer analog (Germany). This application to a simple, yet insightful, case demonstrates the potential use of this general and physically based vulnerability assessment method to complex aquifers.