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.

Dealing with Errors in QCA

This paper discusses five strategies to deal with five types of errors in Qualitative Comparative Analysis (QCA): condition errors, systematic errors, random errors, calibration errors, and deviant case errors. These strategies are the comparative inspection of complex, intermediary, and parsimonious solutions; the use of an adjustment factor, the use of probabilistic criteria, the test of the robustness of calibration parameters, and the use of a frequency threshold for observed combinations of conditions. The strategies are systematically reviewed, assessed, and evaluated as regards their applicability, advantages, limitations, and complementarities.

Free Flows, Limited Diversification: Explaining the Fall and Rise of Stock Market Correlations, 1890-2001

Using a new dataset on capital account openness, we investigate why equity return correlations changed over the last century. Based on a new, long-run dataset on capital account regulations in a group of 16 countries over the period 1890-2001, we show that correlations increase as financial markets are liberalized. These findings are robust to controlling for both the Forbes-Rigobon bias and global averages in equity return correlations. We test the robustness of our conclusions, and show that greater synchronization of fundamentals is not the main cause of increasing correlations. These results imply that the home bias puzzle may be smaller than traditionally claimed.

Are There Classical Business Cycles?

The aim of this paper is to test formally the classical business cycle hypothesis, using data from industrialized countries for the time period since 1960. The hypothesis is characterized by the view that the cyclical structure in GDP is concentrated in the investment series: fixed investment has typically a long cycle, while the cycle in inventory investment is shorter. To check the robustness of our results, we subject the data for 15 OECD countries to a variety of detrending techniques. While the hypothesis is not confirmed uniformly for all countries, there is a considerably high number for which the data display the predicted pattern. None of the countries shows a pattern which can be interpreted as a clear rejection of the classical hypothesis.

Error and attack tolerance of layered complex networks.

Many complex systems may be described by not one but a number of complex networks mapped on each other in a multi-layer structure. Because of the interactions and dependencies between these layers, the state of a single layer does not necessarily reflect well the state of the entire system. In this paper we study the robustness of five examples of two-layer complex systems: three real-life data sets in the fields of communication (the Internet), transportation (the European railway system), and biology (the human brain), and two models based on random graphs. In order to cover the whole range of features specific to these systems, we focus on two extreme policies of system's response to failures, no rerouting and full rerouting. Our main finding is that multi-layer systems are much more vulnerable to errors and intentional attacks than they appear from a single layer perspective.

Fat signal suppression for coronary MRA at 3T using a water-selective adiabatic T2 -preparation technique.

PURPOSE: To improve fat saturation in coronary MRA at 3T by using a spectrally selective adiabatic T2 -Prep (WSA-T2 -Prep). METHODS: A conventional adiabatic T2 -Prep (CA-T2 -Prep) was modified, such that the excitation and restoration pulses were of differing bandwidths. On-resonance spins are T2 -Prepared, whereas off-resonance spins, such as fat, are spoiled. This approach was combined with a CHEmically Selective Saturation (CHESS) pulse to achieve even greater fat suppression. Numerical simulations were followed by phantom validation and in vivo coronary MRA. RESULTS: Numerical simulations demonstrated that augmenting a CHESS pulse with a WSA-T2 -Prep improved robustness to B1 inhomogeneities and that this combined fat suppression was effective over a broader spectral range than that of a CHESS pulse in a conventional T2 -Prepared sequence. Phantom studies also demonstrated that the WSA-T2 -Prep+CHESS combination produced greater fat suppression across a range of B1 values than did a CA-T2 -Prep+CHESS combination. Lastly, in vivo measurements demonstrated that the contrast-to-noise ratio between blood and myocardium was not adversely affected by using a WSA-T2 -Prep, despite the improved abdominal and epicardial fat suppression. Additionally, vessel sharpness improved. CONCLUSION: The proposed WSA-T2 -Prep method was shown to improve fat suppression and vessel sharpness as compared to a CA-T2 -Prep technique, and to also increase fat suppression when combined with a CHESS pulse.

Role of myosin V in actin cable organization in fission yeast

