Turning sunlight into stone: the oxalate-carbonate pathway in a tropical tree ecosystem.

An African oxalogenic tree, the iroko tree (Milicia excelsa), has the property to enhance carbonate precipitation in tropical oxisols, where such accumulations are not expected due to the acidic conditions in these types of soils. This uncommon process is linked to the oxalate-carbonate pathway, which increases soil pH through oxalate oxidation. In order to investigate the oxalate-carbonate pathway in the iroko system, fluxes of matter have been identified, described, and evaluated from field to microscopic scales. In the first centimeters of the soil profile, decaying of the organic matter allows the release of whewellite crystals, mainly due to the action of termites and saprophytic fungi. In addition, a concomitant flux of carbonate formed in wood tissues contributes to the carbonate flux and is identified as a direct consequence of wood feeding by termites. Nevertheless, calcite biomineralization of the tree is not a consequence of in situ oxalate consumption, but rather related to the oxalate oxidation inside the upper part of the soil. The consequence of this oxidation is the presence of carbonate ions in the soil solution pumped through the roots, leading to preferential mineralization of the roots and the trunk base. An ideal scenario for the iroko biomineralization and soil carbonate accumulation starts with oxalatization: as the iroko tree grows, the organic matter flux to the soil constitutes the litter, and an oxalate pool is formed on the forest ground. Then, wood rotting agents (mainly termites, saprophytic fungi, and bacteria) release significant amounts of oxalate crystals from decaying plant tissues. In addition, some of these agents are themselves producers of oxalate (e.g. fungi). Both processes contribute to a soil pool of "available" oxalate crystals. Oxalate consumption by oxalotrophic bacteria can then start. Carbonate and calcium ions present in the soil solution represent the end products of the oxalate-carbonate pathway. The solution is pumped through the roots, leading to carbonate precipitation. The main pools of carbon are clearly identified as the organic matter (the tree and its organic products), the oxalate crystals, and the various carbonate features. A functional model based on field observations and diagenetic investigations with δ13C signatures of the various compartments involved in the local carbon cycle is proposed. It suggests that the iroko ecosystem can act as a long-term carbon sink, as long as the calcium source is related to non-carbonate rocks. Consequently, this carbon sink, driven by the oxalate carbonate pathway around an iroko tree, constitutes a true carbon trapping ecosystem as defined by ecological theory.

Transitions from reproductive systems governed by two self-incompatible Loci to one in fungi.

Self-incompatibility (SI), a reproductive system broadly present in plants, chordates, fungi, and protists, might be controlled by one or several multiallelic loci. How a transition in the number of SI loci can occur and the consequences of such events for the population's genetics and dynamics have not been studied theoretically. Here, we provide analytical descriptions of two transition mechanisms: linkage of the two SI loci (scenario 1) and the loss of function of one incompatibility gene within a mating type of a population with two SI loci (scenario 2). We show that invasion of populations by the new mating type form depends on whether the fitness of the new type is lowered, and on the allelic diversity of the SI loci and the recombination between the two SI loci in the starting population. Moreover, under scenario 1, it also depends on the frequency of the SI alleles that became linked. We demonstrate that, following invasion, complete transitions in the reproductive system occurs under scenario 2 and is predicted only for small populations under scenario 1. Interestingly, such events are associated with a drastic reduction in mating type number.

Effect of changes in the deuterium content of drinking water on the hydrogen isotope ratio of urinary steroids in the context of sports drug testing.

