Inversion prepared coronary MR angiography: direct visualization of coronary blood flow.

PURPOSE: Visualization of coronary blood flow by means of a slice-selective inversion pre-pulse in concert with bright-blood coronary MRA. MATERIALS AND METHODS: Coronary magnetic resonance angiography (MRA) of the right coronary artery (RCA) was performed in eight healthy adult subjects on a 1.5 Tesla MR system (Gyroscan ACS-NT, Philips Medical Systems, Best, NL) using a free-breathing navigator-gated and cardiac-triggered 3D steady-state free-precession (SSFP) sequence with radial k-space sampling. Imaging was performed with and without a slice-selective inversion pre-pulse, which was positioned along the main axis of the coronary artery but perpendicular to the imaging volume. Objective image quality parameters such as SNR, CNR, maximal visible vessel length, and vessel border definition were analyzed. RESULTS: In contrast to conventional bright-blood 3D coronary MRA, the selective inversion pre-pulse provided a direct measure of coronary blood flow. In addition, CNR between the RCA and right ventricular blood pool was increased and the vessels had a tendency towards better delineation. Blood SNR and CNR between right coronary blood and epicardial fat were comparable in both sequences. CONCLUSION: The combination of a free-breathing navigator-gated and cardiac-triggered 3D SSFP sequence with a slice-selective inversion pre-pulse allows for direct and directional visualization of coronary blood flow with the additional benefit of improved contrast between coronary and right ventricular blood pool.

Flow cytometry in the study of cell death

In this report we present a concise review concerning the use of flow cytometric methods to characterize and differentiate between two different mechanisms of cell death, apoptosis and necrosis. The applications of these techniques to clinical and basic research are also considered. The following cell features are useful to characterize the mode of cell death: (1) activation of an endonuclease in apoptotic cells results in extraction of the low molecular weight DNA following cell permeabilization, which, in turn, leads to their decreased stainability with DNA-specific fluorochromes. Measurements of DNA content make it possible to identify apoptotic cells and to recognize the cell cycle phase specificity of apoptotic process; (2) plasma membrane integrity, which is lost in necrotic but not in apoptotic cells; (3) the decrease in forward light scatter, paralleled either by no change or an increase in side scatter, represent early changes during apoptosis. The data presented indicate that flow cytometry can be applied to basic research of the molecular and biochemical mechanisms of apoptosis, as well as in the clinical situations, where the ability to monitor early signs of apoptosis in some systems may be predictive for the outcome of some treatment protocols.

Genetic clusters and sex-biased gene flow in a unicolonial Formica ant.

BACKGROUND: Animal societies are diverse, ranging from small family-based groups to extraordinarily large social networks in which many unrelated individuals interact. At the extreme of this continuum, some ant species form unicolonial populations in which workers and queens can move among multiple interconnected nests without eliciting aggression. Although unicoloniality has been mostly studied in invasive ants, it also occurs in some native non-invasive species. Unicoloniality is commonly associated with very high queen number, which may result in levels of relatedness among nestmates being so low as to raise the question of the maintenance of altruism by kin selection in such systems. However, the actual relatedness among cooperating individuals critically depends on effective dispersal and the ensuing pattern of genetic structuring. In order to better understand the evolution of unicoloniality in native non-invasive ants, we investigated the fine-scale population genetic structure and gene flow in three unicolonial populations of the wood ant F. paralugubris. RESULTS: The analysis of geo-referenced microsatellite genotypes and mitochondrial haplotypes revealed the presence of cryptic clusters of genetically-differentiated nests in the three populations of F. paralugubris. Because of this spatial genetic heterogeneity, members of the same clusters were moderately but significantly related. The comparison of nuclear (microsatellite) and mitochondrial differentiation indicated that effective gene flow was male-biased in all populations. CONCLUSION: The three unicolonial populations exhibited male-biased and mostly local gene flow. The high number of queens per nest, exchanges among neighbouring nests and restricted long-distance gene flow resulted in large clusters of genetically similar nests. The positive relatedness among clustermates suggests that kin selection may still contribute to the maintenance of altruism in unicolonial populations if competition occurs among clusters.

