A fully adaptive numerical approximation for a two-dimensional epidemic model with nonlinear cross-diffusion

An epidemic model is formulated by a reactionâeuro"diffusion system where the spatial pattern formation is driven by cross-diffusion. The reaction terms describe the local dynamics of susceptible and infected species, whereas the diffusion terms account for the spatial distribution dynamics. For both self-diffusion and cross-diffusion, nonlinear constitutive assumptions are suggested. To simulate the pattern formation two finite volume formulations are proposed, which employ a conservative and a non-conservative discretization, respectively. An efficient simulation is obtained by a fully adaptive multiresolution strategy. Numerical examples illustrate the impact of the cross-diffusion on the pattern formation.

On some numerical aspects of an active strain model in cardiac mechanics

We are interested in the development, implementation and testing of an orthotropic model for cardiac contraction based on an active strain decomposition. Our model addresses the coupling of a transversely isotropic mechanical description at the cell level, with an orthotropic constitutive law for incompressible tissue at the macroscopic level. The main differences with the active stress model are addressed in detail, and a finite element discretization using Taylor-Hood and MINI elements is proposed and illustrated with numerical examples.

The mycorrhizal contribution to plant productivity, plant nutrition and soil structure in experimental grassland.

Previous studies have shown that arbuscular mycorrhizal fungi (AMF) can influence plant diversity and ecosystem productivity. However, little is known about the effects of AMF and different AMF taxa on other important community properties such as nutrient acquisition, plant survival and soil structure. We established experimental grassland microcosms and tested the impact of AMF and of different AMF taxa on a number of grassland characteristics. We also tested whether plant species benefited from the same or different AMF taxa in subsequent growing seasons. AMF enhanced phosphorus acquisition, soil aggregation and survival of several plant species, but AMF did not increase total plant productivity. Moreover, AMF increased nitrogen acquisition by some plant species, but AMF had no effect on total N uptake by the plant community. Plant growth responses to AMF were temporally variable and some plant species obtained the highest biomass with different AMF in different years. Hence the results indicate that it may be beneficial for a plant to be colonized by different AMF taxa in different seasons. This study shows that AMF play a key role in grassland by improving plant nutrition and soil structure, and by regulating the make-up of the plant community.

Numerical analysis of wave-induced fluid flow effects on seismic data: Application to monitoring of CO2 storage at the Sleipner Field

In this work we analyze how patchy distributions of CO2 and brine within sand reservoirs may lead to significant attenuation and velocity dispersion effects, which in turn may have a profound impact on surface seismic data. The ultimate goal of this paper is to contribute to the understanding of these processes within the framework of the seismic monitoring of CO2 sequestration, a key strategy to mitigate global warming. We first carry out a Monte Carlo analysis to study the statistical behavior of attenuation and velocity dispersion of compressional waves traveling through rocks with properties similar to those at the Utsira Sand, Sleipner field, containing quasi-fractal patchy distributions of CO2 and brine. These results show that the mean patch size and CO2 saturation play key roles in the observed wave-induced fluid flow effects. The latter can be remarkably important when CO2 concentrations are low and mean patch sizes are relatively large. To analyze these effects on the corresponding surface seismic data, we perform numerical simulations of wave propagation considering reservoir models and CO2 accumulation patterns similar to the CO2 injection site in the Sleipner field. These numerical experiments suggest that wave-induced fluid flow effects may produce changes in the reservoir's seismic response, modifying significantly the main seismic attributes usually employed in the characterization of these environments. Consequently, the determination of the nature of the fluid distributions as well as the proper modeling of the seismic data constitute important aspects that should not be ignored in the seismic monitoring of CO2 sequestration problems.

Persistence of a biocontrol Pseudomonas inoculant as high populations of culturable and non-culturable cells in 200-cm-deep soil profiles

