Quantitative integration of geophysical and hydraulic data : from the local towards the regional scale range

Des progrès significatifs ont été réalisés dans le domaine de l'intégration quantitative des données géophysique et hydrologique l'échelle locale. Cependant, l'extension à de plus grandes échelles des approches correspondantes constitue encore un défi majeur. Il est néanmoins extrêmement important de relever ce défi pour développer des modèles fiables de flux des eaux souterraines et de transport de contaminant. Pour résoudre ce problème, j'ai développé une technique d'intégration des données hydrogéophysiques basée sur une procédure bayésienne de simulation séquentielle en deux étapes. Cette procédure vise des problèmes à plus grande échelle. L'objectif est de simuler la distribution d'un paramètre hydraulique cible à partir, d'une part, de mesures d'un paramètre géophysique pertinent qui couvrent l'espace de manière exhaustive, mais avec une faible résolution (spatiale) et, d'autre part, de mesures locales de très haute résolution des mêmes paramètres géophysique et hydraulique. Pour cela, mon algorithme lie dans un premier temps les données géophysiques de faible et de haute résolution à travers une procédure de réduction déchelle. Les données géophysiques régionales réduites sont ensuite reliées au champ du paramètre hydraulique à haute résolution. J'illustre d'abord l'application de cette nouvelle approche dintégration des données à une base de données synthétiques réaliste. Celle-ci est constituée de mesures de conductivité hydraulique et électrique de haute résolution réalisées dans les mêmes forages ainsi que destimations des conductivités électriques obtenues à partir de mesures de tomographic de résistivité électrique (ERT) sur l'ensemble de l'espace. Ces dernières mesures ont une faible résolution spatiale. La viabilité globale de cette méthode est testée en effectuant les simulations de flux et de transport au travers du modèle original du champ de conductivité hydraulique ainsi que du modèle simulé. Les simulations sont alors comparées. Les résultats obtenus indiquent que la procédure dintégration des données proposée permet d'obtenir des estimations de la conductivité en adéquation avec la structure à grande échelle ainsi que des predictions fiables des caractéristiques de transports sur des distances de moyenne à grande échelle. Les résultats correspondant au scénario de terrain indiquent que l'approche d'intégration des données nouvellement mise au point est capable d'appréhender correctement les hétérogénéitées à petite échelle aussi bien que les tendances à gande échelle du champ hydraulique prévalent. Les résultats montrent également une flexibilté remarquable et une robustesse de cette nouvelle approche dintégration des données. De ce fait, elle est susceptible d'être appliquée à un large éventail de données géophysiques et hydrologiques, à toutes les gammes déchelles. Dans la deuxième partie de ma thèse, j'évalue en détail la viabilité du réechantillonnage geostatique séquentiel comme mécanisme de proposition pour les méthodes Markov Chain Monte Carlo (MCMC) appliquées à des probmes inverses géophysiques et hydrologiques de grande dimension . L'objectif est de permettre une quantification plus précise et plus réaliste des incertitudes associées aux modèles obtenus. En considérant une série dexemples de tomographic radar puits à puits, j'étudie deux classes de stratégies de rééchantillonnage spatial en considérant leur habilité à générer efficacement et précisément des réalisations de la distribution postérieure bayésienne. Les résultats obtenus montrent que, malgré sa popularité, le réechantillonnage séquentiel est plutôt inefficace à générer des échantillons postérieurs indépendants pour des études de cas synthétiques réalistes, notamment pour le cas assez communs et importants où il existe de fortes corrélations spatiales entre le modèle et les paramètres. Pour résoudre ce problème, j'ai développé un nouvelle approche de perturbation basée sur une déformation progressive. Cette approche est flexible en ce qui concerne le nombre de paramètres du modèle et lintensité de la perturbation. Par rapport au rééchantillonage séquentiel, cette nouvelle approche s'avère être très efficace pour diminuer le nombre requis d'itérations pour générer des échantillons indépendants à partir de la distribution postérieure bayésienne. - Significant progress has been made with regard to the quantitative integration of geophysical and hydrological data at the local scale. However, extending corresponding approaches beyond the local scale still represents a major challenge, yet is critically important for the development of reliable groundwater flow and contaminant transport models. To address this issue, I have developed a hydrogeophysical data integration technique based on a two-step Bayesian sequential simulation procedure that is specifically targeted towards larger-scale problems. The objective is to simulate the distribution of a target hydraulic parameter based on spatially exhaustive, but poorly resolved, measurements of a pertinent geophysical parameter and locally highly resolved, but spatially sparse, measurements of the considered geophysical and hydraulic parameters. To this end, my algorithm links the low- and high-resolution geophysical data via a downscaling procedure before relating the downscaled regional-scale geophysical data to the high-resolution hydraulic parameter field. I first illustrate the application of this novel data integration approach to a realistic synthetic database consisting of collocated high-resolution borehole measurements of the hydraulic and electrical conductivities and spatially exhaustive, low-resolution electrical conductivity estimates obtained from electrical resistivity tomography (ERT). The overall viability of this method is tested and verified by performing and comparing flow and transport simulations through the original and simulated hydraulic conductivity fields. The corresponding results indicate that the proposed data integration procedure does indeed allow for obtaining faithful estimates of the larger-scale hydraulic conductivity structure and reliable predictions of the transport characteristics over medium- to regional-scale distances. The approach is then applied to a corresponding field scenario consisting of collocated high- resolution measurements of the electrical conductivity, as measured using a cone penetrometer testing (CPT) system, and the hydraulic conductivity, as estimated from electromagnetic flowmeter and slug test measurements, in combination with spatially exhaustive low-resolution electrical conductivity estimates obtained from surface-based electrical resistivity tomography (ERT). The corresponding results indicate that the newly developed data integration approach is indeed capable of adequately capturing both the small-scale heterogeneity as well as the larger-scale trend of the prevailing hydraulic conductivity field. The results also indicate that this novel data integration approach is remarkably flexible and robust and hence can be expected to be applicable to a wide range of geophysical and hydrological data at all scale ranges. In the second part of my thesis, I evaluate in detail the viability of sequential geostatistical resampling as a proposal mechanism for Markov Chain Monte Carlo (MCMC) methods applied to high-dimensional geophysical and hydrological inverse problems in order to allow for a more accurate and realistic quantification of the uncertainty associated with the thus inferred models. Focusing on a series of pertinent crosshole georadar tomographic examples, I investigated two classes of geostatistical resampling strategies with regard to their ability to efficiently and accurately generate independent realizations from the Bayesian posterior distribution. The corresponding results indicate that, despite its popularity, sequential resampling is rather inefficient at drawing independent posterior samples for realistic synthetic case studies, notably for the practically common and important scenario of pronounced spatial correlation between model parameters. To address this issue, I have developed a new gradual-deformation-based perturbation approach, which is flexible with regard to the number of model parameters as well as the perturbation strength. Compared to sequential resampling, this newly proposed approach was proven to be highly effective in decreasing the number of iterations required for drawing independent samples from the Bayesian posterior distribution.

