Quantitative integration of high-resolution hydrogeophysical data through Monte-Carlo-type conditional simulations

Geophysical techniques can help to bridge the inherent gap with regard to spatial resolution and the range of coverage that plagues classical hydrological methods. This has lead to the emergence of the new and rapidly growing field of hydrogeophysics. Given the differing sensitivities of various geophysical techniques to hydrologically relevant parameters and their inherent trade-off between resolution and range the fundamental usefulness of multi-method hydrogeophysical surveys for reducing uncertainties in data analysis and interpretation is widely accepted. A major challenge arising from such endeavors is the quantitative integration of the resulting vast and diverse database in order to obtain a unified model of the probed subsurface region that is internally consistent with all available data. To address this problem, we have developed a strategy towards hydrogeophysical data integration based on Monte-Carlo-type conditional stochastic simulation that we consider to be particularly suitable for local-scale studies characterized by high-resolution and high-quality datasets. Monte-Carlo-based optimization techniques are flexible and versatile, allow for accounting for a wide variety of data and constraints of differing resolution and hardness and thus have the potential of providing, in a geostatistical sense, highly detailed and realistic models of the pertinent target parameter distributions. Compared to more conventional approaches of this kind, our approach provides significant advancements in the way that the larger-scale deterministic information resolved by the hydrogeophysical data can be accounted for, which represents an inherently problematic, and as of yet unresolved, aspect of Monte-Carlo-type conditional simulation techniques. We present the results of applying our algorithm to the integration of porosity log and tomographic crosshole georadar data to generate stochastic realizations of the local-scale porosity structure. Our procedure is first tested on pertinent synthetic data and then applied to corresponding field data collected at the Boise Hydrogeophysical Research Site near Boise, Idaho, USA.

Exposure to radio-frequency electromagnetic fields from broadcast transmitters and risk of childhood cancer: a census-based cohort study.

We investigated the association between exposure to radio-frequency electromagnetic fields (RF-EMFs) from broadcast transmitters and childhood cancer. First, we conducted a time-to-event analysis including children under age 16 years living in Switzerland on December 5, 2000. Follow-up lasted until December 31, 2008. Second, all children living in Switzerland for some time between 1985 and 2008 were included in an incidence density cohort. RF-EMF exposure from broadcast transmitters was modeled. Based on 997 cancer cases, adjusted hazard ratios in the time-to-event analysis for the highest exposure category (>0.2 V/m) as compared with the reference category (<0.05 V/m) were 1.03 (95% confidence interval (CI): 0.74, 1.43) for all cancers, 0.55 (95% CI: 0.26, 1.19) for childhood leukemia, and 1.68 (95% CI: 0.98, 2.91) for childhood central nervous system (CNS) tumors. Results of the incidence density analysis, based on 4,246 cancer cases, were similar for all types of cancer and leukemia but did not indicate a CNS tumor risk (incidence rate ratio = 1.03, 95% CI: 0.73, 1.46). This large census-based cohort study did not suggest an association between predicted RF-EMF exposure from broadcasting and childhood leukemia. Results for CNS tumors were less consistent, but the most comprehensive analysis did not suggest an association.

Simulated-annealing-based conditional simulation for the local-scale characterization of heterogeneous aquifers

Simulated-annealing-based conditional simulations provide a flexible means of quantitatively integrating diverse types of subsurface data. Although such techniques are being increasingly used in hydrocarbon reservoir characterization studies, their potential in environmental, engineering and hydrological investigations is still largely unexploited. Here, we introduce a novel simulated annealing (SA) algorithm geared towards the integration of high-resolution geophysical and hydrological data which, compared to more conventional approaches, provides significant advancements in the way that large-scale structural information in the geophysical data is accounted for. Model perturbations in the annealing procedure are made by drawing from a probability distribution for the target parameter conditioned to the geophysical data. This is the only place where geophysical information is utilized in our algorithm, which is in marked contrast to other approaches where model perturbations are made through the swapping of values in the simulation grid and agreement with soft data is enforced through a correlation coefficient constraint. Another major feature of our algorithm is the way in which available geostatistical information is utilized. Instead of constraining realizations to match a parametric target covariance model over a wide range of spatial lags, we constrain the realizations only at smaller lags where the available geophysical data cannot provide enough information. Thus we allow the larger-scale subsurface features resolved by the geophysical data to have much more due control on the output realizations. Further, since the only component of the SA objective function required in our approach is a covariance constraint at small lags, our method has improved convergence and computational efficiency over more traditional methods. Here, we present the results of applying our algorithm to the integration of porosity log and tomographic crosshole georadar data to generate stochastic realizations of the local-scale porosity structure. Our procedure is first tested on a synthetic data set, and then applied to data collected at the Boise Hydrogeophysical Research Site.

Active deformation fields: dense deformation field estimation for atlas-based segmentation using the active contour framework.

