Multi-modal Change Detection, Application to the Detection of Flooded Areas: Outcome of the 2009-2010 data fusion contest

The 2009-2010 Data Fusion Contest organized by the Data Fusion Technical Committee of the IEEE Geoscience and Remote Sensing Society was focused on the detection of flooded areas using multi-temporal and multi-modal images. Both high spatial resolution optical and synthetic aperture radar data were provided. The goal was not only to identify the best algorithms (in terms of accuracy), but also to investigate the further improvement derived from decision fusion. This paper presents the four awarded algorithms and the conclusions of the contest, investigating both supervised and unsupervised methods and the use of multi-modal data for flood detection. Interestingly, a simple unsupervised change detection method provided similar accuracy as supervised approaches, and a digital elevation model-based predictive method yielded a comparable projected change detection map without using post-event data.

Mass conservative three-dimensional water tracer distribution from Markov chain Monte Carlo inversion of time-lapse ground-penetrating radar data

Time-lapse geophysical measurements are widely used to monitor the movement of water and solutes through the subsurface. Yet commonly used deterministic least squares inversions typically suffer from relatively poor mass recovery, spread overestimation, and limited ability to appropriately estimate nonlinear model uncertainty. We describe herein a novel inversion methodology designed to reconstruct the three-dimensional distribution of a tracer anomaly from geophysical data and provide consistent uncertainty estimates using Markov chain Monte Carlo simulation. Posterior sampling is made tractable by using a lower-dimensional model space related both to the Legendre moments of the plume and to predefined morphological constraints. Benchmark results using cross-hole ground-penetrating radar travel times measurements during two synthetic water tracer application experiments involving increasingly complex plume geometries show that the proposed method not only conserves mass but also provides better estimates of plume morphology and posterior model uncertainty than deterministic inversion results.

Estimating vadose zone hydraulic properties using ground penetrating radar: The impact of prior information

A number of geophysical methods, such as ground-penetrating radar (GPR), have the potential to provide valuable information on hydrological properties in the unsaturated zone. In particular, the stochastic inversion of such data within a coupled geophysical-hydrological framework may allow for the effective estimation of vadose zone hydraulic parameters and their corresponding uncertainties. A critical issue in stochastic inversion is choosing prior parameter probability distributions from which potential model configurations are drawn and tested against observed data. A well chosen prior should reflect as honestly as possible the initial state of knowledge regarding the parameters and be neither overly specific nor too conservative. In a Bayesian context, combining the prior with available data yields a posterior state of knowledge about the parameters, which can then be used statistically for predictions and risk assessment. Here we investigate the influence of prior information regarding the van Genuchten-Mualem (VGM) parameters, which describe vadose zone hydraulic properties, on the stochastic inversion of crosshole GPR data collected under steady state, natural-loading conditions. We do this using a Bayesian Markov chain Monte Carlo (MCMC) inversion approach, considering first noninformative uniform prior distributions and then more informative priors derived from soil property databases. For the informative priors, we further explore the effect of including information regarding parameter correlation. Analysis of both synthetic and field data indicates that the geophysical data alone contain valuable information regarding the VGM parameters. However, significantly better results are obtained when we combine these data with a realistic, informative prior.

Bayesian Markov-chain-Monte-Carlo inversion of time-lapse cross hole ground-penetrating radar data to characterize the vadose zone at the Arrenaes field site, Denmark

The ground-penetrating radar (GPR) geophysical method has the potential to provide valuable information on the hydraulic properties of the vadose zone because of its strong sensitivity to soil water content. In particular, recent evidence has suggested that the stochastic inversion of crosshole GPR traveltime data can allow for a significant reduction in uncertainty regarding subsurface van Genuchten-Mualem (VGM) parameters. Much of the previous work on the stochastic estimation of VGM parameters from crosshole GPR data has considered the case of steady-state infiltration conditions, which represent only a small fraction of practically relevant scenarios. We explored in detail the dynamic infiltration case, specifically examining to what extent time-lapse crosshole GPR traveltimes, measured during a forced infiltration experiment at the Arreneas field site in Denmark, could help to quantify VGM parameters and their uncertainties in a layered medium, as well as the corresponding soil hydraulic properties. We used a Bayesian Markov-chain-Monte-Carlo inversion approach. We first explored the advantages and limitations of this approach with regard to a realistic synthetic example before applying it to field measurements. In our analysis, we also considered different degrees of prior information. Our findings indicate that the stochastic inversion of the time-lapse GPR data does indeed allow for a substantial refinement in the inferred posterior VGM parameter distributions compared with the corresponding priors, which in turn significantly improves knowledge of soil hydraulic properties. Overall, the results obtained clearly demonstrate the value of the information contained in time-lapse GPR data for characterizing vadose zone dynamics.