The hydrogen isotope ratio (HIR) of body water and, therefore, of all endogenously synthesized compounds in humans, is mainly affected by the HIR of ingested drinking water. As a consequence, the entire organism and all of its synthesized substrates will reflect alterations in the isotope ratio of drinking water, which depends on the duration of exposure. To investigate the effect of this change on endogenous urinary steroids relevant to doping-control analysis the hydrogen isotope composition of potable water was suddenly enriched from -50 to 200 0/00 and maintained at this level for two weeks for two individuals. The steroids under investigation were 5β-pregnane-3α,20α-diol, 5α-androst-16-en-3α-ol, 3α-hydroxy-5α-androstan-17-one (ANDRO), 3α-hydroxy-5β-androstan-17-one (ETIO), 5α-androstane-3α,17β-diol, and 5β-androstane-3α,17β-diol (excreted as glucuronides) and ETIO, ANDRO and 3β-hydroxyandrost-5-en-17-one (excreted as sulfates). The HIR of body water was estimated by determination of the HIR of total native urine, to trace the induced changes. The hydrogen in steroids is partly derived from the total amount of body water and cholesterol-enrichment could be calculated by use of these data. Although the sum of changes in the isotopic composition of body water was 150 0/00, shifts of approximately 30 0/00 were observed for urinary steroids. Parallel enrichment in their HIR was observed for most of the steroids, and none of the differences between the HIR of individual steroids was elevated beyond recently established thresholds. This finding is important to sports drug testing because it supports the intended use of this novel and complementary methodology even in cases where athletes have drunk water of different HIR, a plausible and, presumably, inevitable scenario while traveling.

Reconstructing the origins of high-alpine niches and cushion life form in the genus Androsace S.L. (Primulaceae).

Relatively, few species have been able to colonize extremely cold alpine environments. We investigate the role played by the cushion life form in the evolution of climatic niches in the plant genus Androsace s.l., which spreads across the mountain ranges of the Northern Hemisphere. Using robust methods that account for phylogenetic uncertainty, intraspecific variability of climatic requirements and different life-history evolution scenarios, we show that climatic niches of Androsace s.l. exhibit low phylogenetic signal and that they evolved relatively recently and punctually. Models of niche evolution fitted onto phylogenies show that the cushion life form has been a key innovation providing the opportunity to occupy extremely cold environments, thus contributing to rapid climatic niche diversification in the genus Androsace s.l. We then propose a plausible scenario for the adaptation of plants to alpine habitats.

Quantitative integration of geophysical and hydraulic data : from the local towards the regional scale range

