Impact of Terrain Heterogeneity on Coherent Structure Properties: Numerical Approach

A three-dimensional large-eddy simulation (LES) model, which includes the effects of plant–atmosphere interactions, is used to study the effects of surface inhomogeneities on near-surface coherent structures over an open field and behind a forest canopy. These simulated conditions are representative of two wind sectors of the Site Instrumental de Recherche par Télédétection Atmosphérique (SIRTA) experimental site at the Institut Pierre Simon Laplace, Palaiseau, France. Coherent structure properties deduced from wavelet transforms of the simulated near-surface vertical velocity time series are not modified by upstream terrain heterogeneities, in agreement with site measurements. This feature is related to the nature of structures detected from the vertical velocity time series. The turbulence close to the surface seems composed of both local coherent structures and large coherent structures reflecting outer-layer properties, which depend on the overall surface heterogeneity or upstream heterogeneity. It is argued that the streamwise velocity is representative of these large outer-layer structures that impinge onto the ground through a top-down mechanism as identified through the space–time correlation of the wind velocity components. In contrast, the vertical velocity is more representative of small structures resulting from the impingement of the large outer-layer structures. These small structures represent locally-generated, active turbulence, which adjusts rapidly to local surface conditions, and consequently they are only weakly dependent on upstream heterogeneities.     

Mapping the mean monthly precipitation of a small island using kriging with external drifts

This study focuses on the spatial distribution of mean annual and monthly precipitation in a small island (1128 km²) named Martinique, located in the Lesser Antilles. Only 35 meteorological stations are available on the territory, which has a complex topography. With a digital elevation model (DEM), 17 covariates that are likely to explain precipitation were built. Several interpolation methods, such as regression-kriging (????????, ????????, and ????????) and external drift kriging (??????) were tested using a cross-validation procedure. For the regression methods, predictors were chosen by established techniques whereas a new approach is proposed to select external drifts in a kriging which is based on a stepwise model selection by the Akaike Information Criterion (AIC). The prediction accuracy was assessed at validation sites with three different skill scores. Results show that using methods with no predictors such as inverse distance weighting (??????) or universal kriging (????) is inappropriate in such a territory. ?????? appears to outperform regression methods for any criteria, and selecting predictors by our approach improves the prediction of mean annual precipitation compared to kriging with only elevation as drift. Finally, the predicting performance was also studied by varying the size of the training set leading to less conclusive results for ?????? and its performance. Nevertheless, the proposed method seems to be a good way to improve the mapping of climatic variables in a small island.     

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.     

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.     

Inversion of interference hydraulic pumping tests in both homogeneous and fractal dual media

This work proposes a complete method for automatic inversion of data from hydraulic interference pumping tests based on both homogeneous and fractal dual-medium approaches. The aim is to seek a new alternative concept able to interpret field data, identify macroscopic hydraulic parameters and therefore enhance the understanding of flow in porous fractured reservoirs. Because of its much contrasted sensitivities to parameters, the dual-medium approach yields an ill-posed inverse problem that requires a specific optimization procedure including the calculation of analytical sensitivities and their possible re-scaling. Once these constraints are fulfilled, the inversion proves accurate, provides unambiguous and reliable results. In the fractal context inverting several drawdown curves from different locations at the same time reveals more accurate. Finally, hydraulic parameters drawn from inversion should be taken into account to improve in various situations the conditioning of up-scaled flow in fractured rocks.     

Towards Improved Remedial Outcomes in Categorical Aquifers with an Iterative Ensemble Smoother

Categorical parameter distributions consisting of geologic facies with distinct properties, for example, high-permeability channels embedded in a low-permeability matrix, are common at contaminated sites. At these sites, low-permeability facies store solute mass, acting as secondary sources to higher-permeability facies, sustaining concentrations for decades while increasing risk and cleanup costs. Parameter estimation is difficult in such systems because the discontinuities in the parameter space hinder the inverse problem. This paper presents a novel approach based on Traveling Pilot Points (TRIPS) and an iterative ensemble smoother (IES) to solve the categorical inverse problem. Groundwater flow and solute transport in a hypothetical aquifer with a categorical parameter distribution are simulated using MODFLOW 6. Heads and concentrations are recorded at multiple monitoring locations. IES is used to generate posterior ensembles assuming a TRIPS prior and an approximate multi-Gaussian prior. The ensembles are used to predict solute concentrations and mass into the future. The evaluation also includes an assessment of how the number of measurements and the choice of the geological prior determine the characteristics of the posterior ensemble and the resulting predictions. The results indicate that IES was able to efficiently sample the posterior distribution and showed that even with an approximate geological prior, a high degree of parameterization and history matching could lead to parameter ensembles that can be useful for making certain types of predictions (heads, concentrations). However, the approximate geological prior was insufficient for predicting mass. The analysis demonstrates how decision-makers can quantify uncertainty and make informed decisions with an ensemble-based approach.     

Cretaceous length of day perturbation by mantle avalanche

These last 10 years, numerical models of mantle convection have emphasized the role of the 670 km endothermic phase change in generating avalanches that trigger catastrophic mass transfers between upper and lower mantle. On the other hand, scientists have emphasized the concomitance of large-scale worldwide geophysical and tectonic events, which could find their deep thermal roots in the huge mass transfers induced by the avalanches. In particular, the paleontological records show two periods of length of day (l.o.d.) shortening between 420 and 360, and 200 and 80 Myr BP. This last event is synchronous with a strong true polar wander and a global warming of the upper mantle. In order to study the potential effects of the avalanche on the main component of the Earth’s rotation, the Liouville equation has been solved and the l.o.d. evolution has been calculated from the perturbations of the inertia tensor. The results show that the inertia tensor of the Earth’s is mainly sensitive to the global transfers through the 670 km discontinuity. The l.o.d. perturbations will be synchronous with the global thermal effects of the avalanche. These theoretical results allow proposing a self-consistent physical mechanism to explain periods of the Earth’s rotation acceleration. Within this context, the l.o.d. shortening during the Cenozoic and Cretaceous brings one more clue to the possible participation of a mantle avalanche in generating the concomitant large scale events which have occurred during this very particular period of the Earth’s history.