Source wavelet estimation for waveform inversion of crosshole georadar data

AbstractFor a wide range of environmental, hydrological, and engineering applications there is a fast growing need for high-resolution imaging. In this context, waveform tomographic imaging of crosshole georadar data is a powerful method able to provide images of pertinent electrical properties in near-surface environments with unprecedented spatial resolution. In contrast, conventional ray-based tomographic methods, which consider only a very limited part of the recorded signal (first-arrival traveltimes and maximum first-cycle amplitudes), suffer from inherent limitations in resolution and may prove to be inadequate in complex environments. For a typical crosshole georadar survey the potential improvement in resolution when using waveform-based approaches instead of ray-based approaches is in the range of one order-of- magnitude. Moreover, the spatial resolution of waveform-based inversions is comparable to that of common logging methods. While in exploration seismology waveform tomographic imaging has become well established over the past two decades, it is comparably still underdeveloped in the georadar domain despite corresponding needs. Recently, different groups have presented finite-difference time-domain waveform inversion schemes for crosshole georadar data, which are adaptations and extensions of Tarantola's seminal nonlinear generalized least-squares approach developed for the seismic case. First applications of these new crosshole georadar waveform inversion schemes on synthetic and field data have shown promising results. However, there is little known about the limits and performance of such schemes in complex environments. To this end, the general motivation of my thesis is the evaluation of the robustness and limitations of waveform inversion algorithms for crosshole georadar data in order to apply such schemes to a wide range of real world problems.One crucial issue to making applicable and effective any waveform scheme to real-world crosshole georadar problems is the accurate estimation of the source wavelet, which is unknown in reality. Waveform inversion schemes for crosshole georadar data require forward simulations of the wavefield in order to iteratively solve the inverse problem. Therefore, accurate knowledge of the source wavelet is critically important for successful application of such schemes. Relatively small differences in the estimated source wavelet shape can lead to large differences in the resulting tomograms. In the first part of my thesis, I explore the viability and robustness of a relatively simple iterative deconvolution technique that incorporates the estimation of the source wavelet into the waveform inversion procedure rather than adding additional model parameters into the inversion problem. Extensive tests indicate that this source wavelet estimation technique is simple yet effective, and is able to provide remarkably accurate and robust estimates of the source wavelet in the presence of strong heterogeneity in both the dielectric permittivity and electrical conductivity as well as significant ambient noise in the recorded data. Furthermore, our tests also indicate that the approach is insensitive to the phase characteristics of the starting wavelet, which is not the case when directly incorporating the wavelet estimation into the inverse problem.Another critical issue with crosshole georadar waveform inversion schemes which clearly needs to be investigated is the consequence of the common assumption of frequency- independent electromagnetic constitutive parameters. This is crucial since in reality, these parameters are known to be frequency-dependent and complex and thus recorded georadar data may show significant dispersive behaviour. In particular, in the presence of water, there is a wide body of evidence showing that the dielectric permittivity can be significantly frequency dependent over the GPR frequency range, due to a variety of relaxation processes. The second part of my thesis is therefore dedicated to the evaluation of the reconstruction limits of a non-dispersive crosshole georadar waveform inversion scheme in the presence of varying degrees of dielectric dispersion. I show that the inversion algorithm, combined with the iterative deconvolution-based source wavelet estimation procedure that is partially able to account for the frequency-dependent effects through an "effective" wavelet, performs remarkably well in weakly to moderately dispersive environments and has the ability to provide adequate tomographic reconstructions.

Adaptation en français et en allemand d'une échelle de pression des pairs pour jeunes adultes : le Peer Pressure Inventory [Adaptation of a peer pressure scale in French and German: the Peer Pressure Inventory].

BACKGROUND: Peer pressure is regarded as an important determinant of substance use, sexual behavior and juvenile delinquency. However, few peer pressure scales are validated, especially in French or German. Little is known about the factor structure of such scales or the kind of scale needed: some scales takes into account both peer pressure to do and peer pressure not to do, while others consider only peer pressure to do. The aim of the present study was to adapt French and German versions of the Peer Pressure Inventory, which is one of the most widely used scales in this field. We considered its factor structure and concurrent validity. METHODS: Five thousand eight hundred and sixty-seven young Swiss men filled in a questionnaire on peer pressure, substance use, and other variables (conformity, involvement) in a cohort study. RESULTS: We identified a four-factor structure, with the three factors of the initial Peer Pressure Inventory (involvement, conformity, misconduct) and adding a new one (relationship with girls). A non-valued scale (from no peer pressure to peer pressure to do only) showed stronger psychometric qualities than a valued scale (from peer pressure not to do to peer pressure to do). Concurrent validity was also good. Each behavior or attitude was significantly associated with peer pressure. CONCLUSION: Peer pressure seems to be a multidimensional concept. In this study, peer pressure to do showed the strongest influence on participants. Indeed, peer pressure not to do did not add anything useful. Only peer pressure to do affected young Swiss men's behaviors and attitudes and was reliable.