Des progrès significatifs ont été réalisés dans le domaine de l'intégration quantitative des données géophysique et hydrologique l'échelle locale. Cependant, l'extension à de plus grandes échelles des approches correspondantes constitue encore un défi majeur. Il est néanmoins extrêmement important de relever ce défi pour développer des modèles fiables de flux des eaux souterraines et de transport de contaminant. Pour résoudre ce problème, j'ai développé une technique d'intégration des données hydrogéophysiques basée sur une procédure bayésienne de simulation séquentielle en deux étapes. Cette procédure vise des problèmes à plus grande échelle. L'objectif est de simuler la distribution d'un paramètre hydraulique cible à partir, d'une part, de mesures d'un paramètre géophysique pertinent qui couvrent l'espace de manière exhaustive, mais avec une faible résolution (spatiale) et, d'autre part, de mesures locales de très haute résolution des mêmes paramètres géophysique et hydraulique. Pour cela, mon algorithme lie dans un premier temps les données géophysiques de faible et de haute résolution à travers une procédure de réduction déchelle. Les données géophysiques régionales réduites sont ensuite reliées au champ du paramètre hydraulique à haute résolution. J'illustre d'abord l'application de cette nouvelle approche dintégration des données à une base de données synthétiques réaliste. Celle-ci est constituée de mesures de conductivité hydraulique et électrique de haute résolution réalisées dans les mêmes forages ainsi que destimations des conductivités électriques obtenues à partir de mesures de tomographic de résistivité électrique (ERT) sur l'ensemble de l'espace. Ces dernières mesures ont une faible résolution spatiale. La viabilité globale de cette méthode est testée en effectuant les simulations de flux et de transport au travers du modèle original du champ de conductivité hydraulique ainsi que du modèle simulé. Les simulations sont alors comparées. Les résultats obtenus indiquent que la procédure dintégration des données proposée permet d'obtenir des estimations de la conductivité en adéquation avec la structure à grande échelle ainsi que des predictions fiables des caractéristiques de transports sur des distances de moyenne à grande échelle. Les résultats correspondant au scénario de terrain indiquent que l'approche d'intégration des données nouvellement mise au point est capable d'appréhender correctement les hétérogénéitées à petite échelle aussi bien que les tendances à gande échelle du champ hydraulique prévalent. Les résultats montrent également une flexibilté remarquable et une robustesse de cette nouvelle approche dintégration des données. De ce fait, elle est susceptible d'être appliquée à un large éventail de données géophysiques et hydrologiques, à toutes les gammes déchelles. Dans la deuxième partie de ma thèse, j'évalue en détail la viabilité du réechantillonnage geostatique séquentiel comme mécanisme de proposition pour les méthodes Markov Chain Monte Carlo (MCMC) appliquées à des probmes inverses géophysiques et hydrologiques de grande dimension . L'objectif est de permettre une quantification plus précise et plus réaliste des incertitudes associées aux modèles obtenus. En considérant une série dexemples de tomographic radar puits à puits, j'étudie deux classes de stratégies de rééchantillonnage spatial en considérant leur habilité à générer efficacement et précisément des réalisations de la distribution postérieure bayésienne. Les résultats obtenus montrent que, malgré sa popularité, le réechantillonnage séquentiel est plutôt inefficace à générer des échantillons postérieurs indépendants pour des études de cas synthétiques réalistes, notamment pour le cas assez communs et importants où il existe de fortes corrélations spatiales entre le modèle et les paramètres. Pour résoudre ce problème, j'ai développé un nouvelle approche de perturbation basée sur une déformation progressive. Cette approche est flexible en ce qui concerne le nombre de paramètres du modèle et lintensité de la perturbation. Par rapport au rééchantillonage séquentiel, cette nouvelle approche s'avère être très efficace pour diminuer le nombre requis d'itérations pour générer des échantillons indépendants à partir de la distribution postérieure bayésienne. - Significant progress has been made with regard to the quantitative integration of geophysical and hydrological data at the local scale. However, extending corresponding approaches beyond the local scale still represents a major challenge, yet is critically important for the development of reliable groundwater flow and contaminant transport models. To address this issue, I have developed a hydrogeophysical data integration technique based on a two-step Bayesian sequential simulation procedure that is specifically targeted towards larger-scale problems. The objective is to simulate the distribution of a target hydraulic parameter based on spatially exhaustive, but poorly resolved, measurements of a pertinent geophysical parameter and locally highly resolved, but spatially sparse, measurements of the considered geophysical and hydraulic parameters. To this end, my algorithm links the low- and high-resolution geophysical data via a downscaling procedure before relating the downscaled regional-scale geophysical data to the high-resolution hydraulic parameter field. I first illustrate the application of this novel data integration approach to a realistic synthetic database consisting of collocated high-resolution borehole measurements of the hydraulic and electrical conductivities and spatially exhaustive, low-resolution electrical conductivity estimates obtained from electrical resistivity tomography (ERT). The overall viability of this method is tested and verified by performing and comparing flow and transport simulations through the original and simulated hydraulic conductivity fields. The corresponding results indicate that the proposed data integration procedure does indeed allow for obtaining faithful estimates of the larger-scale hydraulic conductivity structure and reliable predictions of the transport characteristics over medium- to regional-scale distances. The approach is then applied to a corresponding field scenario consisting of collocated high- resolution measurements of the electrical conductivity, as measured using a cone penetrometer testing (CPT) system, and the hydraulic conductivity, as estimated from electromagnetic flowmeter and slug test measurements, in combination with spatially exhaustive low-resolution electrical conductivity estimates obtained from surface-based electrical resistivity tomography (ERT). The corresponding results indicate that the newly developed data integration approach is indeed capable of adequately capturing both the small-scale heterogeneity as well as the larger-scale trend of the prevailing hydraulic conductivity field. The results also indicate that this novel data integration approach is remarkably flexible and robust and hence can be expected to be applicable to a wide range of geophysical and hydrological data at all scale ranges. In the second part of my thesis, I evaluate in detail the viability of sequential geostatistical resampling as a proposal mechanism for Markov Chain Monte Carlo (MCMC) methods applied to high-dimensional geophysical and hydrological inverse problems in order to allow for a more accurate and realistic quantification of the uncertainty associated with the thus inferred models. Focusing on a series of pertinent crosshole georadar tomographic examples, I investigated two classes of geostatistical resampling strategies with regard to their ability to efficiently and accurately generate independent realizations from the Bayesian posterior distribution. The corresponding results indicate that, despite its popularity, sequential resampling is rather inefficient at drawing independent posterior samples for realistic synthetic case studies, notably for the practically common and important scenario of pronounced spatial correlation between model parameters. To address this issue, I have developed a new gradual-deformation-based perturbation approach, which is flexible with regard to the number of model parameters as well as the perturbation strength. Compared to sequential resampling, this newly proposed approach was proven to be highly effective in decreasing the number of iterations required for drawing independent samples from the Bayesian posterior distribution.