Lack of Correlation between Seminal and Plasma HIV-1 Viral Loads is Associated with CD4T Cell Depletion in Therapy-naïve HIV-1+ Patients

The present study was conducted to investigate a possible correlation between plasma (PVL) and seminal viral load (SVL) on treatment-naïve HIV-1-infected patients in Vitória, ES, Brazil. We also evaluated whether the progressive immunosuppression associated with HIV disease (as evidenced by declining CD4 T cell counts) has any impact on the correlation between PVL and SVL HIV-1. Viral load on paired blood and semen samples from 56 consecutive treatment-naïve patients were evaluated and compared to CD4 cell counts. Viral load and T cell counts (cells/µl) were determined by NASBA and by flow cytometry, respectively. Overall, a strong positive correlation between PVL and SVL (rho = 0.438, p = 0.001) was observed. However, when patients were grouped according to their CD4 counts, this correlation was only significant among patients with CD4 counts > 200 cells/µl. Results presented here demonstrate the existence of a strong correlation between PVL and SVL on patients with CD4 cell counts > 200 cells/µl, suggesting that this association may correlate with disease progression.

Correlation of male genital filaments and female spermathecal ducts in New World sand flies of the Lutzomyia intermedia species complex (Diptera: Psychodidae, Phlebotominae)

