Estimation of the correlation structure of crustal velocity heterogeneity from seismic reflection data

Numerous sources of evidence point to the fact that heterogeneity within the Earth's deep crystalline crust is complex and hence may be best described through stochastic rather than deterministic approaches. As seismic reflection imaging arguably offers the best means of sampling deep crustal rocks in situ, much interest has been expressed in using such data to characterize the stochastic nature of crustal heterogeneity. Previous work on this problem has shown that the spatial statistics of seismic reflection data are indeed related to those of the underlying heterogeneous seismic velocity distribution. As of yet, however, the nature of this relationship has remained elusive due to the fact that most of the work was either strictly empirical or based on incorrect methodological approaches. Here, we introduce a conceptual model, based on the assumption of weak scattering, that allows us to quantitatively link the second-order statistics of a 2-D seismic velocity distribution with those of the corresponding processed and depth-migrated seismic reflection image. We then perform a sensitivity study in order to investigate what information regarding the stochastic model parameters describing crustal velocity heterogeneity might potentially be recovered from the statistics of a seismic reflection image using this model. Finally, we present a Monte Carlo inversion strategy to estimate these parameters and we show examples of its application at two different source frequencies and using two different sets of prior information. Our results indicate that the inverse problem is inherently non-unique and that many different combinations of the vertical and lateral correlation lengths describing the velocity heterogeneity can yield seismic images with the same 2-D autocorrelation structure. The ratio of all of these possible combinations of vertical and lateral correlation lengths, however, remains roughly constant which indicates that, without additional prior information, the aspect ratio is the only parameter describing the stochastic seismic velocity structure that can be reliably recovered.

Unraveling dynamics of human physical activity patterns in chronic pain conditions.

Chronic pain is a complex disabling experience that negatively affects the cognitive, affective and physical functions as well as behavior. Although the interaction between chronic pain and physical functioning is a well-accepted paradigm in clinical research, the understanding of how pain affects individuals' daily life behavior remains a challenging task. Here we develop a methodological framework allowing to objectively document disruptive pain related interferences on real-life physical activity. The results reveal that meaningful information is contained in the temporal dynamics of activity patterns and an analytical model based on the theory of bivariate point processes can be used to describe physical activity behavior. The model parameters capture the dynamic interdependence between periods and events and determine a 'signature' of activity pattern. The study is likely to contribute to the clinical understanding of complex pain/disease-related behaviors and establish a unified mathematical framework to quantify the complex dynamics of various human activities.

The four- vs. alternative six-factor structure of the French WISC-IV: Comparison using confirmatory factor analyses

Exploratory and confirmatory factor analyses reported in the French technical manual of the WISC-IV provides evidence supporting a structure with four indices: Verbal Comprehension (VCI), Perceptual Reasoning (PRI), Working Memory (WMI), and Processing Speed (PSI). Although the WISC-IV is more attuned to contemporary theory, it is still not in total accordance with the dominant theory: the Cattell-Horn-Carroll (CHC) theory of cognitive ability. This study was designed to determine whether the French WISC-IV is better described with the four-factor solution or whether an alternative model based on the CHC theory is more appropriate. The intercorrelations matrix reported in the French technical manual was submitted to confirmatory factor analysis. A comparison of competing models suggests that a model based on the CHC theory fits the data better than the current WISC-IV structure. It appears that the French WISC-IV in fact measures six factors: crystallized intelligence (Gc), fluid intelligence (Gf), short-term memory (Gsm), processing speed (Gs), quantitative knowledge (Gq), and visual processing (Gv). We recommend that clinicians interpret the subtests of the French WISC-IV in relation to this CHC model in addition to the four indices.

Syndrome d'hyperstimulation ovarienne [The syndrome of ovarian hyperstimulation].

The ovarian hyperstimulation syndrome (SHO) can be defined as an iatrogenic pathology induced by active substances administered for controlling follicular maturation and ovulation. The etiology, the physiopathology, the diagnostic and therapeutic methods available are discussed. A theoretical model, based on clinical data, allows identification of a set of criteria which should help determining prospectively the chances of development of such a pathology.

