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.

Improved multiexponential transient spectroscopy by iterative deconvolution

The analysis of multiexponential decays is challenging because of their complex nature. When analyzing these signals, not only the parameters, but also the orders of the models, have to be estimated. We present an improved spectroscopic technique specially suited for this purpose. The proposed algorithm combines an iterative linear filter with an iterative deconvolution method. A thorough analysis of the noise effect is presented. The performance is tested with synthetic and experimental data.

Discrete analysis of dividend payments in a non-life insurance portfolio

The process of free reserves in a non-life insurance portfolio as defined in the classical model of risk theory is modified by the introduction of dividend policies that set maximum levels for the accumulation of reserves. The first part of the work formulates the quantification of the dividend payments via the expectation of their current value under diferent hypotheses. The second part presents a solution based on a system of linear equations for discrete dividend payments in the case of a constant dividend barrier, illustrated by solving a specific case.

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.

Discrete contribution to psi ' -> J/psi plus gamma gamma

HadronPhysics2 (Grant Agreement No. 227431) (EU)

Dynamical properties of discrete breathers in curved chains with first and second neighbor interaction

We present the study of discrete breather dynamics in curved polymerlike chains consisting of masses connected via nonlinear springs. The polymer chains are one dimensional but not rectilinear and their motion takes place on a plane. After constructing breathers following numerically accurate procedures, we launch them in the chains and investigate properties of their propagation dynamics. We find that breather motion is strongly affected by the presence of curved regions of polymers, while the breathers themselves show a very strong resilience and remarkable stability in the presence of geometrical changes. For chains with strong angular rigidity we find that breathers either pass through bent regions or get reflected while retaining their frequency. Their motion is practically lossless and seems to be determined through local energy conservation. For less rigid chains modeled via second neighbor interactions, we find similarly that chain geometry typically does not destroy the localized breather states but, contrary to the angularly rigid chains, it induces some small but constant energy loss. Furthermore, we find that a curved segment acts as an active gate reflecting or refracting the incident breather and transforming its velocity to a value that depends on the discrete breathers frequency. We analyze the physical reasoning behind these seemingly general breather properties.

Interaction of a discrete breather with a lattice junction

We study the scattering of a moving discrete breather (DB) on a junction in a Fermi-Pasta-Ulam chain consisting of two segments with different masses of the particles. We consider four distinct cases: (i) a light-heavy (abrupt) junction in which the DB impinges on the junction from the segment with lighter mass, (ii) a heavy-light junction, (iii) an up mass ramp in which the mass in the heavier segment increases continuously as one moves away from the junction point, and (iv) a down mass ramp. Depending on the mass difference and DB characteristics (frequency and velocity), the DB can either reflect from, or transmit through, or get trapped at the junction or on the ramp. For the heavy-light junction, the DB can even split at the junction into a reflected and a transmitted DB. The latter is found to subsequently split into two or more DBs. For the down mass ramp the DB gets accelerated in several stages, with accompanying radiation (phonons). These results are rationalized by calculating the Peierls-Nabarro barrier for the various cases. We also point out implications of our results in realistic situations such as electron-phonon coupled chains.

Confinement of discrete breathers in inhomogeneously profiled nonlinear chains

We investigate numerically the scattering of a moving discrete breather on a pair of junctions in a Fermi-Pasta-Ulam chain. These junctions delimit an extended region with different masses of the particles. We consider (i) a rectangular trap, (ii) a wedge shaped trap, and (iii) a smoothly varying convex or concave mass profile. All three cases lead to DB confinement, with the ease of trapping depending on the profile of the trap. We also study the collision and trapping of two DBs within the profile as a function of trap width, shape, and approach time at the two junctions. The latter controls whether one or both DBs are trapped.

Multifractal wavelet filter of natural images

Natural images are characterized by the multiscaling properties of their contrast gradient, in addition to their power spectrum. In this Letter we show that those properties uniquely define an intrinsic wavelet and present a suitable technique to obtain it from an ensemble of images. Once this wavelet is known, images can be represented as expansions in the associated wavelet basis. The resulting code has the remarkable properties that it separates independent features at different resolution level, reducing the redundancy, and remains essentially unchanged under changes in the power spectrum. The possible generalization of this representation to other systems is discussed.

Airway epithelial IL-15 transforms monocytes into dendritic cells.

IL-15 has recently been shown to induce the differentiation of functional dendritic cells (DCs) from human peripheral blood monocytes. Since DCs lay in close proximity to epithelial cells in the airway mucosa, we investigated whether airway epithelial cells release IL-15 in response to inflammatory stimuli and thereby induce differentiation and maturation of DCs. Alveolar (A549) and bronchial (BEAS-2B) epithelial cells produced IL-15 spontaneously and in a time- and dose-dependent manner after stimulation with IL-1beta, IFN-gamma, or TNF-alpha. Airway epithelial cell supernatants induced an increase of IL-15Ralpha gene expression in ex vivo monocytes, and stimulated DCs enhanced their IL-15Ralpha gene expression up to 300-fold. Airway epithelial cell-conditioned media induced the differentiation of ex vivo monocytes into partially mature DCs (HLA-DR+, DC-SIGN+, CD14+, CD80-, CD83+, CD86+, CCR3+, CCR6(+), CCR7-). Based on their phenotypic (CD123+, BDCA2+, BDCA4+, BDCA1(-), CD1a-) and functional properties (limited maturation upon stimulation with LPS and limited capacity to induce T cell proliferation), these DCs resembled plasmacytoid DCs. The effects of airway epithelial cell supernatants were largely blocked by a neutralizing monoclonal antibody to IL-15. Thus, our results demonstrate that airway epithelial cell-conditioned media have the capacity to differentiate monocytes into functional DCs, a process substantially mediated by epithelial-derived IL-15.

Discrete analysis of dividend payments in a non-life insurance portfolio

The process of free reserves in a non-life insurance portfolio as defined in the classical model of risk theory is modified by the introduction of dividend policies that set maximum levels for the accumulation of reserves. The first part of the work formulates the quantification of the dividend payments via the expectation of their current value under diferent hypotheses. The second part presents a solution based on a system of linear equations for discrete dividend payments in the case of a constant dividend barrier, illustrated by solving a specific case.

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.

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.

Retaining principal components for discrete variables

The present study discusses retention criteria for principal components analysis (PCA) applied to Likert scale items typical in psychological questionnaires. The main aim is to recommend applied researchers to restrain from relying only on the eigenvalue-than-one criterion; alternative procedures are suggested for adjusting for sampling error. An additional objective is to add evidence on the consequences of applying this rule when PCA is used with discrete variables. The experimental conditions were studied by means of Monte Carlo sampling including several sample sizes, different number of variables and answer alternatives, and four non-normal distributions. The results suggest that even when all the items and thus the underlying dimensions are independent, eigenvalues greater than one are frequent and they can explain up to 80% of the variance in data, meeting the empirical criterion. The consequences of using Kaiser"s rule are illustrated with a clinical psychology example. The size of the eigenvalues resulted to be a function of the sample size and the number of variables, which is also the case for parallel analysis as previous research shows. To enhance the application of alternative criteria, an R package was developed for deciding the number of principal components to retain by means of confidence intervals constructed about the eigenvalues corresponding to lack of relationship between discrete variables.