Functional specialization is tightly linked to the ability of eukaryotic cells to acquire a particular shape. Cell morphogenesis, in turn, relies on the capacity to establish and maintain cell "polarity", which is achieved by orienting the trafficking of signaling molecules and organelles towards specific cellular locations and/or membrane domains. The "oriented" transport is based upon cytoskeletal polymers, microtubules and actin filaments, which serve as tracks for molecular motors. These latter generate motion that is translated either into pulling forces or directed transport. Fission yeast, a rod-like unicellular eukaryote, shapes itself by restricting growth at cell tips through the concerted activity of microtubules and actin cables. Microtubules, which assemble into 2-6 bundles and run parallel to the long axis of the cell, serve to orient growth to the tips. Growth is supported by the actin cytoskeleton, which provides tracks, the cables, for motor-based transport of secretory vesicles. The molecular motors, which bind cargos and deliver them to the tips along cables, are also known as type V myosins (hereafter indicated as myosin V). How the bundles of parallel actin filaments, i.e. the cables, extend from the tips through the cell and whether they serve any other purpose, besides providing tracks, is poorly understood. It is also unclear how the crosstalk between the two cytoskeletal systems is achieved. These are the basic questions I addressed during my PhD. The first part of the thesis work (Chapter two) suggests that the sole function of actin cables in polarized growth is to serve as tracks for motors. The data indicate that cells may have evolved two cytoskeletal systems to provide robustness to the polarization process but in principle a unique cytoskeleton might have been able to direct and support polarized growth. How actin cables are organized within the cell to optimize cargo transport is addressed later on (Chapter three). The major finding, based on the actin cable defect of cells lacking myosin Vs, is that actin filaments self-organize through the activity of the transport motors. In fact, by delivering cargos to cell tips and exerting physical pulling forces on actin filaments, Myosin Vs contribute not only to polarize cargo transport but also actin tracks. Among the cargos transported by Myosin V, which may be relevant to its function in organizing cables, there is likely the endoplasmic reticulum (ER). Actin cables, which run parallel to cortical ER, may serve as tracks for Myosin V. Myosin V-driven displacement, in turn, may account for the dynamic expansion and organization of ER during polarized growth as suggested in Chapter four. The last part of the work (Chapter five) highlights the existence of a crosstalk between actin and microtubules. In absence of myosin V, indeed, microtubules contribute to actin cable organization, likely playing a scaffolding/tethering function. Whether or not the kinesin 1, Klp3, plays any role in such process has to be demonstrated. In conclusion the work proposes a novel role for myosin Vs in actin organization, besides its transport function, and provides molecular tools to further dissect the role of this type of myosin in fission yeast. - La spécialisation fonctionnelle est étroitement connectée à la capacité des cellules eucaryotes d'acquérir une forme particulière. La morphogenèse cellulaire à son tour, est basée sur la capacité d'établir et de maintenir la polarité cellulaire, polarité réalisée en orientant le trafic des molécules signales et des organelles vers des zones cellulaires spécifiques. Ce transport directionnel dépend des polymères du cytosquelette, microtubules et microfilaments, qui servent comme des voies pour les moteurs moléculaires. Ces derniers engendrent du mouvement, traduit soit en force de traction soit en transport directionnel. La levure fissipare, un eucaryote unicellulaire en forme de bâtonnet, acquière sa forme en limitant sa croissance aux extrémités par l'action concertée des microtubules et de l'actine. Les microtubules, qui s'assemblent de façon antiparallèle et parcourent la cellule parallèlement à l'axe longitudinal, servent à orienter la croissance aux extrémités. Cette croissance est permise par le cytosquelette d'actine, fournissant des voies, les câbles, pour le transport actif des vésicules de sécrétion. Les moteurs moléculaires, responsables de ce transport actif sont aussi appelés myosines de type V (par la suite appelés myosines V). La manière dont ces câbles s'étendent depuis l'extrémité jusqu'à l'intérieur de la cellule est peu connue. De plus, on ignore également si ces câbles présentent une fonction autre que le transport. L'interaction entre les deux cytosquelettes est également obscure. Ce sont ces questions de base auxquelles j'ai tenté de répondre lors de ma thèse. La première partie de cette thèse (chapitre II) suggère que les câbles d'actine, pendant la croissance polarisée, fonctionnent uniquement comme des voies pour les moteurs moléculaires. Les données indiqueraient que les cellules ont fait évoluer deux systèmes de cytosquelette pour assurer plus de robustesse au processus de polarisation, bien que, comme nous le verrons, un système unique est suffisant. Au chapitre III, nous verrons comment les câbles d'actine sont organisés à l'intérieur de la cellule afin d'optimiser le transport des cargo. La découverte majeure, réalisée en observant des cellules dont la myosine V fait défaut, est que ces filaments d'actine s'auto organisent grâce au passage des moteurs moléculaires le long de ces voies. En réalité, en délivrant les cargos aux extrémités de la cellule et en exerçant des forces de traction sur les câbles, les myosines V contribuent non seulement à polariser le transport mais également à polariser les voies elles mêmes. Nous verrons également au chapitre IV, que parmi les cargos importants pour l'organisation des câbles, il y aurait le réticulum endoplasmique (RE). En effet, les câbles d'actine, qui s'étalent parallèlement au RE cortical, pourraient servir comme voie pour la myosine V. Cette dernière en retour pourrait être responsable de l'expansion dynamique et de l'organisation du RE pendant la croissance polarisée.

High-resolution 3D whole-heart coronary MRA: a study on the combination of data acquisition in multiple breath-holds and 1D residual respiratory motion compensation.

