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.

Focused time-lapse inversion of radio and audio magnetotelluric data

Geoelectrical techniques are widely used to monitor groundwater processes, while surprisingly few studies have considered audio (AMT) and radio (RMT) magnetotellurics for such purposes. In this numerical investigation, we analyze to what extent inversion results based on AMT and RMT monitoring data can be improved by (1) time-lapse difference inversion; (2) incorporation of statistical information about the expected model update (i.e., the model regularization is based on a geostatistical model); (3) using alternative model norms to quantify temporal changes (i.e., approximations of l(1) and Cauchy norms using iteratively reweighted least-squares), (4) constraining model updates to predefined ranges (i.e., using Lagrange Multipliers to only allow either increases or decreases of electrical resistivity with respect to background conditions). To do so, we consider a simple illustrative model and a more realistic test case related to seawater intrusion. The results are encouraging and show significant improvements when using time-lapse difference inversion with non l(2) model norms. Artifacts that may arise when imposing compactness of regions with temporal changes can be suppressed through inequality constraints to yield models without oscillations outside the true region of temporal changes. Based on these results, we recommend approximate l(1)-norm solutions as they can resolve both sharp and smooth interfaces within the same model. (C) 2012 Elsevier B.V. All rights reserved.

Selective coronary artery plaque visualization and differentiation by contrast-enhanced inversion prepared MRI.

AIMS: We sought to evaluate the utility of contrast-enhanced coronary magnetic resonance imaging (CE-MRI) for selective visualization and non-invasive differentiation of atherosclerotic coronary plaque in humans. METHODS AND RESULTS: Nine patients with coronary artery disease (CAD) as confirmed by X-ray angiography and multidetector computed tomography (MDCT) were studied by T1-weighted black blood inversion recovery coronary MRI before (N-IR) and after administration of Gd-DTPA (CE-IR). Plaques were categorized as calcified, non-calcified, and mixed based on their Hounsfield number derived from MDCT. With MDCT, a total of 29 plaques were identified, including calcified (n=6), non-calcified (n=6), and mixed calcified/non-calcified (n=17). On N-IR MRI, 26 plaques (90%) were dark, whereas three plaques (two non-calcified and one mixed) appeared bright. On CE-MRI, 13/29 (45%) plaques, 11 of which were mixed, one non-calcified, and one calcified showed contrast uptake. All others remained dark. CONCLUSION: In this preliminary study, we demonstrate the potential utility of CE-IR MRI for selective plaque visualization and differentiation of plaque types. The observed contrast uptake may be associated with endothelial dysfunction, neovascularization, inflammation, and/or fibrosis.

Landmark detection for fusion of fundus and MRI toward a patient-specific multimodal eye model.

Ophthalmologists typically acquire different image modalities to diagnose eye pathologies. They comprise, e.g., Fundus photography, optical coherence tomography, computed tomography, and magnetic resonance imaging (MRI). Yet, these images are often complementary and do express the same pathologies in a different way. Some pathologies are only visible in a particular modality. Thus, it is beneficial for the ophthalmologist to have these modalities fused into a single patient-specific model. The goal of this paper is a fusion of Fundus photography with segmented MRI volumes. This adds information to MRI that was not visible before like vessels and the macula. This paper contributions include automatic detection of the optic disc, the fovea, the optic axis, and an automatic segmentation of the vitreous humor of the eye.

The SM protein Vps33 and the t-SNARE H(abc) domain promote fusion pore opening.

Intracellular membrane fusion proceeds via distinct stages of membrane docking, hemifusion and fusion pore opening and depends on interacting families of Rab, SNARE and SM proteins. Trans-SNARE complexes dock the membranes in close apposition. Efficient fusion requires further SNARE-associated proteins. They might increase the number of trans-SNARE complexes or the fusogenic potential of a single SNARE complex. We investigated the contributions of the SM protein Vps33 to hemifusion and pore opening between yeast vacuoles. Mutations in Vps33 that weaken its interactions with the SNARE complex allowed normal trans-SNARE pairing and lipid mixing but retarded content mixing. Deleting the H(abc) domain of the vacuolar t-SNARE Vam3, which interacts with Vps33, had the same effect. This suggests that SM proteins promote fusion pore opening by enhancing the fusogenic activity of a SNARE complex. They should thus be considered integral parts of the fusion machinery.

Simple method for detection of MYH11 DNA rearrangements in patients with inv(16)(p13q22) and acute myeloid leukemia.

The pericentric inversion on chromosome 16 [inv(16)(p13q22)] and related t(16;16)(p13;q22) are recurrent aberrations associated with acute myeloid leukemia (AML) M4 Eo. Both abberations result in a fusion of the core binding factor beta (CBFB) and smooth muscle myosin heavy chain gene (MYH11). A selected genomic 6.9-kb BamHl probe detects MYH11 DNA rearrangements in 18 of 19 inv(16)/t(16;16) patients tested using HindIII digested DNA. The rearranged fragments were not detectable after remission in two cases tested, while they were present after relapse in one of these two cases tested.

