MBIS: multivariate Bayesian image segmentation tool.

We present MBIS (Multivariate Bayesian Image Segmentation tool), a clustering tool based on the mixture of multivariate normal distributions model. MBIS supports multichannel bias field correction based on a B-spline model. A second methodological novelty is the inclusion of graph-cuts optimization for the stationary anisotropic hidden Markov random field model. Along with MBIS, we release an evaluation framework that contains three different experiments on multi-site data. We first validate the accuracy of segmentation and the estimated bias field for each channel. MBIS outperforms a widely used segmentation tool in a cross-comparison evaluation. The second experiment demonstrates the robustness of results on atlas-free segmentation of two image sets from scan-rescan protocols on 21 healthy subjects. Multivariate segmentation is more replicable than the monospectral counterpart on T1-weighted images. Finally, we provide a third experiment to illustrate how MBIS can be used in a large-scale study of tissue volume change with increasing age in 584 healthy subjects. This last result is meaningful as multivariate segmentation performs robustly without the need for prior knowledge.

Production of nanometric particles in radio frequency glow discharges in mixtures of silane and methane

The formation of silicon particles in rf glow discharges has attracted attention due to their effect as a contaminant during film deposition or etching. However, silicon and silicon alloy powders produced by plasma¿enhanced chemical vapor deposition (PECVD) are promising new materials for sintering ceramics, for making nanoscale filters, or for supporting catalytic surfaces. Common characteristics of these powders are their high purity and the easy control of their stoichiometry through the composition of the precursor gas mixture. Plasma parameters also influence their structure. Nanometric powders of silicon¿carbon alloys exhibiting microstructural properties such as large hydrogen content and high surface/volume ratio have been produced in a PECVD reactor using mixtures of silane and methane at low pressure (-1 Torr) and low frequency square¿wave modulated rf power (13.56 MHz). The a¿Si1¿xCx:H powders were obtained from different precursor gas mixtures, from R=0.05 to R=9, where R=[SiH4]/([SiH4]+[CH4]). The structure of the a¿Si1¿xCx:H powder was analyzed by several techniques. The particles appeared agglomerated, with a wide size distribution between 5 and 100 nm. The silane/methane gas mixture determined the vibrational features of these powders in the infrared. Silicon-hydrogen groups were present for every gas composition, whereas carbon¿hydrogen and silicon¿carbon bonds appeared in methane¿rich mixtures (R-0.6). The thermal desorption of hydrogen revealed two main evolutions at about 375 and 660¿°C that were ascribed to hydrogen bonded to silicon and carbon, respectively. The estimated hydrogen atom concentration in the sample was about 50%.

Weakly linked binary mixtures of F = 1 87Rb Bose-Einstein condensates

We present a study of binary mixtures of Bose-Einstein condensates confined in a double-well potential within the framework of the mean field Gross-Pitaevskii (GP) equation. We re-examine both the single component and the binary mixture cases for such a potential, and we investigate what are the situations in which a simpler two-mode approach leads to an accurate description of their dynamics. We also estimate the validity of the most usual dimensionality reductions used to solve the GP equations. To this end, we compare both the semi-analytical two-mode approaches and the numerical simulations of the one-dimensional (1D) reductions with the full 3D numerical solutions of the GP equation. Our analysis provides a guide to clarify the validity of several simplified models that describe mean-field nonlinear dynamics, using an experimentally feasible binary mixture of an F = 1 spinor condensate with two of its Zeeman manifolds populated, m = ±1.

Weakly linked binary mixtures of F = 1 87Rb Bose-Einstein condensates

We present a study of binary mixtures of Bose-Einstein condensates confined in a double-well potential within the framework of the mean field Gross-Pitaevskii (GP) equation. We re-examine both the single component and the binary mixture cases for such a potential, and we investigate what are the situations in which a simpler two-mode approach leads to an accurate description of their dynamics. We also estimate the validity of the most usual dimensionality reductions used to solve the GP equations. To this end, we compare both the semi-analytical two-mode approaches and the numerical simulations of the one-dimensional (1D) reductions with the full 3D numerical solutions of the GP equation. Our analysis provides a guide to clarify the validity of several simplified models that describe mean-field nonlinear dynamics, using an experimentally feasible binary mixture of an F = 1 spinor condensate with two of its Zeeman manifolds populated, m = ±1.

Simulated-annealing-based conditional simulation for the local-scale characterization of heterogeneous aquifers

