Bayesian-MCMC inversion of time-lapse geophysical data to characterize vadose zone hydraulic behaviour

Geophysical methods have the potential to provide valuable information on hydrological properties in the unsaturated zone. In particular, time-lapse geophysical data, when coupled with a hydrological model and inverted stochastically, may allow for the effective estimation of subsurface hydraulic parameters and their corresponding uncertainties. In this study, we use a Bayesian Markov-chain-Monte-Carlo (MCMC) inversion approach to investigate how much information regarding vadose zone hydraulic properties can be retrieved from time-lapse crosshole GPR data collected at the Arrenaes field site in Denmark during a forced infiltration experiment.

Inversion of multiple intersecting high-resolution crosshole GPR profiles for hydrological characterization at the Boise Hydrogeophysical Research Site

The integration of geophysical data into the subsurface characterization problem has been shown in many cases to significantly improve hydrological knowledge by providing information at spatial scales and locations that is unattainable using conventional hydrological measurement techniques. In particular, crosshole ground-penetrating radar (GPR) tomography has shown much promise in hydrology because of its ability to provide highly detailed images of subsurface radar wave velocity, which is strongly linked to soil water content. Here, we develop and demonstrate a procedure for inverting together multiple crosshole GPR data sets in order to characterize the spatial distribution of radar wave velocity below the water table at the Boise Hydrogeophysical Research Site (BHRS) near Boise, Idaho, USA. Specifically, we jointly invert 31 intersecting crosshole GPR profiles to obtain a highly resolved and consistent radar velocity model along the various profile directions. The model is found to be strongly correlated with complementary neutron porosity-log data and is further corroborated by larger-scale structural information at the BHRS. This work is an important prerequisite to using crosshole GPR data together with existing hydrological measurements for improved groundwater flow and contaminant transport modeling.

Data-driven topo-climatic mapping with machine learning methods

Automatic environmental monitoring networks enforced by wireless communication technologies provide large and ever increasing volumes of data nowadays. The use of this information in natural hazard research is an important issue. Particularly useful for risk assessment and decision making are the spatial maps of hazard-related parameters produced from point observations and available auxiliary information. The purpose of this article is to present and explore the appropriate tools to process large amounts of available data and produce predictions at fine spatial scales. These are the algorithms of machine learning, which are aimed at non-parametric robust modelling of non-linear dependencies from empirical data. The computational efficiency of the data-driven methods allows producing the prediction maps in real time which makes them superior to physical models for the operational use in risk assessment and mitigation. Particularly, this situation encounters in spatial prediction of climatic variables (topo-climatic mapping). In complex topographies of the mountainous regions, the meteorological processes are highly influenced by the relief. The article shows how these relations, possibly regionalized and non-linear, can be modelled from data using the information from digital elevation models. The particular illustration of the developed methodology concerns the mapping of temperatures (including the situations of Föhn and temperature inversion) given the measurements taken from the Swiss meteorological monitoring network. The range of the methods used in the study includes data-driven feature selection, support vector algorithms and artificial neural networks.

Short- and long-term changes in chromosomal inversion polymorphism and global warming: Drosophila subobscura from the Balkans

The chromosomal inversion polymorphism of Drosophila subobscura is adaptive to environmental changes. The population of Petnica, Serbia, was chosen to analyze short- and long-term changes in this polymorphism. Short-term changes were studied in the samples collected in May, June, and August of 1995. The inversion polymorphism varied over these months, although various interpretations are possible. To analyze long-term changes, samples obtained in May 1995 and May 2010 were compared. The frequency of the 'cold' adapted inversions (Ast, Jst, Ust, Est, and Ost) decreased and that of the 'warm' adapted inversions (A2, J1, U1+2, and O3+4) increased, from 1995 to 2010. These changes are consistent with the general increase in temperature recorded in Petnica for the same period. Finally, the possible response of chromosomal polymorphism to global warming was analyzed at the regional level (Balkan peninsula). This polymorphism depends on the ecological conditions of the populations, and the changes observed appear to be consistent with global warming expectations. Natural selection seems to be the main mechanism responsible for the evolution of this chromosomal polymorphism.

Short- and long-term changes in chromosomal inversion polymorphism and global warming: Drosophila subobscura from the Balkans

The chromosomal inversion polymorphism of Drosophila subobscura is adaptive to environmental changes. The population of Petnica, Serbia, was chosen to analyze short- and long-term changes in this polymorphism. Short-term changes were studied in the samples collected in May, June, and August of 1995. The inversion polymorphism varied over these months, although various interpretations are possible. To analyze long-term changes, samples obtained in May 1995 and May 2010 were compared. The frequency of the 'cold' adapted inversions (Ast, Jst, Ust, Est, and Ost) decreased and that of the 'warm' adapted inversions (A2, J1, U1+2, and O3+4) increased, from 1995 to 2010. These changes are consistent with the general increase in temperature recorded in Petnica for the same period. Finally, the possible response of chromosomal polymorphism to global warming was analyzed at the regional level (Balkan peninsula). This polymorphism depends on the ecological conditions of the populations, and the changes observed appear to be consistent with global warming expectations. Natural selection seems to be the main mechanism responsible for the evolution of this chromosomal polymorphism.

