Spatiotemporal adaptive multiscale multiphysics simulations of two-phase flow

We present a spatiotemporal adaptive multiscale algorithm, which is based on the Multiscale Finite Volume method. The algorithm offers a very efficient framework to deal with multiphysics problems and to couple regions with different spatial resolution. We employ the method to simulate two-phase flow through porous media. At the fine scale, we consider a pore-scale description of the flow based on the Volume Of Fluid method. In order to construct a global problem that describes the coarse-scale behavior, the equations are averaged numerically with respect to auxiliary control volumes, and a Darcy-like coarse-scale model is obtained. The space adaptivity is based on the idea that a fine-scale description is only required in the front region, whereas the resolution can be coarsened elsewhere. Temporal adaptivity relies on the fact that the fine-scale and the coarse-scale problems can be solved with different temporal resolution (longer time steps can be used at the coarse scale). By simulating drainage under unstable flow conditions, we show that the method is able to capture the coarse-scale behavior outside the front region and to reproduce complex fluid patterns in the front region.

Elemental (C/N ratios) and isotopic (δ15N(org), δ13C(org)) compositions of sedimentary organic matter from a high-altitude mountain lake (Meidsee, 2661 m a.s.l., Switzerland): Implications for Lateglacial and Holocene Alpine landscape evolution

The deposition of Late Pleistocene and Holocene sediments in the high-altitude lake Meidsee (located at an altitude of 2661 m a.s.l. in the Southwestern Alps) strikingly coincided with global ice-sheet and mountain-glacier decay in the Alpine forelands and the formation of perialpine lakes. Radiocarbon ages of bottom-core sediments point out (pre-) Holocene ice retreat below 2700 m a.s.l., at about 16, 13, 10, and 9 cal. kyr BP. The Meidsee sedimentary record therefore provides information about the high-altitude Alpine landscape evolution since the Late Pleistocene/Holocene deglaciation in the Swiss Southwestern Alps. Prior to 5 cal. kyr BP, the C/N ratio and the isotopic composition of sedimentary organic matter (delta N-15(org), delta C-13(org)) indicate the deposition of algal-derived organic matter with limited input of terrestrial organic matter. The early Holocene and the Holocene climatic optimum (between 7.0 and 5.5 cal. kyr BP) were characterized by low erosion (decreasing magnetic susceptibility, chi) and high content of organic matter (C-org > 13 wt.%), enriched in C-13(org) (>-18 parts per thousand) with a low C/N (similar to 10) ratio, typical of modern algal matter derived from in situ production. During the late Holocene, there was a long-term increasing contribution of terrestrial organic matter into the lake (C/N > 11), with maxima between 2.4 and 0.9 cal. kyr BP. A major environmental change took place 800 years ago, with an abrupt decrease in the relative contribution of terrestrial organic material into the lake compared with aquatic organic material which subsequently largely dominated (C/N drop from 16 to 10). Nonetheless, this event was marked by a rise in soil erosion (chi), in nutrients input (N and P contents) and in anthropogenic lead deposition, suggesting a human disturbance of Alpine ecosystems 800 years ago. Indeed, this time period coincided with the migration of the Walser Alemannic people in the region, who settled at relatively high altitude in the Southwestern Alps for farming and maintaining Alpine passes.

A matter of time: delayed mate encounter postpones mating window initiation and reduces the strength of female choosiness

Reproductive success is determined by the presence and timing of encounter of mates. The latter depends on species-specific reproductive characteristics (e.g. initiation/duration of the mating window), season, and reproductive strategies (e.g. intensity of choosiness) that may potentially mitigate constraints imposed by mating windows. Despite their potentially crucial role for fitness and population dynamics, limited evidence exists about mating window initiation, duration and reproductive strategies. Here, we experimentally tested the mechanisms of initiation and the duration of the common lizard's Zootoca vivipara mating window, by manipulating the timing of mate encounter and analyzing its effect on (re-)mating probability. We furthermore tested treatment effects on female reproductive strategies, by measuring female choosiness. The timing of mate encounter and season did not significantly affect mating probability. However, a longer delay until mate encounter reduced female choosiness. Re-mating probability decreased with re-mating delay and was independent of mating delay. This indicates that mating window initiation depends on mate encounter, that its duration is fixed, and that plastic reproductive strategies exist. These findings contrast with previous beliefs and shows that mating windows per se may not necessarily constrain reproductive success, which is congruent with rapid range expansion and absence of positive density-effects on reproductive success (Allee effects). In summary, our results show that predicting the effect of mating windows on reproduction is complex and that experimental evidence is essential for evaluating their effect on reproduction and reproductive strategies, both being important determinants of population dynamics and the colonization of new habitats.

