Niche, distribution and global change : modeling insights into biogeography and conservation biology Olivier Broennimann

RESUME Les évidences montrant que les changements globaux affectent la biodiversité s'accumulent. Les facteurs les plus influant dans ce processus sont les changements et destructions d'habitat, l'expansion des espèces envahissantes et l'impact des changements climatiques. Une évaluation pertinente de la réponse des espèces face à ces changements est essentielle pour proposer des mesures permettant de réduire le déclin actuel de la biodiversité. La modélisation de la répartition d'espèces basée sur la niche (NBM) est l'un des rares outils permettant cette évaluation. Néanmoins, leur application dans le contexte des changements globaux repose sur des hypothèses restrictives et demande une interprétation critique. Ce travail présente une série d'études de cas investiguant les possibilités et limitations de cette approche pour prédire l'impact des changements globaux. Deux études traitant des menaces sur les espèces rares et en danger d'extinction sont présentées. Les caractéristiques éco-géographiques de 118 plantes avec un haut degré de priorité de conservation sont revues. La prévalence des types de rareté sont analysées en relation avec leur risque d'extinction UICN. La revue souligne l'importance de la conservation à l'échelle régionale. Une évaluation de la rareté à échelle globale peut être trompeuse pour certaine espèces car elle ne tient pas en compte des différents degrés de rareté que présente une espèce à différentes échelles spatiales. La deuxième étude test une approche pour améliorer l'échantillonnage d'espèces rares en incluant des phases itératives de modélisation et d'échantillonnage sur le terrain. L'application de l'approche en biologie de la conservation (illustrée ici par le cas du chardon bleu, Eryngium alpinum), permettrait de réduire le temps et les coûts d'échantillonnage. Deux études sur l'impact des changements climatiques sur la faune et la flore africaine sont présentées. La première étude évalue la sensibilité de 227 mammifères africains face aux climatiques d'ici 2050. Elle montre qu'un nombre important d'espèces pourrait être bientôt en danger d'extinction et que les parcs nationaux africains (principalement ceux situé en milieux xériques) pourraient ne pas remplir leur mandat de protection de la biodiversité dans le futur. La seconde étude modélise l'aire de répartition en 2050 de 975 espèces de plantes endémiques du sud de l'Afrique. L'étude propose l'inclusion de méthodes améliorant la prédiction des risques liés aux changements climatiques. Elle propose également une méthode pour estimer a priori la sensibilité d'une espèce aux changements climatiques à partir de ses propriétés écologiques et des caractéristiques de son aire de répartition. Trois études illustrent l'utilisation des modèles dans l'étude des invasions biologiques. Une première étude relate l'expansion de la laitue sáuvage (Lactuca serriola) vers le nord de l'Europe en lien avec les changements du climat depuis 250 ans. La deuxième étude analyse le potentiel d'invasion de la centaurée tachetée (Centaures maculosa), une mauvaise herbe importée en Amérique du nord vers 1890. L'étude apporte la preuve qu'une espèce envahissante peut occuper une niche climatique différente après introduction sur un autre continent. Les modèles basés sur l'aire native prédisent de manière incorrecte l'entier de l'aire envahie mais permettent de prévoir les aires d'introductions potentielles. Une méthode alternative, incluant la calibration du modèle à partir des deux aires où l'espèce est présente, est proposée pour améliorer les prédictions de l'invasion en Amérique du nord. Je présente finalement une revue de la littérature sur la dynamique de la niche écologique dans le temps et l'espace. Elle synthétise les récents développements théoriques concernant le conservatisme de la niche et propose des solutions pour améliorer la pertinence des prédictions d'impact des changements climatiques et des invasions biologiques. SUMMARY Evidences are accumulating that biodiversity is facing the effects of global change. The most influential drivers of change in ecosystems are land-use change, alien species invasions and climate change impacts. Accurate projections of species' responses to these changes are needed to propose mitigation measures to slow down the on-going erosion of biodiversity. Niche-based models (NBM) currently represent one of the only tools for such projections. However, their application in the context of global changes relies on restrictive assumptions, calling for cautious interpretations. In this thesis I aim to assess the effectiveness and shortcomings of niche-based models for the study of global change impacts on biodiversity through the investigation of specific, unsolved limitations and suggestion of new approaches. Two studies investigating threats to rare and endangered plants are presented. I review the ecogeographic characteristic of 118 endangered plants with high conservation priority in Switzerland. The prevalence of rarity types among plant species is analyzed in relation to IUCN extinction risks. The review underlines the importance of regional vs. global conservation and shows that a global assessment of rarity might be misleading for some species because it can fail to account for different degrees of rarity at a variety of spatial scales. The second study tests a modeling framework including iterative steps of modeling and field surveys to improve the sampling of rare species. The approach is illustrated with a rare alpine plant, Eryngium alpinum and shows promise for complementing conservation practices and reducing sampling costs. Two studies illustrate the impacts of climate change on African taxa. The first one assesses the sensitivity of 277 mammals at African scale to climate change by 2050 in terms of species richness and turnover. It shows that a substantial number of species could be critically endangered in the future. National parks situated in xeric ecosystems are not expected to meet their mandate of protecting current species diversity in the future. The second study model the distribution in 2050 of 975 endemic plant species in southern Africa. The study proposes the inclusion of new methodological insights improving the accuracy and ecological realism of predictions of global changes studies. It also investigates the possibility to estimate a priori the sensitivity of a species to climate change from the geographical distribution and ecological proprieties of the species. Three studies illustrate the application of NBM in the study of biological invasions. The first one investigates the Northwards expansion of Lactuca serriola L. in Europe during the last 250 years in relation with climate changes. In the last two decades, the species could not track climate change due to non climatic influences. A second study analyses the potential invasion extent of spotted knapweed, a European weed first introduced into North America in the 1890s. The study provides one of the first empirical evidence that an invasive species can occupy climatically distinct niche spaces following its introduction into a new area. Models fail to predict the current full extent of the invasion, but correctly predict areas of introduction. An alternative approach, involving the calibration of models with pooled data from both ranges, is proposed to improve predictions of the extent of invasion on models based solely on the native range. I finally present a review on the dynamic nature of ecological niches in space and time. It synthesizes the recent theoretical developments to the niche conservatism issues and proposes solutions to improve confidence in NBM predictions of the impacts of climate change and species invasions on species distributions.

