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.

Modelling the influence of change in fire regime on the local distribution of a Mediterranean pyrophytic plant species (Cistus salviifolius) at its northern range limit

Aims: To assess the potential distribution of an obligate seeder and active pyrophyte, Cistus salviifolius, a vulnerable species in the Swiss Red List; to derive scenarios by changing the fire return interval; and to discuss the results from a conservation perspective. A more general aim is to assess the impact of fire as a natural factor influencing the vegetation of the southern slopes of the Alps. Locations: Alps, southern Switzerland. Methods: Presence-absence data to fit the model were obtained from the most recent field mapping of C. salviifolius. The quantitative environmental predictors used in this study include topographic, climatic and disturbance (fire) predictors. Models were fitted by logistic regression and evaluated by jackknife and bootstrap approaches. Changes in fire regime were simulated by increasing the time-return interval of fire (simulating longer periods without fire). Two scenarios were considered: no fire in the past 15 years; or in the past 35 years. Results: Rock cover, slope, topographic position, potential evapotranspiration and time elapsed since the last fire were selected in the final model. The Nagelkerke R-2 of the model for C. salviifolius was 0.57 and the Jackknife area under the curve evaluation was 0.89. The bootstrap evaluation revealed model robustness. By increasing the return interval of fire by either up to 15 years, or 35 years, the modelled C. salviifolius population declined by 30-40%, respectively. Main conclusions: Although fire plays a significant role, topography and rock cover appear to be the most important predictors, suggesting that the distribution of C. salviifolius in the southern Swiss Alps is closely related to the availability of supposedly competition-free sites, such as emerging bedrock, ridge locations or steep slopes. Fire is more likely to play a secondary role in allowing C. salviifolius to extend its occurrence temporarily, by increasing germination rates and reducing the competition from surrounding vegetation. To maintain a viable dormant seed bank for C. salviifolius, conservation managers should consider carrying out vegetation clearing and managing wild fire propagation to reduce competition and ensure sufficient recruitment for this species.

Growth decrease in Alpine whitefish: investigating the relative contribution of fishing induced selection and environmental change and its implication for conservation measures

La pression exercée par les activités humaines menace pratiquement tous les écosystèmes aquatiques du globe. Ainsi, sous l'effet de divers facteurs tels que la pollution, le réchauffement climatique ou encore la pêche industrielle, de nombreuses populations de poissons ont vu leurs effectifs chuter et divers changements morphologiques ont été observés. Dans cette étude, nous nous sommes intéressés à une menace particulière: la sélection induite par la pêche sur la croissance des poissons. En effet, la génétique des populations prédit que la soustraction régulière des individus les plus gros peut entraîner des modifications rapides de certains traits physiques comme la croissance individuelle. Cela a par ailleurs été observé dans de nombreuses populations marines ou lacustres, dont les populations de féras, bondelles et autres corégones des lacs suisses. Toutefois, malgré un nombre croissant d'études décrivant ce phénomène, peu de plans de gestion en tiennent compte, car l'importance des effets génétiques liés à la pêche est le plus souvent négligée par rapport à l'impact des changements environnementaux. Le but premier de cette étude a donc été de quantifier l'importance des facteurs génétiques et environnementaux. Dans le premier chapitre, nous avons étudié la population de palée du lac de Joux (Coregonus palaea). Nous avons déterminé les différentiels de sélection dus à la pêche, c'est-à-dire l'intensité de la sélection sur le taux de croissance, ainsi que les changements nets de croissance au cours du temps. Nous avons observé une baisse marquée de croissance et un différentiel de sélection important indiquant qu'au moins 30% de la diminution de croissance observée était due à la pression de sélection induite par la pêche. Dans le deuxième chapitre, nous avons effectué les mêmes analyses sur deux espèces proches du lac de Brienz (C. albellus et C. fatioi) et avons observé des effets similaires dont l'intensité était spécifique à chaque espèce. Dans le troisième chapitre, nous avons analysé deux autres espèces : C. palaea et C. confusus du lac de Bienne, et avons constaté que le lien entre la pression de sélection et la diminution de croissance était influencé par des facteurs environnementaux. Finalement, dans le dernier chapitre, nous avons étudié les effets potentiels de différentes modifications de la taille des mailles des filets utilisés pour la pêche à l'aide de modèles mathématiques. Nous concluons que la pêche a un effet génétique non négligeable (et donc peu réversible) sur la croissance individuelle dans les populations observée, que cet effet est lié à la compétition pour la nourriture et à la qualité de l'environnement, et que certaines modifications simples de la taille des mailles des filets de pêche pourraient nettement diminuer l'effet de sélection et ainsi ralentir, voir même renverser la diminution de croissance observée.