Estimation of the noisy component of anatomical backgrounds.

The knowledge of the relationship that links radiation dose and image quality is a prerequisite to any optimization of medical diagnostic radiology. Image quality depends, on the one hand, on the physical parameters such as contrast, resolution, and noise, and on the other hand, on characteristics of the observer that assesses the image. While the role of contrast and resolution is precisely defined and recognized, the influence of image noise is not yet fully understood. Its measurement is often based on imaging uniform test objects, even though real images contain anatomical backgrounds whose statistical nature is much different from test objects used to assess system noise. The goal of this study was to demonstrate the importance of variations in background anatomy by quantifying its effect on a series of detection tasks. Several types of mammographic backgrounds and signals were examined by psychophysical experiments in a two-alternative forced-choice detection task. According to hypotheses concerning the strategy used by the human observers, their signal to noise ratio was determined. This variable was also computed for a mathematical model based on the statistical decision theory. By comparing theoretical model and experimental results, the way that anatomical structure is perceived has been analyzed. Experiments showed that the observer's behavior was highly dependent upon both system noise and the anatomical background. The anatomy partly acts as a signal recognizable as such and partly as a pure noise that disturbs the detection process. This dual nature of the anatomy is quantified. It is shown that its effect varies according to its amplitude and the profile of the object being detected. The importance of the noisy part of the anatomy is, in some situations, much greater than the system noise. Hence, reducing the system noise by increasing the dose will not improve task performance. This observation indicates that the tradeoff between dose and image quality might be optimized by accepting a higher system noise. This could lead to a better resolution, more contrast, or less dose.

Dosage optimization of treatments using population pharmacokinetic modeling and simulation.

Pharmacokinetic variability in drug levels represent for some drugs a major determinant of treatment success, since sub-therapeutic concentrations might lead to toxic reactions, treatment discontinuation or inefficacy. This is true for most antiretroviral drugs, which exhibit high inter-patient variability in their pharmacokinetics that has been partially explained by some genetic and non-genetic factors. The population pharmacokinetic approach represents a very useful tool for the description of the dose-concentration relationship, the quantification of variability in the target population of patients and the identification of influencing factors. It can thus be used to make predictions and dosage adjustment optimization based on Bayesian therapeutic drug monitoring (TDM). This approach has been used to characterize the pharmacokinetics of nevirapine (NVP) in 137 HIV-positive patients followed within the frame of a TDM program. Among tested covariates, body weight, co-administration of a cytochrome (CYP) 3A4 inducer or boosted atazanavir as well as elevated aspartate transaminases showed an effect on NVP elimination. In addition, genetic polymorphism in the CYP2B6 was associated with reduced NVP clearance. Altogether, these factors could explain 26% in NVP variability. Model-based simulations were used to compare the adequacy of different dosage regimens in relation to the therapeutic target associated with treatment efficacy. In conclusion, the population approach is very useful to characterize the pharmacokinetic profile of drugs in a population of interest. The quantification and the identification of the sources of variability is a rational approach to making optimal dosage decision for certain drugs administered chronically.

Shielding requirements in helical tomotherapy.

Helical tomotherapy is a relatively new intensity-modulated radiation therapy (IMRT) treatment for which room shielding has to be reassessed for the following reasons. The beam-on-time needed to deliver a given target dose is increased and leads to a weekly workload of typically one order of magnitude higher than that for conventional radiation therapy. The special configuration of tomotherapy units does not allow the use of standard shielding calculation methods. A conventional linear accelerator must be shielded for primary, leakage and scatter photon radiations. For tomotherapy, primary radiation is no longer the main shielding issue since a beam stop is mounted on the gantry directly opposite the source. On the other hand, due to the longer irradiation time, the accelerator head leakage becomes a major concern. An analytical model based on geometric considerations has been developed to determine leakage radiation levels throughout the room for continuous gantry rotation. Compared to leakage radiation, scatter radiation is a minor contribution. Since tomotherapy units operate at a nominal energy of 6 MV, neutron production is negligible. This work proposes a synthetic and conservative model for calculating shielding requirements for the Hi-Art II TomoTherapy unit. Finally, the required concrete shielding thickness is given for different positions of interest.