3D coronary vessel wall imaging utilizing a local inversion technique with spiral image acquisition.

Current 2D black blood coronary vessel wall imaging suffers from a relatively limited coverage of the coronary artery tree. Hence, a 3D approach facilitating more extensive coverage would be desirable. The straightforward combination of a 3D-acquisition technique together with a dual inversion prepulse can decrease the effectiveness of the black blood preparation. To minimize artifacts from insufficiently suppressed blood signal of the nearby blood pools, and to reduce residual respiratory motion artifacts from the chest wall, a novel local inversion technique was implemented. The combination of a nonselective inversion prepulse with a 2D selective local inversion prepulse allowed for suppression of unwanted signal outside a user-defined region of interest. Among 10 subjects evaluated using a 3D-spiral readout, the local inversion pulse effectively suppressed signal from ventricular blood, myocardium, and chest wall tissue in all cases. The coronary vessel wall could be visualized within the entire imaging volume.

Parametric Approach to Blind Deconvolution of Nonlinear Channels

A parametric procedure for the blind inversion of nonlinear channels is proposed, based on a recent method of blind source separation in nonlinear mixtures. Experiments show that the proposed algorithms perform efficiently, even in the presence of hard distortion. The method, based on the minimization of the output mutual information, needs the knowledge of log-derivative of input distribution (the so-called score function). Each algorithm consists of three adaptive blocks: one devoted to adaptive estimation of the score function, and two other blocks estimating the inverses of the linear and nonlinear parts of the channel, (quasi-)optimally adapted using the estimated score functions. This paper is mainly concerned by the nonlinear part, for which we propose two parametric models, the first based on a polynomial model and the second on a neural network, while [14, 15] proposed non-parametric approaches.

Quasi-Nonparametric Blind Inversion of Wiener Systems

An e cient procedure for the blind inversion of a nonlinear Wiener system is proposed. We proved that the problem can be expressed as a problem of blind source separation in nonlinear mixtures, for which a solution has been recently proposed. Based on a quasi-nonparametric relative gradient descent, the proposed algorithm can perform e ciently even in the presence of hard distortions.

A Fast Gradient Approximation for Nonlinear Blind Signal Processing

When dealing with nonlinear blind processing algorithms (deconvolution or post-nonlinear source separation), complex mathematical estimations must be done giving as a result very slow algorithms. This is the case, for example, in speech processing, spike signals deconvolution or microarray data analysis. In this paper, we propose a simple method to reduce computational time for the inversion of Wiener systems or the separation of post-nonlinear mixtures, by using a linear approximation in a minimum mutual information algorithm. Simulation results demonstrate that linear spline interpolation is fast and accurate, obtaining very good results (similar to those obtained without approximation) while computational time is dramatically decreased. On the other hand, cubic spline interpolation also obtains similar good results, but due to its intrinsic complexity, the global algorithm is much more slow and hence not useful for our purpose.

A Canonical Genetic Algorithm for Blind Inversion of Linear Channels

It is well known the relationship between source separation and blind deconvolution: If a filtered version of an unknown i.i.d. signal is observed, temporal independence between samples can be used to retrieve the original signal, in the same manner as spatial independence is used for source separation. In this paper we propose the use of a Genetic Algorithm (GA) to blindly invert linear channels. The use of GA is justified in the case of small number of samples, where other gradient-like methods fails because of poor estimation of statistics.

Initialisation of Nonlinearities for PNL and Wiener systems Inversion

This paper proposes a very fast method for blindly initial- izing a nonlinear mapping which transforms a sum of random variables. The method provides a surprisingly good approximation even when the basic assumption is not fully satis¯ed. The method can been used success- fully for initializing nonlinearity in post-nonlinear mixtures or in Wiener system inversion, for improving algorithm speed and convergence.

Coherency and sharpness measures by using ICA algorithms. An investigation for Alzheimer’s disease discrimination

In this paper, we present a comprehensive study of different Independent Component Analysis (ICA) algorithms for the calculation of coherency and sharpness of electroencephalogram (EEG) signals, in order to investigate the possibility of early detection of Alzheimer’s disease (AD). We found that ICA algorithms can help in the artifact rejection and noise reduction, improving the discriminative property of features in high frequency bands (specially in high alpha and beta ranges). In addition to different ICA algorithms, the optimum number of selected components is investigated, in order to help decision processes for future works.

Investigation of ICA algorithms for feature extraction of EEG signals in discrimination of Alzheimer disease

In this paper we present a quantitative comparisons of different independent component analysis (ICA) algorithms in order to investigate their potential use in preprocessing (such as noise reduction and feature extraction) the electroencephalogram (EEG) data for early detection of Alzhemier disease (AD) or discrimination between AD (or mild cognitive impairment, MCI) and age-match control subjects.

Blind Inversion of Wiener Systems

A system in which a linear dynamic part is followed by a non linear memoryless distortion a Wiener system is blindly inverted This kind of systems can be modelised as a postnonlinear mixture and using some results about these mixtures an e cient algorithm is proposed Results in a hard situation are presented and illustrate the e ciency of this algorithm