The effect of continuous positive airway pressure on total cerebral blood flow in healthy awake volunteers.

PURPOSE: Continuous positive airway pressure (CPAP) is the gold standard treatment for obstructive sleep apnea. However, the physiologic impact of CPAP on cerebral blood flow (CBF) is not well established. Ultrasound can be used to estimate CBF, but there is no widespread accepted protocol. We studied the physiologic influence of CPAP on CBF using a method integrating arterial diameter and flow velocity (FV) measurements obtained for each vessel supplying blood to the brain. METHODS: FV and lumen diameter of the left and right internal carotid, vertebral, and middle cerebral arteries were measured using duplex Doppler ultrasound with and without CPAP at 15 cm H(2)O, applied in a random order. Transcutaneous carbon dioxide (PtcCO(2)), heart rate (HR), blood pressure (BP), and oxygen saturation were monitored. Results were compared with a theoretical prediction of CBF change based on the effect of partial pressure of carbon dioxide on CBF. RESULTS: Data were obtained from 23 healthy volunteers (mean ± SD; 12 male, age 25.1 ± 2.6 years, body mass index 21.8 ± 2.0 kg/m(2)). The mean experimental and theoretical CBF decrease under CPAP was 12.5 % (p < 0.001) and 11.9 % (p < 0.001), respectively. The difference between experimental and theoretical CBF reduction was not statistically significant (3.84 ± 79 ml/min, p = 0.40). There was a significant reduction in PtcCO(2) with CPAP (p = <0.001) and a significant increase in mean BP (p = 0.0017). No significant change was observed in SaO(2) (p = 0.21) and HR (p = 0.62). CONCLUSION: Duplex Doppler ultrasound measurements of arterial diameter and FV allow for a noninvasive bedside estimation of CBF. CPAP at 15 cm H(2)O significantly decreased CBF in healthy awake volunteers. This effect appeared to be mediated predominately through the hypocapnic vasoconstriction coinciding with PCO(2) level reduction. The results suggest that CPAP should be used cautiously in patients with unstable cerebral hemodynamics.

Estimation of the correlation structure of crustal velocity heterogeneity from seismic reflection data

Numerous sources of evidence point to the fact that heterogeneity within the Earth's deep crystalline crust is complex and hence may be best described through stochastic rather than deterministic approaches. As seismic reflection imaging arguably offers the best means of sampling deep crustal rocks in situ, much interest has been expressed in using such data to characterize the stochastic nature of crustal heterogeneity. Previous work on this problem has shown that the spatial statistics of seismic reflection data are indeed related to those of the underlying heterogeneous seismic velocity distribution. As of yet, however, the nature of this relationship has remained elusive due to the fact that most of the work was either strictly empirical or based on incorrect methodological approaches. Here, we introduce a conceptual model, based on the assumption of weak scattering, that allows us to quantitatively link the second-order statistics of a 2-D seismic velocity distribution with those of the corresponding processed and depth-migrated seismic reflection image. We then perform a sensitivity study in order to investigate what information regarding the stochastic model parameters describing crustal velocity heterogeneity might potentially be recovered from the statistics of a seismic reflection image using this model. Finally, we present a Monte Carlo inversion strategy to estimate these parameters and we show examples of its application at two different source frequencies and using two different sets of prior information. Our results indicate that the inverse problem is inherently non-unique and that many different combinations of the vertical and lateral correlation lengths describing the velocity heterogeneity can yield seismic images with the same 2-D autocorrelation structure. The ratio of all of these possible combinations of vertical and lateral correlation lengths, however, remains roughly constant which indicates that, without additional prior information, the aspect ratio is the only parameter describing the stochastic seismic velocity structure that can be reliably recovered.

Extinction debt of high-mountain plants under 21st-century climate change

Quantitative estimates of the range loss of mountain plants under climate change have so far mostly relied on static geographical projections of species' habitat shifts(1-3). Here, we use a hybrid model(4) that combines such projections with simulations of demography and seed dispersal to forecast the climate-driven spatio-temporal dynamics of 150 high-mountain plant species across the European Alps. This model predicts average range size reductions of 44-50% by the end of the twenty-first century, which is similar to projections from the most 'optimistic' static model (49%). However, the hybrid model also indicates that population dynamics will lag behind climatic trends and that an average of 40% of the range still occupied at the end of the twenty-first century will have become climatically unsuitable for the respective species, creating an extinction debt(5,6). Alarmingly, species endemic to the Alps seem to face the highest range losses. These results caution against optimistic conclusions from moderate range size reductions observed during the twenty-first century as they are likely to belie more severe longer-term effects of climate warming on mountain plants.