Little is known about the ecology of soil inoculants used for pathogen biocontrol, biofertilization and bioremediation under field conditions. We investigated the persistence and the physiological states of soil-inoculated Pseudomonas protegens (previously Pseudomonas fluorescens) CHA0 (108 CFU g−1 surface soil) in different soil microbial habitats in a planted ley (Medicago sativa L.) and an uncovered field plot. At 72 days, colony counts of the inoculant were low in surface soil (uncovered plot) and earthworm guts (ley plot), whereas soil above the plow pan (uncovered plot), and the rhizosphere and worm burrows present until 1.2 m depth (ley plot) were survival hot spots (105-106 CFU g−1 soil). Interestingly, strain CHA0 was also detected in the subsoil of both plots, at 102-105 CFU g−1 soil between 1.8 and 2 m depth. However, non-cultured CHA0 cells were also evidenced based on immunofluorescence microscopy. Kogure's direct viable counts of nutrient-responsive cells showed that many more CHA0 cells were in a viable but non-culturable (VBNC) or a non-responsive (dormant) state than in a culturable state, and the proportion of cells in those non-cultured states depended on soil microbial habitat. At the most, cells in a VBNC state amounted to 34% (above the plow pan) and those in a dormant state to 89% (in bulk soil between 0.6 and 2 m) of all CHA0 cells. The results indicate that field-released Pseudomonas inoculants may persist at high cell numbers, even in deeper soil layers, and display a combination of different physiological states whose prevalence fluctuates according to soil microbial habitats.

Identification of opsA, a gene involved in solute stress mitigation and survival in soil, in the polycyclic aromatic hydrocarbon-degrading bacterium Novosphingobium sp. strain LH128.

The aim of this study was to identify genes involved in solute and matric stress mitigation in the polycyclic aromatic hydrocarbon (PAH)-degrading Novosphingobium sp. strain LH128. The genes were identified using plasposon mutagenesis and by selection of mutants that showed impaired growth in a medium containing 450 mM NaCl as a solute stress or 10% (wt/vol) polyethylene glycol (PEG) 6000 as a matric stress. Eleven and 14 mutants showed growth impairment when exposed to solute and matric stresses, respectively. The disrupted sequences were mapped on a draft genome sequence of strain LH128, and the corresponding gene functions were predicted. None of them were shared between solute and matric stress-impacted mutants. One NaCl-affected mutant (i.e., NA7E1) with a disruption in a gene encoding a putative outer membrane protein (OpsA) was susceptible to lower NaCl concentrations than the other mutants. The growth of NA7E1 was impacted by other ions and nonionic solutes and by sodium dodecyl sulfate (SDS), suggesting that opsA is involved in osmotic stress mitigation and/or outer membrane stability in strain LH128. NA7E1 was also the only mutant that showed reduced growth and less-efficient phenanthrene degradation in soil compared to the wild type. Moreover, the survival of NA7E1 in soil decreased significantly when the moisture content was decreased but was unaffected when soluble solutes from sandy soil were removed by washing. opsA appears to be important for the survival of strain LH128 in soil, especially in the case of reduced moisture content, probably by mitigating the effects of solute stress and retaining membrane stability.

Characterization of gravitational rock slope deformations at different spatial scales based on field, remote sensing and numerical approaches