The lengths of the male genital filaments and female spermathecal ducts were measured in phlebotomine sand flies of the Lutzomyia intermedia species complex and the ratios between these characters calculated. Ratios for L. intermedia s. s. from Northeast vs Southeast Brazil (Espírito Santo and Minas Gerais), Espírito Santo/Minas Gerais vs Rio de Janeiro/São Paulo and L. intermedia vs L. neivai were significantly different at P < 0.1, 0.05 and 0.01 respectively when compared using ANOVA. The spermathecal ducts and genital filaments of L. intermedia were significantly longer than those of L. neivai (P < 0.01) and could be used to differentiate these species. The taxonomic and biological significance of these differences is discussed.

Antibiotic-dependent correlation between drug-induced killing and loss of luminescence in Streptococcus gordonii expressing luciferase.

Measuring antibiotic-induced killing relies on time-consuming biological tests. The firefly luciferase gene (luc) was successfully used as a reporter gene to assess antibiotic efficacy rapidly in slow-growing Mycobacterium tuberculosis. We tested whether luc expression could also provide a rapid evaluation of bactericidal drugs in Streptococcus gordonii. The suicide vectors pFW5luc and a modified version of pJDC9 carrying a promoterless luc gene were used to construct transcriptional-fusion mutants. One mutant susceptible to penicillin-induced killing (LMI2) and three penicillin-tolerant derivatives (LMI103, LMI104, and LMI105) producing luciferase under independent streptococcal promoters were tested. The correlation between antibiotic-induced killing and luminescence was determined with mechanistically unrelated drugs. Chloramphenicol (20 times the MIC) inhibited bacterial growth. In parallel, luciferase stopped increasing and remained stable, as determined by luminescence and Western blots. Ciprofloxacin (200 times the MIC) rapidly killed 1.5 log10 CFU/ml in 2-4 hr. Luminescence decreased simultaneously by 10-fold. In contrast, penicillin (200 times the MIC) gave discordant results. Although killing was slow (&lt; or = 0.5 log10 CFU/ml in 2 hr), luminescence dropped abruptly by 50-100-times in the same time. Inactivating penicillin with penicillinase restored luminescence, irrespective of viable counts. This was not due to altered luciferase expression or stability, suggesting some kind of post-translational modification. Luciferase shares homology with aminoacyl-tRNA synthetase and acyl-CoA ligase, which might be regulated by macromolecule synthesis and hence affected in penicillin-inhibited cells. Because of resemblance, luciferase might be down-regulated simultaneously. Luminescence cannot be universally used to predict antibiotic-induced killing. Thus, introducing reporter enzymes sharing mechanistic similarities with normal metabolic reactions might reveal other effects than those expected.

Telomerase activity and human telomerase reverse transcriptase mRNA expression in soft tissue tumors: correlation with grade, histology, and proliferative activity.