N-dimensional dynamical systems exploiting instabilities in full

We report experimental and numerical results showing how certain N-dimensional dynamical systems are able to exhibit complex time evolutions based on the nonlinear combination of N-1 oscillation modes. The experiments have been done with a family of thermo-optical systems of effective dynamical dimension varying from 1 to 6. The corresponding mathematical model is an N-dimensional vector field based on a scalar-valued nonlinear function of a single variable that is a linear combination of all the dynamic variables. We show how the complex evolutions appear associated with the occurrence of successive Hopf bifurcations in a saddle-node pair of fixed points up to exhaust their instability capabilities in N dimensions. For this reason the observed phenomenon is denoted as the full instability behavior of the dynamical system. The process through which the attractor responsible for the observed time evolution is formed may be rather complex and difficult to characterize. Nevertheless, the well-organized structure of the time signals suggests some generic mechanism of nonlinear mode mixing that we associate with the cluster of invariant sets emerging from the pair of fixed points and with the influence of the neighboring saddle sets on the flow nearby the attractor. The generation of invariant tori is likely during the full instability development and the global process may be considered as a generalized Landau scenario for the emergence of irregular and complex behavior through the nonlinear superposition of oscillatory motions

Red blood cell microparticles and blood group antigens: an analysis by flow cytometry.

BACKGROUND: The storage of blood induces the formation of erythrocytes-derived microparticles. Their pathogenic role in blood transfusion is not known so far, especially the risk to trigger alloantibody production in the recipient. This work aims to study the expression of clinically significant blood group antigens on the surface of red blood cells microparticles. MATERIAL AND METHODS: Red blood cells contained in erythrocyte concentrates were stained with specific antibodies directed against blood group antigens and routinely used in immunohematology practice. After inducing erythrocytes vesiculation with calcium ionophore, the presence of blood group antigens was analysed by flow cytometry. RESULTS: The expression of several blood group antigens from the RH, KEL, JK, FY, MNS, LE and LU systems was detected on erythrocyte microparticles. The presence of M (MNS1), N (MNS2) and s (MNS4) antigens could not be demonstrated by flow cytometry, despite that glycophorin A and B were identified on microparticles using anti-CD235a and anti-MNS3. DISCUSSION: We conclude that blood group antigens are localized on erythrocytes-derived microparticles and probably keep their immunogenicity because of their capacity to bind specific antibody. Selective segregation process during vesiculation or their ability to elicit an immune response in vivo has to be tested by further studies.

S-system parameter estimation for noisy metabolic profiles using newton-flow analysis.

Biochemical systems are commonly modelled by systems of ordinary differential equations (ODEs). A particular class of such models called S-systems have recently gained popularity in biochemical system modelling. The parameters of an S-system are usually estimated from time-course profiles. However, finding these estimates is a difficult computational problem. Moreover, although several methods have been recently proposed to solve this problem for ideal profiles, relatively little progress has been reported for noisy profiles. We describe a special feature of a Newton-flow optimisation problem associated with S-system parameter estimation. This enables us to significantly reduce the search space, and also lends itself to parameter estimation for noisy data. We illustrate the applicability of our method by applying it to noisy time-course data synthetically produced from previously published 4- and 30-dimensional S-systems. In addition, we propose an extension of our method that allows the detection of network topologies for small S-systems. We introduce a new method for estimating S-system parameters from time-course profiles. We show that the performance of this method compares favorably with competing methods for ideal profiles, and that it also allows the determination of parameters for noisy profiles.

Computational Representation of Collaborative Learning Flow Patterns Using IMS Learning Design