A mixed approach for proving non-inferiority in clinical trials with binary endpoints

When a new treatment is compared to an established one in a randomized clinical trial, it is standard practice to statistically test for non-inferiority rather than for superiority. When the endpoint is binary, one usually compares two treatments using either an odds-ratio or a difference of proportions. In this paper, we propose a mixed approach which uses both concepts. One first defines the non-inferiority margin using an odds-ratio and one ultimately proves non-inferiority statistically using a difference of proportions. The mixed approach is shown to be more powerful than the conventional odds-ratio approach when the efficacy of the established treatment is known (with good precision) and high (e.g. with more than 56% of success). The gain of power achieved may lead in turn to a substantial reduction in the sample size needed to prove non-inferiority. The mixed approach can be generalized to ordinal endpoints.

Numerical simulation of ambient seismic wavefield modification caused by pore-fluid effects in an oil reservoir

We have modeled numerically the seismic response of a poroelastic inclusion with properties applicable to an oil reservoir that interacts with an ambient wavefield. The model includes wave-induced fluid flow caused by pressure differences between mesoscopic-scale (i.e., in the order of centimeters to meters) heterogeneities. We used a viscoelastic approximation on the macroscopic scale to implement the attenuation and dispersion resulting from this mesoscopic-scale theory in numerical simulations of wave propagation on the kilometer scale. This upscaling method includes finite-element modeling of wave-induced fluid flow to determine effective seismic properties of the poroelastic media, such as attenuation of P- and S-waves. The fitted, equivalent, viscoelastic behavior is implemented in finite-difference wave propagation simulations. With this two-stage process, we model numerically the quasi-poroelastic wave-propagation on the kilometer scale and study the impact of fluid properties and fluid saturation on the modeled seismic amplitudes. In particular, we addressed the question of whether poroelastic effects within an oil reservoir may be a plausible explanation for low-frequency ambient wavefield modifications observed at oil fields in recent years. Our results indicate that ambient wavefield modification is expected to occur for oil reservoirs exhibiting high attenuation. Whether or not such modifications can be detected in surface recordings, however, will depend on acquisition design and noise mitigation processing as well as site-specific conditions, such as the geologic complexity of the subsurface, the nature of the ambient wavefield, and the amount of surface noise.

Analysis of an iterative deconvolution approach for estimating the source wavelet during waveform inversion of crosshole ground-penetrating radar data

A major issue in the application of waveform inversion methods to crosshole georadar data is the accurate estimation of the source wavelet. Here, we explore the viability and robustness of incorporating this step into a time-domain waveform inversion procedure through an iterative deconvolution approach. Our results indicate that, at least in non-dispersive electrical environments, such an approach provides remarkably accurate and robust estimates of the source wavelet even in the presence of strong heterogeneity in both the dielectric permittivity and electrical conductivity. Our results also indicate that the proposed source wavelet estimation approach is relatively insensitive to ambient noise and to the phase characteristics of the starting wavelet. Finally, there appears to be little-to-no trade-off between the wavelet estimation and the tomographic imaging procedures.

Termination of major ductile strike-slip shear and differential cooling along the Insubric line (Central Alps): U-Pb, Rb-Sr and 40Ar/39Ar ages of cross-cutting pegmatites