Sperm velocity in an Alpine whitefish: effects of age, size, condition, fluctuating asymmetry and gonad abnormalities

The relationship between sperm velocity and individual age, size, body condition and fluctuating asymmetry was investigated in Alpine whitefish Coregonus fatioi. The fish analysed belonged to one among several sympatric whitefish populations of Lake Thun, Switzerland, which are characterized by a high prevalence of gonad alterations. Therefore, sperm velocity data were also tested for a link between gonad deformation and sperm swimming speed. Sperm velocity was significantly lower in larger-grown individuals and in individuals of higher body condition. As expected, sperm velocity was higher in males with higher levels of fluctuating asymmetry, but it did not significantly vary with male age. Moreover, variation in sperm velocity was found to be significantly higher in individuals showing some types of gonad alterations but it did not significantly correlate with the presence of other types of alterations. (C) 2007 The Authors Journal compilation (C) 2007 The Fisheries Society of the British Isles.

Early change in maladaptive defense style and development of the therapeutic alliance

Résumé Cette étude examine les changements précoces dans le Style Défensif Maladaptatif (SDM), le développement de l'alliance thérapeutique et la relation entre le SDM et l'alliance au cours d'une psychothérapie psychodynamique ultra-brève. Soixante-huit patients ambulatoires du centre de consultation psychiatrique et psychothérapique ont bénéficié d'une intervention psychodynamique en quatre séances. Les mesures des défenses et de l'alliance étaient effectuées à la première et à la dernière séance. Les patients qui ont débuté l'intervention avec une alliance faible et qui l'ont terminée avec une alliance haute (groupe de patients avec une alliance de croissance linéaire) ont diminué leur utilisation de défenses maladaptatives de manière significative au cours de la thérapie, alors que ce n'a pas été le cas pour les patients des groupes à alliances haute-stable et basse-stable. Les résultats ont montré qu'à la fin de l'intervention, le SDM et l'alliance étaient corrélés pour tous les patients. Cette corrélation intéressait plus particulièrement le groupe avec une alliance de croissance linéaire. Ces résultats suggèrent, que le développement de l'alliance thérapeutique reflètent le travail de collaboration entre le patient et son thérapeute alors qu'ils essayent de mieux comprendre les causes de la crise du patient. Cette compréhension peut aider à réduire les défenses initialement activées pour permettre au patient de se défendre de l'anxiété et d'un sentiment de détresse. Abstract This study examined the early change in Maladaptive Defense Style (MDS), the development of the Therapeutic Alliance, and the relationship between MDS and alliance, in a short psychodynamic intervention. Sixty-eight outpatients from a psychiatric clinic completed a four-session psychodynamic intervention. Defense and alliance measures were collected at the intake and the final session. Patients who began the intervention with a poor alliance but ended with a good alliance (linear growth therapeutic alliance group) significantly decreased their use of maladaptive defenses over the course of therapy, while patients in the high and low alliance groups did not. Results showed that at the end of the intervention, MDS and alliance were related across all patients. This relation concerned particularly the linear growth therapeutic alliance profile. These results suggest that the developing therapeutic alliance might reflect the collaborative work between the patient and the therapist as they try to understand the causes of the crisis. This understanding might help reduce maladaptive defenses that were initially activated to ward off anxiety and distress.

Stages of change, barriers, benefits, and preferences for exercise in RA patients: a cross-sectional study.

OBJECTIVES: To determine the distribution of exercise stages of change in a rheumatoid arthritis (RA) cohort, and to examine patients' perceptions of exercise benefits, barriers, and their preferences for exercise. METHODS: One hundred and twenty RA patients who attended the Rheumatology Unit of a University Hospital were asked to participate in the study. Those who agreed were administered a questionnaire to determine their exercise stage of change, their perceived benefits and barriers to exercise, and their preferences for various features of exercise. RESULTS: Eighty-nine (74%) patients were finally included in the analyses. Their mean age was 58.4 years, mean RA duration 10.1 years, and mean disease activity score 2.8. The distribution of exercise stages of change was as follows: precontemplation (n = 30, 34%), contemplation (n = 11, 13%), preparation (n = 5, 6%), action (n = 2, 2%), and maintenance (n = 39, 45%). Compared to patients in the maintenance stage of change, precontemplators exhibited different demographic and functional characteristics and reported less exercise benefits and more barriers to exercise. Most participants preferred exercising alone (40%), at home (29%), at a moderate intensity (64%), with advice provided by a rheumatologist (34%) or a specialist in exercise and RA (34%). Walking was by far the preferred type of exercise, in both the summer (86%) and the winter (51%). CONCLUSIONS: Our cohort of patients with RA was essentially distributed across the precontemplation and maintenance exercise stages of change. These subgroups of patients exhibit psychological and functional differences that make their needs different in terms of exercise counselling.