New pool of cortical interneuron precursors in the early postnatal dorsal white matter.

The migration of cortical γ-aminobutyric acidergic interneurons has been extensively studied in rodent embryos, whereas few studies have documented their postnatal migration. Combining in vivo analysis together with time-lapse imaging on cortical slices, we explored the origin and migration of cortical interneurons during the first weeks of postnatal life. Strikingly, we observed that a large pool of GAD65-GFP-positive cells accumulate in the dorsal white matter region during the first postnatal week. Part of these cells divides and expresses the transcription factor paired box 6 indicating the presence of local transient amplifying precursors. The vast majority of these cells are immature interneurons expressing the neuronal marker doublecortin and partly the calcium-binding protein calretinin. Time-lapse imaging reveals that GAD65-GFP-positive neurons migrate from the white matter pool into the overlying anterior cingulate cortex (aCC). Some interneurons in the postnatal aCC express the same immature neuronal markers suggesting ongoing migration of calretinin-positive interneurons. Finally, bromodeoxyuridine incorporation experiments confirm that a small fraction of interneurons located in the aCC are generated during the early postnatal period. These results altogether reveal that at postnatal ages, the dorsal white matter contains a pool of interneuron precursors that divide and migrate into the aCC.

Short-term increase in particulate matter blunts nocturnal blood pressure dipping and daytime urinary sodium excretion

Short-term exposure to ambient particulate matter with aerodynamic diameters<10 µm were found to be positively associated with blood pressure. Yet, little information exists regarding the association between particles and circadian rhythm of blood pressure. Hence, we analyzed the association of exposure to particulate matter with aerodynamic diameters<10 µm on the day of examination and ≤7 days before with ambulatory blood pressure and with sodium excretion in 359 adults from the general population using multiple linear regression. After controlling for potential confounders, a 10-µg/m3 increase in particulate matter with aerodynamic diameters<10 µm levels was associated with nighttime systolic blood pressure (β=1.32 mm Hg 95% CI, 0.06-2.58 mm Hg at lag 0; P=0.04), nighttime diastolic blood pressure (0.72 mm Hg 95% CI, 0.03-1.42 mm Hg at lag 2; P=0.04), nocturnal systolic blood pressure dipping (-0.96 mm Hg 95% CI, -1.89 to -0.03 mm Hg at lag 0; P=0.044), and daytime urinary sodium excretion (-0.05 log-mmol/min 95% CI, -0.10 to -0.01 log-mmol/min at lag 0; P=0.027) but not with nighttime sodium excretion. The associations with blood pressure rapidly diminished with increasing lag days, and the associations with daytime sodium excretion were maximal with particulate matter with aerodynamic diameters<10 µm in exposures 2 to 5 days before. The associations of short-term increases in particulate matter with aerodynamic diameters<10 µm with higher nighttime blood pressure and blunted systolic blood pressure dipping were preceded by associations with reduced ability of the kidney to excrete sodium during daytime. The underlying mechanism linking air pollution to increased cardiovascular risk may include disturbed circadian rhythms of renal sodium handling and blood pressure.

Glucose metabolism, gray matter structure, and memory decline in subjective memory impairment.

OBJECTIVE: To identify biological evidence for Alzheimer disease (AD) in individuals with subjective memory impairment (SMI) and unimpaired cognitive performance and to investigate the longitudinal cognitive course in these subjects. METHOD: [¹⁸F]fluoro-2-deoxyglucose PET (FDG-PET) and structural MRI were acquired in 31 subjects with SMI and 56 controls. Cognitive follow-up testing was performed (average follow-up time: 35 months). Differences in baseline brain imaging data and in memory decline were assessed between both groups. Associations of memory decline with brain imaging data were tested. RESULTS: The SMI group showed hypometabolism in the right precuneus and hypermetabolism in the right medial temporal lobe. Gray matter volume was reduced in the right hippocampus in the SMI group. At follow-up, subjects with SMI showed a poorer performance than controls on measures of episodic memory. Longitudinal memory decline in the SMI group was associated with reduced glucose metabolism in the right precuneus at baseline. CONCLUSION: The cross-sectional difference in 2 independent neuroimaging modalities indicates early AD pathology in SMI. The poorer memory performance at follow-up and the association of reduced longitudinal memory performance with hypometabolism in the precuneus at baseline support the concept of SMI as the earliest manifestation of AD.