Native-like, long synthetic peptides as components of sub-unit vaccines: practical and theoretical considerations for their use in humans.

Vaccines have been used as a successful tool in medicine by way of controlling many major diseases. In spite of this, vaccines today represent only a handful of all infectious diseases. Therefore, there is a pressing demand for improvements of existing vaccines with particular reference to higher efficacy and undisputed safety profiles. To this effect, as an alternative to available vaccine technologies, there has been a drive to develop vaccine candidate polypeptides by chemical synthesis. In our laboratory, we have recently developed a technology to manufacture long synthetic peptides of up to 130 residues, which are correctly folded and biologically active. This paper discusses the advantages of the molecularly defined, long synthetic peptide approach in the context of vaccine design, development and use in human vaccination.

Malaria vaccine development using synthetic peptides as a technical platform.

The review covers the development of synthetic peptides as vaccine candidates for Plasmodium falciparum- and Plasmodium vivax-induced malaria from its beginning up to date and the concomitant progress of solid phase peptide synthesis (SPPS) that enables the production of long peptides in a routine fashion. The review also stresses the development of other complementary tools and actions in order to achieve the long sought goal of an efficacious malaria vaccine.

The antibody response in mice to carrier-free synthetic polymers of Plasmodium falciparum circumsporozoite repetitive epitope is I-Ab-restricted: possible implications for malaria vaccines.

The immunogenicity of a novel synthetic peptide consisting of an average of 40 (Asn-Ala-Asn-Pro) repeats of the circumsporozoite protein of Plasmodium falciparum, (NANP)40, was studied in mice without using any carrier proteins. First, high titers of anti-(NANP)40 antibodies could be obtained after immunization of C57BL/6 mice. These antibodies also reacted with an extract of mosquitoes infected with P. falciparum sporozoites. C57BL/6 nu/nu mice did not produce antibodies against (NANP)40. Secondly, when 14 strains of mice with nine different H-2 haplotypes were immunized with (NANP)40 without carrier, only H-2b mice were found to produce anti-(NANP)40 antibodies, whereas all non-H-2b mice were consistently unresponsive. This response was demonstrated to be I-A-linked by using recombinant and mutant mice. I-Ab [B10.A(5R)] mice produced anti-(NANP)40 antibodies as well as H-2b inbred mice. B6CH-2bm12 I-Ab-mutant mice showed only a very low response. Third, the antibody response against (NANP)40 could be induced in nonresponder mice by immunization with the peptide coupled to a carrier protein. In view of the existence of such an exceptional H-2b restriction in the response to sporozoite synthetic peptides in mice, the triggering of peptide-specific T cell responses in humans receiving sporozoite malaria vaccines might be difficult to achieve.

Attenuated poxviruses expressing a synthetic HIV protein stimulate HLA-A2-restricted cytotoxic T-cell responses.

Efficient HIV vaccines have to trigger cell-mediated immunity directed against various viral antigens. However little is known about the breadth of the response induced by vaccines carrying multiple proteins. Here, we report on the immunogenicity of a construct harbouring a fusion of the HIV-1 IIIB gag, pol and nef genes (gpn) designed for optimal safety and equimolar expression of the HIV proteins. The attenuated poxviruses, MVA and NYVAC, harbouring the gpn construct, induced potent immune responses in conventional mice characterised by stimulation of Gpn-specific IFN-gamma-producing cells and cytotoxic T cells. In HLA-A2 transgenic mice, recombinant MVA elicited cytotoxic responses against epitopes recognised in most HLA-A2+ HIV-1-infected individuals. We also found that the MVA vaccine triggered the in vitro expansion of peripheral blood cells isolated from a HIV-1-seropositive patient and with similar specificity as found in immunised HLA-A2 transgenic mice. In conclusion, the synthetic HIV polyantigen Gpn delivered by MVA is immunogenic, efficiently processed and presented by human MHC class I molecules.