Rock slope instabilities such as rock slides, rock avalanche or deep-seated gravitational slope deformations are widespread in Alpine valleys. These phenomena represent at the same time a main factor that control the mountain belts erosion and also a significant natural hazard that creates important losses to the mountain communities. However, the potential geometrical and dynamic connections linking outcrop and slope-scale instabilities are often unknown. A more detailed definition of the potential links will be essential to improve the comprehension of the destabilization processes and to dispose of a more complete hazard characterization of the rock instabilities at different spatial scales. In order to propose an integrated approach in the study of the rock slope instabilities, three main themes were analysed in this PhD thesis: (1) the inventory and the spatial distribution of rock slope deformations at regional scale and their influence on the landscape evolution, (2) the influence of brittle and ductile tectonic structures on rock slope instabilities development and (3) the characterization of hazard posed by potential rock slope instabilities through the development of conceptual instability models. To prose and integrated approach for the analyses of these topics, several techniques were adopted. In particular, high resolution digital elevation models revealed to be fundamental tools that were employed during the different stages of the rock slope instability assessment. A special attention was spent in the application of digital elevation model for detailed geometrical modelling of past and potential instabilities and for the rock slope monitoring at different spatial scales. Detailed field analyses and numerical models were performed to complete and verify the remote sensing approach. In the first part of this thesis, large slope instabilities in Rhone valley (Switzerland) were mapped in order to dispose of a first overview of tectonic and climatic factors influencing their distribution and their characteristics. Our analyses demonstrate the key influence of neotectonic activity and the glacial conditioning on the spatial distribution of the rock slope deformations. Besides, the volumes of rock instabilities identified along the main Rhone valley, were then used to propose the first estimate of the postglacial denudation and filling of the Rhone valley associated to large gravitational movements. In the second part of the thesis, detailed structural analyses of the Frank slide and the Sierre rock avalanche were performed to characterize the influence of brittle and ductile tectonic structures on the geometry and on the failure mechanism of large instabilities. Our observations indicated that the geometric characteristics and the variation of the rock mass quality associated to ductile tectonic structures, that are often ignored landslide study, represent important factors that can drastically influence the extension and the failure mechanism of rock slope instabilities. In the last part of the thesis, the failure mechanisms and the hazard associated to five potential instabilities were analysed in detail. These case studies clearly highlighted the importance to incorporate different analyses and monitoring techniques to dispose of reliable and hazard scenarios. This information associated to the development of a conceptual instability model represents the primary data for an integrated risk management of rock slope instabilities. - Les mouvements de versant tels que les chutes de blocs, les éboulements ou encore les phénomènes plus lents comme les déformations gravitaires profondes de versant représentent des manifestations courantes en régions montagneuses. Les mouvements de versant sont à la fois un des facteurs principaux contrôlant la destruction progressive des chaines orogéniques mais aussi un danger naturel concret qui peut provoquer des dommages importants. Pourtant, les phénomènes gravitaires sont rarement analysés dans leur globalité et les rapports géométriques et mécaniques qui lient les instabilités à l'échelle du versant aux instabilités locales restent encore mal définis. Une meilleure caractérisation de ces liens pourrait pourtant représenter un apport substantiel dans la compréhension des processus de déstabilisation des versants et améliorer la caractérisation des dangers gravitaires à toutes les échelles spatiales. Dans le but de proposer un approche plus globale à la problématique des mouvements gravitaires, ce travail de thèse propose trois axes de recherche principaux: (1) l'inventaire et l'analyse de la distribution spatiale des grandes instabilités rocheuses à l'échelle régionale, (2) l'analyse des structures tectoniques cassantes et ductiles en relation avec les mécanismes de rupture des grandes instabilités rocheuses et (3) la caractérisation des aléas rocheux par une approche multidisciplinaire visant à développer un modèle conceptuel de l'instabilité et une meilleure appréciation du danger . Pour analyser les différentes problématiques traitées dans cette thèse, différentes techniques ont été utilisées. En particulier, le modèle numérique de terrain s'est révélé être un outil indispensable pour la majorité des analyses effectuées, en partant de l'identification de l'instabilité jusqu'au suivi des mouvements. Les analyses de terrain et des modélisations numériques ont ensuite permis de compléter les informations issues du modèle numérique de terrain. Dans la première partie de cette thèse, les mouvements gravitaires rocheux dans la vallée du Rhône (Suisse) ont été cartographiés pour étudier leur répartition en fonction des variables géologiques et morphologiques régionales. En particulier, les analyses ont mis en évidence l'influence de l'activité néotectonique et des phases glaciaires sur la distribution des zones à forte densité d'instabilités rocheuses. Les volumes des instabilités rocheuses identifiées le long de la vallée principale ont été ensuite utilisés pour estimer le taux de dénudations postglaciaire et le remplissage de la vallée du Rhône lié aux grands mouvements gravitaires. Dans la deuxième partie, l'étude de l'agencement structural des avalanches rocheuses de Sierre (Suisse) et de Frank (Canada) a permis de mieux caractériser l'influence passive des structures tectoniques sur la géométrie des instabilités. En particulier, les structures issues d'une tectonique ductile, souvent ignorées dans l'étude des instabilités gravitaires, ont été identifiées comme des structures très importantes qui contrôlent les mécanismes de rupture des instabilités à différentes échelles. Dans la dernière partie de la thèse, cinq instabilités rocheuses différentes ont été étudiées par une approche multidisciplinaire visant à mieux caractériser l'aléa et à développer un modèle conceptuel trois dimensionnel de ces instabilités. A l'aide de ces analyses on a pu mettre en évidence la nécessité d'incorporer différentes techniques d'analyses et de surveillance pour une gestion plus objective du risque associée aux grandes instabilités rocheuses.

