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.

Estimating overidentified, nonrecursive, time-varying coefficients structural VARs

This paper provides a method to estimate time varying coefficients structuralVARs which are non-recursive and potentially overidentified. The procedureallows for linear and non-linear restrictions on the parameters, maintainsthe multi-move structure of standard algorithms and can be used toestimate structural models with different identification restrictions. We studythe transmission of monetary policy shocks and compare the results with thoseobtained with traditional methods.

Study of the combined radial post-feeding dispersion of the blowflies Chrysomya megacephala (Fabricius) and C. albiceps (Wiedemann) (Diptera, Calliphoridae)

Blowflies use discrete and ephemeral substrates to feed their larvae. After they run out of food, the larvae begin to disperse in order to find adequate places for pupation or additional food sources, a process named post-feeding larval dispersion. Some important aspects of this process were studied in a circular arena allowing the combined radial post-feeding dispersion from the center of the arena of C. albiceps and C. megacephala larvae. To determine the location of each pupa, the arena was divided in 72 identical sections starting from the center. The distance from the center, the depth and weight of each pupa were evaluated. Statistical tests were done to verify the relation between weight, depth and distance for pupation. From the total an average of 976 larvae released (488 for each species) were collected considering both experiments 456 C. megacephala pupae and 488 of C. albiceps. This demonstrates that C. albiceps probably preyed on 32 C. megacephala larvae during post-feeding dispersion. The study of this dispersion process can be used to estimate the postmortem interval (PMI) of human cadavers in legal medicine.

Measuring time-varying economic fears with consumption-based stochastic discount factors

This paper analyzes empirically the volatility of consumption-based stochastic discount factors as a measure of implicit economic fears by studying its relationship with future economic and stock market cycles. Time-varying economic fears seem to be well captured by the volatility of stochastic discount factors. In particular, the volatility of recursive utility-based stochastic discount factor with contemporaneous growth explains between 9 and 34 percent of future changes in industrial production at short and long horizons respectively. They also explain ex-ante uncertainty and risk aversion. However, future stock market cycles are better explained by a similar stochastic discount factor with long-run consumption growth. This specification of the stochastic discount factor presents higher volatility and lower pricing errors than the specification with contemporaneous consumption growth.

Speculative dynamics in the term structure of interest rates

When long maturity bonds are traded frequently and traders have non-nestedinformation sets, speculative behavior in the sense of Harrison and Kreps (1978) arises.Using a term structure model displaying such speculative behavior, this paper proposesa conceptually and observationally distinct new mechanism generating time varying predictableexcess returns. It is demonstrated that (i) dispersion of expectations about futureshort rates is sufficient for individual traders to systematically predict excess returns and(ii) the new term structure dynamics driven by speculative trade is orthogonal to publicinformation in real time, but (iii) can nevertheless be quantified using only publicly availableyield data. The model is estimated using monthly data on US short to medium termTreasuries from 1964 to 2007 and it provides a good fit of the data. Speculative dynamicsare found to be quantitatively important, potentially accounting for a substantial fractionof the variation of bond yields and appears to be more important at long maturities.

Dispersion of natural arsenic in the Malcantone watershed, Southern Switzerland: field evidence for repeated sorption-desorption and oxidation-reduction processes

In recent years, elevated arsenic concentrations have been found in waters and soils of many, countries, often resulting in a health threat for the local population. Switzerland is not an exception and this paper deals with the release and subsequent fate of arsenic in a 200-km(2) mountainous watershed, characterized by crystalline silicate rocks (gneisses, schists, amphibolites) that contain abundant As-bearing sulfide ore deposits, some of which have been mined for iron and gold in the past. Using analytical methods common for mineralogical, ground water and soil studies (XRD, XRF, XAS-XANES and -EXAFS, electron microprobe, extraction, ICP, AAS with hydride generator, ion chromatography), seven different field situations and related dispersion processes of natural arsenic have been studied: (1) release by rock weathering, (2) transport and deposition by water and ice; (3) release of As to the ground and surface water due to increasing pH; (4) accumulation in humic soil horizons; (5) remobilization by reduction in water-saturated soils and stagnant ground waters; (6) remobilization by using P-rich fertilizers or dung and (7) oxidation, precipitation and dilution in surface waters. Comparison of the results with experimental adsorption studies and speciation diagrams from the literature allows us to reconstruct and identify the typical behavior of arsenic in a natural environment under temperate climatic conditions. The main parameters identified are: (a) once liberated from the primary minerals, sorption processes on Fe-oxy-hydroxides dominate over Al-phases, such as Al-hydroxides or clay minerals and limit the As concentrations in the spring and well waters between 20 and 300 mug/l. (b) Precipitation as secondary minerals is limited to the weathering domain, where the As concentrations are still high and not yet too diluted by rain and soils waters. (c) Although neutral and alkaline pH conditions clearly increase the mobility of As, the main factor to mobilize As is a low redox potential (Eh close or below 0 mV), which favors the dissolution of the Fe-oxy-hydroxides on which the As is sorbed. (d) X-ray absorption spectroscopy (XAS) of As in water-logged humic forest soils indicates that the reduction to As III only occurs at the solid-water interface and that the solid contains As as As V (e) A and Bh horizons of humic cambisols can effectively capture As when As-rich waters flow through them. Complex spatial and temporal variation of the various parameters in a watershed results in repeated mobilization and immobilization of As, which continuously transports As from the upper to the lower part of a watershed and ultimately to the ocean. (C) 2004 Elsevier B.V. All rights reserved.

Comprendre la dispersion des espèces dans l'espace et dans le temps : un défi pour les biogéographes