Simulated-annealing-based conditional simulations provide a flexible means of quantitatively integrating diverse types of subsurface data. Although such techniques are being increasingly used in hydrocarbon reservoir characterization studies, their potential in environmental, engineering and hydrological investigations is still largely unexploited. Here, we introduce a novel simulated annealing (SA) algorithm geared towards the integration of high-resolution geophysical and hydrological data which, compared to more conventional approaches, provides significant advancements in the way that large-scale structural information in the geophysical data is accounted for. Model perturbations in the annealing procedure are made by drawing from a probability distribution for the target parameter conditioned to the geophysical data. This is the only place where geophysical information is utilized in our algorithm, which is in marked contrast to other approaches where model perturbations are made through the swapping of values in the simulation grid and agreement with soft data is enforced through a correlation coefficient constraint. Another major feature of our algorithm is the way in which available geostatistical information is utilized. Instead of constraining realizations to match a parametric target covariance model over a wide range of spatial lags, we constrain the realizations only at smaller lags where the available geophysical data cannot provide enough information. Thus we allow the larger-scale subsurface features resolved by the geophysical data to have much more due control on the output realizations. Further, since the only component of the SA objective function required in our approach is a covariance constraint at small lags, our method has improved convergence and computational efficiency over more traditional methods. Here, we present the results of applying our algorithm to the integration of porosity log and tomographic crosshole georadar data to generate stochastic realizations of the local-scale porosity structure. Our procedure is first tested on a synthetic data set, and then applied to data collected at the Boise Hydrogeophysical Research Site.

Pore-scale modeling of two-phase flow instabilities in porous media