Functionally and phylogenetically diverse plant communities key to soil biota

Recent studies assessing the role of biological diversity for ecosystem functioning indicate that the diversity of functional traits and the evolutionary history of species in a community, not the number of taxonomic units, ultimately drives the biodiversity-ecosystem-function relationship. Here, we simultaneously assessed the importance of plant functional trait and phylogenetic diversity as predictors of major trophic groups of soil biota (abundance and diversity), six years from the onset of a grassland biodiversity experiment. Plant functional and phylogenetic diversity were generally better predictors of soil biota than the traditionally used species or functional group richness. Functional diversity was a reliable predictor for most biota, with the exception of soil microorganisms, which were better predicted by phylogenetic diversity. These results provide empirical support for the idea that the diversity of plant functional traits and the diversity of evolutionary lineages in a community are important for maintaining higher abundances and diversity of soil communities.

Numerical simulation of gel electrophoresis of DNA knots in weak and strong electric fields.

Gel electrophoresis allows one to separate knotted DNA (nicked circular) of equal length according to the knot type. At low electric fields, complex knots, being more compact, drift faster than simpler knots. Recent experiments have shown that the drift velocity dependence on the knot type is inverted when changing from low to high electric fields. We present a computer simulation on a lattice of a closed, knotted, charged DNA chain drifting in an external electric field in a topologically restricted medium. Using a Monte Carlo algorithm, the dependence of the electrophoretic migration of the DNA molecules on the knot type and on the electric field intensity is investigated. The results are in qualitative and quantitative agreement with electrophoretic experiments done under conditions of low and high electric fields.

Flow structures at an idealized bifurcation : a numerical experiment

River bifurcations are key nodes within braided river systems controlling the flow and sediment partitioning and therefore the dynamics of the river braiding process. Recent research has shown that certain geometrical configurations induce instabilities that lead to downstream mid-channel bar formation and the formation of bifurcations. However, we currently have a poor understanding of the flow division process within bifurcations and the flow dynamics in the downstream bifurcates, both of which are needed to understand bifurcation stability. This paper presents results of a numerical sensitivity experiment undertaken using computational fluid dynamics (CFD) with the purpose of understanding the flow dynamics of a series of idealized bifurcations. A geometric sensitivity analysis is undertaken for a range of channel slopes (0.005 to 0.03), bifurcation angles (22 degrees to 42 degrees) and a restricted set of inflow conditions based upon simulating flow through meander bends with different curvature on the flow field dynamics through the bifurcation. The results demonstrate that the overall slope of the bifurcation affects the velocity of flow through the bifurcation and when slope asymmetry is introduced, the flow structures in the bifurcation are modified. In terms of bifurcation evolution the most important observation appears to be that once slope asymmetry is greater than 0.2 the flow within the steep bifurcate shows potential instability and the potential for alternate channel bar formation. Bifurcation angle also defines the flow structures within the bifurcation with an increase in bifurcation angle increasing the flow velocity down both bifurcates. However, redistributive effects of secondary circulation caused by upstream curvature can very easily counter the effects of local bifurcation characteristics. Copyright (C) 2011 John Wiley & Sons, Ltd.

Density-based hierarchical clustering of pyro-sequences on a large scale--the case of fungal ITS1.

MOTIVATION: Analysis of millions of pyro-sequences is currently playing a crucial role in the advance of environmental microbiology. Taxonomy-independent, i.e. unsupervised, clustering of these sequences is essential for the definition of Operational Taxonomic Units. For this application, reproducibility and robustness should be the most sought after qualities, but have thus far largely been overlooked. RESULTS: More than 1 million hyper-variable internal transcribed spacer 1 (ITS1) sequences of fungal origin have been analyzed. The ITS1 sequences were first properly extracted from 454 reads using generalized profiles. Then, otupipe, cd-hit-454, ESPRIT-Tree and DBC454, a new algorithm presented here, were used to analyze the sequences. A numerical assay was developed to measure the reproducibility and robustness of these algorithms. DBC454 was the most robust, closely followed by ESPRIT-Tree. DBC454 features density-based hierarchical clustering, which complements the other methods by providing insights into the structure of the data. AVAILABILITY: An executable is freely available for non-commercial users at ftp://ftp.vital-it.ch/tools/dbc454. It is designed to run under MPI on a cluster of 64-bit Linux machines running Red Hat 4.x, or on a multi-core OSX system. CONTACT: dbc454@vital-it.ch or nicolas.guex@isb-sib.ch.