Simulation of interface dynamics across pore discontinuities

We study the dynamics of a water-oil meniscus moving from a smaller to a larger pore. The process is characterised by an abrupt change in the configuration, yielding a sudden energy release. A theoretic study for static conditions provides analytical solutions of the surface energy content of the system. Although the configuration after the sudden energy release is energetically more convenient, an energy barrier must be overcome before the process can happen spontaneously. The energy barrier depends on the system geometry and on the flow parameters. The analytical results are compared to numerical simulations that solve the full Navier-Stokes equation in the pore space and employ the Volume Of Fluid (VOF) method to track the evolution of the interface. First, the numerical simulations of a quasi-static process are validated by comparison with the analytical solutions for a static meniscus, then numerical simulations with varying injection velocity are used to investigate dynamic effects on the configuration change. During the sudden energy jump the system exhibits an oscillatory behaviour. Extension to more complex geometries might elucidate the mechanisms leading to a dynamic capillary pressure and to bifurcations in final distributions of fluid phases in porous

Se former à l'entretien motivationnel pour aider ses patients à changer: Les étudiants en médecine le confirment [Learning motivational interviewing to help patients change their health-related behaviors: medical students confirm it].

At the Lausanne University, 5th year medical students were trained in Motivational interviewing (MI). Eight hours of training improved their competence in the use of this approach. This experience supports the implementation of MI training in medical schools. Motivational interviewing allows the health professional to actively involve the patient in this behavior change process (drinking, smoking, diet, exercise, medication adherence, etc.), by encouraging reflection and reinforcing personal motivation and resources.

Early parenteral lipids and growth velocity in extremely-low-birth-weight infants.

BACKGROUND & AIMS: Whether early parenteral lipids improve postnatal growth of preterm neonates remains unclear. We aimed to assess the effects of parenteral lipids on growth velocity in extremely-low-birth-weight infants. METHODS: This retrospective cohort study included 121 extremely-low-birth-weight infants. The associations between parenteral lipids (cumulative intakes during the first week and delays in their introduction) and growth velocities (weight, head circumference and length) up to 28 days of life and to 36 weeks of corrected age were analysed using uni- and multivariate linear regression. RESULTS: Univariate analyses showed a significant positive association between the cumulative intakes of parenteral lipids during the first week and i) weight gain up to day 28; ii) weight gain up to 36 weeks of corrected age; iii) head circumference growth up to day 28. There was a negative correlation between the delay in parenteral lipid introduction and weight gain up to day 28. In multivariate analyses, the association between the cumulative intakes of parenteral lipids and weight gain up to 28 days was independent of gestational age at birth, birth weight, sex, smallness for gestational age, and enteral intakes (regression coefficient: 0.19; 95% CI: 0.01-0.38) and, up to 36 weeks, independent of gestational age, birth weight, sex, smallness for gestational age and parenteral glucose and amino acids (0.16; 95% CI: 0.04-0.27). CONCLUSIONS: Parenteral lipids during the first week were positively associated with weight gain in extremely-low-birth-weight infants and could improve early nutritional support of preterm neonates.

Predicting future distributions of mountain plants under climate change: does dispersal capacity matter?

Many studies have investigated the impacts that climate change could potentially have on the distribution of plant species, but few have attempted to constrain projections through plant dispersal limitations. Instead, most studies published so far have been using the simplification of considering dispersal as either unlimited or null. However, depending on a species' dispersal capacity, landscape fragmentation, and the rate of climatic change, these assumptions can lead to serious over- or underestimation of a species' future distribution. To quantify the discrepancies between unlimited, realistic, and no dispersal scenarios, we carried out projections of future distribution over the 21st century for 287 mountain plant species in a study area of the Western Swiss Alps. For each species, simulations were run for four dispersal scenarios (unlimited dispersal, no dispersal, realistic dispersal and realistic dispersal with long-distance dispersal events) and under four climate change scenarios. Although simulations accounting for realistic dispersal limitations did significantly differ from those considering dispersal as unlimited or null in terms of projected future distribution, using the unlimited dispersal simplification nevertheless provided good approximations for species extinctions under more moderate climate change scenarios. Overall, simulations accounting for dispersal limitations produced, for our mountainous study area, results that were significantly closer to unlimited dispersal than to no dispersal. Finally, analyzing the temporal pattern of species extinctions over the entire 21st century showed that, due to the possibility of a large number of species shifting their distribution to higher elevation, important species extinctions for our study area might not occur before the 2080-2100 time periods.