Les problèmes d'écoulements multiphasiques en média poreux sont d'un grand intérêt pour de nombreuses applications scientifiques et techniques ; comme la séquestration de C02, l'extraction de pétrole et la dépollution des aquifères. La complexité intrinsèque des systèmes multiphasiques et l'hétérogénéité des formations géologiques sur des échelles multiples représentent un challenge majeur pour comprendre et modéliser les déplacements immiscibles dans les milieux poreux. Les descriptions à l'échelle supérieure basées sur la généralisation de l'équation de Darcy sont largement utilisées, mais ces méthodes sont sujettes à limitations pour les écoulements présentant de l'hystérèse. Les avancées récentes en terme de performances computationnelles et le développement de méthodes précises pour caractériser l'espace interstitiel ainsi que la distribution des phases ont favorisé l'utilisation de modèles qui permettent une résolution fine à l'échelle du pore. Ces modèles offrent un aperçu des caractéristiques de l'écoulement qui ne peuvent pas être facilement observées en laboratoire et peuvent être utilisé pour expliquer la différence entre les processus physiques et les modèles à l'échelle macroscopique existants. L'objet premier de la thèse se porte sur la simulation numérique directe : les équations de Navier-Stokes sont résolues dans l'espace interstitiel et la méthode du volume de fluide (VOF) est employée pour suivre l'évolution de l'interface. Dans VOF, la distribution des phases est décrite par une fonction fluide pour l'ensemble du domaine et des conditions aux bords particulières permettent la prise en compte des propriétés de mouillage du milieu poreux. Dans la première partie de la thèse, nous simulons le drainage dans une cellule Hele-Shaw 2D avec des obstacles cylindriques. Nous montrons que l'approche proposée est applicable même pour des ratios de densité et de viscosité très importants et permet de modéliser la transition entre déplacement stable et digitation visqueuse. Nous intéressons ensuite à l'interprétation de la pression capillaire à l'échelle macroscopique. Nous montrons que les techniques basées sur la moyenne spatiale de la pression présentent plusieurs limitations et sont imprécises en présence d'effets visqueux et de piégeage. Au contraire, une définition basée sur l'énergie permet de séparer les contributions capillaires des effets visqueux. La seconde partie de la thèse est consacrée à l'investigation des effets d'inertie associés aux reconfigurations irréversibles du ménisque causé par l'interface des instabilités. Comme prototype pour ces phénomènes, nous étudions d'abord la dynamique d'un ménisque dans un pore angulaire. Nous montrons que, dans un réseau de pores cubiques, les sauts et reconfigurations sont si fréquents que les effets d'inertie mènent à différentes configurations des fluides. A cause de la non-linéarité du problème, la distribution des fluides influence le travail des forces de pression, qui, à son tour, provoque une chute de pression dans la loi de Darcy. Cela suggère que ces phénomènes devraient être pris en compte lorsque que l'on décrit l'écoulement multiphasique en média poreux à l'échelle macroscopique. La dernière partie de la thèse s'attache à démontrer la validité de notre approche par une comparaison avec des expériences en laboratoire : un drainage instable dans un milieu poreux quasi 2D (une cellule Hele-Shaw avec des obstacles cylindriques). Plusieurs simulations sont tournées sous différentes conditions aux bords et en utilisant différents modèles (modèle intégré 2D et modèle 3D) afin de comparer certaines quantités macroscopiques avec les observations au laboratoire correspondantes. Malgré le challenge de modéliser des déplacements instables, où, par définition, de petites perturbations peuvent grandir sans fin, notre approche numérique apporte de résultats satisfaisants pour tous les cas étudiés. - Problems involving multiphase flow in porous media are of great interest in many scientific and engineering applications including Carbon Capture and Storage, oil recovery and groundwater remediation. The intrinsic complexity of multiphase systems and the multi scale heterogeneity of geological formations represent the major challenges to understand and model immiscible displacement in porous media. Upscaled descriptions based on generalization of Darcy's law are widely used, but they are subject to several limitations for flow that exhibit hysteric and history- dependent behaviors. Recent advances in high performance computing and the development of accurate methods to characterize pore space and phase distribution have fostered the use of models that allow sub-pore resolution. These models provide an insight on flow characteristics that cannot be easily achieved by laboratory experiments and can be used to explain the gap between physical processes and existing macro-scale models. We focus on direct numerical simulations: we solve the Navier-Stokes equations for mass and momentum conservation in the pore space and employ the Volume Of Fluid (VOF) method to track the evolution of the interface. In the VOF the distribution of the phases is described by a fluid function (whole-domain formulation) and special boundary conditions account for the wetting properties of the porous medium. In the first part of this thesis we simulate drainage in a 2-D Hele-Shaw cell filled with cylindrical obstacles. We show that the proposed approach can handle very large density and viscosity ratios and it is able to model the transition from stable displacement to viscous fingering. We then focus on the interpretation of the macroscopic capillary pressure showing that pressure average techniques are subject to several limitations and they are not accurate in presence of viscous effects and trapping. On the contrary an energy-based definition allows separating viscous and capillary contributions. In the second part of the thesis we investigate inertia effects associated with abrupt and irreversible reconfigurations of the menisci caused by interface instabilities. As a prototype of these phenomena we first consider the dynamics of a meniscus in an angular pore. We show that in a network of cubic pores, jumps and reconfigurations are so frequent that inertia effects lead to different fluid configurations. Due to the non-linearity of the problem, the distribution of the fluids influences the work done by pressure forces, which is in turn related to the pressure drop in Darcy's law. This suggests that these phenomena should be taken into account when upscaling multiphase flow in porous media. The last part of the thesis is devoted to proving the accuracy of the numerical approach by validation with experiments of unstable primary drainage in a quasi-2D porous medium (i.e., Hele-Shaw cell filled with cylindrical obstacles). We perform simulations under different boundary conditions and using different models (2-D integrated and full 3-D) and we compare several macroscopic quantities with the corresponding experiment. Despite the intrinsic challenges of modeling unstable displacement, where by definition small perturbations can grow without bounds, the numerical method gives satisfactory results for all the cases studied.

Short-term storage in vitro and large-scale propagation of grapevine genotypes

The objective of this work was to evaluate the large-scale propagation of grapevine genotypes after short-term storage in vitro. Microshoots from ten grapevine genotypes were used. The following storage temperatures were evaluated: 10, 20, and 25°C. After short-term storage, the shoots were propagated in up to five successive subcultures, to assess the large-scale propagation of the germplasm maintained under conditions of minimal growth. The propagated shoots were rooted in different concentrations of indolbutiric acid (IBA) and acclimatized in greenhouse. The best temperature for short-term storage in vitro and survival of the genotypes was 20°C. In the propagation phase, the highest number of shoots per explant was found in the subcultures 4 and 5, with averages of 4.9 and 4.8 shoots per explant, respectively. In the rooting phase, the best results for number of roots were obtained using a culture medium supplemented with 0.4 µmol L-1 of IBA, with an average of three roots per shoot. During the acclimation phase, a survival rate higher than 95% was achieved after 30 days in the greenhouse. Grapevine genotypes maintained for six months in vitro, at 20ºC, can be micropropagated in large scale.

Local and Regional Scale Determinants of Biodiversity Patterns in Boreal Agricultural Landscapes