Prothèses totales de genoux. Analyse clinique et numérique des micromouvements de l'implant tibial [Total knee prosthesis. Clinical and numerical study of micromovements of the tibial implant].

INTRODUCTION: The importance of the micromovements in the mechanism of aseptic loosening is clinically difficult to evaluate. To complete the analysis of a series of total knee arthroplasties (TKA), we used a tridimensional numerical model to study the micromovements of the tibial implant. MATERIAL AND METHODS: Fifty one patients (with 57 cemented Porous Coated Anatomic TKAs) were reviewed (mean follow-up 4.5 year). Radiolucency at the tibial bone-cement interface was sought on the AP radiographs and divided in 7 areas. The distribution of the radiolucency was then correlated with the axis of the lower limb as measured on the orthoradiograms. The tridimensional numerical model is based on the finite element method. It allowed the measurement of the cemented prosthetic tibial implant's displacements and the micromovements generated at bone-ciment interface. A total load (2000 Newton) was applied at first vertically and asymetrically on the tibial plateau, thereby simulating an axial deviation of the lower limbs. The vector's posterior inclination then permitted the addition of a tangential component to the axial load. This type of effort is generated by complex biomechanical phenomena such as knee flexion. RESULTS: 81 per cent of the 57 knees had a radiolucent line of at least 1 mm, at one or more of the tibial cement-epiphysis jonctional areas. The distribution of these lucent lines showed that they came out more frequently at the periphery of the implant. The lucent lines appeared most often under the unloaded margin of the tibial plateau, when axial deviation of lower limbs was present. Numerical simulations showed that asymetrical loading on the tibial plateau induced a subsidence of the loaded margin (0-100 microns) and lifting off at the opposite border (0-70 microns). The postero-anterior tangential component induced an anterior displacement of the tibial implant (160-220 microns), and horizontal micromovements with non homogenous distribution at the bone-ciment interface (28-54 microns). DISCUSSION: Comparison of clinical and numerical results showed a relation between the development of radiolucent lines and the unloading of the tibial implant's margin. The deleterious effect of lower limbs' axial deviation is thereby proven. The irregular distribution of lucent lines under the tibial plateau was similar of the micromovements' repartition at the bone-cement interface when tangential forces were present. A causative relation between the two phenomenaes could not however be established. Numerical simulation is a truly useful method of study; it permits to calculate micromovements which are relative, non homogenous and of very low amplitude. However, comparative clinical studies remain as essential to ensure the credibility of results.

Heated Optical Fiber for Distributed Soil-Moisture Measurements: A Lysimeter Experiment

An Actively Heated Fiber Optics (AHFO) method to estimate soil moisture is tested and the analysis technique improved on. The measurements were performed in a lysimeter uniformly packed with loam soil with variable water content profiles. In the first meter of the soil profi le, 30 m of fiber optic cable were installed in a 12 loops coil. The metal sheath armoring the fiber cable was used as an electrical resistance heater to generate a heat pulse, and the soil response was monitored with a Distributed Temperature Sensing (DTS) system. We study the cooling following three continuous heat pulses of 120 s at 36 W m(-1) by means of long-time approximation of radial heat conduction. The soil volumetric water contents were then inferred from the estimated thermal conductivities through a specifically calibrated model relating thermal conductivity and volumetric water content. To use the pre-asymptotic data we employed a time correction that allowed the volumetric water content to be estimated with a precision of 0.01-0.035 (m(3) m(-3)). A comparison of the AHFO measurements with soil-moisture measurements obtained with calibrated capacitance-based probes gave good agreement for wetter soils [discrepancy between the two methods was less than 0.04 (m(3) m(-3))]. In the shallow drier soils, the AHFO method underestimated the volumetric water content due to the longertime required for the temperature increment to become asymptotic in less thermally conductive media [discrepancy between the two methods was larger than 0.1 (m(3) m(-3))]. The present work suggests that future applications of the AHFO method should include longer heat pulses, that longer heating and cooling events are analyzed, and, temperature increments ideally be measured with higher frequency.