Structural joint inversion of time-lapse crosshole ERT and GPR traveltime data

Time-lapse geophysical monitoring and inversion are valuable tools in hydrogeology for monitoring changes in the subsurface due to natural and forced (tracer) dynamics. However, the resulting models may suffer from insufficient resolution, which leads to underestimated variability and poor mass recovery. Structural joint inversion using cross-gradient constraints can provide higher-resolution models compared with individual inversions and we present the first application to time-lapse data. The results from a synthetic and field vadose zone water tracer injection experiment show that joint 3-D time-lapse inversion of crosshole electrical resistance tomography (ERT) and ground penetrating radar (GPR) traveltime data significantly improve the imaged characteristics of the point injected plume, such as lateral spreading and center of mass, as well as the overall consistency between models. The joint inversion method appears to work well for cases when one hydrological state variable (in this case moisture content) controls the time-lapse response of both geophysical methods. Citation: Doetsch, J., N. Linde, and A. Binley (2010), Structural joint inversion of time-lapse crosshole ERT and GPR traveltime data, Geophys. Res. Lett., 37, L24404, doi: 10.1029/2010GL045482.

Waveform inversion of crosshole georadar data: influence of source wavelet variability and the suitability of a single wavelet assumption

Waveform-based tomographic imaging of crosshole georadar data is a powerful method to investigate the shallow subsurface because of its ability to provide images of electrical properties in near-surface environments with unprecedented spatial resolution. A critical issue with waveform inversion is the a priori unknown source signal. Indeed, the estimation of the source pulse is notoriously difficult but essential for the effective application of this method. Here, we explore the viability and robustness of a recently proposed deconvolution-based procedure to estimate the source pulse during waveform inversion of crosshole georadar data, where changes in wavelet shape with location as a result of varying near-field conditions and differences in antenna coupling may be significant. Specifically, we examine whether a single, average estimated source current function can adequately represent the pulses radiated at all transmitter locations during a crosshole georadar survey, or whether a separate source wavelet estimation should be performed for each transmitter gather. Tests with synthetic and field data indicate that remarkably good tomographic reconstructions can be obtained using a single estimated source pulse when moderate to strong variability exists in the true source signal with antenna location. Only in the case of very strong variability in the true source pulse are tomographic reconstructions clearly improved by estimating a different source wavelet for each transmitter location.

Model-based Segmentation and Fusion of 3D Computed Tomography and 3D Ultrasound of the Eye for Radiotherapy Planning

Computed Tomography (CT) represents the standard imaging modality for tumor volume delineation for radiotherapy treatment planning of retinoblastoma despite some inherent limitations. CT scan is very useful in providing information on physical density for dose calculation and morphological volumetric information but presents a low sensitivity in assessing the tumor viability. On the other hand, 3D ultrasound (US) allows a highly accurate definition of the tumor volume thanks to its high spatial resolution but it is not currently integrated in the treatment planning but used only for diagnosis and follow-up. Our ultimate goal is an automatic segmentation of gross tumor volume (GTV) in the 3D US, the segmentation of the organs at risk (OAR) in the CT and the registration of both modalities. In this paper, we present some preliminary results in this direction. We present 3D active contour-based segmentation of the eye ball and the lens in CT images; the presented approach incorporates the prior knowledge of the anatomy by using a 3D geometrical eye model. The automated segmentation results are validated by comparing with manual segmentations. Then, we present two approaches for the fusion of 3D CT and US images: (i) landmark-based transformation, and (ii) object-based transformation that makes use of eye ball contour information on CT and US images.

Crocidura cossyrensis Contoli, 1989 (Mammalia, Soridiae): karyotype, biochemical genetics and hybridization experiments

Crocidura cossyrensis Contoli, 1989 (Mammalia, Soricidae): karyotype, biochemical genetics and hybridization experiments. - The shrew Crocidura cossyrensis Contoli, 1989 from Pantelleria (I), a Mediterranean island 100 km south of Sicily and 70 km west from Tunisia, was investigated in order to understand its origin and its relationship with C. russula from Tunisia, Morocco and Switzerland. With the exception of a single heterozygote centric fusion, C. cossyrensis had a karyotype identical with that of C russula from Tunisia (2N = 42, NF = 70 to 72), but it differed from C russula from Morocco and Switzerland (2N = 42, NF = 60). The former have 5-6 pairs of chromosomes with small arms that are acrocentric in the latter. Genetic comparisons with allozyme data revealed small genetic distance (0.04) between C cossyrensis and C russula from Tunisia. In contrast, this eastern clade (Tunisia and Pantelleria) is separated from the western clade (Switzerland and Morocco) by a genetic distance of 0.14. A hybridization experiment between shrews from Pantelleria and Switzerland lead rapidly to an F1 generation. From 12 F1 hybrids that were backcrossed, females reproduced normally, but none of the males did so. Concluding from the results, C. cossyrensis from Pantelleria and C. russula cf. agilis from Tunisia belong to the same taxon that may have reached the differentiation of a biological species within the C. russula group. More geographic samples are needed to determine the definitive taxonomic positions of these shrews.

