Source-wavelet estimation during full-waveform inversion of ground-penetrating radar data

A major issue in the application of waveform inversion methods to crosshole ground-penetrating radar (GPR) data is the accurate estimation of the source wavelet. Here, we explore the viability and robustness of incorporating this step into a recently published time-domain 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 of 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.

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.

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.

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.

Full-waveform inversion of cross-hole ground-penetrating radar data to characterize a gravel aquifer close to the Thur River, Switzerland

Cross-hole radar tomography is a useful tool for mapping shallow subsurface electrical properties viz. dielectric permittivity and electrical conductivity. Common practice is to invert cross-hole radar data with ray-based tomographic algorithms using first arrival traveltimes and first cycle amplitudes. However, the resolution of conventional standard ray-based inversion schemes for cross-hole ground-penetrating radar (GPR) is limited because only a fraction of the information contained in the radar data is used. The resolution can be improved significantly by using a full-waveform inversion that considers the entire waveform, or significant parts thereof. A recently developed 2D time-domain vectorial full-waveform crosshole radar inversion code has been modified in the present study by allowing optimized acquisition setups that reduce the acquisition time and computational costs significantly. This is achieved by minimizing the number of transmitter points and maximizing the number of receiver positions. The improved algorithm was employed to invert cross-hole GPR data acquired within a gravel aquifer (4-10 m depth) in the Thur valley, Switzerland. The simulated traces of the final model obtained by the full-waveform inversion fit the observed traces very well in the lower part of the section and reasonably well in the upper part of the section. Compared to the ray-based inversion, the results from the full-waveform inversion show significantly higher resolution images. At either side, 2.5 m distance away from the cross-hole plane, borehole logs were acquired. There is a good correspondence between the conductivity tomograms and the natural gamma logs at the boundary of the gravel layer and the underlying lacustrine clay deposits. Using existing petrophysical models, the inversion results and neutron-neutron logs are converted to porosity. Without any additional calibration, the values obtained for the converted neutron-neutron logs and permittivity results are very close and similar vertical variations can be observed. The full-waveform inversion provides in both cases additional information about the subsurface. Due to the presence of the water table and associated refracted/reflected waves, the upper traces are not well fitted and the upper 2 m in the permittivity and conductivity tomograms are not reliably reconstructed because the unsaturated zone is not incorporated into the inversion domain.

A glycine-arginine domain in control of the human MRE11 DNA repair protein.

Human MRE11 is a key enzyme in DNA double-strand break repair and genome stability. Human MRE11 bears a glycine-arginine-rich (GAR) motif that is conserved among multicellular eukaryotic species. We investigated how this motif influences MRE11 function. Human MRE11 alone or a complex of MRE11, RAD50, and NBS1 (MRN) was methylated in insect cells, suggesting that this modification is conserved during evolution. We demonstrate that PRMT1 interacts with MRE11 but not with the MRN complex, suggesting that MRE11 arginine methylation occurs prior to the binding of NBS1 and RAD50. Moreover, the first six methylated arginines are essential for the regulation of MRE11 DNA binding and nuclease activity. The inhibition of arginine methylation leads to a reduction in MRE11 and RAD51 focus formation on a unique double-strand break in vivo. Furthermore, the MRE11-methylated GAR domain is sufficient for its targeting to DNA damage foci and colocalization with gamma-H2AX. These studies highlight an important role for the GAR domain in regulating MRE11 function at the biochemical and cellular levels during DNA double-strand break repair.

Regulation of Nedd4-2 self-ubiquitination and stability by a PY motif located within its HECT-domain.

Ubiquitin ligases play a pivotal role in substrate recognition and ubiquitin transfer, yet little is known about the regulation of their catalytic activity. Nedd4 (neural-precursor-cell-expressed, developmentally down-regulated 4)-2 is an E3 ubiquitin ligase composed of a C2 domain, four WW domains (protein-protein interaction domains containing two conserved tryptophan residues) that bind PY motifs (L/PPXY) and a ubiquitin ligase HECT (homologous with E6-associated protein C-terminus) domain. In the present paper we show that the WW domains of Nedd4-2 bind (weakly) to a PY motif (LPXY) located within its own HECT domain and inhibit auto-ubiquitination. Pulse-chase experiments demonstrated that mutation of the HECT PY-motif decreases the stability of Nedd4-2, suggesting that it is involved in stabilization of this E3 ligase. Interestingly, the HECT PY-motif mutation does not affect ubiquitination or down-regulation of a known Nedd4-2 substrate, ENaC (epithelial sodium channel). ENaC ubiquitination, in turn, appears to promote Nedd4-2 self-ubiquitination. These results support a model in which the inter- or intra-molecular WW-domain-HECT PY-motif interaction stabilizes Nedd4-2 by preventing self-ubiquitination. Substrate binding disrupts this interaction, allowing self-ubiquitination of Nedd4-2 and subsequent degradation, resulting in down-regulation of Nedd4-2 once it has ubiquitinated its target. These findings also point to a novel mechanism employed by a ubiquitin ligase to regulate itself differentially compared with substrate ubiquitination and stability.