Effects of alternative sets of climatic predictors on species distribution models and associated estimates of extinction risk: a test with plants in an arid environment

Aim To evaluate the effects of using distinct alternative sets of climatic predictor variables on the performance, spatial predictions and future projections of species distribution models (SDMs) for rare plants in an arid environment. . Location Atacama and Peruvian Deserts, South America (18º30'S - 31º30'S, 0 - 3 000 m) Methods We modelled the present and future potential distributions of 13 species of Heliotropium sect. Cochranea, a plant group with a centre of diversity in the Atacama Desert. We developed and applied a sequential procedure, starting from climate monthly variables, to derive six alternative sets of climatic predictor variables. We used them to fit models with eight modelling techniques within an ensemble forecasting framework, and derived climate change projections for each of them. We evaluated the effects of using these alternative sets of predictor variables on performance, spatial predictions and projections of SDMs using Generalised Linear Mixed Models (GLMM). Results The use of distinct sets of climatic predictor variables did not have a significant effect on overall metrics of model performance, but had significant effects on present and future spatial predictions. Main conclusion Using different sets of climatic predictors can yield the same model fits but different spatial predictions of current and future species distributions. This represents a new form of uncertainty in model-based estimates of extinction risk that may need to be better acknowledged and quantified in future SDM studies.

The estimation of PMP and PMF on alpine basins in Switzerland

This paper presents a very fine grid hydrological model based on the spatiotemporal repartition of precipitation and on the topography. The goal is to estimate the flood on a catchment area, using a Probable Maximum Precipitation (PMP) leading to a Probable Maximum Flood (PMF). The spatiotemporal distribution of the precipitation was realized using six clouds modeled by the advection-diffusion equation. The equation shows the movement of the clouds over the terrain and also gives the evolution of the rain intensity in time. This hydrological modeling is followed by a hydraulic modeling of the surface and subterranean flows, done considering the factors that contribute to the hydrological cycle, such as the infiltration, the exfiltration and the snowmelt. This model was applied to several Swiss basins using measured rain, with results showing a good correlation between the simulated and observed flows. This good correlation proves that the model is valid and gives us the confidence that the results can be extrapolated to phenomena of extreme rainfall of PMP type. In this article we present some results obtained using a PMP rainfall and the developed model.

Outcome evaluation of ankle osteoarthritis treatments using spatio-temporal gait parameters and plantar pressure during unconstrained long distance walking : P25

Introduction: Ankle arthrodesis (AD) and total ankle replacement (TAR) are typical treatments for ankle osteoarthritis (AO). Despite clinical interest, there is a lack of their outcome evaluation using objective criteria. Gait analysis and plantar pressure assessment are appropriate to detect pathologies in orthopaedics but they are mostly used in lab with few gait cycles. In this study, we propose an ambulatory device based on inertial and plantar pressure sensors to compare the gait during long-distance trials between healthy subjects (H) and patients with AO or treated by AD and TAR. Methods: Our study included four groups: 11 patients with AO, 9 treated by TAR, 7 treated by AD and 6 control subjects. An ambulatory system (Physilog®, CH) was used for gait analysis; plantar pressure measurements were done using a portable insole (Pedar®-X, DE). The subjects were asked to walk 50 meters in two trials. Mean value and coefficient of variation of spatio-temporal gait parameters were calculated for each trial. Pressure distribution was analyzed in ten subregions of foot. All parameters were compared among the four groups using multi-level model-based statistical analysis. Results: Significant difference (p <0.05) with control was noticed for AO patients in maximum force in medial hindfoot and forefoot and in central forefoot. These differences were no longer significant in TAR and AD groups. Cadence and speed of all pathologic groups showed significant difference with control. Both treatments showed a significant improvement in double support and stance. TAR decreased variability in speed, stride length and knee ROM. Conclusions: In spite of a small sample size, this study showed that ankle function after AO treatments can be evaluated objectively based on plantar pressure and spatio-temporal gait parameters measured during unconstrained walking outside the lab. The combination of these two ambulatory techniques provides a promising way to evaluate foot function in clinics.