The identification and integration of reusable and customizable CSCL (Computer Supported Collaborative Learning) may benefit from the capture of best practices in collaborative learning structuring. The authors have proposed CLFPs (Collaborative Learning Flow Patterns) as a way of collecting these best practices. To facilitate the process of CLFPs by software systems, the paper proposes to specify these patterns using IMS Learning Design (IMS-LD). Thus, teachers without technical knowledge can particularize and integrate CSCL tools. Nevertheless, the support of IMS-LD for describing collaborative learning activities has some deficiencies: the collaborative tools that can be defined in these activities are limited. Thus, this paper proposes and discusses an extension to IMS-LD that enables to specify several characteristics of the use of tools that mediate collaboration. In order to obtain a Unit of Learning based on a CLFP, a three stage process is also proposed. A CLFP-based Unit of Learning example is used to illustrate the process and the need of the proposed extension.

Towards reusability and tailorability in collaborative learning systems using IMS-LD and grid services

CSCL applications are complex distributed systems that posespecial requirements towards achieving success in educationalsettings. Flexible and efficient design of collaborative activitiesby educators is a key precondition in order to provide CSCL tailorable systems, capable of adapting to the needs of eachparticular learning environment. Furthermore, some parts ofthose CSCL systems should be reused as often as possible inorder to reduce development costs. In addition, it may be necessary to employ special hardware devices, computational resources that reside in other organizations, or even exceed thepossibilities of one specific organization. Therefore, theproposal of this paper is twofold: collecting collaborativelearning designs (scripting) provided by educators, based onwell-known best practices (collaborative learning flow patterns) in a standard way (IMS-LD) in order to guide the tailoring of CSCL systems by selecting and integrating reusable CSCL software units; and, implementing those units in the form of grid services offered by third party providers. More specifically, this paper outlines a grid-based CSCL system having these features and illustrates its potential scope and applicability by means of a sample collaborative learning scenario.

Inbreeding and sex-biased gene flow in the ant Formica exsecta.

The objective of this study was to assess breeding and dispersal patterns of both males and females in a monogyne (a single queen per colony) population of ants. Monogyny is commonly associated with extensive nuptial flights, presumably leading to considerable gene flow over large areas. Opposite to these expectations we found evidence of both inbreeding and sex-biased gene flow in a monogyne population of Formica exsecta. We found a significant degree of population subdivision at a local scale (within islands) for queens (females heading established colonies) and workers, but not for colony fathers (the males mated to the colony queens). However, we found little evidence of population subdivision at a larger scale (among islands). More conclusive support for sex-biased gene flow comes from the analysis of isolation by distance on the largest island, and from assignment tests revealing differences in female and male philopatry. The genetic similarity between pairs of queens decreased significantly when geographical distance increased, demonstrating limited dispersal and isolation by distance in queens. By contrast, we found no such pattern for colony fathers. Furthermore, a significantly greater fraction of colony queens were assigned as having originated from the population of residence, as compared to colony fathers. Inbreeding coefficients were significantly positive for workers, but not for mother queens. The queen-male relatedness coefficient of 0.23 (regression relatedness) indicates that mating occurs between fairly close relatives. These results suggest that some monogyne species of ants have complex dispersal and mating systems that can result in genetic isolation by distance over small geographical scales. More generally, this study also highlights the importance of identifying the relevant scale in analyses of population structure and dispersal.

Immunohistochemical and flow cytometric analysis of long-term label-retaining cells in the adult heart.