De tout temps, l'homme s'est interrogé sur les facteurs qui régissent la distribution des espèces. Comment se dispersent-elles? Comment occupent-elles un milieu? Pourquoi certaines caractéristiques biologiques d'une même entité spécifique varient-elles en différents points de l'espace? Pour répondre à ces questions, il faut faire appel à de nombreuses disciplines scientifiques, telles que l'écologie, la géologie, la climatologie ou la paléontologie. Comprendre la distribution spatiale de la vie, quel que soit le niveau taxonomique, tel est l'objectif de la biogéographie

The stationary distribution of a continuously varying strategy in a class-structured population under mutation-selection-drift balance.

Many traits and/or strategies expressed by organisms are quantitative phenotypes. Because populations are of finite size and genomes are subject to mutations, these continuously varying phenotypes are under the joint pressure of mutation, natural selection and random genetic drift. This article derives the stationary distribution for such a phenotype under a mutation-selection-drift balance in a class-structured population allowing for demographically varying class sizes and/or changing environmental conditions. The salient feature of the stationary distribution is that it can be entirely characterized in terms of the average size of the gene pool and Hamilton's inclusive fitness effect. The exploration of the phenotypic space varies exponentially with the cumulative inclusive fitness effect over state space, which determines an adaptive landscape. The peaks of the landscapes are those phenotypes that are candidate evolutionary stable strategies and can be determined by standard phenotypic selection gradient methods (e.g. evolutionary game theory, kin selection theory, adaptive dynamics). The curvature of the stationary distribution provides a measure of the stability by convergence of candidate evolutionary stable strategies, and it is evaluated explicitly for two biological scenarios: first, a coordination game, which illustrates that, for a multipeaked adaptive landscape, stochastically stable strategies can be singled out by letting the size of the gene pool grow large; second, a sex-allocation game for diploids and haplo-diploids, which suggests that the equilibrium sex ratio follows a Beta distribution with parameters depending on the features of the genetic system.

Predictors of developmental dyslexia in European orthographies with varying complexity.

Background:  The relationship between phoneme awareness, rapid automatized naming (RAN), verbal short-term/working memory (ST/WM) and diagnostic category is investigated in control and dyslexic children, and the extent to which this depends on orthographic complexity. Methods:  General cognitive, phonological and literacy skills were tested in 1,138 control and 1,114 dyslexic children speaking six different languages spanning a large range of orthographic complexity (Finnish, Hungarian, German, Dutch, French, English). Results:  Phoneme deletion and RAN were strong concurrent predictors of developmental dyslexia, while verbal ST/WM and general verbal abilities played a comparatively minor role. In logistic regression models, more participants were classified correctly when orthography was more complex. The impact of phoneme deletion and RAN-digits was stronger in complex than in less complex orthographies. Conclusions:  Findings are largely consistent with the literature on predictors of dyslexia and literacy skills, while uniquely demonstrating how orthographic complexity exacerbates some symptoms of dyslexia.

Accurate calculations of intermolecular interaction energies using explicitly correlated coupled cluster wave functions and a dispersion-weighted MP2 method.

Explicitly correlated coupled-cluster calculations of intermolecular interaction energies for the S22 benchmark set of Jurecka, Sponer, Cerny, and Hobza (Chem. Phys. Phys. Chem. 2006, 8, 1985) are presented. Results obtained with the recently proposed CCSD(T)-F12a method and augmented double-zeta basis sets are found to be in very close agreement with basis set extrapolated conventional CCSD(T) results. Furthermore, we propose a dispersion-weighted MP2 (DW-MP2) approximation that combines the good accuracy of MP2 for complexes with predominately electrostatic bonding and SCS-MP2 for dispersion-dominated ones. The MP2-F12 and SCS-MP2-F12 correlation energies are weighted by a switching function that depends on the relative HF and correlation contributions to the interaction energy. For the S22 set, this yields a mean absolute deviation of 0.2 kcal/mol from the CCSD(T)-F12a results. The method, which allows obtaining accurate results at low cost, is also tested for a number of dimers that are not in the training set.

Phosphorous forms in cultivated indian black earth (anthrosols) of varying texture in the brazilian Amazon

Despite the agricultural importance of Indian Black Earth (IBE) in the Amazon region, there are few studies that report on the relation between soil texture and chemical fertility of IBE. These soils of pre-Colombian origin, with high contents of P, Ca and other nutrients are found across the Amazon valley. IBE profiles were studied to evaluate the total contents of P, its primary chemical forms and the P transformation phases in areas with IBE soils of variable texture and in adjacent reference soils. The soil texture strongly influenced soil fertility, changing in terms of transformation of the primary P forms and, consequently, predominant P forms in IBE. Soils with texture varying between clay and heavy clay had higher total P contents and primary Ca-P forms. Highest P-Al and lowest total P amounts were observed at the site Rio Preto da Eva, where texture varies from sandy loam to sandy clay loam. In the IBE with clay texture the amounts of soluble P, extracted with NH4Cl were highest, although different from Mehlich 1-extractable amounts.

Bulk-plasmon dispersion relations in metals

Using an extended-random-phase-approximation sum-rule technique, we have investigated the bulk-plasmon dispersion relation, incorporating in a simple way exchange and correlation effects within the jellium model. The results obtained are compared with recent experimental findings. The key role played by exchange and correlation effects in improving the agreement between theory and experiment is stressed. The static polarizability has also been calculated as a function of q. The formulas can be easily modified to incorporate band-structure effects (through an intraband electron effective mass) and core-polarization effects (through a static dielectric constant).