Systemic exposure to tobramycin after local antibiotic treatment with calcium sulphate as carrier material.

INTRODUCTION: Osteoset(®) T is a calcium sulphate void filler containing 4% tobramycin sulphate, used to treat bone and soft tissue infections. Despite systemic exposure to the antibiotic, there are no pharmacokinetic studies in humans published so far. Based on the observations made in our patients, a model predicting tobramycin serum levels and evaluating their toxicity potential is presented. METHODS: Following implantation of Osteoset(®) T, tobramycin serum concentrations were monitored systematically. A pharmacokinetic analysis was performed using a non-linear mixed effects model based on a one compartment model with first-degree absorption. RESULTS: Data from 12 patients treated between October 2006 and March 2008 were analysed. Concentration profiles were consistent with the first-order slow release and single-compartment kinetics, whilst showing important variability. Predicted tobramycin serum concentrations depended clearly on both implanted drug amount and renal function. DISCUSSION AND CONCLUSION: Despite the popularity of aminoglycosides for local antibiotic therapy, pharmacokinetic data for this indication are scarce, and not available for calcium sulphate as carrier material. Systemic exposure to tobramycin after implantation of Osteoset(®) T appears reassuring regarding toxicity potential, except in case of markedly impaired renal function. We recommend in adapting the dosage to the estimated creatinine clearance rather than solely to the patient's weight.

Ink dating using thermal desorption and gas chromatography / mass spectrometry: comparison of results obtained in two laboratories

Recent ink dating methods focused mainly on changes in solvent amounts occurring over time. A promising method was developed at the Landeskriminalamt of Munich using thermal desorption (TD) followed by gas chromatography / mass spectrometry (GC/MS) analysis. Sequential extractions of the phenoxyethanol present in ballpoint pen ink entries were carried out at two different temperatures. This method is applied in forensic practice and is currently implemented in several laboratories participating to the InCID group (International Collaboration on Ink Dating). However, harmonization of the method between the laboratories proved to be a particularly sensitive and time consuming task. The main aim of this work was therefore to implement the TD-GC/MS method at the Bundeskriminalamt (Wiesbaden, Germany) in order to evaluate if results were comparable to those obtained in Munich. At first validation criteria such as limits of reliable measurements, linearity and repeatability were determined. Samples were prepared in three different laboratories using the same inks and analyzed using two TDS-GC/MS instruments (one in Munich and one in Wiesbaden). The inter- and intra-laboratory variability of the ageing parameter was determined and ageing curves were compared. While inks stored in similar conditions yielded comparable ageing curves, it was observed that significantly different storage conditions had an influence on the resulting ageing curves. Finally, interpretation models, such as thresholds and trend tests, were evaluated and discussed in view of the obtained results. Trend tests were considered more suitable than threshold models. As both approaches showed limitations, an alternative model, based on the slopes of the ageing curves, was also proposed.

Meiotic drive favors Robertsonian metacentric chromosomes in the common shrew (Sorex araneus, Insectivora, mammalia).

Meiotic drive has attracted much interest because it concerns the robustness of Mendelian segregation and its genetic and evolutionary stability. We studied chromosomal meiotic drive in the common shrew (Sorex araneus, Insectivora, Mammalia), which exhibits one of the most remarkable chromosomal polymorphisms within mammalian species. The open question of the evolutionary success of metacentric chromosomes (Robertsonian fusions) versus acrocentrics in the common shrew prompted us to test whether a segregation distortion in favor of metacentrics is present in female and/or male meiosis. Performing crosses under controlled laboratory conditions with animals from natural populations, we found a clear trend toward a segregation distortion in favor of metacentrics during male meiosis, two chromosome combinations (gm and jl) being significantly preferred over their acrocentric homologs. Apart for one Robertsonian fusion (hi), this trend was absent in female meiosis. We propose a model based on recombination events between twin acrocentrics to explain the difference in transmission ratios of the same metacentric in different sexes and unequal drive of particular metacentrics in the same sex. Pooled data for female and male meiosis revealed a trend toward stronger segregation distortion for larger metacentrics. This is partially in agreement with the frequency of metacentrics occurring in natural populations of a chromosome race showing a high degree of chromosomal polymorphism.