Robust accelerated failure time regression

Robust estimators for accelerated failure time models with asymmetric (or symmetric) error distribution and censored observations are proposed. It is assumed that the error model belongs to a log-location-scale family of distributions and that the mean response is the parameter of interest. Since scale is a main component of mean, scale is not treated as a nuisance parameter. A three steps procedure is proposed. In the first step, an initial high breakdown point S estimate is computed. In the second step, observations that are unlikely under the estimated model are rejected or down weighted. Finally, a weighted maximum likelihood estimate is computed. To define the estimates, functions of censored residuals are replaced by their estimated conditional expectation given that the response is larger than the observed censored value. The rejection rule in the second step is based on an adaptive cut-off that, asymptotically, does not reject any observation when the data are generat ed according to the model. Therefore, the final estimate attains full efficiency at the model, with respect to the maximum likelihood estimate, while maintaining the breakdown point of the initial estimator. Asymptotic results are provided. The new procedure is evaluated with the help of Monte Carlo simulations. Two examples with real data are discussed.

The relationship between Lamb weather types and long-term changes in flood frequency, River Eden, UK

Research has found that both flood magnitude and frequency in the UK may have increased over the last five decades. However, evaluating whether or not this is a systematic trend is difficult because of the lack of longer records. Here we compile and consider an extreme flood record that extends back to 1770. Since 1770, there have been 137 recorded extreme floods. However, over this period, there is not a unidirectional trend of rising extreme flood risk over time. Instead, there are clear flood-rich and flood-poor periods. Three main flood-rich periods were identified: 18731904, 19231933, and 1994 onwards. To provide a first analysis of what is driving these periods, and given the paucity of more sophisticated datasets that extend back to the 18th century, objective Lamb weather types were used. Of the 27 objective Lamb weather types, only 11 could be associated with the extreme floods during the gauged period, and only 5 of these accounted for > 80% of recorded extreme floods The importance of these five weather types over a longer timescale for flood risk in Carlisle was assessed, through calculating the proportion of each hydrological year classified as being associated with these flood-generating weather types. Two periods clearly had more than the average proportions of the year classified as one of the flood causing weather types; 19001940 and 19832007; and these two periods both contained flood-rich hydrological records. Thus, the analysis suggests that systematic organisation of the North Atlantic climate system may be manifest as periods of elevated and reduced flood risk, an observation that has major implications for analyses that assume that climatic drivers of flood risk can be either statistically stationary or are following a simple trend. Copyright (c) 2011 Royal Meteorological Society