Cardiac-resident stem/progenitor cells have been identified based on expression of stem cell-associated antigens. However, no single surface marker allows to identify a definite cardiac stem/progenitor cell entity. Hence, functional stem cell markers have been extensively searched for. In homeostatic systems, stem cells divide infrequently and therefore retain DNA labels such as 5-bromo-2'-deoxyuridine, which are diluted with division. We used this method to analyze long-term label-retaining cells in the mouse heart after 14 days of 5-bromo-2'-deoxyuridine administration. Labeled cells were detected using immunohistochemical and flow-cytometric methods after varying chasing periods up to 12 months. Using mathematical models, the observed label dilution could consistently be described in the context of a 2-population model, whereby a population of rapidly dividing cells accounted for an accelerated early decline, and a population of slowly dividing cells accounted for decelerated dilution on longer time scales. Label-retaining cells were preferentially localized in the atria and apical region and stained negative for markers of the major cell lineages present in the heart. Most cells with long-term label-retention expressed stem cell antigen-1 (Sca-1). Sca-1(+)CD31(-) cells formed cell aggregates in culture, out of which lineage-negative (Lin(-))Sca-1(+)CD31(-) cells emerged, which could be cultured for many passages. These cells formed cardiospheres and showed differentiation potential into mesenchymal cell lineages. When cultured in cardiomyogenic differentiation medium, they expressed cardiac-specific genes. In conclusion, recognition of slow-cycling cells provides functional evidence of stem/progenitor cells in the heart. Lin(-)Sca-1(+)CD31(-) cardiac-derived progenitors have a potential for differentiation into cardiomyogenic and mesenchymal cell lineages.

Nutrient supply by mass flow and diffusion to maize plants in response to soil aggregate size and water potential

Nutrients are basically transported to the roots by mass flow and diffusion. The aim of this study was to quantify the contribution of these two mechanisms to the acquisition of macronutrients (N, P, K, Ca, Mg, and S) and cationic micronutrients (Fe, Mn, Zn, and Cu) by maize plants as well as xylem exudate volume and composition in response to soil aggregate size and water availability. The experiment was conducted in a greenhouse with samples of an Oxisol, from under two management systems: a region of natural savanna-like vegetation (Cerradão, CER) and continuous maize under conventional management for over 30 years (CCM). The treatments were arranged in a factorial [2 x (1 + 2) x 2] design, with two management systems (CER and CCM), (1 + 2) soil sifted through a 4 mm sieve and two aggregate classes (< 0.5 mm and 0.5 - 4.0 mm) and two soil matric potentials (-40 and -10 kPa). These were evaluated in a randomized block design with four replications. The experiment was conducted for 70 days after sowing. The influence of soil aggregate size and water potential on the nutrient transport mechanisms was highest in soil samples with higher nutrient concentrations in solution, in the CER system; diffusion became more relevant when water availability was higher and in aggregates < 0.5 mm. The volume of xylem exudate collected from maize plants increased with the decrease in aggregate size and the increased availability of soil water in the CER system. The highest Ca and Mg concentrations in the xylem exudate of plants grown on samples from the CER system were related to the high concentrations of these nutrients in the soil solution of this management system.

Contribution of macroporosity to water flux of a soil under different tillage systems

In view of the importance of the macroporosity for the water transport properties of soils, its quantitative assessment is a challenging task. Measurements of hydraulic conductivity (K) at different soil water tensions and the quantification of water-conducting macropores (θM) of a soil under different tillage systems could help understand the effects on the soil porous system and related hydraulic properties. The purpose of this study was to assess the effects of Conventional Tillage (CT), Chisel Plow (CP) and No Tillage (NT) on θM and on K; and to quantify the contribution of macroporosity to total water flux in a loam soil. A tension disc infiltrometer was used at two soil water pressure heads (-5 cm, and 0) to infer θM and K, during fallow. Macroporosity was determined based on the flow contribution between 0 and -5 cm water potentials (K0, K5, respectively), according to the Hagen-Poiseuille equation. The K0 values were statistically higher for CT than for NT and CP. The K5 values did not differ statistically among treatments. The mean K values varied between 0.20 and 3.70 cm/h. For CT, θM was significantly greater than for CP and NT, following the same trend as K0. No differences in θM were detected between CP and NT. With CT, the formation of water-conducting macropores with persistence until post-harvest was possible, while under CP preparation, the water-conducting macropores were not persistent. These results support the idea that tillage affects the soil water movement mainly by the resulting water-conducting macropores. Future studies on tillage effects on water movement should focus on macroporosity.