Computational quest for understanding the role of astrocyte signaling in synaptic transmission and plasticity.

The complexity of the signaling network that underlies astrocyte-synapse interactions may seem discouraging when tackled from a theoretical perspective. Computational modeling is challenged by the fact that many details remain hitherto unknown and conventional approaches to describe synaptic function are unsuitable to explain experimental observations when astrocytic signaling is taken into account. Supported by experimental evidence is the possibility that astrocytes perform genuine information processing by means of their calcium signaling and are players in the physiological setting of the basal tone of synaptic transmission. Here we consider the plausibility of this scenario from a theoretical perspective, focusing on the modulation of synaptic release probability by the astrocyte and its implications on synaptic plasticity. The analysis of the signaling pathways underlying such modulation refines our notion of tripartite synapse and has profound implications on our understanding of brain function.

Future hotspots of terrestrial mammal loss.

Current levels of endangerment and historical trends of species and habitats are the main criteria used to direct conservation efforts globally. Estimates of future declines, which might indicate different priorities than past declines, have been limited by the lack of appropriate data and models. Given that much of conservation is about anticipating and responding to future threats, our inability to look forward at a global scale has been a major constraint on effective action. Here, we assess the geography and extent of projected future changes in suitable habitat for terrestrial mammals within their present ranges. We used a global earth-system model, IMAGE, coupled with fine-scale habitat suitability models and parametrized according to four global scenarios of human development. We identified the most affected countries by 2050 for each scenario, assuming that no additional conservation actions other than those described in the scenarios take place. We found that, with some exceptions, most of the countries with the largest predicted losses of suitable habitat for mammals are in Africa and the Americas. African and North American countries were also predicted to host the most species with large proportional global declines. Most of the countries we identified as future hotspots of terrestrial mammal loss have little or no overlap with the present global conservation priorities, thus confirming the need for forward-looking analyses in conservation priority setting. The expected growth in human populations and consumption in hotspots of future mammal loss mean that local conservation actions such as protected areas might not be sufficient to mitigate losses. Other policies, directed towards the root causes of biodiversity loss, are required, both in Africa and other parts of the world.

Genetic structure of Gymnures (genus Hylomys; Erinaceidae) on continental islands of Southeast Asia: Historical effects of fragmentation