Telomerase activity (TA) is detected in most human cancers but, with few exceptions, not in normal somatic cells. Little is known about TA in soft tissue tumors. We have examined a series of benign and malignant soft tissue tumors for TA using the telomerase repeat amplification protocol assay. Analysis of the expression of the human telomerase reverse transcriptase was also carried out using RT-PCR. TA was undetectable in benign lesions (15 of 15) and low-grade sarcomas (6 of 6) and was detectable in 50% (19 of 38) of intermediate-/high-grade sarcomas. Although the presence of TA in soft tissue tumors is synonymous with malignancy, it is neither a reliable method in making the distinction between reactive/benign and malignant (especially low-grade) lesions nor a reliable marker of tumor aggressiveness. Leiomyosarcomas and storiform/pleomorphic malignant fibrous histiocytomas rarely showed TA, irrespective of their grade. A strong correlation between human telomerase reverse transcriptase mRNA expression and TA was observed, supporting the close relationship between both parameters. No significant relationship was observed between proliferative activity (as assessed by MIB-1 immunolabeling) and TA. We verified that the absence of telomerase expression was not due to the presence of telomerase inhibitors and therefore alternative mechanism(s) for cell immortalization, yet to be determined, seem to be involved in the development and/or maintenance of some soft tissue sarcomas.

CD30 in PTCLS: correlation between RNA and protein levels in a large european series from the LYSA

Introduction: CD22 is expressed on most B-non-Hodgkin lymphomas (NHL); inotuzumab ozogamicin (INO) is an anti-CD22 antibody conjugated to calicheamicin. This study evaluated the safety and tolerability of INO plus R-CVP in patients (pts) with relapsed/refractory CD22+ B-NHL. Efficacy data were also collected. Methods: Part 1 of this open-label study identified a maximum tolerated dose (MTD) of INO 0.8mg/m,2 on day 2 plus R-CVP (rituximab 375mg/m,2 cyclophosphamide 750mg/m,2 and vincristine 1.4mg/m,2 on day 1; prednisone 40mg/m,2 on days 1-5) every 21 days. Subsequently, pts were enrolled in the MTD confirmation cohort (part 2, n = 10), which required a dose-limiting toxicity rate of <33% in cycle 1 and <4 pts discontinuing prior to cycle 3 due to an adverse event (AE) in the MTD expansion cohort (part 3, n = 22), which explored preliminary activity. Results: Parts 2 and 3 enrolled 32 pts: 16 pts with diffuse large B-cell lymphoma, 15 with follicular lymphoma and one with mantle cell lymphoma. Median age was 64.5 years (range 44-81 years); 34% of pts had 1 prior regimen, 34% had 2, 28% had ≥3 and 3% had none (median 2; range 0-6).Median treatment duration was five cycles (range 1-6). Part 2 confirmed the MTD as standard dose R-CVP plus INO 0.8mg/m,2; 2/10 pts had a dose-limiting toxicity (grade 3 increased ALT/AST, grade 4 neutropenia requiring G-CSF). One pt discontinued because of an AE prior to cycle 3. Common treatment-related AEs were thrombocytopenia (78%), neutropenia (66%), fatigue (50%), leukopenia (50%), nausea (41%) and lymphopenia (38%); common grade 3/4 AEs were neutropenia (63%), thrombocytopenia (53%), leukopenia (38%) and lymphopenia (31%). There was one case of treatment-related fatal pneumonia with grade 4 neutropenia. Ten pts discontinued treatment due to AEs; thrombocytopenia/delayed platelet recovery was the leading cause (grade 1/2, n = 6; grade 3/4, n = 3). Objective response rate (ORR) was 77% (n = 24/31 evaluable pts), including 26% (n=8/31) with complete response (CR); three pts had stable disease. Of the pts with follicular lymphoma, ORR was 100% (n = 15/15), including seven pts with CR. Of the pts with diffuse large B-cell lymphoma, ORR was 60% (n = 9/16), including one pt with CR. Conclusions: Results suggest that INOplus R-CVP has acceptable toxicity and promising activity in relapsed/refractory CD22+ B-NHL. The most common grade 3/4 AEs were hematologic. Follow-up for progression-free and overall survival is ongoing.

Measuring an effect size from dichotomized data : contrasted results whether using a correlation or an odds-ratio