This thesis examines the local and regional scale determinants of biodiversity patterns using existing species and environmental data. The research focuses on agricultural environments that have experienced rapid declines of biodiversity during past decades. Existing digital databases provide vast opportunities for habitat mapping, predictive mapping of species occurrences and richness and understanding the speciesenvironment relationships. The applicability of these databases depends on the required accuracy and quality of the data needed to answer the landscape ecological and biogeographical questions in hand. Patterns of biodiversity arise from confounded effects of different factors, such as climate, land cover and geographical location. Complementary statistical approaches that can show the relative effects of different factors are needed in biodiversity analyses in addition to classical multivariate models. Better understanding of the key factors underlying the variation in diversity requires the analyses of multiple taxonomic groups from different perspectives, such as richness, occurrence, threat status and population trends. The geographical coincidence of species richness of different taxonomic groups can be rather limited. This implies that multiple geographical regions should be taken into account in order to preserve various groups of species. Boreal agricultural biodiversity and in particular, distribution and richness of threatened species is strongly associated with various grasslands. Further, heterogeneous agricultural landscapes characterized by moderate field size, forest patches and non-crop agricultural habitats enhance the biodiversity of rural environments. From the landscape ecological perspective, the major threats to Finnish agricultural biodiversity are the decline of connected grassland habitat networks, and general homogenization of landscape structure resulting from both intensification and marginalization of agriculture. The maintenance of key habitats, such as meadows and pastures is an essential task in conservation of agricultural biodiversity. Furthermore, a larger landscape context should be incorporated in conservation planning and decision making processes in order to respond to the needs of different species and to maintain heterogeneous rural landscapes and viable agricultural diversity in the future.

Prediction of plant species distributions across six millennia.

The usefulness of species distribution models (SDMs) in predicting impacts of climate change on biodiversity is difficult to assess because changes in species ranges may take decades or centuries to occur. One alternative way to evaluate the predictive ability of SDMs across time is to compare their predictions with data on past species distributions. We use data on plant distributions, fossil pollen and current and mid-Holocene climate to test the ability of SDMs to predict past climate-change impacts. We find that species showing little change in the estimated position of their realized niche, with resulting good model performance, tend to be dominant competitors for light. Different mechanisms appear to be responsible for among-species differences in model performance. Confidence in predictions of the impacts of climate change could be improved by selecting species with characteristics that suggest little change is expected in the relationships between species occurrence and climate patterns.

Cross-scale analysis of the region effect on vascular plant species diversity in southern and northern European mountain ranges.

Background: The divergent glacial histories of southern and northern Europe affect present-day species diversity at coarse-grained scales in these two regions, but do these effects also penetrate to the more fine-grained scales of local communities?Methodology/Principal Findings: We carried out a cross-scale analysis to address this question for vascular plants in two mountain regions, the Alps in southern Europe and the Scandes in northern Europe, using environmentally paired vegetation plots in the two regions (n = 403 in each region) to quantify four diversity components: (i) total number of species occurring in a region (total gamma-diversity), (ii) number of species that could occur in a target plot after environmental filtering (habitat-specific gamma-diversity), (iii) pair-wise species compositional turnover between plots (plot-to-plot beta-diversity) and (iv) number of species present per plot (plot gamma-diversity). We found strong region effects on total gamma-diversity, habitat-specific gamma-diversity and plot-to-plot beta-diversity, with a greater diversity in the Alps even towards distances smaller than 50 m between plots. In contrast, there was a slightly greater plot alpha-diversity in the Scandes, but with a tendency towards contrasting region effects on high and low soil-acidity plots.Conclusions/Significance: We conclude that there are strong regional differences between coarse-grained (landscape- to regional-scale) diversity components of the flora in the Alps and the Scandes mountain ranges,but that these differences do not necessarily penetrate to the finest-grained (plot-scale) diversity component, at least not on acidic soils. Because different processes can lead to a similar pattern, we discuss the consistency of our results with Quaternary history and other divergent features between the two regions such as habitat connectivity, selection for vagility and environmental differences not accounted for in our analyses

Multi-scale effects of nestling diet on breeding performance in a terrestrial top predator inferred from stable isotope analysis