Eustatic sea level changes during Pleistocene climatic fluctuations produced several cycles of connection-isolation among continental islands of the Sunda shelf. To explore the potential effects of these fluctuations, we reconstructed a model of the vicariant events that separated these islands, based on bathymetric information. Among many possible scenarios, two opposite phylogenetic patterns of evolution were predicted for terrestrial organisms living in this region: one is based on the classical allopatric speciation mode of evolution, while the other is the outcome of a sequential dispersal colonization of the archipelago. We tested the applicability of these predictions with an analysis of sequence variation of the cytochrome b gene from several taxa of Hylomys. They were sampled throughout SE-Asia and the Sunda islands. High levels of haplotype differentiation characterize the different island taxa. Such levels of differentiation support the existence of several allopatric species, as was suggested by previous allozyme and morphological data. Also in accordance with previous results, the occurrence of two sympatric species from Sumatra is suggested by their strongly divergent haplotypes. One species, Hylomys suillus maxi, is found both on Sumatra and in Peninsular Malaysia, while the other, H. parvus, is endemic to Sumatra. Its closest relative is H. suillus dorsalis from Borneo. Phylogenetic reconstructions also demonstrate the existence of a Sundaic clade composed of all island taxa, as opposed to those from the continent. Although there is no statistical support for either proposed biogeographic model of evolution, we argue that the sequential dispersal scenario is more appropriate to describe the genetic variation found among the Hylomys taxa. However, despite strong differentiation among island haplotypes, the cladistic relationships between some island taxa could not be resolved. We argue that this is evidence of a rapid radiation, suggesting that the separation of the islands may have been perceived as a simultaneous event rather than as a succession of vicariant events. Furthermore, the estimates of divergence times between the haplotypes of these taxa suggest that this radiation may actually have predated the climatic fluctuations of the Pleistocene. Further refinement of the initial palaeogeographic models of evolution are therefore needed to account for these results.

The robustness of the weak selection approximation for the evolution of altruism against strong selection.

The weak selection approximation of population genetics has made possible the analysis of social evolution under a considerable variety of biological scenarios. Despite its extensive usage, the accuracy of weak selection in predicting the emergence of altruism under limited dispersal when selection intensity increases remains unclear. Here, we derive the condition for the spread of an altruistic mutant in the infinite island model of dispersal under a Moran reproductive process and arbitrary strength of selection. The simplicity of the model allows us to compare weak and strong selection regimes analytically. Our results demonstrate that the weak selection approximation is robust to moderate increases in selection intensity and therefore provides a good approximation to understand the invasion of altruism in spatially structured population. In particular, we find that the weak selection approximation is excellent even if selection is very strong, when either migration is much stronger than selection or when patches are large. Importantly, we emphasize that the weak selection approximation provides the ideal condition for the invasion of altruism, and increasing selection intensity will impede the emergence of altruism. We discuss that this should also hold for more complicated life cycles and for culturally transmitted altruism. Using the weak selection approximation is therefore unlikely to miss out on any demographic scenario that lead to the evolution of altruism under limited dispersal.

10

Research on individual social policy preferences has highlighted a number of socio-structural cleavages as determinants. Studies investigating public opinion on the various redistributive schemes that make up today's welfare states have shown the relevance of class-related factors such as income or education as key explanatory variables (Ferrera 1993; Taylor-Gooby 1995, 1998; and Svallfors 1997). More recent studies, however, have suggested that other factors are also likely to play a role. Among these, the most important are age, gender, and individual values (Armingeon 2006; Deitch 2004; and Roller 2000, 2002). The scenario that emerges from the existing literature is one of multiple intersecting cleavages, but it remains unclear as to what today is the relative weight and specific impact of each of these cleavages.

A deterministic model to assess the impact of AIDS interventions in countries with generalized HIV epidemic. An application to Botswana