OBJECT: To study a scan protocol for coronary magnetic resonance angiography based on multiple breath-holds featuring 1D motion compensation and to compare the resulting image quality to a navigator-gated free-breathing acquisition. Image reconstruction was performed using L1 regularized iterative SENSE. MATERIALS AND METHODS: The effects of respiratory motion on the Cartesian sampling scheme were minimized by performing data acquisition in multiple breath-holds. During the scan, repetitive readouts through a k-space center were used to detect and correct the respiratory displacement of the heart by exploiting the self-navigation principle in image reconstruction. In vivo experiments were performed in nine healthy volunteers and the resulting image quality was compared to a navigator-gated reference in terms of vessel length and sharpness. RESULTS: Acquisition in breath-hold is an effective method to reduce the scan time by more than 30 % compared to the navigator-gated reference. Although an equivalent mean image quality with respect to the reference was achieved with the proposed method, the 1D motion compensation did not work equally well in all cases. CONCLUSION: In general, the image quality scaled with the robustness of the motion compensation. Nevertheless, the featured setup provides a positive basis for future extension with more advanced motion compensation methods.

Fusion of data sets in multivariate linear regression with errors-in-variables

We consider the application of normal theory methods to the estimation and testing of a general type of multivariate regressionmodels with errors--in--variables, in the case where various data setsare merged into a single analysis and the observable variables deviatepossibly from normality. The various samples to be merged can differ on the set of observable variables available. We show that there is a convenient way to parameterize the model so that, despite the possiblenon--normality of the data, normal--theory methods yield correct inferencesfor the parameters of interest and for the goodness--of--fit test. Thetheory described encompasses both the functional and structural modelcases, and can be implemented using standard software for structuralequations models, such as LISREL, EQS, LISCOMP, among others. An illustration with Monte Carlo data is presented.

International business cycles, financial markets and household production

This paper investigates the properties of an international real business cycle model with household production. We show that a model with disturbances to both market and household technologies reproduces the main regularities of the data and improves existing models in matching international consumption, investment and output correlations without irrealistic assumptions on the structure of international financial markets. Sensitivity analysis shows the robustness of the results to alternative specifications of the stochastic processes for the disturbances and to variations of unmeasured parameters within a reasonable range.

Multi-sample analysis of moment-structures: Asymptotic validity of inferences based on second-order moments

In moment structure analysis with nonnormal data, asymptotic valid inferences require the computation of a consistent (under general distributional assumptions) estimate of the matrix $\Gamma$ of asymptotic variances of sample second--order moments. Such a consistent estimate involves the fourth--order sample moments of the data. In practice, the use of fourth--order moments leads to computational burden and lack of robustness against small samples. In this paper we show that, under certain assumptions, correct asymptotic inferences can be attained when $\Gamma$ is replaced by a matrix $\Omega$ that involves only the second--order moments of the data. The present paper extends to the context of multi--sample analysis of second--order moment structures, results derived in the context of (simple--sample) covariance structure analysis (Satorra and Bentler, 1990). The results apply to a variety of estimation methods and general type of statistics. An example involving a test of equality of means under covariance restrictions illustrates theoretical aspects of the paper.

Are these classical business cycles?

The aim of this paper is to test formally the classical business cyclehypothesis, using data from industrialized countries for the timeperiod since 1960. The hypothesis is characterized by the view that the cyclical structure in GDP is concentrated in the investment series: fixed investment has typically a long cycle, while the cycle in inventory investment is shorter. To check the robustness of our results, we subject the data for 15 OECD countries to a variety of detrending techniques. While the hypothesis is not confirmed uniformly for all countries, there is a considerably high number for which the data display the predicted pattern. None of the countries shows a pattern which can be interpreted as a clear rejection of the classical hypothesis.

Label-free cytotoxicity screening assay by digital holographic microscopy

Abstract We introduce a label-free technology based on digital holographic microscopy (DHM) with applicability for screening by imaging, and we demonstrate its capability for cytotoxicity assessment using mammalian living cells. For this first high content screening compatible application, we automatized a digital holographic microscope for image acquisition of cells using commercially available 96-well plates. Data generated through both label-free DHM imaging and fluorescence-based methods were in good agreement for cell viability identification and a Z'-factor close to 0.9 was determined, validating the robustness of DHM assay for phenotypic screening. Further, an excellent correlation was obtained between experimental cytotoxicity dose-response curves and known IC values for different toxic compounds. For comparable results, DHM has the major advantages of being label free and close to an order of magnitude faster than automated standard fluorescence microscopy.

On the robust effects of technology shocks on hours worked and output

We analyze the effects of neutral and investment-specific technology shockson hours and output. Long cycles in hours are captured in a variety of ways.Hours robustly fall in response to neutral shocks and robustly increase inresponse to investment specific shocks. The percentage of the variance ofhours (output) explained by neutral shocks is small (large); the opposite istrue for investment specific shocks. News shocks are uncorrelated with theestimated technology shocks.

On the justice of voting systems

What are the best voting systems in terms of utilitarianism? Or in terms of maximin, or maximax? We study these questions for the case of three alternatives and a class of structurally equivalent voting rules. We show that plurality, arguably the most widely used voting system, performs very poorly in terms of remarkable ideals of justice, such as utilitarianism or maximin, and yet is optimal in terms of maximax. Utilitarianism is bestapproached by a voting system converging to the Borda count, while the best way to achieve maximin is by means of a voting system converging to negative voting. We study the robustness of our results across different social cultures, measures of performance, and population sizes.