Inter-individual diet variation within populations is likely to have important ecological and evolutionary implications. The diet-fitness relationships at the individual level and the emerging population processes are, however, poorly understood for most avian predators inhabiting complex terrestrial ecosystems. In this study, we use an isotopic approach to assess the trophic ecology of nestlings in a long-lived raptor, the Bonelli"s eagle Aquila fasciata, and investigate whether nestling dietary breath and main prey consumption can affect the species" reproductive performance at two spatial scales: territories within populations and populations over a large geographic area. At the territory level, those breeding pairs whose nestlings consumed similar diets to the overall population (i.e. moderate consumption of preferred prey, but complemented by alternative prey categories) or those disproportionally consuming preferred prey were more likely to fledge two chicks. An increase in the diet diversity, however, related negatively with productivity. The age and replacements of breeding pair members had also an influence on productivity, with more fledglings associated to adult pairs with few replacements, as expected in long-lived species. At the population level, mean productivity was higher in those population-years with lower dietary breadth and higher diet similarity among territories, which was related to an overall higher consumption of preferred prey. Thus, we revealed a correspondence in diet-fitness relationships at two spatial scales: territories and populations. We suggest that stable isotope analyses may be a powerful tool to monitor the diet of terrestrial avian predators on large spatio-temporal scales, which could serve to detect potential changes in the availability of those prey on which predators depend for breeding. We encourage ecologists and evolutionary and conservation biologists concerned with the multi-scale fitness consequences of inter-individual variation in resource use to employ similar stable isotope-based approaches, which can be successfully applied to complex ecosystems such as the Mediterranean.

Sidelobe elimination for generalized synthetic discriminant functions by a two-filter correlation and subsequent postprocessing of the intensity distributions

One of the most important problems in optical pattern recognition by correlation is the appearance of sidelobes in the correlation plane, which causes false alarms. We present a method that eliminate sidelobes of up to a given height if certain conditions are satisfied. The method can be applied to any generalized synthetic discriminant function filter and is capable of rejecting lateral peaks that are even higher than the central correlation. Satisfactory results were obtained in both computer simulations and optical implementation.

On Option Price Information Content and Extraction of Implied Probability Distributions

In this study we used market settlement prices of European call options on stock index futures to extract implied probability distribution function (PDF). The method used produces a PDF of returns of an underlying asset at expiration date from implied volatility smile. With this method, the assumption of lognormal distribution (Black-Scholes model) is tested. The market view of the asset price dynamics can then be used for various purposes (hedging, speculation). We used the so called smoothing approach for implied PDF extraction presented by Shimko (1993). In our analysis we obtained implied volatility smiles from index futures markets (S&P 500 and DAX indices) and standardized them. The method introduced by Breeden and Litzenberger (1978) was then used on PDF extraction. The results show significant deviations from the assumption of lognormal returns for S&P500 options while DAX options mostly fit the lognormal distribution. A deviant subjective view of PDF can be used to form a strategy as discussed in the last section.

Imaging the time-integrated cerebral metabolic activity with subcellular resolution through nanometer-scale detection of biosynthetic products deriving from (13)C-glucose.

Glucose is the primary source of energy for the brain but also an important source of building blocks for proteins, lipids, and nucleic acids. Little is known about the use of glucose for biosynthesis in tissues at the cellular level. We demonstrate that local cerebral metabolic activity can be mapped in mouse brain tissue by quantitatively imaging the biosynthetic products deriving from [U-(13)C]glucose metabolism using a combination of in situ electron microscopy and secondary ion mass-spectroscopy (NanoSIMS). Images of the (13)C-label incorporated into cerebral ultrastructure with ca. 100nm resolution allowed us to determine the timescale on which the metabolic products of glucose are incorporated into different cells, their sub-compartments and organelles. These were mapped in astrocytes and neurons in the different layers of the motor cortex. We see evidence for high metabolic activity in neurons via the nucleus (13)C enrichment. We observe that in all the major cell compartments, such as e.g. nucleus and Golgi apparatus, neurons incorporate substantially higher concentrations of (13)C-label than astrocytes.

Environmental behavior of soils and mixtures of soil-whitewash mud

The present study sought to observe the behavior of soils in natural state and in mixtures, in different ratios, with the industrial solid residue called whitewash mud. The work was conducted with samples of typical soils from the region of Alagoinhas, Bahia-Brazil. Wet chemical analysis and atomic absorption spectrophotometry were used in order to obtain the classification of the industrial solid residue. Solubilization and leaching tests were performed and X-ray diffraction and electron microscopy techniques were carried out. The results showed that the whitewash mud was classified as non-inert, but with great capacity of heavy metal retention largely owed to the kaolinite and goethite presence in the clay fraction of the soils, making it difficult to have heavy metals readily available for exchange.