To constrain the age of strike-slip shear, related granitic magmatism, and cooling along the Insubric line, 29 size fractions of monazite and xenotime were dated by the U-Pb method, and a series of 25 Rb-Sr and Ar-40/Ar-39 ages were measured on different size fractions of muscovite and biotite. The three pegmatitic intrusions analyzed truncate high-grade metamorphic mylonite gneisses of the Simplon shear zone, a major Alpine structure produced in association with dextral strike-slip movements along the southern edge of the European plate, after collision with its Adriatic indenter. Pegmatites and aplites were produced between 29 and 25 Ma in direct relation to right-lateral shear along the Insubric line, by melting of continental crust having Sr-87/Sr-86 between 0.7199 and 0.7244 at the time of melting. High-temperature dextral strike-slip shear was active at 29.2 +/- 0.2 (2 sigma) Ma, and it terminated before 26.4 +/- 0.1 Ma. During dike injection, temperatures in the country rocks of the Isorno-Orselina and Monte Rosa structural units did not exceed approximate to 500 degrees C, leading to fast initial cooling, followed by slower cooling to approximate to 350 degrees C within several million years. In one case, initial cooling to approximate to 500 degrees C was significantly delayed by about 4 m.y., with final cooling to approximate to 300 degrees C at 20-19 Ma in all units. For the period between 29 and 19 Ma, cooling of the three sample localities was non-uniform in space and time, with significant variations on the kilometre scale. These differences are most likely due to strongly varying heat flow, and/or heterogeneous distribution of unroofing rates within the continuously deforming Insubric line. If entirely ascribed to differences in unroofing, corresponding rates would vary between 0.5 and 2.5 mm/y, for a thermal gradient of 30 degrees/km.

From Regional Landslide Detection to Site-Specific Slope Deformation Monitoring and Modelling Based on Active Remote Sensors

Landslide processes can have direct and indirect consequences affecting human lives and activities. In order to improve landslide risk management procedures, this PhD thesis aims to investigate capabilities of active LiDAR and RaDAR sensors for landslides detection and characterization at regional scales, spatial risk assessment over large areas and slope instabilities monitoring and modelling at site-specific scales. At regional scales, we first demonstrated recent boat-based mobile LiDAR capabilities to model topography of the Normand coastal cliffs. By comparing annual acquisitions, we validated as well our approach to detect surface changes and thus map rock collapses, landslides and toe erosions affecting the shoreline at a county scale. Then, we applied a spaceborne InSAR approach to detect large slope instabilities in Argentina. Based on both phase and amplitude RaDAR signals, we extracted decisive information to detect, characterize and monitor two unknown extremely slow landslides, and to quantify water level variations of an involved close dam reservoir. Finally, advanced investigations on fragmental rockfall risk assessment were conducted along roads of the Val de Bagnes, by improving approaches of the Slope Angle Distribution and the FlowR software. Therefore, both rock-mass-failure susceptibilities and relative frequencies of block propagations were assessed and rockfall hazard and risk maps could be established at the valley scale. At slope-specific scales, in the Swiss Alps, we first integrated ground-based InSAR and terrestrial LiDAR acquisitions to map, monitor and model the Perraire rock slope deformation. By interpreting both methods individually and originally integrated as well, we therefore delimited the rockslide borders, computed volumes and highlighted non-uniform translational displacements along a wedge failure surface. Finally, we studied specific requirements and practical issues experimented on early warning systems of some of the most studied landslides worldwide. As a result, we highlighted valuable key recommendations to design new reliable systems; in addition, we also underlined conceptual issues that must be solved to improve current procedures. To sum up, the diversity of experimented situations brought an extensive experience that revealed the potential and limitations of both methods and highlighted as well the necessity of their complementary and integrated uses.

Mate and fuse: how yeast cells do it.

Many cells are able to orient themselves in a non-uniform environment by responding to localized cues. This leads to a polarized cellular response, where the cell can either grow or move towards the cue source. Fungal haploid cells secrete pheromones to signal mating, and respond by growing a mating projection towards a potential mate. Upon contact of the two partner cells, these fuse to form a diploid zygote. In this review, we present our current knowledge on the processes of mating signalling, pheromone-dependent polarized growth and cell fusion in Saccharomyces cerevisiae and Schizosaccharomyces pombe, two highly divergent ascomycete yeast models. While the global architecture of the mating response is very similar between these two species, they differ significantly both in their mating physiologies and in the molecular connections between pheromone perception and downstream responses. The use of both yeast models helps enlighten both conserved solutions and species-specific adaptations to a general biological problem.

Update of the FANTOM web resource: from mammalian transcriptional landscape to its dynamic regulation.