Développement et validation d'une méthode de calcul indépendant des unités moniteur basée sur les simulations Monte Carlo pour la vérification des traitements IMRT

Résumé : La radiothérapie par modulation d'intensité (IMRT) est une technique de traitement qui utilise des faisceaux dont la fluence de rayonnement est modulée. L'IMRT, largement utilisée dans les pays industrialisés, permet d'atteindre une meilleure homogénéité de la dose à l'intérieur du volume cible et de réduire la dose aux organes à risque. Une méthode usuelle pour réaliser pratiquement la modulation des faisceaux est de sommer de petits faisceaux (segments) qui ont la même incidence. Cette technique est appelée IMRT step-and-shoot. Dans le contexte clinique, il est nécessaire de vérifier les plans de traitement des patients avant la première irradiation. Cette question n'est toujours pas résolue de manière satisfaisante. En effet, un calcul indépendant des unités moniteur (représentatif de la pondération des chaque segment) ne peut pas être réalisé pour les traitements IMRT step-and-shoot, car les poids des segments ne sont pas connus à priori, mais calculés au moment de la planification inverse. Par ailleurs, la vérification des plans de traitement par comparaison avec des mesures prend du temps et ne restitue pas la géométrie exacte du traitement. Dans ce travail, une méthode indépendante de calcul des plans de traitement IMRT step-and-shoot est décrite. Cette méthode est basée sur le code Monte Carlo EGSnrc/BEAMnrc, dont la modélisation de la tête de l'accélérateur linéaire a été validée dans une large gamme de situations. Les segments d'un plan de traitement IMRT sont simulés individuellement dans la géométrie exacte du traitement. Ensuite, les distributions de dose sont converties en dose absorbée dans l'eau par unité moniteur. La dose totale du traitement dans chaque élément de volume du patient (voxel) peut être exprimée comme une équation matricielle linéaire des unités moniteur et de la dose par unité moniteur de chacun des faisceaux. La résolution de cette équation est effectuée par l'inversion d'une matrice à l'aide de l'algorithme dit Non-Negative Least Square fit (NNLS). L'ensemble des voxels contenus dans le volume patient ne pouvant être utilisés dans le calcul pour des raisons de limitations informatiques, plusieurs possibilités de sélection ont été testées. Le meilleur choix consiste à utiliser les voxels contenus dans le Volume Cible de Planification (PTV). La méthode proposée dans ce travail a été testée avec huit cas cliniques représentatifs des traitements habituels de radiothérapie. Les unités moniteur obtenues conduisent à des distributions de dose globale cliniquement équivalentes à celles issues du logiciel de planification des traitements. Ainsi, cette méthode indépendante de calcul des unités moniteur pour l'IMRT step-andshootest validée pour une utilisation clinique. Par analogie, il serait possible d'envisager d'appliquer une méthode similaire pour d'autres modalités de traitement comme par exemple la tomothérapie. Abstract : Intensity Modulated RadioTherapy (IMRT) is a treatment technique that uses modulated beam fluence. IMRT is now widespread in more advanced countries, due to its improvement of dose conformation around target volume, and its ability to lower doses to organs at risk in complex clinical cases. One way to carry out beam modulation is to sum smaller beams (beamlets) with the same incidence. This technique is called step-and-shoot IMRT. In a clinical context, it is necessary to verify treatment plans before the first irradiation. IMRT Plan verification is still an issue for this technique. Independent monitor unit calculation (representative of the weight of each beamlet) can indeed not be performed for IMRT step-and-shoot, because beamlet weights are not known a priori, but calculated by inverse planning. Besides, treatment plan verification by comparison with measured data is time consuming and performed in a simple geometry, usually in a cubic water phantom with all machine angles set to zero. In this work, an independent method for monitor unit calculation for step-and-shoot IMRT is described. This method is based on the Monte Carlo code EGSnrc/BEAMnrc. The Monte Carlo model of the head of the linear accelerator is validated by comparison of simulated and measured dose distributions in a large range of situations. The beamlets of an IMRT treatment plan are calculated individually by Monte Carlo, in the exact geometry of the treatment. Then, the dose distributions of the beamlets are converted in absorbed dose to water per monitor unit. The dose of the whole treatment in each volume element (voxel) can be expressed through a linear matrix equation of the monitor units and dose per monitor unit of every beamlets. This equation is solved by a Non-Negative Least Sqvare fif algorithm (NNLS). However, not every voxels inside the patient volume can be used in order to solve this equation, because of computer limitations. Several ways of voxel selection have been tested and the best choice consists in using voxels inside the Planning Target Volume (PTV). The method presented in this work was tested with eight clinical cases, which were representative of usual radiotherapy treatments. The monitor units obtained lead to clinically equivalent global dose distributions. Thus, this independent monitor unit calculation method for step-and-shoot IMRT is validated and can therefore be used in a clinical routine. It would be possible to consider applying a similar method for other treatment modalities, such as for instance tomotherapy or volumetric modulated arc therapy.