It is well known that dichotomizing continuous data has the effect to decrease statistical power when the goal is to test for a statistical association between two variables. Modern researchers however are focusing not only on statistical significance but also on an estimation of the "effect size" (i.e., the strength of association between the variables) to judge whether a significant association is also clinically relevant. In this article, we are interested in the consequences of dichotomizing continuous data on the value of an effect size in some classical settings. It turns out that the conclusions will not be the same whether using a correlation or an odds ratio to summarize the strength of association between the variables: Whereas the value of a correlation is typically decreased by a factor pi/2 after each dichotomization, the value of an odds ratio is at the same time raised to the power 2. From a descriptive statistical point of view, it is thus not clear whether dichotomizing continuous data leads to a decrease or to an increase in the effect size, as illustrated using a data set to investigate the relationship between motor and intellectual functions in children and adolescents

Trypanosoma cruzi: blood parasitism kinetics and their correlation with heart parasitism intensity during long-term infection of Beagle dogs

The goals of the present study were to evaluate the kinetics of blood parasitism by examination of fresh blood, blood culture (BC) and PCR assays and their correlation with heart parasitism during two years of infection in Beagle dogs inoculated with the Be-78, Y and ABC Trypanosoma cruzi strains. Our results showed that the parasite or its kDNA is easily detected during the acute phase in all infected animals. On the other hand, a reduced number of positive tests were verified during the chronic phase of the infection. The frequency of positive tests was correlated with T. cruzi strain. The percentage of positive BC and blood PCR performed in samples from animals inoculated with Be-78 and ABC strains were similar and significantly larger in relation to animals infected with the Y strain.Comparison of the positivity of PCR tests performed using blood and heart tissue samples obtained two years after infection showed two different patterns associated with the inoculated T. cruzi strain: (1) high PCR positivity for both blood and tissue was observed in animals infected with Be-78 or ABC strains; (2) lower and higher PCR positivity for the blood and tissue, respectively, was detected in animals infected with Y strains. These data suggest that the sensitivity of BC and blood PCR was T. cruzi strain dependent and, in contrast, the heart tissue PCR revealed higher sensitivity regardless of the parasite stock.

Early and late skin reactions to radiotherapy for breast cancer and their correlation with radiation-induced DNA damage in lymphocytes

INTRODUCTION Radiotherapy outcomes might be further improved by a greater understanding of the individual variations in normal tissue reactions that determine tolerance. Most published studies on radiation toxicity have been performed retrospectively. Our prospective study was launched in 1996 to measure the in vitro radiosensitivity of peripheral blood lymphocytes before treatment with radical radiotherapy in patients with breast cancer, and to assess the early and the late radiation skin side effects in the same group of patients. We prospectively recruited consecutive breast cancer patients receiving radiation therapy after breast surgery. To evaluate whether early and late side effects of radiotherapy can be predicted by the assay, a study was conducted of the association between the results of in vitro radiosensitivity tests and acute and late adverse radiation effects. METHODS Intrinsic molecular radiosensitivity was measured by using an initial radiation-induced DNA damage assay on lymphocytes obtained from breast cancer patients before radiotherapy. Acute reactions were assessed in 108 of these patients on the last treatment day. Late morbidity was assessed after 7 years of follow-up in some of these patients. The Radiation Therapy Oncology Group (RTOG) morbidity score system was used for both assessments. RESULTS Radiosensitivity values obtained using the in vitro test showed no relation with the acute or late adverse skin reactions observed. There was no evidence of a relation between acute and late normal tissue reactions assessed in the same patients. A positive relation was found between the treatment volume and both early and late side effects. CONCLUSION After radiation treatment, a number of cells containing major changes can have a long survival and disappear very slowly, becoming a chronic focus of immunological system stimulation. This stimulation can produce, in a stochastic manner, late radiation-related adverse effects of varying severity. Further research is warranted to identify the major determinants of normal tissue radiation response to make it possible to individualize treatments and improve the outcome of radiotherapy in cancer patients.

An iterative Multiscale Finite-Volume algorithm converging to the exact solution