Dynamical Systems approach to Saffman-Taylor fingering. A Dynamical Solvability Scenario

A dynamical systems approach to competition of Saffman-Taylor fingers in a Hele-Shaw channel is developed. This is based on global analysis of the phase space flow of the low-dimensional ordinary-differential-equation sets associated with the classes of exact solutions of the problem without surface tension. Some simple examples are studied in detail. A general proof of the existence of finite-time singularities for broad classes of solutions is given. Solutions leading to finite-time interface pinchoff are also identified. The existence of a continuum of multifinger fixed points and its dynamical implications are discussed. We conclude that exact zero-surface tension solutions taken in a global sense as families of trajectories in phase space are unphysical because the multifinger fixed points are nonhyperbolic, and an unfolding does not exist within the same class of solutions. Hyperbolicity (saddle-point structure) of the multifinger fixed points is argued to be essential to the physically correct qualitative description of finger competition. The restoring of hyperbolicity by surface tension is proposed as the key point to formulate a generic dynamical solvability scenario for interfacial pattern selection.

Hydrological processes obtained on the plot scale under four simulated rainfall tests during the cycle of different crop systems

The cropping system influences the interception of water by plants, water storage in depressions on the soil surface, water infiltration into the soil and runoff. The aim of this study was to quantify some hydrological processes under no tillage cropping systems at the edge of a slope, in 2009 and 2010, in a Humic Dystrudept soil, with the following treatments: corn, soybeans, and common beans alone; and intercropped corn and common bean. Treatments consisted of four simulated rainfall tests at different times, with a planned intensity of 64 mm h-1 and 90 min duration. The first test was applied 18 days after sowing, and the others at 39, 75 and 120 days after the first test. Different times of the simulated rainfall and stages of the crop cycle affected soil water content prior to the rain, and the time runoff began and its peak flow and, thus, the surface hydrological processes. The depth of the runoff and the depth of the water intercepted by the crop + soil infiltration + soil surface storage were affected by the crop systems and the rainfall applied at different times. The corn crop was the most effective treatment for controlling runoff, with a water loss ratio of 0.38, equivalent to 75 % of the water loss ratio exhibited by common bean (0.51), the least effective treatment in relation to the others. Total water loss by runoff decreased linearly with an increase in the time that runoff began, regardless of the treatment; however, soil water content on the gravimetric basis increased linearly from the beginning to the end of the rainfall.

Multiscale descriptions of density-driven flow instabilities in porous media