This paper presents a new and original variational framework for atlas-based segmentation. The proposed framework integrates both the active contour framework, and the dense deformation fields of optical flow framework. This framework is quite general and encompasses many of the state-of-the-art atlas-based segmentation methods. It also allows to perform the registration of atlas and target images based on only selected structures of interest. The versatility and potentiality of the proposed framework are demonstrated by presenting three diverse applications: In the first application, we show how the proposed framework can be used to simulate the growth of inconsistent structures like a tumor in an atlas. In the second application, we estimate the position of nonvisible brain structures based on the surrounding structures and validate the results by comparing with other methods. In the final application, we present the segmentation of lymph nodes in the Head and Neck CT images, and demonstrate how multiple registration forces can be used in this framework in an hierarchical manner.

Minimization of Nyquist ghosting for echo-planar imaging at ultra-high fields based on a "negative readout gradient" strategy.

PURPOSE: To improve the traditional Nyquist ghost correction approach in echo planar imaging (EPI) at high fields, via schemes based on the reversal of the EPI readout gradient polarity for every other volume throughout a functional magnetic resonance imaging (fMRI) acquisition train. MATERIALS AND METHODS: An EPI sequence in which the readout gradient was inverted every other volume was implemented on two ultrahigh-field systems. Phantom images and fMRI data were acquired to evaluate ghost intensities and the presence of false-positive blood oxygenation level-dependent (BOLD) signal with and without ghost correction. Three different algorithms for ghost correction of alternating readout EPI were compared. RESULTS: Irrespective of the chosen processing approach, ghosting was significantly reduced (up to 70% lower intensity) in both rat brain images acquired on a 9.4T animal scanner and human brain images acquired at 7T, resulting in a reduction of sources of false-positive activation in fMRI data. CONCLUSION: It is concluded that at high B(0) fields, substantial gains in Nyquist ghost correction of echo planar time series are possible by alternating the readout gradient every other volume.

A random linear functional approach to efficiency bounds

This paper suggests a method for obtaining efficiency bounds in models containing either only infinite-dimensional parameters or both finite- and infinite-dimensional parameters (semiparametric models). The method is based on a theory of random linear functionals applied to the gradient of the log-likelihood functional and is illustrated by computing the lower bound for Cox's regression model

In vivo Structural Imaging of the Cerebellum, the Contribution of Ultra-High Fields.

This review covers some of the contributions to date from cerebellar imaging studies performed at ultra-high magnetic fields. A short overview of the general advantages and drawbacks of the use of such high field systems for imaging is given. One of the biggest advantages of imaging at high magnetic fields is the improved spatial resolution, achievable thanks to the increased available signal-to-noise ratio. This high spatial resolution better matches the dimensions of the cerebellar substructures, allowing a better definition of such structures in the images. The implications of the use of high field systems is discussed for several imaging sequences and image contrast mechanisms. This review covers studies which were performed in vivo in both rodents and humans, with a special focus on studies that were directed towards the observation of the different cerebellar layers.

Self-avoiding random walks and Olbers' paradox

In this paper, we prove that a self-avoiding walk of infinite length provides a structure that would resolve Olbers' paradox. That is, if the stars of a universe were distributed like the vertices of an infinite random walk with each segment length of about a parsec, then the night sky could be as dark as actually observed on the Earth. Self-avoiding random walk structure can therefore resolve the Olbers' paradox even in a static universe.

Linking genetic and species diversities in a community of freshwater gastropods