Paclitaxel-eluting biodegradable synthetic vascular prostheses: a step towards reduction of neointima formation?

BACKGROUND: Clinical small-caliber vascular prostheses are unsatisfactory. Reasons for failure are early thrombosis and late intimal hyperplasia. We thus prepared biodegradable small-caliber vascular prostheses using electrospun polycaprolactone (PCL) with slow-releasing paclitaxel (PTX), an antiproliferative drug. METHODS AND RESULTS: PCL solutions containing PTX were used to prepare nonwoven nanofibre-based 2-mm ID prostheses. Mechanical morphological properties and drug loading, distribution, and release were studied in vitro. Infrarenal abdominal aortic replacement was carried out with nondrug-loaded and drug-loaded prostheses in 18 rats and followed for 6 months. Patency, stenosis, tissue reaction, and drug effect on endothelialization, vascular remodeling, and neointima formation were studied in vivo. In vitro prostheses showed controlled morphology mimicking extracellular matrix with mechanical properties similar to those of native vessels. PTX-loaded grafts with suitable mechanical properties and controlled drug-release were obtained by factorial design. In vivo, both groups showed 100% patency, no stenosis, and no aneurysmal dilatation. Endothelial coverage and cell ingrowth were significantly reduced at 3 weeks and delayed at 12 and 24 weeks in PTX grafts, but as envisioned, neointima formation was significantly reduced in these grafts at 12 weeks and delayed at 6 months. CONCLUSIONS: Biodegradable, electrospun, nanofibre, polycaprolactone prostheses are promising because in vitro they maintain their mechanical properties (regardless of PTX loading), and in vivo show good patency, reendothelialize, and remodel with autologous cells. PTX loading delays endothelialization and cellular ingrowth. Conversely, it reduces neointima formation until the end point of our study and thus may be an interesting option for small caliber vascular grafts.

Synthetic Plasmodium falciparum circumsporozoite peptides elicit heterogenous L3T4+ T cell proliferative responses in H-2b mice.

The ability of synthetic P. falciparum (NANP)n circumsporozoite peptides to elicit murine T cell proliferative responses was studied. When C57BL/6, C3H, and DBA/2 mice were injected with (NANP)40, only C57BL/6 (H-2b)-immune lymph node cells proliferated on restimulation in vitro with the same peptide. By using anti-I-A monoclonal antibodies or spleen cells from congenic H-2b mice as a source of antigen-presenting cells, the T cell proliferative response was shown to be restricted to the I-Ab region of the C57BL/6 haplotype. These results are in agreement with previous experiments which demonstrated that the anti-(NANP)40 antibody response was uniquely restricted to C57BL/6 (H-2b) mice. Several C57BL/6 long-term (NANP)n-specific T cell lines and clones were derived. All of the clones exhibited the L3T4 helper T cell phenotype. A considerable heterogeneity of T cell responses was observed when the lines and clones were stimulated with different concentrations of the various peptides studied. The results, together with the observed genetic restriction for both antibody and T cell responses, suggest that perhaps not all individuals who receive a similar repetitive tetrapeptide sporozoite malaria vaccine will develop T cell and or antibody responses.

Pharmacokinetics of synthetic atrial natriuretic peptides in normal men.

The kinetics of atrial natriuretic peptides (ANP) and the kinetic profile of their effect on blood pressure and renal hemodynamic and electrolyte excretion were investigated in 20 salt-loaded healthy volunteers during and after constant rate infusion. At steady state, mean plasma concentrations of ANP were measured at 210, 430, and 2990 pg/ml and mean systemic clearance was 2.6, 2.5, and 1.7 L/min for ANP infusion rates of 0.5, 1, and 5 micrograms/min, respectively, which corresponds to the clearance rate of other vasoactive peptide hormones. The apparent volume of distribution averaged 17 L and the mean half-life was 4.5 minutes. ANP induced dose-related effects on systemic and renal hemodynamic, as well as urinary electrolyte excretion, albeit with a time lag between onset and full effect.

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.