The specificity of TRIM5 alpha-mediated restriction is influenced by its coiled-coil domain.

Retroviruses are both powerful evolutionary forces and dangerous threats to genome integrity. As such, they have imposed strong selective pressure on their hosts, notably triggering the emergence of restriction factors, such as TRIM5 alpha, that act as potent barriers to their cross-species transmission. TRIM5 alpha orthologues from different primates have distinct retroviral restriction patterns, largely dictated by the sequence of their C-terminal PRYSPRY domain, which binds the capsid protein of incoming virions. Here, by combining genetic and functional analyses of human and squirrel monkey TRIM5 alpha, we demonstrate that the coiled-coil domain of this protein, thus far essentially known for mediating oligomerization, also conditions the spectrum of antiretroviral activity. Furthermore, we identify three coiled-coil residues responsible for this effect, one of which has been under positive selection during primate evolution, notably in New World monkeys. These results indicate that the PRYSPRY and coiled-coil domains cooperate to determine the specificity of TRIM5 alpha-mediated capture of retroviral capsids, shedding new light on this complex event.

Evaluation of the reconstruction limits of a frequency-independent crosshole georadar waveform inversion scheme in the presence of dispersion

Waveform tomographic imaging of crosshole georadar data is a powerful method to investigate the shallow subsurface because of its ability to provide images of pertinent petrophysical parameters with extremely high spatial resolution. All current crosshole georadar waveform inversion strategies are based on the assumption of frequency-independent electromagnetic constitutive parameters. However, in reality, these parameters are known to be frequency-dependent and complex and thus recorded georadar data may show significant dispersive behavior. In this paper, we evaluate synthetically the reconstruction limits of a recently published crosshole georadar waveform inversion scheme in the presence of varying degrees of dielectric dispersion. Our results indicate that, when combined with a source wavelet estimation procedure that provides a means of partially accounting for the frequency-dependent effects through an "effective" wavelet, the inversion algorithm performs remarkably well in weakly to moderately dispersive environments and has the ability to provide adequate tomographic reconstructions.

Dirigent domain-containing protein is part of the machinery required for formation of the lignin-based Casparian strip in the root.

The endodermis acts as a "second skin" in plant roots by providing the cellular control necessary for the selective entry of water and solutes into the vascular system. To enable such control, Casparian strips span the cell wall of adjacent endodermal cells to form a tight junction that blocks extracellular diffusion across the endodermis. This junction is composed of lignin that is polymerized by oxidative coupling of monolignols through the action of a NADPH oxidase and peroxidases. Casparian strip domain proteins (CASPs) correctly position this biosynthetic machinery by forming a protein scaffold in the plasma membrane at the site where the Casparian strip forms. Here, we show that the dirigent-domain containing protein, enhanced suberin1 (ESB1), is part of this machinery, playing an essential role in the correct formation of Casparian strips. ESB1 is localized to Casparian strips in a CASP-dependent manner, and in the absence of ESB1, disordered and defective Casparian strips are formed. In addition, loss of ESB1 disrupts the localization of the CASP1 protein at the casparian strip domain, suggesting a reciprocal requirement for both ESB1 and CASPs in forming the casparian strip domain.

Contribution of TIR domain-containing adapter inducing IFN-beta-mediated IL-18 release to LPS-induced liver injury in mice.

BACKGROUND/AIMS: After treatment with heat-killed Propionibacterium acnes mice show dense hepatic granuloma formation. Such mice develop liver injury in an interleukin (IL)-18-dependent manner after challenge with a sublethal dose LPS. As previously shown, LPS-stimulated Kupffer cells secrete IL-18 depending on caspase-1 and Toll-like receptor (TLR)-4 but independently of its signal adaptor myeloid differentiation factor 88 (MyD88), suggesting importance of another signal adaptor TIR domain-containing adapter inducing IFN-beta (TRIF). Nalp3 inflammasome reportedly controls caspase-1 activation. Here we investigated the roles of MyD88 and TRIF in P. acnes-induced hepatic granuloma formation and LPS-induced caspase-1 activation for IL-18 release. METHODS: Mice were sequentially treated with P. acnes and LPS, and their serum IL-18 levels and liver injuries were determined by ELISA and ALT/AST measurement, respectively. Active caspase-1 in LPS-stimulated Kupffer cells was determined by Western blotting. RESULTS: Macrophage-ablated mice lacked P. acnes-induced hepatic granuloma formation and LPS-induced serum IL-18 elevation and liver injury. Myd88(-/-) Kupffer cells, but not Trif(-/-) cells, exhibited normal caspase-1 activation upon TLR4 engagement in vitro. Myd88(-/-) mice failed to develop hepatic granulomas after P. acnes treatment and liver injury induced by LPS challenge. In contrast, Trif(-/-) mice normally formed the hepatic granulomas, but could not release IL-18 or develop the liver injury. Nalp3(-/-) mice showed the same phenotypes of Trif(-/-) mice. CONCLUSIONS: Propionibacterium acnes treatment MyD88-dependently induced hepatic granuloma formation. Subsequent LPS TRIF-dependently activated caspase-1 via Nalp3 inflammasome and induced IL-18 release, eventually leading to the liver injury.