The international Functional Annotation Of the Mammalian Genomes 4 (FANTOM4) research collaboration set out to better understand the transcriptional network that regulates macrophage differentiation and to uncover novel components of the transcriptome employing a series of high-throughput experiments. The primary and unique technique is cap analysis of gene expression (CAGE), sequencing mRNA 5'-ends with a second-generation sequencer to quantify promoter activities even in the absence of gene annotation. Additional genome-wide experiments complement the setup including short RNA sequencing, microarray gene expression profiling on large-scale perturbation experiments and ChIP-chip for epigenetic marks and transcription factors. All the experiments are performed in a differentiation time course of the THP-1 human leukemic cell line. Furthermore, we performed a large-scale mammalian two-hybrid (M2H) assay between transcription factors and monitored their expression profile across human and mouse tissues with qRT-PCR to address combinatorial effects of regulation by transcription factors. These interdependent data have been analyzed individually and in combination with each other and are published in related but distinct papers. We provide all data together with systematic annotation in an integrated view as resource for the scientific community (http://fantom.gsc.riken.jp/4/). Additionally, we assembled a rich set of derived analysis results including published predicted and validated regulatory interactions. Here we introduce the resource and its update after the initial release.

Inversion of crosshole seismic data in heterogenous environments : Comparison of waveform and ray-based approaches

High-resolution tomographic imaging of the shallow subsurface is becoming increasingly important for a wide range of environmental, hydrological and engineering applications. Because of their superior resolution power, their sensitivity to pertinent petrophysical parameters, and their far reaching complementarities, both seismic and georadar crosshole imaging are of particular importance. To date, corresponding approaches have largely relied on asymptotic, ray-based approaches, which only account for a very small part of the observed wavefields, inherently suffer from a limited resolution, and in complex environments may prove to be inadequate. These problems can potentially be alleviated through waveform inversion. We have developed an acoustic waveform inversion approach for crosshole seismic data whose kernel is based on a finite-difference time-domain (FDTD) solution of the 2-D acoustic wave equations. This algorithm is tested on and applied to synthetic data from seismic velocity models of increasing complexity and realism and the results are compared to those obtained using state-of-the-art ray-based traveltime tomography. Regardless of the heterogeneity of the underlying models, the waveform inversion approach has the potential of reliably resolving both the geometry and the acoustic properties of features of the size of less than half a dominant wavelength. Our results do, however, also indicate that, within their inherent resolution limits, ray-based approaches provide an effective and efficient means to obtain satisfactory tomographic reconstructions of the seismic velocity structure in the presence of mild to moderate heterogeneity and in absence of strong scattering. Conversely, the excess effort of waveform inversion provides the greatest benefits for the most heterogeneous, and arguably most realistic, environments where multiple scattering effects tend to be prevalent and ray-based methods lose most of their effectiveness.

General practitioners can evaluate the material, social and health dimensions of patient social status

OBJECTIVE: To identify which physician and patient characteristics are associated with physicians' estimation of their patient social status.DESIGN: Cross-sectional ulticentric survey. SETTING: Fourty-seven primary care private offices in Western Switzerland. PARTICIPANTS: Random sample of 2030 patients ≥ 16, who encountered a general practitioner (GP) between September 2010 and February 2011. MAIN MEASURES: PRIMARY OUTCOME: patient social status perceived by GPs, using the MacArthur Scale of Subjective Social Status, ranging from the bottom (0) to the top (10) of the social scale.Secondary outcome: Difference between GP's evaluation and patient's own evaluation of their social status. Potential patient correlates: material and social deprivation using the DiPCare-Q, health status using the EQ-5D, sources of income, and level of education. GP characteristics: opinion regarding patients' deprivation and its influence on health and care. RESULTS: To evaluate patient social status, GPs considered the material, social, and health aspects of deprivation, along with education level, and amount and type of income. GPs declaring a frequent reflexive consideration of their own prejudice towards deprived patients, gave a higher estimation of patients' social status (+1.0, p = 0.002). Choosing a less costly treatment for deprived patients was associated with a lower estimation (-0.7, p = 0.002). GP's evaluation of patient social status was 0.5 point higher than the patient's own estimate (p<0.0001). CONCLUSIONS: GPs can perceive the various dimensions of patient social status, although heterogeneously, according partly to their own characteristics. Compared to patients' own evaluation, GPs overestimate patient social status.

Hierarchical management of carbon sources is regulated similarly by the CbrA/B systems in Pseudomonas aeruginosa and Pseudomonas putida.