Sequence analysis and expression of the attachment and fusion proteins of canine distemper virus wild-type strain A75/17.

The biological properties of wild-type A75/17 and cell culture-adapted Onderstepoort canine distemper virus differ markedly. To learn more about the molecular basis for these differences, we have isolated and sequenced the protein-coding regions of the attachment and fusion proteins of wild-type canine distemper virus strain A75/17. In the attachment protein, a total of 57 amino acid differences were observed between the Onderstepoort strain and strain A75/17, and these were distributed evenly over the entire protein. Interestingly, the attachment protein of strain A75/17 contained an extension of three amino acids at the C terminus. Expression studies showed that the attachment protein of strain A75/17 had a higher apparent molecular mass than the attachment protein of the Onderstepoort strain, in both the presence and absence of tunicamycin. In the fusion protein, 60 amino acid differences were observed between the two strains, of which 44 were clustered in the much smaller F2 portion of the molecule. Significantly, the AUG that has been proposed as a translation initiation codon in the Onderstepoort strain is an AUA codon in strain A75/17. Detailed mutation analyses showed that both the first and second AUGs of strain A75/17 are the major translation initiation sites of the fusion protein. Similar analyses demonstrated that, also in the Onderstepoort strain, the first two AUGs are the translation initiation codons which contribute most to the generation of precursor molecules yielding the mature form of the fusion protein.

Inference of chromosomal inversion dynamics from Pool-Seq data in natural and laboratory populations of Drosophila melanogaster.

Sequencing of pools of individuals (Pool-Seq) represents a reliable and cost-effective approach for estimating genome-wide SNP and transposable element insertion frequencies. However, Pool-Seq does not provide direct information on haplotypes so that, for example, obtaining inversion frequencies has not been possible until now. Here, we have developed a new set of diagnostic marker SNPs for seven cosmopolitan inversions in Drosophila melanogaster that can be used to infer inversion frequencies from Pool-Seq data. We applied our novel marker set to Pool-Seq data from an experimental evolution study and from North American and Australian latitudinal clines. In the experimental evolution data, we find evidence that positive selection has driven the frequencies of In(3R)C and In(3R)Mo to increase over time. In the clinal data, we confirm the existence of frequency clines for In(2L)t, In(3L)P and In(3R)Payne in both North America and Australia and detect a previously unknown latitudinal cline for In(3R)Mo in North America. The inversion markers developed here provide a versatile and robust tool for characterizing inversion frequencies and their dynamics in Pool-Seq data from diverse D. melanogaster populations.

Free-breathing renal MR angiography with steady-state free-precession (SSFP) and slab-selective spin inversion: initial results.

BACKGROUND: The aim of our study was the investigation of a novel navigator-gated three-dimensional (3D) steady-state free-precession (SSFP) sequence for free-breathing renal magnetic resonance angiography (MRA) without contrast medium, and to examine the advantage of an additional inversion prepulse for improved contrast. METHODS: Eight healthy volunteers (mean age 29 years) and eight patients (mean age 53 years) were investigated on a 1.5 Tesla MR system (ACS-NT, Philips, Best, The Netherlands). Renal MRA was performed using three navigator-gated free-breathing cardiac-triggered 3D SSFP sequences [repetition time (TR) = 4.4 ms, echo time (TE) = 2.2 ms, flip angle 85 degrees, spatial resolution 1.25 x 1.25 x 4.0 mm(3), scanning time approximately 1 minute 30 seconds]. The same sequence was performed without magnetization preparation, with a non-slab selective and a slab-selective inversion prepulse. Signal-to-noise ratio (SNR), contrast-to-noise (CNR) vessel length, and subjective image quality were compared. RESULTS: Three-dimensional SSFP imaging combined with a slab-selective inversion prepulse enabled selective and high contrast visualization of the renal arteries, including the more distal branches. Standard SSFP imaging without magnetization preparation demonstrated overlay by veins and renal parenchyma. A non-slab-selective prepulse abolished vessel visualization. CNR in SSFP with slab-selective inversion was 43.6 versus 10.6 (SSFP without magnetization preparation) and 0.4 (SSFP with non-slab-selective inversion), P < 0.008. CONCLUSION: Navigator-gated free-breathing cardiac-triggered 3D SSFP imaging combined with a slab-selective inversion prepulse is a novel, fast renal MRA technique without the need for contrast media.