Multimodal MRI analysis of structural and functional networks in the human brain : application to two clinical studies : driving abilities under THC intoxication and early white matter changes in FXTAS

Les approches multimodales dans l'imagerie cérébrale non invasive sont de plus en plus considérées comme un outil indispensable pour la compréhension des différents aspects de la structure et de la fonction cérébrale. Grâce aux progrès des techniques d'acquisition des images de Resonance Magnetique et aux nouveaux outils pour le traitement des données, il est désormais possible de mesurer plusieurs paramètres sensibles aux différentes caractéristiques des tissues cérébraux. Ces progrès permettent, par exemple, d'étudier les substrats anatomiques qui sont à la base des processus cognitifs ou de discerner au niveau purement structurel les phénomènes dégénératifs et développementaux. Cette thèse met en évidence l'importance de l'utilisation d'une approche multimodale pour étudier les différents aspects de la dynamique cérébrale grâce à l'application de cette approche à deux études cliniques: l'évaluation structurelle et fonctionnelle des effets aigus du cannabis fumé chez des consommateurs réguliers et occasionnels, et l'évaluation de l'intégrité de la substance grise et blanche chez des jeunes porteurs de la prémutations du gène FMR1 à risque de développer le FXTAS (Fragile-X Tremor Ataxia Syndrome). Nous avons montré que chez les fumeurs occasionnels de cannabis, même à faible concentration du principal composant psychoactif (THC) dans le sang, la performance lors d'une tâche visuo-motrice est fortement diminuée, et qu'il y a des changements dans l'activité des trois réseaux cérébraux impliqués dans les processus cognitifs: le réseau de saillance, le réseau du contrôle exécutif, et le réseau actif par défaut (Default Mode). Les sujets ne sont pas en mesure de saisir les saillances dans l'environnement et de focaliser leur attention sur la tâche. L'augmentation de la réponse hémodynamique dans le cortex cingulaire antérieur suggère une augmentation de l'activité introspective. Une investigation des ef¬fets au niveau cérébral d'une exposition prolongée au cannabis, montre des changements persistants de la substance grise dans les régions associées à la mémoire et au traitement des émotions. Le niveau d'atrophie dans ces structures corrèle avec la consommation de cannabis au cours des trois mois précédant l'étude. Dans la deuxième étude, nous démontrons des altérations structurelles des décennies avant l'apparition du syndrome FXTAS chez des sujets jeunes, asymptomatiques, et porteurs de la prémutation du gène FMR1. Les modifications trouvées peuvent être liées à deux mécanismes différents. Les altérations dans le réseau moteur du cervelet et dans la fimbria de l'hippocampe, suggèrent un effet développemental de la prémutation. Elles incluent aussi une atrophie de la substance grise du lobule VI du cervelet et l'altération des propriétés tissulaires de la substance blanche des projections afférentes correspondantes aux pédoncules cérébelleux moyens. Les lésions diffuses de la substance blanche cérébrale peu¬vent être un marquer précoce du développement de la maladie, car elles sont liées à un phénomène dégénératif qui précède l'apparition des symptômes du FXTAS. - Multimodal brain imaging is becoming a leading tool for understanding different aspects of brain structure and function. Thanks to the advances in Magnetic Resonance imaging (MRI) acquisition schemes and data processing techniques, it is now possible to measure different parameters sensitive to different tissue characteristics. This allows for example to investigate anatomical substrates underlying cognitive processing, or to disentangle, at a pure structural level degeneration and developmental processes. This thesis highlights the importance of using a multimodal approach for investigating different aspects of brain dynamics by applying this approach to two clinical studies: functional and structural assessment of the acute effects of cannabis smoking in regular and occasional users, and grey and white matter assessment in young FMR1 premutation carriers at risk of developing FXTAS. We demonstrate that in occasional smokers cannabis smoking, even at low concentration of the main psychoactive component (THC) in the blood, strongly decrease subjects' performance on a visuo-motor tracking task, and globally alters the activity of the three brain networks involved in cognitive processing: the Salience, the Control Executive, and the Default Mode networks. Subjects are unable to capture saliences in the environment and to orient attention to the task; the increase in Hemodynamic Response in the Anterior Cingulate Cortex suggests an increase in self-oriented mental activity. A further investigation on long term exposure to cannabis, shows a persistent grey matter modification in brain regions associated with memory and affective processing. The degree of atrophy in these structures also correlates with the estimation of drug use in the three months prior the participation to the study. In the second study we demonstrate structural changes in young asymptomatic premutation carriers decades before the onset of FXTAS that might be related to two different mechanisms. Alteration of the cerebellar motor network and of the hippocampal fimbria/ fornix, may reflect a potential neurodevelopmental effect of the premutation. These include grey matter atrophy in lobule VI and modification of white matter tissue property in the corresponding afferent projections through the Middle Cerebellar Peduncles. Diffuse hemispheric white matter lesions that seem to appear closer to the onset of FXTAS and be related to a neurodegenerative phenomenon may mark the imminent onset of FXTAS.