The CbrA/B system in pseudomonads is involved in the utilization of carbon sources and carbon catabolite repression (CCR) through the activation of the small RNAs crcZ in Pseudomonas aeruginosa, and crcZ and crcY in Pseudomonas putida. Interestingly, previous works reported that the CbrA/B system activity in P. aeruginosa PAO1 and P. putida KT2442 responded differently to the presence of different carbon sources, thus raising the question of the exact nature of the signal(s) detected by CbrA. Here, we demonstrated that the CbrA/B/CrcZ(Y) signal transduction pathway is similarly activated in the two Pseudomonas species. We show that the CbrA sensor kinase is fully interchangeable between the two species and, moreover, responds similarly to the presence of different carbon sources. In addition, a metabolomics analysis supported the hypothesis that CCR responds to the internal energy status of the cell, as the internal carbon/nitrogen ratio seems to determine CCR and non-CCR conditions. The strong difference found in the 2-oxoglutarate/glutamine ratio between CCR and non-CCR conditions points to the close relationship between carbon and nitrogen availability, or the relationship between the CbrA/B and NtrB/C systems, suggesting that both regulatory systems sense the same sort or interrelated signal.

Psychogenic seizures and frontal disconnection: EEG synchronisation study.

Objective Psychogenic non-epileptic seizures (PNES) are paroxysmal events that, in contrast to epileptic seizures, are related to psychological causes without the presence of epileptiform EEG changes. Recent models suggest a multifactorial basis for PNES. A potentially paramount, but currently poorly understood factor is the interplay between psychiatric features and a specific vulnerability of the brain leading to a clinical picture that resembles epilepsy. Hypothesising that functional cerebral network abnormalities may predispose to the clinical phenotype, the authors undertook a characterisation of the functional connectivity in PNES patients. Methods The authors analysed the whole-head surface topography of multivariate phase synchronisation (MPS) in interictal high-density EEG of 13 PNES patients as compared with 13 age- and sex-matched controls. MPS mapping reduces the wealth of dynamic data obtained from high-density EEG to easily readable synchronisation maps, which provide an unbiased overview of any changes in functional connectivity associated with distributed cortical abnormalities. The authors computed MPS maps for both Laplacian and common-average-reference EEGs. Results In a between-group comparison, only patchy, non-uniform changes in MPS survived conservative statistical testing. However, against the background of these unimpressive group results, the authors found widespread inverse correlations between individual PNES frequency and MPS within the prefrontal and parietal cortices. Interpretation PNES appears to be associated with decreased prefrontal and parietal synchronisation, possibly reflecting dysfunction of networks within these regions.

On the use of spatially discrete data to compute energy and mass balance

In many practical applications the state of field soils is monitored by recording the evolution of temperature and soil moisture at discrete depths. We theoretically investigate the systematic errors that arise when mass and energy balances are computed directly from these measurements. We show that, even with no measurement or model errors, large residuals might result when finite difference approximations are used to compute fluxes and storage term. To calculate the limits set by the use of spatially discrete measurements on the accuracy of balance closure, we derive an analytical solution to estimate the residual on the basis of the two key parameters: the penetration depth and the distance between the measurements. When the thickness of the control layer for which the balance is computed is comparable to the penetration depth of the forcing (which depends on the thermal diffusivity and on the forcing period) large residuals arise. The residual is also very sensitive to the distance between the measurements, which requires accurately controlling the position of the sensors in field experiments. We also demonstrate that, for the same experimental setup, mass residuals are sensitively larger than the energy residuals due to the nonlinearity of the moisture transport equation. Our analysis suggests that a careful assessment of the systematic mass error introduced by the use of spatially discrete data is required before using fluxes and residuals computed directly from field measurements.

Evaluating effects of spectral training data distribution on continuous field mapping performance

Continuous field mapping has to address two conflicting remote sensing requirements when collecting training data. On one hand, continuous field mapping trains fractional land cover and thus favours mixed training pixels. On the other hand, the spectral signature has to be preferably distinct and thus favours pure training pixels. The aim of this study was to evaluate the sensitivity of training data distribution along fractional and spectral gradients on the resulting mapping performance. We derived four continuous fields (tree, shrubherb, bare, water) from aerial photographs as response variables and processed corresponding spectral signatures from multitemporal Landsat 5 TM data as explanatory variables. Subsequent controlled experiments along fractional cover gradients were then based on generalised linear models. Resulting fractional and spectral distribution differed between single continuous fields, but could be satisfactorily trained and mapped. Pixels with fractional or without respective cover were much more critical than pure full cover pixels. Error distribution of continuous field models was non-uniform with respect to horizontal and vertical spatial distribution of target fields. We conclude that a sampling for continuous field training data should be based on extent and densities in the fractional and spectral, rather than the real spatial space. Consequently, adequate training plots are most probably not systematically distributed in the real spatial space, but cover the gradient and covariate structure of the fractional and spectral space well. (C) 2009 International Society for Photogrammetry and Remote Sensing, Inc. (ISPRS). Published by Elsevier B.V. All rights reserved.