General introductionThe Human Immunodeficiency/Acquired Immunodeficiency Syndrome (HIV/AIDS) epidemic, despite recent encouraging announcements by the World Health Organization (WHO) is still today one of the world's major health care challenges.The present work lies in the field of health care management, in particular, we aim to evaluate the behavioural and non-behavioural interventions against HIV/AIDS in developing countries through a deterministic simulation model, both in human and economic terms. We will focus on assessing the effectiveness of the antiretroviral therapies (ART) in heterosexual populations living in lesser developed countries where the epidemic has generalized (formerly defined by the WHO as type II countries). The model is calibrated using Botswana as a case study, however our model can be adapted to other countries with similar transmission dynamics.The first part of this thesis consists of reviewing the main mathematical concepts describing the transmission of infectious agents in general but with a focus on human immunodeficiency virus (HIV) transmission. We also review deterministic models assessing HIV interventions with a focus on models aimed at African countries. This review helps us to recognize the need for a generic model and allows us to define a typical structure of such a generic deterministic model.The second part describes the main feed-back loops underlying the dynamics of HIV transmission. These loops represent the foundation of our model. This part also provides a detailed description of the model, including the various infected and non-infected population groups, the type of sexual relationships, the infection matrices, important factors impacting HIV transmission such as condom use, other sexually transmitted diseases (STD) and male circumcision. We also included in the model a dynamic life expectancy calculator which, to our knowledge, is a unique feature allowing more realistic cost-efficiency calculations. Various intervention scenarios are evaluated using the model, each of them including ART in combination with other interventions, namely: circumcision, campaigns aimed at behavioral change (Abstain, Be faithful or use Condoms also named ABC campaigns), and treatment of other STD. A cost efficiency analysis (CEA) is performed for each scenario. The CEA consists of measuring the cost per disability-adjusted life year (DALY) averted. This part also describes the model calibration and validation, including a sensitivity analysis.The third part reports the results and discusses the model limitations. In particular, we argue that the combination of ART and ABC campaigns and ART and treatment of other STDs are the most cost-efficient interventions through 2020. The main model limitations include modeling the complexity of sexual relationships, omission of international migration and ignoring variability in infectiousness according to the AIDS stage.The fourth part reviews the major contributions of the thesis and discusses model generalizability and flexibility. Finally, we conclude that by selecting the adequate interventions mix, policy makers can significantly reduce the adult prevalence in Botswana in the coming twenty years providing the country and its donors can bear the cost involved.Part I: Context and literature reviewIn this section, after a brief introduction to the general literature we focus in section two on the key mathematical concepts describing the transmission of infectious agents in general with a focus on HIV transmission. Section three provides a description of HIV policy models, with a focus on deterministic models. This leads us in section four to envision the need for a generic deterministic HIV policy model and briefly describe the structure of such a generic model applicable to countries with generalized HIV/AIDS epidemic, also defined as pattern II countries by the WHO.

The balanced lethal system of crested newts: a ghost of sex chromosomes past?

Balanced lethal systems are more than biological curiosities: as theory predicts, they should quickly be eliminated through the joint forces of recombination and selection. That such systems might become fixed in natural populations poses a challenge to evolutionary theory. Here we address the case of a balanced lethal system fixed in crested newts and related species, which makes 50% of offspring die early in development. All adults are heteromorphic for chromosome pair 1. The two homologues (1A and 1B) have different recessive deleterious alleles fixed on a nonrecombining segment, so that heterozygotes are viable, while homozygotes are lethal. Given such a strong segregation load, how could autosomes stop recombining? We propose a role for a sex-chromosome turnover from pair 1 (putative ancestral sex chromosome) to pair 4 (currently active sex chromosome). Accordingly, 1A and 1B represent two variants (Y(A) and Y(B)) of the Y chromosome from an ancestral male-heterogametic system. We formalize a scenario in which turnovers are driven by sex ratio selection stemming from gene-environment interactions on sex determination. Individual-based simulations show that a balanced lethal system can be fixed with significant likelihood, provided the masculinizing allele on chromosome 4 appears after the elimination of the feminizing allele on chromosome 1. Our study illustrates how strikingly maladaptive traits might evolve through natural selection.

Lung infections: the radiologist's perspective.