Being Small for Gestational Age: Does it Matter for the Neurodevelopment of Premature Infants? A Cohort Study.

BACKGROUND: Whether being small for gestational age (SGA) increases the risk of adverse neurodevelopmental outcome in premature infants remains controversial. OBJECTIVE: to study the impact of SGA (birthweight < percentile 10) on cognition, behavior, neurodevelopmental impairment and use of therapy at 5 years old. METHODS: This population-based prospective cohort included infants born before 32 weeks of gestation. Cognition was evaluated with the K-ABC, and behavior with the Strengths and Difficulties Questionnaire (SDQ). Primary outcomes were cognitive and behavioral scores, as well as neurodevelopmental impairment (cognitive score < 2SD, hearing loss, blindness, or cerebral palsy). The need of therapy, an indirect indicator of neurodevelopmental impairment, was a secondary outcome. Linear and logistic regression models were used to analyze the association of SGA with neurodevelopment. RESULTS: 342/515 (76%) premature infants were assessed. SGA was significantly associated with hyperactivity scores of the SDQ (coefficient 0.81, p < 0.04), but not with cognitive scores, neurodevelopmental impairment or the need of therapy. Gestational age, socio-economic status, and major brain lesions were associated with cognitive outcome in the univariate and multivariate model, whereas asphyxia, sepsis and bronchopulmonary dysplasia were associated in the univariate model only. Severe impairment was associated with fetal tobacco exposition, asphyxia, gestational age and major brain lesions. Different neonatal factors were associated with the use of single or multiple therapies: children with one therapy were more likely to have suffered birth asphyxia or necrotizing enterocolitis, whereas the need for several therapies was predicted by major brain lesions. DISCUSSION: In this large cohort of premature infants, assessed at 5 years old with a complete panel of tests, SGA was associated with hyperactive behavior, but not with cognition, neurodevelopmental impairment or use of therapy. Birthweight <10th percentile alone does not appear to be an independent risk factor of neurodevelopmental adverse outcome in preterm children.

Conditioning of multiple-point statistics facies simulations to tomographic images

Geophysical tomography captures the spatial distribution of the underlying geophysical property at a relatively high resolution, but the tomographic images tend to be blurred representations of reality and generally fail to reproduce sharp interfaces. Such models may cause significant bias when taken as a basis for predictive flow and transport modeling and are unsuitable for uncertainty assessment. We present a methodology in which tomograms are used to condition multiple-point statistics (MPS) simulations. A large set of geologically reasonable facies realizations and their corresponding synthetically calculated cross-hole radar tomograms are used as a training image. The training image is scanned with a direct sampling algorithm for patterns in the conditioning tomogram, while accounting for the spatially varying resolution of the tomograms. In a post-processing step, only those conditional simulations that predicted the radar traveltimes within the expected data error levels are accepted. The methodology is demonstrated on a two-facies example featuring channels and an aquifer analog of alluvial sedimentary structures with five facies. For both cases, MPS simulations exhibit the sharp interfaces and the geological patterns found in the training image. Compared to unconditioned MPS simulations, the uncertainty in transport predictions is markedly decreased for simulations conditioned to tomograms. As an improvement to other approaches relying on classical smoothness-constrained geophysical tomography, the proposed method allows for: (1) reproduction of sharp interfaces, (2) incorporation of realistic geological constraints and (3) generation of multiple realizations that enables uncertainty assessment.