Time for change: a roadmap to guide the implementation of the World Anti-Doping Code 2015

A medical and scientific multidisciplinary consensus meeting was held from 29 to 30 November 2013 on Anti-Doping in Sport at the Home of FIFA in Zurich, Switzerland, to create a roadmap for the implementation of the 2015 World Anti-Doping Code. The consensus statement and accompanying papers set out the priorities for the antidoping community in research, science and medicine. The participants achieved consensus on a strategy for the implementation of the 2015 World Anti-Doping Code. Key components of this strategy include: (1) sport-specific risk assessment, (2) prevalence measurement, (3) sport-specific test distribution plans, (4) storage and reanalysis, (5) analytical challenges, (6) forensic intelligence, (7) psychological approach to optimise the most deterrent effect, (8) the Athlete Biological Passport (ABP) and confounding factors, (9) data management system (Anti-Doping Administration & Management System (ADAMS), (10) education, (11) research needs and necessary advances, (12) inadvertent doping and (13) management and ethics: biological data. True implementation of the 2015 World Anti-Doping Code will depend largely on the ability to align thinking around these core concepts and strategies. FIFA, jointly with all other engaged International Federations of sports (Ifs), the International Olympic Committee (IOC) and World Anti-Doping Agency (WADA), are ideally placed to lead transformational change with the unwavering support of the wider antidoping community. The outcome of the consensus meeting was the creation of the ad hoc Working Group charged with the responsibility of moving this agenda forward.

Anatomical distribution of areas of preserved cartilage in advanced femorotibial osteoarthritis using CT arthrography (Part 1).

OBJECTIVE: To determine subregions of normal and abnormal cartilage in advanced stages of femorotibial osteoarthritis (OA) by mapping the entire femorotibial joint in a cohort of pre-total knee replacement (TKR) OA knees. DESIGN: We defined an areal subdivision of the femorotibial articular cartilage surface on CT arthrography (CTA), allowing the division of the femorotibial articular surface into multiple (up to n = 204 per knee) subregions and the comparison of the same areas between different knees. Two readers independently classified each cartilage area as normal, abnormal or non-assessable in 41 consecutive pre-TKR OA knees. RESULTS: A total of 6447 cartilage areas (from 41 knees) were considered assessable by both readers. The average proportion of preserved cartilage was lower in the medial femorotibial joint than in the lateral femorotibial joint for both readers (32.0/69.8% and 33.9/68.5% (medial/lateral) for reader 1 and 2 respectively, all P < 0.001). High frequencies of normal cartilage were observed at the posterior aspect of the medial condyle (up to 89%), and the anterior aspect of the lateral femorotibial compartment (up to 100%). The posterior aspect of the medial condyle was the area that most frequently exhibited preserved cartilage in the medial femorotibial joint, contrasting with the high frequency of cartilage lesions in the rest of that compartment. CONCLUSIONS: Cartilage at the posterior aspect of the medial condyle, and at the anterior aspect of the lateral femorotibial compartment, may be frequently preserved in advanced grades of OA.

Debris flows as a factor of hillslope evolution controlled by a continuous or a pulse process?

Flood effectiveness observations imply that two families of processes describe the formation of debris flow volume. One is related to the rainfall?erosion relationship, and can be seen as a gradual process, and one is related to additional geological/geotechnical events, those named hereafter extraordinary events. In order to discuss the hypothesis of coexistence of two modes of volume formation, some methodologies are applied. Firstly, classical approaches consisting in relating volume to catchments characteristics are considered. These approaches raise questions about the quality of the data rather than providing answers concerning the controlling processes. Secondly, we consider statistical approaches (cumulative number of events distribution and cluster analysis) and these suggest the possibility of having two distinct families of processes. However the quantitative evaluation of the threshold differs from the one that could be obtained from the first approach, but they all agree in the sense of the coexistence of two families of events. Thirdly, a conceptual model is built exploring how and why debris flow volume in alpine catchments changes with time. Depending on the initial condition (sediment production), the model shows that large debris flows (i.e. with important volume) are observed in the beginning period, before a steady-state is reached. During this second period debris flow volume such as is observed in the beginning period is not observed again. Integrating the results of the three approaches, two case studies are presented showing: (1) the possibility to observe in a catchment large volumes that will never happen again due to a drastic decrease in the sediment availability, supporting its difference from gradual erosion processes; (2) that following a rejuvenation of the sediment storage (by a rock avalanche) the magnitude?frequency relationship of a torrent can be differentiated into two phases, the beginning one with large and frequent debris flow and a later one with debris flow less intense and frequent, supporting the results of the conceptual model. Although the results obtained cannot identify a clear threshold between the two families of processes, they show that some debris flows can be seen as pulse of sediment differing from that expected from gradual erosion.

Is it the grain size or the characteristic pore size that controls the induced polarization relaxation time of clean sands and sandstones?

There is a wide range of evidence to suggest that permeability can be constrained through of induced polarization measurements. For clean sands and sandstones, current mechanistic models of induced polarization predict a relationship between the low-frequency time constant inferred from induced polarization measurements and the grain diameter. A number of observations do, however, disagree with this and indicate that the observed relaxation behavior is rather governed by the so-called dynamic pore radius L. To test this hypothesis, we have developed a set of new scaling relationships, which allow the relaxation time to be computed from the pore size and the permeability to be computed from both the Cole-Cole time constant and the formation factor. Moreover, these new scaling relationships can be also used to predict the dependence of the Cole-Cole time constant as a function of the water saturation under unsaturated conditions. Comparative tests of the proposed new relationships with regard to various published experimental results for saturated clean sands and sandstones as well as for partially saturated clean sandstones, do indeed confirm that the dynamic pore radius L is a much more reliable indicator of the observed relaxation behavior than grain-size-based models.

Tuning In to Sound: Frequency-Selective Attentional Filter in Human Primary Auditory Cortex.

Cocktail parties, busy streets, and other noisy environments pose a difficult challenge to the auditory system: how to focus attention on selected sounds while ignoring others? Neurons of primary auditory cortex, many of which are sharply tuned to sound frequency, could help solve this problem by filtering selected sound information based on frequency-content. To investigate whether this occurs, we used high-resolution fMRI at 7 tesla to map the fine-scale frequency-tuning (1.5 mm isotropic resolution) of primary auditory areas A1 and R in six human participants. Then, in a selective attention experiment, participants heard low (250 Hz)- and high (4000 Hz)-frequency streams of tones presented at the same time (dual-stream) and were instructed to focus attention onto one stream versus the other, switching back and forth every 30 s. Attention to low-frequency tones enhanced neural responses within low-frequency-tuned voxels relative to high, and when attention switched the pattern quickly reversed. Thus, like a radio, human primary auditory cortex is able to tune into attended frequency channels and can switch channels on demand.

Ambulatory assessment of 3D ground reaction force using plantar pressure distribution.

This study aimed to use the plantar pressure insole for estimating the three-dimensional ground reaction force (GRF) as well as the frictional torque (T(F)) during walking. Eleven subjects, six healthy and five patients with ankle disease participated in the study while wearing pressure insoles during several walking trials on a force-plate. The plantar pressure distribution was analyzed and 10 principal components of 24 regional pressure values with the stance time percentage (STP) were considered for GRF and T(F) estimation. Both linear and non-linear approximators were used for estimating the GRF and T(F) based on two learning strategies using intra-subject and inter-subjects data. The RMS error and the correlation coefficient between the approximators and the actual patterns obtained from force-plate were calculated. Our results showed better performance for non-linear approximation especially when the STP was considered as input. The least errors were observed for vertical force (4%) and anterior-posterior force (7.3%), while the medial-lateral force (11.3%) and frictional torque (14.7%) had higher errors. The result obtained for the patients showed higher error; nevertheless, when the data of the same patient were used for learning, the results were improved and in general slight differences with healthy subjects were observed. In conclusion, this study showed that ambulatory pressure insole with data normalization, an optimal choice of inputs and a well-trained nonlinear mapping function can estimate efficiently the three-dimensional ground reaction force and frictional torque in consecutive gait cycle without requiring a force-plate.

Distribution of metals exposure and associations with cardiometabolic risk factors in the "Modeling the Epidemiologic Transition Study".

BACKGROUND: Metals are known endocrine disruptors and have been linked to cardiometabolic diseases via multiple potential mechanisms, yet few human studies have both the exposure variability and biologically-relevant phenotype data available. We sought to examine the distribution of metals exposure and potential associations with cardiometabolic risk factors in the "Modeling the Epidemiologic Transition Study" (METS), a prospective cohort study designed to assess energy balance and change in body weight, diabetes and cardiovascular disease risk in five countries at different stages of social and economic development. METHODS: Young adults (25-45 years) of African descent were enrolled (N = 500 from each site) in: Ghana, South Africa, Seychelles, Jamaica and the U.S.A. We randomly selected 150 blood samples (N = 30 from each site) to determine concentrations of selected metals (arsenic, cadmium, lead, mercury) in a subset of participants at baseline and to examine associations with cardiometabolic risk factors. RESULTS: Median (interquartile range) metal concentrations (μg/L) were: arsenic 8.5 (7.7); cadmium 0.01 (0.8); lead 16.6 (16.1); and mercury 1.5 (5.0). There were significant differences in metals concentrations by: site location, paid employment status, education, marital status, smoking, alcohol use, and fish intake. After adjusting for these covariates plus age and sex, arsenic (OR 4.1, 95% C.I. 1.2, 14.6) and lead (OR 4.0, 95% C.I. 1.6, 9.6) above the median values were significantly associated with elevated fasting glucose. These associations increased when models were further adjusted for percent body fat: arsenic (OR 5.6, 95% C.I. 1.5, 21.2) and lead (OR 5.0, 95% C.I. 2.0, 12.7). Cadmium and mercury were also related with increased odds of elevated fasting glucose, but the associations were not statistically significant. Arsenic was significantly associated with increased odds of low HDL cholesterol both with (OR 8.0, 95% C.I. 1.8, 35.0) and without (OR 5.9, 95% C.I. 1.5, 23.1) adjustment for percent body fat. CONCLUSIONS: While not consistent for all cardiometabolic disease markers, these results are suggestive of potentially important associations between metals exposure and cardiometabolic risk. Future studies will examine these associations in the larger cohort over time.

Building the niche through time: using 13,000 years of data to predict the effects of climate change on three tree species in Europe

Aim Species distribution models (SDMs) based on current species ranges underestimate the potential distribution when projected in time and/or space. A multi-temporal model calibration approach has been suggested as an alternative, and we evaluate this using 13,000 years of data. Location Europe. Methods We used fossil-based records of presence for Picea abies, Abies alba and Fagus sylvatica and six climatic variables for the period 13,000 to 1000yr bp. To measure the contribution of each 1000-year time step to the total niche of each species (the niche measured by pooling all the data), we employed a principal components analysis (PCA) calibrated with data over the entire range of possible climates. Then we projected both the total niche and the partial niches from single time frames into the PCA space, and tested if the partial niches were more similar to the total niche than random. Using an ensemble forecasting approach, we calibrated SDMs for each time frame and for the pooled database. We projected each model to current climate and evaluated the results against current pollen data. We also projected all models into the future. Results Niche similarity between the partial and the total-SDMs was almost always statistically significant and increased through time. SDMs calibrated from single time frames gave different results when projected to current climate, providing evidence of a change in the species realized niches through time. Moreover, they predicted limited climate suitability when compared with the total-SDMs. The same results were obtained when projected to future climates. Main conclusions The realized climatic niche of species differed for current and future climates when SDMs were calibrated considering different past climates. Building the niche as an ensemble through time represents a way forward to a better understanding of a species' range and its ecology in a changing climate.

Systèmes de délivrance des médicaments pour le segment antérieur: bases fondamentales et applications cliniques [Drug delivery systems to target the anterior segment of the eye: fundamental bases and clinical applications].

The development of new drug delivery systems to target the anterior segment of the eye may offer many advantages: to increase the biodisponibility of the drug, to allow the penetration of drug that cannot be formulated as solutions, to obtain constant and sustained drug release, to achieve higher local concentrations without systemic effects, to target more specifically one tissue or cell type, to reduce the frequency of instillation and therefore increase the observance and comfort of the patient while reducing side effects of frequent instillation. Several approaches are developed, aiming to increase the corneal contact time by modified formulation or reservoir systems, or by increasing the tissue permeability using iontophoresis. To date, no ocular drug delivery system is ideal for all purposes. To maximize treatment efficacy, careful evaluation of the specific pathological condition, the targeted Intraocular tissue and the location of the most severe pathology must be made before selecting the method of delivery most suitable for each individual patient.

NR2A at CA1 synapses is obligatory for the susceptibility of hippocampal plasticity to sleep loss.

A loss in the necessary amount of sleep alters expression of genes and proteins implicated in brain plasticity, but key proteins that render neuronal circuits sensitive to sleep disturbance are unknown. We show that mild (4-6 h) sleep deprivation (SD) selectively augmented the number of NR2A subunits of NMDA receptors on postsynaptic densities of adult mouse CA1 synapses. The greater synaptic NR2A content facilitated induction of CA3-CA1 long-term depression in the theta frequency stimulation range and augmented the synaptic modification threshold. NR2A-knock-out mice maintained behavioral response to SD, including compensatory increase in post-deprivation resting time, but hippocampal synaptic plasticity was insensitive to sleep loss. After SD, the balance between synaptically activated and slowly recruited NMDA receptor pools during temporal summation was disrupted. Together, these results indicate that NR2A is obligatory for the consequences of sleep loss on hippocampal synaptic plasticity. These findings could advance pharmacological strategies aiming to sustain hippocampal function during sleep restriction.

Frequency and chronology of extraintestinal manifestations in the Swiss Inflammatory Bowel Disease Cohort

Background: Data on the frequency of extraintestinal manifestations (EIM) in inflammatory bowel disease (IBD) is scarce, especially the one evaluating the time of occurrence of the EIM relative to IBD diagnosis. Aim: To assess the type and frequency of EIM in IBD patients and to evaluate when EIMs occur relative to IBD diagnosis. Methods: Analysis of data from the Swiss Inflammatory Bowel Disease Cohort (SIBDCS) which collects data on a large sample of IBD patients from hospitals and private practices across Switzerland starting in 2005. While parametric data are shown as mean ± SD, non-parametric data are presented as median and interquartile range (IQR). Results: A total of 1143 patients were analyzed (572 (50%) female, mean age 42.1 ± 14.4 years): 629 (55%) with Crohn's disease (CD), 501 (44%) with ulcerative colitis (UC), and 13 (1%) with indeterminate colitis (IC). Of 1143 patients, 374 (32.7%) presented with EIM (65% with CD, 33% with UC, 2% with IC). Of those patients suffering from EIMs, 37.4% presented with one, 41.7% with two, 12.4% with three, 5.3% with four, and 3.2% with five EIM during their lifetime. The IBD patients initially presented with the following EIMs: peripheral arthritis (PA) 63.4%, ankylosing spondylitis (AS) 8.1%, primary sclerosing cholangitis (PSC) 6.0%, uveitis 5.7%, oral aphthosis 5.7%, erythema nodosum (EN) 5.0%, pyoderma gangrenosum 1.8%, psoriasis 0.7%. While 92.9% of EIM occurred once IBD diagnosis was established (median 72 months, IQR 9-147 months, p < 0.001), 7.1% of EIMs preceded IBD diagnosis (median time 28 months before IBD diagnosis, IQR 7-60 months). Over a course of a lifetime, IBD patients presented with the following EIM (total exceeds 100 due to potential presence of multiple EIM): peripheral arthritis 69.3%, oral aphthosis 23%, ankylosing spondylitis 19.4%, uveitis 15.5%, erythema nodosum 14.5%, PSC 7.8%, pyoderma gangrenosum 6%, psoriasis 2.8%. Conclusion: EIMs frequently occur in a lifetime of IBD patients. The vast majority of patients present with EIMs once IBD diagnosis has been established. IBD patients most often present with peripheral arthritis, ankylosing spondylitis and PSC as their first EIM. However, peripheral arthritis, oral aphthosis, and ankylosing spondylitis are the most common EIMs in a lifetime of IBD patients.

Phosphorus burial in the ocean over glacial-interglacial time scales

The role of nutrients, such as phosphorus (P), and their impact on primary productivity and the fluctuations in atmospheric CO2 over glacial-interglacial periods are intensely debated. Suggestions as to the importance of P evolved from an earlier proposal that P actively participated in changing productivity rates and therefore climate change, to most recent ones that changes in the glacial ocean inventory of phosphorus were important but not influential if compared to other macronutrients, such as nitrate. Using new data coming from a selection of ODP sites, we analyzed the distribution of oceanic P sedimentary phases and calculate reactive P burial fluxes, and we show how P burial fluxes changed over the last glacial-interglacial period at these sites. Concentrations of reactive P are generally lower during glacial times, while mass accumulation rates (MAR) of reactive P show higher variability. If we extrapolate for the analyzed sites, we may assume that in general glacial burial fluxes of reactive P are lower than those during interglacial periods by about 8%, because the lack of burial of reactive P on the glacial shelf reduced in size, was apparently not compensated by burial in other regions of the ocean. Using the calculated changes in P burial, we evaluate their possible impact on the phosphate inventory in the world oceans. Using a simple mathematical approach, we find that these changes alone could have increased the phosphate inventory of glacial ocean waters by 17-40% compared to interglacial stages. Variations in the distribution of sedimentary P phases at the investigated sites seem to indicate that at the onset of interglacial stages, shallower sites experienced an increase in reactive P concentrations, which seems to point to P-richer waters at glacial terminations. All these findings would support the Shelf-Nutrient Hypothesis, which assumes that during glacial low stands nutrients are transferred from shallow sites to deep sea with possible feedback on the carbon cycle.