The multiscale finite volume (MsFV) method has been developed to efficiently solve large heterogeneous problems (elliptic or parabolic); it is usually employed for pressure equations and delivers conservative flux fields to be used in transport problems. The method essentially relies on the hypothesis that the (fine-scale) problem can be reasonably described by a set of local solutions coupled by a conservative global (coarse-scale) problem. In most cases, the boundary conditions assigned for the local problems are satisfactory and the approximate conservative fluxes provided by the method are accurate. In numerically challenging cases, however, a more accurate localization is required to obtain a good approximation of the fine-scale solution. In this paper we develop a procedure to iteratively improve the boundary conditions of the local problems. The algorithm relies on the data structure of the MsFV method and employs a Krylov-subspace projection method to obtain an unconditionally stable scheme and accelerate convergence. Two variants are considered: in the first, only the MsFV operator is used; in the second, the MsFV operator is combined in a two-step method with an operator derived from the problem solved to construct the conservative flux field. The resulting iterative MsFV algorithms allow arbitrary reduction of the solution error without compromising the construction of a conservative flux field, which is guaranteed at any iteration. Since it converges to the exact solution, the method can be regarded as a linear solver. In this context, the schemes proposed here can be viewed as preconditioned versions of the Generalized Minimal Residual method (GMRES), with a very peculiar characteristic that the residual on the coarse grid is zero at any iteration (thus conservative fluxes can be obtained).

An approach to vision-based station keeping for an unmanned underwater vehicle

This paper presents an automatic vision-based system for UUV station keeping. The vehicle is equipped with a down-looking camera, which provides images of the sea-floor. The station keeping system is based on a feature-based motion detection algorithm, which exploits standard correlation and explicit textural analysis to solve the correspondence problem. A visual map of the area surveyed by the vehicle is constructed to increase the flexibility of the system, allowing the vehicle to position itself when it has lost the reference image. The testing platform is the URIS underwater vehicle. Experimental results demonstrating the behavior of the system on a real environment are presented

Correlation of biological serum markers with the degree of hepatic fibrosis and necroinflammatory activity in hepatitis C and schistosomiasis patients

Liver biopsy is the gold-standard method to stage fibrosis; however, it is an invasive procedure and is potentially dangerous. The main objective of this study was to evaluate biological markers, such as cytokines IL-13, IFN-γ, TNF-α and TGF-β, platelets, bilirubins (Bil), alanine aminotransferase (ALT) and aspartate aminotransferase (AST), total proteins, γ-glutamil transferase (γ-GT) and alkaline phosphatase (AP), that could be used to predict the severity of hepatic fibrosis in schistosomiasis and hepatitis C (HC) as isolated diseases or co-infections. The following patient groups were selected: HC (n = 39), HC/hepatosplenic schistosomiasis (HSS) (n = 19), HSS (n = 22) and a control group (n = 13). ANOVA and ROC curves were used for statistical analysis. P < 0.05 was considered significant. With HC patients we showed that TNF-α (p = 0.020) and AP (p = 0.005) could differentiate mild and severe fibrosis. With regard to necroinflammatory activity, AST (p = 0.002), γ-GT (p = 0.034) and AP (p = 0.001) were the best markers to differentiate mild and severe activity. In HC + HSS patients, total Bil (p = 0.008) was capable of differentiating between mild and severe fibrosis. In conclusion, our study was able to suggest biological markers that are non-invasive candidates to evaluate fibrosis and necroinflammatory activity in HC and HC + HSS.

Speculative parallelization of multipath radiosity algorithm

In computer graphics, global illumination algorithms take into account not only the light that comes directly from the sources, but also the light interreflections. This kind of algorithms produce very realistic images, but at a high computational cost, especially when dealing with complex environments. Parallel computation has been successfully applied to such algorithms in order to make it possible to compute highly-realistic images in a reasonable time. We introduce here a speculation-based parallel solution for a global illumination algorithm in the context of radiosity, in which we have taken advantage of the hierarchical nature of such an algorithm