Summary Biodiversity is usually studied through species or genetic diversities. To date, these two levels of diversity have remained the independent .fields of investigations of community ecologists and population geneticists. However, recent joint analyses of species and genetic diversities have suggested that common processes may underlie the two levels. Positive correlations between species diversity and genetic diversity may arise when the effects of drift and migration overwhelm selective effects. The first goal of this thesis was to make a joint investigation of the patterns of species and genetic diversity in a community of freshwater gastropods living in a floodplain habitat. The second goal was to determine, as far as possible, the relative influences of the processes underlying the patterns observed at each level. In chapter 2 we investigate the relative influences of different evolutionary forces in shaping the genetic structure of Radix balthica populations. Results revealed that the structure inferred using quantitative traits was lower or equal to the one inferred using neutral molecular markers. Consequently, the pattern of structure observed could be only due to random drift, possibly to uniform selection, but definitely not to selection for local optima. In chapter 3, we analyze the temporal variation of species and genetic diversities in five localities. An extended period of drought occurred at the end of the study period leading to decay of both species and genetic diversities. This parallel loss of diversity following a natural perturbation highlighted the role sometimes predominant of random drift over selection on patterns of biodiversity in a floodplain habitat. In chapter 4, we compare the spatial genetic structures of two sympatric species: Radix balthica and Planorbis carinatus. We found that R. balthica populations are weakly structured and have moderate to high values of gene diversity. In contrast, P. carinatus populations are highly structured and poorly diverse. Then we measured correlations between various indices of species and genetic diversity using genetic data .from the two species. We found only one significant correlation: between species richness and gene diversity of P. carinatus. This result highlights the .need to use genetic date from more than one species to infer correlations between species and genetic diversities. Overall, this thesis provided new insights into the common processes underlying patterns of species and genetic diversity. Résumé La biodiversité est généralement étudiée au niveau de la diversité génétique ou spécifique. Ces deux niveaux sont restés jusqu'à maintenant les domaines d'investigation séparés des généticiens des populations et des écologistes des communautés. Cependant, des analyses conjointes des diversités génétique et spécifique ont récemment suggéré que des processus similaires pouvaient influencer ces deux niveaux. Des corrélations positives entre les diversités génétique et spécifique pourraient être dues aux effets de migration et de dérive qui dominent les effets sélectifs. Le premier but de cette thèse était de faire une étude conjointe des diversités génétique et spécifique dans une communauté de gastéropodes d'eau douce. Le second objectif était de déterminer les influences relatives des différents processus liés à chaque niveau de diversité. Dans le chapitre 2 nous cherchons à déterminer quelles forces évolutives influencent la structure génétique de quatre populations de Radix balthica. La structure mesurée sur des traits quantitatifs s'est révélée être plus faible ou égale à celle mesurée avec des marqueurs moléculaires neutres. La structure observée pourrait ainsi être due uniquement à la dérive génétique, potentiellement à la sélection uniforme, mais en aucun cas à la sélection locale pour différents optima. Dans le chapitre 3 nous analysons la variation temporelle des diversités génétique et spécifique dans cinq localités. Une récente période de sécheresse a causé une diminution parallèle des deux niveaux de diversité. Cette perturbation à mis en évidence le rôle parfois prépondérant de la dérive par rapport à celui de la sélection dans le déterminisme de la biodiversité dans un écosytème alluvial. Dans le chapitre 4, nous comparons la structure génétique spatiale de deux espèces vivant en sympatrie : Radix balthica et Planorbis carinatus. Les populations de R. balthica sont peu structurées et présentent un niveau de diversité relativement élevé alors que celles de P. carinatus sont fortement structurées et peu diversifiées. Nous avons ensuite mesuré différentes corrélations entre les diversités génétique et spécifique, mais la seule relation significative a été trouvée entre la richesse spécifique et la diversité génétique de P. carinatus. Ainsi, cette thèse a permis de découvrir de nouveaux aspects des processus qui influencent en parallèle la diversité aux niveaux génétique et spécifique.

Scaling behavior of random knots.

Using numerical simulations we investigate how overall dimensions of random knots scale with their length. We demonstrate that when closed non-self-avoiding random trajectories are divided into groups consisting of individual knot types, then each such group shows the scaling exponent of approximately 0.588 that is typical for self-avoiding walks. However, when all generated knots are grouped together, their scaling exponent becomes equal to 0.5 (as in non-self-avoiding random walks). We explain here this apparent paradox. We introduce the notion of the equilibrium length of individual types of knots and show its correlation with the length of ideal geometric representations of knots. We also demonstrate that overall dimensions of random knots with a given chain length follow the same order as dimensions of ideal geometric representations of knots.

Geometric stopping of a random walk and its applications to valuing equity-linked death benefits

We study discrete-time models in which death benefits can depend on a stock price index, the logarithm of which is modeled as a random walk. Examples of such benefit payments include put and call options, barrier options, and lookback options. Because the distribution of the curtate-future-lifetime can be approximated by a linear combination of geometric distributions, it suffices to consider curtate-future-lifetimes with a geometric distribution. In binomial and trinomial tree models, closed-form expressions for the expectations of the discounted benefit payment are obtained for a series of options. They are based on results concerning geometric stopping of a random walk, in particular also on a version of the Wiener-Hopf factorization.

Priming psychic and conjuring abilities of a magic demonstration influences event interpretation and random number generation biases

Magical ideation and belief in the paranormal is considered to represent a trait-like character; people either believe in it or not. Yet, anecdotes indicate that exposure to an anomalous event can turn skeptics into believers. This transformation is likely to be accompanied by altered cognitive functioning such as impaired judgments of event likelihood. Here, we investigated whether the exposure to an anomalous event changes individuals' explicit traditional (religious) and non-traditional (e.g., paranormal) beliefs as well as cognitive biases that have previously been associated with non-traditional beliefs, e.g., repetition avoidance when producing random numbers in a mental dice task. In a classroom, 91 students saw a magic demonstration after their psychology lecture. Before the demonstration, half of the students were told that the performance was done respectively by a conjuror (magician group) or a psychic (psychic group). The instruction influenced participants' explanations of the anomalous event. Participants in the magician, as compared to the psychic group, were more likely to explain the event through conjuring abilities while the reverse was true for psychic abilities. Moreover, these explanations correlated positively with their prior traditional and non-traditional beliefs. Finally, we observed that the psychic group showed more repetition avoidance than the magician group, and this effect remained the same regardless of whether assessed before or after the magic demonstration. We conclude that pre-existing beliefs and contextual suggestions both influence people's interpretations of anomalous events and associated cognitive biases. Beliefs and associated cognitive biases are likely flexible well into adulthood and change with actual life events.