Imaging plays a key role in lung infections. A CT scan must be carried out when there is a strong clinical suspicion of pneumonia that is accompanied by normal, ambiguous, or nonspecific radiography, a scenario that occurs most commonly in immunocompromised patients. CT allows clinicians to detect associated abnormalities or an underlying condition and it can guide bronchoalveolar lavage or a percutaneous or transbronchial lung biopsy. An organism can vary in how it is expressed depending on the extent to which the patient is immunocompromised. This is seen in tuberculosis in patients with AIDS. The infective agents vary with the type of immune deficiency and some infections can quickly become life-threatening. Clinicians should be aware of the complex radiological spectrum of pulmonary aspergillosis, given that this diagnosis must be considered in specific settings.

The demographic benefits of belligerence and bravery: defeated group repopulation or victorious group size expansion?

Intraspecific coalitional aggression between groups of individuals is a widespread trait in the animal world. It occurs in invertebrates and vertebrates, and is prevalent in humans. What are the conditions under which coalitional aggression evolves in natural populations? In this article, I develop a mathematical model delineating conditions where natural selection can favor the coevolution of belligerence and bravery between small-scale societies. Belligerence increases an actor's group probability of trying to conquer another group and bravery increase the actors's group probability of defeating an attacked group. The model takes into account two different types of demographic scenarios that may lead to the coevolution of belligerence and bravery. Under the first, the fitness benefits driving the coevolution of belligerence and bravery come through the repopulation of defeated groups by fission of victorious ones. Under the second demographic scenario, the fitness benefits come through a temporary increase in the local carrying capacity of victorious groups, after transfer of resources from defeated groups to victorious ones. The analysis of the model suggests that the selective pressures on belligerence and bravery are stronger when defeated groups can be repopulated by victorious ones. The analysis also suggests that, depending on the shape of the contest success function, costly bravery can evolve in groups of any size.

Magma generation at the easternmost section of the Hellenic arc: Hf, Nd, Pb and Sr isotope geochemistry of Nisyros and Yali volcanoes (Greece)

Geochemical and petrographical studies of lavas and ignimbrites from the Quaternary Nisyros-Yali volcanic system in the easternmost part of the Hellenic arc (Greece) reveal insight into magma generating processes. A compositional gap between 61 and 68 wt.% SiO2 is recognized that coincides with the stratigraphic distinction between pre-caldera and postcaldera volcanic units. Trace element systematics support the subdivision of Nisyros and Yali volcanic units into two distinct suites of rocks. The variation of Nd and Hf present day isotope data and the fact that they are distinct from the isotope compositions of MORB rule out an origin by pure differentiation and require assimilation of a crustal component. Lead isotope ratios of Nisyros and Yali volcanic rocks support mixing of mantle material with a lower crust equivalent. However, Sr-87/Sr-86 ratios of 0.7036-0.7048 are incompatible with a simple binary mixing scenario and give low depleted mantle extraction ages (< 0.1 Ga), in contrast with Pb model ages of 0.3 Ga and Hf and Nd model ages of ca. 0.8 Ga. The budget of fluid-mobile elements Sr and Pb is likely to be dominated by abundant hydrous fluids characterised by mantle-like Sr isotope ratios. Late stage fluids probably were enriched in CO2, needed to explain the high Th concentrations. The occurrence of hydrated minerals (e.g., amphibole) in the first post-caldera unit with the lowermost Sr-87/Sr-86 ratio of 0.7036 +/- 2 can be interpreted as the result of the increased water activity in the source. The presence of two different plagioclase phenocryst generations in the first lava subsequent to the caldera-causing event is indicative for a longer storage time of this magma at a shallower level. A model capable of explaining these observations involves three evolutionary stages. First stage, assimilation of lower crustal material by a primitive magma of mantle origin (as modelled by Nd-Hf isotope systematics). This stage ended by an interruption in replenishment that led to an increase of crystallization and, hence, an increase in viscosity, suppressing eruption. During this time gap, differentiation by fractional crystallization led to enrichment of incompatible species, especially aqueous fluids, to silica depolymerisation and to a decrease in viscosity, finally enabling eruption again in the third stage. (c) 2005 Elsevier B.V. All rights reserved.