Crystal size distribution of periclase in contact metamorphic dolomite marbles from the southern Adamello Massif, Italy

Crystal size distributions (CSD) of periclase in contact metamorphic dolomite marbles are presented for two profiles near the Cima Uzza summit in the southern Adamello Massif (Italy). The database was combined with geochemical and petrological information to deduce the controls on the periclase-forming reaction. The contact metamorphic dolomite marbles are exposed at the contact of mafic intrusive rocks and are partially surrounded by them. Brucite is retrograde and pseudomorphs spherical periclase crystals. Prograde periclase growth is the consequence of limited infiltration of water-rich fluid at T near 605C. Stable isotope data show depletion in (13)C and (18)O over a narrow region (40 cm) near the magmatic contact, whereas the periclase-forming reaction front extends up to 4 m from the contact. CSD analyses along the two profiles show that the median grain size of the periclase crystals does not change, but that there is a progressively greater distribution of grain sizes, including a greater proportion of larger grains, with increasing distance from the contact. A qualitative model, based on the textural and geochemical data, attributes these variations in grain size to changing reaction affinities along a kinetically dispersed infiltration front. This study highlights the need to invoke disequilibrium processes for metamorphic mineral growth and expands the use of CSDs to systems of mineral formation driven by fluid infiltration.

Robust Nonlinear Mapping of Soil Contamination Using Support Regression

Spatial data analysis mapping and visualization is of great importance in various fields: environment, pollution, natural hazards and risks, epidemiology, spatial econometrics, etc. A basic task of spatial mapping is to make predictions based on some empirical data (measurements). A number of state-of-the-art methods can be used for the task: deterministic interpolations, methods of geostatistics: the family of kriging estimators (Deutsch and Journel, 1997), machine learning algorithms such as artificial neural networks (ANN) of different architectures, hybrid ANN-geostatistics models (Kanevski and Maignan, 2004; Kanevski et al., 1996), etc. All the methods mentioned above can be used for solving the problem of spatial data mapping. Environmental empirical data are always contaminated/corrupted by noise, and often with noise of unknown nature. That's one of the reasons why deterministic models can be inconsistent, since they treat the measurements as values of some unknown function that should be interpolated. Kriging estimators treat the measurements as the realization of some spatial randomn process. To obtain the estimation with kriging one has to model the spatial structure of the data: spatial correlation function or (semi-)variogram. This task can be complicated if there is not sufficient number of measurements and variogram is sensitive to outliers and extremes. ANN is a powerful tool, but it also suffers from the number of reasons. of a special type ? multiplayer perceptrons ? are often used as a detrending tool in hybrid (ANN+geostatistics) models (Kanevski and Maignank, 2004). Therefore, development and adaptation of the method that would be nonlinear and robust to noise in measurements, would deal with the small empirical datasets and which has solid mathematical background is of great importance. The present paper deals with such model, based on Statistical Learning Theory (SLT) - Support Vector Regression. SLT is a general mathematical framework devoted to the problem of estimation of the dependencies from empirical data (Hastie et al, 2004; Vapnik, 1998). SLT models for classification - Support Vector Machines - have shown good results on different machine learning tasks. The results of SVM classification of spatial data are also promising (Kanevski et al, 2002). The properties of SVM for regression - Support Vector Regression (SVR) are less studied. First results of the application of SVR for spatial mapping of physical quantities were obtained by the authorsin for mapping of medium porosity (Kanevski et al, 1999), and for mapping of radioactively contaminated territories (Kanevski and Canu, 2000). The present paper is devoted to further understanding of the properties of SVR model for spatial data analysis and mapping. Detailed description of the SVR theory can be found in (Cristianini and Shawe-Taylor, 2000; Smola, 1996) and basic equations for the nonlinear modeling are given in section 2. Section 3 discusses the application of SVR for spatial data mapping on the real case study - soil pollution by Cs137 radionuclide. Section 4 discusses the properties of the modelapplied to noised data or data with outliers.