Les instabilités engendrées par des gradients de densité interviennent dans une variété d'écoulements. Un exemple est celui de la séquestration géologique du dioxyde de carbone en milieux poreux. Ce gaz est injecté à haute pression dans des aquifères salines et profondes. La différence de densité entre la saumure saturée en CO2 dissous et la saumure environnante induit des courants favorables qui le transportent vers les couches géologiques profondes. Les gradients de densité peuvent aussi être la cause du transport indésirable de matières toxiques, ce qui peut éventuellement conduire à la pollution des sols et des eaux. La gamme d'échelles intervenant dans ce type de phénomènes est très large. Elle s'étend de l'échelle poreuse où les phénomènes de croissance des instabilités s'opèrent, jusqu'à l'échelle des aquifères à laquelle interviennent les phénomènes à temps long. Une reproduction fiable de la physique par la simulation numérique demeure donc un défi en raison du caractère multi-échelles aussi bien au niveau spatial et temporel de ces phénomènes. Il requiert donc le développement d'algorithmes performants et l'utilisation d'outils de calculs modernes. En conjugaison avec les méthodes de résolution itératives, les méthodes multi-échelles permettent de résoudre les grands systèmes d'équations algébriques de manière efficace. Ces méthodes ont été introduites comme méthodes d'upscaling et de downscaling pour la simulation d'écoulements en milieux poreux afin de traiter de fortes hétérogénéités du champ de perméabilité. Le principe repose sur l'utilisation parallèle de deux maillages, le premier est choisi en fonction de la résolution du champ de perméabilité (grille fine), alors que le second (grille grossière) est utilisé pour approximer le problème fin à moindre coût. La qualité de la solution multi-échelles peut être améliorée de manière itérative pour empêcher des erreurs trop importantes si le champ de perméabilité est complexe. Les méthodes adaptatives qui restreignent les procédures de mise à jour aux régions à forts gradients permettent de limiter les coûts de calculs additionnels. Dans le cas d'instabilités induites par des gradients de densité, l'échelle des phénomènes varie au cours du temps. En conséquence, des méthodes multi-échelles adaptatives sont requises pour tenir compte de cette dynamique. L'objectif de cette thèse est de développer des algorithmes multi-échelles adaptatifs et efficaces pour la simulation des instabilités induites par des gradients de densité. Pour cela, nous nous basons sur la méthode des volumes finis multi-échelles (MsFV) qui offre l'avantage de résoudre les phénomènes de transport tout en conservant la masse de manière exacte. Dans la première partie, nous pouvons démontrer que les approximations de la méthode MsFV engendrent des phénomènes de digitation non-physiques dont la suppression requiert des opérations de correction itératives. Les coûts de calculs additionnels de ces opérations peuvent toutefois être compensés par des méthodes adaptatives. Nous proposons aussi l'utilisation de la méthode MsFV comme méthode de downscaling: la grille grossière étant utilisée dans les zones où l'écoulement est relativement homogène alors que la grille plus fine est utilisée pour résoudre les forts gradients. Dans la seconde partie, la méthode multi-échelle est étendue à un nombre arbitraire de niveaux. Nous prouvons que la méthode généralisée est performante pour la résolution de grands systèmes d'équations algébriques. Dans la dernière partie, nous focalisons notre étude sur les échelles qui déterminent l'évolution des instabilités engendrées par des gradients de densité. L'identification de la structure locale ainsi que globale de l'écoulement permet de procéder à un upscaling des instabilités à temps long alors que les structures à petite échelle sont conservées lors du déclenchement de l'instabilité. Les résultats présentés dans ce travail permettent d'étendre les connaissances des méthodes MsFV et offrent des formulations multi-échelles efficaces pour la simulation des instabilités engendrées par des gradients de densité. - Density-driven instabilities in porous media are of interest for a wide range of applications, for instance, for geological sequestration of CO2, during which CO2 is injected at high pressure into deep saline aquifers. Due to the density difference between the C02-saturated brine and the surrounding brine, a downward migration of CO2 into deeper regions, where the risk of leakage is reduced, takes place. Similarly, undesired spontaneous mobilization of potentially hazardous substances that might endanger groundwater quality can be triggered by density differences. Over the last years, these effects have been investigated with the help of numerical groundwater models. Major challenges in simulating density-driven instabilities arise from the different scales of interest involved, i.e., the scale at which instabilities are triggered and the aquifer scale over which long-term processes take place. An accurate numerical reproduction is possible, only if the finest scale is captured. For large aquifers, this leads to problems with a large number of unknowns. Advanced numerical methods are required to efficiently solve these problems with today's available computational resources. Beside efficient iterative solvers, multiscale methods are available to solve large numerical systems. Originally, multiscale methods have been developed as upscaling-downscaling techniques to resolve strong permeability contrasts. In this case, two static grids are used: one is chosen with respect to the resolution of the permeability field (fine grid); the other (coarse grid) is used to approximate the fine-scale problem at low computational costs. The quality of the multiscale solution can be iteratively improved to avoid large errors in case of complex permeability structures. Adaptive formulations, which restrict the iterative update to domains with large gradients, enable limiting the additional computational costs of the iterations. In case of density-driven instabilities, additional spatial scales appear which change with time. Flexible adaptive methods are required to account for these emerging dynamic scales. The objective of this work is to develop an adaptive multiscale formulation for the efficient and accurate simulation of density-driven instabilities. We consider the Multiscale Finite-Volume (MsFV) method, which is well suited for simulations including the solution of transport problems as it guarantees a conservative velocity field. In the first part of this thesis, we investigate the applicability of the standard MsFV method to density- driven flow problems. We demonstrate that approximations in MsFV may trigger unphysical fingers and iterative corrections are necessary. Adaptive formulations (e.g., limiting a refined solution to domains with large concentration gradients where fingers form) can be used to balance the extra costs. We also propose to use the MsFV method as downscaling technique: the coarse discretization is used in areas without significant change in the flow field whereas the problem is refined in the zones of interest. This enables accounting for the dynamic change in scales of density-driven instabilities. In the second part of the thesis the MsFV algorithm, which originally employs one coarse level, is extended to an arbitrary number of coarse levels. We prove that this keeps the MsFV method efficient for problems with a large number of unknowns. In the last part of this thesis, we focus on the scales that control the evolution of density fingers. The identification of local and global flow patterns allows a coarse description at late times while conserving fine-scale details during onset stage. Results presented in this work advance the understanding of the Multiscale Finite-Volume method and offer efficient dynamic multiscale formulations to simulate density-driven instabilities. - Les nappes phréatiques caractérisées par des structures poreuses et des fractures très perméables représentent un intérêt particulier pour les hydrogéologues et ingénieurs environnementaux. Dans ces milieux, une large variété d'écoulements peut être observée. Les plus communs sont le transport de contaminants par les eaux souterraines, le transport réactif ou l'écoulement simultané de plusieurs phases non miscibles, comme le pétrole et l'eau. L'échelle qui caractérise ces écoulements est définie par l'interaction de l'hétérogénéité géologique et des processus physiques. Un fluide au repos dans l'espace interstitiel d'un milieu poreux peut être déstabilisé par des gradients de densité. Ils peuvent être induits par des changements locaux de température ou par dissolution d'un composé chimique. Les instabilités engendrées par des gradients de densité revêtent un intérêt particulier puisque qu'elles peuvent éventuellement compromettre la qualité des eaux. Un exemple frappant est la salinisation de l'eau douce dans les nappes phréatiques par pénétration d'eau salée plus dense dans les régions profondes. Dans le cas des écoulements gouvernés par les gradients de densité, les échelles caractéristiques de l'écoulement s'étendent de l'échelle poreuse où les phénomènes de croissance des instabilités s'opèrent, jusqu'à l'échelle des aquifères sur laquelle interviennent les phénomènes à temps long. Etant donné que les investigations in-situ sont pratiquement impossibles, les modèles numériques sont utilisés pour prédire et évaluer les risques liés aux instabilités engendrées par les gradients de densité. Une description correcte de ces phénomènes repose sur la description de toutes les échelles de l'écoulement dont la gamme peut s'étendre sur huit à dix ordres de grandeur dans le cas de grands aquifères. Il en résulte des problèmes numériques de grande taille qui sont très couteux à résoudre. Des schémas numériques sophistiqués sont donc nécessaires pour effectuer des simulations précises d'instabilités hydro-dynamiques à grande échelle. Dans ce travail, nous présentons différentes méthodes numériques qui permettent de simuler efficacement et avec précision les instabilités dues aux gradients de densité. Ces nouvelles méthodes sont basées sur les volumes finis multi-échelles. L'idée est de projeter le problème original à une échelle plus grande où il est moins coûteux à résoudre puis de relever la solution grossière vers l'échelle de départ. Cette technique est particulièrement adaptée pour résoudre des problèmes où une large gamme d'échelle intervient et évolue de manière spatio-temporelle. Ceci permet de réduire les coûts de calculs en limitant la description détaillée du problème aux régions qui contiennent un front de concentration mobile. Les aboutissements sont illustrés par la simulation de phénomènes tels que l'intrusion d'eau salée ou la séquestration de dioxyde de carbone.